Understanding solid nitrogen through molecular dynamics simulations with a machine-learning potential

We construct a fast, transferable, general purpose, machine-learning interatomic potential suitable for large-scale simulations of <a:math xmlns:a="http://www.w3.org/1998/Math/MathML"><a:msub><a:mi mathvariant="normal">N</a:mi><a:mn>2</a:mn></a:msub></a:math>. The potential is trained only on high quality quantum chemical molecule-molecule interactions; no condensed phase information is used. Although there are no explicit or implicit many-molecule interaction terms, the potential reproduces the experimental phase diagram including the melt curve and the molecular solid phases of nitrogen up to <c:math xmlns:c="http://www.w3.org/1998/Math/MathML"><c:mrow><c:mn>10</c:mn><c:mspace width="4pt"/><c:mi>GPa</c:mi></c:mrow></c:math>. This demonstrates that many-molecule interactions are unnecessary to explain the condensed phases of <e:math xmlns:e="http://www.w3.org/1998/Math/MathML"><e:msub><e:mi mathvariant="normal">N</e:mi><e:mn>2</e:mn></e:msub></e:math>. With increased pressure, transitions are observed from cubic <g:math xmlns:g="http://www.w3.org/1998/Math/MathML"><g:mo>(</g:mo><g:mi>α</g:mi><g:mo>)</g:mo></g:math>, which optimizes quadrupole-quadrupole interactions, through tetragonal <h:math xmlns:h="http://www.w3.org/1998/Math/MathML"><h:mo>(</h:mo><h:mi>γ</h:mi><h:mo>)</h:mo></h:math>, which allows more efficient packing, to monoclinic <i:math xmlns:i="http://www.w3.org/1998/Math/MathML"><i:mo>(</i:mo><i:mi>λ</i:mi><i:mo>)</i:mo></i:math>, which packs still more efficiently. On heating, we obtain the hcp three-dimensional (3D) rotor phase <j:math xmlns:j="http://www.w3.org/1998/Math/MathML"><j:mo>(</j:mo><j:mi>β</j:mi><j:mo>)</j:mo></j:math> and, at pressure, the cubic <k:math xmlns:k="http://www.w3.org/1998/Math/MathML"><k:mi>δ</k:mi></k:math> phase which contains both 3D and 2D rotors, tetragonal <l:math xmlns:l="http://www.w3.org/1998/Math/MathML"><l:msup><l:mi>δ</l:mi><l:mo>*</l:mo></l:msup></l:math> phase with 2D rotors, and the rhombohedral <m:math xmlns:m="http://www.w3.org/1998/Math/MathML"><m:mi>ε</m:mi></m:math>. Molecular dynamics demonstrates where these phases are indeed rotors, rather than frustrated order. The model supports the metastability of the complex <n:math xmlns:n="http://www.w3.org/1998/Math/MathML"><n:mi>ι</n:mi></n:math> phase, but not the reported existence of the wide range of bond lengths. The thermodynamic transitions involve both shifts of molecular centers and rotations of molecules: the onset of rotation is rapid, whereas motion of molecular centers is inhibited and we suggest that this is the cause of the experimentally observed sluggishness of transitions. Routine density functional theory calculations give a similar picture to the potential. Published by the American Physical Society 2024