Stable all-nitrogen metallic salt at terapascal pressures

One-Dimensional Non-coplanar Nitrogen Chains in Manganese Tetranitride under High Pressure

Transition metal nitrides have great potential applications as incompressible and high energy density materials. Various polymeric nitrogen structures significantly affect their properties, contributing to their complex bonding modes and coordination conditions. Herein, we first report a new manganese polynitride MnN4 with bifacial trans-cis [N4]n chains by treating with high-pressure and high-temperature conditions in a diamond anvil cell. Our experiments reveal that MnN4 has a P-1 symmetry and could stabilize in the pressure range of 56-127 GPa. Detailed pressure-volume data and calculations of this phase indicate that MnN4 is a potential hard (255 GPa) and high energy density (2.97 kJ/g) material. The asymmetric interactions impel N1 and N4 atoms to hybridize to sp2-3, which causes distortions of [N4]n chains. This work discovers a new polynitride material, fills the gap for the study of manganese polynitride under high pressure, and offers some new insights into the formation of polymeric nitrogen structures.

Quantum tunneling instability of the mythical hexazine and pentazine

Through computational analysis we found that pentazine and hexazine, two hypothetical high-energy density materials, exhibit inherent instability due to quantum tunneling effects. This instability remains even near the absolute zero, and therefore they can be deemed as unsynthesizable. We propose substituents that could potentially stabilize pentazine, especially dimethylamine.

Purely single-bonded spiral nitrogen chains stabilized by trivalent lanthanum ions

Inspired by the single-bonded nitrogen chains stabilized by tetravalent cerium, pentavalent tantalum, and hexavalent tungsten atoms, we explored the possibility of single-bonded nitrogen polymorphs stabilized by trivalent lanthanum ions. To achieve this, we utilized the crystal structure search method on the phase diagram of binary La-N compounds. We identified three novel thermodynamically stable phases, the C2/c LaN3, P-1 LaN4, and P-1 LaN8. Among them, the C2/c phase with infinite helical poly-N6 chains becomes thermodynamically stable above 50 GPa. Each nitrogen atom in the poly-N6 chain acquires one extra electron, and the spiral chain is purely single-bonded. The C2/c phase has an indirect band gap of ∼1.6 eV at 60 GPa. Notably, the band gap exhibits non-monotonic behavior, decreases first and then increases with increasing pressure. This abnormal behavior is attributed to the significant bonding of two La-N bonds at around 35 GPa. Phonon spectrum calculations and AIMD simulations have confirmed that the C2/c phase can be quenched to ambient conditions with slight distortion, and it exhibits excellent detonation properties. Additionally, we also discovered armchair-like nitrogen chains in LaN4 and the armchair and zigzag-like mixed nitrogen chains in LaN8. These results provide valuable insights into the electronic and bonding properties of nitrides under high pressure and may have important implications for the design and development of novel functional materials.

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure's fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features and interface this method with the XtalOpt evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g, the coordination number or chemical environment), and symmetry (e.g, Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN8, with a two-dimensional polymeric nitrogen sublattice, brookite TiO2, and a high pressure BaH4 phase, which was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1̅ WC3N6, was found to possess an energy that is lower than that of two phases proposed in a recent computational study. The method presented here could help in identifying the structures of compounds that have already been synthesized, and in developing new synthesis targets with desired properties.

Ambient-Pressure Recoverable Polynitrogen Solids Assembled by Pentazolate Rings with High Energy Density

Crystal structure predictions and first-principles calculations were used to predict three polynitrogen solids (aP8-N, aP12-N, and oP24-N) that possess competitive enthalpies as compared to the synthesized open-chain N₈ phase at pressures in the range of 0-60 GPa. aP8-N, aP12-N, and oP24-N contain edge-shared, N₂-linked, and N-bridged pentazolate rings and form molecular N₈, molecular N₁₂, and quasi-one-dimensional N_∞ ribbons, respectively. The calculations of formation enthalpies show that the three polynitrogen solids can be synthesized by compressing cyclo-N₅ salts in hydrogen-saturated environments. Molecular simulations suggest that the three polynitrogen solids have the ability of quench recoverability under ambient conditions once being synthesized at high pressure. With estimated energy densities in the range of 5.6-6.5 kJ/g, these three polynitrogen phases show notable promise for applications as high-energy-density materials.

Evidence for Dissociation and Ionization in Shock Compressed Nitrogen to 800 GPa

Triple bonding in the nitrogen molecule (N_{2}) is among the strongest chemical bonds with a dissociation enthalpy of 9.8  eV/molecule. Nitrogen is therefore an excellent test bed for theoretical and numerical methods aimed at understanding how bonding evolves under the influence of the extreme pressures and temperatures of the warm dense matter regime. Here, we report laser-driven shock experiments on fluid molecular nitrogen up to 800 GPa and 4.0  g/cm^{3}. Line-imaging velocimetry measurements and impedance matching method with a quartz reference yield shock equation of state data of initially precompressed nitrogen. Comparison with numerical simulations using path integral Monte Carlo and density functional theory molecular dynamics reveals clear signatures of chemical dissociation and the onset of L-shell ionization. Combining data along multiple shock Hugoniot curves starting from densities between 0.76 and 1.29  g/cm^{3}, our study documents how pressure and density affect these changes in chemical bonding and provides benchmarks for future theoretical developments in this regime, with applications for planetary interior modeling, high energy density science, and inertial confinement fusion research.

Theoretical methods for structural phase transitions in elemental solids at extreme conditions: statics and dynamics

In recent years, theoretical studies have moved from a traditionally supporting role to a more proactive role in the research of phase transitions at high pressures. In many cases, theoretical prediction leads the experimental exploration. This is largely owing to the rapid progress of computer power and theoretical methods, particularly the structure prediction methods tailored for high-pressure applications. This review introduces commonly used structure searching techniques based on static and dynamic approaches, their applicability in studying phase transitions at high pressure, and new developments made toward predicting complex crystalline phases. Successful landmark studies for each method are discussed, with an emphasis on elemental solids and their behaviors under high pressure. The review concludes with a perspective on outstanding challenges and opportunities in the field.

Low pressure metastable single-bonded solid nitrogen phases

Within the framework of the density functional theory, the possibility of the formation of single-bonded solid atomic nitrogen phases as a result of adiabatic compression of molecular and cluster nitrogen structures at zero temperature has been studied. It has been demonstrated that nitrogen clusters N₈(C_2v)-B, which are theoretically predicted as one of the promising candidates for high energy density materials, can transform under compression into a solid atomic phase with crystal lattice symmetry P2₁. The P2₁ phase is dynamically stable under decompression to zero pressure. It is shown that the ε-N₂ molecular phase transforms under compression into a solid atomic phase with R3̄c symmetry, and retains a vibrationally stable crystal structure when the pressure is reduced to 30 GPa, transforming into a stable cluster form at lower pressures. The atoms in the P2₁ and R3̄c solid atomic phases are linked by single bonds; therefore, these structures can store a large amount of energy ≈1.4 eV per atom. A detailed comparison of the properties of new P2₁ and R3̄c solid atomic phases with other nitrogen crystal structures that are dynamically stable at low pressures has been carried out.

Pressure-stabilized high-energy-density material YN10

Polynitrogen compounds have been intensively studied for potential applications as high energy density materials, especially in energy and military fields. Here, using the swarm intelligence algorithm in combination with first-principles calculations, we systematically explored the variable stoichiometries of yttrium-nitrogen compounds on the nitrogen-rich regime at high pressure, where a new stable phase of YN₁₀adoptingI4/msymmetry was discovered at the pressure of 35 GPa and showed metallic character from the analysis of electronic properties. In YN₁₀, all the nitrogen atoms weresp²-hybridized in the form of N₅ring. Furthermore, the gravimetric and volumetric energy densities were estimated to be 3.05 kJ g^-1and 9.27 kJ cm^-1respectively. Particularly, the calculated detonation velocity and pressure of YN₁₀(12.0 km s^-1, 82.7 GPa) was higher than that of TNT (6.9 km s^-1, 19.0 GPa) and HMX (9.1 km s^-1, 39.3 GPa), making it a potential candidate as a high-energy-density material.

Materials by design at high pressures

Pressure, a fundamental thermodynamic variable, can generate two essential effects on materials. First, pressure can create new high-pressure phases via modification of the potential energy surface. Second, pressure can produce new compounds with unconventional stoichiometries via modification of the compositional landscape. These new phases or compounds often exhibit exotic physical and chemical properties that are inaccessible at ambient pressure. Recent studies have established a broad scope for developing materials with specific desired properties under high pressure. Crystal structure prediction methods and first-principles calculations can be used to design materials and thus guide subsequent synthesis plans prior to any experimental work. A key example is the recent theory-initiated discovery of the record-breaking high-temperature superhydride superconductors H3S and LaH10 with critical temperatures of 200 K and 260 K, respectively. This work summarizes and discusses recent progress in the theory-oriented discovery of new materials under high pressure, including hydrogen-rich superconductors, high-energy-density materials, inorganic electrides, and noble gas compounds. The discovery of the considered compounds involved substantial theoretical contributions. We address future challenges facing the design of materials at high pressure and provide perspectives on research directions with significant potential for future discoveries.

Polymerization of Nitrogen in Nitrogen-Fluorine Compounds under Pressure

A wide range of polynitrogen species have attracted much attention because of their potential applications as high-energy-density materials. Until now, predicted polynitrogen was found to be negatively charged, with charge transfer from introduced atoms to nitrogen in nitrogen-bearing compounds. Using an evolutionary algorithm combined with first-principles calculations, stoichiometries and structures in nitrogen-fluorine compounds at pressures ranging from 0 to 200 GPa are investigated. In addition to two fluorine-rich compounds NF3 and NF5, two other compounds, NF and N6F, emerge with increasing pressure. N6F, as a nitrogen-rich compound, will become stable at pressures greater than 180 GPa with a positively charged nitrogen network. Above 120 GPa, the NF compound with polymeric zigzag nitrogen chains is discovered, and it is quenchable to the ambient conditions, acquiring the highest energy density of 5.38 kJ/g among reported binary covalent nitrogen compounds. These newly predicted N-F compounds are useful in understanding the chemistry of polynitrogen.

Multistep Dissociation of Fluorine Molecules under Extreme Compression

All elements that form diatomic molecules, such as H_{2}, N_{2}, O_{2}, Cl_{2}, Br_{2}, and I_{2}, are destined to become atomic solids under sufficiently high pressure. However, as revealed by many experimental and theoretical studies, these elements show very different propensity and transition paths due to the balance of reduced volume, lone pair electrons, and interatomic bonds. The study of F under pressure can illuminate this intricate behavior since F, owing to its unique position on the periodic table, can be compared with H, with N and O, and also with other halogens. Nevertheless, F remains the only element whose solid structure evolution under pressure has not been thoroughly studied. Using a large-scale crystal structure search method based on first principles calculations, we find that, before reaching an atomic phase, F solid transforms first into a structure consisting of F_{2} molecules and F polymer chains and then into a structure consisting of F polymer chains and F atoms, a distinctive evolution with pressure that has not been seen in any other elements. Both intermediate structures are found to be metallic and become superconducting, a result that adds F to the elemental superconductors.

Rules of formation of H-C-N-O compounds at high pressure and the fates of planetary ices

The solar system's outer planets, and many of their moons, are dominated by matter from the H-C-N-O chemical space, based on solar system abundances of hydrogen and the planetary ices [Formula: see text]O, [Formula: see text], and [Formula: see text] In the planetary interiors, these ices will experience extreme pressure conditions, around 5 Mbar at the Neptune mantle-core boundary, and it is expected that they undergo phase transitions, decompose, and form entirely new compounds. While temperature will dictate the formation of compounds, ground-state density functional theory allows us to probe the chemical effects resulting from pressure alone. These structural developments in turn determine the planets' interior structures, thermal evolution, and magnetic field generation, among others. Despite its importance, the H-C-N-O system has not been surveyed systematically to explore which compounds emerge at high-pressure conditions, and what governs their stability. Here, we report on and analyze an unbiased crystal structure search among H-C-N-O compounds between 1 and 5 Mbar. We demonstrate that simple chemical rules drive stability in this composition space, which explains why the simplest possible quaternary mixture HCNO-isoelectronic to diamond-emerges as a stable compound and discuss dominant decomposition products of planetary ice mixtures.

Properties of Alkaline-Earth-Metal Polynitrogen Ternary Materials at High Pressure

We report the formation of cubic and tetragonal BaSrN3 at 100 GPa using an ab initio structure search method. Pressure ramping to 0 GPa reveals a reaction pressure threshold of 4.92 and 7.23 GPa for the cubic and tetragonal BaSrN3, respectively. The cubic phase is stabilized by coulombic interaction between the ions. Meanwhile, tetragonal BaSrN3 is stabilized through an expansion of the d-orbital in Ba and Sr atoms that is compensated by delocalization of π-electrons in N through reduction of π overlap. Elastic properties analysis suggests that both phases are mechanically stable. The structures also have high melting points as predicted using an empirical model, and all imaginary modes vanishes at about 2000 K. These results have significant implication for the design of cleaner and environmentally friendly high energy density materials.

High-Pressure Polymeric Nitrogen Allotrope with the Black Phosphorus Structure

Studies of polynitrogen phases are of great interest for fundamental science and for the design of novel high energy density materials. Laser heating of pure nitrogen at 140 GPa in a diamond anvil cell led to the synthesis of a polymeric nitrogen allotrope with the black phosphorus structure, bp-N. The structure was identified in situ using synchrotron single-crystal x-ray diffraction and further studied by Raman spectroscopy and density functional theory calculations. The discovery of bp-N brings nitrogen in line with heavier pnictogen elements, resolves incongruities regarding polymeric nitrogen phases and provides insights into polynitrogen arrangements at extreme densities.

Predictions on High-Power Trivalent Metal Pentazolate Salts

High-energy-density materials (HEDMs) have been intensively studied for their significance in fundamental sciences and practical applications. Here, using the molecular crystal structure search method based on first-principles calculations, we have predicted a series of metastable energetic trivalent metal pentazolate salts MN₁₅ (M= Al, Ga, Sc, and Y). These compounds have high energy densities, with the highest nitrogen content among the studied nitrides so far. Pentazolate N₅^- molecules stack up face-to-face and form wave-like patterns in the C222₁ and Cc symmetries. The strong covalent bonding and very weak noncovalent interactions with nonbonded overlaps coexist in these ionic-like structures. We find MN₁₅ molecular structures are mechanically stable up to high temperature (∼1000 K) and ambient pressure. More importantly, these trivalent metal pentazolate salts have high detonation pressure (∼80 GPa) and velocity (∼12 km/s). Their detonation pressures exceeding that of TNT and HMX make them good candidates for high-brisance green energetic materials.

Metallic and anti-metallic properties of strongly covalently bonded energetic AlN5 nitrides

High pressure can stimulate numerous novel physical effects which are not observed under ambient conditions, such as the electronic redistribution and delocalization phenomenon in strongly covalently bonded nitrides. Through first principles simulations, we report a new N-rich aluminum nitride AlN5, which crystallizes with the space group P1[combining macron] at 20 GPa and then transforms into the I4[combining macron]2d phase at 60 GPa. We have identified and proved the delocalization effects of π electrons in the strongly covalent Lewis poly-nitrogen structure via the one-dimensional particle in a box mechanism, which contributes to the metallization and stability of the system. This implies that not all strongly covalently bonded systems with highly localized electrons exhibit nonmetallic properties in III-V main group nitrides. Furthermore, pressure results in the hybridization configuration mutation from sp2 in the P1[combining macron] phase to a mixture of sp2 and sp3 hybridization in the I4[combining macron]2d phase, which leads to phase transition from metal to insulator. With increasing pressure, the band gap increases abnormally, exhibiting anti-metallization induced by the strong hybridization. Interestingly, the P1[combining macron] and I4[combining macron]2d structures are simultaneously accompanied by a high energy density and hardness, which enable them to have a greater ability to resist elasticity, plastic deformation and external force destruction in potential applications. Their energy density and hardness are up to 3.29 kJ g-1 and 15.2 GPa in the P1[combining macron] phase but especially 6.14 kJ g-1 and 31.7 GPa in the I4[combining macron]2d phase.

Hexagonal Layered Polymeric Nitrogen Phase Synthesized near 250 GPa

The nitrogen triple bond dissociates in the 100 GPa pressure range and a rich variety of single-bonded polymeric nitrogen structures unique to this element have been predicted up to the terapascal pressure range. The nonmolecular cubic-gauche (cg-N) structure was first observed above 110 GPa, coupled to high temperature (>2000  K) to overcome the kinetic barrier. A mixture of cg-N with a layered phase was afterwards reported between 120 and 180 GPa. Here, by laser heating pure nitrogen from 180 GPa, a sole crystalline phase is characterized above 240 GPa while an amorphous transparent phase is obtained at pressures below. X-ray diffraction and Raman vibrational data reveal a tetragonal lattice (P4_{2}bc) that matches the predicted hexagonal layered polymeric nitrogen (HLP-N) structure. Density-functional theory calculations which include the thermal and dispersive interaction contributions are performed to discuss the stability of the HLP-N structure.

Origin of the considerably high thermal stability of cyclo-N5− containing salts at ambient conditions

Ionization, conjugation, hydrogen bonding, coordination bonding and π–π stacking consolidate the cyclo-N₅⁻ caged in salt crystals.

Structural prediction of host-guest structure in lithium at high pressure

Ab initio random structure searching (AIRSS) technique is used to identify the high-pressure phases of lithium (Li). We proposed the transition mechanism from the fcc to host-guest (HG) structures at finite temperature and high pressure. This complex structural phase transformation has been calculated using ab initio lattice dynamics with finite displacement method which confirms the dynamical harmonic stabilization of the HG structure. The electron distribution between the host-host atoms has also been investigated by electron localization function (ELF). The strongly localized electron of p bond has led to the stability of the HG structure. This remarkable result put the HG structure to be a common high-pressure structure among alkali metals.

Route to high-energy density polymeric nitrogen t-N via He-N compounds

Polymeric nitrogen, stabilized by compressing pure molecular nitrogen, has yet to be recovered to ambient conditions, precluding its application as a high-energy density material. Here we suggest a route for synthesis of a tetragonal polymeric nitrogen, denoted t-N, via He-N compounds at high pressures. Using first-principles calculations with structure searching, we predict a class of nitrides with stoichiometry HeN₄ that are energetically stable (relative to a mixture of solid He and N₂) above 8.5 GPa. At high pressure, HeN₄ comprises a polymeric channel-like nitrogen framework filled with linearly arranged helium atoms. The nitrogen framework persists to ambient pressure on decompression after removal of helium, forming pure polymeric nitrogen, t-N. t-N is dynamically and mechanically stable at ambient pressure with an estimated energy density of ~11.31 kJ/g, marking it out as a remarkable high-energy density material. This expands the known polymeric forms of nitrogen and indicates a route to its synthesis.

Melting curve of SiO2 at multimegabar pressures: implications for gas giants and super-Earths

Ultrahigh-pressure phase boundary between solid and liquid SiO₂ is still quite unclear. Here we present predictions of silica melting curve for the multimegabar pressure regime, as obtained from first principles molecular dynamics simulations. We calculate the melting temperatures from three high pressure phases of silica (pyrite-, cotunnite-, and Fe₂P-type SiO₂) at different pressures using the Z method. The computed melting curve is found to rise abruptly around 330 GPa, an increase not previously reported by any melting simulations. This is in close agreement with recent experiments reporting the α-PbO₂–pyrite transition around this pressure. The predicted phase diagram indicates that silica could be one of the dominant components of the rocky cores of gas giants, as it remains solid at the core of our Solar System’s gas giants. These results are also relevant to model the interior structure and evolution of massive super-Earths.

Diverse Chemistry of Stable Hydronitrogens, and Implications for Planetary and Materials Sciences

Nitrogen hydrides, e.g., ammonia (NH3), hydrazine (N2H4) and hydrazoic acid (HN3), are compounds of great fundamental and applied importance. Their high-pressure behavior is important because of their abundance in giant planets and because of the hopes of discovering high-energy-density materials. Here, we have performed a systematic investigation on the structural stability of N-H system in a pressure range up to 800 GPa through evolutionary structure prediction. Surprisingly, we found that high pressure stabilizes a series of previously unreported compounds with peculiar structural and electronic properties, such as the N4H, N3H, N2H and NH phases composed of nitrogen backbones, the N9H4 phase containing two-dimensional metallic nitrogen planes and novel N8H, NH2, N3H7, NH4 and NH5 molecular phases. Another surprise is that NH3 becomes thermodynamically unstable above ~460 GPa. We found that high-pressure chemistry of hydronitrogens is much more diverse than hydrocarbon chemistry at normal conditions, leading to expectations that N-H-O and N-H-O-S systems under pressure are likely to possess richer chemistry than the known organic chemistry. This, in turn, opens a possibility of nitrogen-based life at high pressure. The predicted phase diagram of the N-H system also provides a reference for synthesis of high-energy-density materials.

Pressure-induced localisation of the hydrogen-bond network in KOH-VI

Using a combination of ab initio crystal structure prediction and neutron diffraction techniques, we have solved the full structure of KOH-VI at 7 GPa. Rather than being orthorhombic and proton-ordered as had previously be proposed, we find that this high-pressure phase of potassium hydroxide is tetragonal (space group I4/mmm) and proton disordered. It has an unusual hydrogen bond topology, where the hydroxyl groups form isolated hydrogen-bonded square planar (OH)4 units. This structure is stable above 6.5 GPa and, despite being macroscopically proton-disordered, local ice rules enforce microscopic order of the hydrogen bonds. We suggest the use of this novel type of structure to study concerted proton tunneling in the solid state, while the topology of the hydrogen bond network could conceivably be exploited in data storage applications based solely on the manipulations of hydrogen bonds. The unusual localisation of the hydrogen bond network under applied pressure is found to be favored by a more compact packing of the constituents in a distorted cesium chloride structure.

Transformation pathways in high-pressure solid nitrogen: from molecular N2 to polymeric cg-N

The transformation pathway in high-pressure solid nitrogen from N2 molecular state to polymeric cg-N phase was investigated by means of ab initio molecular dynamics and metadynamics simulations. In our study, we observed a transformation mechanism starting from molecular Immm phase that initiated with formation of trans-cis chains. These chains further connected within layers and formed a chain-planar state, which we describe as a mixture of two crystalline structures--trans-cis chain phase and planar phase, both with Pnma symmetry. This mixed state appeared in molecular dynamics performed at 120 GPa and 1500 K and in the metadynamics run at 110 GPa and 1500 K, where the chains continued to reorganize further and eventually formed cg-N. During separate simulations, we also found two new phases--molecular P2(1)/c and two-three-coordinated chain-like Cm. The transformation mechanism heading towards cg-N can be characterized as a progressive polymerization process passing through several intermediate states of variously connected trans-cis chains. In the final stage of the transformation chains in the layered form rearrange collectively and develop new intraplanar as well as interplanar bonds leading to the geometry of cg-N. Chains with alternating trans and cis conformation were found to be the key entity--structural pattern governing the dynamics of the simulated molecular-polymeric transformation in compressed nitrogen.

Pressure-induced symmetry-lowering transition in dense nitrogen to layered polymeric nitrogen (LP-N) with colossal Raman intensity

We present the discovery of a novel nitrogen phase synthesized using laser-heated diamond anvil cells at pressures between 120-180 GPa well above the stability field of cubic gauche (cg)-N. This new phase is characterized by its singly bonded, layered polymeric (LP) structure similar to the predicted Pba2 and two colossal Raman bands (at ∼1000 and 1300  cm^{-1} at 150 GPa), arising from two groups of highly polarized nitrogen atoms in the bulk and surface of the layer, respectively. The present result also provides a new constraint for the nitrogen phase diagram, highlighting an unusual symmetry-lowering 3D cg-N to 2D LP-N transition and thereby the enhanced electrostatic contribution to the stabilization of this densely packed LP-N (ρ=4.85  g/cm^{3} at 120 GPa).