Exchange-correlation effects in interatomic energies for pure density functionals and their application to the molecular energy prediction

In this proof-of-concept paper, we show how exchange-correlation effects can be simply recovered for interatomic energies within the interacting quantum atoms decomposition when local, gradient generalized, or meta-gradient generalized approximations are used in density functional theory (DFT) calculations. We also demonstrate how inhomogeneity and non-local effects can be introduced even from a pure local scheme, without resorting to any orbital information. Finally, we provide numerical evidence on a database of selected energetic molecules that this decomposition scheme can be efficiently used to build accurate models for the prediction of molecular energies from an initial "cheap" DFT calculation.

Computational insight into the crystal structures of cubane and azacubanes

CONTEXT

Using quantum chemistry and atom-atom potential methods, the molecular and crystal structures of cubane 1 and all types of unsubstituted azacubanes 2-22 were calculated. Alternative possible polymorphs of cubane 1 have been proposed. The thermochemical properties of azacubanes in the gas and solid phases were assessed. Thermodynamic aspects of stability are considered, and a significant decrease in stability is revealed upon transition from cubane 1 to octaazacubane 22. It has been shown that the density and energetic properties of azacubanes depend nonlinearly on the number of nitrogen atoms in the structure and the density of octaazacubane 22 at room temperature is 1.546 g cm-3, which is significantly lower than the previously given estimate.

METHODS

In this work, DFT calculations were conducted through the software Gaussian 09 using B3LYP functional with basis set aug-cc-PVDZ and the Grimme dispersion correction D2. For crystal structure optimization, the atom-atom potential methods with PMC (packing of molecules in crystal) program were used. Charges for molecular electrostatic potential were fitted by FitMEP, and enthalpies of formation in gas phase were assessed by G3B3.

A polymeric nitrogen N[Formula: see text]-N[Formula: see text] system with enhanced stability at low pressure

Postulated in 1992 and synthesized in 2004 above 2000 K and 110 GPa, the singly-bonded nitrogen cubic gauche crystal (cg-PN) is still considered to be the ultimate high energy density material (HEDM). The search however has continued for a method to synthesize cg-PN at more ambient conditions or find HEDMs which can be synthesized at lower pressure and temperature. Here, using ab initio evolutionary crystal prediction techniques, a simpler nitrogen-based molecular crystal consisting of N[Formula: see text] and N[Formula: see text] molecules is revealed to be a more favorable polynitrogen at lower pressures. The energetic gain of  534 meV/atom over cg-PN and  138 meV/atom over the N[Formula: see text] molecular crystal at zero pressure makes the N[Formula: see text]-N[Formula: see text] system more appealing. Dynamical and mechanical stabilities are investigated at 5 and 0 GPa, and vibrational frequencies are assessed for its Raman and IR spectra. The prospects of an experimental synthesis of the N[Formula: see text]-N[Formula: see text] polymeric system compared to cg-PN is higher because the C[Formula: see text] symmetry of N[Formula: see text] within this crystal would be easier to target from the readily available N[Formula: see text] azides and the observed N[Formula: see text] and N[Formula: see text] radicals.

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.

Molecular and Electronic Structures of Neutral Polynitrogens: Review on the Theory and Experiment in 21st Century

The data on the existence and physicochemical characteristics of uncharged single element chemical compounds formed by nitrogen atoms and containing more than two nuclides of this element (from N4 to N120, oligomeric and polymeric polynitrogens) have been systematized and generalized. It has been noticed that these data have a predominantly predictive character and were obtained mainly using quantum chemical calculations of various levels (HF, DFT, MP, CCSD etc.). The possibility of the practical application of these single element compounds has been considered. The review mainly covers articles published in the last 25 years. The bibliography contains 128 references.

Insight into energetic performance from structures: A density functional theory study on N6

In this study, in order to compare the energetic performance of poly-nitrogen compounds with different structures, N6 was taken as an example, and three bi-rings structures of N6 were presented...

First-principles study of the structural phase transition process of solid nitrogen under pressure

Due to the diversity of solid nitrogen structure, its phase transition has been a hot topic for many scientists. Herein, we first studied the structural softening of rhombohedral solid nitrogen under pressure using first-principles calculations. Then, a new criterion, Egret criterion, was proposed to predict the whole process from beginning to end of structural phase transition of solid nitrogen. Based on the discussion of acoustic phonons, we concluded that the phase transition of rhombohedral solid nitrogen starts from k-point F along the [- 1, - 1, 0] direction in a-axis, and the structural phase transition velocity is slow. Also, we use the Egret criterion proposed by us to predict the emergence of ξ-N₂ and the stability of ξ-N₂ at 17 GPa and 22 GPa, respectively, and this result is in good agreement with the phase diagram of nitrogen.

All-Nitrogen Cages and Molecular Crystals: Topological Rules, Stability, and Pyrolysis Paths

We combined ab initio molecular dynamics with the intrinsic reaction coordinate in order to investigate the mechanisms of stability and pyrolysis of N4 ÷ N120 fullerene-like nitrogen cages. The stability of the cages was evaluated in terms of the activation barriers and the activation Gibbs energies of their thermal-induced breaking. We found that binding energies, bond lengths, and quantum-mechanical descriptors failed to predict the stability of the cages. However, we derived a simple topological rule that adjacent hexagons on the cage surface resulted in its instability. For this reason, the number of stable nitrogen cages is significantly restricted in comparison with their carbon counterparts. As a rule, smaller clusters are more stable, whereas the earlier proposed large cages collapse at room temperature. The most stable all-nitrogen cages are the N4 and N6 clusters, which can form the van der Waals crystals with densities of 1.23 and 1.36 g/cm3, respectively. The examination of their band structures and densities of electronic states shows that they are both insulators. Their power and sensitivity are not inferior to the modern advanced high-energy nanosystems.

Crystalline Structures and Energetic Properties of Lithium Pentazolate under Ambient Conditions

Recently, it has been reported that high-pressure synthesized lithium pentazolates could be quenched down to ambient conditions. However, the crystalline structures of LiN5 under ambient conditions are still ambiguous. In this work, the structures of LiN5 compound were directly explored at atmospheric pressure by using a new constrain structure search method. By using this method, three new allotropes were confirmed, and they show lower energy than the previous reported LiN5 phases. Both their thermodynamic and dynamic stability were confirmed through formation enthalpies, phonon spectrum, and ab initio molecular dynamics simulations under ambient conditions. Moreover, these three allotropes show similar formation enthalpies and properties, which suggests that it is hard to obtain a single LiN5 phase, which is well consistent with the experimental phenomenon. Furthermore, because of their low formation energy, all of them possess low energy density when they directly decompose to Li3N and nitrogen (0.52 kJ/g). Instead, the decomposed energy could be further improved to 3.78 kJ/g when they decompose under an oxygen-rich environment.

Bipentazole (N10): A Low-Energy Molecular Nitrogen Allotrope with High Intrinsic Stability

In this Letter, we report a crystal structure prediction and characterization of a molecular nitrogen allotrope N₁₀ (bipentazole) using state-of-the-art computational methods. To date, in the form of a P2₁ space group crystal, this allotrope is the most stable predicted form of nitrogen, other than N₂, in the pressure range 0-42 GPa. Its metastability at ambient conditions was justified using phonon dispersion and mechanical properties calculations as well as ab initio molecular dynamics simulations. Due to a high intrinsic stability caused by aromaticity, bipentazole may appear to be the first nitrogen allotrope stable enough for a large-scale synthesis at ambient conditions. The calculations of propulsive characteristics revealed that bipentazole is an excellent "green" energetic material. A potential strategy for the synthesis of this compound is offered and rationalized. The unique electronic structure of bipentazole makes it a strongly electrophilic all-nitrogen reagent, which can exhibit unusual chemistry.

Beyond molecular nitrogen: revelation of two ambient-pressure metastable single- and double-bonded nitrogen allotropes built from three-membered rings

In this study, we present a crystal structure prediction and characterization for two novel single- and double-bonded nitrogen allotropes (denoted as CubN and DobN) bearing three-membered nitrogen cycles. These materials are ambient-pressure metastable robust solids with strong cohesive energies. Due to the presence of strained nitrogen cycles, these allotropes are high-lying phases, which retain high energy densities at least up to 200 GPa, when CubN becomes lower in energy than ζ-N2. The studied allotropes are indirect wide-bandgap semiconductors. A detailed spectral characterization of these materials is also presented herein. Due to their extremely high heats of formation and crystal densities, CubN and DobN demonstrate excellent detonation properties that are comparable to those of previously reported phases (cg-N and TrigN). Propulsive properties (including specific impulse and characteristic velocity) of CubN and DobN as solid monopropellants were also estimated. Thus, if synthesized, the present nitrogen allotropes will have great potential for practical applications and achieving fundamental knowledge about nitrogen as a chemical element.

A unified model of impact sensitivity of metal azides

A comprehensive theoretical study of 18 metal azides is reported.

Spontaneous disproportionation of lithium biphenyl in solution: a combined experimental and theoretical study

Spontaneous disproportionation of lithium biphenyl radical anion into a lithium biphenyl dianion plus neutral biphenyl, according to 2LiBiph ⇄ Li₂Biph + Biph, has been evidenced experimentally by UV-vis spectroscopy and backed up theoretically using TDDFT methods.

Two isomeric solid carbon nitrides with 1 : 1 stoichiometry which exhibit strong mechanical anisotropy

Two isomeric layered carbon nitride polymorphs are characterized using reliable theoretical methods. The NCNC phase, which is predicted for the first time, has a number of differences with the isomeric NCCN polymorph in its electronic, spectral and mechanical properties.