Benchmark ab initio thermochemistry of the isomers of diimide, N2H2, using accurate computed structures and anharmonic force fields

Many-body analysis and spectroscopic characterization of diazene oligomers: A theoretical study

The present study reports the many-body analysis and spectroscopic characterization of linear and cyclic diazene oligomers in gas and water solvent states. The oligomers of diazene from monomer to pentamer have been considered for the study. The spectroscopic studies such as geometrical parameters, infrared spectra, electronic absorption spectra, and natural transition orbitals (NTOs) were reported. Many-body analysis techniques have been implemented to study the interactions among the diazene oligomers. These calculations have been performed using exchange and correlation functional (B3LYP) and 6-311++G (d,p) basis set. The geometrical parameters and infrared modes of monomer diazene in the gas state are well-matched with the available experimental determinations at this level of theory. A significant change in vibrational modes of linear and cyclic diazene oligomers has been observed in the gas phase-to-water solvent state. The time-dependent density functional theory (TD-DFT) has been used to calculate the electronic absorption spectra of diazene oligomers. The Wavelength of electronic transitions, oscillator strength, and HOMO to LUMO gap has been reported. Many-body analysis shows that two-, three-, four-, and five-body energies have a remarkable contribution to the binding energy in addition to relaxation energies. All these calculations have been performed using Gaussian 16 program package.

Theoretical-Computational Modeling of Gas-State Thermodynamics in Flexible Molecular Systems: Ionic Liquids in the Gas Phase as a Case Study

A theoretical-computational procedure based on the quasi-Gaussian entropy (QGE) theory and molecular dynamics (MD) simulations is proposed for the calculation of thermodynamic properties for molecular and supra-molecular species in the gas phase. The peculiarity of the methodology reported in this study is its ability to construct an analytical model of all the most relevant thermodynamic properties, even within a wide temperature range, based on a practically automatic sampling of the entire conformational repertoire of highly flexible systems, thereby bypassing the need for an explicit search for all possible conformers/rotamers deemed relevant. In this respect, the reliability of the presented method mainly depends on the quality of the force field used in the MD simulations and on the ability to discriminate in a physically coherent way between semi-classical and quantum degrees of freedom. The method was tested on six model systems (n-butane, n-butane, n-octanol, octadecane, 1-butyl-3-methylimidazolium hexafluorophosphate and 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic pairs), which, being experimentally characterized and already addressed by other theoretical-computational methods, were considered as particularly suitable to allow us to evaluate the method's accuracy and efficiency, bringing out advantages and possible drawbacks. The results demonstrate that such a physically coherent yet relatively simple method can represent a further valid computational tool that is alternative and complementary to other extremely efficient computational methods, as it is particularly suited for addressing the thermodynamics of gaseous systems with a high conformational complexity over a large range of temperature.

Theoretical investigation of the isomerization pathways of diazenes: torsion vs. inversion

Diazenes are an important family of organic compounds used widely in synthetic and materials chemistry. These molecules have a planar geometry and exhibit cis-trans isomerization. The simplest of all these molecules - diazene (N2H2) - has been subjected to several experimental and theoretical studies. Two mechanisms have been proposed for the cis-trans isomerization of diazene, which are an in-plane inversion and an out-of-plane torsion. The activation energies for these pathways are similar and the competition between these two mechanisms has been discussed in the literature based on electronic structure theory calculations. Three decades ago, a classical dynamics investigation of diazene isomerization was carried out using a model Hamiltonian and it was indicated that the in-plane inversion is forbidden classically because of a centrifugal barrier and the out-of-plane torsion is the only isomerization pathway. In the present work, we investigated the cis-trans isomerization dynamics of diazene using ab initio classical trajectory simulations at the CASSCF(2,2)/aug-cc-pVDZ level of electronic structure theory. The simulation results confirmed the presence of the aforementioned centrifugal barrier for the inversion and torsion was the only observed pathway. The calculations were repeated for a similar system (difluorodiazene, N2F2) and again the centrifugal barrier prevented the inversion pathway.

Combined Electron-Diffraction, Spectroscopic, and Theoretical Determination of the Structure of N-Deuterio- trans-Methyldiazene, CH3N═ND. Conformational Effects of the N═N Double Bond

Gas phase electron-diffraction (GED) data obtained at a nozzle-tip temperature of 273 K have been combined with spectroscopic vibrational-rotational constants to determine the structure of trans-methyldiazene, an important prototype for the N═N double bond. The N-deuterio form CH₃N═ND was used in the study since it is appreciably more stable than CH₃N═NH. Both the theoretical and experimental results are consistent with a planar C _s trans-CNND framework. The experimental results ( r_α⁰/ r_g²⁷³) are 1.465(2)/1.467(2) Å for the CN bond, 1.248(1)/1.251(1) Å for the N═N double bond, and 1.037(17)/1.048(17) Å for the ND bond. The NND angle is 105.9(20)/105.6(20)° and the CNN angle is 112.4(5)/112.2(5)°, where the uncertainties in parentheses are twice the standard deviation from a simultaneous least-squares fit of the GED and microwave data. For the methyl group, both theory and experiment indicate that two CH bonds are symmetrically arranged out of the molecular plane while the third CH' lies in the plane in an eclipsed (not staggered) cis-H'CNN arrangement. Theoretical calculations (B3LYP/cc-PVnZ and CCSD(T)/cc-PVnZ) suggest a slight distortion of the methyl group, with a tilt of the methyl top axis about 5° away from the N═N bond. The experimental data are consistent with this picture but are equally consistent with an undistorted methyl group. Inclusion of distortions predicted by theory in a complete basis set limit (CBS) lead to a preferred analysis with average values of 1.086(5)/1.106(5) Å for the CH bond length and an average HCH angle of 108.3(8)/107.8(8)°. Features of the structure of methyldiazene and related compounds are discussed. It is found that the short N═N bond length in the diazenes produces much greater steric repulsion than in analogous ethylene compounds and this effect leads to some interesting conformational and distortion differences for attached CH₃ groups.

Bias-Free Chemically Diverse Test Sets from Machine Learning

Current benchmarking methods in quantum chemistry rely on databases that are built using a chemist's intuition. It is not fully understood how diverse or representative these databases truly are. Multivariate statistical techniques like archetypal analysis and K-means clustering have previously been used to summarize large sets of nanoparticles however molecules are more diverse and not as easily characterized by descriptors. In this work, we compare three sets of descriptors based on the one-, two-, and three-dimensional structure of a molecule. Using data from the NIST Computational Chemistry Comparison and Benchmark Database and machine learning techniques, we demonstrate the functional relationship between these structural descriptors and the electronic energy of molecules. Archetypes and prototypes found with topological or Coulomb matrix descriptors can be used to identify smaller, statistically significant test sets that better capture the diversity of chemical space. We apply this same method to find a diverse subset of organic molecules to demonstrate how the methods can easily be reapplied to individual research projects. Finally, we use our bias-free test sets to assess the performance of density functional theory and quantum Monte Carlo methods.

A Database of Formation Enthalpies of Nitrogen Species by Compound Methods (CBS-QB3, CBS-APNO, G3, G4)

Accurate thermochemical data for compounds containing C/H/N/O are required to underpin kinetics simulation and modeling of the reactions of these species in different environments. There is a dearth of experimental data so computational quantum chemistry has stepped in to fill this breach and to verify whether particular experiments are in need of revision. A number of composite model chemistries (CBS-QB3, CBS-APNO, G3, and G4) are used to compute theoretical atomization energies and hence enthalpies of formation at 0 and 298.15 K, and these are benchmarked against the best available compendium of values, the Active Thermochemical Tables or ATcT. In general the agreement is very good for some 28 species with the only discrepancy being for hydrazine. It is shown that, although individually the methods do not perform that well, collectively the mean unsigned error is <1.7 kJ mol(-1); hence, this approach provides a useful tool to screen published values and validate new experimental results. Using multiple model chemistries does have some drawbacks but can produce good results even for challenging molecules like HOON and CN2O2. The results for these smaller validated molecules are then used as anchors for determining the formation enthalpies of larger species such as methylated hydrazines and diazenes, five- and six-membered heterocyclics via carefully chosen isodesmic working reactions with the aim of resolving some discrepancies in the literature and establishing a properly validated database. This expanded database could be useful in testing the performance of computationally less-demanding density function methods with newer functionals that have the capacity to treat much larger systems than those tested here.

The heats of formation of diazene, hydrazine, N2H3+, N2H5+, N2H, and N2H3 and the Methyl Derivatives CH3NNH, CH3NNCH3, and CH3HNNHCH3

The heats of formation of N(2)H, diazene (cis- and trans-N(2)H(2)), N(2)H(3), and hydrazine (N(2)H(4)), as well as their protonated species (diazenium, N(2)H(3)(+), and hydrazinium, N(2)H(5)(+)), have been calculated by using high level electronic structure theory. Energies were calculated by using coupled cluster theory with a perturbative treatment of the triple excitations (CCSD(T)) and employing augmented correlation consistent basis sets (aug-cc-pVnZ) up to quintuple-zeta, to perform a complete basis set extrapolation for the energy. Geometries were optimized at the CCSD(T) level with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. Core-valence and scalar relativistic corrections were included, as well as scaled zero point energies. We find the following heats of formation (kcal/mol) at 0 (298) K: DeltaH(f)(N(2)H) = 60.8 (60.1); DeltaH(f)(cis-N(2)H(2)) = 54.9 (53.2); DeltaH(f)(trans-N(2)H(2)) = 49.9 (48.1) versus >/=48.8 +/- 0.5 (exptl, 0 K); DeltaH(f)(N(2)H(4)) = 26.6 (23.1) versus 22.8 +/- 0.2 (exptl, 298 K); DeltaH(f)(N(2)H(3)) = 56.2 (53.6); DeltaH(f)(N(2)H(3)(+)) = 231.6 (228.9); and DeltaH(f)(N(2)H(5)(+)) = 187.1 (182.7). In addition, we calculated the heats of formation of CH(3)NH(2), CH(3)NNH, and CH(3)HNNHCH(3) by using isodesmic reactions and at the G3(MP2) level. The calculated results for the hydrogenation reaction RNNR + H(2) --> RHNNHR show that substitution of an organic substituent for H improved the energetics, suggesting that these types of compounds may be possible to use in a chemical hydrogen storage system.

Predicting Catalysis: Understanding Ammonia Synthesis from First-Principles Calculations

Here, we give a full account of a large collaborative effort toward an atomic-scale understanding of modern industrial ammonia production over ruthenium catalysts. We show that overall rates of ammonia production can be determined by applying various levels of theory (including transition state theory with or without tunneling corrections, and quantum dynamics) to a range of relevant elementary reaction steps, such as N(2) dissociation, H(2) dissociation, and hydrogenation of the intermediate reactants. A complete kinetic model based on the most relevant elementary steps can be established for any given point along an industrial reactor, and the kinetic results can be integrated over the catalyst bed to determine the industrial reactor yield. We find that, given the present uncertainties, the rate of ammonia production is well-determined directly from our atomic-scale calculations. Furthermore, our studies provide new insight into several related fields, for instance, gas-phase and electrochemical ammonia synthesis. The success of predicting the outcome of a catalytic reaction from first-principles calculations supports our point of view that, in the future, theory will be a fully integrated tool in the search for the next generation of catalysts.

Development of density functionals for thermochemical kinetics

A density functional theory exchange-correlation functional for the exploration of reaction mechanisms is proposed. This functional, denoted BMK (Boese-Martin for Kinetics), has an accuracy in the 2 kcal/mol range for transition state barriers but, unlike previous attempts at such a functional, this improved accuracy does not come at the expense of equilibrium properties. This makes it a general-purpose functional whose domain of applicability has been extended to transition states, rather than a specialized functional for kinetics. The improvement in BMK rests on the inclusion of the kinetic energy density together with a large value of the exact exchange mixing coefficient. For this functional, the kinetic energy density appears to correct "back" the excess exact exchange mixing for ground-state properties, possibly simulating variable exchange.