Dynamics of (H−,H2) collisions: A time-dependent quantum mechanical investigation on a newab initiopotential energy surface

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Quantum state-to-state study for (H−(D−),HD) collisions on two potential energy surfaces

Revealing the reaction mechanisms of the H⁻/D⁻ + HD reaction – an exact quantum dynamics study on two potential energy surfaces.

Low frequency vibrational anharmonicity and nuclear spin effects of Cl-(H2) and Cl-(D2)

Low frequency combination bands of ³⁵Cl^-(H₂) and ³⁵Cl^-(D₂) have been measured in the region between 600 and 1100 cm^-1 by infrared predissociation spectroscopy in a cryogenic 22-pole ion trap using a free electron laser at the FELIX Laboratory as a tunable light source. The ³⁵Cl^-(H₂) (³⁵Cl^-(D₂)) spectrum contains three bands at 773 cm^-1 (620 cm^-1), 889 cm^-1 (692 cm^-1), and 978 cm^-1 (750 cm^-1) with decreasing intensity toward higher photon energies. Comparison of the experimentally determined transition frequencies with anharmonic vibrational self-consistent field and vibrational configuration interaction calculations suggests the assignment of the combination bands v₁ + v₂, 2v₁ + v₂, and 3v₁ + v₂ for ³⁵Cl^-(H₂) and 2v₁ + v₂, 3v₁ + v₂, and 4v₁ + v₂ for ³⁵Cl^-(D₂), where v₁ is the ³⁵Cl^-⋯H₂ stretching fundamental and v₂ is the Cl^-(H₂) bend. The observed asymmetric temperature dependent line shape of the v₁ + v₂ transition can be modeled by a series of ∑⁺-∏ ro-vibrational transitions, when substantially decreasing the rotational constant in the vibrationally excited state by 35%. The spectrum of ³⁵Cl^-(D₂) shows a splitting of 7 cm^-1 for the strongest band which can be attributed to the tunneling of the ortho/para states of D₂.

Cation-specific conformational behavior of polyelectrolyte brushes: from aqueous to nonaqueous solvent

We have investigated changes in the cation-specific conformational behavior of poly(sodium styrenesulfonate) (PSS) brushes as the solvent changes from water to methanol using a quartz crystal microbalance with dissipation (QCM-D). A solvation to desolvation transition of the grafted chains accompanied by swelling to the collapse transition of the brushes is observed for Na(+). In the case of Cs(+), the brushes undergo solvation to desolvation to resolvation accompanied by swelling to collapse to reswelling transitions. The resolvation and reswelling transitions for Cs(+) are induced by the charge inversion of the brushes via van der Waals interactions between Cs(+) and the brushes. All of the transitions for monovalent cations become less obvious as the methanol content increases. For divalent Ca(2+) and trivalent La(3+), a solvation to desolvation to resolvation transition of the grafted chains accompanied by a swelling to collapse to reswelling transition of the brushes can be observed. The resolvation and reswelling of the brushes for the multivalent cations are induced by the charge inversion of the brushes via charge-image charge interactions. The extent of the transitions for the PSS brushes in the presence of multivalent cations is only slightly influenced by the methanol content.

Potential energy and dipole moment surfaces of H(3) (-) molecule

A new potential energy surface for the electronic ground state of the simplest triatomic anion H(3) (-) is determined for a large number of geometries. Its accuracy is improved at short and large distances compared to previous studies. The permanent dipole moment surface of the state is also computed for the first time. Nine vibrational levels of H(3) (-) and 14 levels of D(3) (-) are obtained, bound by at most approximately 70 and approximately 126 cm(-1), respectively. These results should guide the spectroscopic search of the H(3) (-) ion in cold gases (below 100K) of molecular hydrogen in the presence of H(-) ions.

Effect of Coriolis coupling in chemical reaction dynamics

It is essential to evaluate the role of Coriolis coupling effect in molecular reaction dynamics. Here we consider Coriolis coupling effect in quantum reactive scattering calculations in the context of both adiabaticity and nonadiabaticity, with particular emphasis on examining the role of Coriolis coupling effect in reaction dynamics of triatomic molecular systems. We present the results of our own calculations by the time-dependent quantum wave packet approach for H + D2 and F(2P3/2,2P1/2) + H2 as well as for the ion-molecule collisions of He + H2 +, D(-) + H2, H(-) + D2, and D+ + H2, after reviewing in detail other related research efforts on this issue.

The rates of S(N)2 reactions and their relation to molecular and solvent properties

The energy barriers of symmetrical methyl exchanges in the gas phase have been calculated with the reaction path of the intersecting/interacting-state model (ISM). Reactive bond lengths increase down a column of the Periodic Table and compensate for the decrease in the force constants, which explains the near constancy of the intrinsic barriers in the following series of nucleophiles: F(-) approximately Cl(-) approximately Br(-) approximately I(-). This compensation is absent along the rows of the Periodic Table and the trend in the reactivity is dominated by the increase in the electrophilicity index of the nucleophile in the series C<N<O<F. Solvent effects have been quantitatively incorporated into the ISM model through a correlation between electrophilicity and the solvent acceptor number. This correlation is transferable between nucleophiles and solvents and allows the methyl transfer rate constants in solution to be calculated with remarkable simplicity and accuracy. The relationship between the S(N)2 and electron-transfer mechanisms is clarified and it is shown that smaller solvent static effects should be expected for electron transfer in the absence of a thermodynamic driving force.

Influence of Reagent Rotation on (H-, D2) and (D-, H2) Collisions: A Quantum Mechanical Study

Time-independent quantum mechanical (TIQM) approach (helicity basis truncated at k = 2) has been used for computing differential and integral cross sections for the exchange reaction H- + D2 (v = 0, j = 0-4) --> HD + D- and D- + H2 (v = 0, j = 0-3) --> HD + H- in three dimensions on an accurate ab initio potential energy surface. It is shown that the j-weighted differential reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder at four different relative translational energies (Etrans = 0.55, 0.93, 1.16 and 1.48 eV) for (H-, D2) and at one relative translational energy (Etrans = 0.6 eV) by Haufler et al. for both (H-, D2) and (D-, H2) collisions. The j-weighted integral reaction cross section values are in good agreement with the crossed beam measurements by Zimmer and Linder in the Etrans range 0.5-1.5 eV and close to the guided ion beam results by Haufler et al. for (H-, D2) in the range 0.8-1.2 eV. Time-dependent quantum mechanical (TDQM) results obtained using centrifugal sudden approximation are reported in the form of integral reaction cross section values as a function of Etrans in the range 0.3-3.0 eV for both reactions in three dimensions on the same potential energy surface. The TDQM reaction cross section values decline more sharply than the TIQM results with increase in the initial rotational quantum number (j) for the D2 molecules in their ground vibrational state (v = 0) for (H-, D2) collisions. The computed j-weighted reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder for (H-, D2) collisions and guided ion beam results by Haufler et al. for both (H-, D2) and (D-, H2) collisions for energies below the threshold for electron detachment channel.

Three Dimensional Quantum Dynamics of (H-, H2) and Its Isotopic Variants

We present the results of a time-dependent quantum mechanical investigation using centrifugal sudden approximation in the form of reaction probability as a function of collision energy (E(trans)) in the range 0.3-3.0 eV for a range of total angular momentum (J) values and the excitation function sigma(E(trans)) for the exchange reaction H(-) + H(2) (v = 0, j = 0) --> H(2) + H(-) and its isotopic variants in three dimensions on an accurate ab initio potential energy surface published recently (J. Chem. Phys. 2004, 121, 9343). The excitation function results are shown to be in excellent agreement with those obtained from crossed beam measurements by Zimmer and Linder for H(-) + D(2) collisions for energies below the threshold for electron detachment channel and somewhat larger than the most recent results of Haufler et al. for (H(-), D(2)) and (D(-), H(2)) collisions.