Exploration of NVE classical trajectories as a tool for molecular crystal structure prediction, with tests on ice polymorphs.

Following an initial Communication [Buch et al., J. Chem. Phys. 123, 051108 (2005)], a new molecular-dynamics-based approach is explored to search for candidate crystal structures of molecular solids corresponding to minima of the enthalpy. The approach is based on the observation of phase transitions in an artificial periodic system with a small unit cell and relies on the existence of an optimal energy range for observing freezing to low-lying minima in the course of classical trajectories. Tests are carried out for O structures of nine H2O-ice polymorphs. NVE trajectories for a range of preimposed box shapes display freezing to the different crystal polymorphs whenever the box dimensions approximate roughly the appropriate unit cell; the exception is ice II for which freezing requires unit cell dimensions close to the correct ones. In an alternate version of the algorithm, an initial box shape is picked at random and subsequently readjusted at short trajectory intervals by enthalpy minimization. Tests reveal the existence of ice forms which are "difficult" and "easy" to locate in this way. The former include ice IV, which is also difficult to crystallize experimentally from the liquid, and ice II, which does not interface with the liquid in the phase diagram. On the other hand, the latter crystal search procedure located successfully the remaining seven ice polymorphs, including ice V, which corresponds to the most complicated structure of all ice phases, with a monoclinic cell of 28 molecules.