Prediction of densities for solid energetic molecules with molecular surface electrostatic potentials

The densities of high energetic molecules in the solid state were calculated with a simplified scheme based on molecular surface electrostatic potentials (MSEP). The MSEP scheme for density estimation, originally developed by Politzer et al., was further modified to calculate electrostatic potential on a simpler van der Waals surface. Forty-one energetic molecules containing at least one nitro group were selected from among a variety of molecular types and density values, and were used to test the suitability of the MSEP scheme for predicting the densities of solid energetic molecules. For comparison purposes, we utilized the group additivity method (GAM) incorporating the parameter sets developed by Stine (Stine-81) and by Ammon (Ammon-98 and -00). The absolute average error in densities from our MSEP scheme was 0.039 g/cc. The results based on our MSEP scheme were slightly better than the GAM results. In addition, the errors in densities generated by the MSEP scheme were almost the same for various molecule types, while those predicted by GAM were somewhat dependent upon the molecule types.     

Comparative theoretical studies of energetic azo s-triazines

In this work, the properties of the synthesized high-nitrogen compounds 4,4',6,6'-tetra(azido)azo-1,3,5-triazine (TAAT) and 4,4',6,6'-tetra(azido)hydrazo-1,3,5-triazine (TAHT), and a set of designed bridged triazines with similar bridges were studied theoretically to facilitate further developments for the molecules of interests. The gas-phase heats of formation were predicted based on the isodesmic reactions by using the DFT-B3LYP/AUG-cc-PVDZ method. The estimates of the condensed-phase heats of formation and heats of sublimation were estimated in the framework of the Politzer approach. Calculation results show that the method gives a good estimation for enthalpies, in comparison with available experimental data for TAAT and TAHT. The crystal density has been computed using molecular packing calculations. The calculated detonation velocities and detonation pressures indicate that -NF(2), -NO(2), -N═N-, and -N═N(O)- groups are effective structural units for improving the detonation performance of the bridged triazines. The synthesized TAAT and TAHT are not preferred energetic materials due to their inferior detonation performance. The p→π conjugation effect between the triazine rings and bridges makes the molecule stable as a whole. The electrostatic behavior of the bridged triazines is characterized by an anomalous surface potential imbalance when incorporating the strongly electron-withdrawing -NF(2) and -NO(2) groups into the molecule. An analysis of the bond dissociation energies shows that all these derivatives have good thermal stability over RDX and HMX, and the -NH-NH- bridge is more helpful for improving the stability than -N═N(O)- and -N═N- bridges. Considering the detonation performance and thermal stability, three bridged triazines may be considered as the potential candidates of high-energy density materials (HEDMs).     

Importance of van der Waals Interactions in QM/MM Simulations

The importance of accurately treating van der Waals interactions between the quantum mechanical (QM) and molecular mechanical (MM) atoms in hybrid QM/MM simulations has been investigated systematically. First, a set of van der Waals (vdW) parameters was optimized for an approximate density functional method, the self-consistent charge-tight binding density functional (SCC-DFTB) approach, based on small hydrogen-bonding clusters. The sensitivity of condensed phase observables to the SCC-DFTB vdW parameters was then quantitatively investigated by SCC-DFTB/MM simulations of several model systems using the optimized set and two sets of extreme vdW parameters selected from the CHARMM22 forcefield. The model systems include a model FAD molecule in solution and a solvated enediolate, and the properties studied include the radial distribution functions of water molecules around the solute (model FAD and enediolate), the reduction potential of the model FAD and the potential of mean force for an intramolecular proton transfer in the enediolate. Although there are noticeable differences between parameter sets for gas-phase clusters and solvent structures around the solute, thermodynamic quantities in the condensed phase (e.g., reduction potential and potential of mean force) were found to be less sensitive to the numerical values of vdW parameters. The differences between SCC-DFTB/MM results with the three vdW parameter sets for SCC-DFTB atoms were explained in terms of the effects of the parameter set on solvation. The current study has made it clear that efforts in improving the reliability of QM/MM methods for energetical properties in the condensed phase should focus on components other than van der Waals interactions between QM and MM atoms.     

Molecular design of prismane-based potential energetic materials with high detonation performance and low impact sensitivity

To develop new energetic materials, the eleven nitroester substitution derivatives of prismane were investigated at the B3LYP/6-311G** level of density functional theory (DFT). The gas phase heats of formation were calculated by isodesmic reactions and the solid-state heats of formation were obtained by the Politzer approach using the heats of sublimation for the designed compounds. The detonation velocities and pressures of all molecules were calculated by Kamlet–Jacobs equations based on molecular density and heat of detonation. The results show that the nitroester group in prismane is helpful for enhancing molecular detonation properties and power index. Among all molecules, 1,2,3,4-tetrnitroesterprismane has excellent detonation properties (detonation pressure = 40.05 GPa, detonation velocity = 9.28 km/s) and large power index value. The molecular stabilities were evaluated by calculating bond dissociation energies and characteristic heights (H₅₀). The results indicate that the bond dissociation energies of all molecules are above 80 kJ/mol, and all molecules have a larger H₅₀ value than hexanitrohexaazaisowurtzitane (CL-20, 12 cm). The obtained structure–property relationships may provide basic information for the molecular design of novel high-energy materials.     

Application of Molecular Dynamics Simulations in Molecular Property Prediction I: Density and Heat of Vaporization

Molecular mechanical force field (FF) methods are useful in studying condensed phase properties. They are complementary to experiment and can often go beyond experiment in atomic details. Even a FF is specific for studying structures, dynamics and functions of biomolecules, it is still important for the FF to accurately reproduce the experimental liquid properties of small molecules that represent the chemical moieties of biomolecules. Otherwise, the force field may not describe the structures and energies of macromolecules in aqueous solutions properly. In this work, we have carried out a systematic study to evaluate the General AMBER Force Field (GAFF) in studying densities and heats of vaporization for a large set of organic molecules that covers the most common chemical functional groups. The latest techniques, such as the particle mesh Ewald (PME) for calculating electrostatic energies, and Langevin dynamics for scaling temperatures, have been applied in the molecular dynamics (MD) simulations. For density, the average percent error (APE) of 71 organic compounds is 4.43% when compared to the experimental values. More encouragingly, the APE drops to 3.43% after the exclusion of two outliers and four other compounds for which the experimental densities have been measured with pressures higher than 1.0 atm. For heat of vaporization, several protocols have been investigated and the best one, P4/ntt0, achieves an average unsigned error (AUE) and a root-mean-square error (RMSE) of 0.93 and 1.20 kcal/mol, respectively. How to reduce the prediction errors through proper van der Waals (vdW) parameterization has been discussed. An encouraging finding in vdW parameterization is that both densities and heats of vaporization approach their "ideal" values in a synchronous fashion when vdW parameters are tuned. The following hydration free energy calculation using thermodynamic integration further justifies the vdW refinement. We conclude that simple vdW parameterization can significantly reduce the prediction errors. We believe that GAFF can greatly improve its performance in predicting liquid properties of organic molecules after a systematic vdW parameterization, which will be reported in a separate paper.     

Predicting solid-state heats of formation of newly synthesized polynitrogen materials by using quantum mechanical calculations

We present density functional theory level predictions and analysis of the basic properties of newly synthesized high-nitrogen compounds together with 3,6-bis(2H-tetrazol-5-yl)-1,2,4,5-tetrazine (BTT) and 3,3'-azobis(6-amino-1,2,4,5-tetrazine) (DAAT), for which experimental data are available. The newly synthesized high-nitrogen compounds are based on tricycle fused 1,2,4-triazine and 1,2,4,5-tetrazine heterocycles. In this work, the molecules BTT and DAAT have been studied in order to validate the theoretical approach and to facilitate further progress developments for the molecules of interest. Molecular structural properties are clarified, and IR spectra predictions are provided to help detection of those compounds in the experiment. The energy content of the molecules in the gas phase is evaluated by calculating standard enthalpies of formation, by using a special selection of isodesmic reaction paths. We also include estimates of the condensed-phase heats of formation and heats of sublimation in the framework of the Politzer approach. The obtained properties are consistent with those new high-nitrogen compounds being a promising set of advanced energetic materials.     

Ab initio characterization of the Ne-I2 van der Waals complex: intermolecular potentials and vibrational bound states

A theoretical study of the potential energy surface and bound states is performed for the ground state of the NeI(2) van der Waals (vdW) complex. The three-dimensional interaction energies are obtained from ab initio coupled-cluster, coupled-cluster single double (triple)/complete basis set, calculations using large basis sets, of quadruple- through quintuple-zeta quality, in conjunction with relativistic effective core potentials for the heavy iodine atoms. For the analytical representation of the surface two different schemes, based on fitting and interpolation surface generation techniques, are employed. The surface shows a double-minimum topology for linear and T-shaped configurations. Full variational quantum mechanical calculations are carried out using the model surfaces, and the vibrationally averaged structures and energetics for the NeI(2) isomers are determined. The accuracy of the potential energy surfaces is validated by a comparison between the present results and the corresponding experimental data available. In lieu of more experimental measurements, we also report our results/predictions on higher bound vibrational vdW levels, and the influence of the employed surface on them is discussed.     

Formation of self-assembled chains of tetrathiafulvalene on a Cu(100) surface

Formation of self-assembled chains of tetrathiafulvalene (TTF) on the Cu(100) surface has been investigated by scanning tunneling microscopy and density functional theory calculations that include semiempirical van der Waals (vdW) interaction corrections. The calculations show that the chain structures observed in the experiments can only be explained by including the vdW interactions. The molecules are tilted along the chain in order to achieve maximal intermolecular interaction. The chains are metastable on the surface, which is consistent with the experimental observation that they disappear after annealing. The fact that all TTF chains observed in the experiment are short might be possibly explained by the interplay between the stabilizing vdW molecule-molecule interaction and the destabilizing rearrangement of surface atoms due to the strong molecule-substrate interaction.     

The van der Waals one-fluid model for viscosity in Lennard–Jones fluids: Influence of size and energy parameters

This work aims to estimate the limitations of the van der Waals one-fluid (vdW1) approximation in the prediction of the viscosity of Lennard–Jones (LJ) mixtures. To do so, non-equilibrium molecular dynamics simulations have been performed. Results on mixtures have been compared to those of their equivalent pure fluids (in the sense of the vdW1 model). Several systems (146 configurations) are studied, which are composed of binary and ternary mixtures in various thermodynamic states and for different combining rules. In a first step, deviations induced separately by LJ molecular parameters (size or energy) have been analyzed. It is shown that the vdW1 model is well designed for the energy parameter in every configuration. On the contrary, for the size parameter, deviations induced by this one-fluid approach are shown to be large in dense state and low temperature systems. In a second step, the coupling effects of the LJ size and energy parameters with the mass are studied. It appears that an accurate one-fluid approximation for viscosity should involve a coupling between the mass and the size parameters in its formulation (which is not the case with the vdW1 model) but not between the mass and the energy.     

Monte Carlo Simulations of Adsorbed Solutions in Heterogeneous Porous Materials

We present results of a Monte Carlo simulation study of binary mixtures of ethane and methane in silica gel. The molecular model treats the adsorbent as a matrix of silica microspheres. The adsorption isotherms, adsorption selectivities and isosteric heats of adsorption have been determined for these systems. The results are compared with predictions from the ideal adsorbed solution (IAS) theory and with experiment. The heats of adsorption are accurately described by the IAS theory. The adsorption isotherms are accurately described by the IAS theory at low bulk pressure but the IAS theory overpredicts the density at high bulk pressure. This latter effect is opposite to that observed in bulk mixtures of this type where nonidealities generally lead to a density increase on mixing. The pressure dependence of the selectivity does not exhibit a maximum at low pressure. We discuss this effect in terms of the adsorbent microstructure.     

On the Energetics of nh3 Adsorption and Decomposition on nb(100) Surface : a ubi-qep Study

The energetics of the adsorption and decomposition of ammonia on Nb(100) surface was investigated by the method of unity bond index-quadratic exponential potential (UBI-QEP). The experimental as well as theoretically derived atomic heats of adsorptions were used as the input data and the results were compared with the findings of local density functional theory (LDFT) and non-local density functional theory (NLDFT). The method was capable of correctly predicting the fragmentation of ammonia on 4-fold hollow sites on the surface of Nb(100) and magnitudes of the heats of adsorptions and activation energies were more realistic.     

Vibrational and thermodynamic properties of β-HMX: a first-principles investigation

Thermodynamic properties of β-HMX crystal are investigated using the quasi-harmonic approximation and density functional theory within the local density approximation (LDA), generalized gradient approximation (GGA), and GGA + empirical van der Waals (vdW) correction. It is found that GGA well describes the thermal expansion coefficient and heat capacity but fails to produce correct bulk modulus and equilibrium volume. The vdW correction improves the bulk modulus and volume, but worsens the thermal expansion coefficient and heat capacity. In contrast, LDA describes all thermodynamic properties with reasonable accuracy, and overall is a good exchange-correlation functional for β-HMX molecular crystal. The results also demonstrate significant contributions of phonons to the equation of state. The static calculation of equilibrium volume for β-HMX differs from the room-temperature value incorporating lattice vibrations by over 5%. Therefore, for molecular crystals, it is essential to include phonon contributions when calculated equation of state is compared with experimental data at ambient condition.     

Water monomer interaction with gold nanoclusters from van der Waals density functional theory

We investigate the interaction between water molecules and gold nanoclusters Au(n) through a systematic density functional theory study within both the generalized gradient approximation and the nonlocal van der Waals (vdW) density functional theory. Both planar (n = 6-12) and three-dimensional (3D) clusters (n = 17-20) are studied. We find that applying vdW density functional theory leads to an increase in the Au-Au bond length and a decrease in the cohesive energy for all clusters studied. We classify water adsorption on nanoclusters according to the corner, edge, and surface adsorption geometries. In both corner and edge adsorptions, water molecule approaches the cluster through the O atom. For planar clusters, surface adsorption occurs in a O-up/H-down geometry with water plane oriented nearly perpendicular to the cluster. For 3D clusters, water instead favors a near-flat surface adsorption geometry with the water O atom sitting nearly atop a surface Au atom, in agreement with previous study on bulk surfaces. Including vdW interaction increases the adsorption energy for the weak surface adsorption but reduces the adsorption energy for the strong corner adsorption due to increased water-cluster bond length. By analyzing the adsorption induced charge rearrangement through Bader's charge partitioning and electron density difference and the orbital interaction through the projected density of states, we conclude that the bonding between water and gold nanocluster is determined by an interplay between electrostatic interaction and covalent interaction involving both the water lone-pair and in-plane orbitals and the gold 5d and 6s orbitals. Including vdW interaction does not change qualitatively the physical picture but does change quantitatively the adsorption structure due to the fluxionality of gold nanoclusters.     

Sensitivity of PSA process performance to input variables

Mathematical models for pressure swing adsorption (PSA) processes essentially require the simultaneous solutions of mass, heat and momentum balance equations for each step of the process using appropriate boundary conditions for the steps. The key model input variables needed for estimating the separation performance of the process are the multicomponent adsorption equilibria, kinetics and heats of adsorption for the system of interest. A very detailed model of an adiabatic Skarstrom PSA cycle for production of high purity methane from a ethylene-methane bulk mixture is developed to study the sensitivity of the process performance to the input variables. The adsorption equilibria are described by the heterogeneous Toth model which accounts for variations of isosteric heats of adsorption of the components with adsorbate loading. A linear driving force model is used to describe the kinetics. The study shows that small errors in the heats of adsorption of the components can severely alter the overall performance of the process (methane recovery and productivity). The adsorptive mass transfer coefficients of the components also must be known fairly accurately in order to obtain precise separation performance.     

Thermodynamics of CO2 adsorption on functionalized SBA-15 silica. NLDFT analysis of surface energetic heterogeneity

Siliceous SBA-15 mesoporous molecular sieves were functionalized with different amounts of 3-aminopropyl-trimethoxysilane. To obtain a more detailed insight into the material properties of the prepared samples, their textural parameters were combined with results of thermal analysis. Adsorption isotherms of carbon dioxide on parent and functionalized SBA-15 were measured in the temperature range from 273 to 333 K. From the temperature dependence of CO(2) isotherms the isosteric adsorption heats of CO(2) were determined and discussed. Information about the surface energetic heterogeneity caused by tethered 3-aminopropyl groups were obtained from CO(2) adsorption energy distributions calculated using the theoretical CO(2) adsorption isotherms derived from the non-local density functional theory. The values of isosteric heats and the energy distributions of CO(2) adsorption detect highly energetic sites and enabled quantification of their concentrations.     

Computational investigation of the heat of formation, detonation properties of furazan-based energetic materials

The heats of formation (HOFs) for a series of furazan-based energetic materials were calculated by density functional theory. The isodesmic reaction method was employed to estimate the HOFs. The result shows that the introductions of azo and azoxy groups can increase the HOF, but the introduction of azo group can increase the more HOF, when compared with azoxy group. The detonation velocities and detonation pressures of the furazan-based energetic materials are further evaluated at B3LYP/6-31G* level. Dioxoazotetrafurazan and azoxytetrafurazan may be regarded as the potential candidates of high-energy density materials because of good detonation performance. In addition, there are good linear correlations between OB and detonation velocities, and OB and detonation pressures. The energy gaps between the HOMO and LUMO of the studied compounds are also investigated. These results provide basic information for the molecular design of novel high-energy density materials.