Tunneling Splittings in Water Clusters from Path Integral Molecular Dynamics.

We present calculations of tunneling splittings in selected small water clusters, based on a recently developed path integral molecular dynamics (PIMD) method. The ground-rotational-state tunneling motions associated with the largest splittings in the water dimer, trimer, and hexamer are considered, and we show that the PIMD predictions are in very good agreement with benchmark quantum and experimental results. As the tunneling spectra are highly sensitive to both the details of the quantum dynamics and the potential energy surface, our calculations are a validation of the MB-Pol surface as well as the accuracy of PIMD. The favorable scaling of PIMD with system size paves the way for calculations of tunneling splittings in large, nonrigid molecular systems with motions that cannot be treated accurately by other methods, such as the semiclassical instanton.

Tunneling splittings from path-integral molecular dynamics using a Langevin thermostat.

We report an improved method for the calculation of tunneling splittings between degenerate configurations in molecules and clusters using path-integral molecular dynamics (PIMD). Starting from an expression involving a ratio of thermodynamic density matrices at the bottom of the symmetric wells, we use thermodynamic integration with molecular dynamics simulations and a Langevin thermostat to compute the splittings stochastically. The thermodynamic integration is performed by sampling along the semiclassical instanton path, which provides an efficient reaction coordinate as well as being physically well-motivated. This approach allows us to carry out PIMD calculations of the multi-well tunneling splitting pattern in the water dimer and to refine previous PIMD calculations for one-dimensional models and malonaldehyde. The large (acceptor) splitting in the water dimer agrees to within 20% of benchmark variational results, and the smaller splittings agree to within 10%.

Nearside-farside and local angular momentum analyses of time-independent scattering amplitudes for the H + D2 (vi = 0, ji = 0) --> HD (vf = 3, jf = 0) + D reaction.

The scattering dynamics of the state-to-state reaction H + D2 (v(i) = 0, j(i) = 0, m(i) = 0) --> HD (v(f) = 3, j(f) = 0, m(f) = 0) + D is investigated, where vi, ji, mi and vf, jf, mf are initial and final vibrational, rotational, and helicity quantum numbers, respectively. We use accurate quantum scattering matrix elements for total energies in the range 1.52-2.50 eV (calculated stepwise in 0.01 eV increments). The theoretical tools used are a nearside-farside (NF) analysis of the partial wave series (PWS) for the scattering amplitude, together with NF local angular momentum (LAM) theory. We find that the backward scattering, which is the energy-domain analog of the time-direct reaction mechanism, is N dominated, whereas the forward scattering (time-delayed analog) is a result of NF interference between the more slowly varying N and F subamplitudes. The LAM analysis reveals the existence of a "trench-ridge" structure. We also resum the PWS up to three times prior to making the NF decomposition. We show that such resummations usually provide an improved physical interpretation of the NF differential cross sections (DCSs) and NF LAMs. We analyze two resummed scattering amplitudes in more detail, where particular values of the resummation parameters give rise to unexpected unphysical behavior in the N and F DCSs over a small angular range. We analyze the cause of this unphysical behavior and describe viable workarounds to the problem. The energy-domain calculations in this paper complement the time-domain results reported earlier by Monks, P. D. D.; Connor, J. N. L.; Althorpe, S. C. J. Phys. Chem. A 2006, 110, 741.

Theory of time-dependent reactive scattering: cumulative time-evolving differential cross sections and nearside-farside analyses of time-dependent scattering amplitudes for the H + D2 --> HD + D reaction.

Nearside-farside (NF) theory, originally developed in the energy domain for the time-independent description of molecular collisions and chemical reactions, is applied to the plane wave packet (PWP) formulation of time-dependent scattering. The NF theory decomposes the partial wave series representation for the time-dependent PWP scattering amplitude into two time-dependent subamplitudes: one N, the other F. In addition, NF local angular momentum (LAM) theory is applied to the PWP scattering amplitude. The novel concept of a cumulative time-evolving differential cross section is introduced, in which the upper infinite time limit of a half-Fourier transform is replaced by a finite time. In a similar way, a cumulative energy-evolving angular distribution is defined. Application is made to the state-to-state reaction, H + D2(v(i) = 0, j(i) = 0) --> HD(v(f) = 3, j(f) = 0) + D, where v(i), j(i) and v(f), j(f) are vibrational and rotational quantum numbers for the initial and final states, respectively. This reaction exhibits time-direct and time-delayed (by about 25 fs) collision mechanisms. It is shown that the direct-time mechanism is N dominant scattering, whereas the time-delayed mechanism exhibits characteristics of NF interference. The NF and LAM theories provide valuable insights into the time-dependent properties of a reaction, as do snapshots from a movie of the cumulative time-evolving differential cross section.