Meredith J. T. Jordan

Polyatomic molecular potential energy surfaces by interpolation in local internal coordinates

We present a method for expressing a potential energy surface (PES) for polyatomic molecules as an interpolation of local Taylor expansions in internal coordinates. This approach extends and replaces an earlier method which was only directly applicable to molecules of no more than four atoms. In general, the local Taylor expansions are derived from ab initio quantum calculations. Here, the methodology is evaluated by comparison with an analytic surface for the reactions H+CH4⇌H2+CH3. Approximately 1000–1300 data points are required for an accurate 12-dimensional surface which describes both forward and backward reactions, at the energy studied.

Convergence of molecular potential energy surfaces by interpolation: Application to the OH+H2→H2O+H reaction

A recently proposed scheme for interpolating and iteratively improving molecular potential energy surfaces [Ischtwan and Collins, J. Chem. Phys. 100, 8080 (1994)] is evaluated by comparison with an analytic surface for the OH+H2→H2O+H reaction. An improvement in the procedure for constructing the potential surface is suggested and implemented. The most efficient means of converging the surface is determined. It is found that the probability of reaction, for example, may be accurately calculated using of the order of 200–400 data points to define the potential energy surface.

Roaming is the dominant mechanism for molecular products in acetaldehyde photodissociation

Reaction pathways that bypass the conventional saddle-point transition state (TS) are of considerable interest and importance. An example of such a pathway, termed “roaming,” has been described in the photodissociation of H2CO. In a combined experimental and theoretical study, we show that roaming pathways are important in the 308-nm photodissociation of CH3CHO to CH4 + CO. The CH4 product is found to have extreme vibrational excitation, with the vibrational distribution peaked at ≈95% of the total available energy. Quasiclassical trajectory calculations on full-dimensional potential energy surfaces reproduce these results and are used to infer that the major route to CH4 + CO products is via a roaming pathway where a CH3 fragment abstracts an H from HCO. The conventional saddle-point TS pathway to CH4 + CO formation plays only a minor role. H-atom roaming is also observed, but this is also a minor pathway. The dominance of the CH3 roaming mechanism is attributed to the fact that the CH3 + HCO radical asymptote and the TS saddle-point barrier to CH4 + CO are nearly isoenergetic. Roaming dynamics are therefore not restricted to small molecules such as H2CO, nor are they limited to H atoms being the roaming fragment. The observed dominance of the roaming mechanism over the conventional TS mechanism presents a significant challenge to current reaction rate theory.

Classical trajectory studies of the reaction CH4+H→CH3+H2

Trajectory data are reported for the reaction CH4+H→CH3+H2, designed to provide information that can be used to test approximate quantitative theories for the dynamics of abstraction reactions. A potential function was devised which properly reflects the nuclear permutation symmetry of the process. Microscopic reaction rate coefficients were obtained as functions of fixed rotational and vibrational energy, and of the angular momentum. The data indicated significant uncoupling between the various modes although, at a minimum, the symmetric stretch is directly coupled to the reaction coordinate at the transition state. The data were used to test the assumption that the total angular momentum, J, may be approximated by the orbital angular momentum, L. L is approximately conserved from the reactant to the saddle point configuration in reactive and nonreactive collisions and may be well approximated by J. The angular momentum about the long axis of the reacting system (equivalent to the K quantum number) is no...

The utility of higher order derivatives in constructing molecular potential energy surfaces by interpolation

In this paper we evaluate the use of higher order derivatives in the construction of an interpolated potential energy surface for the OH+H2→H2O+H reaction. The surface construction involves interpolating between local Taylor expansions about a set of known data points. We examine the use of first, second, third, and fourth order Taylor expansions in the interpolation scheme. The convergence of the various interpolated surfaces is evaluated in terms of the probability of reaction. We conclude that first order Taylor expansions (and by implication zeroth order expansions) are not suitable for constructing potential energy surfaces for reactive systems. We also conclude that it is inefficient to use fourth order derivatives. The factors differentiating between second and third order Taylor expansions are less clear. Although third order surfaces require substantially fewer data points to converge than second order surfaces, this faster convergence does not offset the large cost incurred in calculating numeri...

Molecular potential energy surfaces by interpolation in Cartesian coordinates

We present a new method for expressing a molecular potential energy surface (PES) as an interpolation of local Taylor expansions. By using only Cartesian coordinates for the atomic positions, this method avoids redundancy problems associated with the use of internal coordinates. The correct translation, rotation, inversion, and permutation invariance are incorporated in the PES via the interpolation method itself. The method is most readily employed for bound molecules or clusters and is demonstrated by application to the vibrational motion of acetylene.

An interpolated unrestricted Hartree–Fock potential energy surface for the OH+H2→H2O+H reaction

In this paper we demonstrate, at the UHF/6‐311++G(d,p) level of theory, the practical feasibility of using ab initio quantum chemical calculations to generate a molecular potential energy surface (PES) for the OH+H2→H2O+H reaction using our previously suggested interpolation and iteration schemes. The successful, and almost completely automated, merger of the PES algorithm and quantum chemical calculations involves a number of significant practical problems, the solutions of which are presented in detail. The convergence of the interpolated potential surface was monitored in terms of reaction probability and we find that the surface converges once the energy, gradient and Hessian have been calculated at approximately 350 geometries. We also find that, although the initial geometries used consisted only of points along a reaction path for the OH+H2→H2O+H reaction, the potential energy surface iteration process rapidly adds information about other, energetically accessible, reaction channels.

What a difference a decade makes: progress in ab initio studies of the hydrogen bond

Abstract This article provides a summary of our studies of hydrogen-bonded complexes during the decade of the 90s. These studies began with systematic investigations of the methodological dependence of the computed structures and binding energies of these complexes. The MP2/6-31+G(d,p) level of theory was identified as the minimum level required to obtain reliable structures, while reliable energetics required larger polarized split-valence basis sets that include diffuse functions. While the experimental frequency shift of the A–H stretching band upon formation of an A–H–B hydrogen bond could also be reproduced at MP2/6-31+G(d,p) for a variety of hydrogen-bonded complexes, significant discrepancies were observed for others, including complexes of HCl and HBr with ammonia, trimethylamine, and 4-substituted pyridines. Resolving these discrepancies became the primary focus of our work, and redefined our research efforts. We solved a model two-dimensional nuclear Schrodinger equation to obtain anharmonic dimer- and proton-stretching frequencies, modeled matrix effects with external electric fields, and characterized hydrogen bond types as traditional, proton-shared, and ion-pair. We were able to resolve the observed discrepancies between theory and experiment, and explain the rather disparate effects of matrices on the IR spectra of closely related complexes. We also initiated studies of the NMR properties of the chemical shift of the hydrogen-bonded proton, and the A–B spin–spin coupling constant across the A–H–B hydrogen bond. We demonstrated the dominance of the Fermi-contact term for determining coupling constants in complexes with N–H–N, N–H–O, O–H–O, and Cl–H–N hydrogen bonds, and the distance dependence of this term. We also showed that the IR anharmonic proton-stretching frequency and the NMR spin–spin coupling constant are spectroscopic fingerprints of hydrogen bond type, which provide information about intermolecular distances in hydrogen-bonded complexes.

Photo-tautomerization of Acetaldehyde to Vinyl Alcohol: A Potential Route to Tropospheric Acids

Enols in the Atmosphere? Keto/enol tautomerization (HC−C=O→C=C−OH) plays a central role in the chemistry of carbonyl compounds in a solution in which solvent and catalytic acids or bases can facilitate the proton transfer from C to O and back again. In contrast, analyses of atmospheric chemistry tend to exclude enol structure, on the assumption that tautomerization does not proceed regularly in gas phase. Andrews et al. (p. 1203, published online 16 August) used isotopic labeling to probe the photoisomerization pathway of gaseous acetaldehyde in the lab and discovered evidence for an enol. Subsequent modeling indicates that photogenerated enols could build up sufficiently in the troposphere to account for previously puzzling observations of organic acids in the atmosphere. Enol tautomers may play a bigger role in atmospheric chemistry than previously suspected. Current atmospheric models underestimate the production of organic acids in the troposphere. We report a detailed kinetic model of the photochemistry of acetaldehyde (ethanal) under tropospheric conditions. The rate constants are benchmarked to collision-free experiments, where extensive photo-isomerization is observed upon irradiation with actinic ultraviolet radiation (310 to 330 nanometers). The model quantitatively reproduces the experiments and shows unequivocally that keto-enol photo-tautomerization, forming vinyl alcohol (ethenol), is the crucial first step. When collisions at atmospheric pressure are included, the model quantitatively reproduces previously reported quantum yields for photodissociation at all pressures and wavelengths. The model also predicts that 21 ± 4% of the initially excited acetaldehyde forms stable vinyl alcohol, a known precursor to organic acid formation, which may help to account for the production of organic acids in the troposphere.

Vibrational spectroscopy of the hydrogen bond: An ab initio quantum-chemical perspective

The hydrogen bond has long been recognized as an important type of intermolecular interaction. Its infrared (IR) spectroscopic signature is the shift to lower frequency and the increase in intensity of the A-H stretching band upon formation of the A-H…B hydrogen bond. Ab initio calculations carried out with an appropriate wavefunction model and basis set, and using the harmonic approximation, can reasonably reproduce the shift of the A-H stretching band upon hydrogen bonding, if the equilibrium structure exists in a relatively deep potential well on the surface, so that both the V=0 and the V=1 vibrational states of the proton-stretching mode are confined within this well. However, if the equilibrium structure is found in a region of the surface which is broad and relatively flat, or if a second region of the surface can be accessed in either the V=0 or the V=1 vibrational state of the proton-stretching mode, then the harmonic approximation fails to describe the anharmonicity inherent in the surface. For ...