Assessing the Performances of Dispersion-Corrected Density Functional Methods for Predicting the Crystallographic Properties of High Nitrogen Energetic Salts

Assessing Exchange-Correlation Functionals for Heterogeneous Catalysis of Nitrogen Species

Increasing interest in the sustainable synthesis of ammonia, nitrates, and urea has led to an increase in studies of catalytic conversion between nitrogen-containing compounds using heterogeneous catalysts. Density functional theory (DFT) is commonly employed to obtain molecular-scale insight into these reactions, but there have been relatively few assessments of the exchange-correlation functionals that are best suited for heterogeneous catalysis of nitrogen compounds. Here, we assess a range of functionals ranging from the generalized gradient approximation (GGA) to the random phase approximation (RPA) for the formation energies of gas-phase nitrogen species, the lattice constants of representative solids from several common classes of catalysts (metals, oxides, and metal-organic frameworks (MOFs)), and the adsorption energies of a range of nitrogen-containing intermediates on these materials. The results reveal that the choice of exchange-correlation functional and van der Waals correction can have a surprisingly large effect and that increasing the level of theory does not always improve the accuracy for nitrogen-containing compounds. This suggests that the selection of functionals should be carefully evaluated on the basis of the specific reaction and material being studied.

Performance of Dispersion-Inclusive Density Functional Theory Methods for Energetic Materials

Molecular crystals of energetic materials (EMs) are denser than typical molecular crystals and are characterized by distinct intermolecular interactions between nitrogen-containing moieties. To assess the performance of dispersion-inclusive density functional theory (DFT) methods, we have compiled a data set of experimental sublimation enthalpies of 31 energetic materials. We evaluate the performance of three methods: the semilocal Perdew-Burke-Ernzerhof (PBE) functional coupled with the pairwise Tkatchenko-Scheffler (TS) dispersion correction, PBE with the many-body dispersion (MBD) method, and the PBE-based hybrid functional (PBE0) with MBD. Zero-point energy contributions and thermal effects are described using the quasi-harmonic approximation (QHA), including explicit treatment of thermal expansion, which we find to be non-negligible for EMs. The lattice energies obtained with PBE0+MBD are the closest to experimental sublimation enthalpies with a mean absolute error of 9.89 kJ/mol. However, the state-of-the-art treatment of vibrational and thermal contributions makes the agreement with experiment worse. Pressure-volume curves are also examined for six representative materials. For pressure-volume curves, all three methods provide reasonable agreement with experimental data with mean absolute relative errors of 3% or less. Most of the intermolecular interactions typical of EMs, namely nitro-amine, nitro-nitro, and nitro-hydrogen interactions, are more sensitive to the choice of the dispersion method than to the choice of the exchange-correlation functional. The exception is π-π stacking interactions, which are also very sensitive to the choice of the functional. Overall, we find that PBE+TS, PBE+MBD, and PBE0+MBD do not perform as well for energetic materials as previously reported for other classes of molecular crystals. This highlights the importance of testing dispersion-inclusive DFT methods for diverse classes of materials and the need for further method development.

High pressure structural behaviour of 5,5'-bitetrazole-1,1'-diolate based energetic materials: a comparative study from first principles calculations

Pressure on the scale of gigapascals can cause incredible variations in the physicochemical and detonation characteristics of energetic materials. As a continuation of our earlier work (B. Moses Abraham, et al., Phys. Chem. Chem. Phys., 2018, 20, 29693-29707), here we report the high pressure structural and vibrational properties of 5,5'-bitetrazole-1,1'-diolate based energetic ionic salts via dispersion-corrected density functional theory calculations. Remarkably, these energetic materials exhibit anisotropic behavior along three crystallographic directions with progressing pressure; especially, the maximum and minimum reduction in volume is observed for HA-BTO and TKX-50, respectively. The large bulk modulus of TKX-50 (28.64) indicates its hard nature when compared to other BTO-based energetic salts. The effect of pressure on hydrogen bonded D-H⋯A energetic materials induces spectral shift (lengthening/shortening) in the donor group (D-H) of the stretching vibrations and is widely recognized as the signature of hydrogen bonding. We observed unusual contraction of the D-H bond under compression due to the short range repulsive forces encountered by the H atom while the molecule attempts to stabilize. The Hirshfeld surface analysis highlights the pressure induced stabilization of HA-BTO due to increased N⋯H/H⋯N and O⋯H/H⋯O close contact of hydrogen bond acceptors and donors. These studies provide theoretical guidance as a function of pressure, on how the micro-structures and intermolecular interactions can tune macroscopic properties to enhance the energetic performance.

A density functional theory based approach for predicting melting points of ionic liquids

Accurate prediction of melting points of ILs is important both from the fundamental point of view and from the practical perspective for screening ILs with low melting points and broadening their utilization in a wider temperature range. In this work, we present an ab initio approach to calculate melting points of ILs with known crystal structures and illustrate its application for a series of 11 ILs containing imidazolium/pyrrolidinium cations and halide/polyatomic fluoro-containing anions. The melting point is determined as a temperature at which the Gibbs free energy of fusion is zero. The Gibbs free energy of fusion can be expressed through the use of the Born-Fajans-Haber cycle via the lattice free energy of forming a solid IL from gaseous phase ions and the sum of the solvation free energies of ions comprising IL. Dispersion-corrected density functional theory (DFT) involving (semi)local (PBE-D3) and hybrid exchange-correlation (HSE06-D3) functionals is applied to estimate the lattice enthalpy, entropy, and free energy. The ions solvation free energies are calculated with the SMD-generic-IL solvation model at the M06-2X/6-31+G(d) level of theory under standard conditions. The melting points of ILs computed with the HSE06-D3 functional are in good agreement with the experimental data, with a mean absolute error of 30.5 K and a mean relative error of 8.5%. The model is capable of accurately reproducing the trends in melting points upon variation of alkyl substituents in organic cations and replacement one anion by another. The results verify that the lattice energies of ILs containing polyatomic fluoro-containing anions can be approximated reasonably well using the volume-based thermodynamic approach. However, there is no correlation of the computed lattice energies with molecular volume for ILs containing halide anions. Moreover, entropies of solid ILs follow two different linear relationships with molecular volume for halides and polyatomic fluoro-containing anions. Continuous progress in predicting crystal structures of organic salts with halide anions will be a key factor for successful prediction of melting points with no prior knowledge of the crystal structure.

First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

Vapor Liquid Equilibria of Hydrofluorocarbons Using Dispersion-Corrected and Nonlocal Density Functionals

Recent developments in dispersion corrected and nonlocal density functionals are aimed at accurately capturing dispersion interactions, a key shortcoming of local and semilocal approximations of density functional theory. These functionals have shown significant promise for dimers and small clusters of molecules as well as crystalline materials. However, their efficacy for predicting vapor liquid equilibria is largely unexplored. In this work, we examine the accuracy of dispersion-corrected and nonlocal van der Waals functionals by computing the vapor liquid coexistence curves (VLCCs) of hydrofluoromethanes. Our results indicate that the PBE-D3 functional performs significantly better in predicting saturated liquid densities than the rVV10 functional. With the PBE-D3 functional, we also find that as the number of fluorine atoms increase in the molecule, the accuracy of saturated liquid density prediction improves as well. All the functionals significantly underpredict the saturated vapor densities, which also result in an underprediction of saturated vapor pressure of all compounds. Despite the differences in the bulk liquid densities, the local microstructures of the liquid CFH3 and CF2H2 are relatively insensitive to the density functional employed. For CF3H, however, rVV10 predicts slightly more structured liquid than the PBE-D3 functional.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Structural, electronic and optical properties of well-known primary explosive: Mercury fulminate

Mercury Fulminate (MF) is one of the well-known primary explosives since 17th century and it has rendered invaluable service over many years. However, the correct molecular and crystal structures are determined recently after 300 years of its discovery. In the present study, we report pressure dependent structural, elastic, electronic and optical properties of MF. Non-local correction methods have been employed to capture the weak van der Waals interactions in layered and molecular energetic MF. Among the non-local correction methods tested, optB88-vdW method works well for the investigated compound. The obtained equilibrium bulk modulus reveals that MF is softer than the well known primary explosives Silver Fulminate (SF), silver azide and lead azide. MF exhibits anisotropic compressibility (b > a > c) under pressure, consequently the corresponding elastic moduli decrease in the following order: C22 > C11 > C33. The structural and mechanical properties suggest that MF is more sensitive to detonate along c-axis (similar to RDX) due to high compressibility of Hg⋯O non-bonded interactions along that axis. Electronic structure and optical properties were calculated including spin-orbit (SO) interactions using full potential linearized augmented plane wave method within recently developed Tran-Blaha modified Becke-Johnson (TB-mBJ) potential. The calculated TB-mBJ electronic structures of SF and MF show that these compounds are indirect bandgap insulators. Also, SO coupling is found to be more pronounced for 4d and 5d-states of Ag and Hg atoms of SF and MF, respectively. Partial density of states and electron charge density maps were used to describe the nature of chemical bonding. Ag-C bond is more directional than Hg-C bond which makes SF to be more unstable than MF. The effect of SO coupling on optical properties has also been studied and found to be significant for both (SF and MF) of the compounds.