Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic 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 379 records · Page 21

Spectroscopy of the fractal Hofstadter energy spectrum

Hofstadter’s butterfly, the predicted energy spectrum for non-interacting electrons confined to a two-dimensional lattice in a magnetic field, is one of the most remarkable fractal structures in nature. At rational ratios of magnetic flux quanta per lattice unit cell, this spectrum shows self-similar distributions of energy levels that reflect its recursive construction. For most materials, Hofstadter’s butterfly is predicted under experimental conditions that are unachievable using laboratory-scale magnetic fields. More recently, electrical transport studies have provided evidence for Hofstadter’s butterfly in materials engineered to have artificially large lattice constants, such as those with moiré superlattices. Yet, so far, direct spectroscopy of the fractal energy spectrum predicted by Hofstadter nearly 50 years ago has remained out of reach. Here, in this work, we use high-resolution scanning tunnelling microscopy/spectroscopy (STM/STS) to investigate the flat electronic bands in twisted bilayer graphene (TBG) near the predicted second magic angle, an ideal setting for spectroscopic studies of Hofstadter’s spectrum. Our study shows the fractionalization of flat moiré bands into discrete Hofstadter subbands and discerns experimental signatures of self-similarity of this spectrum. Moreover, our measurements uncover a spectrum that evolves dynamically with electron density, showing phenomena beyond that of Hofstadter’s original model owing to the combined effects of strong correlations, Coulomb interactions and the quantum degeneracy of electrons in TBG.

electronic properties and materials↗

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗

Solvent oligomerization pathways facilitated by electrolyte additives during solid-electrolyte interphase formation

The solid–electrolyte interphase (SEI) layer formation is known to play an important role in determining the lifetime of lithium-ion batteries. A thin, stable SEI layer is linked to overall improved battery performance and longevity, however, the factors and mechanisms that lead to optimal SEI morphology and composition are not well understood. Inclusion of electrolyte additives (fluoroethylene carbonate, FEC; and vinylene carbonate, VC) is often necessary for improving SEI characteristics. To understand how these electrolyte additives impact SEI formation, we employed molecular dynamics (MD) and density functional theory (DFT) simulations to study the reaction networks and oligomerization pathways, respectively, for three systems containing ethylene carbonate (EC), a lithium ion, and FEC or VC. MD simulations suggest radical oligomerization pathways analogous to traditional oligomerization with nucleophilic alkoxide species via SN1 reaction mechanisms. Both S N 1 and S N 2 mechanisms were studied for all three systems using DFT. Oligomerization reactions were studied with both a standard alkoxide species and a ring-opened EC radical as the nucleophiles and EC, FEC, and VC as the electrophiles. For all cases, FEC and VC exhibited lower free energy barriers and more stable adducts when compared with EC. We conclude that one of the role of additives is to modify the oligomerization process of EC by introducing branching points (FEC) or termination points (VC).

25 ENERGY STORAGE↗

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗

New interaction potentials for borate glasses with mixed network formers

Here, we adapt and apply a recently developed optimization scheme used to obtain effective potentials for aluminosilicate glasses to include the network former boron into the interaction parameter set. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by ab initio molecular dynamics simulations, and density, coordination, and elastic modulus of glass at room temperature from experiments. The new interaction potentials are shown to reliably reproduce the structure, coordination, and mechanical properties over a wide range of compositions for binary alkali borates. Furthermore, the transferability of these new interaction parameters allows mixing to reliably reproduce the properties of various boroaluminate and borosilicate glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A windowed mean trajectory approximation for condensed phase dynamics

We propose a trajectory-based quasi-classical method for approximating dynamics in condensed phase systems. Building upon the previously developed optimized mean trajectory approximation that has been used to compute linear and nonlinear spectra, we borrow some ideas from filtering trajectory methods to obtain a novel semiclassical method for the dynamical propagation of density matrices. This new approximation is tested rigorously against standard multistate electronic models, spin-boson models, and models of the Fenna–Matthews–Olson complex. For dissipative systems, the current method is significantly better or as good as many other semiclassical methods available, especially at low temperatures and for off-diagonal density matrix elements, whereas for scattering models, the current method bears similar limitations as mean-field propagation schemes. All results are tested against the numerically exact hierarchical equations of motion method. In conclusion, the new method shows excellent agreement across various parameter regimes with numerically exact results, highlighting the robustness and accuracy of our approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗

High-pressure reduction of carbon dioxide in reactive liquid mixtures

Density-functional theory based molecular-dynamics simulations were used to investigate high-pressure chemical reactions in liquid mixtures of CO 2 with several elements (Si, Mn, and Fe) at high temperatures of 2000-3000 K. Our ab initio simulations indicate that these reactant elements can reduce CO 2 to C at high pressures (20 GPa) leading to the formation of C-C chains, with Si by far the most effective carbon-reducing agent. A combined chemical analysis using Bader charge analysis and Crystal Orbital Hamilton Population (COHP) on simulation snapshots shows that significant charge transfer from the reducing element to the C atoms creates instability in the C-O covalent bonds. COHP analysis further shows that Mn/Fe-O and Mn/Fe-C bonding interactions are weaker compared to the Si counterparts. These results further our understanding of the redox chemistry of CO 2 at conditions relevant to planetary mantle interiors and demonstrate the effectiveness of high pressure in the reduction of CO 2 directly to solid carbon.

Ab initio molecular dynamics↗

Comparative many-body study of Pr 4 Ni 3 O 8 and NdNiO 2

We study the many-body electronic structure of the stoichiometric and electron-doped trilayer nickelate Pr4Ni 3 O 8 in comparison to that of the stoichiometric and hole-doped infinite layer nickelate NdNiO 2 within the framework of density functional plus dynamical mean field theory, noting that Pr4Ni 3 O 8 has the same nominal carrier concentration as NdNiO2 doped to a level of 1/3 holes/Ni. In this work, we find that the correlated Ni-3$\textit{d}$ shells of both of these low valence nickelates have similar many-body configurations with correlations dominated by the $d_{x^2 - y^2}$ orbital. Additionally, when compared at the same nominal carrier concentration, the materials exhibit similar many-body electronic structures, self energies, and correlation strengths, but differ in Fermiology. Compared to cuprates, these materials are closer to the Mott-Hubbard regime due to their larger charge transfer energies. Moreover, doping involves the charge reservoir provided by the rare earth 5$\textit{d}$ electrons, as opposed to cuprates where it is realized via the oxygen 2$\textit{p}$ electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electric fields and substrates dramatically accelerate spin relaxation in graphene

Electrons in graphene are theoretically expected to retain spin states much longer than most materials, making graphene a promising platform for spintronics and quantum information technologies. Here, in this work, we use first-principles density-matrix (FPDM) dynamics simulations to show that interaction with electric fields and substrates strongly enhances spin relaxation through scattering with phonons. Consequently, the relaxation time at room temperature reduces from microseconds in free-standing graphene to nanoseconds in graphene on the hexagonal boron nitride (hBN) substrate, which is the order of magnitude typically measured in experiments. Further, inversion symmetry breaking by hBN introduces a stronger asymmetry in electron and hole spin lifetimes than predicted by the conventional D'yakonov-Perel' (DP) model for spin relaxation. Deviations from the conventional DP model are stronger for in-plane spin relaxation, resulting in out-of-plane to in-plane lifetime ratios much greater than 1/2 with a maximum close to the Dirac point. These FPDM results, independent of symmetry-specific assumptions or material-dependent parameters, also validate recent modifications of the DP model to explain such deviations. Overall, our results indicate that spin-phonon relaxation in the presence of substrates may be more important in graphene than typically assumed, requiring consideration for graphene-based spin technologies at room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Valence, charge transfer, and orbital-dependent correlation in bilayer nickelates Nd 3 Ni 2 O 7

We examine the bulk electronic structure of Nd 3 Ni 2 O 7 using Ni 2 p core-level hard x-ray photoemission spectroscopy combined with density functional theory + dynamical mean-field theory. Our results reveal a large deviation of the Ni 3 d occupation from the formal Ni 2.5 + valency, highlighting the importance of the charge transfer from oxygen ligands. We find that the dominant d 8 configuration is accompanied by nearly equal contributions from d 7 and d 9 states, exhibiting an unusual valence state among Ni-based oxides. Finally, we discuss the Ni d x 2 − y 2 and d z 2 orbital-dependent hybridization, correlation and local spin dynamics. Published by the American Physical Society 2025

Takegami, Daisuke (ORCID:0000000338321849)↗

Incommensurate Spin Fluctuations in the Spin-Triplet Superconductor Candidate UTe 2

Spin-triplet superconductors are of extensive current interest because they can host topological state and Majorana fermions important for quantum computation. The uranium-based heavy-fermion superconductor UTe 2 has been argued as a spin-triplet superconductor similar to UGe2, URhGe, and UCoGe, where the superconducting phase is near (or coexists with) a ferromagnetic (FM) instability and spin-triplet electron pairing is driven by FM spin fluctuations. Here we use neutron scattering to show that, although UTe 2 exhibits no static magnetic order down to 0.3 K, its magnetism in the [0,K,L] plane is dominated by incommensurate spin fluctuations near an antiferromagnetic ordering wave vector and extends to at least 2.6 meV. We are able to understand the dominant incommensurate spin fluctuations of UTe 2 in terms of its electronic structure calculated using a combined density-functional and dynamic mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Understanding the lattice thermal conductivity of SrTiO 3 from an ab initio perspective

Here, we present a detailed analysis of the structure dependence of the lattice thermal conductivity of SrTiO 3 . We have used both ab initio molecular dynamic simulations and density functional theory calculations to decouple the effect of different structural distortions on the thermal conductivity. We have identified two main mechanisms for tuning the thermal conductivity when a distortion is applied. First, the modification of the acoustic-modes energy dispersion when a change in the lattice parameters is imposed and second, the low energy polar modes. In particular and counterintuitively, we have found that an increase in the angle of the oxygen octahedral rotations increases the thermal conductivity due to its coupling to these polar modes.

36 MATERIALS SCIENCE↗

Analytically Derived Piezoelectric Ejecta Data Transformations, with Implications for Claims of Experimentally Observed Self-Similarity

Building upon a previously documented analytic formalism, we derive a set of data transformations that relate piezoelectrically inferred time-dependent ejecta masses, dynamical ejecta cloud densities, and "radiographic" ejecta cloud density distributions. These transformation rules are derived rigorously, with careful attention to the assumption of instant ejecta production which is inherent in the piezoelectric mass inference procedure. Some of these relationships have not appeared previously in the literature, including solutions for the unique ejecta velocity distribution and radiographic ejecta cloud density distribution consistent with the instant-production assumption for a given piezoelectrically inferred mass dataset. In deriving the latter, we prove that—again subject to the ubiquitous instant-creation assumption—measurements at a single sensor are sufficient to determine properties of the ejecta cloud everywhere at all times. Other data transformation relationships have appeared in the literature without a rigorous derivation, or have been stated in forms that obscure the key mathematical relationships. We apply these results to arguments for self-similar ejecta cloud evolution, thereby deriving mathematical expressions by which claims of self-similarity may be examined. We prove that one analysis used in the literature to claim experimental evidence of a short ejecta production timescale—and consequently of self-similar cloud evolution—in fact places no constraint the ejecta production interval and therefore does not support the associated conclusion. We also provide a potentially useful uniqueness proof.

30 DIRECT ENERGY CONVERSION↗

Briefing Book: Office of Experimental Sciences

Key to OES science is the advancement of experimental capabilities and vibrant partnerships with other elements of NNSA’s Defense Programs. OES fosters an experimental portfolio that spans the space from small scale studies utilizing tabletop instruments through complex and dynamic high energy density; hydrodynamic; and subcritical experiments on flagship national facilities like Los Alamos Neutron Science Center (LANSCE), National Ignition Facility (NIF), Z machine at Sandia National Laboratories, Dual Axis Radiographic Hydrodynamic Test (DARHT), and the U1a Complex in Nevada. Essential to meeting the deterrence challenges that drive OES is a functional strategy for integration of OES with Defense Programs pursuits in modeling and simulation, engineering and technology maturation, and stockpile and production modernization.

42 ENGINEERING↗

Detection and Association of Operational Events using DAS and Seismometers (FY 2025 Mid-Year Report)

This mid-year report summarizes ongoing work to identify anomalous vibration signals indicative of potential containment breaches. This work includes compiling continuous seismic datasets and testing and refining underground detection and geolocation techniques. In the first two quarters of FY25, we have completed two project work plan tasks: (1) creating a database of continuous waveforms and ground truth event data from multiple modalities and (2) refining and implementing a detection and association algorithm to create a catalog of anomalous underground activities. This report contains a summary of the seismic database including the continuous seismic data collected by a dense array of surface seismic stations above Pleasant Gap Mine, and continuous seismic data collected using subsurface distributed acoustic sensing (DAS) in the subsurface at Sanford Underground Research Facility (SURF) and the ground truth information gathered from both sites. This report also includes results from refining and applying a dynamic power spectral density detector to both continuous seismic datasets. Finally, the report provides an initial catalog of subsurface operational events from both sensing modalities.

58 GEOSCIENCES↗

The GREEN ‘omics of Nutrient Feedbacks to Soil Warming

The GREEN ‘omics of Nutrient Feedbacks in Soil project advanced the DOE Biological and Environmental Research (BER) mission by developing and applying isotope-enabled ’omics tools to understand how soil microbes regulate carbon and nutrient cycling. Guided by the Growth Rate, growth Efficiency, and stoichiometry of Essential Nutrients (GREEN ’omics) framework, the project aimed to build a predictive, systems-level understanding of microbial traits that control ecosystem biogeochemistry. In a collaboration among Northern Arizona University (lead), West Virginia University, Lawrence Livermore National Laboratory, and Pacific Northwest National Laboratory, we combined quantitative stable isotope probing (qSIP), Chip-SIP, NanoSIMS, and genome-resolved metagenomics across long-term experiments in Arctic, boreal, temperate, and tropical ecosystems. The project produced three key outcomes: 1) We showed that community-weighted temperature sensitivities of bacterial growth (Q10) can predict ecosystem-scale soil respiration responses across diverse soils. 2) We provided the first in situ evidence for density-dependent population dynamics in soil bacteria and demonstrated that nutrient additions intensify competition, concentrating carbon use into fewer taxa. 3) We improved and extended isotope-enabled ’omics methods by quantifying qSIP measurement error to guide experimental design and coupling SIP with genome-resolved metagenomics to reveal cross-kingdom interactions among bacteria, fungi, and viruses. Together, these results show that a small number of microbial traits and taxa exert disproportionate control over soil carbon and nutrient cycling, providing critical data and methods to improve representation of microbial processes in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

FY 2026 Midyear Report: Seismic Monitoring of Underground Vibration Sources Using Distributed Acoustic Sensing and Seismometers

Safeguards-relevant temporal changes in underground facilities can be observed using geophysical monitoring techniques. Seismic waves, in particular, provide valuable insights into subsurface activities and can serve as an important tool for detecting anomalous events that may indicate containment breaches at geological repositories. This midyear report summarizes ongoing efforts to automatically and rapidly detect and locate anomalous vibration signals that could be indicative of potential containment breaches. Previous work during FY25 focused on compiling continuous seismic datasets from two underground sites and developing a database of continuous waveforms and ground-truth event data derived from multiple sensing modalities. Building on this foundation, we are adapting anomaly detection and geolocation algorithms to explore methods for monitoring underground activities using two relatively low-maintenance sensing technologies: a dense surface geophone array deployed at the Pleasant Gap mine in Pennsylvania, and a three-dimensional fiber-optic cable array for distributed acoustic sensing (DAS) installed in the subsurface at the Sanford Underground Research Facility (SURF) in South Dakota. This report summarizes work conducted during the first two quarters of FY26, during which we refined a dynamic power spectral density (PSD)-based detector, applied it independently to each geophone station, and then combined the per‑station detections with density-based spatial clustering of applications with noise (DBSCAN) to cluster events and produce spatial maps over a nine‑day interval. In addition, we outline plans for a field trial at the Waste Isolation Pilot Plant (WIPP) in New Mexico to compare traditional seismic monitoring approaches with DAS techniques and to evaluate the benefits of combined data analysis. Activities during the past two quarters have included the preparation and submission of a Field Test Plan to WIPP for approval, as well as submission to headquarters for review and feedback.

58 GEOSCIENCES↗