Engineering topics

Gonzalez, Joseph M.

Publications and source records attributed to Gonzalez, Joseph M..

3 records

DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗

DOE OSTI2024

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. Here we propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Crystal structure of silver pentazolates AgN 5 and AgN 6

Silver pentazolate, a high energy density compound containing the cyclo-N 5 – anion, has recently been synthesized under ambient conditions. However, due to high sensitivity to irradiation, its crystal structure has not been determined. In this work, silver–nitrogen crystalline compounds under ambient conditions and at high pressures, up to 100 GPa, are predicted and characterized by performing first-principles evolutionary crystal structure searching with variable stoichiometry. It is found that newly discovered AgN 5 and AgN 6 are the only thermodynamically stable silver–nitrogen compounds at pressures between 42 and 80 GPa. In contrast to AgN 5 , the pentazolate AgN 6 compound contains N 2 diatomic molecules in addition to cyclo-N 5 – . These AgN 5 and AgN 6 crystals are metastable under ambient conditions with positive formation enthalpies of 54.95 kJ mol –1 and 46.24 kJ mol –1 , respectively. Here, the underlying cause of the stability of cyclo-N 5 – silver pentazolates is the enhanced aromaticity enabled by the charge transfer from silver atoms to nitrogen rings. To aid in the experimental identification of these materials, calculated Raman spectra are reported at ambient pressure: the frequencies of N 5 – vibrational modes of AgN 5 are in good agreement with those measured in the experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗