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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Multireference Density Functional Theory with Generalized Auxiliary Systems for Ground and Excited States

To describe static correlation, we create a new approach to density functional theory (DFT), which uses a generalized auxiliary system that is of a different symmetry, such as particle number or spin, from that of the physical system. The total energy of the physical system consists of two parts: the energy of the auxiliary system, which is determined with a chosen density functional approximation (DFA), and the excitation energy from an approximate linear response theory that restores the symmetry to that of the physical system, thus rigorously leading to a multideterminant description of the physical system. The electron density of the physical system is different from that of the auxiliary system and is uniquely determined from the functional derivative of the total energy with respect to the external potential. Our energy functional is thus an implicit functional of the physical system density, but an explicit functional of the auxiliary system density. We demonstrate that the total energy minimum and stationary states, describing the ground and excited states of the physical system, can be obtained by a self-consistent optimization with respect to the explicit variable, the generalized Kohn–Sham noninteracting density matrix. We have developed the generalized optimized effective potential method for the self-consistent optimization. Among options of the auxiliary system and the associated linear response theory, reformulated versions of the particle–particle random phase approximation (pp-RPA) and the spin-flip time-dependent density functional theory (SF-TDDFT) are selected for illustration of principle. Numerical findings show that our multireference DFT successfully describes static correlation in bond dissociation and double bond rotation.

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

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗
DOE OSTI2023

Ultrafast Molecular Frame Quantum Tomography

Here, we develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH 3 ). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH 3 , we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Competing phases, phase separation, and coexistence in the extended one-dimensional bosonic Hubbard model

Here, we study the phase diagram of the one-dimensional bosonic Hubbard model with contact (U) and near neighbor (V ) interactions focusing on the gapped Haldane insulating (HI) phase which is characterized by an exotic nonlocal order parameter. The parameter regime (U, V and μ) where this phase exists and how it competes with other phases such as the supersolid (SS) phase, is incompletely understood. We use the Stochastic Green Function quantum Monte Carlo algorithm as well as the density matrix renormalization group to map out the phase diagram. The HI exists only at = 1, the SS phase exists for a very wide range of parameters (including commensurate fillings) and displays power law decay in the one body Green function were our main conclusions. Additionally, we show that at fixed integer density, the system exhibits phase separation in the (U, V ) plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Superconductivity in the doped t – J model: Results for four-leg cylinders

Here, we report a density-matrix renormalization group study of the lightly doped t–J model on a four-leg cylinder with doped hole concentrations per site δ = 5%–12.5%. By keeping an unusually large number of states and long system sizes, we are able to accurately document the interplay between superconductivity, spin, and charge-density-wave orders. The long-distance behavior is consistent with that of a Luther-Emery liquid with a spin gap and power-law charge-density-wave and superconducting correlations. This is the widest t–J or Hubbard system in which power-law superconducting correlations have been established.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Nonequilibrium Steady-State Transport in Quantum Impurity Models: A Thermofield and Quantum Quench Approach Using Matrix Product States

Here, the numerical renormalization group (NRG) is tailored to describe interacting impurity models in equilibrium, but it faces limitations for steady-state nonequilibrium, arising, e.g., due to an applied bias voltage. We show that these limitations can be overcome by describing the thermal leads using a thermofield approach, integrating out high energy modes using NRG, and then treating the nonequilibrium dynamics at low energies using a quench protocol, implemented using the time-dependent density matrix renormalization group. This yields quantitatively reliable results for the current (with errors ≲3%) down to the exponentially small energy scales characteristic of impurity models. We present results of benchmark quality for the temperature and magnetic field dependence of the zero-bias conductance peak for the single-impurity Anderson model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Intertwined order in a frustrated four-leg t - J cylinder

Here, we report a density-matrix renormalization group study of the t–J model with nearest (t 1 and J 1 ) and next-nearest (t 2 and J 2 ) interactions on a four-leg cylinder with concentration δ=1/8 of doped holes. We observe an astonishingly complex interplay between uniform d-wave superconductivity (SC) and strong spin and charge-density wave ordering tendencies (SDW and CDW). Depending on parameters, the CDWs can be commensurate with period 4 or 8. By comparing the charge ordering vectors with 2k F , we rule out Fermi surface nesting-induced density wave order in our model. Magnetic frustration (i.e., J 2 /J 1 ~1/2) significantly quenches SDW correlations with little effect on the CDW. Typically, the SC order is strongly modulated at the CDW ordering vector and exhibits d-wave symmetry around the cylinder. There is no evidence of a near-degenerate tendency to pair-density wave (PDW) ordering, charge 4e SC, or orbital current order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Linear scaling pseudo Fermi-operator expansion for fractional occupation

Recursive Fermi-operator expansion methods for the calculation of the idempotent density matrix are valid only at zero electronic temperature with integer occupation numbers. We show how such methods can be modified to include fractional occupation numbers of an approximate or pseudo Fermi-Dirac distribution and how the corresponding entropy term of the free energy is calculated. The proposed methodology is demonstrated and evaluated for different electronic structure methods including density functional tight-binding theory, Kohn-Sham density functional theory using numerical orbitals, and quantum chemistry Hartree-Fock theory using Gaussian basis functions.

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

Pair density wave in the doped three-band Hubbard model on two-leg square cylinders

A pair density wave (PDW) is a superconducting (SC) state with a spatially modulated order parameter. Although much is known about the properties of the PDW state, its realization in microscopic models with divergent susceptibility has been challenging. Here, in this work, we report a density-matrix renormalization group study of a three-band Hubbard model (also known as the Emery model) for cuprates on long two-leg square cylinders. Upon light doping, we find that the ground state of the system is consistent with that of a PDW state with mutually commensurate and power-law SC, charge (CDW), and spin (SDW) density wave correlations. The SC correlations are dominant between neighboring Cu sites with d-wave pairing symmetry. Interestingly, we find that the near-neighbor interactions, especially the near-neighbor attractive V pd interaction between neighboring Cu and oxygen sites, can notably enhance the SC correlations while simultaneously suppressing the CDW correlations. For a modestly strong attractive V pd , the SC correlations become quasi-long-ranged with a divergent PDW susceptibility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Improved Kelbg Potentials for Z > 1 and Application to Carbon Plasmas

In this work, we present a general form for the electron‐ion diffractive potential derived from the quantum pair density matrix and fit to the improved Kelbg potential for atomic numbers up to $Z = 54$. We apply classical molecular dynamics using the improved Kelbg potential for carbon with various forms of the Pauli potential to compute internal energies and pressures for hot, dense plasma conditions. Our results are compared to an equation of state model based on path integral Monte Carlo and density functional theory simulations to examine the extent to which the improved Kelbg potential reproduces the internal energy and pressure of carbon plasmas. The regions of validity for carbon agree generally with those derived previously for hydrogen once pressure ionization effects are incorporated. Based on our carbon results and previously published hydrogen studies, we discuss the general applicability and limitations of these potentials for equation of state studies in warm dense matter and high energy density plasmas.

general physics↗
DOE OSTI2018

Large-Scale Electron Correlation Calculations: Rank-Reduced Full Configuration Interaction

Here, we present the rank-reduced full configuration interaction (RR-FCI) method, a variational approach for the calculation of extremely large full configuration interaction (FCI) wave functions. In this report, we show that RR-FCI can provide ground state singlet and triplet energies within kcal/mol accuracy of full CI (FCI) with computational effort scaling as the square root of the number of determinants in the CI space (compared to conventional FCI methods which scale linearly with the number of determinants). Fast graphical processing unit (GPU) accelerated projected σ = Hc matrix–vector product formation enables calculations with configuration spaces as large as 30 electrons in 30 orbitals, corresponding to an FCI calculation with over 2.4 × 10 16 configurations. We apply this method in the context of complete active space configuration interaction calculations to acenes with 2–5 aromatic rings, comparing absolute energies against FCI when possible and singlet/triplet excitation energies against both density matrix renormalization group (DMRG) and experimental results. The dissociation of molecular nitrogen was also examined using both FCI and RR-FCI. In each case, we found that RR-FCI provides a low cost alternative to FCI, with particular advantages when relative energies are desired.

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

Tuning topological phase and quantum anomalous Hall effect by interaction in quadratic band touching systems

Interaction-driven topological phases significantly enrich the class of topological materials and thus are of great importance. Here in this paper, we study the phase diagram of interacting spinless fermions filling the two-dimensional checkerboard lattice with a quadratic band touching (QBT) point. By developing new diagnosis based on the state-of-the-art density-matrix renormalization group and exact diagonalization, we determine accurate quantum phase diagram for such a system at half-filling with three distinct phases. For weak nearest-neighboring interactions, we demonstrate the instability of the QBT towards an interaction-driven spontaneous quantum anomalous Hall (QAH) effect. For strong interactions, the system breaks the rotational symmetry realizing a nematic charge-density-wave (CDW) phase. Interestingly, for intermediate interactions we discover a symmetry-broken bond-ordered critical phase sandwiched in between the QAH and CDW phases, which splits the QBT into two Dirac points driven by interaction. Instead of the direct transition between QAH and CDW phases, our identification of an intermediate phase sheds new light on the theoretical understanding of the interaction-driven phases in QBT systems.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Nanoporous Cu-C composites based on carbon-nanotube aerogels

Current synthesis methods of nanoporous Cu–C composites offer limited control of the material composition, structure, and properties, particularly for large Cu loadings of ≳20 wt%. Here, we describe two related approaches to realize novel nanoporous Cu–C composites based on the templating of recently developed carbon-nanotube aerogels (CNT-CAs). Our first approach involves the trapping of Cu nanoparticles while CNT-CAs undergo gelation. This method yields nanofoams with relatively high densities of ≳65 mg cm -3 for Cu loadings of ≳10 wt%. Our second approach overcomes this limitation by filling the pores of undoped CNT-CA monoliths with an aqueous solution of CuSO 4 followed by (i) freeze-drying to remove water and (ii) thermal decomposition of CuSO 4 . With this approach, we demonstrate Cu–C composites with a C matrix density of -25 mg cm -3 and Cu loadings of up to 70 wt%. These versatile methods could be extended to fabricate other nanoporous metal–carbon composite materials geared for specific applications.

Charnvanichborikarn, S.↗
DOE OSTI2019

Energy window stochastic density functional theory

We report that linear scaling density functional theory is important for understanding electronic structure properties of nanometer scale systems. Recently developed stochastic density functional theory can achieve linear or even sublinear scaling for various electronic properties without relying on the sparsity of the density matrix. The basic idea relies on projecting stochastic orbitals onto the occupied space by expanding the Fermi-Dirac operator and repeating this for N χ stochastic orbitals. Often, a large number of stochastic orbitals are required to reduce the statistical fluctuations (which scale as N$-1/2\atop{χ}$) below a tolerable threshold. In this work, we introduce a new stochastic density functional theory that can efficiently reduce the statistical fluctuations for certain observable and can also be integrated with an embedded fragmentation scheme. The approach is based on dividing the occupied space into energy windows and projecting the stochastic orbitals with a single expansion onto all windows simultaneously. This allows for a significant reduction of the noise as illustrated for bulk silicon with a large supercell. Finally, we also provide theoretical analysis to rationalize why the noise can be reduced only for a certain class of ground state properties, such as the forces and electron density.

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

Ground-state phase diagram and superconductivity of the doped Hubbard model on six-leg square cylinders

Here, we have studied the ground state properties of Hubbard model on long six-leg square cylinders with doped hole concentration per site 5.55%≤δ≤12.5% using density-matrix renormalization group. By keeping a large number of states for long system sizes, we find that the nature of the ground state is remarkably sensitive to the presence of next-nearest-neighbor electron hopping t'. In the positive t' side, we find a robust d-wave superconducting (SC) phase characterized by coexisting quasi-long-range SC and charge density wave (CDW) correlations. Without t' the ground state forms an insulating stripe phase with long-range CDW order but short-range spin-spin and SC correlations. In stark contrast to four-leg cylinders, our results show that the lightly doped Hubbard model on six-leg cylinders remains insulating in the negative t' side where the SC correlations decay exponentially with short correlation lengths. In the larger negative t' side, the doped holes form a novel holon Wigner crystal with one doped hole per emergent unit cell and short-range spin-spin correlations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various density functional theory packages

We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gapless spin liquid and pair density wave of the Hubbard model on three-leg triangular cylinders

We study the ground state properties of the Hubbard model on three-leg triangular cylinders using large-scale density-matrix renormalization group simulations. At half-filling, we identify an intermediate gapless spin liquid phase, which has one gapless spin mode and algebraic spin–spin correlations but exponential decay scalar chiral–chiral correlations, between a metallic phase at weak coupling and Mott insulating dimer phase at strong interaction. Upon light doping the gapless spin liquid, the system exhibits power-law charge-density-wave (CDW) correlations but short-range single-particle, spin–spin, and chiral–chiral correlations. Similar to CDW correlations, the superconducting correlations also decay in power-law but oscillate in sign as a function of distance, which is consistent with the striped pair-density wave. When further doping the gapless spin liquid phase or doping the dimer order phase, another phase takes over, which has similar CDW correlations but all other correlations decay exponentially.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗