Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Matrix”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

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

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2026

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

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

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

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

Photoinduced Charge Transfer versus Fragmentation Pathways in Lanthanum Cyclopentadienyl Complexes

This study compares two competing pathways of photoexcitations in gas-phase metal–organic complexes: first, a sequence of phonon-assisted electronic transitions leading to dissipation of the energy of photoexcitations and, second, a sequence of light-driven electronic transitions leading to photolysis. Phonon-assisted charge carrier dynamics is investigated by combination of the density matrix formalism and on-the-fly nonadiabatic couplings. Light-driven fragmentation is modeled by a time-dependent excited-state molecular-dynamics (TDESMD) algorithm based on Rabi theory and principles similar to the trajectory surface hopping approximation. Numerical results indicate that, under the medium intensity of the laser field, light-driven electronic transitions are more probable than phonon-assisted ones. The formation of multiple products is observed in TDESMD trajectories. Simulated mass spectra are extracted from TDESMD simulations and compared to experimental photoionization time-of-flight (PI-TOF) mass spectra. It is found that several features in the experimental mass spectra are reproduced by the simulations.

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

Can TDDFT Describe Excited Electronic States of Naphthol Photoacids? A Closer Look with EOM-CCSD

The 1 L b and 1 L a excited states of naphthols are characterized by using time-dependent density functional theory (TDDFT), configuration interaction with singles (CIS), and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) methods. TDDFT fails dramatically at predicting the energy and ordering of the 1 L a and 1 L b excited states as observed experimentally, while EOM-CCSD accurately predicts the excited states as characterized by natural transition orbital analysis. Here, the limitations of TDDFT are attributed to the absence of correlation from doubly excited configurations as well as the inconsistent description of excited electronic states of naphthol photoacids revealed by excitation analysis based on the one-electron transition density matrix.

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

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

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

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

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

Ultrafast Electronic Relaxation through a Conical Intersection: Nonadiabatic Dynamics Disentangled through an Oscillator Strength-Based Diabatization Framework

Here, we employ surface hopping trajectories to model the short-time dynamics of gas-phase and partially solvated 4-(N,N-dimethylamino)benzonitrile (DMABN), a dual fluorescent molecule that is known to undergo a nonadiabatic transition through a conical intersection. To compare theory vs time-resolved fluorescence measurements, we calculate the mixed quantum–classical density matrix and the ensemble averaged transition dipole moment. We introduce a diabatization scheme based on the oscillator strength to convert the TDDFT adiabatic states into diabatic states of L a and L b character. Somewhat surprisingly, we find that the rate of relaxation reported by emission to the ground state is almost 50% slower than the adiabatic population relaxation. Although our calculated adiabatic rates are largely consistent with previous theoretical calculations and no obvious effects of decoherence are seen, the diabatization procedure introduced here enables an explicit picture of dynamics in the branching plane, raising tantalizing questions about geometric phase effects in systems with dozens of atoms.

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

Excitonic Energy Landscape of the Y16F Mutant of the Chlorobium tepidum Fenna–Matthews–Olson (FMO) Complex: High Resolution Spectroscopic and Modeling Studies

We report high-resolution (low-temperature) absorption, emission, and nonresonant/resonant hole-burned (HB) spectra and results of excitonic calculations using a non-Markovian reduced density matrix theory (with an improved algorithm for parameter optimization in heterogeneous samples) obtained for the Y16F mutant of the Fenna–Matthews–Olson (FMO) trimer from the green sulfur bacterium Chlorobium tepidum. We show that the Y16F mutant is a mixture of FMO complexes with three independent low-energy traps (located near 817, 821, and 826 nm), in agreement with measured composite emission and HB spectra. Two of these traps belong to mutated FMO subpopulations characterized by significantly modified low-energy excitonic states. Hamiltonians for the two major subpopulations (Sub 821 and Sub 817 ) provide new insight into extensive changes induced by the single-point mutation in the vicinity of BChl 3 (where tyrosine Y16 was replaced with phenylalanine F16). The average decay time(s) from the higher exciton state(s) in the Y16F mutant depends on frequency and occurs on a picosecond time scale.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Mixed Upper Exciton State of the Special Pair in Bacterial Reaction Centers

Excitonic interactions between two closely separated bacteriochlorophyll a molecules (BChls) in the special pair of the reaction center (RC) of purple bacteria determine the positions and relative oscillator strengths of its two excitonic components. While the absorption of the lower excitonic band is well defined, the position and the intensity of the upper excitonic band (P Y+ ) are still under debate. Recent 77 K two-dimensional electronic spectroscopy data on Rba. capsulatus suggested that the P Y+ component absorbs at ~840 nm i.e. at a significantly lower energy than previously suggested. In the present work we argue that the P Y+ state is mixed with the excited states of the accessory BChls (B*/P Y+ ) leading to excitons contributing to the 785-825 nm spectral region which is consistent with previously published data. Our conclusion is based on hole-burning/linear dichroism data and modeling studies of the excitonic structure of the RC using non-Markovian reduced density matrix approach.

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

How the Size and Density of Charge-Transfer Excitons Depend on Heterojunction’s Architecture

We have characterized the size, intensity, density, and distribution of charge-transfer (CT) excitons as a function of the acceptor–donor architecture of prototypical organic interfaces. This characterization was done by computational analysis of 17 models of varying numbers, positions, and orientations of the donor and acceptor molecules. The models’ building blocks were phenyl-C61-butyric acid methyl ester (PCBM) fullerene acceptors and dual-band donor polymers composed of thiophene, benzothiadiazole, and benzotriazole subunits. The electronic structure of the donor–acceptor complexes was computed with the time-dependent long-range-corrected density-functional tight-binding method and analyzed with the fragment-based one-electron transition density matrix. In all models, the complexes with edge-on orientation have denser spectra of low-energy CT states lying below the absorption bands compared to the complexes with face-on orientation. This CT-state distribution in edge-on complexes provides a gate to efficiently populate cold CT excitons. Moreover, the cold CT excitons have a higher degree of charge separation in the edge-on than in the face-on complexes. The CT amount and the CT exciton size generally increase with the energy of the CT states, although the electron remains localized on a single molecule in cold CT states. Delocalization over two PCBM molecules was observed for high-energy CT states. The exciton size also depends on the orientation. Larger excitons are produced by the delocalization of the electrons perpendicularly to the interface. When the delocalization is parallel, the smaller electron–hole distances yield moderately sized CT excitons.

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

Spin Unrestricted Nonradiative Relaxation Dynamics of Cobalt-Doped Anatase Nanowire

Spin-resolved charge transfer dynamics at the interface of a Co(NH3)2-doped (001) anatase TiO2 nanowire and liquid water calculations based on density functional theory and density matrix formalism are considered. Three models with the same stoichiometry but different electronic structure are explored. While one model had no change to electron count and spin count (neutral model), the other two models were assigned a charge of 2+, one in a doublet and the other a quartet spin configuration. Co2+ is the most probable state for dopant in all models and Co acts as an electron acceptor. The optical absorption spectra show a rather unique pattern post 400 nm where the α and β absorptions happen independently at different frequency ranges. Essentially, the relaxation dynamics can be controlled as when an α electron is excited, there is a near zero probability of a β electron also undergoing an excitation and vice versa. The isolated exemption is between 400 and 650 nm in the neutral model. All the models absorb light in the visible range, while the electrons and holes have drastic differences in relaxation rates. The spatial charge separation occurs upon excitation and subsequent trapping, designating the considered system a prospect for electrochemical cell applications. (2+) models trap states more effectively slow down electron relaxation, making the charged models slightly better options for application compared to the neutral model.

Jensen, Stephanie J.↗
DOE OSTI2017

Hole Transfer in Dye-Sensitized Cesium Lead Halide Perovskite Photovoltaics: Effect of Interfacial Bonding

Lead halide perovskites have gained attention as an active material in solid-state dye-sensitized photovoltaics due to their high absorption of visible light and long charge-transport lengths. In perovskite-based dye-sensitized photovoltaic architectures the perovskite material is typically paired with a hole-transport material, such as spiro-OMeTAD, which extracts a hole from the photoexcited perovskite to generate free charge carriers. In this study, we explored two competing charge-transfer pathways at the interface between lead halide perovskite and spiro-OMeTAD: “through-bond” and “through-space”. For the through-bond pathway we use a segment of spiro-OMeTAD that contains methoxy linker groups, which will be referred to as “dye with methoxy linker groups” (DML). For the through-space pathway we use a segment of spiro-OMeTAD with removed linker groups, triphenylamine, which will be referred to as “dye”. Four atomistic models were studied: (I) a periodic cesium lead iodide (CsPbI3) perovskite nanowire (NW) that is paired with the dye molecule, (II) a periodic CsPbI3 perovskite NW paired with the DML molecule where the linker groups form coordination bond to the surface of the nanowire, (III) a CsPbI 3 perovskite thin film (TF) paired with the dye molecule, and (IV) a CsPbI 3 perovskite TF paired with the DML molecule. Charge-transfer dynamics, providing rates of electron/hole relaxation and relaxation pathways, are calculated using reduced density matrix formalism using Redfield theory. Finally, it was found that the terminal surface of the perovskite (Pb–I vs Cs–I) has important implications for energetic alignment at the perovskite–dye interface due to band bending. Finally, computed charge-transfer rates match well with upper and lower bounds of reported experimental results where “fast” picosecond rates correspond to through-bond pathway and “slow” nanosecond rates correspond to through-space.

14 SOLAR ENERGY↗
DOE OSTI2019

Photoexcited Electron Lifetimes Influenced by Momentum Dispersion in Silicon Nanowires

Silicon nanowires (SiNWs) exhibit unique optoelectronic properties originating from one-dimensional confinement effect. Technologically relevant properties—including band gap, radiative and nonradiative transitions, and charge carrier relaxation rates—are influenced by sampling of momentum along the growth direction. In this work, we quantify the influence of momentum dispersion on hot-electron relaxation rates and nonradiative lifetimes for SiNWs with ⟨100⟩ and ⟨111⟩ growth directions. Specifically, we account for nonradiative transitions between electronic states with different momentum, where changes in energy and momentum are achieved via interaction with a thermal bath of nuclear degrees of freedom. Photoexcited dynamic processes in reference SiNWs are computed via “on-the-fly” nonadiabatic couplings between electronic and nuclear degrees of freedom based on density functional theory (DFT). The dynamics of electronic degrees of freedom is propagated by a Redfield equation of motion for the reduced density matrix. Our results show that transitions allowing change of momentum prompt electron relaxation faster than those not allowing change of momentum. Our study also indicates that the electron relaxation time in ⟨100⟩ SiNW is longer than that in ⟨111⟩ SiNW.

Fatima, N/A↗
DOE OSTI2019

Absence of Coulomb Blockade in the Anderson Impurity Model at the Symmetric Point

In this study, we investigate the characteristics of the electric current in the so-called symmetric Anderson impurity model. We study the nonequilibrium model using two complementary approximate methods, the perturbative quantum master equation approach to the reduced density matrix and a self-consistent equation of motion approach to the nonequilibrium Green's function. We find that, at a particular symmetry point, an interacting Anderson impurity model recovers the same steady-state current as an equivalent noninteracting model, akin a two-band resonant level model. We show this in the Coulomb blockade regime for both high and low temperatures, where either the approximate master equation approach or the Green's function method provides accurate results for the current. We conclude that the steady-state current in the symmetric Anderson model at this regime does not encode characteristics of a many-body interacting system.

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

Entangled Electrons Drive a Non-superexchange Mechanism in a Cobalt Quinoid Dimer Complex

A central theme in chemistry is the understanding of the mechanisms that drive chemical transformations. A well-known, highly cited mechanism in organometallic chemistry is the superexchange mechanism in which unpaired electrons on two or more metal centers interact through an electron pair of the bridging ligand. We use a combination of novel synthesis and computation to show that such interactions may in fact occur by a more direct mechanism than superexchange that is based on direct quantum entanglement of the two metal centers. Specifically, we synthesize and experimentally characterize a novel cobalt dimer complex with benzoquinoid bridging ligands and investigate its electronic structure with the variational two-electron reduced density matrix method using large active spaces. The result draws novel connections between inorganic mechanisms and quantum entanglement, thereby opening new possibilities for the design of strongly correlated organometallic compounds whose magnetic and spin properties have applications in superconductors, energy storage, thermoelectrics, and spintronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗