First-principles calculations (Chen et al., Matter and Radiation at Extremes 10, 015801) predict that substituting 12.5–25% of the nitrogen atoms in black-phosphorus-structured polymeric nitrogen (BP-N) with phosphorus — uniformly within the puckered layers and with perpendicular interlayer arrangement — makes the phase structurally and dynamically stable at 0 GPa, with a predicted gravimetric energy density of 8.07 kJ/g (about twice TNT). However, this 'realization' is purely computational: it rests on PBE+vdW energetics, zero-temperature phonon analysis, exhaustive sampling only at 12.5% doping, and the empirical Kamlet–Jacobs detonation model. The open problem is whether ambient-pressure-stable P-doped BP-N can actually be synthesized, recovered, and validated: identifying a synthesis route that produces the required dopant uniformity and interlayer registry, establishing stability under finite temperature, disorder, defects, and atmospheric exposure, and confirming the energetics with higher-accuracy methods — culminating in experimental realization or refutation of the prediction.

Background

Nitrogen gas owes its inertness to the N≡N triple bond (about 954 kJ/mol), far stronger than N=N double (418 kJ/mol) or N−N single (160 kJ/mol) bonds. Polymeric nitrogen solids built from single- or double-bonded nitrogen networks are therefore metastable high-energy-density materials (HEDMs): decomposing them back into N₂ releases large amounts of energy while producing only harmless N₂. High-pressure experiments have produced several such phases — cubic gauche nitrogen (cg-N, a three-dimensional network of single-bonded N with a distorted cubic lattice) above 110 GPa (Eremets et al., 2004), hexagonal layered polymeric nitrogen HLP-N near 250 GPa (Laniel et al., 2019), and others — but all of them decompose back toward molecular N₂ or chain-like structures when pressure is released (cg-N at about 42 GPa, HLP-N at about 66 GPa), so none can be recovered to ambient conditions. In 2020, BP-N, a polymeric nitrogen isostructural to black phosphorus (the layered allotrope of phosphorus built from puckered, honeycomb-like N sheets), was synthesized experimentally at roughly 130–150 GPa (Ji et al., Science Advances 2020; companion Liu et al., Chinese Physics B 2020), but it too decomposes on decompression, in its case into chain-like structures. Phosphorus sits directly below nitrogen in group V of the periodic table, and black phosphorus itself is perfectly stable at ambient conditions. This motivated Chen, Zhang, Zhu, Cao, Zhang, and Wang (Matter and Radiation at Extremes 10, 015801; arXiv:2405.06212) to ask whether phosphorus doping could stabilize BP-N at ambient pressure. Their density-functional-theory (DFT) study, using the PBE (Perdew–Burke–Ernzerhof) exchange-correlation functional with van der Waals (vdW) corrections — corrections needed for the weak attraction between the puckered N layers — examined 16-atom supercells at 12.5% and 25% doping. They found that a uniform distribution of P atoms within each layer maintains structural stability at 0 GPa — uniformity is quantified in the source by a Gaussian-based metric over pairwise P–P distances, and configurations with uniformity below 0.80 are unstable because of P clustering — and that when P atoms in adjacent layers sit directly above one another's interlayer positions in a perpendicular arrangement, the resulting N–P electric dipoles lock the layers against interlayer sliding, giving at 12.5% doping the only configuration with no imaginary frequencies in its phonon dispersion (the standard zero-kelvin criterion for dynamic stability); at 25% doping, four of the sampled configurations — all with perpendicular interlayer P arrangement and interlayer N–P dipoles — are free of imaginary frequencies. For the 12.5% configuration they computed a gravimetric energy density of 8.07 kJ/g and volumetric energy density of 24.71 kJ/cm³ — about two (three) times TNT — using the empirical Kamlet–Jacobs detonation equations, a 1960s engineering correlation for detonation velocity and pressure. Three limitations define the open problem. First, synthesis accessibility: no experimental route is known that could incorporate phosphorus into BP-N at gigapascal pressures while achieving the required statistical rarity of configurations — intralayer uniformity and perpendicular interlayer registry — since dopant placement is not directly controllable during high-pressure synthesis, and it is unknown whether any route (high-pressure N–P reaction, controlled diffusion, or other protocols) can produce and quench the doped phase. Second, stability beyond zero kelvin: the analysis used 0 K phonons, exhaustive sampling only at 12.5% (50 random configurations at 25%), and no molecular dynamics, so nothing is known about finite-temperature stability, sensitivity to defects and realistic configurational disorder, mechanical/thermal stability, or reactivity with air and containers. Third, accuracy: PBE+vdW energetics and the empirical Kamlet–Jacobs model carry unquantified errors for a novel doped polynitride; higher-level energetics (hybrid functionals, the random-phase approximation (RPA), or quantum Monte Carlo) and, ultimately, experimental measurement are needed.

Problem Statement

Establish whether phosphorus-doped BP-N — black-phosphorus-structured polymeric nitrogen with roughly 12.5–25% of the N atoms replaced by P, uniformly distributed within the puckered layers and arranged perpendicularly between adjacent layers — can be physically realized and remain stable at ambient pressure and temperature, as predicted by PBE+vdW density-functional calculations (Chen et al., Matter Radiat. Extremes 10, 015801). Stability here is meant in the operational sense of the source calculations: structural stability — the doped BP-N lattice is maintained under full relaxation at 0 GPa without decomposing into molecular N₂ or chain-like nitrogen — and dynamic stability — no imaginary frequencies in the phonon dispersion. Because doped BP-N is, like all polymeric nitrogens, metastable by construction rather than a ground state (the source computes positive decomposition enthalpies of 1.34 eV/atom at 12.5% and 1.18 eV/atom at 25% relative to α′-P₃N₅ + N₂), stability of a recovered sample at ambient pressure and temperature means persistence as the doped BP-N phase rather than decomposition. The single overarching objective is validation of the predicted ambient stability of doped BP-N as a practical high-energy-density material; its three requested components jointly serve that objective: (i) synthesis — identify an experimental route (high-pressure synthesis from N and P precursors, controlled diffusion, or another protocol) that produces the doped BP-N phase with the required intralayer dopant uniformity and interlayer registry, and determine whether it can be recovered (quenched) to ambient conditions; (ii) stability — determine whether the predicted stability survives finite temperature, realistic substitutional disorder and defects (beyond the complete enumeration at 12.5% and 50-configuration random sampling at 25%), mechanical and thermal loading, and exposure to atmosphere or containers, using for example ab initio molecular dynamics and explicit configurational sampling; and (iii) energetics — determine the decomposition enthalpy and energy density with higher-accuracy methods than PBE+vdW and beyond the empirical Kamlet–Jacobs detonation model, and validate them against measurement. Completion boundary: the task is resolved either affirmatively — by proposing a synthesis route that could produce the intralayer-uniform, interlayer-perpendicular doped phase and could be recovered (quenched) to ambient conditions, together with computational evidence that the predicted stability and energy density survive beyond the source's zero-kelvin, single-configuration PBE+vdW treatment — or negatively — by a decisive demonstration that the doped phase cannot be so realized or is not stable at ambient pressure and temperature. The three components (i)–(iii) are the requested joint scope of the investigation rather than independent pass/fail gates, and the conditions listed under (ii) define the class of effects to be examined; a decisive refutation established through any of them resolves the overarching question. An affirmative resolution may be analytical rather than experimental — a proposed synthesis protocol supported by such computational evidence, not a demonstrated synthesis or a measured energy release, suffices — and component (iii)'s comparison with measurement is a longer-term goal rather than a condition for resolution.

References

  1. Guo Chen, Chengfeng Zhang, Yuanqin Zhu, Bingqing Cao, Jie Zhang, Xianlong Wang, Realized stable BP-N at ambient pressure by phosphorus doping, Matter and Radiation at Extremes 10, 015801 (2025), DOI: 10.1063/5.0239841, arXiv:2405.06212, https://arxiv.org/abs/2405.06212
  2. Cheng Ji, Adebayo A. Adeleke, Liuxiang Yang, Biao Wan, Huiyang Gou, Yansun Yao, Bing Li, Yue Meng, Jesse S. Smith, Vitali B. Prakapenka, Wenjun Liu, Guoyin Shen, Wendy L. Mao, Ho-kwang Mao, Nitrogen in black phosphorus structure, Science Advances 6, eaba9206 (2020), DOI: 10.1126/sciadv.aba9206, https://doi.org/10.1126/sciadv.aba9206
  3. Ying Liu, Haipeng Su, Caoping Niu, Xianlong Wang, Junran Zhang, Zhongxue Ge, Yanchun Li, Synthesis of black phosphorus structured polymeric nitrogen, Chinese Physics B 29, 106201 (2020), DOI: 10.1088/1674-1056/aba9bd, https://doi.org/10.1088/1674-1056/aba9bd
  4. Mikhail I. Eremets, Alexander G. Gavriliuk, Ivan A. Trojan, Dymitro A. Dzivenko, Reinhard Boehler, Single-bonded cubic form of nitrogen, Nature Materials 3, 558–563 (2004), DOI: 10.1038/nmat1146, https://doi.org/10.1038/nmat1146
  5. D. Laniel, G. Geneste, G. Weck, M. Mezouar, P. Loubeyre, Hexagonal Layered Polymeric Nitrogen Phase Synthesized near 250 GPa, Physical Review Letters 122, 066001 (2019), DOI: 10.1103/PhysRevLett.122.066001, https://doi.org/10.1103/PhysRevLett.122.066001
  6. Jingyi Liu, Binbin Wu, Leilei Zhang, Xue Chang, Xinyu Zhao, Chunmei Fan, Meifang Pu, Chunyin Zhou, Xiaoli Wang, Li Lei, Observation of the Phonon Splitting Associated with Strong Interlayered Interaction in Black Phosphorus Structured Polymeric Nitrogen, Journal of Physical Chemistry C 129, 9755–9763 (2025), DOI: 10.1021/acs.jpcc.5c00788, https://doi.org/10.1021/acs.jpcc.5c00788
  7. Matteo Ceppatelli, Manuel Serrano-Ruiz, Marta Morana, Kamil Dziubek, Demetrio Scelta, Roberto Bini, Maurizio Peruzzini, High-Pressure Synthesis of Pnictogen Nitrides, Accounts of Chemical Research 58, 2927–2938 (2025), DOI: 10.1021/acs.accounts.5c00309, https://doi.org/10.1021/acs.accounts.5c00309
  8. Yanfeng Zhang, Guo Chen, Chengfeng Zhang, Jie Zhang, Xianlong Wang, Exfoliating from high-pressure arsenic–nitrogen compounds: an efficient way to obtain 2D-AsN, Physical Chemistry Chemical Physics 27, 7532–7537 (2025), DOI: 10.1039/d5cp00461f, https://doi.org/10.1039/d5cp00461f
  9. Hui Liang, Di Wang, Rui Xu, Hao Chen, Dan Zhou, Yunwei Zhang, Quan Li, Ultrahigh strength of cage-like polymeric nitrogen surpassing diamond under high pressure, Matter and Radiation at Extremes 11, 017802 (2026), DOI: 10.1063/5.0296196, https://doi.org/10.1063/5.0296196
  10. Ran Liu, Miao Yu, Haodi Zhang, Peng Wang, Shuang Liu, Zhen Yao, Bingbing Liu, Shuangchen Lu, Capturing High-Energy Polymeric Nitrogen Structures Stabilized by Lutetium at Ambient Conditions, Inorganic Chemistry 64, 2748–2757 (2025), DOI: 10.1021/acs.inorgchem.4c04710, https://doi.org/10.1021/acs.inorgchem.4c04710