Rotational excitations in para-H2+para-H2 collisions: full- and reduced-dimensional quantum wave packet studies comparing different potential energy surfaces

We study the process of rotational excitation in the collisions of para-H(2) with para-H(2) by propagating wave packets with the multiconfiguration time-dependent Hartree (MCTDH) algorithm. Transition probabilities are then calculated by the method of Tannor and Weeks based on time-correlation functions. Calculations were carried out up to a total angular momentum of J=70 to compute integral cross sections up to 1.2 eV in collision energy and thermal rate coefficients from 100 to 3000 K. The process is studied on the full-dimensional potential energy surface of Boothroyd-Martin-Keogh-Peterson (BMKP) as well as on the rigid rotor surface of Diep and Johnson. We test the validity of the rigid rotor approximation by also considering two rigid rotor restrictions of the BMKP potential energy surface (PES). Additionally, we investigate a variant of the BMKP PES suggested by Pogrebnya and Clary [Chem. Phys. Lett. 363, 523 (2002)] with reduced anisotropy. We compare our results with previous theoretical data for the cross sections and with experimental data for the rate coefficients at low temperatures.     

Rovibrational energy transfer in ortho-H2+para-H2 collisions

We present the results of a full-dimensional quantum mechanical study of the rovibrational energy transfer in the collision between ortho-H2 and para-H2 in the energy range of 0.1-1.0 eV. The multiconfiguration time-dependent Hartree algorithm has been used to propagate the wave packets on the global potential energy surface by Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)] and on a modification of this surface where the short range anisotropy is reduced. State-to-state attributes such as probabilities or integral cross sections are obtained using the formalism of Tannor and Weeks [J. Chem. Phys. 98, 3884 (1993)] by Fourier transforming the correlation functions. The effect of initial rotation of the diatoms on the inelastic and de-excitation processes is investigated.     

Rotational excitation of H2O by para-H2 from an adiabatically reduced dimensional potential

Cross sections and rate coefficients for low lying rotational transitions in H(2)O colliding with para-hydrogen pH(2) are computed using an adiabatic approximation which reduces the dimensional dynamics from a 5D to a 3D problem. Calculations have been performed at the close-coupling level using the recent potential of Valiron et al. [J. Chem. Phys. 129, 134306 (2008)]. A good agreement is found between the reduced adiabatic calculations and the 5D exact calculations, with an impressive time saving and memory gain. This adiabatic reduction of dimensionality seems very promising for scattering studies involving the excitation of a heavy target molecule by a light molecular projectile.     

Inelastic collisions in para-H2: translation-rotation state-to-state rate coefficients and cross sections at low temperature and energy

We report an experimental determination of the k(00-->02) rate coefficient for inelastic H(2):H(2) collisions in the temperature range from 2 to 110 K based on Raman spectroscopy data in supersonic expansions of para-H(2). For this purpose a more accurate method for inverting the master equation of rotational populations is presented. The procedure permits us to reduce the measured k(00-->02) rate coefficient to the corresponding sigma(00-->02) cross section in the range of precollisional energy from 360 to 600 cm(-1). Numerical calculations of sigma(00-->02) carried out in the frame of the coupled channel method are also reported for different intermolecular potentials of H(2). A good agreement is found between the experimental cross section and the numerical one derived from Diep and Johnson's potential [J. Chem. Phys. 112, 4465 (2000)].     

Quantum dynamics study of H+NH3-->H2+NH2 reaction

We report in this paper a quantum dynamics study for the reaction H+NH3-->NH2+H2 on the potential energy surface of Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The quantum dynamics calculation employs the semirigid vibrating rotor target model [J. Z. H. Zhang, J. Chem. Phys. 111, 3929 (1999)] and time-dependent wave packet method to propagate the wave function. Initial state-specific reaction probabilities are obtained, and an energy correction scheme is employed to account for zero point energy changes for the neglected degrees of freedom in the dynamics treatment. Tunneling effect is observed in the energy dependency of reaction probability, similar to those found in H+CH4 reaction. The influence of rovibrational excitation on reaction probability and stereodynamical effect are investigated. Reaction rate constants from the initial ground state are calculated and are compared to those from the transition state theory and experimental measurement.     

Theoretical study of the rovibrational spectrum of H2O-H2

In this paper we report transition frequencies and line strengths computed for H(2)O-H(2) and compare with the experimental observations of [M. J. Weida and D. J. Nesbitt, J. Chem. Phys. 110, 156 (1999)]. To compute the spectra we use a symmetry adapted Lanczos algorithm and an uncoupled product basis set. Our results corroborate the assignments of Weida and Nesbitt and there is good agreement between calculated and observed transitions. Possible candidates for lines that Weida and Nesbitt were not able to assign are presented. Several other bands that may be observable are also discovered. Although all the observed bands are associated with states localized near the global potential minimum, at which H(2)O acts as proton acceptor, a state with significant amplitude near the T-shape secondary potential minimum at which H(2)O acts as proton donor is identified by examining many different probability density plots.     

Theoretical study of the complex-forming CH + H2 --> CH2 + H reaction

The complex-forming CH + H2 --> CH2 + H reaction is studied employing a recently developed global potential energy function. The reaction probability in the total angular momentum J = 0 limit is estimated with a four-atom quantum wave packet method and compared with classical trajectory and statistical theory results. The formation of complexes from different reactant internal states is also determined with wave packet calculations. While there is no barrier to reaction along the minimum energy path, we find that there are angular constraints to complex formation. Trajectory-based estimates of the low-pressure rate constants are made and compared with experimental results. We find that zero-point energy violation in the trajectories is a particularly severe problem for this reaction.     

Rovibrational states of the H2O-H2 complex: an ab initio calculation

All bound rovibrational levels of the H(2)O-H(2) dimer are calculated for total angular momentum J = 0-5 on two recent intermolecular potential surfaces reported by Valiron et al. [J. Chem. Phys. 129, 134306 (2008)] and Hodges et al. [J. Chem. Phys. 120, 710 (2004)] obtained through ab initio calculations. The method used handles correctly the large amplitude internal motions in this complex; it involves a discrete variable representation of the intermolecular distance coordinate R and a basis of coupled free rotor wave functions for the hindered internal rotations and the overall rotation of the dimer. The basis is adapted to the permutation symmetry associated with the para/ortho (p/o) nature of both H(2)O and H(2) as well as to inversion symmetry. Dimers containing oH(2) are more strongly bound than dimers with pH(2), as expected, with dissociation energies D(0) of 33.57, 36.63, 53.60, and 59.04 cm(-1)for pH(2)O-pH(2), oH(2)O-pH(2), pH(2)O-oH(2), and oH(2)O-oH(2), respectively, on the potential of Valiron et al. that corresponds to a binding energy D(e) of 235.14 cm(-1). Rovibrational wave functions are computed as well and the nature of the bound states in the four different dimer species is discussed. Converged rovibrational levels on both potentials agree well with the high-resolution spectrum reported by Weida and Nesbitt [J. Chem. Phys. 110, 156 (1999)]; the hindered internal rotor model that was used to interpret this spectrum is qualitatively correct.     

Quantum wave packet dynamics of N(2D)+H2 reaction

The quantum wave packet dynamics of the title reaction within the coupled state approximation is examined here and initial state-selected reaction probabilities, integral reaction cross sections, and thermal rate constants are reported. The ab initio potential energy surface of the electronic ground state (1(2)A(")) of the system recently reported by Ho et al. [J. Chem. Phys., 119, 3063 (2003)] is employed in this investigation. All partial wave contributions up to the total angular momentum J=55 were necessary to obtain converged integral reaction cross sections up to a collision energy of 1.0 eV. Thermal rate constants are calculated from the reaction cross sections and compared with the available theoretical and experimental results. Typical resonances formed during the course of the reaction and elucidating the insertion type mechanism for the product formation are calculated. Vibrational energy levels supported by the deep well (approximately 5.5 eV) of the 1(2)A(") potential energy surface of NH(2) are also calculated for the total angular momentum J=0. A statistical analysis of the spacing between the adjacent levels of this energy spectrum is performed and the extent of irregularity in the spectral sequence is assessed.     

Variational quantum mechanical and active database approaches to the rotational-vibrational spectroscopy of ketene, H2CCO

A variational quantum mechanical protocol is presented for the computation of rovibrational energy levels of semirigid molecules using discrete variable representation of the Eckart-Watson Hamiltonian, a complete, "exact" inclusion of the potential energy surface, and selection of a vibrational subspace. Molecular symmetry is exploited via a symmetry-adapted Lanczos algorithm. Besides symmetry labels, zeroth-order rigid-rotor and harmonic-oscillator quantum numbers are employed to characterize the computed rovibrational states. Using the computational molecular spectroscopy algorithm presented, a large number of rovibrational states, up to J = 50, of the ground electronic state of the parent isotopologue of ketene, H(2) (12)C=(12)C=(16)O, were computed and characterized. Based on 12 references, altogether 3982 measured and assigned rovibrational transitions of H(2) (12)C=(12)C=(16)O have been collected, from which 3194 were validated. These transitions form two spectroscopic networks (SN). The ortho and the para SNs contain 2489 and 705 validated transitions and 1251 and 471 validated energy levels, respectively. The computed energy levels are compared with energy levels obtained, up to J = 41, via an inversion protocol based on this collection of validated measured rovibrational transitions. The accurate inverted energy levels allow new assignments to be proposed. Some regularities and irregularities in the rovibrational spectrum of ketene are elucidated.     

Geometric phase effects in resonance-mediated scattering: H + H2(+) on its lowest triplet electronic state

We report a quantum dynamics study of H + H(2)(+) (v(0) = 0, j(0) = 0) scattering on its lowest triplet state, for J = 0 total angular momentum and total energies up to 1.85 eV. This provides a benchmark example of indirect resonance-mediated reaction in presence of a conical intersection (CI). Visualization of time-dependent wave packets shows significant "looping" around the CI, which is facilitated by long-lived H(3)(+) scattering resonances, predominant at low energies. State-to-state inelastic transition probabilities exhibit a highly oscillatory structure and pronounced geometric phase effects, which, in contrast to direct reactions, are more strongly marked at lower energies.     

Water line lists close to experimental accuracy using a spectroscopically determined potential energy surface for H2(16)O, H2(17)O, and H2(18)O

Line lists of vibration-rotation transitions for the H(2) (16)O, H(2) (17)O, and H(2) (18)O isotopologues of the water molecule are calculated, which cover the frequency region of 0-20 000 cm(-1) and with rotational states up to J=20 (J=30 for H(2) (16)O). These variational calculations are based on a new semitheoretical potential energy surface obtained by morphing a high accuracy ab initio potential using experimental energy levels. This potential reproduces the energy levels with J=0, 2, and 5 used in the fit with a standard deviation of 0.025 cm(-1). Linestrengths are obtained using an ab initio dipole moment surface. That these line lists make an excellent starting point for spectroscopic modeling and analysis of rotation-vibration spectra is demonstrated by comparison with recent measurements of Lisak and Hodges [J. Mol. Spectrosc. (unpublished)]: assignments are given for the seven unassigned transitions and the intensity of the strong lines are reproduced to with 3%. It is suggested that the present procedure may be a better route to reliable line intensities than laboratory measurements.     

Spatial delocalization in para-H2 clusters

We applied the quantum path integral Monte Carlo method for the study of (para-H)N (N = 5-33) clusters at T = 2 K, exploring static and dynamic order, which originates from the effects of zero-point energy, kinetic energy, and thermal fluctuations in quantum clusters. Information on dynamic structure was inferred from the asymptotic tails of the cage correlation function calculated from the centroid Monte Carlo trajectory. The centroid cage correlation function decays to zero for large clusters (N = 15-33), manifesting the interchange of molecules between different solvation shells, with statistically diminishing back interchange. Further evidence for the floppiness of para-hydrogen clusters emerges from the Monte Carlo evolution of the centroid of a tagged molecule, which exhibits significant changes in the list of its first and second solvation shells due to the interchange of molecules between these shells.     

Time-dependent quantum study of the kinetics of the H(2S) + FO(2II) OH(2II) + F(2P) reaction

Quantum mechanical wave packet calculations are carried out for the H((2)S) + FO((2)II) --> OH((2)II) + F((2)P) reaction on the adiabatic potential energy surface of the ground (3)A' triplet state. The state-to-state and state-to-all reaction probabilities for total angular momentum J = 0 have been calculated. The probabilities for J > 0 have been estimated from the J = 0 results by using J-shifting approximation based on a capture model. Then, the integral cross sections and initial state-selected rate constants have been calculated. The calculations show that the initial state-selected reaction probabilities are dominated by many sharp peaks. The reaction cross section does not manifest any sharp oscillations and the initial state-selected rate constants are sensitive to the temperature.     

Femtosecond study of Cu(H(2)O) dynamics

The short-time nuclear dynamics of Cu(H(2)O) is investigated using femtosecond photodetachment-photoionization spectroscopy and time-dependent quantum wave packet calculations. The Cu(H(2)O) dynamics is initiated in the electronic ground state of the complex by electron photodetachment from the Cu(-)(H(2)O) complex, where hydrogen atoms are oriented toward Cu. Several time-resolved resonant multiphoton ionization schemes are used to probe the ensuing reorientation and dissociation. Immediately following photodetachment, the neutral complex is far from its minimum energy geometry and possesses an internal energy comparable to the Cu-H(2)O dissociation energy and undergoes both large-amplitude H(2)O motion and dissociation. Dissociation is observed to occur on three distinct time scales: 0.6, 8, and 100 ps. These results are compared to the results of time-dependent J=0 wave packet calculations, propagating the initial anion vibrational wave functions on the ground-state potential of the neutral complex. An excellent agreement is obtained between the experimental results and the ionization signals derived from the calculated probability amplitudes. Related experiments and calculations are carried out on the Cu(D(2)O) complex, with results very similar to those of Cu(H(2)O).     

Quantum approaches for the insertion dynamics of the H+ + D2 and D+ + H2 reactive collisions

The H(+)+D(2) and D(+)+H(2) reactive collisions are studied using a recently proposed adiabatic potential energy surface of spectroscopic accuracy. The dynamics is studied using an exact wave packet method on the adiabatic surface at energies below the curve crossing occurring at approximately 1.5 eV above the threshold. It is found that the reaction is very well described by a statistical quantum method for a zero total angular momentum (J) as compared with the exact ones, while for higher J some discrepancies are found. For J >0 different centrifugal sudden approximations are proposed and compared with the exact and statistical quantum treatments. The usual centrifugal sudden approach fails by considering too high reaction barriers and too low reaction probabilities. A new statistically modified centrifugal sudden approach is considered which corrects these two failures to a rather good extent. It is also found that an adiabatic approximation for the helicities provides results in very good agreement with the statistical method, placing the reaction barrier properly. However, both statistical and adiabatic centrifugal treatments overestimate the reaction probabilities. The reaction cross sections thus obtained with the new approaches are in rather good agreement with the exact results. In spite of these deficiencies, the quantum statistical method is well adapted for describing the insertion dynamics, and it is then used to evaluate the differential cross sections.     

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.     

NMR studies on derivatives of heteroaromatic compounds: exceptionally small carbon-carbon couplings, 1J(CC), in 2-lithiothiophene, 2-lithio-N-methylpyrrole and 2-lithiofuran

Spin-spin carbon-carbon coupling constants across one bond and carbon proton coupling constants across one, two and three bonds have been measured for a large series of derivatives of five-membered heterocyclic compounds. This included 2-methyl and 2-lithio derivatives of furan, N-methyl pyrrole and thiophene and a series of 2-R-substituted thiophenes where R = O-t-Bu, Cl, Br, I, Si(CH3)3, MgBr and MgTh. For the long-range C-H couplings their signs have been determined in several compounds by the use of modern 2D NMR techniques, and in all the cases they have been found to be positive. A good linear dependence upon electronegativity of the substituent has been observed for 1J(CC), 2J(C2H3), and 3J(C2H4). Very small 1J(CC) couplings of ca. 30 Hz only have been found in all three lithio compounds; they belong to the smallest couplings of this type reported up to now. They are accompanied by very large and positive two-bond couplings, 2J(C2H3) of ca. +20 Hz, which in turn are the largest couplings of this type reported so far. In both cases the changes observed are interpreted in terms of the changes in the Fermi contact contribution.