Vacuum‐ultraviolet photolysis of N2H2: Generation of NH fragments

Extrapolation in quantum chemistry: Insights on energetics and reaction dynamics

Since there is no exact solution for problems in physics and chemistry, extrapolation methods may assume a key role in quantitative quantum chemistry. Two topics where it bears considerable impact are addressed, both at the heart of computational quantum chemistry: electronic structure and reaction dynamics. In the first, the problem of extrapolating the energy obtained by solving the electronic Schrödinger equation to the limit of the complete one-electron basis set is addressed. With the uniform-singlet-and-triplet-extrapolation (USTE) scheme at the focal point, the emphasis is on recent updates covering from the energy itself to other molecular properties. The second topic refers to extrapolation of quantum mechanical reactive scattering probabilities from zero total angular momentum to any of the values that it may assume when running quasiclassical trajectories, QCT/QM-[Formula: see text]J. With the extrapolation guided in both cases by physically motivated asymptotic theories, realism is seeked by avoiding unsecure jumps into the unknown. Although, mostly review oriented, a few issues are addressed for the first time here and there. Prospects for future work conclude the overview.

Optimal basis sets for CBS extrapolation of the correlation energy: oVxZ and oV(x+d)Z

We seek correlation consistent double- and triple-zeta basis sets that perform optimally for extrapolating the correlation energy to the one-electron complete basis set limit. Since the methods used are approximate, the novel basis sets become method specific in the sense of performing best for the chosen level of theory. Such basis sets are also shown to perform accurately for tensorial properties and do not significantly alter the Hartree-Fock energy. Quantitatively, the extrapolated correlation energies from (oVdZ, oVtZ) outperform typically by three- to fivefold those obtained from traditional ansatzes with similar flexibility, thus being (VtZ, VqZ) type or even better. They may even outperform explicitly correlated ones. Not surprisingly, the outperformance in relative energies (e.g., atomization and dissociation energies, and ionization potential) is somewhat downscaled, albeit consistently better than with traditional basis sets. As a case study, we also consider the polarizability of p-nitroaniline, a sizeable system for which complete basis set (CBS)(oVdZ, oVtZ) calculations are shown to outperform equally expensive CBS(VdZ, VtZ) results.

Accurate ab initio based DMBE potential energy surface for the ground electronic state of N2H2

A global single-sheeted double many-body expansion potential energy surface is reported for the ground electronic state of N(2)H(2). Starting from an approximate cluster expansion of the molecular potential that utilizes previously reported functions of the same family for the triatomic fragments, four-body energy terms have been calibrated from extensive accurate ab initio data so as to reproduce the main features of the title system. The switching function formalism previously suggested for three-body systems [A. J. C. Varandas and L. Poveda, Theor. Chem. Acc. 116, 404 (2006)] has been generalized to approximate the true multisheeted nature of N(2)H(2) potential energy surface, thus allowing the correct behavior at the N((2)D) + NH(2)((2)A(")) and N((4)S) + NH(2)((4)A(")) dissociation limits. The resulting fully six-dimensional potential energy function reproduces the correct symmetry under permutation of identical atoms and predicts the main stationary points of the molecule in the valence and long-range regions in good agreement with available experimental and theoretical data on the diazene molecule.

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.

Density Functional Theory Study of the Properties of N−H···N, Noncooperativities, and Intermolecular Interactions in Linear trans-Diazene Clusters up to Ten Molecules

We investigate aspects of N-H...N hydrogen bonding in the linear trans-diazene clusters (n=2-10) such as the N...H and N-H lengths, n(N) --> sigma(N-H) interactions, N...H strengths, and frequencies of the N-H stretching vibrations utilizing the DFT/B3LYP theory, the natural bond orbital (NBO) method, and the theory of atoms in molecules (AIM). Our calculations indicate that the structure and energetics are qualitatively different from the conventional H-bonded systems, which usually exhibit distinct cooperative effects, as cluster size increases. First, a shortening rather than lengthening of the N-H bond is found and thus a blue rather than red shift is predicted. Second, for the title clusters, any sizable cooperative changes in the N-H and N...H lengths, n(N) --> sigma(N-H) charge transfers, N...H strengths, and frequencies of the N-H stretching vibrations for the linear H-bonded trans-diazene clusters do not exist. Because the n(N) --> sigma(N-H) interaction hardly exhibits cooperative effects, the capability of the linear trans-diazene cluster to localize electrons at the N...H bond critical point is almost independent of cluster size and thereby leads to the noncooperative changes in the N...H lengths and strengths and the N-H stretching frequencies. Third, the dispersion energy is sizable and important; more than 30% of short-range dispersion energy not being reproduced by the DFT leads to the underestimation of the interaction energies by DFT/B3LYP. The calculated nonadditive interaction energies show that, unlike the conventional H-boned systems, the trans-diazene clusters indeed exhibit very weak nonadditive interactions.

Besides N(2), What Is the Most Stable Molecule Composed Only of Nitrogen Atoms?

Polynitrogen molecules have been studied systematically at high levels of ab initio and density functional theory (DFT). Besides N(2), the thermodynamically most stable N(n)() molecules, located with the help of a newly developed energy increment system, are all based on pentazole units. The geometric, energetic, and magnetic criteria establish pentazole (2) and its anion (3) to be as aromatic as their isoelectronic analogues, e.g., furan, pyrrole, and the cyclopentadienyl anion. The bond lengths in 2 and 3 are equalized; both have large aromatic stabilization energies (ASE) and also substantial magnetic susceptibility exaltations (Lambda). The C(s)() symmetric azidopentazole (14), a candidate for experimental investigation, is the lowest energy N(8) isomer but is still 196.7 kcal/mol higher in energy than four N(2) molecules. Octaazapentalene (12) with 10 pi electrons also is aromatic. The D(2)(d)() symmetric bispentazole (21) is the lowest energy N(10) minimum but is 260 kcal/mol higher in energy than five N(2) molecules. For strain-free molecules, the average deviation is +/-2.6 kcal/mol between the DFT energies and those based on the increment scheme. The increment scheme also provides estimates of the strain energies of polynitrogen compounds, e.g., tetraazatetrahedrane (8, 48.2 kcal/mol), octaazacubane (11, 192.6 kcal/mol), and N(20) (27, 294.6 kcal/mol), and is useful in searching for new high-energy-high-density materials.