A 2016 computational study proposed fifty phosphorus-doped nitrogen clusters NnPm\mathrm{N}_n\mathrm{P}_m (n+m=4,8,10,14,16,20n+m = 4, 8, 10, 14, 16, 20) as candidate high-energy-density materials, based solely on density-functional-theory formation and binding energies computed for mostly single, hand-picked isomers per stoichiometry, and concluded that the mixed clusters are more stable than the pure Nn+m\mathrm{N}_{n+m}/Pn+m\mathrm{P}_{n+m} clusters of the same size. That conclusion is a positive claim, not a posed question; whether it holds — whether any of these compositions is a practicable energetic material — remains open, and settling it requires: (i) establishing whether the reported structures are the true global minima of their stoichiometries, and whether the claimed thermodynamic advantage of mixed clusters over pure Nn+m\mathrm{N}_{n+m}/Pn+m\mathrm{P}_{n+m} survives systematic, converged global-minimum sampling; (ii) determining the decomposition pathways, activation barriers, and lifetimes — the kinetic stability — of the proposed clusters, which were never computed; and (iii) determining by what realistic route the clusters could form and persist. No published work has yet addressed these points: adjacent experimental work has reached only six-atom or smaller N–P species, and global-optimization or barrier methodology has so far been applied only to pure phosphorus, pure nitrogen, and boron-doped nitrogen clusters.

Background

High-energy-density materials (HEDMs) store large amounts of chemical energy per unit mass and release it by decomposing into very stable small molecules. Nitrogen-rich compounds are the canonical 'green' HEDM candidates because energy can be stored in N–N single and double bonds and recovered on relaxation into the N≡N triple bond of N2\mathrm{N}_2; the central obstacle is that such compounds are typically only metastable, so both the thermodynamic energy content and the kinetic stability (the height of the barriers protecting the metastable state, and hence the lifetime) must be established before any composition is practicable.

Valadbeigi (2016) computationally screened phosphorus-doped nitrogen clusters NnPm\mathrm{N}_n\mathrm{P}_m: fifty clusters with total atom counts n+m=4,8,10,14,16,20n+m = 4, 8, 10, 14, 16, 20. Geometries were optimized at two Kohn–Sham density-functional theory (DFT) levels — the M06-2X and B3LYP exchange–correlation functionals with triple-zeta Gaussian basis sets 6-311++G(d,p) and 6-311++G(2df,p) — testing spin multiplicities 1, 3, and 5 (the multiplicity counts the orientations of the total electronic spin and distinguishes competing electronic ground states). For each cluster the study reported the enthalpy and Gibbs free energy of formation from atomic N and P, and the binding energy per atom Eb=(nEN+mEP−ET)/(n+m)E_b = (nE_N + mE_P - E_T)/(n+m), positive for a bound cluster — the convention consistent with the values the paper reports — where ETE_T is the cluster's total electronic energy and ENE_N, EPE_P the energies of the free atoms. The computed EbE_b values (roughly 2.3–3.6 eV per atom) led to the claim that many mixed stoichiometries are more cohesive than the pure clusters Nn+m\mathrm{N}_{n+m} and Pn+m\mathrm{P}_{n+m} of the same total size — for instance, the paper's abstract cites Eb(N4P6)=3.48E_b(\mathrm{N}_4\mathrm{P}_6) = 3.48 eV, about 50% larger than Eb(N10)=2.49E_b(\mathrm{N}_{10}) = 2.49 eV and Eb(P10)=2.34E_b(\mathrm{P}_{10}) = 2.34 eV — motivating the HEDM candidacy. At the time, the largest well-characterized molecular N–P species were small: the aromatic diphosphatriazolate anion P2N3−\mathrm{P}_2\mathrm{N}_3^- (isolated as a room-temperature-stable salt in 2015) and unassigned PmNn+\mathrm{P}_m\mathrm{N}_n^+ ion series observed by laser ablation of α\alpha-P3N5\mathrm{P}_3\mathrm{N}_5.

Three limitations define the open problem. First, isomer sampling: an isomer is a distinct bonding arrangement of the same set of atoms, and the study considered mostly one chosen structure per stoichiometry (for example a dodecahedral N20\mathrm{N}_{20} and adamantane-type N10\mathrm{N}_{10}/P10\mathrm{P}_{10}), so the reported minima are not established as global minima (the lowest-energy structure over all isomers), and whether the mixed-superiority claim survives proper sampling remains open. Second, kinetics: no transition states (stationary points with exactly one imaginary vibrational frequency, i.e. saddle points on the potential-energy surface separating minima along a reaction coordinate), activation barriers, or decomposition pathways were computed, so the lifetime of every proposed cluster is unknown. Third, formation: all calculations were gas-phase electronic-structure results, with no modeling of temperature, pressure, condensed-phase environments, or any concrete synthesis scenario. A further reliability caveat is that the two DFT levels were not validated against higher-level wavefunction methods or experiment.

Problem Statement

For the phosphorus-doped nitrogen clusters NnPm\mathrm{N}_n\mathrm{P}_m with both elements present (n,m≥1n, m \geq 1) and n+m=4,8,10,14,16,20n+m = 4, 8, 10, 14, 16, 20 proposed as stable high-energy-density materials by Valadbeigi (2016), establish whether the source's stability and HEDM-candidacy claim holds, with quantified and independently reproducible results on: (1) Structures — the isomer landscape of each stoichiometry, determined by systematic unbiased global-optimization searches (for example basin-hopping or genetic-algorithm sampling — stochastic global-optimization methods that repeatedly perturb and locally re-relax candidate structures — followed by DFT refinement), identifying the global minimum of each NnPm\mathrm{N}_n\mathrm{P}_m and testing whether the claimed thermodynamic advantage of mixed clusters over the pure Nn+m\mathrm{N}_{n+m} and Pn+m\mathrm{P}_{n+m} clusters of the same size (for example the paper's Eb(N4P6)=3.48E_b(\mathrm{N}_4\mathrm{P}_6) = 3.48 eV versus Eb(N10)=2.49E_b(\mathrm{N}_{10}) = 2.49 eV and Eb(P10)=2.34E_b(\mathrm{P}_{10}) = 2.34 eV) survives systematic, converged global-minimum sampling — where the sampling counts as converged when its method, search range, and number of candidate isomers are stated and repeated or extended searches no longer find a lower-energy isomer of that stoichiometry. (2) Kinetics — the decomposition pathways and activation barriers of the proposed clusters, via transition-state searches and rate theory where needed, yielding defensible lifetimes under realistic conditions. (3) Formation — whether realistic formation and stabilization routes exist for those compositions that prove thermodynamically and kinetically viable, such as gas-phase cluster assembly, matrix isolation (trapping the clusters in a frozen inert-gas matrix), laser ablation, or high-pressure synthesis, including how such conditions modify the thermodynamics and kinetics. The single overarching objective is a viability verdict for the proposed NnPm\mathrm{N}_n\mathrm{P}_m HEDM family. A demonstrated refutation of the source's stability claim on any one of (1)–(3) — for example global-minimum searches overturning the reported thermodynamic advantage, or computed decomposition barriers implying short lifetimes — resolves the task and must rest on such demonstrated results, not on the absence of supporting analysis; when (1) or (2) already refutes the claim, dimension (3) need not be completed. A positive viability verdict instead requires quantified results on all three dimensions and must state the temperature, pressure, and lifetime conditions under which each surviving composition is judged viable.

References

  1. Younes Valadbeigi, 'Phosphorus-doped nitrogen clusters (NnPm): Stable high energy density materials', Chemical Physics Letters 645 (2016) 195–199, DOI 10.1016/j.cplett.2016.01.007, https://doi.org/10.1016/j.cplett.2016.01.007
  2. Chaojiang Zhang, Cheng Zhu, André K. Eckhardt, Ralf I. Kaiser, 'Formation of the elusive tetrahedral P3N molecule', Science Advances 8 (2022) eabo5792, DOI 10.1126/sciadv.abo5792, https://doi.org/10.1126/sciadv.abo5792
  3. Shyam Sharan Tripathi, Srikanth Korutla, Rishabh Kumar Pandey, Tammineni Rajagopala Rao, 'Theoretical Elucidation of Photoelectron Spectrum of Diphosphatriazolate Anion, P2N3−', The Journal of Physical Chemistry A 129 (2025) 7617–7629, DOI 10.1021/acs.jpca.5c02662, https://doi.org/10.1021/acs.jpca.5c02662
  4. Wen-Jie Wu, Wei-Jie Chi, Quan-Song Li, Ze-Sheng Li, 'Strategy for designing stable and powerful nitrogen-rich high-energy materials by introducing boron atoms', Journal of Molecular Modeling 23 (2017) 191, DOI 10.1007/s00894-017-3360-6, https://doi.org/10.1007/s00894-017-3360-6
  5. S. Mahtout, N. Amatousse, F. Rabilloud, 'Structural, electronic and magnetic properties of Pn+1 and FePn (n = 1–14) clusters', Computational and Theoretical Chemistry 1122 (2017) 16–26, DOI 10.1016/j.comptc.2017.10.010, https://doi.org/10.1016/j.comptc.2017.10.010
  6. M. A. Roumili, S. Mahtout, K. Baddari, M. Lasmi, F. Rabilloud, 'Phosphorus, cobalt-phosphorus, and nickel-phosphorus clusters: Growth behavior, electronic, and magnetic properties', International Journal of Quantum Chemistry 123 (2023) e27189, DOI 10.1002/qua.27189, https://doi.org/10.1002/qua.27189
  7. Alexandra Velian, Christopher C. Cummins, 'Synthesis and characterization of P2N3−: An aromatic ion composed of phosphorus and nitrogen', Science 348 (2015) 1001–1004, DOI 10.1126/science.aab0204, https://doi.org/10.1126/science.aab0204
  8. Yifan Jin, Ajith Perera, Rodney J. Bartlett, 'Spectroscopic analysis of diphosphatriazolate anion (P2N3−) by coupled-cluster methods as a step toward N5−', Chemical Physics Letters 640 (2015) 68–71, DOI 10.1016/j.cplett.2015.10.014, https://doi.org/10.1016/j.cplett.2015.10.014
  9. Cheng Zhu, André K. Eckhardt, Sankhabrata Chandra, Andrew M. Turner, Peter R. Schreiner, Ralf I. Kaiser, 'Identification of a prismatic P3N3 molecule formed from electron irradiated phosphine-nitrogen ices', Nature Communications 12 (2021) 5467, DOI 10.1038/s41467-021-25775-1, https://doi.org/10.1038/s41467-021-25775-1
  10. Sachinkumar Dagurao Pangavhane, Lucie Hebedová, Milan Alberti, Josef Havel, 'Laser ablation synthesis of new phosphorus nitride clusters from α-P3N5 via laser desorption ionization and matrix assisted laser desorption ionization time-of-flight mass spectrometry', Rapid Communications in Mass Spectrometry 25 (2011) 917–924, DOI 10.1002/rcm.4937, https://doi.org/10.1002/rcm.4937
  11. A. C. R. Gomes, M. X. Silva, B. R. L. Galvão, 'Stability of neutral molecular polynitrogens: energy content and decomposition mechanisms', RSC Advances 11 (2021) 21567–21578, DOI 10.1039/d1ra03259c, https://doi.org/10.1039/d1ra03259c
  12. Alyona A. Starikova, Natalia M. Boldyreva, Ruslan M. Minyaev, Alexander I. Boldyrev, Vladimir I. Minkin, 'Computational Assessment of an Elusive Aromatic N3P3 Molecule', ACS Omega 3 (2018) 286–291, DOI 10.1021/acsomega.7b01719, https://doi.org/10.1021/acsomega.7b01719
  13. Jiale Li, Meicheng Chen, Kaiwen Liu, Kewei Ding, Hongguang Xu, Weijun Zheng, Cheng Lu, 'High-nitrogen-content energetic BNn+ (n = 4–16) clusters', Physical Chemistry Chemical Physics 27 (2025) 3567–3573, DOI 10.1039/d4cp03566f, https://doi.org/10.1039/d4cp03566f
  14. Linwei Sai, Xiaoming Huang, Xiaoqing Liang, Xue Wu, Ruili Shi, Di Wu, 'Structural Evolution of Medium-Sized Phosphorus Clusters (P20–P36) from Ab Initio Global Search', Journal of Cluster Science 31 (2019) 567–574, DOI 10.1007/s10876-019-01754-x, https://doi.org/10.1007/s10876-019-01754-x
  15. Min Zhou, Yicheng Xu, Yongliang Cui, Xianyi Zhang, Xianglei Kong, 'Search for Global Minimum Structures of P2n+1+ (n = 1–15) Using xTB-Based Basin-Hopping Algorithm', Frontiers in Chemistry 9 (2021) 694156, DOI 10.3389/fchem.2021.694156, https://doi.org/10.3389/fchem.2021.694156