| 83 |
|
Xi(4, 4) = Xirr; |
| 84 |
|
Xi(5, 5) = Xirr; |
| 85 |
|
|
| 86 |
< |
const RealType convertConstant = 6.023; //convert poise.angstrom to amu/fs |
| 87 |
< |
Xi *= convertConstant; |
| 86 |
> |
Xi *= OOPSEConstant::viscoConvert; |
| 87 |
|
XiCopy = Xi; |
| 88 |
|
|
| 89 |
|
invertMatrix(XiCopy, D); |
| 90 |
< |
RealType kt = OOPSEConstant::kB * temperature; |
| 91 |
< |
D *= kt; |
| 93 |
< |
Xi *= OOPSEConstant::kb * temperature; |
| 90 |
> |
RealType kt = OOPSEConstant::kb * temperature; // in kcal mol^-1 |
| 91 |
> |
D *= kt; // now in angstroms^2 fs^-1 (at least for Trans-trans) |
| 92 |
|
|
| 93 |
|
HydroProp* hprop = new HydroProp(V3Zero, Xi, D); |
| 94 |
|
|
