Abstract:We describe the first principle‐based electron force field (eFF) methodology for modeling the simultaneous dynamics of electrons and nuclei (eMD) evolving nonadiabatically under transient extreme conditions. We introduce the parallel implementation of eFF (pEFF) that makes it practical to perform simulations of the nonadiabatic dynamics of materials in extreme environments involving millions of nuclei and electrons, over multi‐picoseconds time scales, and demonstrate its application to: (i) accurately determine density and predict percent ionization of hydrogen at high pressure (∼61 GPa) and temperatures up to 15,300 K and (ii) determine, the single shock Hugoniot for lithium metal directly from the shock wave kinematics, i.e., mass velocities (Up) and shock wave velocities (Us), and shock density data. For (i), the density and ionization fractions of hydrogen atoms were calculated using the isobaric–isothermal ensemble at an isotropic pressure of 61.4 GPa and for temperatures between 300 K and 15,300 K. The results at 15,300 K describe a molecular fluid with density ρ = 0.36 g/cm3, in close agreement with existing experiments and theory, and ∼0.5% ionization. This result provides no indication of the existence of a critical plasma phase‐transition point at this particular temperature and pressure, as previously predicted by others. For (ii), the relationship between Up and Us was characterized to be linear and plastic in the range 1–20 km/s, and the single shock Hugoniot was determined in close agreement with published results for experimentally reported Ups. In addition to this, we provide a description of the materials' behavior for large Ups in terms of the appearance of a weak metallic plasma phase by Up = 10 km/s, with ≃ 8% ionization, gradually transitioning to a denser plasma with an estimated ≃ 35% ionization by Up = 15 km/s. Last but not least, we confirm the computational efficiency and scalability of pEFF by comparing its single processor performance against the fastest existing serial code, which results in a linear speedup ∼10× faster for every 16,000 particles in favor of pEFF, and by evaluating its parallel performance in terms of its strong and weak scaling capabilities. Our results, on Los Alamos's Lobo supercomputer (a 38TFLOPSs Linux HPC with Quad‐core AMD Opteron nodes interconnected with an Infiniband), show strong scaling with near ideal speedups for loads >62 particles per processor. Weak scaling is shown to be close to linear under the same per‐processor load range. As an absolute reference, an NVT run with 2 million particle lithium bulk system (0.5 M nuclei and 1.5 M electrons) on Lobo takes ∼0.44 s/timestep on 1024 processors (∼1 day/ps using an integration timestep of 0.005 fs). © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011

High-Performance Scalable Molecular Dynamics Simulations of a Polarizable Force Field Based on Classical Drude Oscillators in NAMD

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Molecular dynamics simulations using the drude polarizable force field on GPUs with OpenMM: Implementation, validation, and benchmarks

Balancing Group I Monoatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Super-Scalable Molecular Dynamics Algorithm

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

Scalable molecular dynamics on CPU and GPU architectures with NAMD

Convergence of stochastic-extended Lagrangian molecular dynamics method for polarizable force field simulation

Large-scale, Long-Term Nonadiabatic Electron Molecular Dynamics for Describing Material Properties and Phenomena in Extreme Environments.

Incorporating surface polarization effects into large-scale coarse-grained Molecular Dynamics simulation

Interfacing the Core-Shell or the Drude Polarizable Force Field With Car-Parrinello Molecular Dynamics for QM/MM Simulations

Advanced Potential Energy Surfaces for Molecular Simulation

A fast, convenient, polarizable electrostatic model for molecular dynamics

Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Electronic coarse graining enhances the predictive power of molecular simulation allowing challenges in water physics to be addressed

Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field.

CHARMM‐GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Accurate interatomic force field for molecular dynamics simulation by hybridizing classical and machine learning potentials

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques