PubMed Health⌕ Search

Biomedical subjects

Joel M Bowman

Publications and source records attributed to Joel M Bowman.

At least 19 recordsLinked to original sources
PubMed2006

Potential energy surface and MULTIMODE vibrational analysis of C2H3+.

A full dimensional, ab initio-based semiglobal potential energy surface for C(2)H(3) (+) is reported. The ab initio electronic energies for this molecule are calculated using the spin-restricted, coupled cluster method restricted to single and double excitations with triples corrections [RCCSD(T)]. The RCCSD(T) method is used with the correlation-consistent polarized valence triple-zeta basis augmented with diffuse functions (aug-cc-pVTZ). The ab initio potential energy surface is represented by a many-body (cluster) expansion, each term of which uses functions that are fully invariant under permutations of like nuclei. The fitted potential energy surface is validated by comparing normal mode frequencies at the global minimum and secondary minimum with previous and new direct ab initio frequencies. The potential surface is used in vibrational analysis using the "single-reference" and "reaction-path" versions of the code MULTIMODE.

Journal Article↗
PubMed2006

The calculated infrared spectrum of Cl-H2O using a new full dimensional ab initio potential surface and dipole moment surface.

We report a full dimensional, ab initio-based global potential energy surface (PES) and dipole moment surface for Cl-H2O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method [CCSD(T)] with a moderately large basis (aug-cc-pVTZ). Vibrational energies and wave functions are accurately obtained using MULTIMODE. The wave function and dipole moment surface are used to calculate and analyze the pure infrared spectrum at 0 K which is compared with experiment. Vibrational energies and the infrared spectra for DOD and HOD/DOH are also presented.

Journal Article↗
PubMed2006

Quasiclassical trajectory study of the reaction H+CH4(nu3 = 0,1)-->CH3+H2 using a new ab initio potential energy surface.

Detailed quasiclassical trajectory calculations of the reaction H+CH4(nu3 = 0,1)-->CH3 + H2 using a slightly updated version of a recent ab initio-based CH5 potential energy surface [X. Zhang et al., J. Chem. Phys. 124, 021104 (2006)] are reported. The reaction cross sections are calculated at initial relative translational energies of 1.52, 1.85, and 2.20 eV in order to make direct comparison with experiment. The relative reaction cross section enhancement ratio due to the excitation of the C-H antisymmetric stretch varies from 2.2 to 3.0 over this energy range, in good agreement with the experimental result of 3.0 +/- 1.5 [J. P. Camden et al., J. Chem. Phys. 123, 134301 (2005)]. The laboratory-frame speed and center-of-mass angular distributions of CH3 are calculated as are the vibrational and rotational distributions of H2 and CH3. We confirm that this reaction occurs with a combination of stripping and rebound mechanisms by presenting the impact parameter dependence of these distributions and also by direct examination of trajectories.

Journal Article↗
PubMed2006

Quasiclassical trajectory calculations of the reaction C+C2H2-->l-C3H, c-C3H+H, C3+H2 using full-dimensional triplet and singlet potential energy surfaces.

Full-dimensional, density functional theory (B3LYP/6-311g(d,p))-based potential energy surfaces (PESs) are reported and used in quasi-classical calculations of the reaction of C with C(2)H(2). For the triplet case, the PES spans the region of the reactants, the complex region (with numerous minima and saddle points) and the products, linear(l)-C(3)H+H, cyclic(c)-C(3)H+H and c-(3)C(3)+H(2). For the singlet case, the PES describes the complex region and products l-C(3)H+H, c-C(3)H+H and l-(1)C(3)+H(2). The PESs are invariant under permutation of like nuclei and are fit to tens of thousands of electronic energies. Energies and harmonic frequencies of the PESs agree well the DFT ones for all stationary points and for the reactant and the products. Dynamics calculations on the triplet PES find both l-C(3)H and c-C(3)H products, with l-C(3)H being dominant at the energies considered. Limited unimolecular reaction dynamics on the singlet PES find both products in comparable amounts as well as the C(3)+H(2) product.

Journal Article↗
PubMed2006

The roaming atom pathway in formaldehyde decomposition.

We present a detailed experimental and theoretical investigation of formaldehyde photodissociation to H(2) and CO following excitation to the 2(1)4(1) and 2(1)4(3) transitions in S(1). The CO velocity distributions were obtained using dc slice imaging of single CO rotational states (v=0, j(CO)=5-45). These high-resolution measurements reveal the correlated internal state distribution in the H(2) cofragments. The results show that rotationally hot CO (j(CO) approximately 45) is produced in conjunction with vibrationally "cold" H(2) fragments (v=0-5): these products are formed through the well-known skewed transition state and described in detail in the accompanying paper. After excitation of formaldehyde above the threshold for the radical channel (H(2)CO-->H+HCO) we also find formation of rotationally cold CO (j(CO)=5-28) correlated to highly vibrationally excited H(2) (v=6-8). These products are formed through a novel mechanism that involves near dissociation followed by intramolecular H abstraction [D. Townsend et al., Science 306, 1158 (2004)], and that avoids the region of the transition state entirely. The dynamics of this "roaming" mechanism are the focus of this paper. The correlations between the vibrational states of H(2) and rotational states of CO formed following excitation on the 2(1)4(3) transition allow us to determine the relative contribution to molecular products from the roaming atom channel versus the conventional molecular channel.

Carbon Monoxide↗
PubMed2006

Termolecular kinetics for the Mu + CO + M recombination reaction: A unique test of quantum rate theory.

The room-temperature termolecular rate constants, k0, for the Mu + CO + M<==>MuCO + M (M = He, N2, Ar) recombination reaction have been measured by the muSR technique, and are reported for moderator gas pressures of up to approximately 200 bar (densities less, similar 0.4 x 10(22) molec cm(-3)). The experimental relaxation rates reveal an unusual signature, in being dominated by the electron spin-rotation interaction in the MuCO radical that is formed in the addition step. In N2 moderator, k0 = 1.2+/-0.1 x 10(-34) cm(6) s(-1), only about 30% higher than found in Ar or He. The experimental results are compared with theoretical calculations carried out on the Werner-Keller-Schinke (WKS) surface [Keller et al., J. Chem. Phys. 105, 4983 (1996)], within the framework of the isolated resonance model (IRM). The positions and lifetimes of resonance states are obtained by solving the complex Hamiltonian for the nonrotating MuCO system, using an L2 method, with an absorbing potential in the asymptotic region. Accurate values of the vibrational bound and resonance states of MuCO reveal unprecedented isotope effects in comparisons with HCO, due to the remarkable effect of replacing H by the very light Mu atom (m(Mu) approximately (1/9)m(H)). Due to its pronounced zero-point energy shift, there are only two (J = 0) bound states in MuCO. Contributions from nonzero J states to the termolecular rate constants are evaluated through the J-shifting approximation, with rotational constants evaluated at the potential minimum. The value of the important A constant (181 cm(-1)) used in this approximation was supported by accurate J = K = 1 calculations, from which A = 180 cm(-1) was obtained by numerical evaluation. The calculations presented here, with a "weak collision factor" beta c = 0.001, indicative of the very sparse density of MuCO states, give a very good account of both the magnitude and pressure dependence of the experimental rates, but only when the fact that the two initially bound (J = 0) states become resonances for J > 0 is taken into account. This is the first time in IRM calculations of atom-molecule recombination reactions where J not equal to 0 states have proven to be so important, thus providing a truly unique test of quantum rate theory.

Journal Article↗
PubMed2006

Zero-point energy constraint in quasi-classical trajectory calculations.

A method to constrain the zero-point energy in quasi-classical trajectory calculations is proposed and applied to the Henon-Heiles system. The main idea of this method is to smoothly eliminate the coupling terms in the Hamiltonian as the energy of any mode falls below a specified value.

Journal Article↗
PubMed2006

Ab initio calculation of the low-lying vibrational states of C2H2(A) in full dimensionality.

We report full-dimensional calculations of vibrational energies of trans-C2H2(A) using the code MULTIMODE and with a full-dimensional potential energy surface obtained by fitting singles and doubles coupled-cluster equations-of-motion (EOM-CCSD) energies using a [3s 2p 1d] atomic natural orbital basis. The EOM-CCSD calculations were done with the code "ACES II". We compare the properties of the potential surface to previous calculations at the trans minimum and also compare the vibrational energies to experimental ones.

Journal Article↗
PubMed2006

The calculated infrared spectrum of Cl- H2O using a full dimensional ab initio potential surface and dipole moment surface.

We report a full dimensional ab initio-based global potential energy surface (PES) and dipole moment surface (DMS) for Cl(-)H(2)O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method (CCSD(T)) with a moderately large basis (aug-cc-pVTZ). The infrared spectrum and vibrational dynamics are reported and compared to experiment. These results are analyzed by examination of wave function and the dipole surface.

Journal Article↗
PubMed2006

Deuteration effects on the structure and infrared spectrum of CH5(+).

The CH5+ molecular ion is well-known for its large amplitude motions that lead to complete scrambling of the hydrogen atoms, even in the vibrational ground state. Experiments have been reported that probe the consequences of these large amplitude motions. We recently reported that quantum zero-point effects partially quench the scrambling when CH5+ is partially deuterated. Here, the consequences of this quantum localization are investigated through calculations of the low-resolution spectra of CH4D+, CHD4+, and CD5+. The spectra are obtained by convoluting stick spectra, evaluated for individual stationary points on an ab initio potential surface, multiplying them by Diffusion Monte Carlo ground state density at that stationary point, and taking the sum. The CH/D stretch regions of CH4D+ and CD5+ are red-shifted relative to CH5+, while the overall shape of the envelope remains unaffected. In contrast, for CHD4+, the shape of the spectral envelope in the CH/D stretch region differs from the other three isotopologs. These signatures of the quantum localization of the deuterium on the spectra are discussed.

Journal Article↗
PubMed2006

Vibrational analysis of the H5O2+ infrared spectrum using molecular and driven molecular dynamics.

Standard molecular and driven molecular dynamics are used to analyze prominent spectral features in the H5O2+ infrared spectrum. In the driven method, the molecular Hamiltonian is augmented with a time-dependent term, mu x epsilon(0) sin(omegat), where mu is the dipole moment of H5O2+, epsilon0 is the electric field, and omega is the frequency. The magnitude of the electric field determines whether the driving is mild (the harmonic limit) or strong (anharmonic motion and mode coupling). We analyze the spectrum in the wavenumber range from 600 to 1900 cm(-1), where recent experimental measurements are available for H5O2+. On the basis of the simulations, we have assigned the broad feature around 1000 cm(-1) to the proton transfer coupled with the torsion motion. Intense absorption near 1780 cm(-1) is assigned to the H2O monomer bend coupled with proton transfer.

Models, Chemical↗
PubMed2006

An ab initio based global potential energy surface describing CH5+ --> CH3+ + H2.

A full-dimensional, ab initio based potential energy surface (PES) for CH(5)(+), which can describe dissociation is reported. The PES is a precise fit to 36173 coupled-cluster [CCSD(T)] calculations of electronic energies done using an aug-cc-pVTZ basis. The fit uses a polynomial basis that is invariant with respect to permutation of the five H atoms, and thus describes all 120 equivalent minima. The rms fitting error is 78.1 cm(-1) for the entire data set of energies up to 30,000 cm(-1) and a normal-mode analysis of CH(5)(+) also verifies the accuracy of the fit. Two saddle points have been located on the surface as well and compared with previous theoretical work. The PES dissociates correctly to the fragments CH(3)(+) + H(2) and the equilibrium geometry and normal-mode analyses of these fragments are also presented. Diffusion Monte Carlo calculations are done for the zero-point energies of CH(5)(+) (and some isotopologs) as well as for the separated fragments of CH(5)(+), CH(3)(+) + H(2) and those of CH(4)D(+), CH(3)(+) + HD and CH(2)D(+) + H(2). Values of D(0) are reported for these dissociations. A molecular dynamics calculation of CH(4)D(+) dissociation at one total energy is also performed to both validate the applicability of the PES for dynamics studies as well as to test a simple classical statistical prediction of the branching ratio of the dissociation products.

Journal Article↗
PubMed2006

Ab initio potential energy and dipole moment surfaces of (H2O)2.

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with Møller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).

Journal Article↗
PubMed2006

An ab initio potential surface describing abstraction and exchange for H+CH4.

We report an ab initio-based global potential energy surface for H+CH4 that describes the abstraction and exchange reactions. The PES, which is invariant with respect to any permutation of five H atoms, is a fit to 20,728 electronic energies calculated using the partially spin-restricted coupled-cluster method (RCCSD(T)) with a moderately large basis (aug-cc-pVTZ). A first set of quasiclassical trajectory calculations using this PES are reported for the H+CD4-->HD+CD3 reaction at collision energies of 0.65 and 1.52 eV and are compared to experiment and recent direct dynamics calculations done with density functional theory.

Journal Article↗
PubMed2006

Quantum deconstruction of the infrared spectrum of CH5+.

We present two quantum calculations of the infrared spectrum of protonated methane (CH5+) using full-dimensional, ab initio-based potential energy and dipole moment surfaces. The calculated spectra compare well with a low-resolution experimental spectrum except below 1000 cm(-1), where the experimental spectrum shows no absorption. The present calculations find substantial absorption features below 1000 cm(-1), in qualitative agreement with earlier classical calculations of the spectrum. The major spectral bands are analyzed in terms of the molecular motions. Of particular interest is an intense feature at 200 cm(-1), which is due to an isomerization mode that connects two equivalent minima. Very recent high-resolution jet-cooled spectra in the CH stretch region (2825 to 3050 cm(-1)) are also reported, and assignments of the band origins are made, based on the present quantum calculations.

Journal Article↗
PubMed2005

New insights on reaction dynamics from formaldehyde photodissociation.

We review the photodissociation dynamics of formaldehyde with an emphasis on recent calculations that make use of a global ab initio-based potential energy surface for the S(0) state. These calculations together with recent experiments reveal striking departures from conventional transition state theory for the formation of the molecular products H(2) + CO. The evidence for this departure is reviewed in detail by examining properties of the new potential surface and results of quasiclassical trajectory dynamics calculations using this surface. We also review very recent work on the dynamics governing the formation of radical products, H + HCO. These products can be formed on the T(1) surface as well as the S(0) one, and we present some results contrasting the dynamics on these two surfaces. This work makes use of a new semi-global ab initio-based T(1) potential energy surface.

Formaldehyde↗
PubMed2005

Quantum studies of the vibrations in H3O2- and D3O2-.

The vibrations of H3O2- and D3O2- are investigated using diffusion Monte Carlo (DMC) and vibrational configuration-interaction approaches, as implemented in the program MULTIMODE. These studies use the potential surface recently developed by Huang et al. [ J. Am. Chem. Soc. 126, 5042 (2004)]. The focus of this work is on the vibrational ground state and fundamentals which occur between 100 and 3700 cm(-1). In most cases, excellent agreement is obtained between the fundamental frequencies calculated by the two approaches. This serves to demonstrate the power of both methods for treating this very anharmonic system. Based on the results of the MULTIMODE and DMC treatments, the extent and nature of the couplings in H3O2- and D3O2- are investigated.

Journal Article↗