Sensitivity of Nucleation Thermodynamics and Kinetics to the Treatment of Long-Range Interactions
Abstract
Nucleation rates are exponentially sensitive to the thermodynamic driving force and can therefore depend strongly on the treatment of long-range intermolecular interactions. Here, using the Lennard--Jones (LJ) system as a benchmark, we combine molecular dynamics (MD) simulations, jumpy forward-flux sampling (jFFS), and free-energy calculations to quantify the effect of potential truncation on melting thermodynamics, homogeneous crystal nucleation kinetics, and computational cost. Within the cutoff-radius range 2.5σ≤ rc≤ 6σ, the melting temperature at zero pressure varies by approximately 11%, while the nucleation rate changes by approximately ten orders of magnitude. By invoking classical nucleation theory (CNT), we show that this pronounced kinetic sensitivity originates primarily from cutoff-induced changes in the chemical potential difference between the liquid and crystalline phases. Building on this observation, we develop a CNT-based framework for extrapolating finite-cutoff rates to the full-potential limit and for estimating the expected rate deviations at other cutoff radii and temperatures. These findings also provide a systematic basis for cutoff selection: the optimal cutoff should minimize computational cost while keeping the deviation from the full-potential rate within acceptable bounds. At kT/ε=0.5, rc=4σ provides a reasonable compromise according to these criteria. We further demonstrate that conventional homogeneous tail corrections do not offer a reliable alternative, as they cannot consistently account for the liquid, crystalline, and interfacial environments present during nucleation. Our findings highlight the need to specify and validate the truncation scheme as an integral component of force-field development in simulations of nucleation and other interfacial phase transitions within inhomogeneous environments.
Related
- Substrate-induced crystallisation in polymers: Exceptional impact of prefreezing versus heterogeneous nucleation on kinetics and morphology
- Parameter Recovery and Noise Sensitivity of Physics-Informed Neural Networks for Michaelis–Menten Kinetics: A Benchmark Against Classical Nonlinear Least Squares
- QCD thermodynamics through the crossover as a gas of confining strings with repulsive interactions
- Projected sensitivity of a spinor Bose-Einstein condensate comagnetometer to proton monopole-dipole interactions
- Consistent treatment of flavour-specific neutrino self-interactions in cosmology
BibTeX
@misc{vargas2026sensitivity,
title = {{Sensitivity of Nucleation Thermodynamics and Kinetics to the Treatment of Long-Range Interactions}},
author = {Fernanda Sulantay Vargas and Kimia Sinaeian and Amir Haji-Akbari},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.19646},
eprint = {2609.19646},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.19646},
}
RIS
TY - UNPB
TI - Sensitivity of Nucleation Thermodynamics and Kinetics to the Treatment of Long-Range Interactions
AU - Fernanda Sulantay Vargas
AU - Kimia Sinaeian
AU - Amir Haji-Akbari
PY - 2026
DA - 2026/09/17
PB - arXiv
DO - 10.48550/arxiv.2609.19646
UR - https://arxiv.org/abs/2609.19646
AB - Nucleation rates are exponentially sensitive to the thermodynamic driving force and can therefore depend strongly on the treatment of long-range intermolecular interactions. Here, using the Lennard--Jones (LJ) system as a benchmark, we combine molecular dynamics (MD) simulations, jumpy forward-flux sampling (jFFS), and free-energy calculations to quantify the effect of potential truncation on melting thermodynamics, homogeneous crystal nucleation kinetics, and computational cost. Within the cutoff-radius range 2.5σ≤ r_c≤ 6σ, the melting temperature at zero pressure varies by approximately 11%, while the nucleation rate changes by approximately ten orders of magnitude. By invoking classical nucleation theory (CNT), we show that this pronounced kinetic sensitivity originates primarily from cutoff-induced changes in the chemical potential difference between the liquid and crystalline phases. Building on this observation, we develop a CNT-based framework for extrapolating finite-cutoff rates to the full-potential limit and for estimating the expected rate deviations at other cutoff radii and temperatures. These findings also provide a systematic basis for cutoff selection: the optimal cutoff should minimize computational cost while keeping the deviation from the full-potential rate within acceptable bounds. At kT/ε=0.5, r_c=4σ provides a reasonable compromise according to these criteria. We further demonstrate that conventional homogeneous tail corrections do not offer a reliable alternative, as they cannot consistently account for the liquid, crystalline, and interfacial environments present during nucleation. Our findings highlight the need to specify and validate the truncation scheme as an integral component of force-field development in simulations of nucleation and other interfacial phase transitions within inhomogeneous environments.
ER -