Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Communication: A difference density picture for the self-consistent field ansatz

We formulate self-consistent field (SCF) theory in terms of an interaction picture where the working variable is the difference density matrix between the true system and a corresponding superposition of atomic densities. As the difference density matrix directly represents the electronic deformations inherent in chemical bonding, this “difference self-consistent field (dSCF)” picture provides a number of significant conceptual and computational advantages. We show that this allows for a stable and efficient dSCF iterative procedure with wholly single-precision Coulomb and exchange matrix builds. We also show that the dSCF iterative procedure can be performed with aggressive screening of the pair space. These approximations are tested and found to be accurate for systems with up to 1860 atoms and >10 000 basis functions, providing for immediate overall speedups of up to 70% in the heavily optimized TERACHEM SCF implementation.

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

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2 p and 3 p 1 h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

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

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

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

N -representability violations in truncated equation-of-motion coupled-cluster methods

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. Furthermore, these observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

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

On-the-fly machine learned force fields for the study of warm dense matter: Application to diffusion and viscosity of CH

We develop a framework for on-the-fly machine learned force field (MLFF) molecular dynamics (MD) simulations of warm dense matter (WDM). In particular, we employ an MLFF scheme based on the kernel method and Bayesian linear regression, with the training data generated from the Kohn–Sham density functional theory (DFT) using the Gauss spectral quadrature method, within which we calculate energies, atomic forces, and stresses. We verify the accuracy of the formalism by comparing the predicted properties of warm dense carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that ab initio MD simulations of WDM can be accelerated by up to three orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the diffusion coefficients and shear viscosity of CH at a density of 1 g/cm3 and temperatures in the range of 75 000–750 000 K. We find that the self- and inter-diffusion coefficients and the viscosity obey a power law with temperature, and that the diffusion coefficient results suggest a weak coupling between C and H in CH. In addition, we find agreement within standard deviation with previous results for C and CH but disagreement for H, demonstrating the need for ab initio calculations as presented here.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

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

Magnetic field-induced intermediate quantum spin liquid with a spinon Fermi surface

The Kitaev model with an applied magnetic field in the H ∥ [ 111 ] direction shows two transitions: from a nonabelian gapped quantum spin liquid (QSL) to a gapless QSL at H c 1 ≃ 0.2 K and a second transition at a higher field H c 2 ≃ 0.35 K to a gapped partially polarized phase, where K is the strength of the Kitaev exchange interaction. We identify the intermediate phase to be a gapless U(1) QSL and determine the spin structure function S ( k ) and the Fermi surface ϵ F S ( k ) of the gapless spinons using the density matrix renormalization group (DMRG) method for large honeycomb clusters. Further calculations of static spin-spin correlations, magnetization, spin susceptibility, and finite temperature-specific heat and entropy corroborate the gapped and gapless nature of the different field-dependent phases. In the intermediate phase, the spin-spin correlations decay as a power law with distance, indicative of a gapless phase.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Exact matrix product state representation and convergence of a fully correlated electronic wavefunction in the infinite-basis limit

Here In this paper we present the exact representation of a fully correlated electronic wavefunction as the single-particle basis approaches completeness. It consists of a half-infinite chain of matrices of exponentially increasing size. The complete basis limit is illustrated numerically using the density-matrix renormalization-group method by computing the core-valence entanglement in the C 2 ground state in increasing subsets of cc-pVTZ and pVQZ bases until convergence is reached.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nature of quantum spin liquids of the S = 1 2 Heisenberg antiferromagnet on the triangular lattice: A parallel DMRG study

Here we study the ground-state properties of the quantum spin liquid (QSL) phases of the spin-1/2 antiferromagnetic Heisenberg model on the triangular lattice with nearest- (J 1 ), next-nearest- (J 2 ), and third-neighbor (J 3 ) interactions by using density-matrix renormalization group (DMRG) method. By combining parallel DMRG with SU(2) spin rotational symmetry, we are able to obtain accurate results on large cylinders with length up to L x =48 and circumference L y =6–12. Our results suggest that the QSL phase of the J 1 -J 2 Heisenberg model is gapped which is characterized by the absence of gapless mode, short-range spin-spin and dimer-dimer correlations. In the presence of J 3 interaction, we find that a new critical QSL with a single gapless mode emerges. While both spin-spin and scalar chiral-chiral correlations are short-ranged, dimer-dimer correlations are quasi-long-ranged which decays as a power-law at long distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Influence of extended interactions on spin dynamics in one-dimensional cuprates

Quasi-one-dimensional (1D) materials provide a unique platform for understanding the importance and influence of extended interactions on the physics of strongly correlated systems due to their relative structural simplicity and the existence of powerful theoretical tools well-adapted to one spatial dimension. Recently, this was highlighted by anomalous observations in the single-particle spectral function A(q, ω) of 1D cuprate chain compounds, measured by angle-resolved photoemission spectroscopy (ARPES), which were explained by the presence of a long-range attractive interaction. Such an extended interaction should leave its fingerprints on other observables, notably the dynamical spin structure factor S(q, ω), measured by neutron scattering or resonant inelastic x-ray scattering (RIXS). Starting from a simple Hubbard Hamiltonian in 1D and using time-dependent density matrix renormalization group (tDMRG) methods, we show that the presence of long-range attractive coupling, directly through an instantaneous Coulomb interaction V or retarded electron phonon (el-ph) coupling, can introduce significant spectral weight redistribution in S(q, ω) across a wide range of doping. Here, this underscores the significant impact that extended interactions can have on dynamical correlations among particles, and the importance of properly incorporating this influence in modeling. Our results demonstrate that S(q, ω) can provide a sensitive experimental constraint, which complements ARPES measurements, in identifying key interactions in 1D cuprates, beyond the standard Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Doping evolution of charge and spin excitations in two-leg Hubbard ladders: Comparing DMRG and FLEX results [Doping evolution of charge and spin excitations in two-leg Hubbard ladders: Comparing DMRG and RPA+FLEX results]

Here, we study the magnetic and charge dynamical response of a Hubbard model in a two-leg ladder geometry using the density matrix renormalization group (DMRG) method and the random phase approximation within the fluctuation-exchange approximation (FLEX). Our calculations reveal that FLEX can capture the main features of the magnetic response from weak up to intermediate Hubbard repulsion for doped ladders, when compared with the numerically exact DMRG results. However, while at weak Hubbard repulsion both the spin and charge spectra can be understood in terms of weakly interacting electron-hole excitations across the Fermi surface, at intermediate coupling DMRG shows gapped spin excitations at large momentum transfer that remain gapless within the FLEX approximation. For the charge response, FLEX can only reproduce the main features of the DMRG spectra at weak coupling and high doping levels, while it shows an incoherent character away from this limit. Overall, our analysis shows that FLEX works surprisingly well for spin excitations at weak and intermediate Hubbard U values even in the difficult low-dimensional geometry such as a two-leg ladder. Finally, we discuss the implications of our results for neutron scattering and resonant inelastic x-ray scattering experiments on two-leg ladder cuprate compounds.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Density matrix renormalization group description of the island of inversion isotopes F 28 – 33

Recent experiments have confirmed that the neutron-rich isotopes 28,29 F belong to the so-called island of inversion (IOI), a region of the nuclear chart around Z = 10 and N = 20 where nuclear structure deviates from the standard shell model predictions due to deformation and continuum effects. However, while the general principles leading to the IOI are relatively well understood, the details of the low-lying structure of the exotic fluorine isotopes 28–33 F are basically unknown. In this study, we perform large-scale shell model calculations including continuum states to investigate the properties of the neutron-rich isotopes 25–33 F, from a core of 24 O and using an effective two-body interaction with a small number of adjustable parameters in the central and tensor channels. We develop two models adjusted on experimentally confirmed states in 25,26 O and 25–27 F based on different assumptions concerning the positions of the neutron 0d 3/2 and 1p 3/2 shells, and solve the many-body problem using the density matrix renormalization group (DMRG) method for open quantum systems in an sd–fp model space. We obtain the first detailed spectroscopy of 25–33F in the continuum and show how the interplay between continuum effects and deformation explains the recent data on 28,29 F. Several deformed one- and two-neutron halo states are predicted in 29,31 F, and we provide some information about the possible structure of the heaviest fluorine isotopes. We also suggest several experimental studies of interest to constraint models and test the present predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Room-Temperature Mechanical Resonator with a Single Added or Subtracted Phonon

A room-temperature mechanical oscillator undergoes thermal Brownian motion with an amplitude much larger than the amplitude associated with a single phonon of excitation. This motion can be read out and manipulated using laser light using a cavity-optomechanical approach. By performing a strong quantum measurement (i.e., counting single photons in the sidebands imparted on a laser), we herald the addition and subtraction of single phonons on the 300 K thermal motional state of a 4 GHz mechanical oscillator. Further, to understand the resulting mechanical state, we implement a tomography scheme and observe highly non-Gaussian phase-space distributions. Using a maximum likelihood method, we infer the density matrix of the oscillator, and we confirm the counterintuitive doubling of the mean phonon number resulting from phonon addition and subtraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2015

Kohn–Sham exchange-correlation potentials from second-order reduced density matrices

We describe a practical algorithm for constructing the Kohn–Sham exchange-correlation potential corresponding to a given second-order reduced density matrix. Unlike conventional Kohn–Sham inversion methods in which such potentials are extracted from ground-state electron densities, the proposed technique delivers unambiguous results in finite basis sets. The approach can also be used to separate approximately the exchange and correlation potentials for a many-electron system for which the reduced density matrix is known. The algorithm is implemented for configuration-interaction wave functions and its performance is illustrated with numerical examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Real-time density-matrix coupled-cluster approach for closed and open systems at finite temperature

We extend the coupled-cluster method to correlated quantum dynamics of both closed and open systems at finite temperatures using the thermofield formalism. The approach expresses the time-dependent density matrix in an exponential ansatz and describes time-evolution along the Keldysh path contour. A distinct advantage of the approach is exact trace-preservation as a function of time, ensuring conservation of probability and particle number. Furthermore, the method avoids the computation of correlated bra-states, simplifying the computational implementation. We develop the method in a thermal quasiparticle representation, which allows seamless connection to the projection method and diagrammatic techniques of the traditional coupled-cluster formalism. For comparison, we also apply the thermofield framework to the density-matrix renormalization-group method to obtain reference results for closed and open systems at finite temperature. We test the singles and doubles approximation to the density-matrix coupled-cluster method on the correlated electronic dynamics of the single-impurity Anderson model, demonstrating that the new method successfully captures the correlated dynamics of both closed systems at finite temperature and driven-dissipative open systems. This encouraging performance motivates future applications to nonequilibrium quantum many-body dynamics in realistic systems.

Shushkov, Philip (ORCID:0000000266136117)↗