Non-Markovian dynamics in ice nucleation
Pablo Montero de Hijes, Sebastian Falkner, Christoph Dellago
arXiv:2512.01479·cond-mat.soft·Published 2025-12-01
In simulation studies of crystallisation, the size of the largest crystalline nucleus is often used as a reaction coordinate to monitor the progress of the nucleation process. Here, we investigate, for the case of homogeneous ice nucleation, whether the nucleus size exhibits Markovian dynamics, as assumed in classical nucleation theory. Using 300 independent nucleation trajectories generated by molecular dynamics, we evaluate the mean recurrence time required to reach selected values of the largest nucleus size. Early recurrences consistently take longer than later ones, revealing a clear history dependence and thus non-Markovian dynamics. To identify the slow modes underlying this behaviour, we analyse several structural descriptors of the nucleus, observing subtle but systematic differences between nuclei at early and late recurrences. By training a neural network on 2,700 short trajectories to learn the committor, we identify relevant collective variables. Based on these features, symbolic regression provides a compact approximation of the committor, i.e., an improved reaction coordinate, which we subsequently test for Markovianity.
TopicsQuantum Chemistry & Force Fields
Tagsmolecular-dynamics
arXiv categoriescond-mat.soft
arXiv abstract pagePDF