Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonadiabatic coupling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2025

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Hard–Soft Acid–Base Theory Explains Photoexcited Carrier Dynamics in Porphyrin/CNT Nanohybrids: Time-Domain Atomistic Analysis

We employ the fundamental chemical concepts of hard− soft acid−base to formulate general principles governing excited-state dynamics in zinc porphyrin (ZnP)/carbon nanotube (CNT) hybrids for energy photoconversion. Atomistic quantum dynamics simulations demonstrate that electron-withdrawing and donating substituents at the ZnP β-pyrrolic position strongly influence the dynamics. ZnP photoexcitation produces subpicosecond electron transfer (ET) from ZnP to CNT, in agreement with the experiment. Substitutions of CN by H and tBu accelerate the ET. The trend is directly related to the hard−soft acid−base concept because the soft−soft interaction between the t Bu- ZnP acid and the mild CNT base enhances the donor−acceptor coupling. Longer coherence and more active vibrational modes facilitate the ET in t Bu-ZnP/CNT. Electron−hole recombination in CN-ZnP/CNT occurs on a hundred picosecond time scale, nicely corroborated by the experiment. The exciton lifetime is extended beyond a nanosecond by the substitutions. The soft−soft interaction in t Bu-ZnP/ CNT increases the splitting between the highest occupied orbitals of the two subsystems, reduces their mixing, and decreases nonadiabatic coupling between the ground and excited states. Rapid decoherence and involvement of low-frequency vibrations favor longer lifetimes. Our investigation reveals that larger p K a of the β-pyrrolic acid gives rapid ET and slow recombination and provides detailed mechanistic information, essential for future optoelectronic applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Photoselective isotope fractionation dynamics of N 2 with cosmo and atmospheric chemistry perspectives

Stable isotope ratio measurements provide valuable insights into a broad range of natural processes, from planetary atmospheres and climate to interstellar chemistry. Nitrogen, which has two stable isotopes, exhibits varying isotope ratios across the solar system. To model these observations, the isotope fraction as a function of energy is essential. At the Advanced Light Source (ALS), we measured the photodissociation of molecular nitrogen (N 2 ) with vacuum UV photons where a single photon is sufficiently energetic to dissociate the strong bond. The nitrogen atoms produced are scavenged with H 2 to form ammonia, whose isotopic makeup is determined. Blending the experiments with dynamical computations that include the shielding of light, we examine the isotopic composition and electronic atomic states produced. The measured photodissociation of N 2 at a natural isotopic composition with a frequency broad light beam exceptionally strongly favors the formation of the heavier nitrogen isotope, 15 N. Computations concur and suggest that the maximum in the quantum yield reflects significant variations in the specific electronic quantum states of the product N atoms that have quite different reactivities. Our quantum computations show that at similar energies, photodissociation of 14 N 14 N and 15 N 14 N can lead to different product channels. The computed dynamics include extensive state-selective spin–orbit and nonadiabatic couplings affecting the light absorption and dissociation pathways that proceed via the triplet manifold of states. Our results are relevant for future exploration missions, both in situ and sample-return and for other molecules such as O 2 and CO.

Komarova, Ksenia [Hebrew Univ. of Jerusalem (Israe↗
DOE OSTI2025

Photochemical Au(I)–Au(I) Bond Formation: A Battle between Intersystem Crossing and Internal Conversion

The transition metal complex Au(CN) 2 – has provided experimental evidence of photoinduced bond formation between Au(I) atoms in solution. However, the underlying photochemical driving force for this bond formation reaction remains unclear. In this study, we investigate the ultrafast Au–Au bonding process in the [Au(CN) 2 ] 2 2– dimer using nonadiabatic dynamics simulations that incorporate intersystem crossing and internal conversion pathways. Reaction pathways and transitions among photochemically accessible singlet and triplet excited states are analyzed. Computational results indicate that intersystem crossing is the primary driving force in the early stages of ultrafast photochemical dynamics, while internal conversion among triplet states plays a critical role after the system stabilizes in a higher-lying triplet state. Furthermore, this work provides a mechanistic perspective on modulating photochemical reactions by tuning the relative strengths of spin–orbit coupling and nonadiabatic coupling.

Computational chemistry↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2025

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Roadmap for Molecular Benchmarks in Nonadiabatic Dynamics

Simulating the coupled electronic and nuclear response of a molecule to light excitation requires the application of nonadiabatic molecular dynamics. However, when faced with a specific photophysical or photochemical problem, selecting the most suitable theoretical approach from the wide array of available techniques is not a trivial task. The challenge is further complicated by the lack of systematic method comparisons and rigorous testing on realistic molecular systems. This absence of comprehensive molecular benchmarks remains a major obstacle to advances within the field of nonadiabatic molecular dynamics. A CECAM workshop, Standardizing Nonadiabatic Dynamics: Towards Common Benchmarks, was held in May 2024 to address this issue. This Perspective highlights the key challenges identified during the workshop in defining molecular benchmarks for nonadiabatic dynamics. Specifically, this work outlines some preliminary observations on essential components needed for simulations and proposes a roadmap aiming to establish, as an ultimate goal, a community-driven, standardized molecular benchmark set.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Jahn–Teller distortion controls electron transfer in photoexcited Cu( I ) donor–acceptor systems

The Jahn–Teller distortion (JTD) is a defining structural response to electronic excitation in Cu( I )-based transition metal complexes, yet its role in photoinduced electron transfer (PET) remains largely unexplored. Here, we demonstrate that the JTD governs charge separation (CS) through a vibronically controlled conical intersection in heteroleptic Cu( I ) bisphenanthroline–naphthalene diimide (CuHETPHEN-NDI) donor–acceptor dyads. Using 20-fs broadband transient absorption spectroscopy combined with coherent vibrational wavepacket (CVWP) analysis and quantum chemical calculations, we directly track nuclear motions that steer the system from the metal-to-ligand charge-transfer ( 1 MLCT) state to the CS state. Steric bulkiness of the pendant groups at the 2,9-positions of the phenanthroline ligand systematically slows both JTD and CS. Short-time Fourier transformation and Fourier filtering analyses identify two key vibrational signatures: a low-frequency breathing mode (∼100 cm −1 ) via a bond distance change between Cu and ligated Ns (Cu–N) that modulates the NDI anion absorption and acts as a vibronic coupling coordinate, and a higher-frequency mode (∼313 cm −1 ) that evolves along the PET trajectory. Normal mode analysis and potential energy surface calculations show that the JTD brings the 1 MLCT and CS states into degeneracy, while the Cu–N breathing motion dynamically modulates donor–acceptor electronic coupling to enable ultrafast nonadiabatic electron transfer. Steric hindrance exerted by the groups at the 2,9 positions of the phenanthroline ligands suppresses this vibronic coupling, leading to faster CVWP decoherence for the 313 cm −1 mode and slower CS. These findings unravel JTD-controlled vibronic coupling at conical intersection as a governing factor for CS and provide insight into designing Cu-based photosensitizers by harnessing structural dynamics to control PET.

Kim, Pyosang [Argonne National Laboratory (ANL), A↗
DOE OSTI2026

Transition-state lattice modes and the breakdown of adiabatic tunneling for hydrogen and deuterium in bcc Nb

Interstitial hydrogen and deuterium in body-centered-cubic metals constitute archetypal quantum tunneling systems. Their relevance has been renewed by the connection between hydrogenic tunneling in Nb and defect-induced decoherence in superconducting qubits, motivating a predictive microscopic theory. Existing theoretical treatments invoke an adiabatic separation between the light interstitial and the host lattice, an assumption whose validity has not been rigorously established for hydrogenic species. Here, we show that the experimentally measured tunnel splittings of O-trapped H and D in bcc Nb are quantitatively reproduced only within a five-dimensional (5D) Lattice-Renormalized Born-Oppenheimer (LRBO) framework. This approach treats three interstitial modes and two judiciously selected lattice modes, which includes a transition-state mode, on equal quantum footing. By recasting nested Born-Oppenheimer hierarchies within this same formalism and benchmarking against modern \textit{ab initio} potential energy surfaces, we show that adiabatic separation of the light particle from lattice dynamics is satisfied only in the positive-muon ($μ^{+}$) mass limit. In contrast, tunneling for H and D is fundamentally a collective, nonadiabatic process mediated by anharmonic lattice couplings. Finally, we show that the breakdown of adiabaticity can be anticipated from simple energy estimates involving the ground-state light-particle energy evaluated at a small number of fixed lattice configurations, providing a practical criterion for assessing the validity of adiabatic tunneling theories in other systems.

Pritchard, P. Graham [Northwestern U.]↗
DOE OSTI2026

Many roads to the seam: How conformational flexibility drives nonadiabatic relaxation in a prototypical tetrapyrrolic chromophore

Large and structurally flexible chromophores pose challenges for in silico modeling of photodeactivation due to the many vibrational modes that can funnel the system toward energy degeneracy. In this work, we examine how the multiple degrees of freedom in biliverdin, a prototypical tetrapyrrolic chromophore, cooperate to drive access to the S 1 /S 0 intersection seam in vacuo. We begin by mapping the ground-state potential energy surface to identify representative biliverdin conformers relevant to photoexcitation. We then use a CASSCF-based framework to map the excited-state landscape and characterize the intersection seam, identifying distinct conical-intersection types. Finally, we employ ab initio multiple spawning to resolve the dynamical pathways by which the system accesses these regions. DFT potential-energy and free-energy mappings indicate that, although several conformers are relevant, the “locked-helix” ZsZsZs conformer predominates in the ground state. The intersection seam comprises numerous geometrically distinct regions characterized by varying degrees and combinations of dihedral torsion, pyramidalization, and bond-length alternation. Yet only select regions lie within energetic reach, and moderate barriers separate them from the S 1 minimum. Nonadiabatic dynamics combined with multivariate analyses show that, despite extensive mode coupling during deactivation that guides the system toward multiple regions of the seam, a single dihedral torsion, together with bond-length alternation, predominantly drives energy degeneracy. This work offers new insight into biliverdin’s intrinsic photochemical response and underscores a general feature of flexible chromophores: many modes may participate during photorelaxation, but only a limited subset ultimately dictates seam accessibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Nonadiabatically Driven Quantum Interference Effects in the Ultracold K + KRb → Rb + K 2 Chemical Reaction

The K + KRb → Rb + K 2 chemical reaction is the first ultracold atom–diatom chemical reaction for which experimental results have been reported for temperatures below 1 μK more than a decade ago. The reaction occurs through coupling with an excited electronic state that is accessible even in the ultracold limit. A previous quantum dynamics study, excluding nonadiabatic effects, has reported a rate coefficient that is about 35% below the experimental value. Here, we report the first nonadiabatic quantum dynamics study of this reaction and obtain rate coefficients in better agreement with experiments. Our results show that short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Signatures of Antisymmetric Vibrations in the Ultrafast Dynamics of Quadrupolar Dyes

Antisymmetric molecular vibrations are central to ultrafast, nonadiabatic photophysical and photochemical processes such as conical intersection dynamics, Herzberg–Teller couplings and, potentially, singlet fission. Direct spectroscopic identification of such vibrations is, however, challenging, since they are typically Raman inactive and affect optical transitions only weakly. Here, we report experimental signatures of vibronic couplings to a high-frequency antisymmetric vibration in the excited state dynamics of a quasi-quadrupolar molecule by ultrafast two-dimensional electronic spectroscopy (2DES). The early time, sub-50 fs 2DES maps reveal an asymmetric peak pattern with characteristic low-energy cross-peaks. We show that these peaks arise from stimulated emission transitions from an anharmonic, double-minimum excited state potential energy surface formed by vibronic coupling to a high-frequency antisymmetric mode. Simulations based on a phenomenological essential state model support the results. Our findings offer a new approach for identifying antisymmetric vibrations in ultrafast 2DES and track excited state wavepacket motion before relaxation washes out the spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fast methods for multisite charge transfer processes. I. Constrained, state averaged CASSCF(1,n) and CASSCF(2n − 1,n) simulations

We design a dynamically weighted state-averaged constrained complete active space self-consistent field (DW-SA-cCASSCF) algorithm to treat electrons or holes moving between n molecular fragments (where n can be larger than 2). Within such a so-called eDSCn/hDSCn approach, we consider configurations that are mutually single excitations of each other, and we apply a generalized set of constraints to tailor the method for studying charge transfer problems. The constrained optimization problem is efficiently solved using a DIIS-SQP algorithm, thus maintaining computational efficiency. We demonstrate the method for a finite Su–Schrieffer–Heeger chain, successfully reproducing the expected exponential decay of diabatic couplings with distance. When combined with a gradient, the current extension immediately enables efficient nonadiabatic dynamics simulations of complex multi-state charge transfer processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Mechanism of the Non‐Kasha Fluorescence in Pyrene

The high-energy shoulder in the gas-phase fluorescence emission spectrum of pyrene is a well-known example of non-Kasha emission. Here, we comparatively assess two approaches, vibronic perturbation theory and nonadiabatic dynamics, in their ability to predict and explain the gas-phase fluorescence spectrum of pyrene. While both methods qualitatively capture the non-Kasha emission, they differ in their computational requirements, accuracy, and physical interpretation. Vibronic perturbation theory and nonadiabatic dynamics are complementary and can be combined in a two-step approach to non-Kasha fluorescence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Carrier Localization and Spontaneous Formation of Two-Dimensional Polarization Domain in Halide Perovskites

Halide perovskites are known for their rich phase diagram and superior performance in diverse optoelectronics applications. The latter property is often attributed to the long electron-hole recombination time, whose underlying physical mechanism has been a long-standing controversy. In this Letter, we investigate the transport and localization properties of electron and hole carriers in a prototypical halide perovskite (CsPbBr 3 ), through ab initio tight-binding nonadiabatic dynamics approach for large-scale (tens of nm size) supercell calculations, to simulate electron and ion dynamics on the same footing. Here, we found distinct structural, lattice polarization, and electron-phonon coupling properties at low (below 100 K) and high temperatures, consistent with experimental observations. In particular, at low temperature we find spontaneous formation of polar grain boundaries in the nonpolar bulk systems, which result in two-dimensional polarization patterns that serve to localize and separate electrons and holes. We reveal phonon-assisted variable-range hopping mostly responsible for low-temperature transport, and their characteristic frequency correlates with temperature-dependent phonon power spectrum and energy oscillation frequency in nonadiabatic dynamics. We answer the critical questions of long electron-hole recombination lifetime at low temperature and offer the correlation among polarization domains, electron-phonon couplings, and photocarrier dynamics.

Carrier generation & recombination↗