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 523 records · Page 29

Noninvasive measurements of spin transport properties of an antiferromagnetic insulator

Antiferromagnetic insulators (AFIs) are of substantial interest because of their potential in the development of next-generation spintronic devices. One major effort in this emerging field is to harness AFIs for long-range spin information communication and storage. Here, we report a noninvasive method to optically access the intrinsic spin transport properties of an archetypical AFI α-Fe 2 O 3 via nitrogen-vacancy (NV) quantum spin sensors. By NV relaxometry measurements, we successfully detect the frequency-dependent dynamic fluctuations of the spin density of α-Fe 2 O 3 along the Néel order parameter, from which an intrinsic spin diffusion constant of α-Fe 2 O 3 is experimentally measured in the absence of external spin biases. Our results highlight the significant opportunity offered by NV centers in diagnosing the underlying spin transport properties in a broad range of high-frequency magnetic materials such as two-dimensional magnets, spin liquids, and magnetic Weyl semimetals, which are challenging to access by the conventional measurement techniques.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reversible phase transformations between Pb nanocrystals and a viscous liquid-like phase

Phase transformations have been a prominent topic of study for both fundamental and applied science. Solid-liquid reaction–induced phase transformations can be hard to characterize, and the transformation mechanisms are often not fully understood. Here, we report reversible phase transformations between a metal (Pb) nanocrystal and a viscous liquid-like phase unveiled by in situ liquid cell transmission electron microscopy. The reversible phase transformations are obtained by modulating the electron current density (between 1000 and 3000 electrons Å –2 s –1 ). The metal-organic viscous liquid-like phase exhibits short-range ordering with a preferred Pb-Pb distance of 0.5 nm. Assisted by density functional theory and molecular dynamics calculations, we show that the viscous liquid-like phase results from the reactions of Pb with the CH 3 O fragments from the triethylene glycol solution under electron beam irradiation. Such reversible phase transformations may find broad implementations.

36 MATERIALS SCIENCE↗

Di-$\pi^0$ correlations in p+p, p+Al and p+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at STAR

The STAR Collaboration reports measurements of back-to-back azimuthal correlations of di- \pi^0 π 0 produced at forward pseudorapidity ( 2.6<\eta<4.0 2.6 < η < 4.0 ) in p+ p + p p , p+ p + \mathrm{Al} A l and p+ p + \mathrm{Au} A u collisions at a center-of-mass energy per nucleon-nucleon pair of 200 GeV. A clear suppression of the correlated away-side yields is observed in p+ p + \mathrm{Au} A u and for the first time in p+ p + \mathrm{Al} A l collisions, compared with the p p + + p p data. The enhanced suppression found in p+ p + \mathrm{Au} A u with respect to p+ p + \mathrm{Al} A l collisions exhibits the saturation scale ( Q_s^2 Q s 2 ) dependence on the mass number A A . The observed suppression of back-to-back pairs as a function of event activity and transverse momentum points to non-linear gluon dynamics arising at high parton densities.

Chu, Xiaoxuan↗

Interpenetrating Plasma Simulations

This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY20 Proxy App Suite Release: Report for ECP Proxy App Project Milestone ADCD-504-10

Version 4.0 of the ECP Proxy App Suite is practically unchanged from the previous release. The current set of proxies has proven useful for many aspects of benchmarking and co-design and we see little reason to alter the suite. Although there have been few changes to the ECP suite, the team has been hard at work in other areas. In the area of Machine Learning (ML) we have now created a separate proxy suite dedicated to this scientific applications of ML. The suite includes: miniGAN (Generative Adversarial Networks), miniRL (Reinforcement Learning), CRADL (inline inference), Cosmoflow-Benchmark (Convolutional Neural Network), and MLPerf-DeepCam (Climate Segmentation Benchmark). Section 3 contains more information about these proxies as well as the principles that are guiding the development of the suite. We have surveyed available proxies for several application domains including Computational Fluid Dynamics, Quantum Chemistry, Quantum Computing Simulation, Molecular Dynamics, Monte Carlo Transport, and Density Functional Theory to identify gaps in proxy coverage. Several new proxy apps are either already available or will be released soon to fill these gaps. Section 4 provides full details. Finally, section 5 reports on our continued collaboration with the ECP Continuous Integration (CI) effort to use proxy apps to help identify problems and roadblocks to cross-lab CI. We have also assisted the El Capitan Center of Excellence (COE) to stand up a CI system that will be used to test releases of HPE and AMD software stacks. We hope that the COE effort can serve as a model for ECP by showing how the Proxy App Team can work with AD and ST teams identify critical features, kernels, patterns, etc. and incorporate them into a CI system that will help ensure that Frontier and Aurora will provide those needed capabilities.

97 MATHEMATICS AND COMPUTING↗

Fundamental Study of Fatigue Crack Initiation at Grain and Twin Boundaries in Austentic Stainless Steel (Final Technical Report)

The overarching goal of this project was to understand the dislocation-based mechanisms driving fatigue crack initiation at grain and twin boundaries. The hypothesis of the work is that dislocation structures equivalent to persistent slip bands (PSBs) form in the boundaries, leading to local stress concentrations and fatigue crack initiation. The sub-goals associated with this work include: • Understanding the evolution of the dynamic stress state and dislocation density surrounding PSB/grain boundary interactions. This will require the development of new techniques that couple high resolution electron backscatter diffraction (HREBSD) with in situ scanning electron microscopy (SEM) deformation. • Uncover the mechanisms by which dislocations are accommodated in grain and twin boundaries during fatigue. The proposed approach here is to conduct in situ transmission electron microscopy (TEM) deformation experiments on fatigued samples. • Establish correlations between microstructure, deformation accommodation, and crack initiation across time and length scales. This work primarily focuses on understanding fatigue crack initiation in austenitic stainless steels, though promising work by an undergraduate researcher extended the work to high purity Al.

42 ENGINEERING↗

Development of scalable design optimization parameters for bi-component protective systems

Bi-component protection-systems need to withstand a variety of extreme loading conditions such as high-velocity projectile impact to shock and blast mitigation. The large number of parameters at each stage of bi-component protection system design process results in a near-infinite number of possible variables for optimization. The aim of this project is to identify the scaling laws linking material and geometric scales of the problem to the final impact performance, further including the effects of wave propagation/interaction and inter-component bonding. Specifically, the impact dynamics for different optimal projectile/target density ratios and relative component mass fractions are quantified for bonded ceramic/titanium alloy systems. The impact behavior such as wave mechanics and damage propagation within the target system is characterized and quantified via digital image correlation or via x-ray phase-contrast imaging.

36 MATERIALS SCIENCE↗

MaNGA DynPop – I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles, and dark matter

ABSTRACT This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here, we present a catalogue of dynamically determined quantities for over 10 000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically aligned JAMcyl or spherically aligned JAMsph. We also make four assumptions for the models’ dark versus luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results.

Astronomy & Astrophysics↗

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗

Combined, time-resolved, in situ neutron reflectometry and X-ray diffraction analysis of dynamic SEI formation during electrochemical N 2 reduction

One means of improving performance for electrochemical ammonia production through the Li-mediated N 2 reduction reaction (Li-NRR) is by cycling the current driving the reaction between open-circuit conditions and periods of applied current density. Herein, we have investigated the dynamics of the electrode–electrolyte interface under Li-NRR conditions during current cycling using in situ time-resolved neutron reflectometry and grazing-incidence synchrotron X-ray diffraction. During cycling, measured neutron reflectivity curves indicated bilayer formation in which Li-containing species such as LiOH, Li 2 O, and small quantities of Li 3 N and metallic Li primarily appeared in a thin layer at the cathode surface, above which formed a much larger, porous, ‘solid–electrolyte interface’ (SEI) layer. Upon return to open-circuit conditions, Li-containing species quickly moved out of the thin layer, leaving a compact, stable layer of decomposition products underneath the SEI layer. This SEI layer concomitantly filled with electrolyte or dissolved, becoming indistinguishable from the electrolyte via contrast in scattering-length density (SLD). During the second current cycle, Li-containing species again preferentially deposited directly atop the cathode, with the thick SEI-like layer again appearing within a minute. This SEI layer exhibited a lower SLD more quickly than in the first cycle, which might suggest that Li-containing species become distributed within the porous SEI layer. Furthermore, these time-resolved observations of SEI and plated layers during current cycling suggest that benefits associated with return to open-circuit conditions between periods of applied current density may be related to the concomitant loss of Li-containing species from a thin layer at the cathode surface into a porous SEI layer that becomes filled with electrolyte or dissolves.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Trimeron-phonon coupling in magnetite

Using density functional theory, we study the lattice dynamical properties of magnetite (Fe 3 O 4 ) in the high-temperature cubic and low-temperature monoclinic phases. The calculated phonon dispersion curves and density of states are compared with the available experimental data obtained by inelastic neutron, inelastic x-ray, and nuclear inelastic scattering. We find a very good agreement between the theoretical and experimental results for the monoclinic Cc structure revealing the strong coupling between the charge-orbital (trimeron) order and specific phonon modes. For the cubic phase, clear discrepancies arise due to fluctuation effects, which are not included in the calculation method. Despite this shortcoming, we argue that the main spectral features can be understood assuming that the strong trimeron-phonon coupling is extended above the Verwey transition, with lattice dynamics influenced by the short-range order instead of the average cubic structure. Our results indicate the validity of trimerons (and trimeron-phonon coupling) to explain the physics of magnetite much beyond their original formulation.

36 MATERIALS SCIENCE↗

Grid-responsive hydrogen production: Capital utilization and current density vs. efficiency in variable electricity markets

To achieve low-cost hydrogen production from water electrolyzers, grid tied electrolysis may need to operate dynamically to minimize the cost of supplying energy to the electrolyzer stack and produce hydrogen during low-cost hours and turn off/down during high-cost hours. Operating systems in this way can decrease capital utilization (capacity factor) and electricity costs. This strategy would shift the dominant cost drivers away from electricity (and thus efficiency) to the capital costs of the system, due to the underutilized capital when operating at low-capacity factors. Increasing the operational current density of the system could, in effect, reduce the capital cost of the system while producing hydrogen at a lower efficiency on a per unit energy basis. In the variable electricity cost profiles analyzed in this paper, increasing the current density for liquid alkaline from 0.5 A/cm2 to 1.5 Ac/m2 and proton exchange membrane electrolyzers from 2 A/cm2 to 4 A/cm2 resulted in substantial reductions in the levelized cost of hydrogen. Additionally, as capacity factors and electricity costs decrease, the optimal operating current density of the electrolyzer systems analyzed increases. These findings suggest R&D efforts should focus on increasing the operational current densities, reducing the turn down ratios, and understanding the durability implications of those strategies on low-temperature liquid alkaline and proton exchange membrane electrolyzers.

08 HYDROGEN↗

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗