Tomas Zimmermann

Path integral evaluation of equilibrium isotope effects

A general and rigorous methodology to compute the quantum equilibrium isotope effect is described. Unlike standard approaches, ours does not assume separability of rotational and vibrational motions and does not make the harmonic approximation for vibrations or rigid rotor approximation for the rotations. In particular, zero point energy and anharmonicity effects are described correctly quantum mechanically. The approach is based on the thermodynamic integration with respect to the mass of isotopes and on the Feynman path integral representation of the partition function. An efficient estimator for the derivative of free energy is used whose statistical error is independent of the number of imaginary time slices in the path integral, speeding up calculations by a factor of approximately 60 at 500 K and more at room temperature. We describe the implementation of the methodology in the molecular dynamics package AMBER 10. The method is tested on three [1,5] sigmatropic hydrogen shift reactions. Because of the computational expense, we use ab initio potentials to evaluate the equilibrium isotope effects within the harmonic approximation and then the path integral method together with semiempirical potentials to evaluate the anharmonicity corrections. Our calculations show that the anharmonicity effects amount up to 30% of the symmetry reduced reaction free energy. The numerical results are compared with recent experiments of Doering et al., [J. Am. Chem. Soc. 128, 9080 (2006); J. Am. Chem. Soc.129, 2488 (2007)] confirming the accuracy of the most recent measurement on 2,4,6,7,9-pentamethyl-5-(5,5-(2)H(2))methylene-11,11a-dihydro-12H-naphthacene as well as concerns about compromised accuracy, due to side reactions, of another measurement on 2-methyl-10-(10,10-(2)H(2))methylenebicyclo[4.4.0]dec-1-ene.

Efficient on-the-fly ab initio semiclassical method for computing time-resolved nonadiabatic electronic spectra with surface hopping or Ehrenfest dynamics

We derive a somewhat crude, yet very efficient semiclassical approximation for computing nonadiabatic spectra. The resulting method, which is a generalization of the multiple-surface dephasing representation, includes quantum effects through interference of mixed quantum-classical trajectories and through quantum treatment of the collective electronic degree of freedom. The method requires very little computational effort beyond the fewest-switches surface hopping or Ehrenfest locally mean-field dynamics and is very easy to implement. The proposed approximation is tested by computing the absorption and time-resolved stimulated emission spectra of pyrazine using the four-dimensional three-surface model which allows for comparison with the numerically exact quantum spectra. As expected, the multiple-surface dephasing representation is not suitable for high-resolution linear spectra, yet it seems to capture all the important features of pump-probe spectra. Finally, the method is combined with on-the-fly ab initio evaluation of the electronic structure (i.e., energies, forces, electric-dipole, and nonadiabatic couplings) in order to compute fully dimensional nonadiabatic spectra of pyrazine without approximations inherent to analytical, including vibronic-coupling models. The Appendix provides derivations of perturbative expressions for linear and pump-probe spectra of arbitrary mixed states and for arbitrary laser pulse shapes.

Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene

Recent experiments have confirmed the importance of nuclear quantum effects even in large biomolecules at physiological temperature. Here we describe how the path integral formalism can be used to describe rigorously the nuclear quantum effects on equilibrium and kinetic properties of molecules. Specifically, we explain how path integrals can be employed to evaluate the equilibrium (EIE) and kinetic (KIE) isotope effects, and the temperature dependence of the rate constant. The methodology is applied to the [1,5] sigmatropic hydrogen shift in pentadiene. Both the KIE and the temperature dependence of the rate constant confirm the importance of tunneling and other nuclear quantum effects as well as of the anharmonicity of the potential energy surface. Moreover, previous results on the KIE were improved by using a combination of a high level electronic structure calculation within the harmonic approximation with a path integral anharmonicity correction using a lower level method.

Communications: Evaluation of the nondiabaticity of quantum molecular dynamics with the dephasing representation of quantum fidelity.

We propose an approximate method for evaluating the importance of non-Born-Oppenheimer effects on the quantum dynamics of nuclei. The method uses a generalization of the dephasing representation (DR) of quantum fidelity to several diabatic potential energy surfaces and its computational cost is the cost of dynamics of a classical phase space distribution. It can be implemented easily into any molecular dynamics program and also can utilize on-the-fly ab initio electronic structure information. We test the methodology on three model problems introduced by Tully and on the photodissociation of NaI. The results show that for dynamics close to the diabatic limit, the decay of fidelity due to nondiabatic effects is described accurately by the DR. In this regime, unlike the mixed quantum-classical methods such as surface hopping or Ehrenfest dynamics, the DR can capture more subtle quantum effects than the population transfer between potential energy surfaces. Hence we propose using the DR to estimate the dynamical importance of diabatic, spin-orbit, or other couplings between potential energy surfaces. The acquired information can help reduce the complexity of a studied system without affecting the accuracy of the quantum simulation.

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schrödinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI.

Nonadiabatic effects in electronic and nuclear dynamics

Due to their very nature, ultrafast phenomena are often accompanied by the occurrence of nonadiabatic effects. From a theoretical perspective, the treatment of nonadiabatic processes makes it necessary to go beyond the (quasi) static picture provided by the time-independent Schrödinger equation within the Born-Oppenheimer approximation and to find ways to tackle instead the full time-dependent electronic and nuclear quantum problem. In this review, we give an overview of different nonadiabatic processes that manifest themselves in electronic and nuclear dynamics ranging from the nonadiabatic phenomena taking place during tunnel ionization of atoms in strong laser fields to the radiationless relaxation through conical intersections and the nonadiabatic coupling of vibrational modes and discuss the computational approaches that have been developed to describe such phenomena. These methods range from the full solution of the combined nuclear-electronic quantum problem to a hierarchy of semiclassical approaches and even purely classical frameworks. The power of these simulation tools is illustrated by representative applications and the direct confrontation with experimental measurements performed in the National Centre of Competence for Molecular Ultrafast Science and Technology.

Efficient sampling avoids the exponential wall in classical simulations of fidelity

We analyze the efficiency of available algorithms for the simulation of classical fidelity and show that their computational costs increase exponentially with the number of degrees of freedom. Then we present an algorithm for which the number of trajectories needed for convergence is independent of the systems dimensionality and show that, within a continuous family of algorithms, this algorithm is the only one with this property. Simultaneously we propose a general analytical approach to estimate efficiency of trajectory-based methods and suggest how to use it to accelerate calculations of other classical correlation functions. Converged numerical results are provided for systems with phase space volume 2^{1700} times larger than the volume of the initial state.