Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials

Abstract

The effects of the finite size of the simulation box in equilibrium molecular dynamics simulations are investigated for prototypical superionic conductors of different types, namely the fluorite-structure materials PbF2, CaF2, and UO2 (type-II), and the α phase of AgI (type I). Largely validated empirical force-fields are employed to run ns-long simulations and extract general trends for several properties, at increasing size and in a wide temperature range. This work shows that, for the considered type-II superionic conductors, the diffusivity dramatically depends on the system size and that the superionic regime is shifted to larger temperatures in smaller cells. Furthermore, only simulations of several hundred atoms are able to capture the experimentally-observed, characteristic change in the activation energy of the diffusion process, occurring at the order-disorder transition to the superionic regime. Finite-size effects on ion diffusion are instead much weaker in α-AgI. The thermal conductivity is found generally smaller for smaller cells, where the temperature-independent (Allen-Feldman) regime is also reached at significantly lower temperatures. The finite-size effects on the thermal motion of the non-mobile ions composing the solid matrix follow the simple law which holds for solids.

arXiv comment10 pages, 8 figures

References in corpus (20)

• Canonical sampling through velocity-rescaling
• Unified theory of thermal transport in crystals and disordered solids
• Physics-inspired structural representations for molecules and materials
• Materials Cloud, a platform for open computational science
• Modeling heat transport in crystals and glasses from a unified lattice-dynamical approach
• High-throughput computational screening for solid-state Li-ion conductors
• Deep Potential generation scheme and simulation protocol for the Li10GeP2S12-type superionic conductors
• Uranus evolution models with simple thermal boundary layers
• Heat transport in liquid water from first-principles and deep-neural-network simulations
• Accurate thermal conductivities from optimally short molecular dynamics simulations
• Effect of Non-Adiabatic Thermal Profiles on the Inferred Compositions of Uranus and Neptune
• Modeling lithium-ion solid-state electrolytes with a pinball model
• Theory and numerical simulation of heat transport in multi-component systems
• High-precision molecular dynamics simulation of UO2-PuO2: Anion self-diffusion in UO2
• Temperature- and vacancy-concentration-dependence of heat transport in Li${}_3$ClO from multi-method numerical simulations
• Heat and charge transport in H${}_2$O at ice-giant conditions from ab initio molecular dynamics simulations
• SporTran: a code to estimate transport coefficients from the cepstral analysis of (multivariate) current time series
• Invariance principles in the theory and computation of transport coefficients
• Heat transport in insulators from ab initio Green-Kubo theory
• Electrolytes for Li-ion all-solid-state batteries: a first-principles comparative study of Li10GeP2O12 and Li10GeP2S12 in the LISICON and LGPS phases