The influence of nonadiabatic rotational transitions on the line shapes of the rotational Raman spectrum of H2 in liquid argon

Quantum transition probabilities due to overlapping electromagnetic pulses: Persistent differences between Dirac's form and nonadiabatic perturbation theory

The probability of transition to an excited state of a quantum system in a time-dependent electromagnetic field determines the energy uptake from the field. The standard expression for the transition probability has been given by Dirac. Landau and Lifshitz suggested, instead, that the adiabatic effects of a perturbation should be excluded from the transition probability, leaving an expression in terms of the nonadiabatic response. In our previous work, we have found that these two approaches yield different results while a perturbing field is acting on the system. Here, we prove, for the first time, that differences between the two approaches may persist after the perturbing fields have been completely turned off. We have designed a pair of overlapping pulses in order to establish the possibility of lasting differences, in a case with dephasing. Our work goes beyond the analysis presented by Landau and Lifshitz, since they considered only linear response and required that a constant perturbation must remain as t → ∞. First, a "plateau" pulse populates an excited rotational state and produces coherences between the ground and excited states. Then, an infrared pulse acts while the electric field of the first pulse is constant, but after dephasing has occurred. The nonadiabatic perturbation theory permits dephasing, but dephasing of the perturbed part of the wave function cannot occur within Dirac's method. When the frequencies in both pulses are on resonance, the lasting differences in the calculated transition probabilities may exceed 35%. The predicted differences are larger for off-resonant perturbations.

Application of mean-field and surface-hopping approaches for interrogation of the Xe3+ molecular ion photoexcitation dynamics

The dissociation dynamics of the excited Xe(3) (+) molecular ion through the Pi(12)(u) and Pi(12)(g) conical intersection was interrogated by computational simulation in which no adjustable parameters were used. The electronic ground and excited state potential energy surfaces were generated by the diatomics-in-molecules method, and the Ehrenfest mean-field and Tully surface-hopping approaches treated the nonadiabatic interactions. Reproduction of the experimental spectrum of the symmetric photofragmentation as a function of excitation energy was obtained within the region of interest (2.5-3.75 eV), with the exception of a 0.25 eV width on the red side of the spectral apex. Good agreement was obtained with the experimental dissociated photofragment kinetic energy spectra. It was determined that the greatest contribution to the nonadiabatic coupling between the two states originated from the bending vibrational mode of the molecule in the Sigma(12)(u), ground electronic state before excitation.

Progress in the theory of mixed quantum-classical dynamics

Quantum-classical Liouville dynamics can be used to study the properties of open quantum systems that are coupled to bath or environmental degrees of freedom whose dynamics can be approximated by classical equations of motion. In contrast to many open quantum system approaches, quantum-classical dynamics provides a detailed description of the bath molecules. Such a description is especially appropriate for the study of quantum rate processes, such as proton and electron transport, where the detailed dynamics of the bath has a strong influence on the quantum rate. The quantum-classical Liouville equation can also serve as a starting point for the derivation of reduced descriptions where all or some of the bath degrees of freedom are projected out. Quantum-classical Liouville dynamics can be simulated in terms of an ensemble of surface-hopping trajectories whose character differs from that in other surface-hopping schemes. The results of studies of proton transfer in condensed phase and reactive dynamics in a dissipative environment are presented to illustrate applications of the formalism.

Semiclassical calculation of nonadiabatic thermal rate constants: Application to condensed phase reactions

The nonadiabatic transition state theory proposed recently by Zhao et al. [J. Chem. Phys. 121, 8854 (2004)] is extended to calculate rate constants of complex systems by using the Monte Carlo and umbrella sampling methods. Surface hopping molecular dynamics technique is incorporated to take into account the dynamic recrossing effect. A nontrivial benchmark model of the nonadiabatic reaction in the condensed phase is used for the numerical test. It is found that our semiclassical results agree well with those produced by the rigorous quantum mechanical method. Comparing with available analytical approaches, we find that the simple statistical theory proposed by Straub and Berne [J. Chem. Phys. 87, 6111 (1987)] is applicable for a wide friction region although their formula is obtained using Landau-Zener [Phys. Z. Sowjetunion 2, 46 (1932); Proc. R. Soc. London, Ser. A 137, 696 (1932)] nonadiabatic transition probability along a one-dimensional diffusive coordinate. We also investigate how the nuclear tunneling events affect the dependence of the rate constant on the friction.

Electronic relaxation dynamics of Ni2+-ion aqueous solution: Molecular-dynamics simulation

Electronic relaxation dynamics of Ni2+-ion aqueous solution is investigated using molecular-dynamics (MD) simulations with the model-effective Hamiltonian developed previously. The nonadiabatic transition rates from the first three excited states to the ground state are evaluated by the golden rule formula with the adiabatic MD simulations. The MD simulations with the fewest-switch surface-hopping method are also carried out to obtain a more detailed description of the electronic relaxation dynamics among the excited states. We found out that the transitions among the three excited states are very fast, in the order of 10 fs, while the transition between the excited and ground states is slow, about 800 ps. These findings are consistent with the time scales of energy dissipation detected by the transient lens experiment. In both simulations, we explore the effects of the quantum decoherence, where the decoherence functions are derived by the energy-gap dynamics with the displaced harmonic-oscillator model.

Quantum-classical limit of quantum correlation functions

A quantum-classical limit of the canonical equilibrium time correlation function for a quantum system is derived. The quantum-classical limit for the dynamics is obtained for quantum systems comprising a subsystem of light particles in a bath of heavy quantum particles. In this limit the time evolution of operators is determined by a quantum-classical Liouville operator, but the full equilibrium canonical statistical description of the initial condition is retained. The quantum-classical correlation function expressions derived here provide a way to simulate the transport properties of quantum systems using quantum-classical surface-hopping dynamics combined with sampling schemes for the quantum equilibrium structure of both the subsystem of interest and its environment