Nonadiabatic Rate Research Articles

Highly Twisted N,N-Dialkylamines as a Design Strategy to Tune Simple Aromatic Hydrocarbons as Steric Environment-Sensitive Fluorophores.

The steric-environment sensitivity of fluorescence of 9,10-bis(N,N-dialkylamino)anthracenes (BDAAs) was studied experimentally and theoretically. A new design strategy to tune simple aromatic hydrocarbons as efficient aggregation-induced emission (AIE) luminogens and molecular rotors is proposed. For a variety of BDAAs, prominent Stokes shifts and efficient solid-state fluorescence were observed. Calculations on BDAA-methyl suggested that in the ground state (S0) conformations, the pyramidal amine groups are highly twisted, so that their lone-pair orbitals cannot conjugate with the anthracene π orbitals. Fluorescence takes place from the S1 minima, in which one or both amine groups are planarized. The stability of the S1 excited state minima as well as destabilization of the S0 state is the origin of large Stokes shift. Experimental measurement of the nonadiabatic transition rate suggests that para disubstitution by dialkylamino (or strongly electron-donating) groups is a key for fast internal conversion. Minimum energy conical intersection (MECI) between S1 and S0 states was found to have a Dewar-benzene like structure. Although this can be reached efficiently in liquid phase for fast internal conversion, a large amplitude motion is required to reach this MECI, which is prohibited in the solid state and caused efficient AIE. This strategy is used to find experimentally that naphthalene analogues are also efficient AIE luminogens. The flexibility of alkyl chains on amino groups is also found to be important for allowed charge-transfer transition. Thus, three points [(1) highly twisted N,N-dialkylamines, (2) substitution at the para positions, (3) with flexible alkyl groups] were proposed for activation of small aromatic hydrocarbons.

Direct ab Initio (Meta-)Surface-Hopping Dynamics.

Tractable methods for studying the molecular dynamics of chemical processes driven by electronic nonadiabaticity are highly sought after to provide insight into, for example, photochemical reaction mechanisms, molecular collisions, and thermalized electronic band structures. Starting from the time-dependent Schrödinger equation for a many-body system, a direct ab initio trajectory surface-hopping (TSH) method relying on an analytical treatment of nonadiabatic couplings between electronic states is developed in this work. An approach that combines time-dependent perturbation theory and explicit time evolution via TSH to expedite calculation of nonadiabatic transition rates, namely, meta-surface-hopping dynamics, is presented, and an extrapolatory approach using time-dependent perturbation theory for recovering unbiased transition rates is assessed. The meta-surface-hopping method is applied to the problem of estimating nonradiative relaxation rates of a photoexcited iminium ion, CH₂NH₂⁺, and evidence for internal consistency of the combined dynamics/perturbation theory approach is presented.

Mean field ring polymer molecular dynamics for electronically nonadiabatic reaction rates.

We present a mean field ring polymer molecular dynamics method to calculate the rate of electron transfer (ET) in multi-state, multi-electron condensed-phase processes. Our approach involves calculating a transition state theory (TST) estimate to the rate using an exact path integral in discrete electronic states and continuous Cartesian nuclear coordinates. A dynamic recrossing correction to the TST rate is then obtained from real-time dynamics simulations using mean field ring polymer molecular dynamics. We employ two different reaction coordinates in our simulations and show that, despite the use of mean field dynamics, the use of an accurate dividing surface to compute TST rates allows us to achieve remarkable agreement with Fermi's golden rule rates for nonadiabatic ET in the normal regime of Marcus theory. Further, we show that using a reaction coordinate based on electronic state populations allows us to capture the turnover in rates for ET in the Marcus inverted regime.

Nonadiabatic rate constants for proton transfer and proton-coupled electron transfer reactions in solution: Effects of quadratic term in the vibronic coupling expansion.

Rate constant expressions for vibronically nonadiabatic proton transfer and proton-coupled electron transfer reactions are presented and analyzed. The regimes covered include electronically adiabatic and nonadiabatic reactions, as well as high-frequency and low-frequency proton donor-acceptor vibrational modes. These rate constants differ from previous rate constants derived with the cumulant expansion approach in that the logarithmic expansion of the vibronic coupling in terms of the proton donor-acceptor distance includes a quadratic as well as a linear term. The analysis illustrates that inclusion of this quadratic term in the framework of the cumulant expansion framework may significantly impact the rate constants at high temperatures for proton transfer interfaces with soft proton donor-acceptor modes that are associated with small force constants and weak hydrogen bonds. The effects of the quadratic term may also become significant in these regimes when using the vibronic coupling expansion in conjunction with a thermal averaging procedure for calculating the rate constant. In this case, however, the expansion of the coupling can be avoided entirely by calculating the couplings explicitly for the range of proton donor-acceptor distances sampled. The effects of the quadratic term for weak hydrogen-bonding systems are less significant for more physically realistic models that prevent the sampling of unphysical short proton donor-acceptor distances. Additionally, the rigorous relation between the cumulant expansion and thermal averaging approaches is clarified. In particular, the cumulant expansion rate constant includes effects from dynamical interference between the proton donor-acceptor and solvent motions and becomes equivalent to the thermally averaged rate constant when these dynamical effects are neglected. This analysis identifies the regimes in which each rate constant expression is valid and thus will be important for future applications to proton transfer and proton-coupled electron transfer in chemical and biological processes.

Nonadiabatic Electron Dynamics in Orthogonal Two-Color Laser Fields with Comparable Intensities.

We theoretically investigate the nonadiabatic subcycle electron dynamics in orthogonally polarized two-color laser fields with comparable intensities. The photoelectron dynamics is simulated by exact solution to the 3D time-dependent Schrödinger equation, and also by two other semiclassical methods, i.e., the quantum trajectory MonteCarlo simulation and the Coulomb-corrected strong field approximation. Through these methods, we identify the underlying mechanisms of the subcycle electron dynamics and find that both the nonadiabatic effects and the Coulomb potential play very important roles. The contribution of the nonadiabatic effects manifest in two aspects, i.e., the nonadiabatic ionization rate and the nonzero initial velocities at the tunneling exit. The Coulomb potential has a different impact on the electrons' trajectories for different relative phases between the two pulses.

Laser-Driven Hole Trapping in a Ge/Si Core–Shell Nanocrystal: An Atomistic Configuration Interaction Perspective

This contribution reports on the first example of laser-driven charge carrier confinement in a solid state quantum dot investigated using a fully atomistic, correlated many-electron ansatz. Specifically, a Ge/Si model nanocrystal is designed to retain the main structural characteristics and excitonic properties of experimentally observed self-assembled pyramidal heterostructures. State-selective laser excitations yielding hole confinement in the Ge nanostructure are simulated using the reduced density matrix variant of the time-dependent configuration interaction method (ρ-TDCI). The degree of carrier localization in the quantum dot is determined by analyzing the correlated one-electron densities from configuration interaction electronic states at a singles level. Additionally, dissipation and pure dephasing are included to treat the coupling of the local core–shell structure with the vibrations of the surrounding silicon matrix. For this purpose, a new microscopic model for these nonadiabatic coupling rates is derived. The results reveal that, despite the presence of dissipation, charge carriers can be efficiently confined by localized optical excitations in the model Ge/Si quantum dot to create long-lived, large permanent dipoles in the nanocrystal.

Surface hopping, transition state theory, and decoherence. II. Thermal rate constants and detailed balance.

We investigate a simple approach to compute a non-adiabatic thermal rate constant using the fewest switches surface hopping (FSSH) dynamics. We study the effects of both decoherence (using our augmented-FSSH (A-FSSH) algorithm) and forbidden hops over a large range of parameters, including high and low friction regimes, and weak and strong electronic coupling regimes. Furthermore, when possible, we benchmark our results against exact hierarchy equations of motion results, where we usually find a maximum error of roughly a factor of two (at reasonably large temperatures). In agreement with Hammes-Schiffer and Tully, we find that a merger of transition state theory and surface hopping can be both accurate and efficient when performed correctly. We further show that detailed balance is followed approximately by A-FSSH dynamics.

Electron transfer, decoherence, and protein dynamics: insights from atomistic simulations.

Electron transfer in biological systems drives the processes of life. From cellular respiration to photosynthesis and enzymatic catalysis, electron transfers (ET) are chemical processes on which essential biological functions rely. Over the last 40 years, scientists have sought understanding of how these essential processes function in biology. One important breakthrough was the discovery that Marcus theory (MT) of electron transfer is applicable to biological systems. Chemists have experimentally collected both the reorganization energies (λ) and the driving forces (ΔG°), two parameters of Marcus theory, for a large variety of ET processes in proteins. At the same time, theoretical chemists have developed computational approaches that rely on molecular dynamics and quantum chemistry calculations to access numerical estimates of λ and ΔG°. Yet another crucial piece in determining the rate of an electron transfer is the electronic coupling between the initial and final electronic wave functions. This is an important prefactor in the nonadiabatic rate expression, since it reflects the probability that an electron tunnels from the electron donor to the acceptor through the intervening medium. The fact that a protein matrix supports electron tunneling much more efficiently than vacuum is now well documented, both experimentally and theoretically. Meanwhile, many chemists have provided examples of the rich physical chemistry that can be induced by protein dynamics. This Account describes our studies of the dynamical effects on electron tunneling. We present our analysis of two examples of natural biological systems through MD simulations and tunneling pathway analyses. Through these examples, we show that protein dynamics sustain efficient tunneling. Second, we introduce two time scales: τcoh and τFC. The former characterizes how fast the electronic coupling varies with nuclear vibrations (which cause dephasing). The latter reflects the time taken by the system to leave the crossing region. In the framework of open quantum systems, τFC is a short time approximation of the characteristic decoherence time of the electronic subsystem in interaction with its nuclear environment. The comparison of the respective values of τcoh and τFC allows us to probe the occurrence of non-Condon effects. We use ab initio MD simulations to analyze how decoherence appears in several biological cofactors. We conclude that we cannot account for its order of magnitude by considering only the atoms or bonds directly concerned with the transfer. Decoherence results from contributions from all atoms of the system appearing with a time delay that increases with the distance from the primarily concerned atoms or bonds. The delay and magnitude of the contributions depend on the chemical nature of the system. Finally, we present recent developments based on constrained DFT for efficient and accurate evaluations of the electronic coupling in ab initio MD simulations. These are promising methods to study the subtle fluctuations of the electronic coupling and the mechanisms of electronic decoherence in biological systems.

Field-cycle-resolved photoionization in solids.

The Keldysh theory of photoionization in a solid dielectric is generalized to the case of arbitrarily short driving pulses of arbitrary pulse shape. We derive a closed-form solution for the nonadiabatic ionization rate in a transparent solid with a periodic dispersion relation, which reveals ultrafast ionization dynamics within the field cycle and recovers the key results of the Keldysh theory in the appropriate limiting regimes.

Probing Nonadiabaticity in the Proton-Coupled Electron Transfer Reaction Catalyzed by Soybean Lipoxygenase.

Proton-coupled electron transfer (PCET) plays a vital role in many biological and chemical processes. PCET rate constant expressions are available for various well-defined regimes, and determining which expression is appropriate for a given system is essential for reliable modeling. Quantitative diagnostics have been devised to characterize the vibronic nonadiabaticity between the electron–proton quantum subsystem and the classical nuclei, as well as the electron–proton nonadiabaticity between the electrons and proton(s) within the quantum subsystem. Herein these diagnostics are applied to a model of the active site of the enzyme soybean lipoxygenase, which catalyzes a PCET reaction that exhibits unusually high deuterium kinetic isotope effects at room temperature. Both semiclassical and electronic charge density diagnostics illustrate vibronic and electron–proton nonadiabaticity for this PCET reaction, supporting the use of the Golden rule nonadiabatic rate constant expression with a specific form of the vibronic coupling. This type of characterization will be useful for theoretical modeling of a broad range of PCET processes.

Non-oscillatory flux correlation functions for efficient nonadiabatic rate theory.

There is currently much interest in the development of improved trajectory-based methods for the simulation of nonadiabatic processes in complex systems. An important goal for such methods is the accurate calculation of the rate constant over a wide range of electronic coupling strengths and it is often the nonadiabatic, weak-coupling limit, which being far from the Born-Oppenheimer regime, provides the greatest challenge to current methods. We show that in this limit there is an inherent sign problem impeding further development which originates from the use of the usual quantum flux correlation functions, which can be very oscillatory at short times. From linear response theory, we derive a modified flux correlation function for the calculation of nonadiabatic reaction rates, which still rigorously gives the correct result in the long-time limit regardless of electronic coupling strength, but unlike the usual formalism is not oscillatory in the weak-coupling regime. In particular, a trajectory simulation of the modified correlation function is naturally initialized in a region localized about the crossing of the potential energy surfaces. In the weak-coupling limit, a simple link can be found between the dynamics initialized from this transition-state region and an generalized quantum golden-rule transition-state theory, which is equivalent to Marcus theory in the classical harmonic limit. This new correlation function formalism thus provides a platform on which a wide variety of dynamical simulation methods can be built aiding the development of accurate nonadiabatic rate theories applicable to complex systems.

LES Modeling of the Impact of Heat Losses and Differential Diffusion on Turbulent Stratified Flame Propagation: Application to the TU Darmstadt Stratified Flame

Common combustion chambers often exhibit turbulent flames propagating in partially-premixed mixtures. This propagation is generally governed by aerodynamics, unsteady mixing and chemical processes and may also be affected by conductive heat losses when the reactive zone develops close to the burner lips. The Filtered TAbulated Chemistry for Large Eddy Simulation (F-TACLES) model has been recently developed to include tabulated chemistry in Large Eddy Simulation (LES) of adiabatic stratified flames in flamelet regimes. The present article proposes a modeling approach to account for both differential diffusion and non-adiabatic effects on flame consumption speed following the F-TACLES formalism. The adiabatic F-TACLES model is first detailed using a generalized formalism for diffusive fluxes allowing either to account for differential diffusion or not. The F-TACLES model is then extended to non-adiabatic situations. A correction factor based on the non-adiabatic consumption rate is introduced to recover a realistic filtered flame consumption speed. The objective is here to tackle flame stabilization mechanisms when heat losses affect the reaction zone. The proposed approach is validated through the simulation of the unconfined stratified turbulent jet flame TSF-A for which stabilization process is affected by heat losses. Five simulations are performed for both adiabatic and non-adiabatic flow conditions comparing unity Lewis number and complex diffusion assumptions. The adiabatic F-TACLES model predicts a flame anchored at the burner lip disagreeing with experimental data. The non-adiabatic simulation exhibits local extinction due to heat losses near the burner exit. The flame is then lifted improving the comparison with experiments. Results also show a significant impact of molecular diffusion model on both mean flame consumption rate and angle.

Quantum Stochastic Dynamics in the Presence of a Time-Periodic Rapidly Oscillating Potential: Nonadiabatic Escape Rate

Escape from a metastable state in the presence of a high-frequency field (where the driving becomes nonadiabatic) underlies a broad range of phenomena of physics and chemistry, and thus its understanding is of paramount importance. We study the problem of intermediate-to-high-damping escape from a metastable state of a dissipative system driven by a rapidly oscillating field, one of the most important classes of nonequilibrium systems, in a broad range of field driving frequencies (ω) and amplitudes (a). We construct a Langevin equation using quantum gauge transformation in the light of Floquet theorem and exploiting a systematic perturbative expansion in powers of 1/ω using "Kapitza-Landau time window". The quantum dynamics in a high-frequency field are found to be described by an effective time-independent potential. The temperature dependence of escape rate and the change of its form with varying parameters of the field have been analyzed. It may decrease upon increasing the temperature which is contingent on the effects of intricate interplay between external modulation and dissipation. The crossover temperature between tunnelling and thermal hopping increases with an increase in external modulation so that quantum effects in the escape are relevant at higher temperatures. These observations are uncommon and counterintuitive and, therefore, of considerable interest. Our results might be valuable for the exploration of the dynamics of cold atoms in electromagnetic fields.

Two-photon process via internal conversion by correlated photon pairs

We theoretically investigate a two-photon process via molecular internal conversion, driven by correlated photons with time delay. As an example, a molecular system approximated by Morse oscillators is used in order to evaluate the efficiency of the process, in terms of how the efficiency changes with the nonadiabatic coupling rate of vibronic states. We show that an efficient two-photon process requires two different correlated photon pairs with positive and negative energy correlations, depending on the value of the nonadiabatic coupling rate.

Theoretical Analysis of Proton Relays in Electrochemical Proton-Coupled Electron Transfer

The coupling of long-range electron transfer to proton transport over multiple sites plays a vital role in many biological and chemical processes. Recently the concerted proton-coupled electron transfer (PCET) reaction in a molecule with a hydrogen-bond relay inserted between the proton donor and acceptor sites was studied electrochemically. The standard rate constants and kinetic isotope effects (KIEs) were measured experimentally for this double proton transfer system and a related single proton transfer system. In the present paper, these systems are studied theoretically using vibronically nonadiabatic rate constant expressions for electrochemical PCET. Application of this approach to proton relays requires the calculation of multidimensional proton vibrational wave functions and the incorporation of multiple proton donor-acceptor motions. The decrease in proton donor-acceptor distances due to thermal fluctuations and the contributions from excited electron-proton vibronic states play important roles in these systems. The calculated KIEs and the ratio of the standard rate constants for the single and double proton transfer systems are in agreement with the experimental data. The calculations indicate that the standard PCET rate constant is lower for the double proton transfer system because of the smaller overlap integral between the ground state reduced and oxidized proton vibrational wave functions, resulting in greater contributions from excited electron-proton vibronic states with higher free energy barriers. The theory predicts that this rate constant may be increased by modifying the molecule in a manner that decreases the equilibrium proton donor-acceptor distances or alters the molecular thermal motions to facilitate the concurrent decrease of these distances. These insights may guide the design of more efficient catalysts for energy conversion devices.

Persistence time of charge carriers in defect states of molecular semiconductors

Charge carriers in organic crystals are often trapped in point defects. The persistence time of the charge in these defect states is evaluated by computing the escape rate from this state using non-adiabatic rate theory. Two cases are considered (i) the hopping between separate identical defect states and (ii) the hopping between a defect state and the bulk (delocalized) states. We show that only the second process is likely to happen with realistic defect concentrations and highlight that the inclusion of an effective quantum mode of vibration is essential for accurate computation of the rate. The computed persistence time as a function of the trap energy indicates that trap states shallower than ∼0.3 eV cannot be effectively investigated with some slow spectroscopic techniques such as THz spectroscopy or EPR commonly used to study the nature of excess charge in semiconductors.

A multi-dimensional microcanonical Monte Carlo study of S0 → T1 intersystem crossing of isocyanic acid

A general formula for multi-dimensional Monte Carlo microcanonical nonadiabatic rate constant expressed in configuration space is applied to calculate the rate of intersystem crossing (ISC) between ...

Theory of Proton-Coupled Electron Transfer in Energy Conversion Processes

Proton-coupled electron transfer (PCET) reactions play an essential role in a broad range of energy conversion processes, including photosynthesis and respiration. These reactions also form the basis of many types of solar fuel cells and electrochemical devices. Recent advances in the theory of PCET enable the prediction of the impact of system properties on the reaction rates. These predictions may guide the design of more efficient catalysts for energy production, including those based on artificial photosynthesis and solar energy conversion. This Account summarizes the theoretically predicted dependence of PCET rates on system properties and illustrates potential approaches for tuning the reaction rates in chemical systems. A general theoretical formulation for PCET reactions has been developed over the past decade. In this theory, PCET reactions are described in terms of nonadiabatic transitions between the reactant and product electron-proton vibronic states. A series of nonadiabatic rate constant expressions for both homogeneous and electrochemical PCET reactions have been derived in various well-defined limits. Recently this theory has been extended to include the effects of solvent dynamics and to describe ultrafast interfacial PCET. Analysis of the rate constant expressions provides insight into the underlying physical principles of PCET and enables the prediction of the dependence of the rates on the physical properties of the system. Moreover, the kinetic isotope effect, which is the ratio of the rates for hydrogen and deuterium, provides a useful mechanistic probe. Typically the PCET rate will increase as the electronic coupling and temperature increase and as the total reorganization energy and equilibrium proton donor-acceptor distance decrease. The rate constant is predicted to increase as the driving force becomes more negative, rather than exhibit turnover behavior in the inverted region, because excited vibronic product states associated with low free energy barriers and relatively large vibronic couplings become accessible. The physical basis for the experimentally observed pH dependence of PCET reactions has been debated in the literature. When the proton acceptor is a buffer species, the pH dependence may arise from the protonation equilibrium of the buffer. It could also arise from kinetic complexity of competing concerted and sequential PCET reaction pathways. In electrochemical PCET, the heterogeneous rate constants and current densities depend strongly on the overpotential. The change in equilibrium proton donor-acceptor distance upon electron transfer may lead to asymmetries in the Tafel plots and deviations of the transfer coefficient from the standard value of one-half at zero overpotential. Applications of this theory to experimentally studied systems illustrate approaches that can be utilized to tune the PCET rate. For example, the rate can be tuned by changing the pH or using different buffer species as proton acceptors. The rate can also be tuned with site-specific mutagenesis in biological systems or chemical modifications that vary the substituents on the redox species in chemical systems. Understanding the impact of these changes on the PCET rate may assist experimental efforts to enhance energy conversion processes.

Understanding the rate of spin-forbidden thermolysis of HN3 and CH3N3

The pyrolysis of the simplest azides HN(3) and CH(3)N(3) has been studied computationally. Nitrogen extrusion leads to the production of NH or CH(3)N. The azides have singlet ground states but the nitrenes CH(3)N and NH have triplet ground states. The competition between spin-allowed decomposition to the excited state singlet nitrenes and the spin-forbidden N(2) loss is explored using accurate electronic structure methods (CASSCF/cc-pVTZ and MR-AQCC/cc-pVTZ) as well as statistical rate theories. Nonadiabatic rate theories are used for the dissociation leading to the triplet nitrenes. For HN(3), (3)NH formation is predicted to dominate at low energy, and the calculated rate constant agrees very well with energy-resolved experimental measurements. Under thermal conditions, however, the singlet and triplet pathways are predicted to occur competitively, with the spin-allowed product increasingly favored at higher temperatures. For CH(3)N(3) thermolysis, spin-allowed dissociation to form (1)CH(3)N should largely dominate at all temperatures, with spin-forbidden formation of (3)CH(3)N almost negligible. Singlet methyl nitrene is very unstable and should rearrange to CH(2)NH immediately upon formation, and the latter species may lose H(2) competitively with vibrational cooling, depending on temperature and pressure.

Model system-bath Hamiltonian and nonadiabatic rate constants for proton-coupled electron transfer at electrode-solution interfaces

An extension of the Anderson-Newns-Schmickler model for electrochemical proton-coupled electron transfer (PCET) is presented. This model describes reactions in which electron transfer between a solute complex in solution and an electrode is coupled to proton transfer within the solute complex. The model Hamiltonian is derived in a basis of electron-proton vibronic states defined within a double adiabatic approximation for the electrons, transferring proton, and bath modes. The interaction term responsible for electronic transitions between the solute complex and the electrode depends on the proton donor-acceptor vibrational mode within the solute complex. This model Hamiltonian is used to derive the anodic and cathodic rate constants for nonadiabatic electrochemical PCET. The derivation is based on the master equations for the reduced density matrix of the electron-proton subsystem, which includes the electrons of the solute complex and the electrode, as well as the transferring proton. The rate constant expressions differ from analogous expressions for electrochemical electron transfer because of the summation over electron-proton vibronic states and the dependence of the couplings on the proton donor-acceptor vibrational motion. These differences lead to additional contributions to the total reorganization energy, an additional exponential temperature-dependent prefactor, and a temperature-dependent term in the effective activation energy that has different signs for the anodic and cathodic processes. This model can be generalized to describe both nonadiabatic and adiabatic electrochemical PCET reactions and provides the framework for the inclusion of additional effects, such as the breaking and forming of other chemical bonds.