Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-adiabaticity”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ultrafast formation of topological defects in a two-dimensional charge density wave

Topological defects play a central role in dynamical systems undergoing a non-adiabatic transition. In solids, topological defects as a result of femtosecond laser excitation have attracted increasing interest not only because they are key to understanding phase transitions but also because they can generate a variety of hidden orders that are not accessible in thermal equilibrium. Despite the common occurrence of these defects in a non-equilibrium system, the fundamental limit on how fast they can emerge in solids and the generic pathway for defect creation at such fast timescales have remained open questions. Here we apply ultrafast electron diffraction to study the reciprocal-space signatures of transient defects in a two-dimensional charge density wave, where simultaneous measurements of both defect and phonon dynamics yield a microscopic view of defect formation in the femtosecond regime. We find that one-dimensional domain walls are generated well within 1 ps following photoexcitation, during which the defect growth is not dictated by the amplitude of the order parameter but is mediated by a non-thermal population of longitudinal optical phonons. In conclusion, our work provides a framework for the ultrafast engineering of topological defects that are coupled to specific collective modes, which will prove useful for the dynamical control of non-equilibrium phases in correlated materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Excited state dynamics of azanaphthalenes reveal opportunities for the rational design of photoactive molecules

Abstract Various photoactive molecules contain motifs built on aza-aromatic heterocycles, although a detailed understanding of the excited state photophysics and photochemistry in such systems is not fully developed. To help address this issue, the non-adiabatic dynamics operating in azanaphthalenes under hexane solvation was studied following 267 nm excitation using ultrafast transient absorption spectroscopy. Specifically, the species quinoline, isoquinoline, quinazoline, quinoxaline, 1,6-naphthyridine, and 1,8-naphthyridine were investigated, providing a systematic variation in the relative positioning of nitrogen heteroatom centres within a bicyclic aromatic structure. Our results indicate considerable differences in excited state lifetimes, and in the propensity for intersystem crossingvsinternal conversion across the molecular series. The overall pattern of behaviour can be explained in terms of potential energy barriers and spin-orbit coupling effects, as demonstrated by extensive quantum chemistry calculations undertaken at the SCS-ADC(2) level of theory. The fact that quantum chemistry calculations can achieve such detailed and nuanced agreement with experimental data across a full set of six molecules exhibiting subtle variations in their composition provides an excellent example of the current state-of-the-art and is indicative of future opportunities for rational design of photoactive molecules.

Chemistry↗
DOE OSTI2025

High-throughput and data-driven search for stable optoelectronic AMSe 3 materials

The rapid advancement in emerging optoelectronic technologies demands highly efficient, affordable, and ecofriendly materials. In this context, ternary chalcogenides, especially ternary selenides, show early promise as a material class due to their stability and remarkable electronic, optical, and transport properties. In this work, we integrate first-principles-based high-throughput computations with machine learning (ML) techniques to predict the thermodynamic stability and optoelectronic properties of 920 valency-satisfied selenide compounds. Through investigating polymorphism, our study reveals the edge-sharing orthorhombic Pnma phase (NH 4 CdCl 3 -type) as the most stable structure for most ternary selenides. High-fidelity supervised ML models are trained and tested to accelerate stability and band gap predictions. These data-driven models pin down the most influential features that dominantly control key material characteristics. The multistep high-throughput computations identify the ternary selenides with optimal direct band gaps, light carrier masses, and strong optical absorption edges. The extensive materials screening considering phase stability, toxicity, and defect tolerance, finally identifies the seven most suitable candidates for photovoltaic applications. Two of these final compounds, SrZrSe 3 and SrHfSe 3 , have already been synthesized in a single-phase form, with the latter showing an optically suitable band gap, aligning well with our findings. The non-adiabatic molecular dynamics reveal sufficiently long photoexcited charge carrier lifetimes (on the order of nanoseconds) in some of these selected selenide materials, indicating their exciting characteristics. Overall, our study suggests a robust in silico framework that can be extended to screen large datasets of various material classes for identifying promising photoactive candidates.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ultrafast photochemistry and electron-diffraction spectra in n → (3s) Rydberg excited cyclobutanone resolved at the multireference perturbative level

We study the ultrafast time evolution of cyclobutanone excited to the singlet n → Rydberg state through non-adiabatic surface-hopping simulationsperformed at extended multi-state complete active space second-order perturbation (XMS-CASPT2) level of theory. These dynamics predict relaxation to the ground-state with a timescale of 822 ± 45 fs with minimal involvement of the triplets. The major relaxation path to the ground-state involves a three-state degeneracy region and leads to a variety of fragmented photoproducts. We simulate the resulting time-resolved electron-diffraction spectra, which track the relaxation of the excited state and the formation of various photoproducts in the ground state.

Chemistry↗
DOE OSTI2024

Ultrafast electron diffraction of photoexcited gas-phase cyclobutanone predicted by ab initio multiple cloning simulations

We present the result of our calculations of ultrafast electron diffraction (UED) for cyclobutanone excited into the S2 electronic state, which is based on the non-adiabatic dynamics simulations with the Ab Initio Multiple Cloning (AIMC) method with the electronic structure calculated at the SA(3)-CASSCF(12,12)/aug-cc-pVDZ level of theory. The key features in the UED pattern were identified, which can be used to distinguish between the reaction pathways observed in the AIMC dynamics, although there is a significant overlap between representative signals due to the structural similarity of the products. The calculated UED pattern can be compared with the experiment.

Chemistry↗
DOE OSTI2024

Polariton-assisted incoherent to coherent excitation energy transfer between colloidal nanocrystal quantum dots

Here, we explore the dynamics of energy transfer between two nanocrystal quantum dots placed within an optical microcavity. By adjusting the coupling strength between the cavity photon mode and the quantum dots, we have the capacity to fine-tune the effective coupling between the donor and acceptor. Introducing a non-adiabatic parameter, γ, governed by the coupling to the cavity mode, we demonstrate the system’s capability to shift from the overdamped Förster regime (γ $\ll$ 1) to an underdamped coherent regime (γ $\gg$ 1). In the latter regime, characterized by swift energy transfer rates, the dynamics are influenced by decoherence time. To illustrate this, we study the exciton energy transfer dynamics between two closely positioned CdSe/CdS core/shell quantum dots with sizes and separations relevant to experimental conditions. Employing an atomistic approach, we calculate the excitonic level arrangement, exciton–phonon interactions, and transition dipole moments of the quantum dots within the microcavity. These parameters are then utilized to define a model Hamiltonian. Subsequently, we apply a generalized non-Markovian quantum Redfield equation to delineate the dynamics within the polaritonic framework.

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

Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2025

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗
DOE OSTI2026

General framework for quantifying dissipation pathways in open quantum systems. III. Off-diagonal subsystem–bath couplings

This paper extends the previously reported theory of dissipation pathways [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111 (2024)] to incorporate off-diagonal subsystem–bath coupling, which is often required to model molecular systems where the environment directly influences transitions and couplings between subsystem states. We systematically derive master equations for both population transfer and dissipation into individual bath components, for which we also rigorously prove energy conservation and detailed balance. The approach is based on second-order perturbation theory with respect to the subsystem–bath couplings, whose form is not limited to any specific model. The accuracy of the developed method is tested by applying it to diverse model Hamiltonians involving linearly coupled harmonic oscillator baths and comparing the outcomes against the hierarchical equations of motion (HEOM) method. Overall, our method accurately quantifies the contributions of specific bath components to the overall dissipation while significantly reducing the computational cost compared to numerically exact methods such as HEOM, thus offering a path to examine how vibronic interactions steer non-adiabatic processes in realistic chemical systems.

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

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Measuring the molecular Migdal effect with neutron scattering on diatomic gases

The Migdal effect is a key inelastic signal channel which could be used to detect low-mass dark matter, but it has never been observed experimentally using Standard Model probes. Here we propose a conceptual design for an experiment which could detect the Migdal effect in diatomic molecules through low-energy neutron scattering, and we provide the requirements on the beam spectrum and reducible backgrounds such that a detection may be achieved. The enhancement of the Migdal rate through non-adiabatic couplings, which are absent in isolated atoms, combined with the distinctive photon energies of electronic transitions in CO, suggest that a positive detection of the molecular Migdal effect may be possible with modest beam times at existing neutron facilities. Published by the American Physical Society 2024

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Spin-Controllable Dynamics in Defect-Engineered Carbon Nanotubes as Single Photon Emitters: Data-Driven Modeling and Computations

Quantum technologies, such as quantum computing and sensing, require efficient single-photon emission (SPE) sources that operate at room temperature in telecom wavelengths. While several materials can serve as SPE sources, no single platform meets all the criteria for efficiency, ambient operation, and scalability. Single-walled carbon nanotubes (SWCNTs) with covalently attached molecules offer a promising solution. Their SPE can be easily tuned via modifications of the SWCNT's diameter, chirality, and bonded molecules, enabling emission across near-IR to telecom wavelengths at ambient conditions. However, to fully realize the potential of SWCNTs and unlock their quantum capabilities, a deeper understanding of how structural defects from molecular adducts affect their emission and competing photoexcited processes is essential. To address this gap in our knowledge, this project combined quantum chemistry calculations with data-driven methods of cheminformatics (QSAR) and machine learning (ML). The developed computational approaches have provided several design strategies for covalent functionalization of SWCNTs to improve their optical response. The collaboration with Los Alamos National Lab (LANL) enabled direct comparison of computational and experimental data, facilitating method validation. This partnership was enhanced through access to LANL's Center for Integrated Nanotechnologies (CINT) utilizing User Facility Program and summer internships, which provided three NDSU graduate students with hands-on experience at LANL. The outcomes of this project included (1) Advancing the current stage of computational methods in accurate modeling of non-adiabatic spin-dependent photoexcited dynamics and its applicability to nanosystems consisting of thousands of atoms, realized as open-access codes linked to existing DFT-based software; (2) Establishing the relationship between the structure of adducts and SWCNTs and intrinsic excitonic and spin properties of defect states for guiding novel synthetic strategies and experimental probes of chemically functionalized SWCNTs as near-IR emitting materials; (3) Generating virtual libraries of hypothetical functionalized SWCNTs for virtual screening of their chemical structures and optical properties, leveraging new functionalities of SWCNTs; (4) Offering a unique experience for NDSU graduate students that prepared them for future scientific careers related to materials modeling and big data processing. These results were summarized in 12 published journal papers and 3 recently submitted papers. One of a key finding is that the position of defect sites on the SWCNT surface primarily drives the emission redshift (up to 100 meV), while the polarity of the defect-inducing molecules has a much smaller effect (~10 meV). However, the electron-donating or withdrawing properties of a molecule influence selecting reactivity of defect sites. These insights important for optimizing synthetic protocols for desired emissions in SWCNTs. We also revealed that the interaction between two defects at various positions on the SWCNT enhances the redshift and optical activity of states, favoring strong near-IR emission. This suggests that manipulations in defect concentrations is a promising strategy for controlling efficient emission. Mostly important, the defect position was found controllable by the spin states of photoexcited intermediates: Excited aromatic molecules form ortho defects with SWCNTs at their singlet states in the presence of oxygen, while oxygen-free conditions favor para defects via the triplet-state mechanism. Additionally, a heat-activated [2+2] cycloaddition reaction facilitates divalent defect formation with fewer bonding positions that narrows emission bands. These groundbreaking findings have been experimentally validated and significantly advance our understanding of defect chemistry in SWCNTs. Using a novel encoding technique and 3D-MoRSE descriptors, we developed highly accurate ML/QSAR models to predict both the 3D structure and optical properties of SWCNTs with chemical defects. This model enabled the creation of a virtual library of 125,556 structures, providing new insights into the relationship between SWCNT-defect structure and emission.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗
DOE OSTI2024

Ligand Controls Excited Charge Carrier Dynamics in Metal-Rich CdSe Quantum Dots: Computational Insights

Small metal-rich semiconducting quantum dots (QDs) are promising for solid-state lighting and single-photon emission due to their highly tunable yet narrow emission line widths. Nonetheless, the anionic ligands commonly employed to passivate these QDs exert a substantial influence on the optoelectronic characteristics, primarily owing to strong electron–phonon interactions. In this work, we combine time-domain density functional theory and nonadiabatic molecular dynamics to investigate the excited charge carrier dynamics of Cd 28 Se 17 X 22 QDs (X = HCOO – , OH – , Cl – , and SH – ) at ambient conditions. These chemically distinct but regularly used molecular groups influence the dynamic surface-ligand interfacial interactions in Cd-rich QDs, drastically modifying their vibrational characteristics. The strong electron–phonon coupling leads to substantial transient variations at the band edge states. The strength of these interactions closely depends on the physicochemical characteristics of passivating ligands. Consequently, the ligands largely control the nonradiative recombination rates and emission characteristics in these QDs. Our simulations indicate that Cd 28 Se 17 (OH) 22 has the fastest nonradiative recombination rate due to the strongest electron–phonon interactions. Conversely, QDs passivated with thiolate or chloride exhibit considerably longer carrier lifetimes and suppressed nonradiative processes. The ligand-controlled electron–phonon interactions further give rise to the broadest and narrowest intrinsic optical line widths for OH and Cl-passivated single QDs, respectively. Finally, obtained computational insights lay the groundwork for designing appropriate passivating ligands on metal-rich QDs, making them suitable for a wide range of applications, from blue LEDs to quantum emitters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

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

Polariton spectra under the collective coupling regime. II. 2D non-linear spectra

In our previous work [Mondal et al., J. Chem. Phys. 162, 014114 (2025)], we developed several efficient computational approaches to simulate exciton–polariton dynamics described by the Holstein–Tavis–Cummings (HTC) Hamiltonian under the collective coupling regime. Here, we incorporated these strategies into the previously developed Lindblad-partially linearized density matrix (⁠$\mathscr{L}$-PLDM) approach for simulating 2D electronic spectroscopy (2DES) of exciton–polariton under the collective coupling regime. In particular, we apply the efficient quantum dynamics propagation scheme developed in Paper I to both the forward and the backward propagations in the PLDM and develop an efficient importance sampling scheme and graphics processing unit vectorization scheme that allow us to reduce the computational costs from $\mathscr{O}$($\mathscr{K}$ 2 )$\mathscr{O}$(T 3 ) to $\mathscr{O}$($\mathscr{K}$)$\mathscr{O}$(T 0 ) for the 2DES simulation, where $\mathscr{K}$ is the number of states and T is the number of time steps of propagation. As a result, we further simulated the 2DES for an HTC Hamiltonian under the collective coupling regime and analyzed the signal from both rephasing and non-rephasing contributions of the ground state bleaching, excited state emission, and stimulated emission pathways.

2D non-linear spectra↗