Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Position-dependent Voronoi probability distribution functions for matter and halos

Here, we measure the Voronoi density probability distribution function (PDF) for both dark matter and halos in N-body simulations. For the dark matter, Voronoi densities represent the matter density field smoothed on a uniform mass scale, which approximates the Lagrangian density field. For halos, the Voronoi densities contain information about the local environment of each halo. We measure the halo virial masses, the total amount of dark matter within each halo Voronoi cell, and the halo Voronoi cell volumes, and we show how halo abundances depend on these three quantities. We then study the position-dependent Voronoi density PDF, measured within finite subregions of the Universe, using separate universe simulations. We demonstrate that the spatial variation of the position-dependent PDF is due to large-scale density fluctuations, indicating that the position-dependent PDF is a biased tracer of large-scale structure. We measure this bias for the dark matter, and interpret it as the bias of regions of the Lagrangian density field that are selected based on density. For the halos, this bias can be interpreted as a form of assembly bias. We present the mapping from late-time to early-time Voronoi density for each simulation dark matter particle, which is highly stochastic. We compare the median of this stochastic map with spherical collapse calculations and discuss challenges involved in modeling the evolution of the density field on these scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Kohn–Sham Density in a Slater Orbital Basis Set

Finite, atom-centered Slater basis sets are used to determine approximate Kohn–Sham molecular orbitals. This is achieved by minimizing the kinetic energy plus the sum-squared difference between the Kohn–Sham density and the full configuration interaction density. As a result of the finite basis, a weight factor is introduced to balance the two minimization components. Results herein show that this can be done systematically, without sensitive dependence on the choice of scaling factor. In addition, the algorithm is applied to the LiH diatomic for fractional electron counts, where stretching the bond introduces significant reorganization of the electron density. As a result, the analysis will show the correct KS orbital structure and reveal the effects of correlation and electron locality on the KS solutions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Here, this study overviews and extends a recently developed stochastic finite-temperature Kohn-Sham density functional theory to study warm dense matter using Langevin dynamics, specifically under periodic boundary conditions. The method's algorithmic complexity exhibits nearly linear scaling with system size and is inversely proportional to the temperature. Additionally, a linear-scaling stochastic approach is introduced to assess the Kubo-Greenwood conductivity, demonstrating exceptional stability for dc conductivity. Utilizing the developed tools, we investigate the equation of state, radial distribution, and electronic conductivity of hydrogen at a temperature of 30 000 K. As for the radial distribution functions, we reveal a transition of hydrogen from gaslike to liquidlike behavior as its density exceeds 4 g/cm 3 . As for the electronic conductivity as a function of the density, we identified a remarkable isosbestic point at frequencies around 7 eV, which may be an additional signature of a gas-liquid transition in hydrogen at 30 000 K.

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

Reaction-diffusive dynamics of number-conserving dissipative quantum state preparation

The use of dissipation for the controlled creation of nontrivial quantum many-body correlated states is of much fundamental and practical interest. What is the result of imposing number conservation, which, in closed system, gives rise to diffusive spreading? For this work, we investigate this question for a paradigmatic model of a two-band system, with dissipative dynamics aiming to empty one band and to populate the other, which had been introduced before for the dissipative stabilization of topological states. Going beyond the mean-field treatment of the dissipative dynamics, we demonstrate the emergence of a diffusive regime for the particle and hole density modes at intermediate length- and timescales, which, interestingly, can only be excited in nonlinear response to external fields. We also identify processes that limit the diffusive behavior of this mode at the longest length- and timescales. Strikingly, we find that these processes lead to a reaction-diffusion dynamics governed by the Fisher-Kolmogorov-Petrovsky-Piskunov equation, making the designed dark state unstable towards a state with a finite particle and hole density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Cavity electrodynamics of van der Waals heterostructures

Van der Waals heterostructures host many-body quantum phenomena that are tunable in situ using electrostatic gates. Their constituent two-dimensional materials and gates can naturally form plasmonic self-cavities, confining light in standing waves of current density due to finite-size effects. The plasmonic resonances of typical graphite gates fall in the gigahertz to terahertz range, corresponding to the same microelectronvolt to millielectronvolt energy scale as the phenomena in van der Waals heterostructures that they electrically control. This raises the possibility that the built-in cavity modes of graphite gates are relevant for shaping the low-energy physics of these heterostructures. However, probing these cavity-coupled electrodynamics is challenging as devices are notably smaller than the diffraction limit at the relevant wavelengths. Here we report on the intrinsic cavity conductivity of gate-tunable graphene heterostructures. As the carrier density is tuned, we observe coupling and spectral weight transfer between graphene and graphite plasmonic cavity modes in the ultrastrong coupling regime. We present an analytical model to describe the results and provide general principles for cavity design. Our findings show that intrinsic cavity effects are important for understanding the low-energy electrodynamics of van der Waals heterostructures and open a pathway for useful functionality through cavity control.

Electronic properties and materials↗
DOE OSTI2022

Bioinspired mechanically interlocking holey graphene@SiO 2 anode

Mechanically interlocking structures that can enhance adhesion at the interface and regulate the stress distribution have been widely observed in biological systems. Inspired by the biological structures in the wings of beetles, we synthesized a holey graphene@SiO 2 anode with strong mechanical interlocking, characterized it electrochemically, and explained its performance by finite element analysis and density functional calculations. The mechanically interlocking structure enhances lithium-ion (Li + ) storage by transmitting the strain from SiO 2 to the holey graphene and by a mechano-electrochemical coupling effect. The interlocking fit hinders the abscission of SiO 2 and the distinctive structure reduces the stress and strain of SiO 2 during (de)lithiation. The positive mechano-electrochemical coupling effect preserves the amount of electrochemically active phase (Li x Si) during cycles and facilitates Li + diffusion. Therefore, the capacity shows only a slight attenuation after 8000 cycles (cycling stability), and the specific capacity is ~1200 mA h g –1 at 5 A/g (rate-performance). This study furnishes a novel way to design high-performance Li + /Na + /K + /Al 3+ anodes with large volume expansion.

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

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

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

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

Electronic rotons and Wigner crystallites in a two-dimensional dipole liquid

A key concept proposed by Landau to explain superfluid liquid helium is the elementary excitation of quantum particles called rotons1–8. The irregular arrangement of atoms in a liquid leads to the aperiodic dispersion of rotons, which played a pivotal role in understanding fractional quantum Hall liquids (magneto-rotons)9,10 and the supersolidity of Bose–Einstein condensates11–13. Even for a two-dimensional electron or dipole liquid, in the absence of a magnetic field, the repulsive interactions have been predicted to form a roton minimum14–19, which can be used to trace the transition to Wigner crystals20–24 and superconductivity25–27, although this has not yet been observed. Here, we report the observation of such electronic rotons in a two-dimensional dipole liquid of alkali-metal ions donating electrons to surface layers of black phosphorus. Our data reveal the striking aperiodic dispersion of rotons, which is characterized by a local minimum of energy at finite momentum. As the density of dipoles decreases so that interactions dominate over the kinetic energy, the roton gap reduces to 0, as in a crystal, signalling Wigner crystallization. Our model shows the importance of short-range order arising from repulsion between dipoles, which can be viewed as the formation of Wigner crystallites (bubbles or stripes) floating in the sea of a Fermi liquid. Our results reveal that the primary origin of electronic rotons (and the pseudogap) is strong correlations.

Park, Soobin↗
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2025

New state of matter between the hadronic phase and the quark-gluon plasma?

Lattice-quantum chromodynamics (QCD) simulations and theoretical arguments hint at the existence of an intermediate phase of strongly interacting matter between a confined hadron gas and a deconfined quark-gluon plasma (QGP). We qualitatively and semiquantitatively explore and differentiate the phase structures in the temperature window from the QCD pseudocritical temperature 𝑇 c ≃ 160 MeV to the pure gluonic deconfinement temperature 𝑇 d ≃ 285 MeV. We propose a three-regime picture using a hadron resonance gas description augmented with the exponential spectrum of strings, corresponding to highly excited mesons and glueballs, based on the analysis of a large number 𝑁 c of colors. We estimate the entropy density from our model to confirm that the lattice-QCD data are bracketed with three regimes, i.e., a hadron gas, a QGP, and a new phase for 𝑇 c ≲ 𝑇 ≲ 𝑇 d . In this new phase, which we name a spaghetti of quarks with glueballs (SQGBs), thermal degrees of freedom of quarks are liberated, yet gluons remain confined in glueballs. Since the Hagedorn temperature 𝑇 H ∼ 285 MeV is universal in the meson and the glueball sectors, in the infinite-𝑁 c limit, the phase diagram in the plane of the baryon chemical potential and the temperature is reduced to one with the confined and deconfined phases and quarkyonic matter at high density. At large but finite 𝑁 c , an SQGB window may open between these phases. We point out that the SQGB has interesting similarities with quarkyonic matter and that this matter in the large-𝑁 c limit is confined as measured by the interaction between heavy quarks, but behaves in other respects like a quasifree gas of quarks. As a result of the extrapolation to 𝑁 c = 3, we present a revised phase diagram with the SQGB phase bounded by thermal crossovers. Finally, we give a quantitative analysis of chiral-symmetry restoration in the SQGB phase.

Color confinement↗
DOE OSTI2023

High-temperature self-energy corrections to x-ray absorption spectra

Effects of finite-temperature quasiparticle self-energy corrections on x-ray absorption spectra are investigated within the finite-temperature quasiparticle local density GW approximation up to temperatures T of order the Fermi temperature. To facilitate the calculations, we parametrize the quasiparticle self-energy using low-order polynomial fits. Here we show that temperature-driven decrease in the electron lifetime substantially broadens the spectra in the near-edge region with increasing T. However, the quasiparticle shift is most strongly modified near the onset of plasmon excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Vorticity screening by dense layers

We investigate variable density layers in the presence of vortices, by analogy to line charges in the presence of a dielectric slab. We adapt solutions for dielectric layers from the literature to the variable density fluid case, obtaining exact hypergeometric expressions for flow velocities induced everywhere by a vortex in the vicinity of a finite-thickness slab of different densities. The dimensionless Atwood number, which appears naturally in other variable density phenomena, such as the Rayleigh–Taylor instability, reappears here as the natural ordering parameter for the influence of the slab, which acts primarily to screen the influence of vorticity across the slab to a lower effective circulation. The solution takes the form of scaling correction factors that might be applied to computational models or evaluated directly. Extensions, where vorticity is localized within the dense layer itself and amplified outside of it, are considered, as well as some speculative applications such as inertial confinement fusion capsule designs.

42 ENGINEERING↗
DOE OSTI2023

Superconducting valence bond fluid in lightly doped eight-leg $t-J$ cylinders

Superconductivity in doped quantum paramagnets has been a subject of long theoretical inquiry. In this work, we report a density matrix renormalization group study of lightly doped $t-J$ models on the finite-width square lattice (doped hole densities $δ$ = 1/12 and 1/8) with parameters for which previous studies have suggested that the undoped system in 2D is either a quantum spin liquid or a valence bond crystal. Our studies are performed on cylinders with width up to 8. Here, ground-state correlations are found to be nearly identical for the “doped quantum spin liquid” and “doped valence bond crystal.” Upon increasing the cylinder widths from 4 to 8, we observed a significant strengthening of the quasi-long-range superconducting correlations and a dramatic suppression of any “competing” charge density wave order. Extrapolating from the observed behavior of the width eight cylinders, we speculate that the system has a nodeless $d$-wave superconducting ground state in the 2D limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Eulerian finite element simulations of the drop weight impact test with a dislocation Density-based continuum model

During a drop weight impact test, the kinetic energy of the falling weight is transferred into the sample resting on the anvil. The plastic deformation in the sample is an important mechanism for the dissipation of the input kinetic energy. We use Eulerian finite element analysis to simulate the deformation and temperature evolution in a 1,3,5-trinitro-1,3,5-triazine (RDX) sample consisting of multiple crystals. In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, i.e., advection. In an effort to reduce the advection errors, we use a rate form of a dislocation density-based continuum model by Luscher et al. Here, the simulations predict localization of plastic deformation, and plastic dissipation as a significant source of heat generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Reassessment of Birch's Law on hcp‐Fe From Ultrasonic Sound Velocity Measurement and Implications on the Velocity Profiles of Earth's Inner Core

Here, we performed in situ X-ray diffraction and ultrasonic sound velocity measurements on hcp-Fe up to 15 GPa, 873 K in a multi-anvil apparatus. The elastic moduli and their pressure and temperature derivatives were determined by fitting the velocity and density data to the third-order finite strain equations, yielding K S0 = 169.0(57) GPa, K$^{′}_{S0}$ = 5.4(6), (∂K S /∂T) P = −0.031(3) GPa/K, G 0 = 104.5(27) GPa, G$^{′}_{0}$ = 1.7(2), (∂G/∂T) P = −0.060(2) GPa/K. Within the experimental P-T range, we find significant temperature effect on the density-velocity (ρ − V P or S ) relations and caution the use of the temperature-independent Birch's law for extrapolation to core conditions. Furthermore, temperature-induced velocity decrease is more significant in V S than in V P , offering a possible explanation for the high Poisson's ratio in the core. Extrapolations based on our results together with previous experimental data suggest that V P of hcp-Fe aligns with PREM at Earth's core conditions, while Vs and density are approximately 10% and 2.7% higher than PREM, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multi-scale transport in the DIII-D ITER baseline scenario with direct electron heating and projection to ITER

Multi-scale fluctuations measured by turbulence diagnostics spanning long and short wavelength spatial scales impact energy confinement and the scale-lengths of plasma kinetic profiles in the DIII-D ITER baseline scenario with direct electron heating. Contrasting discharge phases with ECH + neutral beam injection (NBI) and NBI only at similar rotation reveal higher energy confinement and lower fluctuations when only NBI heating is used. Modeling of the core transport with TGYRO using the TGLF turbulent transport model and NEO neoclassical transport reproduces the experimental profile changes upon application of direct electron heating and indicates that multi-scale transport mechanisms are responsible for changes in the temperature and density profiles. Intermediate and high-k fluctuations appear responsible for the enhanced electron thermal flux, and intermediate-k electron modes produce an inward particle pinch that increases the inverse density scale length. Projection to ITER is performed with TGLF and indicates a density profile that has a finite scale length due to intermediate-k electron modes at low collisionality and increases the fusion gain. For a range of E×B shear, the dominant mechanism that increases fusion performance is suppression of outward low-k particle flux and increased density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗