A comparison of semiclassical wavepacket,exact quantum and classical trajectory calculations for diffractive atom/surface scattering

Comparison of a time-dependent wavepacket approach with classical trajectory and exact coupled channel calculations is presented for He/LiF(001) diffractive scattering. The range of validity of the wavepacket method is examined in detail and attention is paid to the dependence of the results on the width of the packets. At higher energies, where the wavepacket method will be almost exact, a comparison is made with the results of classical trajectories. It is shown that, aside from a correct incorporation of quantum effects, the wavepacket approach can have a considerable computational advantage over a pure classical calculation.     

Importance of anharmonicity, recrossing effects, and quantum mechanical tunneling in transition state theory with semiclassical tunneling. A test case: the H2 + Cl hydrogen abstraction reaction

The hydrogen abstraction reaction from H2 by the Cl atom is studied by means of the variational transition state theory with semiclassical tunneling coefficients on the BW2 potential energy surface. Vibrational anharmonicity and coupling between the bending modes are taken into account. The occurrence of trajectories that recross the transition state is estimated by means of the canonical unified statistical method and by classical trajectories calculations. Different semiclassical methods for tunneling calculations are tested. Our results show that anharmonicity has a small but nonnegligible effect on the thermal rate constants, recrossing can be neglected, and tunneling is adequately described by the least-action approximation, and less successfully by the large-curvature version 3 approximation. However, the large-curvature version 4 and small-curvature approximations lead to a severe underestimation of tunneling. Thermal rate constants calculated using transition state theory including anharmonicity and tunneling agree very well with accurate quantal thermal rate constants over a wide temperature range, although the improvement over the harmonic transition state theory with the microcanonically optimized semiclassical tunneling approximation (based on version 3 of the large-curvature tunneling method) used in a previous study of this reaction is only marginal.     

Dynamical symmetries of a hydrogen atom in a circularly polarized electromagnetic field

We study dynamical symmetries for a hydrogen atom exposed to a circularly polarized electromagnetic wave in the weak field limit.     

Absorption and fluorescence of PRODAN in phospholipid bilayers: a combined quantum mechanics and classical molecular dynamics study

Absorption and fluorescence spectra of PRODAN (6-propionyl-2-dimethylaminonaphthalene) were studied by means of the time-dependent density functional theory and the algebraic diagrammatic construction method. The influence of environment, a phosphatidylcholine lipid bilayer and water, was taken into account employing a combination of quantum chemical calculations with empirical force-field molecular dynamics simulations. Additionally, experimental absorption and emission spectra of PRODAN were measured in cyclohexane, water, and lipid vesicles. Both planar and twisted configurations of the first excited state of PRODAN were taken into account. The twisted structure is stabilized in both water and a lipid bilayer, and should be considered as an emitting state in polar environments. Orientation of the excited dye in the lipid bilayer significantly depends on configuration. In the bilayer, the fluorescence spectrum can be regarded as a combination of emission from both planar and twisted structures.     

Environmental effects on the spectroscopic properties of gallic acid: a combined classical and quantum mechanical study

The solvation of gallic acid (in water and acetonitrile) is studied by means of its spectroscopic properties. IR, UV, and NMR spectra are predicted by using various solvation models obtained in terms of both purely classical and density functional approaches. Comparison with experiments is used to validate solvation models. Hydrogen-bond and long-range (or bulk) effects are evaluated by comparing different solvation models. A continuum-only approach, a purely discrete, and a mixed continuum/discrete approach based on quantum-mechanical and classical molecular-dynamics solute-solvent clusters are tested.     

Stabilization of boron clusters with classical fullerene structures by combined doping effect: a quantum chemical study

A combined approach (endohedral doping and exohedral environment) to stabilization of boron clusters with classical fullerene structures has been studied. The boron clusters with classical fullerene structures are stable when heteroatomic part of the complex (endohedral atom and exohedral environment) donates in total 18 electrons to the composite system, stability of which depends on the coordination capabilities and donor ability of the endohedral and surrounding atoms. The most effective stabilization is achieved in the case of the endohedral transition metals atoms, whereas the most effective environment is given by the lithium surrounding.

     

Comparison of quantum mechanical and semiclassical cross sections for rotational excitation of hydrogen

Partial and integral cross sections for rotational transitions in He  H₂ have been calculated by a semiclassical coupled states method. The agreement with corresponding quantum mechanical values is excellent.     

Inelastic collisions of molecular hydrogen: a comparison of results from quantum and classical mechanics

Full-dimensional quantum and classical calculations have been carried out for inelastic (nonreactive) energy transfer in H2+H2 on the ab initio potential energy surface of Boothroyd et al. [J. Chem. Phys. 116, 666 (2002)]. State-to-state cross sections are determined and compared for transitions from H2(0,j(ab))+H2(1,j(cd)). While there is excellent agreement for transitions involving small Deltaj, for larger Deltaj and for vibrational relaxation, significant differences are observed which exhibit no systematic trends. Reasons for this disagreement are discussed.     

YinYang atom: a simple combined ab initio quantum mechanical molecular mechanical model

A simple interface is proposed for combined quantum mechanical (QM) molecular mechanical (MM) calculations for the systems where the QM and MM regions are connected through covalent bonds. Within this model, the atom that connects the two regions, called YinYang atom here, serves as an ordinary MM atom to other MM atoms and as a hydrogen-like atom to other QM atoms. Only one new empirical parameter is introduced to adjust the length of the connecting bond and is calibrated with the molecule propanol. This model is tested with the computation of equilibrium geometries and protonation energies for dozens of molecules. Special attention is paid on the influence of MM point charges on optimized geometry and protonation energy, and it is found that it is important to maintain local charge-neutrality in the MM region in order for the accurate calculation of the protonation and deprotonation energies. Overall the simple YinYang atom model yields comparable results to some other QM/MM models.     

氟原子与氢分子共线反应几率的量子散射计算

基于最新的6SEC势能面,用邓从豪等提出的LCAC-SW方法计算得到了共线反应F+H~2(v=0)→HF(v')+H的态-态反应几率,计算结果准确地反映出势能面的特点,进一步证明LCAC-SW方法是一成功的量子散射方法。     

Heavy atom tunneling in chemical reactions: study of H + LiF collisions

The H+LiF(X (1)sigma(+),upsilon=0-2,j=0)-->HF(X (1)sigma(+),upsilon',j')+Li(2S) bimolecular process is investigated by means of quantum scattering calculations on the chemically accurate X 2A' LiHF potential energy surface of Aguado et al. [A. Aguado, M. Paniagua, C. Sanz, and J. Roncero, J. Chem. Phys. 119, 10088 (2003)]. Calculations have been performed for zero total angular momentum for translational energies from 10(-7) to 10(-1) eV. Initial-state selected reaction probabilities and cross sections are characterized by resonances originating from the decay of metastable states of the H...F-Li and Li...F-H van der Waals complexes. Extensive assignment of the resonances has been carried out by performing quasibound states calculations in the entrance and exit channel wells. Chemical reactivity is found to be significantly enhanced by vibrational excitation at low temperatures, although reactivity appears much less favorable than nonreactive processes due to the inefficient tunneling of the relatively heavy fluorine atom strongly bound in van der Waals complexes.     

Intramolecular hydrogen atom tunneling in 2-chlorobenzoic acid studied by low-temperature matrix-isolation infrared spectroscopy

Intramolecular hydrogen atom tunneling in 2-chlorobenzoic acid has been investigated by low-temperature matrix-isolation infrared spectroscopy with the aid of density functional theory calculation. Infrared spectra of two relatively stable syn isomers, SC and ST, were observed in argon and xenon matrixes. When the matrix samples were annealed after deposition, the isomerization from ST to SC occurred around the benzene-carboxyl bond. Two less stable anti isomers, AT, which has an OH...Cl intramolecular hydrogen bond, and AC, which has no OH...Cl bond, were produced from SC and ST upon UV irradiation. When the matrix samples were kept in the dark after UV irradiation, AT and AC changed to ST and SC, respectively, by spontaneous isomerization around the C-O axis in the carboxyl group. The rate constants of isomerization, AT --> ST, in a Xe matrix were estimated from the absorbance changes at various matrix temperatures. The rate constants showed a drastic decrease in deuteration of the hydrogen atom of the carboxyl group. The relationship between the rate constants and the matrix temperature did not follow the Arrhenius law. These findings lead to the conclusion that the isomerization of AT --> ST and AC --> SC in low-temperature rare-gas matrixes proceeds through intramolecular hydrogen atom tunneling.     

Hydrogen atom transfer reactions of a ruthenium imidazole complex: hydrogen tunneling and the applicability of the Marcus cross relation

The reaction of Ru(II)(acac)2(py-imH) (Ru(II)imH) with TEMPO(*) (2,2,6,6-tetramethylpiperidine-1-oxyl radical) in MeCN quantitatively gives Ru(III)(acac)2(py-im) (Ru(III)im) and the hydroxylamine TEMPO-H by transfer of H(*) (H(+) + e(-)) (acac = 2,4-pentanedionato, py-imH = 2-(2'-pyridyl)imidazole). Kinetic measurements of this reaction by UV-vis stopped-flow techniques indicate a bimolecular rate constant k(3H) = 1400 +/- 100 M(-1) s(-1) at 298 K. The reaction proceeds via a concerted hydrogen atom transfer (HAT) mechanism, as shown by ruling out the stepwise pathways of initial proton or electron transfer due to their very unfavorable thermochemistry (Delta G(o)). Deuterium transfer from Ru(II)(acac)2(py-imD) (Ru(II)imD) to TEMPO(*) is surprisingly much slower at k(3D) = 60 +/- 7 M(-1) s(-1), with k(3H)/k(3D) = 23 +/- 3 at 298 K. Temperature-dependent measurements of this deuterium kinetic isotope effect (KIE) show a large difference between the apparent activation energies, E(a3D) - E(a3H) = 1.9 +/- 0.8 kcal mol(-1). The large k(3H)/k(3D) and DeltaE(a) values appear to be greater than the semiclassical limits and thus suggest a tunneling mechanism. The self-exchange HAT reaction between Ru(II)imH and Ru(III)im, measured by (1)H NMR line broadening, occurs with k(4H) = (3.2 +/- 0.3) x 10(5) M(-1) s(-1) at 298 K and k(4H)/k(4D) = 1.5 +/- 0.2. Despite the small KIE, tunneling is suggested by the ratio of Arrhenius pre-exponential factors, log(A(4H)/A(4D)) = -0.5 +/- 0.3. These data provide a test of the applicability of the Marcus cross relation for H and D transfers, over a range of temperatures, for a reaction that involves substantial tunneling. The cross relation calculates rate constants for Ru(II)imH(D) + TEMPO(*) that are greater than those observed: k(3H,calc)/k(3H) = 31 +/- 4 and k(3D,calc)/k(3D) = 140 +/- 20 at 298 K. In these rate constants and in the activation parameters, there is a better agreement with the Marcus cross relation for H than for D transfer, despite the greater prevalence of tunneling for H. The cross relation does not explicitly include tunneling, so close agreement should not be expected. In light of these results, the strengths and weaknesses of applying the cross relation to HAT reactions are discussed.     

Simulation of tunneling in enzyme catalysis by combining a biased propagation approach and the quantum classical path method: application to lipoxygenase

The ability of using wave function propagation approaches to simulate isotope effects in enzymes is explored, focusing on the large H/D kinetic isotope effect of soybean lipoxygenase-1 (SLO-1). The H/D kinetic isotope effect (KIE) is calculated as the ratio of the rate constants for hydrogen and deuterium transfer. The rate constants are calculated from the time course of the H and D nuclear wave functions. The propagations are done using one-dimensional proton potentials generated as sections from the full multidimensional surface of the reacting system in the protein. The sections are obtained during a classical empirical valence bond (EVB) molecular dynamics simulation of SLO-1. Since the propagations require an extremely long time for treating realistic activation barriers, it is essential to use an effective biasing approach. Thus, we develop here an approach that uses the classical quantum path (QCP) method to evaluate the quantum free energy change associated with the biasing potential. This approach provides an interesting alternative to full QCP simulations and to other current approaches for simulating isotope effects in proteins. In particular, this approach can be used to evaluate the quantum mechanical transmission factor or other dynamical effects, while still obtaining reliable quantized activation free energies due to the QCP correction.     

Near-infrared radiation induced conformational change and hydrogen atom tunneling of 2-chloropropionic acid in low-temperature Ar matrix

Former assignments of the matrix-isolation infrared (MI-IR) spectrum of 2-chloropropionic acid are revised with the help of near-infrared (NIR) laser irradiation induced change in conformer ratios. This method allows not only the unambiguous assignment of each band in the MI-IR spectrum to the two trans (Z) and the cis (E) conformers but also the assignment of the spectral bands to different matrix sites. The tunneling decay of the higher-energy cis conformer prepared from both trans conformers in different sites is also investigated. It is shown that the tunneling decay time is very sensitive to the matrix site, especially if the in situ prepared high-energy conformer has a strained geometry in the matrix cage. The analysis shows that the kinetics of some cis → trans back conversion processes cannot be fitted by a single exponential decay. The possible reasons of this observation are examined and discussed. The present and former results clearly show that, in addition to tunneling processes, the decay rates strongly depend on solid-state effects. Therefore, simple theoretical predictions of decay rates, which do not take into account the solid-state effects, can only be compared to experimental observations only if experimentally proven that these effects do not significantly affect the experimentally measured tunneling rates.     

Overtone-induced dissociation of hydrogen peroxide: A classical trajectory study

The overtone-induced dissociation of hydrogen peroxide is studied by classical trajectories. The dissociation lifetimes are found to be on the order of picoseconds and comparable with statistical lifetimes, although energy redistribution is not complete. We conclude that overtone vibrational excitation features are negligibly lifetime-broadened by reaction.     

Intermolecular hydrogen atom abstraction from the fluoranil triplet excited state: a Raman study

Intermolecular hydrogen abstraction reaction mechanisms in photoexcited ketones have traditionally been studied using time resolved absorption spectroscopy. Another approach is presented involving time resolved resonance Raman spectroscopy to study such reactions, using the fluoranil/isopropanol system as an example. It has been shown that vibrational spectra can be recorded starting from the triplet excited state to the product state (radical anion) via the intermediate state, which is the ketyl radical. Thus, it is demonstrated that following the reaction evolution in terms of structural (vibrational modes) details would prove to be useful not only for mechanistic investigation but also for structure-reactivity correlations in photoexcited systems.     

Rotational excitation of water by hydrogen molecules: comparison of results from classical and quantum mechanics

Quasiclassical trajectory calculations are carried out for rotational excitation of water by hydrogen molecules. State-to-state rate coefficients are determined at 100 K and are compared to available quantum results. A good agreement between classical and quantum rates is observed for downward transitions, with an average accuracy of classical results better than a factor of 2. It is thus found that the ambiguities described by Faure and Wiesenfeld [J. Chem. Phys. 121, 6771 (2004)] can be solved in the particular case of waterlike asymmetric-top molecules.     

Spectroscopic and quantum chemical study of cyclopropylmethylphosphine, a candidate for intramolecular hydrogen bonding

The properties of the novel compound cyclopropylmethylphosphine (C3H5CH2PH2) have been investigated by means of Stark-modulation microwave spectroscopy and high-level quantum chemical calculations. Spectra attributable to the three conformers of the molecule with a synclinal arrangement of the H-C-C-P atoms were recorded and assigned. The experimental rotational constants obtained for these conformers were found to be in good agreement with those generated by ab initio geometry optimizations at the MP2/aug-cc-pVTZ level of theory. An estimate of the relative energies of the three conformers with observable spectra, by means of relative intensity measurements, compared favorably with the results of G3 energy calculations performed for the molecule. In addition to the observation of ground-state rotational spectra for three conformers, spectra belonging to a number of vibrationally excited states were also assigned with the aid of radio frequency microwave double-resonance experiments. A tentative assignment of these excited-state spectra was proposed by appealing to the results of density functional theory vibrational frequency calculations performed at the B3LYP/6-311++(3df,2pd) level. The energetically preferred conformer of the molecule allowed a close approach between a hydrogen atom belonging to the phosphino group and the edge of the cyclopropyl ring. The possibility of the formation of an intramolecular hydrogen bond to electron density associated with so-called banana bonds is discussed.