Isothermal bulk modulus and its first pressure derivative of NaCl at high pressure and high temperature

The isothermal bulk modulus and its first pressure derivative of NaCl are investigated using the classical molecular dynamics method and the quasi-harmonic Debye model. To ensure faithful molecular dynamics simulations, two types of potentials, the shell-model (SM) potential and the two-body rigid-ion Born-Mayer-Huggins-Fumi-Tosi (BMHFT) potential, are fully tested. Compared with the SM potential based simulation, the molecular dynamics simulation with the BMHFT potential is very successful in reproducing accurately the measured bulk modulus of NaCl. Particular attention is paid to the prediction of the isothermal bulk modulus and its first pressure derivative using the reliable potential and to the comparison of the SM and the BMHFT potentials based molecular dynamics simulations with the quasi-harmonic Debye model. The properties of NaCl in the pressure range of 0-30 GPa at temperatures up to the melting temperature of 1050 K are investigated.     

Isothermal bulk modulus and its first pressure derivative of NaCl at high pressure and high temperature

The isothermal bulk modulus and its first pressure derivative of NaCl are investigated using the classical molecular dynamics method and the quasi-harmonic Debye model.To ensure faithful molecular dynamics simulations,two types of potentials,the shell-model(SM) potential and the two-body rigid-ion Born-Mayer-Huggins-Fumi-Tosi(BMHFT) potential,are fully tested.Compared with the SM potential based simulation,the molecular dynamics simulation with the BMHFT potential is very successful in reproducing accurately the measured bulk modulus of NaCl.Particular attention is paid to the prediction of the isothermal bulk modulus and its first pressure derivative using the reliable potential and to the comparison of the SM and the BMHFT potentials based molecular dynamics simulations with the quasi-harmonic Debye model.The properties of NaCl in the pressure range of 0-30 GPa at temperatures up to the melting temperature of 1050 K are investigated.     

Comparative investigations of the P-V-T relationship of NaCl at high pressure and temperature

Two different potential models of molecular dynamics (MD) simulations have been applied to investigate the pressure-volume-temperature (P-V-T) relationship and lattice parameter of NaCl under high pressure and temperature. The first one is the shell model (SM) potentials in which due to the short-range interaction pairs of ions are moved together as is the case in polarization of a crystal due to the motion of the positive and negative ions, and the second one is the two-body rigid-ion Born-Mayer-Huggins-Fumi-Tosi (BMHFT) potentials with full treatment of long-range Coulomb forces. The P-V relationship at 300 K, T-V relationship at zero pressure, and lattice parameter a, have been obtained and compared with the available experimental data and other theoretical results. Compared with SM potentials, the MD simulation with BMHFT potentials is very successful in reproducing accurately the measured volumes of NaCl. At an extended pressure and temperature ranges, P-V relationship under different isotherms at selected temperatures, T-V relationship under different pressures, and lattice parameter a have also been predicted. The properties of NaCl are summarized in the pressure range 0-30 GPa and the temperature up to 2000 K.     

Thermoelastic parameter αKT of sodium chloride at high pressure and high temperature

Two different potential models to the molecular dynamics (MD) simulations have been applied to investigate the thermoelastic parameter αK_T of sodium chloride (NaCl) under high pressure and high temperature. The first one is the shell model (SM) potential that due to the short-range interaction when pairs of ions are moved together as is the case in that polarization of a crystal due to the motion of the positive and negative ions, and the second one is the two-body rigid-ion Born–Mayer–Huggins–Fumi–Tosi (BMHFT) potential with full treatment of long-range Coulomb forces. Particular attention is paid to the comparison of the SM- and BMHFT-MD simulations with the Debye model for the first time, and this model combines with ab initio calculations within local density approximation (LDA) and generalized gradient approximation (GGA) using ultrasoft pseudopotentials and a plane-wave basis in the framework of density functional theory (DFT), and it takes into account the phononic effects within the quasi-harmonic approximation. Note that the MD calculated volumes using SM model is somewhat larger than both the DFT and experimental volumes despite not considering the temperature effect. Compared with SM potential, the MD simulated 300 K isotherm of NaCl with BMHFT potential is very successful in reproducing accurately the measured volumes and the GGA calculated volumes. Generally, it is found that there exist minor differences between the LDA and GGA computed the thermoelastic parameter αK_T of NaCl, with both average results giving good agreement with SM-MD simulations. At an extended pressure and temperature ranges, the variation of thermoelastic parameter αK_T which play a central role in the formulation of approximate equations of state has also been predicted. The properties of NaCl are summarized in the pressure range of 0–300 kbar and the temperature up to 2000 K.     

Melting and Grüneisen parameters of NaCl at high pressure

The Buckingham potential has been employed to simulate the melting and thermodynamic parameters of sodium chloride (NaCl) using the molecular dynamics (MD) method. The constant-volume heat capacity and Grüneisen parameters have been obtained in a wide range of temperatures. The calculated thermodynamic parameters are found to be in good agreement with the available experimental data. The NaCl melting simulations appear to validate the interpretation of superheating of the solid in the one-phase MD simulations. The melting curve of NaCl is compared with the experiments and other calculations at pressure 0－30GPa range.     

Effect of pressure on the thermal expansion of MgO up to 200~GPa

Constant temperature and pressure molecular dynamics (MD) simulations are performed to investigate the thermal expansivity of MgO at high pressure, by using effective pair-wise potentials which consist of Coulomb, dispersion, and repulsion interactions that include polarization effects through the shell model (SM). In order to take into account non-central forces in crystals, the breathing shell model (BSM) is also introduced into the MD simulation. We present a comparison between the volume thermal expansion coefficient α dependences of pressure P at 300 and 2000~K that are obtained from the SM and BSM potentials and those derived from other experimental and theoretical methods in the case of MgO. Compared with the results obtained by using the SM potentials, the MD results obtained by using BSM potentials are more compressible. In an extended pressure and temperature range, the α value is also predicted. The properties of MgO in a pressure range of 0--200~GPa at temperatures up to 3500~K are summarized.     

Alchemical molecular dynamics for inverse design

We present a molecular dynamics (MD) implementation of an extended statistical mechanical ensemble that includes ‘alchemical’ degrees of freedom describing particle attributes as thermodynamic variables. We demonstrate the use of this alchemical MD method in inverse design simulations of particles interacting via the Oscillating Pair Potential (OPP) and the Lennard–Jones–Gauss potential (LJG) – two general, previously studied models for which phase diagrams are known. We show that alchemical MD can quickly and efficiently optimise pair potentials for target structures within a specified design space in the low-temperature regime, where internal energy adequately represents the features of the alchemical free energy landscape. We show that alchemical MD can be also used to inversely design pair potentials to achieve target materials properties (here, bulk modulus) directly, without explicit knowledge of the structure–property relationship. Alchemical MD can easily be generalised and applied to any target materials properties or structures and used with any differentiable interaction potential.     

Dynamics of ions in the selectivity filter of the KcsA channel

The statistical and dynamical properties of ions in the selectivity filter of the KcsA ion channel are considered on the basis of molecular dynamics (MD) simulations of the KcsA protein embedded in a lipid membrane surrounded by an ionic solution. A new approach to the derivation of a Brownian dynamics (BD) model of ion permeation through the filter is discussed, based on unbiased MD simulations. It is shown that depending on additional assumptions, ion’s dynamics can be described either by under-damped Langevin equation with constant damping and white noise or by Langevin equation with a fractional memory kernel. A comparison of the potential of the mean force derived from unbiased MD simulations with the potential produced by the umbrella sampling method demonstrates significant differences in these potentials. The origin of these differences is an open question that requires further clarifications.     

Experimental and theoretical study of the far-infrared spectra of HCl dissolved in liquid Ar,Kr, and Xe

New experimental results are reported on the relative absorption intensity distribution in the FIR spectra of HCl dissolved in liquefied Ar, Kr, and Xe at several temperatures along the liquid—vapour coexistence curve. These are treated further by applying a previously developed quantum-statistical spectral theory, which accounts for the line mixing and memory effects. Theoretical spectra are given in terms of the anisotropic potential time autocorrelation functions obtained from classical MD simulations using several empirical analytical potentials with density-adjusted well depths. Globally fair agreement between the theoretical and experimental spectra is demonstrated, except in the high frequency wings, where the theory underestimates the observed intensities. The choice of a particular radial form for the anisotropic HCl/RG potentials is found to be not critical for reproducing the experimental absorption profiles.     

Phase behaviour and interfacial properties of ternary system CO2 + n-butane + n-decane: coarse-grained theoretical modelling and molecular simulations

ABSTRACT

This work illustrates the application of a three-party approach based on theoretical modelling, molecular dynamic (MD) simulations and available experimental data for describing the phase equilibrium and interfacial properties for the ternary system: carbon dioxide + n-butane + n-decane and its corresponding binary sub-systems at 344.3 K. Specifically, a coarse-grained force field is employed for both theoretical predictions and MD. The interfacial region is described by the square gradient theory where the homogenous Helmholtz energy density contribution is provided by the Statistical Associated Fluid Theory equation of state for potentials of variable range for molecules conformed of segments interacting through the Mie potential (SAFT-VR Mie) and MD simulations in the canonical ensemble where the molecules are represented by a coarse-grained Mie force field. The novelty here is that both the theory and the simulations uniquely share the same underlying intermolecular potentials; hence, the experimental data are employed to verify both the theory and simulations. In this schema, the ternary mixture is full predictive as its parameters are only based on pure fluids parameters and binary interactions. It is observed that the phase equilibria and the interfacial properties are equally well represented by the used approach.     

Test of the Gouy-Chapman theory for a charged lipid membrane against explicit-solvent molecular dynamics simulations

A wealth of experimental data has verified the applicability of the Gouy-Chapman (GC) theory to charged lipid membranes. Surprisingly, a validation of GC by molecular dynamics (MD) simulations has been elusive. Here, we report a test of GC against extensive MD simulations of an anionic lipid bilayer solvated by water at different concentrations of NaCl or KCl. We demonstrate that the ion distributions from the simulations agree remarkably well with GC predictions when information on the adsorption of counterions to the bilayer is incorporated.     

高压下液态硝基甲烷的分子动力学模拟

 对高压下液态硝基甲烷的性质进行经典和基于第一性原理计算的Car-Parrinello分子动力学（CPMD）模拟。利用经典势的分子动力学（MD）模拟研究了高压压缩状态下液态硝基甲烷的结构和热力学性质,得到了高达14.2 GPa压力下的理论Hugoniot数据。对于一些热力学函数,如总能和粒子速度,经典势模拟给出了很好的总趋势,基本特征和实验观测一致。但是在给定的密度下,经典模拟预言的Hugoniot压力偏高。在几个选定的密度下,进行了CPMD模拟,得到了二体相关函数、速度自相关函数、振动光谱和其它的热力学性质,并与经典模拟结果进行了比较。对二体相关函数的分析表明经典势的短程部分的刚性可能太强,从而导致了比实验值高的理论压力值。对于某些二体相关函数,CPMD模拟和经典模拟结果差别很大,可以归结为量子效应。当压力增高时,量子模拟得到的振动光谱向高频部分移动的现象与实验观测相符合。     

Comparative Study on Two Melting Simulation Methods: Melting Curve of Gold

Melting simulation methods are of crucial importance to determining melting temperature of materials efficiently.A high-efficiency melting simulation method saves much simulation time and computational resources.To compare the efficiency of our newly developed shock melting(SM)method with that of the well-established two-phase(TP)method,we calculate the high-pressure melting curve of Au using the two methods based on the optimally selected interatomic potentials.Although we only use 640 atoms to determine the melting temperature of Au in the SM method,the resulting melting curve accords very well with the results from the TP method using much more atoms.Thus,this shows that a much smaller system size in SM method can still achieve a fully converged melting curve compared with the TP method,implying the robustness and efficiency of the SM method.     

Molecular dynamics study of nanoparticle evolution in a background gas under laser ablation conditions

Long-time evolution of nanoparticles produced by short laser interactions is investigated for different materials. To better understand the mechanisms of the nanoparticle formation at a microscopic level, we use molecular dynamics (MD) simulations to analyse the evolution of a cluster in the presence of a background gas with different parameters (density and temperature). In particular, we compare the simulation results obtained for materials with different interaction potentials (Morse, Lennard-Jones, and Embedded Atom Model). Attention is focused on the evaporation and condensation processes of a cluster with different size and initial temperature. As a result of the MD calculations, we determinate the influence of both cluster properties and background gas parameters on the nanoparticle evolution. The role of the interaction potential is discussed based on the results of the simulations.     

Reactive force field potential for carbon deposition on silicon surfaces

In this paper a new interatomic potential based on the Kieffer force field and designed to perform molecular dynamics (MD) simulations of carbon deposition on silicon surfaces is implemented. This potential is a third-order reactive force field that includes a dynamic charge transfer and allows for the formation and breaking of bonds. The parameters for Si-C and C-C interactions are optimized using a genetic algorithm. The quality of the potential is tested on its ability to model silicon carbide and diamond physical properties as well as the formation energies of point defects. Furthermore, MD simulations of carbon deposition on reconstructed (100) silicon surfaces are carried out and compared to similar simulations using a Tersoff-like bond order potential. Simulations with both potentials produce similar results showing the ability to extend the use of the Kieffer potential to deposition studies. The investigation reveals the presence of a channelling effect when depositing the carbon at 45°?incidence angle. This effect is due to channels running in directions symmetrically equivalent to the (110) direction. The channelling is observed to a lesser extent for carbon atoms with 30°?and 60°?incidence angles relative to the surface normal. On a pristine silicon surface, sticking coefficients were found to vary between 100 and 73%, depending on deposition conditions.     

Exploration of Entropy Pair Functional Theory

Evaluation of the entropy from molecular dynamics (MD) simulation remains an outstanding challenge. The standard approach requires thermodynamic integration across a series of simulations. Recent work Nicholson et al. demonstrated the ability to construct a functional that returns excess entropy, based on the pair correlation function (PCF); it was capable of providing, with acceptable accuracy, the absolute excess entropy of iron simulated with a pair potential in both fluid and crystalline states. In this work, the general applicability of the Entropy Pair Functional Theory (EPFT) approach is explored by applying it to three many-body interaction potentials. These potentials are state of the art for large scale models for the three materials in this study: Fe modelled with a modified embedded atom method (MEAM) potential, Cu modelled with an MEAM and Si modelled with a Tersoff potential. We demonstrate the robust nature of EPFT in determining excess entropy for diverse systems with many-body interactions. These are steps toward a universal Entropy Pair Functional, EPF, that can be applied with confidence to determine the entropy associated with sophisticated optimized potentials and first principles simulations of liquids, crystals, engineered structures, and defects.     

Brittle and ductile fracture of semiconductor nanowires – molecular dynamics simulations

Fracture of silicon and germanium nanowires in tension at room temperature is studied by molecular dynamics simulations using several interatomic potential models. While some potentials predict brittle fracture initiated by crack nucleation from the surface, most potentials predict ductile fracture initiated by dislocation nucleation and slip. A simple parameter based on the ratio between the ideal tensile strength and the ideal shear strength is found to correlate very well with the observed brittle versus ductile behaviours for all the potentials used in this study. This parameter is then computed by ab initio methods, which predict brittle fracture at room temperature. A brittle-to-ductile transition (BDT) is observed in MD simulations at higher temperature. The BDT mechanism in semiconductor nanowires is different from that in the bulk, due to the lack of a pre-existing macrocrack that is always assumed in bulk BDT models.     

Molecular dynamic studies on materials under laser shocks

Molecular dynamic (MD) simulations offer a powerful means of understanding the microscopic characteristics of shock-propagation through solids and fluids, especially for the short spatial and temporal scales relevant to laser-driven shocks. First-principles molecular dynamics can be directly compared with time-resolved experimental measurements, and methods based on empirical (embedded-atom) potentials fitted to first-principles quantum-mechanical calculations are effective for MD simulations of shock propagation through many millions of atoms. In comparison, thermodynamic approaches based on free-energy considerations do not provide detailed information about mechanical-relaxation or phase-transformation processes within the shock front. We illustrate these ideas by way of embedded-atom simulations of shock-wave propagation through copper crystals of different orientation.     

Predicting XAFS scattering path cumulants and XAFS spectra for metals (Cu,Ni, Fe,Ti, Au) using molecular dynamics simulations

The ability of molecular dynamics (MD) simulations to support the analysis of X‐ray absorption fine‐structure (XAFS) data for metals is evaluated. The low‐order cumulants (ΔR, σ², C₃) for XAFS scattering paths are calculated for the metals Cu, Ni, Fe, Ti and Au at 300 K using 28 interatomic potentials of the embedded‐atom method type. The MD cumulant predictions were evaluated within a cumulant expansion XAFS fitting model, using global (path‐independent) scaling factors. Direct simulations of the corresponding XAFS spectra, χ(R), are also performed using MD configurational data in combination with the FEFFab initio code. The cumulant scaling parameters compensate for differences between the real and effective scattering path distributions, and for any errors that might exist in the MD predictions and in the experimental data. The fitted value of ΔR is susceptible to experimental errors and inadvertent lattice thermal expansion in the simulation crystallites. The unadjusted predictions of σ² vary in accuracy, but do not show a consistent bias for any metal except Au, for which all potentials overestimate σ². The unadjusted C₃ predictions produced by different potentials display only order‐of‐magnitude consistency. The accuracy of direct simulations of χ(R) for a given metal varies among the different potentials. For each of the metals Cu, Ni, Fe and Ti, one or more of the tested potentials was found to provide a reasonable simulation of χ(R). However, none of the potentials tested for Au was sufficiently accurate for this purpose.     

Thermal and mechanical properties of thermosetting polymers using coarse-grained simulation

We developed coarse-grained (CG) molecular representations of mixtures of diglycidyl ether of bisphenol-A (DGEBA) and poly(oxypropylene) diamine (POP-DA) for use in CG molecular dynamics (MD) simulations. In the CG representation, DGEBA is comprised of three beads of two types and POP-DA also by three beads of two types. Atomistic MD of liquid systems was performed to derive intra- and inter-bead potentials via Boltzmann inversion. While the bonded potentials, composed of bond stretching and angle bending, were parameterized directly from the distribution functions of all atomistic molecular dynamics trajectories, the non-bonded potentials were derived from the iterative Boltzmann Inversion with a given set of coarse-grained interactions. CG systems correctly reproduced liquid and crosslinked densities. Under uniaxial tension, the Young's modulus of the CG systems was much lower than the experimental value, and we show this arises from the assumed form of the extrapolated regions of the CG potentials. By stiffening these regions, we increased the CG Young's modulus of the crosslinked systems without sacrificing the correct prediction of density. This suggests that transferrable CG potentials can be optimized for use in non-equilibrium MD for property estimation.