Nonadiabatic Transitions Research Articles

Use of QM/MM Surface Hopping Simulations to Understand Thermally Activated Rare-Event Nonadiabatic Transitions in the Condensed Phase.

We implement a rare-event sampling scheme for quantifying the rate of thermally activated nonadiabatic transitions in the condensed phase. Our Quantum mechanics/molecular mechanics (QM/MM) methodology uses the recently developed Interface for NonAdiabatic QM/MM in Solvent (INAQS) package to interface an elementary electronic structure package and a popular open-source molecular dynamics software (GROMACS) to simulate an electron transfer event between two stationary ions in a solution of acetonitrile solvent molecules. Nonadiabatic effects are implemented through a surface hopping scheme, and our simulations allow further quantitative insight into the participation ratio of a solvent and the effect of ion separation distance as far as facilitating electron transfer. We also demonstrate that the standard gas-phase approaches for treating frustrated hops and velocity reversal must be refined when working in the condensed phase with many degrees of freedom. The code and methodology developed here can be easily expanded upon and modified to incorporate other systems and should provide a great deal of new insight into a wide variety of condensed phase nonadiabatic phenomena.

Beyond Cavity Born-Oppenheimer: On Nonadiabatic Coupling and Effective Ground State Hamiltonians in Vibro-Polaritonic Chemistry.

The emerging field of vibro-polaritonic chemistry studies the impact of light-matter hybrid states known as vibrational polaritons on chemical reactivity and molecular properties. Here, we discuss vibro-polaritonic chemistry from a quantum chemical perspective beyond the cavity Born-Oppenheimer (CBO) approximation and examine the role of electron-photon correlation in effective ground state Hamiltonians. We first quantitatively review ab initio vibro-polaritonic chemistry based on the molecular Pauli-Fierz Hamiltonian in dipole approximation and a vibrational strong coupling (VSC) Born-Huang expansion. We then derive nonadiabatic coupling elements arising from both "slow" nuclei and cavity modes compared to "fast" electrons via the generalized Hellmann-Feynman theorem, discuss their properties, and reevaluate the CBO approximation. In the second part, we introduce a crude VSC Born-Huang expansion based on adiabatic electronic states, which provides a foundation for widely employed effective Pauli-Fierz Hamiltonians in ground state vibro-polaritonic chemistry. Those do not strictly respect the CBO approximation but an alternative scheme, which we name crude CBO approximation. We argue that the crude CBO ground state misses electron-photon correlation relative to the CBO ground state due to neglected cavity-induced nonadiabatic transition dipole couplings to excited states. A perturbative connection between both ground state approximations is proposed, which identifies the crude CBO ground state as a first-order approximation to its CBO counterpart. We provide an illustrative numerical analysis of the cavity Shin-Metiu model with a focus on nonadiabatic coupling under VSC and electron-photon correlation effects on classical activation barriers. We finally discuss the potential shortcomings of the electron-polariton Hamiltonian when employed in the VSC regime.

Singlet and Triplet Pathways Determine the Thermal Z/E Isomerization of an Arylazopyrazole-Based Photoswitch.

Understanding the thermal isomerization mechanism of azobenzene derivatives is essential to designing photoswitches with tunable half-lives. Herein, we employ quantum chemical calculations, nonadiabatic transition state theory, and photosensitized experiments to unravel the thermal Z/E isomerization of a heteroaromatic azoswitch, the phenylazo-1,3,5-trimethylpyrazole. In contrast to the parent azobenzene, we predict two pathways to be operative at room temperature. One is a conventional ground-state reaction occurring via inversion of the aryl group, and the other is a nonadiabatic process involving intersystem crossing to the lowest-lying triplet state and back to the ground state, accompanied by a torsional motion around the azo bond. Our results illustrate that the fastest reaction rate is not controlled by the mechanism involving the lowest activation energy, but the size of the spin-orbit couplings at the crossing between the singlet and the triplet potential energy surfaces is also determinant. It is therefore mandatory to consider all of the multiple reaction pathways in azoswitches in order to predict experimental half-lives.

Multiple Dissociation Pathways in HNCO Decomposition Governed by Potential Energy Surface Topography.

The exquisite features of molecular photochemistry are key to any complete understanding of the chemical processes governed by potential energy surfaces (PESs). It is well established that multiple dissociation pathways relate to nonadiabatic transitions between multiple coupled PESs. However, little detail is known about how the single PES determines reaction outcomes. Here we perform detailed experiments on HNCO photodissociation, acquiring the state-specific correlations of the NH (a1Δ) and CO (X1Σ+) products. The experiments reveal a trimodal CO rotational distribution. Dynamics simulations based on a full-dimensional machine-learning-based PES of HNCO unveil three dissociation pathways exclusively occurring on the S1 excited electronic state. One pathway, following the minimum energy path (MEP) via the transition state, contributes to mild rotational excitation in CO, while the other two pathways deviating substantially from the MEP account for relatively cold and hot CO rotational state populations. These peculiar dynamics are unambiguously governed by the S1 state PES topography, i.e., a narrow acceptance cone in the vicinity of the transition state region. The dynamical picture shown in this work will serve as a textbook example illustrating the importance of the PES topography in molecular photochemistry.

Surface hopping, electron translation factors, electron rotation factors, momentum conservation, and size consistency.

For a system without spin-orbit coupling, the (i) nuclear plus electronic linear momentum and (ii) nuclear plus orbital electronic angular momentum are good quantum numbers. Thus, when a molecular system undergoes a nonadiabatic transition, there should be no change in the total linear or angular momentum. Now, the standard surface hopping algorithm ignores the electronic momentum and indirectly equates the momentum of the nuclear degrees of freedom to the total momentum. However, even with this simplification, the algorithm still does not conserve either the nuclear linear or the nuclear angular momenta. Here, we show that one way to address these failures is to dress the derivative couplings (i.e., the hopping directions) in two ways: (i) we disallow changes in the nuclear linear momentum by working in a translating basis (which is well known and leads to electron translation factors) and (ii) we disallow changes in the nuclear angular momentum by working in a basis that rotates around the center of mass [which is not well-known and leads to a novel, rotationally removable component of the derivative coupling that we will call electron rotation factors below, cf. Eq. (96)]. The present findings should be helpful in the short term as far as interpreting surface hopping calculations for singlet systems (without spin) and then developing the new surface hopping algorithm in the long term for systems where one cannot ignore the electronic orbital and/or spin angular momentum.

Efficiency of charge transfer in changing the dissociation dynamics of OD+ transients formed after the photo-fragmentation of D2O.

We present an investigation of the relaxation dynamics of deuterated water molecules after direct photo-double ionization at 61eV. We focus on the very rare D+ + O+ + D reaction channel in which the sequential fragmentation mechanisms were found to dominate the dynamics. Aided by theory, the state-selective formation and breakup of the transient OD+(a1Δ, b1Σ+) is traced, and the most likely dissociation path-OD+: a1Δ or b1Σ+ → A 3Π → X 3Σ- → B 3Σ--involving a combination of spin-orbit and non-adiabatic charge transfer transitions is determined. The multi-step transition probability of this complex transition sequence in the intermediate fragment ion is directly evaluated as a function of the energy of the transient OD+ above its lowest dissociation limit from the measured ratio of the D+ + O+ + D and competing D+ + D+ + O sequential fragmentation channels, which are measured simultaneously. Our coupled-channel time-dependent dynamics calculations reproduce the general trends of these multi-state relative transition rates toward the three-body fragmentation channels.

Analytic Evaluation of Nonadiabatic Couplings within the Complex Absorbing Potential Equation-of-Motion Coupled-Cluster Method.

We present the theory for the evaluation of nonadiabatic couplings (NACs) involving resonance states within the complex absorbing potential equation-of-motion coupled-cluster (CAP-EOM-CC) framework implemented within the singles and doubles approximation. Resonance states are embedded in the continuum and undergo rapid decay through autodetachment. In addition, nuclear motion can facilitate transitions between different resonances and between resonances and bound states. These nonadiabatic transitions affect the chemical fate of resonances and have distinct spectroscopic signatures. The NAC vector is a central quantity needed to model such effects. In the CAP-EOM-CC framework, resonance states are treated on the same footing as bound states. Using the example of fumaronitrile, which supports a bound radical anion and several anionic resonances, we analyze the NAC between bound states and pseudocontinuum states, between bound states and resonances, and between two resonances. We find that the NAC between a bound state and a resonance is nearly independent of the CAP strength and thus straightforward to evaluate, whereas the NAC between two resonance states or between a bound state and a pseudocontinuum state is more difficult to evaluate.

Electronic-state chaos, intramolecular electronic energy redistribution, and chemical bonding in persisting multidimensional nonadiabatic systems.

We study the chaotic, huge fluctuation of electronic state, resultant intramolecular energy redistribution, and strong chemical bonding surviving the fluctuation with exceedingly long lifetimes of highly excited boron clusters. Those excited states constitute densely quasi-degenerate state manifolds. The huge fluctuation is induced by persisting multidimensional nonadiabatic transitions among the states in the manifold. We clarify the mechanism of their coexistence and its physical significance. In doing so, we concentrate on two theoretical aspects. One is quantum chaos and energy randomization, which are to be directly extracted from the properties of the total electronic wavefunctions. The present dynamical chaos takes place through frequent transitions from adiabatic states to others, thereby making it very rare for the system to find dissociation channels. This phenomenon leads to the concept of what we call intramolecular nonadiabatic electronic-energy redistribution, which is an electronic-state generaliztion of the notion of intramolecular vibrational energy redistribution. The other aspect is about the peculiar chemical bonding. We investigate it with the energy natural orbitals (ENOs) to see what kind of theoretical structures lie behind the huge fluctuation. The ENO energy levels representing the highly excited states under study appear to have four robust layers. We show that the energy layers responsible for chaotic dynamics and those for chemical bonding are widely separated from each other, and only when an event of what we call "inter-layer crossing" happens to burst can the destruction of these robust energy layers occur, resulting in molecular dissociation. This crossing event happens only rarely because of the large energy gaps between the ENO layers. It is shown that the layers of high energy composed of complex-valued ENOs induce the turbulent flow of electrons and electronic-energy in the cluster. In addition, the random and fast time-oscillations of those high energy ENOs serve as a random force on the nuclear dynamics, which can work to prevent a concentration of high nuclear kinetic energy in the dissociation channels.

Zeno and Anti-Zeno Effects in Nonadiabatic Molecular Dynamics.

Decoherence plays an important role in nonadiabatic (NA) molecular dynamics (MD) simulations because it provides a physical mechanism for trajectory hopping and can alter transition rates by orders of magnitude. Generally, decoherence effects slow quantum transitions, as exemplified by the quantum Zeno effect: in the limit of infinitely fast decoherence, the transitions stop. If the measurements are not sufficiently frequent, an opposite quantum anti-Zeno effect occurs, in which the transitions are accelerated with faster decoherence. Using two common NA-MD approaches, fewest switches surface hopping and decoherence-induced surface hopping, combined with analytic examination, we demonstrate that including decoherence into NA-MD slows down NA transitions; however, many realistic systems operate in the anti-Zeno regime. Therefore, it is important that NA-MD methods describe both Zeno and anti-Zeno effects. Numerical simulations of charge trapping and relaxation in graphitic carbon nitride suggest that time-dependent NA Hamiltonians encountered in realistic systems produce robust results with respect to errors in the decoherence time, a favorable feature for NA-MD simulations.

Time rescaling of nonadiabatic transitions

Applying time-dependent driving is a basic way of quantum control. Driven systems show various dynamics as its time scale is changed due to the different amount of nonadiabatic transitions. The fast-forward scaling theory enables us to observe slow (or fast) time-scale dynamics during moderate time by applying additional driving. Here we discuss its application to nonadiabatic transitions. We derive mathematical expression of additional driving and also find a formula for calculating it. Moreover, we point out relation between the fast-forward scaling theory for nonadiabatic transitions and shortcuts to adiabaticity by counterdiabatic driving.

Theoretical study of nonadiabatic hydrogen atom scattering dynamics on metal surfaces using the hierarchical equations of motion method.

Hydrogen atom scattering on metal surfaces is investigated based on a simplified Newns-Anderson model. Both the nuclear and electronic degrees of freedom are treated quantum mechanically. By partitioning all the surface electronic states as the bath, the hierarchical equations of motion method for the fermionic bath is employed to simulate the scattering dynamics. It is found that, with a reasonable set of parameters, the main features of the recent experimental studies of hydrogen atom scattering on metal surfaces can be reproduced. Vibrational states on the chemisorption state whose energies are close to the incident energy are found to play an important role, and the scattering process is dominated by a single-pass electronic transition forth and back between the diabatic physisorption and chemisorption states. Further study on the effects of the atom-surface coupling strength reveals that, upon increasing the atom-surface coupling strength, the scattering mechanism changes from typical nonadiabatic transitions to dynamics in the electronic friction regime.

Probing C-I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy.

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268nm UV light-induced dynamics in 2-iodothiophene (C4H3IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by abinitio molecular dynamics simulations. Gradual absorption changes up to ∼0.5 to 1ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C-I bond, enabling a nonadiabatic transition to a dissociative π→σC-I* state. For the subsequent fragmentation, a relatively narrow bond-length region along the C-I stretch coordinate between 230 and 280pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

Differential Shannon Entropies Characterizing Electron-Nuclear Dynamics and Correlation: Momentum-Space Versus Coordinate-Space Wave Packet Motion.

We calculate differential Shannon entropies derived from time-dependent coordinate-space and momentum-space probability densities. This is performed for a prototype system of a coupled electron-nuclear motion. Two situations are considered, where one is a Born-Oppenheimer adiabatic dynamics, and the other is a diabatic motion involving strong non-adiabatic transitions. The information about coordinate- and momentum-space dynamics derived from the total and single-particle entropies is discussed and interpreted with the help of analytical models. From the entropies, we derive mutual information, which is a measure for the electron-nuclear correlation. In the adiabatic case, it is found that such correlations are manifested differently in coordinate- and momentum space. For the diabatic dynamics, we show that it is possible to decompose the entropies into state-specific contributions.

Efficient White Light Emission of 0D Lead-Free Indium-Based Halide Perovskite and the Intermediate State Promotion Mechanism of the Nonadiabatic Transition to Self-Trapped Exciton via Antimony(III) Cation Doping

Efficient White Light Emission of 0D Lead-Free Indium-Based Halide Perovskite and the Intermediate State Promotion Mechanism of the Nonadiabatic Transition to Self-Trapped Exciton via Antimony(III) Cation Doping

Driving forces for ultrafast laser-induced sp2 to sp3 structural transformation in graphite

Understanding the microscopic mechanism of photoinduced sp2-to-sp3 structural transformation in graphite is a scientific challenge with great importance. Here, the ultrafast dynamics and characteristics of laser-induced structural transformation in graphite are revealed by non-adiabatic quantum dynamic simulations. Under laser irradiation, graphite undergoes an interlayer compression and sliding stage, followed by a key period of intralayer buckling and interlayer bonding to form an intermediate sp2-sp3 hybrid structure, before completing the full transformation to hexagonal diamond. The process is driven by the cooperation of charge carrier multiplication and selective phonon excitations through electron-phonon interactions, in which photoexcited hot electrons scattered into unoccupied high-energy conduction bands play a key role in the introduction of in-plane instability in graphite. This work identifies a photoinduced non-adiabatic transition pathway from graphite to diamond and shows far-reaching implications for designing optically controlled structural phase transition in materials.

Photonic Floquet Landau-Zener tunneling and temporal beam splitters.

Landau-Zener tunneling (LZT), i.e., the nonadiabatic transition under strong parameter driving in multilevel systems, is ubiquitous in physics, providing a powerful tool for coherent wave control both in quantum and classical systems. While previous works mainly focus on LZT between two energy bands in time-invariant crystals, here, we construct synthetic time-periodic temporal lattices from two coupled fiber loops and demonstrate dc- and ac-driven LZTs between periodic Floquet bands. We show that dc- and ac-driven LZTs display distinctive tunneling and interference characteristics, which can be harnessed to realize fully reconfigurable LZT beam splitter arrangements. As a potential application to signal processing, we realize a 4-bit temporal beam encoder for classical light pulses using a reconfigurable LZT beam splitter network. Our work introduces and experimentally demonstrates a new class of reconfigurable linear optics circuits harnessing Floquet LZT, which may find versatile applications in temporal beam control, signal processing, quantum simulations, and information processing.

Lead-Free Zero-Dimensional Zn-Based Metal Halides of Highly Efficient Blue Luminescence from Self-Trapping Exciton

Recently, hybrid metal halides have received great attention in the field of solid-state lighting because of their diverse structures and excellent photoluminescence properties. In this work, we first reported two hybrid zinc-based metal halides with zero-dimensional structures, (BMPP)2ZnBr4 and (TBA)2ZnBr4, which exhibited broadband emission with large Stokes shifts. Notably, the highest photoluminescence quantum yield of 59.76% was observed. Additionally, the luminescence mechanism of metal halides was investigated by using time-resolved femtosecond transient absorption experiments. A broad excited-state absorption platform with the tendency of slowly decaying was shown in the detection range, demonstrating that after the electrons were excited to the excited state, the free excitons underwent a nonadiabatic transition to self-trapped excitons and went through a radiation recombination process to the ground state. A blue-light-emitting diode could be easily obtained by coating (BMPP)2ZnBr4 on a GaN chip, which indicated that it has good competitiveness in the application of solid-state lighting devices.

Nonadiabatic transition probabilities for quantum systems in electromagnetic fields: Dephasing and population relaxation due to contact with a bath.

We contrast Dirac's theory of transition probabilities and the theory of nonadiabatic transition probabilities, applied to a perturbed system that is coupled to a bath. In Dirac's analysis, the presence of an excited state |k0⟩ in the time-dependent wave function constitutes a transition. In the nonadiabatic theory, a transition occurs when the wave function develops a term that is not adiabatically connected to the initial state. Landau and Lifshitz separated Dirac's excited-state coefficients into a term that follows the adiabatic theorem of Born and Fock and a nonadiabatic term that represents excitation across an energy gap. If the system remains coherent, the two approaches are equivalent. However, differences between the two approaches arise when coupling to a bath causes dephasing, a situation that was not treated by Dirac. For two-level model systems in static electric fields, we add relaxation terms to the Liouville equation for the time derivative of the density matrix. We contrast the results obtained from the two theories. In the analysis based on Dirac's transition probabilities, the steady state of the system is not an equilibrium state; also, the steady-state population ρkk,s increases with increasing strength of the perturbation and its value depends on the dephasing time T2. In the nonadiabatic theory, the system evolves to the thermal equilibrium with the bath. The difference is not simply due to the choice of basis because the difference remains when the results are transformed to a common basis.

Majorana’s approach to nonadiabatic transitions validates the adiabatic-impulse approximation

The approach by Ettore Majorana for non-adiabatic transitions between two quasi-crossing levels is revisited and significantly extended. We rederive the transition probability, known as the Landau–Zener–Stückelberg–Majorana formula, and introduce Majorana’s approach to modern readers. This result, typically referred as the Landau–Zener formula, was published by Majorana before Landau, Zener and Stückelberg. Moreover, we go well beyond previous results and we now obtain the full wave function, including its phase, which is important nowadays for quantum control and quantum information. The asymptotic wave function correctly describes the dynamics away from the avoided-level crossing, while it has limited accuracy in that region.

Proximity‐Encirclement of Multiple Exceptional Points in Non‐Hermitian Photonic Waveguides

AbstractConventionally, dynamical encirclement of exceptional points in non‐Hermitian systems is known to manifest a counterintuitive chiral state conversion. However, the prerequisite of such traits enclosing an exceptional point is broken when only encircling its proximity, preserving a still chiral switching. Research on the proximity‐encirclement in multistate systems is lacking. In this paper, a photonic‐waveguide‐array non‐Hermitian system is proposed to investigate the dynamics by encircling two exceptional points or their proximity. A series of encircling trajectories defined by the parametric equations are designed to steer the evolution of photonic modes in waveguides. The wave propagating along the waveguides is also simulated to capture this non‐Hermitian physics. The chiral behavior in proximity‐encirclement contrasts with the familiar encirclement of one exceptional point and exhibits the unexpected occurrence of nonadiabatic transitions. Furthermore, if two exceptional points are sufficiently encircled, the system will evolve to a stable final state earlier, as a symbol of the occurrence of the nonadiabatic transition. Such novel chiral conversion is maintained only if the encircling trajectories are located at adequate proximity.