Molecular dynamics simulations of compressible ions

Abstract

A representation of the short-range repulsion energy in an ionic system is described which allows for the fact that an ion may be compressed by its neighbours. The total energy of the system is expressed in a pairwise additive form, but the interionic interactions have a many-body character. The form of this representation and the parameters required to represent MgO and CaO are obtained from recent ab-initio electronic structure calculations. The fact that the representation is transferable between crystals with different coordination number is demonstrated by direct comparison with ab-initio results on the different crystal types. Comparison with experimental results on the equation of state of different isomorphs and on the location of the pressure of the transition between them confirms the accuracy of the ab-initio results and of the potential derived from them in representing perfect crystal properties. A computationally efficient molecular dynamics (MD) scheme may be derived for this representation. The additional degrees of freedom which represent the varying ionic radii are constrained to their adiabatic values in the course of the simulation by an adaptation of Car and Parrinello’s method. The MD scheme is used to examine whether an ab-initio parameterized potential model which allows for the spherical compression of an oxide ion by its neighbours and for dipole polarization effects is a sufficiently good representation of the interactions in MgO to allow an accurate calculation of the phonon dispersion curves.

Evidence of dynamical local scale distribution heterogeneities in CsF-AF (A=Li, Na, K and Rb) molten salts

“…67 In addition, the accuracy of the molecular dynamics simulations necessitates taking account the ion polarizability. 69,70,71 In molten fluorides, the variations of the cations chemical shift have often provided crucial and more sensitive information to understand the local structure of the melt. 58,40,72 For this purpose, the chemical shift is measured for several crystalline compounds which structure is perfectly defined.…”

Including many-body effects in models for ionic liquids

“…[24][25][26] The potential is best described as the sum of four different components: charge-charge, dispersion, overlap repulsion and polarization.…”

The dispersion interaction between quantum mechanics and effective fragment potential molecules

“…First, numerous theoretical studies suggest that, in order to compare with experimental or high-level ab initio calculations, a damping function must be used with the dispersion energy term. [28][29][30][31][32][33][34][35] The damping function, F v 6 , is employed to account for both exchange-dispersion and charge penetration effects, which predominate at short intermolecular separations. Additionally, as with EFP2-EFP2 dispersion, 6 the total EFP2-AI dispersion energy expansion in Eq.…”