Dissipative dynamics of laser induced nonadiabatic molecular alignment

Nonadiabatic alignment induced by short, moderately intense laser pulses in molecules coupled to dissipative environments is studied within a nonperturbative density matrix theory. We focus primarily on exploring and extending a recently proposed approach [Phys. Rev. Lett. 95, 113001 (2005)], wherein nonadiabatic laser alignment is used as a coherence spectroscopy that probes the dissipative properties of the solvent. To that end we apply the method to several molecular collision systems that exhibit sufficiently varied behavior to represent a broad variety of chemical environments. These include molecules in low temperature gas jets, in room temperature gas cells, and in dense liquids. We examine also the possibility of prolonging the duration of the field free (post-pulse) alignment in dissipative media by a proper choice of the system parameters.     

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.     

Approach to nonadiabatic transitions by density matrix evolution and molecular dynamics simulations

Many biological processes are characterized by an essentially quantum dynamical event, such as electron or proton transfer, in a complex classical environment. To treat such processes properly by computer simulation, allowing nonadiabatic transitions involving excited states, we recently developed a density matrix evolution (DME ) method [H. J. C. Berendsen and J. Mavri, J. Phys. Chem, 97 , 13464 (1993)] which simulates the dynamics of quantum systems embedded in a classical environment. The formalism of the method is presented and an overview of the applications ranging from collisions of a quantum harmonic oscillator with noble gas atoms to proton tunneling in a double-well hydrogen bond is given. The methodology for treatment of proton-transfer processes with inclusion of excited states is presented. Future application of the method on biologically interesting processes, such as proton transfer in enzymatic reactions, is discussed. © 1996 John Wiley & Sons, Inc.     

Muffin-tin orbital Wannier-like functions for insulators and metals.

Herein, we outline a method that is able to generate truly minimal basis sets that accurately describe either a group of bands, a band, or even just the occupied part of a band. These basis sets are the so-called NMTOs, muffin-tin orbitals of order N. For an isolated set of bands, symmetrical orthonormalization of the NMTOs yields a set of Wannier functions that are atom-centered and localized by construction. They are not necessarily maximally localized, but may be transformed into those Wannier functions. For bands that overlap others, Wannier-like functions can be generated. It is shown that NMTOs give a chemical understanding of an extended system. In particular, orbitals for the pi and sigma bands in an insulator, boron nitride, and a semimetal, graphite, will be considered. In addition, we illustrate that it is possible to obtain Wannier-like functions for only the occupied states in a metallic system by generating NMTOs for cesium. Finally, we visualize the pressure-induced s-->d transition.     

Non-Born-Oppenheimer path in anti-Hermitian dynamics for nonadiabatic transitions

A serious difficulty in the semiclassical Ehrenfest theory for nonadiabatic transitions is that a path passing across the avoided crossing is forced to run on a potential averaged over comprising adiabatic potential surfaces that commit the avoided crossing. Therefore once a path passes through the crossing region, it immediately becomes incompatible with the standard view of "classical trajectory" running on an adiabatic surface. This casts a fundamental question to the theoretical structure of chemical dynamics. In this paper, we propose a non-Born-Oppenheimer path that is generated by an anti-Hermitian Hamiltonian, whose complex-valued eigenenergies can cross in their real parts and avoid crossing in the imaginary parts in the vicinity of the nonadiabatic transition region. We discuss the properties of this non-Born-Oppenheimer path and thereby show its compatibility with the Born-Oppenheimer classical trajectories. This theory not only allows the geometrical branching of the paths but gives the nonadiabatic transition amplitudes and quantum phases along the generated paths.     

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.     

Incorporation of nonadiabatic transition into wave-packet dynamics

Nonadiabatic wave-packet dynamics is factorized into purely adiabatic propagation and instantaneous localized nonadiabatic transition. A general formula is derived for the quantum-mechanical local nonadiabatic operator which is implemented within the framework of the R-matrix method. The operator can be used for incorporating the nonadiabatic transition in semiclassical wave-packet dynamics.     

Quantum nonadiabatic dynamics of hydrogen exchange reactions

In continuation of our earlier effort to understand the nonadiabatic coupling effects in the prototypical H + H₂ exchange reaction [Jayachander Rao et al. Chem. Phys. 333 (2007) 135], we present here further quantum dynamical investigations on its isotopic variants. The present work also corrects a technical scaling error occurred in our previous studies on the H + HD reaction. Initial state-selected total reaction cross sections and Boltzmann averaged thermal rate constants are calculated with the aid of a time-dependent wave packet approach employing the double many body expansion potential energy surfaces of the system. The theoretical results are compared with the experimental and other theoretical data whenever available. The results re-establish our earlier conclusion, on a more general perspective, that the electronic nonadiabatic effects are negligible on the important quantum dynamical observables of these reactive systems reported here.     

Ab initio nonadiabatic molecular dynamics of the ultrafast electron injection across the alizarin-TiO2 interface

The observed 6-fs photoinduced electron transfer (ET) from the alizarin chromophore into the TiO2 surface is investigated by ab initio nonadiabatic (NA) molecular dynamics in real time and at the atomistic level of detail. The system derives from the dye-sensitized semiconductor Gr?tzel cell and addresses the problems of an organic/inorganic interface that are commonly encountered in photovoltaics, photochemistry, and molecular electronics. In contrast to the typical Gr?tzel cell systems, where molecular donors are in resonance with a high density of semiconductor acceptor states, TiO2 sensitized with alizarin presents a novel case in which the molecular photoexcited state is at the edge of the conduction band (CB). The high level ab initio analysis of the optical absorption spectrum supports this observation. Thermal fluctuations of atomic coordinates are particularly important both in generating a nonuniform distribution of photoexcited states and in driving the ET process. The NA simulation resolves the controversy regarding the origin of the ultrafast ET by showing that although ultrafast transfer is possible with the NA mechanism, it proceeds mostly adiabatically in the alizarin-TiO2 system. The simulation indicates that the electron is injected into a localized surface state within 8 fs and spreads into the bulk on a 100-fs or longer time scale. The molecular architecture seen in the alizarin-TiO2 system permits efficient electron injection into the edge of the CB by an adiabatic mechanism without the energy loss associated with injection high into the CB by a NA process.     

Theoretical insight into the aggregation induced emission phenomena of diphenyldibenzofulvene: a nonadiabatic molecular dynamics study

The diphenyldibenzofulvene (DPDBF) molecule appears in two forms: ring open and ring closed. The former fluoresces weakly in solution, but it becomes strongly emissive in the solid phase, exhibiting an exotic aggregation-induced emission phenomenon. The latter presents a normal aggregation quenching phenomenon, as is expected. We implement nonadiabatic molecular dynamics based on the combination of time-dependent Kohn-Sham (TDKS) and density functional tight binding (DFTB) methods with Tully's fewest switches surface hopping algorithm to investigate the excited state nonradiative decay processes. From the analysis of the nonadiabatic coupling vectors, it is found that the low frequency twisting motion in the ring open DPDBF couples strongly with the electronic excitation and dissipates the energy efficiently. While in the closed form, such motion is blocked by a chemical bond. This leads to the nonradiative decay rate for the open form (1.4 ps) becoming much faster than the closed form (24.5 ps). It is expected that, in the solid state, the low frequency motion of the open form will be hindered and the energy dissipation pathway by nonradiative decay will be slowed, presenting a remarkable aggregation enhanced emission phenomenon.     

A coherent state approach to semiclassical nonadiabatic dynamics

A semiclassical (SC) approximation to the quantum mechanical propagator for nonadiabatic systems is derived. Our derivation starts with an exact path integral expression that uses canonical coherent states for the nuclear degrees of freedom and spin coherent states for the electronic degrees of freedom. A stationary path approximation (SPA) is then applied to the path integral to obtain the SC approximation. The SPA results in complex classical trajectories of both nuclear and electronic degrees of freedom and a double ended boundary condition. The root search problem is solved using the previously proposed "real trajectory local search" algorithm. The SC approximation is tested on three simple one dimensional two-state systems proposed by Tully [J. Chem. Phys. 93, 1061 (1990)], and the SC results are compared to Ehrenfest and surface hopping predictions. Excellent agreement with quantum results is reached when the SC trajectory is far away from caustics. We discuss the origin of caustics in this SC formalism and the strengths and weaknesses of this approach.     

Manipulating nonadiabatic conical intersection dynamics by optical cavities

Optical cavities hold great promise to manipulate and control the photochemistry of molecules. We demonstrate how molecular photochemical processes can be manipulated by strong light–matter coupling. For a molecule with an inherent conical intersection, optical cavities can induce significant changes in the nonadiabatic dynamics by either splitting the pristine conical intersections into two novel polaritonic conical intersections or by creating light-induced avoided crossings in the polaritonic surfaces. This is demonstrated by exact real-time quantum dynamics simulations of a three-state two-mode model of pyrazine strongly coupled to a single cavity photon mode. We further explore the effects of external environments through dissipative polaritonic dynamics computed using the hierarchical equation of motion method. We find that cavity-controlled photochemistry can be immune to external environments. We also demonstrate that the polariton-induced changes in the dynamics can be monitored by transient absorption spectroscopy.

Optical cavities hold great promise to manipulate and control the photochemistry of molecules.     

Phonon thermal conductivity in nanolaminated composite metals via molecular dynamics

We use nonequilibrium molecular dynamics to characterize the phonon contribution to thermal conduction of Al nanostructures and the role of interfaces in metallic nanocomposites. We characterize the lattice thermal conductivity of pure Al samples as a function of size and temperature from which we obtain, using kinetic theory, the temperature dependence of the phonon mean free path. We also calculated the thermal conductivity of AlAl* and AlNi nanolaminate composites (where Al* differs from Al only in its mass) for various periodic sizes and compositions as well as the associated interfacial thermal resistivities (ITRs). We find that simple, additive models provide good estimates of the thermal conductivities of the nanocomposites in terms of those of the individual components and interfaces if size effects on the behavior of the individual components are considered. The additive models provide important insight to the decrease in thermal conductivity of the nanolaminates as their periodicity (thickness of a bilayer) is reduced to a size comparable with the phonon mean free path and break down when this characteristic size is reduced further. At this point the system can be regarded as homogeneous and the conductivity increases with decreasing periodicity of the laminates. We also observe that the ITR depends on the direction of the heat flux; this is the first molecular level characterization of such thermal diode behavior in a realistic three dimensional material.     

Semiclassical nonadiabatic dynamics using a mixed wave-function representation

Nonadiabatic effects in quantum dynamics are described using a mixed polar/coordinate space representation of the wave function. The polar part evolves on dynamically determined potential surfaces that have diabatic and adiabatic potentials as limiting cases of weak localized and strong extended diabatic couplings. The coordinate space part, generalized to a matrix form, describes transitions between the surfaces. Choice of the effective potentials for the polar part and partitioning of the wave function enables one to represent the total wave function in terms of smooth components that can be accurately propagated semiclassically using the approximate quantum potential and small basis sets. Examples are given for two-state one-dimensional problems that model chemical reactions that demonstrate the capabilities of the method for various regimes of nonadiabatic dynamics.     

Matching-pursuit/split-operator Fourier-transform simulations of nonadiabatic quantum dynamics

A rigorous and practical approach for simulations of nonadiabatic quantum dynamics is introduced. The algorithm involves a natural extension of the matching-pursuitsplit-operator Fourier-transform (MPSOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] recently developed for simulations of adiabatic quantum dynamics in multidimensional systems. The MPSOFT propagation scheme, extended to nonadiabatic dynamics, recursively applies the time-evolution operator as defined by the standard perturbation expansion to first-, or second-order, accuracy. The expansion is implemented in dynamically adaptive coherent-state representations, generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting propagation method are demonstrated as applied to the canonical model systems introduced by Tully for testing simulations of dual curve-crossing nonadiabatic dynamics.     

A diabatic representation including both valence nonadiabatic interactions and spin-orbit effects for reaction dynamics

A diabatic representation is convenient in the study of electronically nonadiabatic chemical reactions because the diabatic energies and couplings are smooth functions of the nuclear coordinates and the couplings are scalar quantities. A method called the fourfold way was devised in our group to generate diabatic representations for spin-free electronic states. One drawback of diabatic states computed from the spin-free Hamiltonian, called a valence diabatic representation, for systems in which spin-orbit coupling cannot be ignored is that the couplings between the states are not zero in asymptotic regions, leading to difficulties in the calculation of reaction probabilities and other properties by semiclassical dynamics methods. Here we report an extension of the fourfold way to construct diabatic representations suitable for spin-coupled systems. In this article we formulate the method for the case of even-electron systems that yield pairs of fragments with doublet spin multiplicity. For this type of system, we introduce the further simplification of calculating the triplet diabatic energies in terms of the singlet diabatic energies via Slater's rules and assuming constant ratios of Coulomb to exchange integrals. Furthermore, the valence diabatic couplings in the triplet manifold are taken equal to the singlet ones. An important feature of the method is the introduction of scaling functions, as they allow one to deal with multibond reactions without having to include high-energy diabatic states. The global transformation matrix to the new diabatic representation, called the spin-valence diabatic representation, is constructed as the product of channel-specific transformation matrices, each one taken as the product of an asymptotic transformation matrix and a scaling function that depends on ratios of the spin-orbit splitting and the valence splittings. Thus the underlying basis functions are recoupled into suitable diabatic basis functions in a manner that provides a multibond generalization of the switch between Hund's cases in diatomic spectroscopy. The spin-orbit matrix elements in this representation are taken equal to their atomic values times a scaling function that depends on the internuclear distances. The spin-valence diabatic potential energy matrix is suitable for semiclassical dynamics simulations. Diagonalization of this matrix produces the spin-coupled adiabatic energies. For the sake of illustration, diabatic potential energy matrices are constructed along bond-fission coordinates for the HBr and the BrCH(2)Cl molecules. Comparison of the spin-coupled adiabatic energies obtained from the spin-valence diabatics with those obtained by ab initio calculations with geometry-dependent spin-orbit matrix elements shows that the new method is sufficiently accurate for practical purposes. The method formulated here should be most useful for systems with a large number of atoms, especially heavy atoms, and/or a large number of spin-coupled electronic states.     

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.