Calculation of the melting point of alkali halides by means of computer simulations

Compositional transferability of deep learning potentials: a case study for LiCl-KCl melt

CONTEXT

One of the crucial issues related to machine learning potentials is the formation of representative dataset. For multicomponent systems, it is a general methodology to scan the composition range with a certain step. However, there is a lack of information on the compositional transferability of machine learning potentials. In this paper, we extend the knowledge in this area by studying the transferability of deep learning potential over the range of compositions of LiCl-KCl molten mixtures. The training dataset was formed using only the near-eutectic composition of 60% LiCl-40% KCl. Then, we tested the ability of the model to predict physicochemical properties of the melts far from the reference composition. It was found that for the composition range from 0 to 100% of LiCl, the calculated properties concur closely with those of other studies and ab initio calculations. Therefore, the model shows prominent non-intuitive compositional transferability. Moreover, the solid states and solid-liquid coexistence were reproduced. The calculated melting temperatures of LiCl and KCl show the errors of 6.6% and 0.4%, respectively. We argue that such good transferability stems from the local structure configurations that are typical both for pure LiCl and for pure KCl which are implicitly presented in the training dataset because of local fluctuations in composition.

METHODS

To collect the data for the initial dataset, density functional theory was employed. Then, the DeePMD package was used to train a neural network potential. To calculate the properties of the melts, standard equilibrium and non-equilibrium molecular dynamic approaches were utilized.

Estimation of thermodynamic and physicochemical properties of the alkali astatides: On the bond strength of molecular astatine (At2 ) and the hydration enthalpy of astatide (At- )

The recent accurate and precise determination of the electron affinity (EA) of the astatine atom At0 warrants a re-investigation of the estimated thermodynamic properties of At0 and astatine containing molecules as this EA was found to be much lower (by 0.4 eV) than previous estimated values. In this contribution we estimate, from available data sources, the following thermodynamic and physicochemical properties of the alkali astatides (MAt, M = Li, Na, K, Rb, Cs): their solid and gaseous heats of formation, lattice and gas-phase binding enthalpies, sublimation energies and melting temperatures. Gas-phase charge-transfer dissociation energies for the alkali astatides (the energy requirement for M+ At- ➔ M0  + At0 ) have been obtained and are compared with those for the other alkali halides. Use of Born-Haber cycles together with the new AE (At0 ) value allows the re-evaluation of ΔHf (At0 )g (=56 ± 5 kJ/mol); it is concluded that (At2 )g is a weakly bonded species (bond strength <50 kJ/mol), significantly weaker bonded than previously estimated (116 kJ/mol) and much weaker bonded than I2 (148 kJ/mol), but in agreement with the finding from theory that spin-orbit coupling considerably reduces the bond strength in At2 . The hydration enthalpy (ΔHaq ) of At- is estimated to be -230 ± 2 kJ/mol (using ΔHaq [H+ ] = -1150.1 kJ/mol), in good agreement with molecular dynamics calculations. Arguments are presented that the largest alkali halide, CsAt, like the smallest, LiF, will be only sparingly soluble in water, following the generalization from hard/soft acid/base principles that "small likes small" and "large likes large."

First-principles molten salt phase diagrams through thermodynamic integration

Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. We show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid-liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange-correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid-liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

Origin of low melting point of ionic liquids: dominant role of entropy

Large structural entropy makes salts liquid at room temperature.

Study of the structure of a multicomponent salt melt using molecular dynamics modeling

The composition of the electrolyte is critical in the electrodeposition of high-purity silicon. In this work, molecular dynamics modeling of the preparation of liquid salt melt KF-KCl-KI and a detailed study of its structure based on the method of statistical geometry have been performed. Partial radial distribution functions reflect the size of the ions under consideration and the averaged structure of the generated ionic subsystems. Halogen subsystems have domed angular distributions of nearest geometric neighbors, a wide range of face types of combined polyhedra, and fifth order rotational symmetry. The shape of the distribution of distances to the nearest neighbors of a given type depends on the amount of these ions in the melt. Small-scale thermal fluctuations in the halogen subsystems are predominantly represented by small triangular faces in combined polyhedra. The electrodeposition of silicon was carried out in a homogeneous salt melt, in which each halogen ion had from one to three close contacts with halogen ions of any other type. The simulations performed provide a fundamental understanding of the structure of the electrolyte molten salts used to produce solar silicon.

Effective mass path integral simulations of quasiparticles in condensed phases

The quantum many-body problem in condensed phases is often simplified using a quasiparticle description, such as effective mass theory for electron motion in a periodic solid. These approaches are often the basis for understanding many fundamental condensed phase processes, including the molecular mechanisms underlying solar energy harvesting and photocatalysis. Despite the importance of these effective particles, there is still a need for computational methods that can explore their behavior on chemically relevant length and time scales. This is especially true when the interactions between the particles and their environment are important. We introduce an approach for studying quasiparticles in condensed phases by combining effective mass theory with the path integral treatment of quantum particles. This framework incorporates the generally anisotropic electronic band structure of materials into path integral simulation schemes to enable modeling of quasiparticles in quantum confinement, for example. We demonstrate the utility of effective mass path integral simulations by modeling an exciton in solid potassium chloride and electron trapping by a sulfur vacancy in monolayer molybdenum disulfide.

Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems

Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.

Computational Assessment of the Molecular Structure and Properties for High Energy Density Fuel

High energy density fuels (HEDFs) that have high volumetric net heat of combustion (NHOC), high stability, and high environmental resistance are greatly important in the fuel field and in military bases and aerospace applications. In this paper, molecular dynamics and quantum chemistry were used to compute the significant physical properties of candidate molecules for HEDFs, such as their enthalpies of combustion, enthalpies of vaporization, densities, and melting points. A computational protocol for evaluating these properties in the fuel field was established, including a new method for estimating the melting point. By using our protocol, we found that, to improve fuel performances such as the density and volumetric NHOC, cyclopropanation is obviously better than hydrogenation. Our protocol was verified to have good accuracy and can be used to compare and assess different target molecules as potential HEDFs.

Systematically improved melting point prediction: a detailed physical simulation model is required

Our detailed physical simulation model allows for an unprecedented and systematically improved prediction of melting points of alkali halides.

Calculation of the water-octanol partition coefficient of cholesterol for SPC, TIP3P, and TIP4P water

We present a numerical study of the relative solubility of cholesterol in octanol and water. Our calculations allow us to compare the accuracy of the computed values of the excess chemical potential of cholesterol for several widely used water models (SPC, TIP3P, and TIP4P). We compute the excess solvation free energies by means of a cavity-based method [L. Li et al., J. Chem. Phys. 146(21), 214110 (2017)] which allows for the calculation of the excess chemical potential of a large molecule in a dense solvent phase. For the calculation of the relative solubility ("partition coefficient," log₁₀ P o / w ) of cholesterol between octanol and water, we use the OPLS/AA force field in combination with the SPC, TIP3P, and TIP4P water models. For all water models studied, our results reproduce the experimental observation that cholesterol is less soluble in water than in octanol. While the experimental value for the partition coefficient is log₁₀ P o / w = 3.7, SPC, TIP3P, and TIP4P give us a value of log₁₀ P o / w = 4.5, 4.6, and 2.9, respectively. Therefore, although the results for the studied water models in combination with the OPLS/AA force field are acceptable, further work to improve the accuracy of current force fields is needed.

Molecular dynamics determination of liquid-vapor coexistence in molten alkali halides

We show by extensive molecular dynamics simulations that accurate predictions of liquid-vapor coexistence in molten alkali halides can be achieved in terms of a rigid ion potential description in which temperature-dependent ionic diameters are employed. The new ionic sizes result from the fitting of the experimental isothermal compressibilities, a condition whose physical implications and consequences are illustrated. The same diameters also allow us to formulate confident predictions for the compressibilities of salts in cases where the experimental data are lacking. The extension of the present approach to molten alkali-halide mixtures and to other classes of molten salts is discussed.

Single-pass incremental force updates for adaptively restrained molecular dynamics

Adaptively restrained molecular dynamics (ARMD) allows users to perform more integration steps in wall-clock time by switching on and off positional degrees of freedoms. This article presents new, single-pass incremental force updates algorithms to efficiently simulate a system using ARMD. We assessed different algorithms for speedup measurements and implemented them in the LAMMPS MD package. We validated the single-pass incremental force update algorithm on four different benchmarks using diverse pair potentials. The proposed algorithm allows us to perform simulation of a system faster than traditional MD in both NVE and NVT ensembles. Moreover, ARMD using the new single-pass algorithm speeds up the convergence of observables in wall-clock time. © 2017 Wiley Periodicals, Inc.

Crystal structures of model lithium halides in bulk phase and in clusters

We employ lattice energy calculations and molecular dynamics simulations to compare the stability of wurtzite and rock salt crystal structures of four lithium halides (LiF, LiCl, LiBr, and LiI) modeled using the Tosi-Fumi and Joung-Cheatham potentials, which are models frequently used in simulation studies. Both infinite crystals and finite clusters are considered. For the Tosi-Fumi model, we find that all four salts prefer the wurtzite structure both at 0 K and at finite temperatures, in disagreement with experiments, where rock salt is the stable structure and wurtzite exists as a metastable state. For Joung-Cheatham potentials, rock salt is more stable for LiF and LiCl, but the wurtzite structure is preferred by LiBr and LiI. It is clear that the available lithium halide force fields need improvement to bring them into better accord with the experiment. Finite-size clusters that are more stable as rock salt in the bulk phase tend to solidify as small rock salt crystals. However, small clusters of salts that prefer the wurtzite structure as bulk crystals tend to form structures that have hexagonal motifs, but are not finite-size wurtzite crystals. We show that small wurtzite structures are unstable due to the presence of a dipole and rearrange into more stable, size-dependent structures. We also show that entropic contributions can act in favor of the wurtzite structure at higher temperatures. The possible relevance of our results for simulation studies of crystal nucleation from melts and/or aqueous solutions is discussed.

Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route

The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.

Seeding approach to crystal nucleation

We present a study of homogeneous crystal nucleation from metastable fluids via the seeding technique for four different systems: mW water, Tosi-Fumi NaCl, Lennard-Jones, and Hard Spheres. Combining simulations of spherical crystal seeds embedded in the metastable fluid with classical nucleation theory, we are able to successfully describe the nucleation rate for all systems in a wide range of metastability. The crystal-fluid interfacial free energy extrapolated to coexistence conditions is also in good agreement with direct calculations of such parameter. Our results show that seeding is a powerful technique to investigate crystal nucleation.

Ligand-Mediated Interactions between Nanoscale Surfaces Depend Sensitively and Nonlinearly on Temperature, Facet Dimensions, and Ligand Coverage

Nanoparticles are often covered in ligand monolayers, which can undergo a temperature-dependent order-disorder transition that switches the particle-particle interaction from repulsive to attractive in solution. In this work, we examine how changes in the ligand surface coverage and facet dimensions affect the ordering of ligands, the arrangement of nearby solvent molecules, and the interaction between ligand monolayers on different particles. In particular, we consider the case of strongly bound octadecyl ligands on the (100) facet of CdS in the presence of an explicit n-hexane solvent. Depending on the facet dimensions and surface coverage, we observe three distinct ordered states that differ in how the ligands are packed together, and which affect the thickness of the ligand shell and the structure of the ligand-solvent interface. The temperature dependence of the order-disorder transition also broadens and shifts to lower temperature in a nonlinear manner as the nanoscale is approached from above. We find that ligands on nanoscale facets can behave very similarly to those on macroscopic surfaces in solution, and that some facet dimensions affect the ligand alignment more strongly than others. As the ligands order, the interaction between opposing monolayers becomes attractive, even well below full surface coverage. The strength of attraction per unit surface area is strongly affected by ligand coverage, but only weakly by facet width. Conversely, we find that bringing two monolayers together just above the order-disorder transition temperature can induce ordering and attraction.

Lattice mold technique for the calculation of crystal nucleation rates

We present a new simulation method for the calculation of crystal nucleation rates by computer simulation. The method is based on the use of molds to induce crystallization in state points where nucleation is a rare event. The mold is a cluster of potential energy wells placed in the lattice positions of the solid. The method has two distinct steps. In the first one the probability per unit volume of forming a sub-critical crystal cluster in the fluid is computed by means of thermodynamic integration. The thermodynamic route consists in gradually switching on an attractive interaction between the wells and the fluid particles. In the second step, the frequency with which such cluster becomes post-critical is computed in Molecular Dynamics simulations with the mold switched on. We validate our method with a continuous version of the hard sphere potential and with the sodium chloride Tosi–Fumi model. In all studied state points we obtain a good agreement with literature data obtained from other rare event simulation techniques. Our method is quite suitable for the study of both crystal nucleation of arbitrarily complex structures and the competition between different polymorphs in the nucleation stage.