Solid–liquid coexistence in ionic systems and the properties of the interface

The effect of dispersion interactions on the properties of LiF in condensed phases

Classical molecular dynamics simulations are performed on LiF in the framework of the polarizable ion model. The overlap repulsion and polarization terms of the interaction potential are derived on a purely non-empirical, first-principles basis. For the dispersion, three cases are considered: a first one in which the dispersion parameters are set to zero and two others in which they are included, with different parametrizations. Various thermodynamic, structural and dynamic properties are calculated for the solid and liquid phases. The melting temperature is also obtained from direct coexistence simulations of the liquid and solid phases. Dispersion interactions appear to have an important effect on the densities of both phases and on the melting point, although the liquid properties are not affected when simulations are performed in the NVT ensemble at the experimental density.

Electrochemical interface between an ionic liquid and a model metallic electrode

A molecular dynamics simulation model for an electroactive interface in which a metallic electrode is maintained at a preset electrical potential is described. The model, based on earlier work of Siepmann and Sprik [J. Chem. Phys. 102, 511 (1995)], uses variable charges whose magnitudes are adjusted on-the-fly according to a variational procedure to maintain the constant potential condition. As such, the model also allows for the polarization of the electrode by the electrolyte, sometimes described by the introduction of image charges. The model has been implemented in a description of an electrochemical cell as a pair of parallel planar electrodes separated by the electrolyte using a two-dimensional Ewald summation method. The method has been applied to examine the interfacial structure in two ionic liquids, consisting of binary mixtures of molten salts, chosen to exemplify the influences of dissimilar cation size and charge. The stronger coordination of the smaller and more highly charged cations by the anions prevents them from approaching even the negatively charged electrode closely. This has consequences for the capacitance of the electrode and will also have an impact on the rates of electron transfer processes. The calculated capacitances exhibit qualitatively the same dependence on the applied potential as has been observed in experimental studies.

Phase diagrams of alkali halides using two interaction models: A molecular dynamics and free energy study

Phase diagrams for potassium and sodium chlorides are determined by molecular dynamics and free energy calculations. Two rigid-ion interaction models, namely, the Born-Mayer-Huggins (BMH) and Michielsen-Woerlee-Graaf (MWG) effective pair potentials, have been used. The critical and triple point properties are discussed and compared with available experimental and simulation data. The MWG model reproduces the experimental liquid-gas equilibria better than the BMH model, being the accordance very good in the lowest temperature region of the coexistent liquids, particularly for NaCl. However, both models underestimate the critical temperatures of KCl and NaCl. Relatively to the solid-gas equilibria, the models do not reproduce well the experimental data. As for the solid-liquid coexistences either the BMH or the MWG models appear unrealistic.

Mechanisms of heterogeneous crystal growth in atomic systems: Insights from computer simulations

In this paper we analyze the atomic-level structure of solid/liquid interfaces of Lennard-Jones fcc systems. The 001, 011, and 111 faces are examined during steady-state growth and melting of these crystals. The mechanisms of crystallization and melting are explored using averaged configurations generated during these steady-state runs, where subsequent tagging and labeling of particles at the interface provide many insights into the detailed atomic behavior at the freezing and melting interfaces. The interfaces are generally found to be rough and we observe the structure of freezing and melting interfaces to be very similar. Large structural fluctuations with solidlike and liquidlike characteristics are apparent in both the freezing and melting interfaces. The behavior at the interface observed under either growth or melting conditions reflects a competition between ordering and disordering processes. In addition, we observe atom hopping that imparts liquidlike characteristics to the solid side of the interfaces for all three crystal faces. Solid order is observed to extend as rough, three-dimensional protuberances through the interface, particularly for the 001 and 011 faces. We are also able to reconcile our different measures for the interfacial width and address the onset of asymmetry in the growth rates at high rates of crystal growth/melting.