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

Influence of globally spin-aligned vector mesons to the measurements of the chiral magnetic effect in heavy-ion collisions

The chiral magnetic effect (CME) in high-energy heavy-ion collisions arises from the interplay between the chirality imbalance and the intense magnetic field and will cause a charge separation along the magnetic field direction. While the CME search is still ongoing in experiments, the non-CME contributions need to be excluded from the CME observables. In this work, we examine the influence of globally spin-aligned ρ mesons on the γ 112 correlator, the R Ψ 2 (ΔS) correlator, and the signed balance functions, via a toy model and a multiphase transport model (AMPT). We find that the CME observables are sensitive to the 00-component of the spin density matrix, ρ 00 : they receive positive (negative) contributions when ρ 00 is larger (smaller) than 1/3.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Time Domain Simulations of Chemical Bonding Effects in Surface-Enhanced Spectroscopy

We employ the atom-centered density-matrix propagation method to illustrate how time-dependent structural and conformational changes affect the electronic structure and derived spectroscopic properties of a prototypical finite metal cluster-bound π-conjugated organic, namely the Ag7-benzenethiol complex. We find that the calculated spectroscopic properties are dictated by large amplitude motion which controls the coupling between the aromatic ring of the molecule and the metal cluster. The simulated vibrational spectra of Ag7-benzenethiol are in accord with previous experiments which probe Raman scattering from benzenethiol adsorbed on silver surfaces.

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

Computational modeling of latent-heat-storage in PCM modified interior plaster

The latent heat storage systems represent a promising way for decrease of buildings energy consumption with respect to the sustainable development principles of building industry. The presented paper is focused on the evaluation of the effect of PCM incorporation on thermal performance of cement-lime plasters. For basic characterization of the developed materials, matrix density, bulk density, and total open porosity are measured. Thermal conductivity is accessed by transient impulse method. DSC analysis is used for the identification of phase change temperature during the heating and cooling process. Using DSC data, the temperature dependent specific heat capacity is calculated. On the basis of the experiments performed, the supposed improvement of the energy efficiency of characteristic building envelope system where the designed plasters are likely to be used is evaluated by a computational analysis. Obtained experimental and computational results show a potential of PCM modified plasters for improvement of thermal stability of buildings and moderation of interior climate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Investigation of the ozone formation reaction pathway: Comparisons of full configuration interaction quantum Monte Carlo and fixed-node diffusion Monte Carlo with contracted and uncontracted MRCI

The association/dissociation reaction path for ozone (O 2 + O ↔ O 3 ) is notoriously difficult to describe accurately using ab initio electronic structure theory, due to the importance of both strong and dynamic electron correlations. Experimentally, spectroscopic studies of the highest lying recorded vibrational states combined with the observed negative temperature dependence of the kinetics of oxygen isotope exchange reactions confirm that the reaction is barrierless, consistent with the latest potential energy surfaces. Previously reported potentials based on Davidson-corrected internally contracted multireference configuration interaction (MRCI) suffer from a spurious reef feature in the entrance channel even when extrapolated towards the complete basis set limit. In this work, we report an analysis of comparisons between a variety of electronic structure methods including internally contracted and uncontracted MRCI (with and without Davidson corrections), as well as full configuration interaction quantum Monte Carlo, fixed-node diffusion Monte Carlo, and density matrix renormalization group.

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

Non-orthogonal multi-Slater determinant expansions in auxiliary field quantum Monte Carlo

We investigate here the use of non-orthogonal multi-Slater determinant (NOMSD) expansions as trial wavefunctions in auxiliary field quantum Monte Carlo simulations of molecular systems. We show that NOMSD trial wavefunctions with as few as twenty determinants are sufficient in order to achieve chemical accuracy across most of the G1 molecular test set. We also show that NOMSD trial wavefunctions are useful for more challenging strongly correlated systems by computing relative energies along the isomerization path of the [Cu 2 O 2 ] 2+ molecule. Our results for [Cu 2 O 2 ] 2+ compare favorably with other accurate quantum many-body methods, including density matrix renormalization group and completely renormalized coupled cluster methods.

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

Numerical tests of coherence-corrected surface hopping methods using a donor-bridge-acceptor model system

Surface hopping (SH) is a popular mixed quantum-classical method for modeling nonadiabatic excited state processes in molecules and condensed phase materials. The method is simple, efficient, and easy to implement, but the use of classical and independent nuclear trajectories introduces an overcoherence in the electronic density matrix which, if ignored, often leads to spurious results, such as overestimated reaction rates. Several methods have been proposed to incorporate decoherence into SH simulations, but a lack of insightful benchmarks makes their relative accuracy unknown. Herein, we run numerical simulations of common coherence-corrected SH methods including Truhlar’s decay-of-mixing (DOM) and Subotnik’s augmented SH using a Donor-bridge-Acceptor (DbA) model system. Numerical simulations are carried out in the superexchange regime, where charge transfer proceeds from a donor to an acceptor as a result of donor-bridge and bridge-acceptor couplings. The computed donor-to-acceptor reaction rates are compared to the reference Marcus theory results. For the DbA model under consideration, augmented SH recovers Marcus theory with quantitative accuracy, whereas DOM is only qualitatively accurate depending on whether predefined parameters in the decoherence rate are chosen wisely. We propose a general method for parameterizing the decoherence rate in the DOM method, which improves the method’s reaction rates and presumably increases its transferability. Altogether, the decoherence method of choice must be chosen with great care and this work provides insight using an exactly solvable model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toward DMRG-tailored coupled cluster method in the 4c-relativistic domain

There are three essential problems in computational relativistic chemistry: Electrons moving at relativistic speeds, close lying states, and dynamical correlation. Currently available quantum-chemical methods are capable of solving systems with one or two of these issues. However, there is a significant class of molecules in which all the three effects are present. These are the heavier transition metal compounds, lanthanides, and actinides with open d or f shells. For such systems, sufficiently accurate numerical methods are not available, which hinders the application of theoretical chemistry in this field. In this paper, we combine two numerical methods in order to address this challenging class of molecules. These are the relativistic versions of coupled cluster methods and the density matrix renormalization group (DMRG) method. To the best of our knowledge, this is the first relativistic implementation of the coupled cluster method externally corrected by DMRG. The method brings a significant reduction of computational costs as we demonstrate on the system of TlH, AsH, and SbH.

Brandejs, Jan (ORCID:0000000221073095)↗
DOE OSTI2021

Two-color, intracavity pump–probe, cavity ringdown spectroscopy

Here, we report a proof-of-principle demonstration of intracavity pump–probe, cavity ringdown (CRD) detection in a three-mirror, traveling-wave cavity. With cavity-enhanced pump power and probe absorption path length, the technique is a generally applicable, high-sensitivity, high-selectivity detection method. In our experiments, the pump radiation is switched off during every other probe ringdown, which allows uncorrelated measurements of analyte and background cavity decay rates. The net, two-color signal from the difference between the pump-on and pump-off decay rates is immune to empty-CRD drifts and spectral overlaps from non-target molecular transitions. The immunity to the ringdown drifts allows longer signal-averaging and, thus, higher detection sensitivity. The ability to compensate for the background absorption enhances the detection selectivity in spectrally congested regions. Our technique is well-suited for trace-detection in the mid-IR region, where pump–probe schemes based on strong rovibrational transitions can be applied. In this work, two-color CRD detection is implemented on a ladder-type, three-level system based on the N 2 O, ν 3 = 1 ← 0, P(19) (pump) and ν 3 = 2 ← 1, R(18) (probe), rovibrational transitions. By frequency-locking two-quantum cascade lasers to the p-polarization (pump, Finesse = 5280) and s-polarization (probe, Finesse = 67 700) cavity modes, we achieve high intracavity pump power (36 W) and high probe ringdown rates (>2 kHz). The observed two-color spectra are simulated by a density-matrix, three-level system model that is solved under the constraints of the cavity resonance conditions. In addition to its background compensation capability, experimental flexibility in the selection of pump–probe schemes and signal insensitivity to intracavity laser power are further features that enhance the utility of our technique for mid-IR trace-detection.

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

A graphene edge-mediated quantum gate

We propose a quantum gate architecture that allows for the systematic control of the effective exchange interactions between magnetic impurities embedded in nanoscale graphene flakes connected by a gated bridge. The entanglement between the magnetic moment and the edge states of the fragments is used to electrostatically tune the exchange interaction from ferro to antiferromagnetic by merely changing the bridge's carrier density. By characterizing the effects of size and coupling parameters, we explore different operation regimes of this device by means of exact calculations with the density matrix renormalization group. We analyze the results utilizing a simplified model that accounts for the main many-body mechanisms. Finally, we discuss how to use arrays of these devices to build quantum simulators for quantum many-body Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

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

Efficient evaluation of the Breit operator in the Pauli spinor basis

The frequency-independent Coulomb–Breit operator gives rise to the most accurate treatment of two-electron interaction in the non-quantum-electrodynamics regime. The Breit interaction in the Coulomb gauge consists of magnetic and gauge contributions. The high computational cost of the gauge term limits the application of the Breit interaction in relativistic molecular calculations. In this work, we apply the Pauli component integral–density matrix contraction scheme for gauge interaction with a maximum spin- and component separation scheme. We also present two different computational algorithms for evaluating gauge integrals. One is the generalized Obara–Saika algorithm, where the Laplace transformation is used to transform the gauge operator into Gaussian functions and the Obara–Saika recursion is used for reducing the angular momentum. The other algorithm is the second derivative of Coulomb interaction evaluated with Rys-quadrature. This work improves the efficiency of performing Dirac–Hartree–Fock with the variational treatment of Breit interaction for molecular systems. In conclusion, we use this formalism to examine relativistic trends in the Periodic Table and analyze the relativistic two-electron interaction contributions in heavy-element complexes.

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

Molecular dipole moment learning via rotationally equivariant derivative kernels in molecular-orbital-based machine learning

This study extends the accurate and transferable molecular-orbital-based machine learning (MOB-ML) approach to modeling the contribution of electron correlation to dipole moments at the cost of Hartree–Fock computations. A MOB pairwise decomposition of the correlation part of the dipole moment is applied, and these pair dipole moments could be further regressed as a universal function of MOs. The dipole MOB features consist of the energy MOB features and their responses to electric fields. An interpretable and rotationally equivariant derivative kernel for Gaussian process regression (GPR) is introduced to learn the dipole moment more efficiently. The proposed problem setup, feature design, and ML algorithm are shown to provide highly accurate models for both dipole moments and energies on water and 14 small molecules. To demonstrate the ability of MOB-ML to function as generalized density-matrix functionals for molecular dipole moments and energies of organic molecules, we further apply the proposed MOB-ML approach to train and test the molecules from the QM9 dataset. The application of local scalable GPR with Gaussian mixture model unsupervised clustering GPR scales up MOB-ML to a large-data regime while retaining the prediction accuracy. In addition, compared with the literature results, MOB-ML provides the best test mean absolute errors of 4.21 mD and 0.045 kcal/mol for dipole moment and energy models, respectively, when training on 110 000 QM9 molecules. The excellent transferability of the resulting QM9 models is also illustrated by the accurate predictions for four different series of peptides.

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

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2022

Magnetic phases of spatially modulated spin-1 chains in Rydberg excitons: Classical and quantum simulations

Here, in this work, we study the magnetic phases of a spatially modulated chain of spin-1 Rydberg excitons. Using the Density Matrix Renormalization Group (DMRG) technique, we study various magnetic and topologically nontrivial phases using both single-particle properties, such as local magnetization and quantum entropy, and many-body ones, such as pair-wise Néel and long-range string correlations. In particular, we investigate the emergence and robustness of the Haldane phase, a topological phase of anti-ferromagnetic spin-1 chains. Furthermore, we devise a hybrid quantum algorithm employing restricted Boltzmann machine to simulate the ground state of such a system that shows very good agreement with the results of exact diagonalization and DMRG.

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

Two-color polarization spectroscopy measurements of Zeeman state-to-state collision induced transitions of nitric oxide in binary gas mixtures

We investigated collision induced transitions in the (0, 0) band of the A2Σ+–X2Π electronic transition of nitric oxide (NO) using two-color polarization spectroscopy (TCPS). Two sets of TCPS spectra for 1% NO, diluted in different buffer gases at 295 K and 1 atm, were obtained with the pump beam tuned to the R11(11.5) and OP12(1.5) transitions. The buffer gases were He, Ar, and N2. The probe was scanned while the pump beam was tuned to the line center. Theoretical TCPS spectra, calculated by solving the density matrix formulation of the time-dependent Schrödinger wave equation, were compared with the experimental spectra. A collision model based on the modified exponential-gap law was used to model the rotational level-to-rotational level collision dynamics. A model for collisional transfer from an initial to a final Zeeman state was developed based on the difference in cosine of the rotational quantum number J projection angle with the z-axis for the two Zeeman states. Rotational energy transfer rates and Zeeman state collisional dynamics were varied to obtain good agreement between theory and experiment for the two different TCPS pump transitions and for the three different buffer gases. One key finding, in agreement with quasi-classical trajectory calculations, is that the spin-rotation changing transition rate in the A2Σ+ level of NO is almost zero for rotational quantum numbers ≥8. It was necessary to set this rate to near zero to obtain agreement with the TCPS spectra.

Chemistry↗
DOE OSTI2023

Generalized spin σ -SCF method

We introduce a generalization of the σ-SCF method to approximate noncollinear spin ground and excited single-reference electronic states by minimizing the Hamiltonian variance. The new method is based on the σ-SCF method, originally proposed by Ye et al. [J. Chem. Phys. 147, 214104 (2017)], and provides a prescription to determine ground and excited noncollinear spin states on an equal footing. Our implementation was carried out utilizing an initial simulated annealing stage followed by a mean-field iterative self-consistent approach to simplify the cumbersome search introduced by generalizing the spin degrees of freedom. The simulated annealing stage ensures a broad exploration of the Hilbert space spanned by the generalized spin single-reference states with random complex element-wise rotations of the generalized density matrix elements in the simulated annealing stage. The mean-field iterative self-consistent stage employs an effective Fockian derived from the variance, which is utilized to converge tightly to the solutions. This process helps us to easily find complex spin structures, avoiding manipulating the initial guess. As proof-of-concept tests, we present results for Hn (n = 3–7) planar rings and polyhedral clusters with geometrical spin frustration. We show that most of these systems have noncollinear spin excited states that can be interpreted in terms of geometric spin frustration. These states are not directly targeted by energy minimization methods, which are meant to converge to the ground state. This stresses the capability of the σ-SCF methodology to find approximate noncollinear spin structures as mean-field excited states.

Chemistry↗
DOE OSTI2024

Noise reduction of stochastic density functional theory for metals

Density Functional Theory (DFT) has become a cornerstone in the modeling of metals. However, accurately simulating metals, particularly under extreme conditions, presents two significant challenges. First, simulating complex metallic systems at low electron temperatures is difficult due to their highly delocalized density matrix. Second, modeling metallic warm-dense materials at very high electron temperatures is challenging because it requires the computation of a large number of partially occupied orbitals. This study demonstrates that both challenges can be effectively addressed using the latest advances in linear-scaling stochastic DFT methodologies. Despite the inherent introduction of noise into all computed properties by stochastic DFT, this research evaluates the efficacy of various noise reduction techniques under different thermal conditions. Our observations indicate that the effectiveness of noise reduction strategies varies significantly with the electron temperature. Furthermore, we provide evidence that the computational cost of stochastic DFT methods scales linearly with system size for metal systems, regardless of the electron temperature regime.

Chemistry↗
DOE OSTI2024

Efficient simulation of open quantum systems coupled to a reservoir through multiple channels

It is challenging to simulate open quantum systems that are connected to a reservoir through multiple channels. For example, vibrations may induce fluctuations in both energy gaps and electronic couplings, which represent two independent channels of system–bath couplings. Systems of this kind are ubiquitous in the processes of excited state radiationless decay. Combined with density matrix renormalization group (DMRG) and matrix product states (MPS) methods, we develop an interaction-picture chain mapping strategy for vibrational reservoirs to simulate the dynamics of these open systems, resulting in time-dependent spatially local system–bath couplings in the chain-mapped Hamiltonian. This transformation causes the entanglement generated by the system–bath interactions to be restricted within a narrow frequency window of vibrational modes, enabling efficient DMRG/MPS dynamical simulations. We demonstrate the utility of this approach by simulating singlet fission dynamics using a generalized spin-boson Hamiltonian with both diagonal and off-diagonal system–bath couplings. This approach generalizes an earlier interaction-picture chain mapping scheme, allowing for efficient and exact simulation of systems with multi-channel system–bath couplings using matrix product states, which may further our understanding of nonlocal exciton–phonon couplings in exciton transport and the non-Condon effect in energy and electron transfer.

Chemistry↗