Search

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. 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., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

Tin sulfides, Sn x S y , are an important class of materials that are actively investigated as novel photovoltaic and water splitting materials. A first-principles evolutionary crystal structure search is performed with the goal of constructing the complete phase diagram of Sn x S y and discovering new phases as well as new compounds of varying stoichiometry at ambient conditions and pressures up to 100 GPa. The ambient phase of SnS 2 with P 3 ¯ m 1 symmetry remains stable up to 28 GPa. Another ambient phase, SnS, experiences a series of phase transformations including α -SnS to β -SnS at 9 GPa, followed by β -SnS to γ -SnS at 40 GPa. γ -SnS is a new high-pressure metallic phase with P m 3 ¯ m space group symmetry stable up to 100 GPa, which becomes a superconductor with a maximum T c = 9.74 K at 40 GPa. Another new metallic compound, Sn 3 S 4 with I 4 ¯ 3 d space group symmetry, is predicted to be stable at pressures above 15 GPa, which also becomes a superconductor with relatively high T c = 21.9 K at 30 GPa.