Energy localization and excess fluctuations from long-range interactions in equilibrium molecular dynamics

Ralph V. Chamberlin, Vladimiro Mujica, Sergei Izvekov and James P. Larentzos

Physica A: Statistical Mechanics and its Applications, 2020, vol. 540, issue C

Abstract: Molecular dynamics (MD) simulations of standard systems of interacting particles (“atoms”) give excellent agreement with the equipartition theorem for the average energy, but we find that these simulations exhibit finite-size effects in the dynamics that cause local fluctuations in energy to deviate significantly from the analogous energy fluctuation relation (EFR). We have made a detailed analysis of Lennard-Jones atoms to track the origin of such unphysical fluctuations, which must be corrected for an appropriate description of the statistical mechanics of small systems, especially at low temperatures (T). Similar behavior is found in a model of nitromethane at higher T. The main conclusion of our study is that systems separated into nanometer-sized “blocks” inside much larger simulations exhibit excess fluctuations in potential energy (pe) that diverge inversely proportional to T in a manner that is strongly dependent on the range of interaction. Specifically, at low T with long-range interactions pe fluctuations exceed the EFR by at least an order of magnitude, dropping abruptly to below the EFR when interactions include only 1st-neighbor atoms. Thus, excess pe fluctuations cannot be due to simple surface effects from the robustly harmonic 1st-neighbor interactions, nor from any details in the simulations or analysis. A simplistic model that includes 2nd-neighbor interactions matches the excess fluctuations, but only if the 2nd-neighbor energy is not included in Boltzmann’s factor, attributable to anharmonic effects and energy localization. Empirically, adding 2nd-neighbor interactions greatly increases the width of the pe distribution, suggesting that Anderson localization may play a role. Characterizing energy correlations as a function of time and distance reveals that excess pe fluctuations in each block coincide with negative pe correlations between neighboring blocks, whereas reduced pe fluctuations coincide with positive pe correlations. Indeed, anomalous pe fluctuations in small systems can be quantified using a net local energy in Boltzmann’s factor that includes pe from the surrounding shell of similarly small systems, or equivalently an effective local temperature. Our results and analysis elucidate the source of non-Boltzmann fluctuations, and the need to include mesoscopic thermal effects from the local environment for a consistent theoretical description of the equilibrium fluctuations in MD simulations of standard models with long-range interactions.

Keywords: Statistical physics; MD simulations; Fluctuations; Nanoscale dynamics; Specific heat; Localization (search for similar items in EconPapers)
Date: 2020
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0378437119318126
Full text for ScienceDirect subscribers only. Journal offers the option of making the article available online on Science direct for a fee of $3,000

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:phsmap:v:540:y:2020:i:c:s0378437119318126

DOI: 10.1016/j.physa.2019.123228

Access Statistics for this article

Physica A: Statistical Mechanics and its Applications is currently edited by K. A. Dawson, J. O. Indekeu, H.E. Stanley and C. Tsallis

More articles in Physica A: Statistical Mechanics and its Applications from Elsevier
Bibliographic data for series maintained by Catherine Liu ().