| 35 |
|
* |
| 36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
| 37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
| 38 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
| 38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
| 39 |
|
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
| 40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
| 41 |
|
*/ |
| 157 |
|
|
| 158 |
|
// Calculate lattice constant (in Angstroms) |
| 159 |
|
|
| 160 |
< |
RealType avgMass = getMolMass(oldInfo->getMoleculeStamp(0), |
| 161 |
< |
oldInfo->getForceField()); |
| 160 |
> |
RealType avgMass = MoLocator::getMolMass(oldInfo->getMoleculeStamp(0), |
| 161 |
> |
oldInfo->getForceField()); |
| 162 |
|
|
| 163 |
|
latticeConstant = pow(rhoConvertConst * nMolPerCell * avgMass / density, |
| 164 |
|
(RealType)(1.0 / 3.0)); |
| 199 |
|
|
| 200 |
|
createMdFile(inputFileName, outputFileName, nSites); |
| 201 |
|
|
| 202 |
< |
if (oldInfo != NULL) |
| 203 |
< |
delete oldInfo; |
| 202 |
> |
delete oldInfo; |
| 203 |
|
|
| 204 |
|
// We need to read in the new SimInfo object, then Parse the |
| 205 |
|
// md file and set up the system |