Niyazi Bulut

A detailed quantum mechanical and quasiclassical trajectory study on the dynamics of the H++H2→H2+H+ exchange reaction

The H+ + H2 exchange reaction has been studied theoretically by means of a different variety of methods as an exact time independent quantum mechanical, approximate quantum wave packet, statistical quantum, and quasiclassical trajectory approaches. Total and state-to-state reaction probabilities in terms of the collision energy for different values of the total angular momentum obtained with these methods are compared. The dynamics of the reaction is extensively studied at the collision energy of E(coll)=0.44 eV. Integral and differential cross sections and opacity functions at this collision energy have been calculated. In particular, the fairly good description of the exact quantum results provided by the statistical quantum method suggests that the dynamics of the process is governed by an insertion mechanism with the formation of a long-lived collision complex.

On the dynamics of the H(+) + D2(v=0, j=0) --> HD + D(+) reaction: A comparison between theory and experiment

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.

Time-dependent wave packet and quasiclassical trajectory study of the C(P3)+OH(X Π2)→CO(X Σ1+)+H(S2) reaction at the state-to-state level

The first calculations of state-to-state reaction probabilities and product state-resolved integral cross sections at selected collision energies (0.05, 0.1, 0.5, and 1.0 eV) for the title reaction on the ab initio potential energy surface of [Zanchet et al. J. Phys. Chem. A 110, 12017 (2006)] with the OH reagent in selected rovibrational states (v = 0-2, j = 0-5) have been carried out by means of the real wave packet (RWP) and quasiclassical trajectory (QCT) methods. State-selected total reaction probabilities have been calculated for total angular momentum J = 0 in a broad range of collision energies. Integral cross sections and state-specific rate coefficients have been obtained from the corresponding J = 0 RWP reaction probabilities for initially selected rovibrational states by means of a capture model. The calculated RWP and QCT state-selected rate coefficients are practically temperature independent. Both RWP and QCT reaction probabilities, integral cross sections, and rate coefficients are almost independent of the initial rotational excitation. The RWP results are found to be in an overall good agreement with the corresponding QCT results. The present results have been compared with earlier wave packet calculations carried out on the same potential energy surface.

Accurate Time-Dependent Wave Packet Study of the Li + H2+ Reaction and Its Isotopic Variants

The dynamics and kinetics of the Li + H₂⁺ reaction and its isotopic variants (D₂⁺ and T₂⁺) have been studied by using a time-dependent wave packet (TDWP) coupled-channel (CC) method on the ab initio potential energy surface (PES) of Martinazzo et al. [J. Chem. Phys. 2003, 119, 21]. Total initial v = 0, j = 0 state-selected reaction probabilities for the Li + H₂⁺ reaction and its isotopic variants have been calculated from the threshold up to 1 eV for total angular momenta J from 0 to 90. Integral cross sections have been evaluated from the reaction probabilities at collision energies from threshold (≈0.2 eV) up to 1.0 eV collision. The calculated rate constants as a function of temperature show an Arrhenius type behavior in the 200 ≤ T ≤ 1000 K temperature interval. It has been found to be a considerable large intermolecular kinetic isotope effect. The TDWP-CC results are in overall good agreement with those obtained applying the TDWP Centrifugal-Sudden (CS) approximation, showing that the CS approximation is rather accurate for the title reaction.

Quantum Mechanical Wave Packet and Quasiclassical Trajectory Calculations for the Li + H2+ Reaction†

The dynamics and kinetics of the Li + H2(+) reaction have been studied by means of quantum mechanical (QM) real wave packet, wave packet with flux operator, and quasiclassical trajectory (QCT) calculations on the ab initio potential energy surface of Martinazzo et al. [J. Chem. Phys., 2003, 119, 21]. Total initial state-selected reaction probabilities for the title reaction have been calculated for total angular momentum J = 0 at collision energies from threshold up to 1 eV. Wave packet reaction probabilities at selected values of the total angular momentum up to J = 60 are obtained using the centrifugal sudden approximation (CSA). Integral cross sections and rate constants have been calculated from the wave packet reactions probabilities by means of a refined J-shifting method and the separable rotation approximation in combination with the CSA for J > 0. The calculated rate constants as function of temperature show an Arrhenius type behavior. The QM results are found to be in overall good agreement with the corresponding QCT data.

Accurate Time-dependent Wave Packet Study of the H++LiH Reaction at Early Universe Conditions

The dynamics and kinetics of the H{sup +} + LiH reaction have been studied using a quantum reactive time-dependent wave packet (TDWP) coupled-channel quantum mechanical method on an ab initio potential energy surface at conditions of the early universe. The total reaction probabilities for the H{sup +} + LiH(v = 0, j = 0) {yields} H{sup +} {sub 2} + Li process have been calculated from 5 Multiplication-Sign 10{sup -3} eV up to 1 eV for total angular momenta J from 0 to 110. Using a Langevin model, integral cross sections have been calculated in that range of collision energies and extrapolated for energies below 5 Multiplication-Sign 10{sup -3} eV. The calculated rate constants are found to be nearly independent of temperature in the 10-1000 K interval with a value of Almost-Equal-To 10{sup -9} cm{sup 3} s{sup -1}, which is in good agreement with estimates used in evolutionary models of the early universe lithium chemistry.

State-to-state Quantum Wave Packet Dynamics of the LiH + H Reaction on Two Ab Initio Potential Energy Surfaces

The dynamics and kinetics of the LiH + H reaction have been studied by using an accurate quantum reactive time-dependent wave packet method on the ab initio ground electronic state potential energy surfaces (PES) developed earlier. Reaction probabilities for the two possible reaction channels, the LiH?+?H? H2?+?Li depletion process and the LiH + H?H + LiH hydrogen exchange reaction, have been calculated from 1?meV up to 1.0?eV collision energies for total angular momenta J from 0 to 80. State-to-state and total integral cross sections for the LiH-depletion and H-exchange channels of the reaction have been calculated over this collision energy range. It is found that the LiH-depletion channel is dominant in the whole range of collision energies for both PESs. Accurate total rate coefficients have been calculated on both surfaces from 100?K to 2000?K and are significantly larger than previous empirical estimates and previous J-shifting results. In addition, the present accurate calculations present noticeable differences with previous calculations using the centrifugal sudden approximation.

State-to-state chemistry and rotational excitation of CH+ in photon-dominated regions

We present a detailed theoretical study of the rotational excitation of CH+ due to reactive and nonreactive collisions involving C+(2P), H2, CH+, H and free electrons. Specifically, the formation of CH+ proceeds through the reaction between C+(2P) and H2(νH2 = 1, 2), while the collisional (de)excitation and destruction of CH+ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000 K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH+ with H atoms at kinetic temperatures below 50 K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H2 has a substantial effect on the excitation of CH+ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH+ toward the Orion Bar and the planetary nebula NGC 7027. Our results further suggest that the population of νH2 = 2 might be significant in the photon-dominated region of NGC 7027.

Quantum mechanical three-dimensional wavepacket study of the O(1D)+ClH→ClO+H reaction

Abstract The state-to-state reaction probabilities at zero total angular momentum have been calculated for the O( 1 D)+ClH→ClO+H reaction. Reaction probabilities from the ground state of ClH and going to all possible open product channels are computed over a wide range of energies. The energy dependence of the reaction probabilities shows that the reaction is strongly dominated by resonances. The product rotational and vibrational distributions are also extracted.

Accurate Time-Dependent Wave Packet Calculations for the O+ + H2 → OH+ + H Ion–Molecule Reaction

Accurate quantum reactive scattering time-dependent wave packet close-coupling calculations have been carried out to determine total reaction probabilities and integral cross sections for the O(+) + H2 → OH(+) + H reaction in a range of collision energies from 10(-3) eV up to 1.0 eV for the H2 rovibrational states (v = 0; j = 0, 1, 2) and (v = 1; j = 0) using the potential energy surface (PES) by Martı́nez et al. As expected for a barrierless reaction, the reaction cross section decays rapidly with collision energy, Ec, following a behavior that nearly corresponds to that predicted by the Langevin model. Rotational excitation of H2 into j = 1, 2 has a very moderate effect on reactivity, similarly to what happens with vibrational excitation below Ec ≈ 0.3 eV. However, at higher collision energies the cross section increases notably when H2 is promoted to v = 1. This effect is explained by resorting to the effective potentials in the entrance channel. The integral cross sections have been used to calculate rate constants in the temperature range 200-1000 K. A good overall agreement has been found with the available experimental data on integral cross sections and rate constants. In addition, time-independent quantum mechanical and quasi-classical trajectory (QCT) calculations have been performed on the same PES aimed to compare the various methodologies and to discern the detailed mechanism of the title reaction. In particular, the analysis of individual trajectories has made it possible to explain, in terms of the coupling between reagent relative velocity and the topography of the PES, the presence of a series of alternating maxima and minima in the collision energy dependence of the QCT reaction probabilities for the reactions with H2(v=0,1,j=0), which are absent in the quantum mechanical calculations.