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.     

Coherent correlation between nonadiabatic rotational excitation and angle-dependent ionization of NO in intense laser fields

We investigate coherent correlation between nonadiabatic rotational excitation and angle-dependent ionization of NO in intense laser fields in the state-resolved manner. When neutral NO molecules are partly ionized in intense laser fields (I(0) > 35 TW/cm(2)), a hole in the rotational wave packet of the remaining neutral NO is created because of the ionization rate depending on the alignment angle of the molecular axis with respect to the laser polarization direction. Rotational state distributions of NO are experimentally observed, and then the characteristic feature that the population at higher J levels is increased by the ionization can be identified. Numerical calculation for solving time-dependent rotational Schro?dinger equations including the effect of the ionization is carried out. The numerical results suggest that NO molecules aligned perpendicular to the laser polarization direction are dominantly ionized at the peak intensity of I(0) = 42 TW/cm(2), where the multiphoton ionization is preferred rather than the tunneling ionization.     

Quantum-classical Liouville dynamics of nonadiabatic proton transfer

A proton transfer reaction in a linear hydrogen-bonded complex dissolved in a polar solvent is studied using mixed quantum-classical Liouville dynamics. In this system, the proton is treated quantum mechanically and the remainder of the degrees of freedom is treated classically. The rates and mechanisms of the reaction are investigated using both adiabatic and nonadiabatic molecular dynamics. We use a nonadiabatic dynamics algorithm which allows the system to evolve on single adiabatic surfaces and on coherently coupled pairs of adiabatic surfaces. Reactive-flux correlation function expressions are used to compute the rate coefficients and the role of the dynamics on the coherently coupled surfaces is elucidated.     

Theory of three-dimensional alignment by intense laser pulses

We introduce a theoretical framework for study of three-dimensional alignment by moderately intense laser pulses and discuss it at an elementary level. Several features of formal interest are noted and clarified. Our approach is nonperturbative, treating the laser field within classical and the material system within quantum mechanics. The theory is implemented numerically using a basis set of rotational eigenstates, transforming the time-dependent Schrodinger equation to a set of coupled differential equations where all matrix elements are analytically soluble. The approach was applied over the past few years to explore different adiabatic and nonadiabatic three-dimensional alignment approaches in conjunction with experiments, but its formal details and numerical implementation were not reported in previous studies. Although we provide simple numerical examples to illustrate the content of the equations, our main goal is to complement previous reports through an introductory discussion of the underlying theory.     

Role of rotational temperature in adiabatic molecular alignment

One-dimensional alignment of molecules in the adiabatic limit, where the pulse duration greatly exceeds the molecular rotational periods, is studied experimentally. Four different asymmetric top molecules (iodobenzene, p-diiodobenzene, 3,4-dibromothiophene, and 4,4'-dibromobiphenyl), rotationally cooled through a high pressure supersonic pulsed valve, are aligned by a 9-ns-long pulse. Their orientations are measured through Coulomb explosion, induced by a 130-fs-long pulse, and by recording the direction of the recoiling ions. The paper focuses on the crucial role of the initial rotational temperature for the degree of alignment. In particular, we show that at molecular temperatures in the 1 K range very strong alignment is obtained already at intensities of a few times 10(11) W/cm2 for all four molecules. At the highest intensities (approximately 10(12) W/cm2) the molecules can tolerate without ionizing >or=0.92 in the case of iodobenzene. This is the strongest degree of alignment ever reported for any molecule.     

Mixed quantum-classical dynamics in the adiabatic representation to simulate molecules driven by strong laser pulses

The dynamics of molecules under strong laser pulses is characterized by large Stark effects that modify and reshape the electronic potentials, known as laser-induced potentials (LIPs). If the time scale of the interaction is slow enough that the nuclear positions can adapt to these externally driven changes, the dynamics proceeds by adiabatic following, where the nuclei gain very little kinetic energy during the process. In this regime we show that the molecular dynamics can be simulated quite accurately by a semiclassical surface-hopping scheme formulated in the adiabatic representation. The nuclear motion is then influenced by the gradients of the laser-modified potentials, and nonadiabatic couplings are seen as transitions between the LIPs. As an example, we simulate the process of adiabatic passage by light induced potentials in Na(2) using the surface-hopping technique both in the diabatic representation based on molecular potentials and in the adiabatic representation based on LIPs, showing how the choice of the representation is crucial in reproducing the results obtained by exact quantum dynamical calculations.     

X-ray diffraction and molecular simulation study of the crystalline and liquid states of succinic anhydride

The crystal structure of succinic anhydride was studied at five temperatures between 100 K and the melting point by single-crystal X-ray diffraction. The temperature dependence of molecular libration tensors was determined. Intermolecular interactions, in particular through unusually close molecule-molecule contacts, are discussed, with a detailed calculation of electrostatic energies. A method for the adaptation of existing crystal force fields to molecular dynamics has been developed; the adapted force field was used to study molecular motion and rotational diffusion with increasing temperature. Equilibration of the crystalline system becomes impossible at a temperature very close to the experimental melting temperature, where a sudden transition to the liquid state occurs, and a partial kinetic picture of the melting process is obtained. After validation of the force field against experimental crystal data, the state equation of the liquid was predicted. Enthalpies of sublimation, melting, and vaporization were calculated. The dynamics of a solution of succinic anhydride in a nonpolar solvent was simulated, for a discussion of the aggregation process leading to demixing and to crystal nucleation.     

Quantum and classical approaches for rotational relaxation and nonresonant laser alignment of linear molecules: a comparison for CO2 gas in the nonadiabatic regime

A quantum approach and classical molecular dynamics simulations (CMDS) are proposed for the modeling of rotational relaxation and of the nonadiabatic alignment of gaseous linear molecules by a nonresonant laser field under dissipative conditions. They are applied to pure CO(2) and compared by looking at state-to-state collisional rates and at the value of induced by a 100 fs laser pulse linearly polarized along z[overhead arrow]. The main results are: (i) When properly requantized, the classical model leads to very satisfactory predictions of the permanent and transient alignments under non-dissipative conditions. (ii) The CMDS calculations of collisional-broadening coefficients and rotational state-to-state rates are in very good agreement with those of a quantum model based on the energy corrected sudden (ECS) approximation. (iii) Both approaches show a strong propensity of collisions, while they change the rotational energy (i.e., J), to conserve the angular momentum orientation (i.e., M/J). (iv) Under dissipative conditions, CMDS and quantum-ECS calculations lead to very consistent decays with time of the "permanent" and transient components of the laser-induced alignment. This result, expected from (i) and (ii), is obtained only if a properly J- and M-dependent ECS model is used. Indeed, rotational state-to-state rates and the decay of the "permanent" alignment demonstrate, for pure CO(2), the limits of a M-independent collisional model proposed previously. Furthermore, computations show that collisions induce a decay of the "permanent" alignment about twice slower than that of the transient revivals amplitudes, a direct consequence of (iii). (v) The analysis of the effects of reorienting and dephasing elastic collisions shows that the latter have a very small influence but that the former play a non-negligible role in the alignment dynamics. (vi) Rotation-translation collisionally induced transfers have also been studied, demonstrating that they only slightly change the alignment dissipation for the considered laser energy conditions.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

Optimal control of molecular alignment in dissipative media

We explore the controllability of nonadiabatic alignment in dissipative media, and the information content of control experiments regarding the bath properties and the bath system interactions. Our approach is based on a solution of the quantum Liouville equation within the multilevel Bloch formalism, assuming Markovian dynamics. We find that the time and energy characteristics of the laser fields that produce desired alignment characteristics at a predetermined instant respond in distinct manners to decoherence and to population relaxation, and are sensitive to both time scales. In particular, the time-evolving spectral composition of the optimal pulse mirrors the time-evolving rotational composition of the wave packet, and points to different mechanisms of rotational excitation in isolated systems, in systems subject to a decoherering bath, and in ones subject to a population relaxing bath.     

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

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

Quantum-classical dynamics of wave fields

An approach to the quantum-classical mechanics of phase space dependent operators, which has been proposed recently, is remodeled as a formalism for wave fields. Such wave fields obey a system of coupled nonlinear equations that can be written by means of a suitable non-Hamiltonian bracket. As an example, the theory is applied to the relaxation dynamics of the spin-boson model. In the adiabatic limit, a good agreement with calculations performed by the operator approach is obtained. Moreover, the theory proposed in this paper can take nonadiabatic effects into account without resorting to surface-hopping approximations. Hence, the results obtained follow qualitatively those of previous surface-hopping calculations and increase by a factor of (at least) 2, the time length over which nonadiabatic dynamics can be propagated with small statistical errors. Moreover, it is worth to note that the dynamics of quantum-classical wave fields proposed here is a straightforward non-Hamiltonian generalization of the formalism for nonlinear quantum mechanics that Weinberg introduced recently.     

On the possibility of aligning paramagnetic molecules or ions in a magnetic field

Strong magnetic fields can hybridize low rotational states of paramagnetic molecules or molecular ions whose electronic angular momentum is coupled to the molecular axis. The hybridization creates pendular states in which the molecular axis is confined to librate over a limited angular range about the field direction. In this way substantial spatial alignment associated with large Zeeman shifts can be attained for many ground-state radicals or ions and electronically excited states of diatomic or linear molecules. The magnetic hybridization is analogous to that recently demonstrated for polar molecules in electric fields. The magnetic version can only provide ensemble alignment rather than orientation, but offers complementary chemical scope by virtue of its applicability to nonpolar molecules and ions.     

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.     

Time-dependent density-functional theory beyond the adiabatic approximation: insights from a two-electron model system

Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local-density approximation (ALDA) for the dynamical exchange-correlation potential V(xc)(r,t). An exact (i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a V(xc)(r,t) with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two-dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with numerical solutions of the two-electron time-dependent Schrodinger equation. It is shown explicitly how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension of the ALDA fails for finite systems but becomes correct in the thermodynamic limit.     

Alignment of molecules in pulsed resonant laser fields

We investigate by numerical simulations the dynamics of alignment of linear molecules in resonant pulsed laser fields and its dependence on pulse length, field strength, and molecular parameters. We propose an analytical short-time approximation for the time-dependent wave packets. We provide a theoretical basis for the occurrence of saturation in the rotational pumping. We present a formula to predict the time at which the maximum alignment occurs. We discuss the magnitude of the laser-induced alignment and we relate it to a theoretical upper limit.     

Inelastic collisions of cold polar molecules in nonparallel electric and magnetic fields

The authors present a detailed study of low-temperature collisions between CaD molecules and He atoms in superimposed electric and magnetic fields with arbitrary orientations. Electric fields do not interact with the electron spin of the molecules directly but modify their rotational structure and, consequently, the spin-rotation interactions. The authors examine molecular Stark and Zeeman energy levels as functions of the angle between the fields and show that rotating fields may induce and shift avoided crossings between the Zeeman levels of the rotationally ground and rotationally excited states of the molecule. The dynamics of molecular collisions are extremely sensitive to external fields near these avoided crossings and it is shown that molecular collisions may be controlled by varying both the strength and the relative orientation of the fields. The effects observed in this study are due to interactions of the isolated molecules with external fields so the conclusions should be relevant for collisions of molecules with other atoms or collisions of molecules with each other. This study demonstrates that electric fields may be used to enhance or suppress spin-rotation interactions in molecules. The spin-rotation interactions induce nonadiabatic couplings between states of different total spins in systems of two open-shell species and it is suggested that electric fields might be used for controlling nonadiabatic spin transitions and spin-forbidden chemical reactions of cold molecules in a magnetic trap.     

Abstractive dissociation of oxygen over Al(111): a nonadiabatic quantum model

The dissociation of oxygen on a clean aluminum surface is studied theoretically. A nonadiabatic quantum dynamical model is used, based on four electronically distinct potential energy surfaces characterized by the extent of charge transfer from the metal to the adsorbate. A flat surface approximation is used to reduce the computation complexity. The conservation of the helicopter angular momentum allows Boltzmann averaging of the outcome of the propagation of a three degrees of freedom wave function. The dissociation event is simulated by solving the time-dependent Schr?dinger equation for a period of 30 femtoseconds. As a function of incident kinetic energy, the dissociation yield follows the experimental trend. An attempt at simulation employing only the lowest adiabatic surface failed, qualitatively disagreeing with both experiment and nonadiabatic calculations. The final products, adsorptive dissociation and abstractive dissociation, are obtained by carrying out a semiclassical molecular dynamics simulation with surface hopping which describes the back charge transfer from an oxygen atom negative ion to the surface. The final adsorbed oxygen pair distribution compares well with experiment. By running the dynamical events backward in time, a correlation is established between the products and the initial conditions which lead to their production. Qualitative agreement is thus obtained with recent experiments that show suppression of abstraction by rotational excitation.     

Molecular dynamics simulation of the nematic liquid crystal phase in the presence of an intense magnetic field

The influence of an intense external field on the dynamics of the nematic liquid crystal phase is investigated using a molecular dynamics simulation for the Gay-Berne nematogen under isobaric-isothermal conditions. The molecular dynamics as a function of the second-rank orientational order parameter P<2> for a system consisting of a nematic liquid crystal in the presence of an intense magnetic field is compared with that of a similar system without the field. The translational motion of molecules is determined as a function of the translational diffusion coefficient tensor and the anisotropy and compared with the values predicted theoretically. The rotational dynamics of molecules is analyzed using the first- and the second-rank orientational time correlation functions. The translational diffusion coefficient parallel with respect to the director is constrained by the intense field, although the perpendicular one is decreased as the P<2> is increased, just as it is in the system without the field. However, no essential effect of the strong magnetic field is observed in the rotational molecular dynamics. Further, the rotational diffusion coefficient parallel with respect to the director obtained from the first-rank orientational time correlation function in the simulation is qualitatively in agreement with that in the real nematic liquid crystalline molecules. The P<2> dependence of the rotational diffusion coefficient for the system with the intense magnetic field shows a tendency similar to that for the system without the field.