Nitrogen oxides under pressure: stability, ionization, polymerization, and superconductivity.

Nitrogen oxides are textbook class of molecular compounds, with extensive industrial applications. Nitrogen and oxygen are also among the most abundant elements in the universe. We explore the N-O system at 0 K and up to 500 GPa though ab initio evolutionary simulations. Results show that two phase transformations of stable molecular NO2 occur at 7 and 64 GPa, and followed by decomposition of NO2 at 91 GPa. All of the NO(+)NO3(-) structures are found to be metastable at T = 0 K, so experimentally reported ionic NO(+)NO3(-) is either metastable or stabilized by temperature. N2O5 becomes stable at 9 GPa, and transforms from P-1 to C2/c structure at 51 GPa. NO becomes thermodynamically stable at 198 GPa. This polymeric phase is superconducting (Tc = 2.0 K) and contains a -N-N- backbone.

Both nitrogen and oxygen have been extensively investigated in experiments and theoretical simulations. Generally, nitrogen is an insulator or a semiconductor. Cubic gauche phase of nitrogen1 is stable in a wide range of pressure2. Other nitrogen structures, such as chain and rings345, have also been reported. All known phases of oxygen are molecular67. Experiments and first-principles calculations for oxygen under high pressure revealed the complex evolution of insulator->semiconductor->metal->semiconductor8. The superconductivity of solid oxygen (T = 0.6 K) was observed at pressure above 96 GPa in experiment9. The known nitrogen oxides are semiconducting (for example, the band gap of Im-3 NO2 calculated is approximately 2.8 eV).

At ambient pressure, nitrogen oxides exist as molecular crystals with many applications in chemical industry and important biological roles. The volumetric behavior of nitrous oxide under pressure has been investigated since 196110. The synthesis and phase transformations of N2O have been analyzed in experimental and theoretical studies111213141516. Different from normal phases containing N2O4 molecules, the ionic NO+NO3− was reported in the range of 1.5 to 3.0 GPa17. The typical N-O stretching frequency of NO+ was characterized at 2234 cm−1, consistent with previous reports1718. In 2001, Somayazulu et al.19 synthesized the ionic NO+NO3− (nitrosonium nitrate) phase from N2O at above 20 GPa and 1000 K, and performed first structural characterization of NO+NO3−. Somayazulu et al.19 proposed an ionic NO+NO3− model based on aragonite with space group of P21cn. Other P21/m20 and Pna2116 models of NO+NO3− were also suggested, and the later one is more stable. However, the simulated XRD data of the Pna21 structure is quite different from that in experiments, indicating that other undiscovered stable NO+NO3− structures might exist.

Results and Discussions

We employed the evolutionary algorithm USPEX21222324 to predict stable N-O compounds and structures under high pressures. Up to 500 GPa, only three stable N-O compounds were found (NO2, N2O5 and NO), as seen in Fig. 1. Most of them retain their molecular structures even under high pressure. Experimentally known “laughing gas” N2O is metastable. The stable phases are discussed as follows.

Figure 1

Phase diagram of the N–O system.

NO2: Besides the known cubic (Im-3) and monoclinic (P21/c) NO2 structures are stable in pressure ranges of 0–7 and 7–64 GPa respectively, another P21/c structure was found to be stable from 64 to 91 GPa (Fig. 2a). Similar to the known phases, this novel NO2 structure also contains N2O4 molecules. Different from known P21/c NO2, the proposed P21/c NO2 is denser and has 8 formula units in the unit cell. NO2 decomposes at 91 GPa.

Figure 2

Structures of (a) P21/c NO2, (b) P-1 and (c) C2/c N2O5, (d) P21/m NO.

N2O5: Molecular N2O5 phases are stable in a wide pressure range (9–446 GPa). At 51 GPa, N2O5 transforms from P-1 (Fig. 2b) to C2/c (Fig. 2c) structure. The N2O5 molecules remain planar. At 446 GPa, N2O5 becomes unstable and decomposes into NO and O. At 0 GPa, the P-1 N2O5 is 0.04 eV/atom more stable than known hexagonal NO2NO3, however, both of them are calculated to be above the convex hull, and therefore metastable.

NO: NO is a metastable compound at low pressures. A polymeric NO structure (Fig. 2d, P21/m) becomes stable at 198 GPa. Nitrogen atoms form a strong covalent backbone (N-N = 1.34 Å) in the shape of a zigzag chain. Indeed, that is right between the typical values of single (1.45 Å) and double (1.25 Å) nitrogen-nitrogen bonds. Each nitrogen atom is also bonded to one oxygen atom (N-O bond length is 1.20 Å). A similar backbone has also been reported for the N-H system25. Distance between neighboring quasi-one-dimensional structures is 1.86 Å. Phonon dispersion curves of this remarkable polymeric phase were calculated (shown in Fig. S6). No imaginary frequencies were observed, implying its dynamical stability.

While most of the stable N-O phases are semiconducting, polymeric NO is metallic. The band structure of NO is shown in Fig. 3. Using the Allen-Dynes modified McMillan equation2627 with value of the Coulomb pseudopotential μ* = 0.13, polymeric NO is superconducting with T = 2.0 K at 200 GPa, which is close to that of oxygen89.

Figure 3

Band structure of P21/m NO at 198 GPa. Z(0,0,0.5), A(0.5,0.5,0.5), M(0.5,0.5,0), G(0,0,0), R(0,0.5,0.5) and X(0,0.5,0).

As mentioned above, ionic NO+NO3- has been observed in several high-pressure experiments192028. However, no stable NO+NO3− structure was found in our variable-composition searches. To find the lowest-enthalpy ionic NO+NO3− structure, we performed (NO)n(NO3)n (n = 6 or 8) calculations at 0–50 GPa, assembling structures from ready-made NO and NO3 units in variable proportion. A novel metastable monoclinic NO+NO3− (P21, Fig. 4) was found to be more stable than orthorhombic phase16 and monoclinic P21/m NO+NO3− 20 at pressures above 1.7 GPa. The main difference between novel P21 and P21/m NO+NO3− models20 is the orientation of the NO+ molecules. Importantly, at all pressures structures made of N2O4 molecules are more stable than ionic NO+NO3− structures (Fig. 4)

Figure 4

Enthalpies of NO2 phases as a function of pressures.

In experiments, the typical Raman frequencies of NO+NO3− are 2234 cm−1 for the N-O stretch in NO+, together with 1345, 1056 and 721 cm−1 for anti-symmetric stretch, symmetric stretch and in-plane deformation for NO3− respectively1719282930. The Raman frequencies and intensities of NO+NO3− and NO2 structures were calculated at 20 GPa. Here, Raman frequencies of NO+ and NO3− were used for comparison. As shown in Fig. 5, the computed Raman spectra of P21/c NO2 and Pna21 NO+NO3− are significantly different from experimental ones.17. The typical Raman frequencies of N-O stretch are 2071 cm−1 of P21 and 2151 cm−1 of P21/m structures. Both of them basically match the experimental data19, but that of P21 NO+NO3− obtains better match in terms of relative intensity. Similar comparison for Raman and XRD data could also be observed under other pressures2031.

Figure 5

Simulated Raman spectra of P21/c NO2, and P21, P21/m20 and Pna2116 NO+NO3− at 20 GPa.

Typical Raman frequencies of NO+ and NO3− in experiment17 were drawn by dotted lines.

In summary, stable NO2, N2O5 and NO phases were found in N-O system up to 500 GPa. The P21/c NO2 becomes stable at 64 GPa and decomposes at 91 GPa. N2O5 with P-1 becomes stable at 9 GPa, transforms to C2/c at 51 GPa and decomposes at 446 GPa. The only metallic structure (P21/m NO) has -N-N- zigzag backbone and possesses superconductivity with T = 2.0 K. Our results show that ionic NO+NO3− is metastable, and we identify a novel P21 structure that matches experimental data better and has lower enthalpy than previously proposed structures.

Methods

An evolutionary algorithm, as implemented in the USPEX code21222324, were utilized to search for the stable compounds and structures. This method has already been successfully applied to study numerous systems, including nitrogen and oxygen under pressure57. Structure relaxations were done using density functional theory (DFT)3233 within the generalized gradient approximation (GGA)34 using the all-electron projector augmented wave (PAW)3536 method as implemented in the VASP code37. The plane-wave kinetic energy cutoff was set to 600 eV and Brillouin zone was sampled at a resolution of 2π×0.06 Å−1. At first, variable-composition were carried out at 0, 10, 20, 30, 50, 100, 150, 200, 250, 300, 350, 400, 450 and 500 GPa. Stability of compounds was judged using the convex hull construction: those compounds which are on the convex hull (i.e. which are more favorable than any isochemical mixture of other phases) are thermodynamically stable at given conditions. The PHONOPY code38 was employed to calculate phonon dispersions for all promising structures, and all the discussed structures were found to be dynamically stable. All Raman frequencies and intensities were calculated according to the method of Porezag and Pederson39. The electron–phonon coupling calculations in Quantum Espresso40 with 180 Ry plane-wave cutoff energy were used to calculate the critical temperature of superconductivity (Tc).

Additional Information

How to cite this article: Li, D. et al. Nitrogen oxides under pressure: stability, ionization, polymerization, and superconductivity. Sci. Rep. 5, 16311; doi: 10.1038/srep16311 (2015).