Active Thermochemical Tables: The Partition Function of Hydroxymethyl (CH2OH) Revisited.

The best currently available set of temperature-dependent nonrigid rotor anharmonic oscillator (NRRAO) thermochemical and thermophysical properties of hydroxymethyl radical is presented. The underlying partition function relies on a critically evaluated complement of accurate experimental and theoretical data and is constructed using a two-pronged strategy that combines contributions from large amplitude motions obtained from direct counts, with contributions from the other internal modes of motion obtained from analytic NRRAO expressions. The contributions from the two strongly coupled large-amplitude motions of CH2OH, OH torsion and CH2 wag, are based on energy levels obtained by solving the appropriate two-dimensional projection of a fully dimensional potential energy surface that was recently obtained at the CCSD(T)/cc-pVTZ level of theory. The contributions of the remaining seven, more rigid, vibrational modes and of the external rotations are captured by NRRAO corrections to the standard rigid rotor harmonic oscillator (RRHO) treatment, which include corrections for vibrational anharmonicities, rotation-vibration interaction, Coriolis effects, and low temperature. The basic spectroscopic constants needed for the construction of the initial RRHO partition function rely on experimental ground-state rotational constants and the best available experimental fundamentals, additionally complemented by fundamentals obtained from the variational solution of the full-dimensional potential energy surface using a recently developed two-component multilayer Lanczos algorithm. The higher-order spectroscopic constants necessary for the NRRAO corrections are extracted from a second-order variational perturbation treatment (VPT2) of the same potential energy surface. The Lanczos solutions of the fully dimensional surface are validated against available experimental data, and the VPT2 results and the solutions of the reduced dimensionality surface are validated both against the Lanczos solutions and available experiments. The NRRAO thermophysical and thermochemical properties, given both in tabular form and as seven- and nine-coefficient NASA polynomials, are compared to previous results. In addition, the latest ATcT values for the enthalpy of formation of CH2OH at 298.15 K (0 K), -16.75 ± 0.27 kJ/mol (-10.45 ± 0.27 kJ/mol), and of other related CH nO m species ( n = 0-4, m = 0,1) are reported, together with a plethora of related bond dissociation enthalpies (BDEs), such as the C-H, O-H, and C-O bond dissociation enthalpies of methanol, 402.16 ± 0.26 kJ/mol (395.61 ± 0.26 kJ/mol), 440.34 ± 0.26 kJ/mol (434.86 ± 0.26 kJ/mol), and 384.85 ± 0.15 kJ/mol (377.14 ± 0.15 kJ/mol), respectively, and analogous BDEs for hydroxymethyl, 343.67 ± 0.37 kJ/mol (339.16 ± 0.37 kJ/mol), 125.54 ± 0.28 kJ/mol (121.11 ± 0.28 kJ/mol), and 445.86 ± 0.29 kJ/mol (438.76 ± 0.29 kJ/mol), respectively. The reasons governing the alternation between strong and weak sequential H atom BDEs of methanol are also discussed.

A rigorous full-dimensional quantum dynamics study of tunneling splitting of rovibrational states of vinyl radical C2H3.

We report a rigorous quantum mechanical study of the rovibrational energy levels of vinyl radical C2H3. The calculations are carried out using a real two-component multi-layer Lanczos algorithm in a set of orthogonal polyspherical coordinates based on a recently developed accurate ab initio potential energy surface of C2H3. All well converged 158 vibrational bands up to 3200 cm-1 are determined, together with a comparison to previous calculations and experimental results. Results show a remarkable multi-dimensional tunneling effect on the vibrational spectra of the radical. The vibrational tunneling splitting is substantially different from that of previous reduced dimensional calculations. The rotational constants of the fundamental vibrational bands of C2H3 are also given. It was found that the rovibrational states are strongly coupled, especially among those bending vibrational modes. In addition, the perturbative iteration approach of Gruebele has been extended to assign the rovibrational energy levels of C2H3 without the requirement of explicit wavefunctions.

A coherent discrete variable representation method on a sphere.

The coherent discrete variable representation (ZDVR) has been extended for constructing a multidimensional potential-optimized DVR basis on a sphere. In order to deal with the non-constant Jacobian in spherical angles, two direct product primitive basis methods are proposed so that the original ZDVR technique can be properly implemented. The method has been demonstrated by computing the lowest states of a two-dimensional (2D) vibrational model. Results show that the extended ZDVR method gives accurate eigenvalues and exponential convergence with increasing ZDVR basis size.

An exact variational method to calculate rovibrational spectra of polyatomic molecules with large amplitude motion.

We report a new full-dimensional variational algorithm to calculate rovibrational spectra of polyatomic molecules using an exact quantum mechanical Hamiltonian. The rovibrational Hamiltonian of system is derived in a set of orthogonal polyspherical coordinates in the body-fixed frame. It is expressed in an explicitly Hermitian form. The Hamiltonian has a universal formulation regardless of the choice of orthogonal polyspherical coordinates and the number of atoms in molecule, which is suitable for developing a general program to study the spectra of many polyatomic systems. An efficient coupled-state approach is also proposed to solve the eigenvalue problem of the Hamiltonian using a multi-layer Lanczos iterative diagonalization approach via a set of direct product basis set in three coordinate groups: radial coordinates, angular variables, and overall rotational angles. A simple set of symmetric top rotational functions is used for the overall rotation whereas a potential-optimized discrete variable representation method is employed in radial coordinates. A set of contracted vibrationally diabatic basis functions is adopted in internal angular variables. Those diabatic functions are first computed using a neural network iterative diagonalization method based on a reduced-dimension Hamiltonian but only once. The final rovibrational energies are computed using a modified Lanczos method for a given total angular momentum J, which is usually fast. Two numerical applications to CH4 and H2CO are given, together with a comparison with previous results.

Full-Dimensional Quantum Calculations of Vibrational Levels of NH4(+) and Isotopomers on An Accurate Ab Initio Potential Energy Surface.

Vibrational energy levels of the ammonium cation (NH4(+)) and its deuterated isotopomers are calculated using a numerically exact kinetic energy operator on a recently developed nine-dimensional permutation invariant semiglobal potential energy surface fitted to a large number of high-level ab initio points. Like CH4, the vibrational levels of NH4(+) and ND4(+) exhibit a polyad structure, characterized by a collective quantum number P = 2(v1 + v3) + v2 + v4. The low-lying vibrational levels of all isotopomers are assigned and the agreement with available experimental data is better than 1 cm(-1).

Rotational and angular distributions of NO products from NO-Rg (Rg = He, Ne, Ar) complex photodissociation.

We present the results of an investigation into the rotational and angular distributions of the NO à state fragment following photodissociation of the NO-He, NO-Ne, and NO-Ar van der Waals complexes excited via the à ← X̃ transition. For each complex, the dissociation is probed for several values of Ea, the available energy above the dissociation threshold. For NO-He, the Ea values probed were 59, 172, and 273 cm(-1); for NO-Ne they were 50 and 166 cm(-1); and for NO-Ar they were 44, 94, 194, and 423 cm(-1). The NO à state rotational distributions arising from NO-He are cold, with most products in low angular momentum states. NO-Ne leads to broader NO rotational distributions but they do not extend to the maximum possible given the energy available. In the case of NO-Ar, the distributions extend to the maximum allowed at that energy and show the unusual shapes associated with the rotational rainbow effect reported in previous studies. This is the only complex for which a rotational rainbow effect is observed at the chosen Ea values. Product angular distributions have also been measured for the NO à photodissociation product for the three complexes. NO-He produces nearly isotropic fragments, but the anisotropy parameter, ß, for NO-Ne and NO-Ar photofragments shows a surprising change in sign from negative to positive as Ea increases within the unstructured excitation profile. Franck-Condon selection of a broader distribution of geometries including more linear geometries at lower excitation energies and more T-shaped geometries at higher energies can account for the changing recoil anisotropy. Two-dimensional wavepacket calculations are reported to model the rotational state distributions and the bound-continuum absorption spectra.

Vibrational energy levels of the simplest Criegee intermediate (CH2OO) from full-dimensional Lanczos, MCTDH, and MULTIMODE calculations.

Accurate vibrational energy levels of the simplest Criegee intermediate (CH2OO) were determined on a recently developed ab initio based nine-dimensional potential energy surface using three quantum mechanical methods. The first is the iterative Lanczos method using a conventional basis expansion with an exact Hamiltonian. The second and more efficient method is the multi-configurational time-dependent Hartree (MCTDH) method in which the potential energy surface is refit to conform to the sums-of-products requirement of MCTDH. Finally, the energy levels were computed with a vibrational self-consistent field/virtual configuration interaction method in MULTIMODE. The low-lying levels obtained from the three methods are found to be within a few wave numbers of each other, although some larger discrepancies exist at higher levels. The calculated vibrational levels are very well represented by an anharmonic effective Hamiltonian.

Accurate quantum dynamics calculations of vibrational spectrum of dideuteromethane CH2D2.

We report a rigorous variational study of the infrared (IR) vibrational spectra of both CH2D2 and (13)CH2D2 isotopomers using an exact molecular Hamiltonian. Calculations are carried out using a recently developed multi-layer Lanczos algorithm based on the accurate refined Wang and Carrington potential energy surface of methane and the low-order truncated ab initio dipole moment surface of Yurchenko et al. [J. Mol. Spectrosc. 291, 69 (2013)]. All well converged 357 vibrational energy levels up to 6100 cm(-1) of CH2D2 are obtained, together with a comparison to previous calculations and 91 experimental bands available. The calculated frequencies are in excellent agreement with the experimental results and give a root-mean-square error of 0.67 cm(-1). In particular, we also compute the transition intensities from the vibrational ground state for both isotopomers. Based on the theoretical results, 20 experimental bands are suggested to be re-assigned. Surprisingly, an anomalous C isotopic effect is discovered in the nν5 modes of CH2D2. The predicted IR spectra provide useful information for understanding those unknown bands.

Neural network iterative diagonalization method to solve eigenvalue problems in quantum mechanics.

We propose a multi-layer feed-forward neural network iterative diagonalization method (NNiDM) to compute some eigenvalues and eigenvectors of large sparse complex symmetric or Hermitian matrices. The NNiDM algorithm is developed by using the complex (or real) guided spectral transform Lanczos (cGSTL) method, thick restart technique, and multi-layered basis contraction scheme. Artificial neurons (or nodes) are defined by a set of formally orthogonal Lanczos polynomials, where the biases and weights are dynamically determined through a series of cGSTL iterations and small matrix diagonalizations. The algorithm starts with one random vector. The last output layer produces wanted eigenvalues and eigenvectors near a given reference value via a linear transform diagonalization approach. Since the algorithm uses the spectral transform technique, it is capable of computing interior eigenstates in dense spectrum regions. The general NNiDM algorithm is applied for calculating energies, widths, and wavefunctions of two typical molecules HO2 and CH4 as examples.

Doppler-Resolved Kinetics of Saturation Recovery.

Frequency-modulated laser transient absorption has been used to monitor the ground-state rotational energy-transfer rates of CN radicals in a double-resonance, depletion recovery experiment. When a pulsed laser is used to burn a hole in the equilibrium ground-state population of one rotational state without velocity selection, the population recovery rate is found to depend strongly on the Doppler detuning of a narrow-band probe laser. Similar effects should be apparent for any relaxation rate process that competes effectively with velocity randomization. Alternative methods of extracting thermal rate constants in the presence of these non-thermal conditions are evaluated. Total recovery rate constants, analogous to total removal rate constants in an experiment preparing a single initial rotational level, are in good agreement with quantum scattering calculations, but are slower than previously reported experiments and show qualitatively different rotational state dependence between Ar and He collision partners. Quasi-classical trajectory studies confirm that the differing rotational state dependence is primarily a kinematic effect.

Multi-layer Lanczos iteration approach to calculations of vibrational energies and dipole transition intensities for polyatomic molecules.

We report a rigorous full dimensional quantum dynamics algorithm, the multi-layer Lanczos method, for computing vibrational energies and dipole transition intensities of polyatomic molecules without any dynamics approximation. The multi-layer Lanczos method is developed by using a few advanced techniques including the guided spectral transform Lanczos method, multi-layer Lanczos iteration approach, recursive residue generation method, and dipole-wavefunction contraction. The quantum molecular Hamiltonian at the total angular momentum J = 0 is represented in a set of orthogonal polyspherical coordinates so that the large amplitude motions of vibrations are naturally described. In particular, the algorithm is general and problem-independent. An application is illustrated by calculating the infrared vibrational dipole transition spectrum of CH4 based on the ab initio T8 potential energy surface of Schwenke and Partridge [Spectrochimica Acta, Part A 57, 887 (2001)] and the low-order truncated ab initio dipole moment surfaces of Yurchenko et al. [J. Mol. Spectrosc. 291, 69 (2013)]. A comparison with experiments is made. The algorithm is also applicable for Raman polarizability active spectra.

Origin of anomalous electronic circular dichroism spectrum of RuPt2(tppz)2Cl2(PF6)4 in acetonitrile.

We report a theoretical study of the structures, energetics, and electronic spectra of the Pt(II)/Ru(II) mixed-metal complex RuPt2(tppz)2Cl2(PF6)4 (tppz = 2,3,5,6-tetra(2-pyridyl)pyrazine) in acetonitrile. The hybrid B3LYP density functional theory and its TDDFT methods were used with a complete basis set (CBS) extrapolation scheme and a conductor polarizable continuum model (C-PCM) for solvation effects. Results showed that the trinuclear complex has four types of stable conformers and/or enantiomers. They are separated by high barriers owing to the repulsive H/H geometrical constraints in tppz. A strong entropy effect was found for the dissociation of RuPt2(tppz)2Cl2(PF6)n in acetonitrile. The UV-visible and emission spectra of the complex were also simulated. They are in good agreement with experiments. In this work we have largely focused on exploring the origin of anomalous electronic circular dichroism (ECD) spectra of the RuPt2(tppz)2Cl2(PF6)4 complex in acetonitrile. As a result, a new mechanism has been proposed together with a clear illustration by using a physical model.

A complex guided spectral transform Lanczos method for studying quantum resonance states.

A complex guided spectral transform Lanczos (cGSTL) algorithm is proposed to compute both bound and resonance states including energies, widths, and wavefunctions. The algorithm comprises of two layers of complex-symmetric Lanczos iterations. A short inner layer iteration produces a set of complex formally orthogonal Lanczos polynomials. They are used to span the guided spectral transform function determined by a retarded Green operator. An outer layer iteration is then carried out with the transform function to compute the eigen-pairs of the system. The guided spectral transform function is designed to have the same wavefunctions as the eigenstates of the original Hamiltonian in the spectral range of interest. Therefore, the energies and/or widths of bound or resonance states can be easily computed with their wavefunctions or by using a root-searching method from the guided spectral transform surface. The new cGSTL algorithm is applied to bound and resonance states of HO2, and compared to previous calculations.

Electronic structure and quantum dynamics of photoinitiated dissociation of O2 on rutile TiO2 nanocluster.

The adsorption and photoinitiated dissociation of molecular oxygen on reduced rutile TiO2 nanocluster have been studied using a hybrid density functional theory (DFT)/time-dependent DFT approach and a time-dependent wavepacket dynamics method. Results show that the most favorable state for O2 at the bridging row O-vacancy site of TiO2 is O2(2-) with an orientation parallel to the surface. We find that its dissociation in the electronic ground state involves a spin forbidden intersystem crossing, and therefore has a large barrier along the reaction pathway. However, time-dependent wavepacket calculations reveal that the photoinitiated O2 dissociation on TiO2 is very fast via a direct mechanism on the excited states. The lifetime of excited O2 molecules is predicted to be about 266 fs. Non-adiabatic effects among the singlet electronic states are found to play an important role in the O2 dissociation whereas the spin-orbit effect is negligible. In addition, adsorption of two O2 molecules at an O-vacancy site shows that the second O2 molecule can stabilize the system by about 0.22 eV.

Exploring the ring-opening pathways in the reaction of morpholinyl radicals with oxygen molecule.

Quantum chemistry calculations using hybrid density functional theory and the coupled-cluster method have been performed to investigate the ring-opening pathways in the oxidation of morpholine (1-oxa-4-aza-cyclohexane). Hydrogen abstraction can form two different carbon-centered radicals, morpholin-2-yl or morpholin-3-yl, or the nitrogen-centered radical, morpholin-4-yl, none of which are found to have low-energy pathways to ring-opening. Extensive exploration of multiple reaction pathways following molecular oxygen addition to these three radicals revealed two competitive low energy pathways to ring-opening. Addition of O(2) to either carbon-centered radical, followed by a 1,4-H shifting mechanism can yield a long-lived cyclic epoxy intermediate, susceptible to ring-opening, following further radical attack. In particular, the second pathway begins with O(2) attack on morpholin-2-yl, followed by a 1,5-H shift and a unimolecular ring-opening without having to overcome a high barrier, releasing a significant amount of heat in the overall ring-opening reaction. The calculations provide valuable context for the development of mechanisms for the low temperature combustion chemistry of nitrogen and oxygen-containing fuels.

The kinetics study of the S + S2 → S3 reaction by the chaperone mechanism.

The recombination of S atoms has been found to be stepwise from the smallest unit, the elemental S atom, to the most abundant molecule S(8). The reaction between S + S(2) → S(3) has not been reported either experimentally or by theory, but may be a key intermediate step in the formation of sulfur aerosols in low-O(2) atmospheres. In this work, the kinetics of this reaction is reported with Ar gas used as the chaperone molecule in the production of S(3) via two complex intermediates: SAr + S(2) and S(2)Ar + S. Quasi-classical and classical trajectory methods are used. The rate constant of the S + S(2) + Ar → S(3) + Ar reaction is determined to be 2.66 × 10(-33) cm(6) mol(-1) s(-1) at 298.15 K. The temperature dependence of the reaction is found to be 2.67 × 10(-33) exp[143.56(1∕T-1∕298.15)]. The second-order rate constant of S + S(2) → S(3) is 6.47 × 10(-14) cm(3) molecule(-1) s(-1) at 298.15 K and the Arrhenius-type rate constant is calculated to be 6.25 × 10(-14) exp[450.15(1∕T-1∕298.15)] cm(3) molecule(-1) s(-1). This work provides a rate coefficient for a key intermediate species in studies of sulfur formation in the modern Venus atmosphere and the primitive Earth atmosphere, for which assumed model rate coefficients have spanned nearly 4 orders of magnitude. Although a symmetry-induced mass-independent isotope effect is not expected for a chaperone mechanism, the present work is an important step toward evaluating whether mass-independence is expected for thiozone formation as is observed for ozone formation.

HOCO radical chemistry.

Free radicals are important species in atmospheric chemistry, combustion, plasma environments, interstellar clouds, and biochemistry. Therefore, researchers would like to understand the formation mechanism, structure, stability, reactivity, spectroscopy, and dynamics of these chemical species. However, due to the presence of one or more unpaired electrons, radicals are often very reactive and have short lifetimes, which makes it difficult to conduct experiments. The HOCO radical appears in the atmosphere as well as in combustion environments and plays an important role in the conversion of CO to CO(2). Through the interplay between theoretical and experimental investigations, researchers have only recently understood the chemical role of the HOCO radical. In this Account, we systematically describe the current state of knowledge of the HOCO radical based on recent theoretical and experimental studies. This radical's two stable conformers, trans- and cis-HOCO, have been identified by high-level ab initio calculations and experimental spectroscopy. trans-HOCO is more stable by approximately 1.8 kcal/mol. The heat of formation of HOCO (298 K) was determined to be -43.0 ± 0.5 kcal/mol, giving a potential well depth of 30.1 ± 0.5 kcal/mol relative to the asymptote of the reactants OH + CO. The HOCO radical is very reactive. In most reactions between the HOCO radical and atoms, the HOCO radical acts as a hydrogen donor to reaction partners. Generally, the hydrogen is transferred through the formation of an association intermediate, which then proceeds through a molecular elimination step to produce the reaction products. The reaction rates of HOCO with some small radicals fall in the range of 10(-11)-10(-10) cm(3) molecule(-1) s(-1). These results clearly illustrate important features in the reactivity of the HOCO radical with other molecules.

Ab initio and RRKM study of the reaction of ClO with HOCO radicals.

The reaction pathways for the ClO + HOCO reaction have been explored using the coupled-cluster method to locate and optimize the critical points on the ground-state potential-energy surface. Results show that the ClO + HOCO reaction can produce Cl + HOC(O)O, HOCl + CO(2), HCl + CO(3), and HClO + CO(2) via an addition or a direct hydrogen abstraction reaction mechanism. The reaction kinetics has been studied using the variational RRKM theory. It is found that the ClO + HOCO reaction is fast and has a negative temperature dependence at low temperatures. At room temperature, the thermal rate coefficient is obtained as 4.26 x 10(-12) cm(3) molecules(-1) s(-1) with product branching fractions of Cl (0.518), HOCl (0.469), HCl (0.01), and HClO (0.003) at zero pressure. The Cl + HOC(O)O products are major, compared to the HOCl + CO(2) products, because of the loose transition state along the dissociation pathway to eliminate Cl. In addition, the RRKM/master equation simulations indicate that the stabilization of the HOC(O)OCl intermediates is noticeable at moderate pressures as its thermal rate constants reach about 6.0 x 10(-13) cm(3) molecules(-1) s(-1). In contrast, the other product branching ratios for the ClO + HOCO reaction are weakly dependent on pressure.

Spherical electron cloud hopping molecular dynamics simulation on dissociative recombination of protonated water.

Dissociative recombination (DR) of H(3)O(+) with electrons at zero collision energy has been studied by a direct ab initio molecular dynamics method on four low-lying electronic states of the system. Initial conditions for trajectories are determined by a spherical electron cloud hopping (SECH) model, while nonadiabatic effects are considered through a surface hopping scheme. The energies, forces, and nonadiabatic coupling strengths (NACS) used in trajectory propagations are calculated on-the-fly via state-average complete active self-consistent field (CASSCF) theory with full valence electrons. Dynamics results show that the H(3)O(+) DR is ultrafast and yields diversity of products. Product branching fractions are predicted to be 0.660 for (OH + 2H), 0.230 for (H(2)O + H), 0.108 for (OH + H(2)), and 0.002 for (O + H + H(2)), which are in excellent agreement with the heavy-ion storage ring experimental results. Kinetic energies of the eliminated hydrogen atoms are large and show a bimodal distribution.

Theoretical study of the reaction of CH3 with HOCO radicals.

The reaction of HOCO radicals with CH3 radicals is examined using the coupled cluster method to locate and optimize the critical points on the ground-state potential energy surface. The results show that the CH3 + HOCO reaction can produce both the H2O + CH2CO and the CH4 + CO2 products through acetic acid and enediol intermediates. Direct ab initio dynamics calculations determine the thermal rate coefficients to be k(T/K) = 3.24 x 10(-11)T0.1024 in cm3 x molec(-1) x s(-1) at T