Pressure-Induced Polymerization of One-Dimensional Nitrogen Chains in K2N2
doi: 10.11858/gywlxb.20240719
-
Abstract: The crystal structure prediction of K2N2 in the pressure range of 0–150 GPa using an advanced particle swarm crystal structure search method was conducted. The results show that the stable ground state phase of K2N2 is a monoclinicC2/mstructure, and three high-pressure structures including Na2N2-type,Cmmm, andC2/care identified at pressures of 1.7, 3.6, 122 GPa, respectively. The volume dependence on pressure shows that the three phase transitions, i. e.,C2/m→Na2N2-type, Na2N2-type→Cmmm, andCmmm→C2/c, are all first order phase transitions, corresponding to volume collapses of 14.4%, 22.5%, and 4.0%, respectively. During the high pressure phase transitions of K2N2, the coordination number of K atom increases from 5 to 10, and a change in the nature of the N-N bonding from N=N dimmer in the ground state ofC2/mstructure to N―N single bond chain in the high-pressureC2/cphase is accompanied. The high-pressureC2/cphase exhibits semiconducting properties with a band gap of 2.0 eV, whileC2/m, Na2N2-type, andCmmmphases have metallic behaviors. Electronic structure calculation and electron-localized function analysis indicate that the high-pressure structural phase transition of K2N2 is due to the K-plone-pair electrons activation and their participation in bonding with N atoms under high pressure.
-
Key words:
- high pressure /
- K2N2 /
- structure predictions /
- phase transitions /
- electronic structures
摘要: 采用先进的粒子群晶体结构搜索方法对K2N2在0~150 GPa压强范围内进行晶体结构预测,结果表明,K2N2的基态稳定相为单斜C2/m结构,且在1.7、3.6和122 GPa压强下的结构分别为Na2N2型、Cmmm、C2/c。体积随压强的变化关系显示C2/m→Na2N2型、Na2N2型→Cmmm和Cmmm→C2/c这3个相变均为一级相变,对应的体积坍塌分别为14.4%、22.5%和4.0%。在K2N2高压相变过程中,K原子的配位数从5增加到10,并伴随着N-N成键性质的变化,即从基态C2/m结构中的准分子N=N双键聚合为高压C2/c相中的N―N单键链。C2/m、Na2N2型、Cmmm相表现出金属性,而高压C2/c相表现出半导体(带隙为2.0 eV)性质。电子结构计算和电子局域函数分析表明,K2N2的高压结构相变来源于高压下K-p孤对电子的激活及其与N原子的成键。 -
Metal nitrides have attracted increasing attention in a wide range application due to their fascinating properties, such as hardness, superconductivity, various types of magnetism, etc.[1–2] Recently, extensive efforts have focused on the synthesis and design of metal polynitrides with single and double nitrogen bonds under high pressure, as metal polynitrides are proposed to be promissing candidates for high-energy density materials (HEDMs)[3–5]. The decomposition of these polynitrides is expected to release huge energy due to the conversion from a single/double-bonded polymeric nitrogen (160/418 kJ/mol) into a triple-bonded nitrogen (954 kJ/mol) molecule without producing pollutants. The atomic polymeric nitrogen with cubic gauche structure (cg-N) at a high pressure and high temperature (HPHT) condition (pressure above 110 GPa and temperature greater than 2 000 K) was successfully synthetized by Eremets et al.[6–8]. Motivated by that synthesis, alkali metal azides AN3, constructed by spherical cations and linear molecular
N−3 anions with N=N double bonds, was proposed to be suitable precursor in the formation of polymeric nitrogen under pressure[9]. Experimentally, the polymeric evolutionary behaviors ofN−3 anions in LiN3[10], NaN3[11–13], KN3[14–16], RbN3[17–18], and CsN3[19–20] have been fully investigated under high pressure. One of the most interesting findings is the reported polymeric nitrogen nets in NaN3[11] at 160 GPa. More recently, a fully sp2-hybridized layered polymeric nitrogen structure, featuring fused 18-membered rings, in potassium supernitride (K2N16) was successfully synthesized under HPHT condition using a laser heating diamond anvil cell (DAC) technique[21]. In our previous works[22–25], theN−3 anions in LiN3, NaN3, and KN3 all transform firstly to “N6” molecular clusters under compression, and then to a polymerized nitrogen phase at high pressures.Besides the well-known alkali metal azides AN3, alkali metal diazenides A2N2 containing the N=N double bonds, were also proposed as precursors of the polymeric atomic nitrogen under pressure. In 2010, two alkali metal diazenides, Li2N2 and Na2N2, were proposed to have orthorhombic Pmmm structure[26]. However, an orthorhombic Immm structure of Li2N2[27] was then synthesized under HPHT condition. Immediately after that, a serial pressure-induced structural phase transitions of Li2N2 were reported in two independent works[28–29]. Moreover, a high-pressure tetragonal I41/acd phase containing the spiral nitrogen chain was uncovered at 242 GPa[29], indicating that the polymerization of nitrogen is realized from Li2N2.
Compared to Li2N2 and Na2N2, K2N2 have not been reported so far, neither for ambient pressure nor for high pressure. Here, we will extensively explore the high-pressure structures of K2N2 up to 150 GPa by using the first principles swarm structure searching method. The structure evolutions and chemical bonding behaviors of K2N2 within different structures were then fully studied to provide an insight into the formation of polymeric nitrogen in alkali metal diazenides.
1. Computational Method
The variable-cell crystal structure predictions for K2N2 (1–4 formula unit (f. u.) in the simulation cell) up to 150 GPa were conducted by CALYPSO code[30–31]. The effectiveness of this method has been fully confirmed in predicting high-pressure structures of various systems[32–35]. During the structure search process, the 60% structures of each generation with lower enthalpies were selected to generate the structures for the next generation by particle swarm optimization (PSO) technique, and the other structures in new generation were randomly generated to increase the structural diversity. Usually, the structure searching simulation stopped at the 20th generation and 40 structures per generation were generated at each pressure point. After finishing the structure search, the most stable structures were used for further structural relaxation and property calculation, a process that is implemented in the Vienna ab Initio Simulation Package[36]. The Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional under the generalized gradient approximation[37] and the projector-augmented wave (PAW)[38] for ion-electron potentials were adopted herein. The van-der-Waals (vdW) correction was applied in the structural relaxation using Grimme’s method (DFT-D3)[39]. The cutoff energy of 550 eV and Monkhorst-Pack k meshes[40] of 2π×0.032 Å−1 were used to ensure the good convergences of total energy and maximum force on each atom. The lattice dynamic stability of each structure was checked by phonon spectrum obtained by PHONOPY code[41]. The charge transfer behavior of K2N2 was performed using Bader atoms-in-molecules (AIM) method[42].
2. Results and Discussion
2.1 Ground-State and High-Pressure Structures of K2N2
The ground-state structure for K2N2 at ambient pressure in our fully structural search (in Fig. 1(a)) is a monoclinic structure with the space group C2/m, which is different from those reported for Li2N2 and Na2N2[26]. In Fig. 1(a), the N2 quasimolecule dimers in C2/m structure are parallel to each other, and are tilted by a certain angle within a-b plane. The N-N bond length is 1.192 Å, which is a typical value for double bond. Also, the calculated Mayer bond order[43] for the N2 anion in the C2/m is 2.18, being consistent with double bond. Due to the weak interaction between the N2 quasimolecules in the C2/m phase, the application of gigapascal pressure leads to the orientation change of N2 quasimolecules and the structure transformation. Indeed, two stable orthorhombic structures, i.e. Pmmm and Cmmm, for K2N2 were uncovered at 2 and 5 GPa, respectively, as shown in Fig. 1(b) and Fig. 1(c). Remarkably, the unveiled Pmmm phase and the previously reported ground-state structure of Na2N2 are exactly isostructrual (the Pmmm phase denoted as the Na2N2-type phase hereafter), and the Cmmm phase of K2N2 and the high-pressure metastable structures of Na2N2 are isostructrual as well[26]. Furthermore, both Na2N2-type and Cmmm orthorhombic phases have the similar atomic arrangement that the parallel N2 dimers are perpendicular to the K atomic layer. In the Na2N2-type and Cmmm phases, the calculated N-N bond lengths are 1.220 and 1.254 Å, while the characteristics of the N=N double bond are as that of the C2/m structure in the ground-state. As the pressures rises to 150 GPa, a stable monoclinic C2/c structure, containing a folded one-dimensional nitrogen chains (see Fig. 1(d)), is identified for K2N2 for the first time. The lengths of two nonequivalent N-N bond in C2/c structure are 1.474 and 1.373 Å at 135 GPa, indicating genuine N―N single bond, and thus the nature of the structure is similar to that of the high-pressure phase of NaN3[24]. Table 1 lists the structural details of the lattice parameters and the atomic positions of various K2N2 structures at selected pressure points.
Table 1. Optimized structural parameters of the C2/m, Na2N2-type, Cmmm, and C2/c phases of K2N2Phase Pressure/GPa Lattice parameters dN-N/Å Atomic fractional coordinates C2/m 0 a = 7.562 Å, b = 3.912 Å, c = 10.923 Å,
α = γ = 90°, β= 134.955°1.192 K 4i (0.613, 0, 0.260)
N 4i (0.519, 0, 0.456)Na2N2-type 2.5 a = 3.326 Å, b = 4.375 Å, c = 5.694 Å,
α = β = γ = 90°1.220 K1 1e (0, 0.500, 0)
K2 1c (0, 0, 0.500)
N 2t (0.500, 0.500, 0.393)Cmmm 20.0 a = 6.896 Å, b = 5.104 Å, c = 2.855 Å,
α = β = γ = 90°1.254 K 4g (0.694, 0, 0)
N 4j (0, 0.877, 0.500)C2/c 135.0 a = 6.971 Å, b = 4.099 Å, c = 4.510 Å,
α = γ = 90°, β= 81.342°1.474, 1.373 K 8f (0.669, 0.909, 0.198)
N 8f (0.041, 0.887, 0.102)The phonon dispersion curves of the four predicted structures are illustrated in Fig. 2 in order to check their dynamic stabilities. Fig. 2 shows that all the K2N2 phases are dynamically stable as the eigen frequencies of their lattice vibrations are positive for all wavevectors in the Brillouin zone. It should be noted that the predicted C2/m, Na2N2-type, and Cmmm phases exhibit higher eigen frequencies of around 50 THz, which correspond to the vibrations of N2 quasimolecules with stronger N=N double bonds. By contrast, the C2/c phase shows relative lower vibration frequencies of around 30−35 THz, corresponding to the vibration frequencies of weaker N―N single bonds.
Herein, we evaluated the energy involved in the decomposition of the C2/c phase with polymerized N―N single bond chains into metal K and gaseous N2. The energy involved in the C2/c-K2N2 decompostion was determined to be 0.344 kJ/g, which is about three times of that of conventional energetic materials (0.104−0.105 kJ/g)[44], but much smaller than those reported for high-energy density materials (3.61 kJ/g for BeN4, 5.30 kJ/g for GaN5, 4.81 kJ/g for GaN10)[45–46]. This may be due to the low nitrogen content in K2N2, so the nitrogen-rich KNx (x> 1) is expected to be potential HEDM under high pressure.
2.2 Pressure-Induced Structural Phase Transitions of K2N2
To determine the phase transition pressure of K2N2, we plotted in Figs. 3(a) and 3(b) the enthalpy difference of the predicted C2/m, Cmmm, and C2/c structures relative to the Na2N2-type phase in the pressure range of 0−150 GPa. The inset of Fig. 3(a) shows that above 1.7 GPa the Na2N2-type phase becomes energetically preferable with respect to the ground-state C2/m structure and is enthalpically stable up to about 3.6 GPa, above which it transforms into the Cmmm structure. From Fig. 3(b), it can be seen that the Cmmm structure persists up to 122 GPa and then transforms to the C2/c structure according to our structural predictions and enthalpy difference calculations. The wide pressure range of stable Cmmm phase may be attributed to the relative smaller ionic radius and atomic mass of K atom,in contrast with the heavier Rb and Cs atoms. It is thus expected that the chemical precompression exerted by heavier atoms, such as Rb and Cs, in A2N2 (as for Cs2N2[43]) can yield lower pressures for the polymerization of one-dimension N chains. From the calculated pressure dependence of the volume per f. u. in Figs. 3(c) and 3(d), the volume collapses of the first-order transformations C2/m→Na2N2-type, Na2N2-type→Cmmm, and Cmmm→C2/c are 14.4%, 22.5%, and 4.0% respectively. Thus, the Cmmm phase is much stiffer than low-pressure C2/m and Na2N2-type phases under compression. This is consistent with the increasing coordination number of the K atom in the different K2N2 phases under compression in Figs. 3(c) and 3(d). The densification of K2N2 under pressure is realized by increasing the cation coordination number from 5 in ground-state C2/m phase to 8 in Na2N2-type (K2 atoms are coordinated by eight N atoms and K1 atoms are fourfold coordinated by four N atoms) and Cmmm phases, and finally to 10 in C2/c phase.
2.3 Electronic Structures of K2N2 Phases
To get deeper insight into the electronic structure evolutions of these K2N2 phases under pressure, we calculated the total density of states (DOS) and the site projected DOS, as shown in Fig. 4, where the Fermi level is the vertical dotted line. In contrast with the semiconducting nature of C2/c phase (Fig. 4(d)), the C2/m, Na2N2-type, and Cmmm phases are all metallic with the N-p electron crossing the Fermi level. Thus, the K2N2 possesses a metal-semiconductor transition, that is similar to that of Li-N system in recent work[28]. From Fig. 4(a)–Fig. 4(c), the C2/m, Na2N2-type, and Cmmm phases exhibit similar projected DOS profiles in the valence band region, the K-s and K-p atomic orbitals are located at the deep energy levels (around −32 and −16 eV, respectively), and the N-s atomic orbitals mainly lie around −22 − −10 eV. Moreover, the interaction between K and N atoms is extremely weak due to the absence of overlapping regions in their atomic projected DOSs. From the inspection of the DOS profiles, a clear broadening of different atomic decomposed states is seen in the high-pressure C2/c phase (Fig. 4(d)), leading to clear hybridizations between K-p and N-p atomic orbitals. In contrast to the first three low-pressure K2N2 phases, this distinction is partially due to the increasing coordination number of K atoms in C2/c phase and to the resulting chemical bonding change that will be illustrated in the following Bader charge transfer analysis.
Fig. 5 presents the PBE-resulted energy band structures of the four K2N2 phases, the projected weights of K-s, K-p, and N-s orbital electrons are not taken into account herein because of their exceptionally tiny contributions around the Fermi level in Fig. 4. Compared to the strong metallic nature of the Cmmm phase, which is mainly dominated by N-px/N-pz at the considered energy range of −2 − 2 eV (Fig. 5(c)), the C2/m ( Na2N2-type) phase exhibits a relative weaker metallic behavior that two bands of N-py/N-pz (N-px/N-py) crossing the Fermi level in Fig. 5(a) (Fig. 5(b)). The indirect semiconducting nature of C2/c phase (band gap of 2.0 eV) is demonstrated in Fig. 5(d), with the valence band maximum (VBM) and conduction band minimum (CBM) is located at A and E points, respectively. It is further found that the VBM and CBM of the C2/c phase are principally composed of N-px and N-pz orbitals, respectively.
2.4 Atomic Chemical Bonding of K2N2 Phases
Fig. 6 shows the projected two-dimensional (2D) electron localization functions (ELF) of the four K2N2 phases to further reveal the atomic chemical bonding. The ELF is used to characterize the localized distribution of electrons of different atoms. The ELF values of 1, 0.5, and close to 0 indicate strong electron localization, electron gas, and nonelectron localization, respectively. From Fig. 6(a)–Fig. 6(c), the high-ELF regions are distributed between two N atoms, indicating the strong covalent N-N bonding in the C2/m, Na2N2-type, and Cmmm phases. Lone-pair electrons of K remain from the C2/m to Cmmm phases under pressure, as illustrated by the ELF value of 0 around K atoms, which is in conformity with the projected atomic orbital DOSs described in Fig. 4(a)–Fig. 4(c). For the compressed C2/c phase, some lone-pair electrons of K still remain, yet some of them are transferred to the regions of the K-N bond. Hence, the pressure compels the K-p lone-pair electrons to participate in the chemical bonding with N-p electrons under compression, realizing the high pressure C2/c phase with increased coordination number of K atoms. The electron localization on N-N bond in the C2/c phase (Fig. 6(d)) is weaker than that in the other three low-pressure K2N2 phases, confirming that the N―N single bond of the polymerized one-dimensional nitrogen chains is weaker than the N=N double bond of the N2 quasimolecules.
Finally, the Bader charges in the four phases are listed in Table 2 by using the Bader AIM method. For the ground-state C2/m phase, each N atom accepts 0.491e which corresponds to the electron loss per K atom. When phase transits into the Na2N2-type and Cmmm phases, a relatively stronger ionic bonding feature of K-N bond is disclosed because that the charge transfer is raised to approximately 0.645e − 0.690e as a result of the electron gain of N. Contrarily, it is noted that the atomic charge transfer from the K atom to N atom decrease to 0.623e in the compressed C2/c phase. This variation may be originated from that partial K atoms participate in covalent bonding with the N atom under compression, as indicated by the orbital hybridizations of K-p and N-p in Fig. 4(d). It then can be concluded that the chemical bonding in the C2/c phase containing one-dimensional polymerized nitrogen chains is a complex mixture of covalent and ionic characteristics.
Table 2. Calculated Bader charges of K and N atoms in C2/m, Na2N2-type, Cmmm, and C2/c phasesPhase Pressure/GPa Atom Charge value/e Charge transfer/e C2/m 0 K
N8.509 (×2)
5.491 (×2)+0.491
−0.491Na2N2-type 2.5 K1
K2
N8.407
8.303
5.645 (×2)+0.593
+0.697
−0.645Cmmm 20.0 K
N8.310 (×2)
5.690 (×2)+0.690
−0.690C2/c 135.0 K
N8.377 (×2)
5.623 (×2)+0.623
−0.6233. Conclusions
In summary, we have performed a systematic crystal structure search for alkali metal diazenide K2N2 under pressure using the advanced swarm intelligence structure simulations. The ground-state phase of K2N2 is determined to be a new monoclinic C2/m structure with N2 quasimolecules, and three high-pressure structures of K2N2, including Na2N2-type (1.7−3.6 GPa), Cmmm (3.6−122.0 GPa), and C2/c (above 122.0 GPa), are identified. Remarkably, the pressure-induced polymerization of one-dimensional nitrogen chains in the high-pressure C2/c semiconducting phase is firstly reported, attributed to the participation of K-p lone-pair electrons in the chemical bonding under pressure. The occurrences of the high-pressure phases of Na2N2-type, Cmmm, and C2/c follow the increased coordination numbers (from 5 to 10) of K atoms with the incresed pressures, accompanied by significant N-N bonding modification from N=N dimers to polymerized N―N single bond one-dimensional chains. Furthermore, a metal-semiconductor transition of K2N2 under pressure is demonstrated by the electronic structures calculations. The present findings provide fundamental insights into the high-pressure structures and the polymerization of N atoms of alkali metal diazenides, supplying useful guidance for the syntheses of K2N2 in further experimental efforts.
-
Table 1. Optimized structural parameters of the C2/m, Na2N2-type, Cmmm, and C2/c phases of K2N2
Phase Pressure/GPa Lattice parameters dN-N/Å Atomic fractional coordinates C2/m 0 a = 7.562 Å, b = 3.912 Å, c = 10.923 Å,
α = γ = 90°, β= 134.955°1.192 K 4i (0.613, 0, 0.260)
N 4i (0.519, 0, 0.456)Na2N2-type 2.5 a = 3.326 Å, b = 4.375 Å, c = 5.694 Å,
α = β = γ = 90°1.220 K1 1e (0, 0.500, 0)
K2 1c (0, 0, 0.500)
N 2t (0.500, 0.500, 0.393)Cmmm 20.0 a = 6.896 Å, b = 5.104 Å, c = 2.855 Å,
α = β = γ = 90°1.254 K 4g (0.694, 0, 0)
N 4j (0, 0.877, 0.500)C2/c 135.0 a = 6.971 Å, b = 4.099 Å, c = 4.510 Å,
α = γ = 90°, β= 81.342°1.474, 1.373 K 8f (0.669, 0.909, 0.198)
N 8f (0.041, 0.887, 0.102)Table 2. Calculated Bader charges of K and N atoms in C2/m, Na2N2-type, Cmmm, and C2/c phases
Phase Pressure/GPa Atom Charge value/e Charge transfer/e C2/m 0 K
N8.509 (×2)
5.491 (×2)+0.491
−0.491Na2N2-type 2.5 K1
K2
N8.407
8.303
5.645 (×2)+0.593
+0.697
−0.645Cmmm 20.0 K
N8.310 (×2)
5.690 (×2)+0.690
−0.690C2/c 135.0 K
N8.377 (×2)
5.623 (×2)+0.623
−0.623 -
[1] PIERSON H O. Handbook of refractory carbides and nitrides: properties, characteristics, processing and applications [M]. Westwood: Noyes, 1996. [2] HORVATH-BORDON E, RIEDEL R, ZERR A, et al. High-pressure chemistry of nitride-based materials [J]. Chemical Society Reviews, 2006, 35(10): 987–1014. doi: 10.1039/b517778m [3] CHRISTE K O. Polynitrogen chemistry enters the ring: a cyclo- N−5 anion has been synthesized as a stable salt and characterized [J]. Science, 2017, 355(6323): 351. doi: 10.1126/science.aal5057[4] LANIEL D, WECK G, LOUBEYRE P. Direct reaction of nitrogen and lithium up to 75 GPa: synthesis of the Li3N, LiN, LiN2, and LiN5 compounds [J]. Inorganic Chemistry, 2018, 57(17): 10685–10693. doi: 10.1021/acs.inorgchem.8b01325 [5] YAO Y S, ADENIYI A O. Solid nitrogen and nitrogen-rich compounds as high-energy-density materials [J]. Physica Status Solidi (B), 2021, 258(6): 2000588. doi: 10.1002/pssb.202000588 [6] EREMETS M I, GAVRILIUK A G, TROJAN I A, et al. Single-bonded cubic form of nitrogen [J]. Nature Materials, 2004, 3(8): 558–563. doi: 10.1038/nmat1146 [7] EREMETS M I, GAVRILIUK A G, TROJAN I A. Single-crystalline polymeric nitrogen [J]. Applied Physics Letters, 2007, 90(17): 171904. doi: 10.1063/1.2731679 [8] GREGORYANZ E, GONCHAROV A F, SANLOUP C, et al. High P-T transformations of nitrogen to 170 GPa [J]. The Journal of Chemical Physics, 2007, 126(18): 184505. doi: 10.1063/1.2723069 [9] ZHANG L J, WANG Y C, LV J, et al. Materials discovery at high pressures [J]. Nature Reviews Materials, 2017, 2(4): 17005. doi: 10.1038/natrevmats.2017.5 [10] MEDVEDEV S A, TROJAN I A, EREMETS M I, et al. Phase stability of lithium azide at pressures up to 60 GPa [J]. Journal of Physics: Condensed Matter, 2009, 21(19): 195404. doi: 10.1088/0953-8984/21/19/195404 [11] EREMETS M I, POPOV M Y, TROJAN I A, et al. Polymerization of nitrogen in sodium azide [J]. The Journal of Chemical Physics, 2004, 120(22): 10618–10623. doi: 10.1063/1.1718250 [12] ZHU H Y, ZHANG F X, JI C, et al. Pressure-induced series of phase transitions in sodium azide [J]. Journal of Applied Physics, 2013, 113(3): 033511. doi: 10.1063/1.4776235 [13] ZHOU M, LIU S J, DU M R, et al. High-pressure-induced structural and chemical transformations in NaN3 [J]. The Journal of Physical Chemistry C, 2020, 124(37): 19904–19910. doi: 10.1021/acs.jpcc.0c04107 [14] JI C, ZHANG F X, HOU D B, et al. High pressure X-ray diffraction study of potassium azide [J]. Journal of Physics and Chemistry of Solids, 2011, 72(6): 736–739. doi: 10.1016/j.jpcs.2011.03.005 [15] JI C, ZHENG R, HOU D B, et al. Pressure-induced phase transition in potassium azide up to 55 GPa [J]. Journal of Applied Physics, 2012, 111(11): 112613. doi: 10.1063/1.4726212 [16] WANG Y, BYKOV M, CHEPKASOV I, et al. Stabilization of hexazine rings in potassium polynitride at high pressure [J]. Nature Chemistry, 2022, 14(7): 794–800. doi: 10.1038/s41557-022-00925-0 [17] LI D M, WU X X, JIANG J R, et al. Pressure-induced phase transitions in rubidium azide: studied by in-situ X-ray diffraction [J]. Applied Physics Letters, 2014, 105(7): 071903. doi: 10.1063/1.4893464 [18] LI D M, LI F F, LI Y, et al. High-pressure studies of rubidium azide by Raman and infrared spectroscopies [J]. The Journal of Physical Chemistry C, 2015, 119(29): 16870–16878. doi: 10.1021/acs.jpcc.5b05208 [19] HOU D B, ZHANG F X, JI C, et al. Series of phase transitions in cesium azide under high pressure studied by in situ X-ray diffraction [J]. Physical Review B, 2011, 84(6): 064127. doi: 10.1103/PhysRevB.84.064127 [20] LI D M, ZHU P F, JIANG J R, et al. High-pressure Raman and infrared spectroscopic studies of cesium azide [J]. The Journal of Physical Chemistry C, 2016, 120(47): 27013–27018. doi: 10.1021/acs.jpcc.6b09811 [21] SUI M H, LIU S, WANG P, et al. High-pressure synthesis of fully sp2-hybridized polymeric nitrogen layer in potassium supernitride [J]. Science Bulletin, 2023, 68(14): 1505–1513. doi: 10.1016/j.scib.2023.06.029 [22] ZHANG M G, YAN H Y, WEI Q, et al. Novel high-pressure phase with pseudo-benzene “N6” molecule of LiN3 [J]. Europhysics Letters, 2013, 101(2): 26004. doi: 10.1209/0295-5075/101/26004 [23] WANG X L, LI J F, BOTANA J, et al. Polymerization of nitrogen in lithium azide [J]. The Journal of Chemical Physics, 2013, 139(16): 164710. doi: 10.1063/1.4826636 [24] ZHANG M G, YIN K T, ZHANG X X, et al. Structural and electronic properties of sodium azide at high pressure: a first principles study [J]. Solid State Communications, 2013, 161: 13–18. doi: 10.1016/j.ssc.2013.01.032 [25] ZHANG M G, YAN H Y, WEI Q, et al. A new high-pressure polymeric nitrogen phase in potassium azide [J]. RSC Advances, 2015, 5(16): 11825–11830. doi: 10.1039/C4RA15699D [26] ZHANG X W, ZUNGER A, TRIMARCHI G. Structure prediction and targeted synthesis: a new Na nN2 diazenide crystalline structure [J]. The Journal of Chemical Physics, 2010, 133(19): 194504. doi: 10.1063/1.3488440 [27] SCHNEIDER S B, FRANKOVSKY R, SCHNICK W. High-pressure synthesis and characterization of the alkali diazenide Li2N2 [J]. Angewandte Chemie International Edition, 2012, 51(8): 1873–1875. doi: 10.1002/anie.201108252 [28] SHEN Y Q, OGANOV A R, QIAN G R, et al. Novel lithium-nitrogen compounds at ambient and high pressures [J]. Scientific Reports, 2015, 5: 14204. doi: 10.1038/srep14204 [29] ZHANG J, WANG X L, YANG K S, et al. The polymerization of nitrogen in Li2N2 at high pressures [J]. Scientific Reports, 2018, 8(1): 13144. doi: 10.1038/s41598-018-31355-z [30] WANG Y C, LV J, ZHU L, et al. Crystal structure prediction via particle-swarm optimization [J]. Physical Review B, 2010, 82(9): 094116. doi: 10.1103/PhysRevB.82.094116 [31] WANG Y C, LV J, ZHU L, et al. CALYPSO: a method for crystal structure prediction [J]. Computer Physics Communications, 2012, 183(10): 2063–2070. doi: 10.1016/j.cpc.2012.05.008 [32] LI Y W, FENG X L, LIU H Y, et al. Route to high energy density polymeric nitrogen t-N via He-N compounds [J]. Nature Communications, 2018, 9(1): 722. doi: 10.1038/s41467-018-03200-4 [33] MA L, WANG K, XIE Y, et al. High-temperature superconducting phase in clathrate calcium hydride CaH6 up to 215 K at a pressure of 172 GPa [J]. Physical Review Letters, 2022, 128(16): 167001. doi: 10.1103/PhysRevLett.128.167001 [34] DUAN Q Z, SHEN J Y, ZHONG X, et al. Structural phase transition and superconductivity of ytterbium under high pressure [J]. Physical Review B, 2022, 105(21): 214503. doi: 10.1103/PhysRevB.105.214503 [35] SUN W G, CHEN B L, LI X F, et al. Ternary Na-P-H superconductor under high pressure [J]. Physical Review B, 2023, 107(21): 214511. doi: 10.1103/PhysRevB.107.214511 [36] KRESSE G, FURTHMÜLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set [J]. Physical Review B, 1996, 54(16): 11169–11186. doi: 10.1103/PhysRevB.54.11169 [37] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865 [38] BLÖCHL P E. Projector augmented-wave method [J]. Physical Review B, 1994, 50(24): 17953–17979. doi: 10.1103/PhysRevB.50.17953 [39] GRIMME S, ANTONY J, EHRLICH S, et al. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu [J]. The Journal of Chemical Physics, 2010, 132(15): 154104. doi: 10.1063/1.3382344 [40] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188 [41] TOGO A, TANAKA I. First principles phonon calculations in materials science [J]. Scripta Materialia, 2015, 108: 1–5. doi: 10.1016/j.scriptamat.2015.07.021 [42] GATTI C. Chemical bonding in crystals: new directions [J]. Zeitschrift für Kristallographie-Crystalline Materials, 2005, 220(5/6): 399–457. [43] PENG F, HAN Y X, LIU H Y, et al. Exotic stable cesium polynitrides at high pressure [J]. Scientific reports, 2015, 5(1): 16902. doi: 10.1038/srep16902 [44] BRILL T B, JAMES K J. Kinetics and mechanisms of thermal decomposition of nitroaromatic explosives [J]. Chemical Reviews, 1993, 93(8): 2667–2692. doi: 10.1021/cr00024a005 [45] ZHANG S T, ZHAO Z Y, LIU L L, et al. Pressure-induced stable BeN4 as a high-energy density material [J]. Journal of Power Sources, 2017, 365: 155–161. doi: 10.1016/j.jpowsour.2017.08.086 [46] ZHAI H, XU R, DAI J H, et al. Stabilized nitrogen framework anions in the Ga-N system [J]. Journal of the American Chemical Society, 2022, 144(47): 21640–21647. doi: 10.1021/jacs.2c09056 -