Unified semiclassical theory for the two-state system: an analytical solution for general nonadiabatic tunneling

Unified semiclasical solution for general nonadiabatic tunneling between two adiabatic potential energy surfaces is established by employing unified semiclassical solution for pure nonadiabatic transition [C. Zhu, J. Chem. Phys. 105, 4159 (1996)] with the certain symmetry transformation. This symmetry comes from a detailed analysis of the reduced scattering matrix for Landau-Zener type of crossing as a special case of nonadiabatic transition and nonadiabatic tunneling. Traditional classification of crossing and noncrossing types of nonadiabatic transition can be quantitatively defined by the rotation angle of adiabatic-to-diabatic transformation, and this rotational angle enters the analytical solution for general nonadiabatic tunneling. The certain two-state exponential potential models are employed for numerical tests, and the calculations from the present general nonadiabatic tunneling formula are demonstrated in very good agreement with the results from exact quantum mechanical calculations. The present general nonadiabatic tunneling formula can be incorporated with various mixed quantum-classical methods for modeling electronically nonadiabatic processes in photochemistry.     

Quantum effect of intramolecular high-frequency vibrational modes on diffusion-controlled electron transfer rate: from the weak to the strong electronic coupling regions

The Sumi-Marcus theory is extended by introducing two approaches to investigate electron transfer reactions from weak-to-strong electronic coupling regime. One of these approaches is the quantum R-matrix theory, useful for dealing with the intramolecular vibrational motions in the whole electronic coupling domain. The other is the split operator approach that is employed to solve the reaction-diffusion equation. The approaches are then applied to electron transfer in the Marcus inverted regime to investigate the nuclear tunneling effect on the long time rate and the survival probabilities. The numerical results illustrate that the adiabatic suppression obtained from the R-matrix approach is much smaller than that from the Landau-Zener theory whereas it cannot be predicted by the perturbation theory. The jointed effects of the electronic coupling and solvent relaxation time on the rates are also explored.     

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.     

Reaction-rate approach to some low temperature processes in alkali halides

The general reaction-rate theory is applied to some processes in alkali halide crystals that occur by lattice tunneling at low temperature. The theoretical equations used refer to the strong coupling case comprising both adiabatic and nonadiabatic transitions. The thermal ionization of excited F centers is considered from this point of view. The theory reproduces rather well the experimental temperature dependence of the ionization efficiency. Similarly the reorientation of OH^? impurity dipoles is treated and a good agreement is again found between theoretical and experimental relaxation rates in several host crystals.     

Dynamics of nonadiabatic chemical reactions

New methods are proposed to treat nonadiabatic chemical dynamics in realistic large molecular systems by using the Zhu-Nakamura (ZN) theory of curve-crossing problems. They include the incorporation of the ZN formulas into the Herman-Kluk type semiclassical wave packet propagation method and the trajectory surface hopping (TSH) method, formulation of the nonadiabatic transition state theory, and its application to the electron transfer problem. Because the nonadiabatic coupling is a vector in multidimensional space, the one-dimensional ZN theory works all right. Even the classically forbidden transitions can be correctly treated by the ZN formulas. In the case of electron transfer, a new formula that can improve the celebrated Marcus theory in the case of normal regime is obtained so that it can work nicely in the intermediate and strong electronic coupling regimes. All these formulations mentioned above are demonstrated to work well in comparison with the exact quantum mechanical numerical solutions and are expected to be applicable to large systems that cannot be treated quantum mechanically numerically exactly. To take into account another quantum mechanical effect, namely, the tunneling effect, an efficient method to detect caustics from which tunneling trajectories emanate is proposed. All the works reported here are the results of recent activities carried out in the author's research group. Finally, the whole set of ZN formulas is presented in Appendix.     

Semiclassical treatment of thermally activated electron transfer in the intermediate to strong electronic coupling regime under the fast dielectric relaxation

The generalized nonadiabatic transition-state theory (NA-TST) (Zhao, Y.; et al. J. Chem. Phys. 2004, 121, 8854) is used to study electron transfer with use of the Zhu-Nakamura (ZN) formulas of nonadiabatic transition in the case of fast dielectric relaxation. The rate constant is expressed as a product of the well-known Marcus formula and a coefficient which represents the correction due to the strong electronic coupling. In the case of general multidimensional systems, the Monte Carlo approach is utilized to evaluate the rate by taking into account the multidimensionality of the crossing seam surface. Numerical demonstration is made by using a model system of a collection of harmonic oscillators in the Marcus normal region. The results are naturally coincident with the perturbation theory in the weak electronic coupling limit; while in the intermediate to strong electronic coupling regime where the perturbation theory breaks down the present results are in good agreement with those from the quantum mechanical flux-flux correlation function within the model of effective one-dimensional mode.     

Semiclassical theory of electronically nonadiabatic chemical dynamics: incorporation of the Zhu-Nakamura theory into the frozen Gaussian propagation method

The title theory is developed by combining the Herman-Kluk semiclassical theory for adiabatic propagation on single potential-energy surface and the semiclassical Zhu-Nakamura theory for nonadiabatic transition. The formulation with use of natural mathematical principles leads to a quite simple expression for the propagator based on classical trajectories and simple formulas are derived for overall adiabatic and nonadiabatic processes. The theory is applied to electronically nonadiabatic photodissociation processes: a one-dimensional problem of H2+ in a cw (continuous wave) laser field and a two-dimensional model problem of H2O in a cw laser field. The theory is found to work well for the propagation duration of several molecular vibrational periods and wide energy range. Although the formulation is made for the case of laser induced nonadiabatic processes, it is straightforwardly applicable to ordinary electronically nonadiabatic chemical dynamics.     

Ultrafast electron-transfer and solvent adiabaticity effects in viologen charge-transfer complexes

We report ultrafast electron transfer (ET) in charge-transfer complexes that shows solvent relaxation effects consistent with adiabatic crossover models of nonadiabatic ET. The complexes of either dimethyl viologen (MV) or diheptyl viologen (HV) with 4,4'-biphenol (BP) (MVBP or HVBP complexes) have identical charge-transfer spectra and kinetics in ethylene glycol with approximately 900 fs ET decay. We assign this decay time as largely due to adiabatic control of a predicted approximately 40 fs nonadiabatic ET. The MVBP decay in methanol of 470 fs is reduced in mixtures having low (2-20%) concentrations of acetonitrile to as short as 330 fs; these effects are associated with faster relaxation time in methanol and its mixtures. In contrast, HVBP has much longer ET decay in methanol (730 fs) and mixture effects that only reduce its decay to 550 fs. We identify the heptyl substituent as creating major perturbations to solvent relaxation times in the methanol solvation shell of HVBP. These charge-transfer systems have reasonably well-defined geometry with weak electronic coupling where the electronic transitions are not dependent on intramolecular motions. We used a nonadiabatic ET model with several models for adiabatic crossover predictions to discuss the small variation of energy gap with solvent and the ET rates derived from adiabatic solvent control. A time correlation model of solvent relaxation was used to define the solvent relaxation times for this case of approximately zero-barrier ET.     

Electronic nonadiabatic effects in the adsorption of hydrogen atoms on metals

The time-dependent, mean-field Newns-Anderson model for a spin-polarized adsorbate approaching a metallic surface is solved in the wide-band limit. Equations for the time evolution of the occupation of the spin dependent adsorbate states and for the nonadiabatic and nearly adiabatic adsorbate-surface energy transfer rates are derived. Numerical solutions are obtained using characteristic parameters derived from density functional theory calculations for the H/Cu(111) system. The time evolution of the model system is shown to be strongly nonadiabatic in the vicinity of the transition point between spin-polarized and nonpolarized ground states. Away from the spin transition the nonadiabatic energy transfer is in close agreement with the nearly adiabatic limit. Near the transition, nonadiabatic effects are large and the nearly adiabatic approximation fails.     

Evaluation of canonical and microcanonical nonadiabatic reaction rate constants by using the Zhu-Nakamura formulas

We consider a problem of calculating both thermal and microcanonical rate constants for nonadiabatic chemical reactions. Instead of using the conventional transition state theory, we use a generalized seam surface and introduce a concept of a coordinate dependent effective nonadiabatic transition probability based on the Zhu-Nakamura theory which can treat the nonadiabatic tunneling properly. The present approach can be combined with Monte Carlo method so as to be applicable to chemical reactions in complicated systems. The method is demonstrated to work well in wide energy and temperature range. Numerical tests also show that it is very essential for accurate evaluation of the thermal rate constant to use the generalized seam surface and take into account the nonadiabatic tunneling effect.     

Vibrational wavefunctionals for the thermal and photochemical syn—anti isomerization of imines

Employing the partitioning scheme for the nuclear coordinates, we recently have discussed the vibration—isomerizational hamiltonian appropriate for the syn—anti rearrangement of imines. Now we present the calculations of the vibrational wavefunctions and energies which were needed for the investigation of the horizontal and vertical isomerization processes. The energy splitting of the adiabatic zero-order states shows clearly the importance of both inversion and rotation, for the rearrangement process in the ground and first triplet state of methylenimine. Using these zero-order states we then have investigated the question as to what extent nonadiabatic coupling influences the rate of the horizontal and vertical reactions. Further we have discussed two mechanisms for the nonadiabatic effects concerning the tunneling frequency, which depend on whether the nonrelevant coordinated are “symmetric” or “antisymmetric” with respect ot the symmetry operator describing the isomerizational states. The use of the vibrational adiabatic states for the calculation of the photoisomerization rates has also been suggested.     

Nonadiabatic transitions between lambda-doubling states in the capture of a diatomic molecule by an ion

The low-energy capture of a dipolar diatomic molecule in an adiabatically isolated electronic state with a good quantum number Omega (Hund's coupling case a) by an ion occurs adiabatically with respect to rotational transitions of the diatom. However, the capture dynamics may be nonadiabatic with respect to transitions between the pair of the Lambda-doubling states belonging to the same value of the intrinsic angular momentum j. In this work, nonadiabatic transition probabilities are calculated which define the Lambda-doubling j-specific capture rate coefficients. It is shown that the transition from linear to quadratic Stark effect in the ion-dipole interaction, which damps the T(-1/2) divergence of the capture rate coefficient calculated with vanishing Lambda-doubling splitting, occurs in the adiabatic regime with respect to transitions between Lambda-doubling adiabatic channel potentials. This allows one to suggest simple analytical expressions for the rate coefficients in the temperature range which covers the region between the sudden and the adiabatic limits with respect to the Lambda-doubling states.     

Sampling the proton transfer reaction coordinate in mixed quantum-classical molecular dynamics simulations

An umbrella sampling approach based on the vibrational energy gap is presented and examined for exploring the reaction coordinate for a proton transfer (PT) reaction. The technique exploits the fact that for a PT reaction the energy gap between the vibrational ground and excited states of the transferring proton reaches a minimum at the transition state. Umbrella sampling is used within mixed quantum-classical simulations to identify the transition state configurations and explore the reaction free energy curve and vibrationally nonadiabatic coupling. The method is illustrated by application to a model phenol-amine proton transfer reaction complex in a nanoconfined solvent. The results from this new umbrella sampling approach are consistent with those obtained from previous umbrella sampling calculations based on a collective solvent coordinate. This sampling approach further provides insight into the vibrationally nonadiabatic coupling for the proton transfer reaction and has potential for simulating vibrational spectra of PT reaction complexes in solution.     

Nonadiabatic transitions of fluctuation origin

The quasi-classical approximation is used to find a solution for nonadiabatic transition by tunneling in a system of two parabolic terms that interact with the surrounding medium. The coupling between the terms is arbitrary.     

Calculation of vibronic couplings for phenoxyl/phenol and benzyl/toluene self-exchange reactions: implications for proton-coupled electron transfer mechanisms

The vibronic couplings for the phenoxyl/phenol and the benzyl/toluene self-exchange reactions are calculated with a semiclassical approach, in which all electrons and the transferring hydrogen nucleus are treated quantum mechanically. In this formulation, the vibronic coupling is the Hamiltonian matrix element between the reactant and product mixed electronic-proton vibrational wavefunctions. The magnitude of the vibronic coupling and its dependence on the proton donor-acceptor distance can significantly impact the rates and kinetic isotope effects, as well as the temperature dependences, of proton-coupled electron transfer reactions. Both of these self-exchange reactions are vibronically nonadiabatic with respect to a solvent environment at room temperature, but the proton tunneling is electronically nonadiabatic for the phenoxyl/phenol reaction and electronically adiabatic for the benzyl/toluene reaction. For the phenoxyl/phenol system, the electrons are unable to rearrange fast enough to follow the proton motion on the electronically adiabatic ground state, and the excited electronic state is involved in the reaction. For the benzyl/toluene system, the electrons can respond virtually instantaneously to the proton motion, and the proton moves on the electronically adiabatic ground state. For both systems, the vibronic coupling decreases exponentially with the proton donor-acceptor distance for the range of distances studied. When the transferring hydrogen is replaced with deuterium, the magnitude of the vibronic coupling decreases and the exponential decay with distance becomes faster. Previous studies designated the phenoxyl/phenol reaction as proton-coupled electron transfer and the benzyl/toluene reaction as hydrogen atom transfer. In addition to providing insights into the fundamental physical differences between these two types of reactions, the present analysis provides a new diagnostic for differentiating between the conventionally defined hydrogen atom transfer and proton-coupled electron transfer reactions.     

利用基于直接动力学的轨线面跳跃方法处理非绝热过程

势能面交叉引起的非绝热过程广泛存在于光化学和光物理中。对这一过程进行描述是理论化学的重要挑战之一。非绝热过程涉及原子核与电子之间的耦合运动,因此量子化学的基本假设之一"玻恩-奥本海默"近似被打破,所以对其进行描述需要发展新的动力学理论方法。在这些方法中,Tully发展的最少轨线面跳跃方法凭借易于程序化、便于计算等优点已经发展成为处理非绝热问题的主要动力学方法之一。其中原子核以经典的方式在单一势能面上进行演化,电子以量子的方式沿着同一轨线进行演化。在整个演化过程中,非绝热跃迁通过轨线在不同势能面间的跃迁来描述,其中跳跃发生的几率与电子的演化有关。如果将该方法与从头算直接动力学相结合,可以在全原子水平上研究实际分子体系的非绝热动力学,给出其激发态寿命、非绝热动力学中分子的主要运动方式、反应通道以及分支比等重要信息。本文旨在讨论最少面跳跃直接动力学方法研究非绝热问题的一些进展,包括动力学基本理论,特别关注将最少面跳跃方法和直接动力学结合的数值实现细节,同时讨论该方法在研究实际体系当中的一些应用,并对轨线面跳跃方法下一步发展的一些方向进行合理的展望。