Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon.

Diamond possesses exceptional physical properties due to its remarkably strong carbon-carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure-temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. We propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

Water in an External Electric Field: Comparing Charge Distribution Methods Using ReaxFF Simulations.

The growing interest in the effects of external electric fields on reactive processes requires predictive methods that can reach longer length and time scales than quantum mechanical simulations. Recently, many studies have included electric fields in ReaxFF, a widely used reactive molecular dynamics method. In the case of modeling an external electric field, the charge distribution method used in ReaxFF is critical. The most common charge distribution method used in previous studies of electric fields is the charge equilibration (QEq) method, which assumes that the system is a contiguous conductor and that charge transfer can occur across any distance. In contrast, many systems of interest are insulators or semiconductors, and long-distance charge transfer should not occur in response to a small difference in potential. This study focuses on the limitations of the QEq method in the context of water in an external electric field. We demonstrate that QEq can predict unphysical charge distributions and exhibits properties that do not converge as a function of system size. Furthermore, we show that electric fields within the recently developed atom-condensed Kohn-Sham density functional theory (DFT) approximated to the second-order (ACKS2) approach address the major limitations of electric fields in QEq. With ACKS2, we observe more physical charge distributions and properties that converge as a function of system size. We do not suggest that ACKS2 is perfect in all circumstances but rather show specific cases where it addresses the major shortcomings of QEq in the context of an external electric field.

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields.

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. In this Perspective, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

Effect of an external field on capillary waves in a dipolar fluid.

The role of an external field on capillary waves at the liquid-vapor interface of a dipolar fluid is investigated using molecular dynamics simulations. For fields parallel to the interface, the interfacial width squared increases linearly with respect to the logarithm of the size of the interface across all field strengths tested. The value of the slope decreases with increasing field strength, indicating that the field dampens the capillary waves. With the inclusion of the parallel field, the surface stiffness increases with increasing field strength faster than the surface tension. For fields perpendicular to the interface, the interfacial width squared is linear with respect to the logarithm of the size of the interface for small field strengths, and the surface stiffness is less than the surface tension. Above a critical field strength that decreases as the size of the interface increases, the interface becomes unstable due to the increased amplitude of the capillary waves.

Liquid-vapor interface of the Stockmayer fluid in a uniform external field.

The effect of a uniform (nonspatially varying) external field on the liquid-vapor interface of the Stockmayer fluid (Lennard-Jones particles embedded with a point dipole) has been investigated by molecular-dynamics simulations. The long-ranged parts of both the dipole and Lennard-Jones interactions are treated using an Ewald summation, which removes the effects of the cutoff. The direction of the field shifts the critical point and interfacial properties in different directions. For an external field parallel to the interface, the critical temperature increases, while for a field applied perpendicular to the interface, it decreases. The effects of the field on surface tension and interfacial width are also investigated. For zero field, dipoles near the liquid-vapor interface show a weak orientation parallel to the interface. For fields parallel to the interface, ordering in the liquid phase is greater than the vapor, while for fields perpendicular to the interface, the opposite is true.

Extension and evaluation of the multilevel summation method for fast long-range electrostatics calculations.

Several extensions and improvements have been made to the multilevel summation method (MSM) of computing long-range electrostatic interactions. These include pressure calculation, an improved error estimator, faster direct part calculation, extension to non-orthogonal (triclinic) systems, and parallelization using the domain decomposition method. MSM also allows fully non-periodic long-range electrostatics calculations which are not possible using traditional Ewald-based methods. In spite of these significant improvements to the MSM algorithm, the particle-particle particle-mesh (PPPM) method was still found to be faster for the periodic systems we tested on a single processor. However, the fast Fourier transforms (FFTs) that PPPM relies on represent a major scaling bottleneck for the method when running on many cores (because the many-to-many communication pattern of the FFT becomes expensive) and MSM scales better than PPPM when using a large core count for two test problems on Sandia's Redsky machine. This FFT bottleneck can be reduced by running PPPM on only a subset of the total processors. MSM is most competitive for relatively low accuracy calculations. On Sandia's Chama machine, however, PPPM is found to scale better than MSM for all core counts that we tested. These results suggest that PPPM is usually more efficient than MSM for typical problems running on current high performance computers. However, further improvements to MSM algorithm could increase its competitiveness for calculation of long-range electrostatic interactions.

Characteristics of thermal conductivity in classical water models.

The thermal conductivities of common water models are compared using equilibrium (EMD) and non-equilibrium molecular dynamics (NEMD) simulation. A complete accounting for electrostatic contributions to the heat flux was found to resolve the previously reported differing results of NEMD and EMD Green-Kubo measurements for the extended simple point-charge (SPC/E) model. Accordingly, we demonstrate the influence of long-range electrostatics on the thermal conductivity with a simple coulomb cutoff, Ewald summation, and by an extended particle-particle particle-mesh method. For each water model, the thermal conductivity is computed and decomposed in terms of frequency-dependent thermodynamic and topological contributions. The rigid, three-site SPC, SPC/E, and transferable intermolecular potential (TIP3P-Ew) water models are shown to have similar thermal conductivity values at standard conditions, whereas models that include bond stretching and angle bending have higher thermal conductivities.

Chemical potential perturbation: extension of the method to lattice sum treatment of intermolecular potentials.

The recently developed chemical potential perturbation (CPP) method [S. G. Moore and D. R. Wheeler, J. Chem. Phys. 134, 114514 (2011)] is extended to the lattice (Ewald) sum treatment of intermolecular potentials. The CPP method predicts chemical potentials for a range of composition points using the local (position-dependent) pressure tensor of an inhomogeneous system. When computing the local pressure tensor, one can use the Irving-Kirkwood (IK) or Harasima (H) contours of distributing the pressure. We compare these two contours and show that for a planar interface, the homogeneous pressure and resulting chemical potential can be approximated with the CPP method using either the IK or the H contour, though with the lattice sum method the H contour has much greater computational efficiency. The proposed methods are validated by calculating the chemical potentials of the Lennard-Jones fluid and extended simple point-charge (SPC/E) water, and results show a high level of agreement with respective equations of state.

Chemical potential perturbation: a method to predict chemical potentials in periodic molecular simulations.

A new method, called chemical potential perturbation (CPP), has been developed to predict the chemical potential as a function of density in periodic molecular simulations. The CPP method applies a spatially varying external force field to the simulation, causing the density to depend upon position in the simulation cell. Following equilibration the homogeneous (uniform or bulk) chemical potential as a function of density can be determined relative to some reference state after correcting for the effects of the inhomogeneity of the system. We compare three different methods of approximating this correction. The first method uses the van der Waals density gradient theory to approximate the inhomogeneous Helmholtz free energy density. The second method uses the local pressure tensor to approximate the homogeneous pressure. The third method uses the Triezenberg-Zwanzig definition of surface tension to approximate the inhomogeneous free energy density. If desired, the homogeneous pressure and Helmholtz free energy can also be predicted by the new method, as well as binodal and spinodal densities of a two-phase fluid region. The CPP method is tested using a Lennard-Jones (LJ) fluid at vapor, liquid, two-phase, and supercritical conditions. Satisfactory agreement is found between the CPP method and an LJ equation of state. The efficiency of the CPP method is compared to that for Widom's method under the tested conditions. In particular, the new method works well for dense fluids where Widom's method starts to fail.

Improved implementation of Kirkwood-Buff solution theory in periodic molecular simulations.

Kirkwood-Buff (KB) solution theory is a means to obtain certain thermodynamic derivatives from knowledge of molecular distributions. In actual practice the required integrals over radial distribution functions suffer inaccuracies due to finite-distance truncation effects and their use in closed systems. In this work we discuss how best to minimize these inaccuracies under traditional KB theory. In addition we implement a method for calculating KB quantities in molecular simulations with periodic boundary conditions and particularly within the canonical ensemble. The method is based on a finite-Fourier-series expansion of molecular concentration fluctuations and leads to more reliable results for a given computational effort. The procedure is validated and compared to the original method for a nonideal liquid mixture of Lennard-Jones particles intended to imitate a real system, carbon tetrafluoride, and methane.