SEARCH · Engineering Papers

Results for “DENSITY OF STATES”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗

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 OSTI2024

What Controls Crystal Diversity and Microphysical Variability in Cirrus Clouds?

Abstract Variability of ice microphysical properties like crystal size and density in cirrus clouds is important for climate through its impact on radiative forcing, but challenging to represent in models. For the first time, recent laboratory experiments of particle growth (tied to crystal morphology via deposition density) are combined with a state‐of‐the‐art Lagrangian particle‐based microphysics model in large‐eddy simulations to examine sources of microphysical variability in cirrus. Simulated particle size distributions compare well against balloon‐borne observations. Overall, microphysical variability is dominated by variability in the particles' thermodynamic histories. However, diversity in crystal morphology notably increases spatial variability of mean particle size and density, especially at mid‐levels in the cloud. Little correlation between instantaneous crystal properties and supersaturation occurs even though the modeled particle morphology is directly tied to supersaturation based on laboratory measurements. Thus, the individual thermodynamic paths of each particle, not the instantaneous conditions, control the evolution of particle properties.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Block2 : A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.

Chemistry↗

DOE OSTI2024

Learning stochastic dynamics and predicting emergent behavior using transformers

We show that a neural network originally designed for language processing can learn the dynamical rules of a stochastic system by observation of a single dynamical trajectory of the system, and can accurately predict its emergent behavior under conditions not observed during training. We consider a lattice model of active matter undergoing continuous-time Monte Carlo dynamics, simulated at a density at which its steady state comprises small, dispersed clusters. We train a neural network called a transformer on a single trajectory of the model. The transformer, which we show has the capacity to represent dynamical rules that are numerous and nonlocal, learns that the dynamics of this model consists of a small number of processes. Forward-propagated trajectories of the trained transformer, at densities not encountered during training, exhibit motility-induced phase separation and so predict the existence of a nonequilibrium phase transition. Transformers have the flexibility to learn dynamical rules from observation without explicit enumeration of rates or coarse-graining of configuration space, and so the procedure used here can be applied to a wide range of physical systems, including those with large and complex dynamical generators.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Tensor-polarized parton density in the $N → Δ$ transition

The generalized parton distributions for transitions between baryon states with different masses have a forward limit in which they behave as parton densities (light-front momentum transfer Δ + , Δ T = 0, energy transfer Δ - ≠ 0). These "transition parton densities" can realize spin/isospin quantum numbers not accessible in the ground-state nucleon. The N → Δ transition gives rise to a new parton density proportional to the 1/2 → 3/2 spin transition tensor. Its properties are derived, and its magnitude is estimated in the chiral quark-soliton model based on the large-N c limit of QCD.

Delta isobar↗

DOE OSTI2021

Stoichiometric irreversibility of aged garnet electrolytes

Solid-state lithium batteries (SSLBs) have been regarded as one of the next-generation energy storage systems. With the adoption of solid-state electrolytes (SSEs) and lithium metal anodes, SSLBs enable higher energy density and more reliable safety than the state-of-the-art lithium-ion batteries. Among potential SSEs, the cation-doped Li 7 La 3 Zr 2 O 12 (LLZO) is promising for its high ionic conductivity (~10 –3 S cm –2 ) at room temperature and high stability with Li metal anode. However, the storage of doped LLZO in the ambient condition suffers the aging effect, including the structural transition (i.e. low-temperature cubic form) and the stoichiometric changes (i.e. Li 2 CO 3 ). These changes are detrimental to LLZO ionic conductivity and interfacial stability in SSLBs. To this end, in this study we are motivated to investigate the structural and stoichiometric reversibility of aged LLZO during thermal treatment. With the help of an in-situ synchrotron-based high-energy X-ray diffraction technique, our experiments revealed that the LLZO powders became a low-temperature cubic phase when exposed to the ambient condition for an extended period of time. A high temperature cubic form can be restored after a thermal treatment of the aged LLZO powder, regardless of the type of dopant. However, the restoration of the stoichiometry remained a challenge, and the degree of the restoration showed a clear dependence on the dopant chemistry.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Strong long-wavelength electron-phonon coupling in Ta 2 ⁢Ni⁡(Se,S) 5

The search for intrinsic excitonic insulators (EI) has long been confounded by coexisting electron–phonon coupling in bulk materials. Although the ground state of an EI may be difficult to differentiate from density-wave orders or other structural instabilities, excited states offer distinctive signatures. One way to provide clarity is to directly inspect the phonon spectral function for long wavelength broadening caused by phonon interaction with the high velocity EI phason. Here, in this study, we report that the quasi-one-dimensional (quasi-1D) EI candidate Ta 2 NiSe 5 shows extremely anisotropic phonon broadening and softening in the semimetallic normal state. In contrast, such behavior is completely absent in the broken symmetry state of Ta 2 NiSe 5 and in the isostructural Ta 2 NiS 5 , where the latter has a fully gapped normal state. By contrasting the expected phonon lifetimes in the BCS and BEC limits of a putative EI, our results suggest that the phase transition in Ta 2 Ni(Se,S) 5 family is closely related to strong interband electron–phonon coupling. We experimentally determine the dimensionless coupling $\frac{g}{ω_0}$ ∼ 10, revealing Ta 2 Ni(Se,S) 5 as a rare “ultrastrong coupling” material.

Kang, Zhibo [Yale University, New Haven, CT (Unite↗

DOE OSTI2023

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree–Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Chemistry↗

DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Practically Accessible All-Solid-State Batteries Enabled by Organosulfide Cathodes and Sulfide Electrolytes

The combination of organic electrode materials and sulfide electrolytes is expected to enable the development of all-solid-state organic batteries featuring high energy density and long-term sustainability. In this work, thiuram hexasulfide is reported as a low-cost and high-capacity organic cathode for solid-state batteries based on sulfide electrolytes, delivering a capacity of ~600 mAh g -1 and retaining 80.8 % capacity after 500 cycles. An electrochemically reversible change of the cathode interface was revealed upon cycling. Full cell displays an oscillating stress change up to 0.6 MPa during cycling, predominated by the anode side. The energy density is 1140 Wh kg -1 at the material level and 376 Wh kg -1 at the electrode level, which are among the best-reported organic cathodes to date. A high areal capacity of 10.4 mAh cm -2 is reached with a high mass loading cathode. A dry-film approach is further explored to manufacture sheet-type cells. Free-standing electrolyte film is merely ~48 μm thick and demonstrates an ultralow areal resistance of 3.9 Ω cm 2 , significantly boosts the cell-level energy density and reduces the cell internal resistance.

25 ENERGY STORAGE↗

DOE OSTI2022

A Semi-Empirical Correlation for the Onset of Density Wave Oscillations in a Helical Coil Steam Generator

NuScale Power has completed full length helical coil steam generator (HCSG) tests at SIET Laboratories in Piacenza, Italy. The test program included an investigation of the onset and evolution of density wave oscillations (DWO). This paper describes the mechanisms leading to the onset of DWO as inferred from tests performed in the full-height TF-2 counter-flow test facility. This research supports the concept that the onset of DWO is caused by liquid bridging inside the boiling length of the tubes. Using a separate cylinders approximation, a DWO onset correlation is obtained by performing a steady-state momentum balance on the density front created by liquid bridging. This yields a momentum balance equation in terms of a set of dimensionless superficial velocities, a corresponding dimensionless pressure drop term and an average void fraction. Envelope theory is then used to obtain a critical void fraction which is substituted into the dimensionless momentum balance to obtain a DWO onset correlation similar to the Wallis counter-current flooding correlation. A comparison of the DWO onset correlation to the TF-2 DWO experimental data is presented and shows excellent agreement when the average equilibrium vapor quality is greater than 0.2.

Density Wave Oscillations, Helical Coil Steam Gene↗

DOE OSTI2025

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Real-time dynamics of the Schwinger model as an open quantum system with Neural Density Operators

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Two-Dimensional Moiré Polaronic Electron Crystals

Two-dimensional moiré materials have emerged as the most versatile platform for realizing quantum phases of electrons. Here, we explore the stability origins of correlated states in WSe 2 /WS 2 moiré superlattices. We find that ultrafast electronic excitation leads to partial melting of the Mott states on timescales 5 times longer than predictions from the charge hopping integrals and that the melting rates are thermally activated, with activation energies of 1⁢8±3 and 1⁢3±2 meV for the one- and two-hole Mott states, respectively, suggesting significant electron-phonon coupling. A density functional theory calculation of the one-hole Mott state confirms polaron formation and yields a hole-polaron binding energy of 16 meV. These findings reveal a close interplay of electron-electron and electron-phonon interactions in stabilizing the polaronic Mott insulators at transition metal dichalcogenide moiré interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Quantum Many-Body Theory from a Solution of the N -Representability Problem

Here, in this study, we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H 8 ring.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Microscopic Origin of Polarization Charges at Ga N / ( Al , Ga ) N Interfaces

Ga⁢N/(Al,Ga)⁢N heterojunctions are at the heart of high-electron-mobility transistors that are being adopted for high-power and high-frequency applications. The strong polarization fields present at this interface significantly enhance the density of the two-dimensional electron gas (2DEG) that is confined on the Ga⁢ N side of the junction. The microscopic origin of these electrons has been debated over the years: after excluding that they would be contributed by bulk donors, a model that identifies surface states on the (Al, Ga)⁢N surface as the source of electrons has become widely adopted. Recently, it has become clear, however, that the measured density of surface states is insufficient to account for the high electron density in the 2DEG. Here we demonstrate, based on state-of-the-art first-principles calculations, that the source of electrons is intrinsic to the overall structure and that the negative charge in the 2DEG is balanced by fixed charge on the surface. We perform a rigorous study of polarization, using our recently developed methodology for quantifying polarization fields within the finite-sized systems that can be addressed with density-functional calculations. The results show that the electrons that appear in the 2DEG originate locally at the interface, and that the net charge at the interface is predominantly compensated by fixed charge on the surface, rather than surface states. We elucidate the source of this fixed charge and associate it with surface reconstructions or the presence of heterovalent impurities (such as oxygen). Our results force a reassessment of the impact of surface states on the density of the 2DEG: rather than serving as a supply of electrons, the surface states mainly act to pin the Fermi level. Our conclusions allow a fresh interpretation of experimental observations and allow devising guidelines for optimizing carrier densities in the 2DEG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗