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

On the Properties of Cosmological Ionization Fronts

We investigate the properties of cosmological ionization fronts during the epoch of reionization using the Cosmic Reionization on Computers simulations. By analyzing reionization timing maps, we characterize ionization front velocities and curvatures and their dependence on the density structure of the intergalactic medium (IGM). The velocity distribution of ionization fronts in the simulations indicates that while the barrier-crossing analytical model captures the overall shape in high-velocity regions, it fails to reproduce the low-velocity tail, highlighting the non-Gaussian nature of the IGM’s density field. Ionization front velocities are inversely correlated with local density, propagating faster in underdense regions and more slowly in overdense environments. Faster ionization fronts also lead to higher postionization temperatures, reaching a plateau at ∼2 × 10$^{4}$ K for velocities exceeding 3000 km s$^{−1}$. Examining curvature statistics further establishes a connection between ionization front structure and the normalized density contrast ν, with trends in overdense regions aligning well with barrier-crossing model predictions, while deviations appear in underdense environments due to model limitations. These results provide a detailed characterization of ionization front dynamics and their interaction with the underlying density field, bridging small-scale reionization physics with large-scale observables such as the 21 cm signal and the IGM’s thermal history.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Coupled Experimental/Computational Investigation of the Dynamics of Interacting Magnetized Plasmas

The interaction, or interpenetration, of magnetized plasmas of different density and/or pressure occurs in a wide variety of natural and man-made systems. Such systems include extragalactic jets propagating into the intergalactic medium, solar coronal mass ejections into background solar wind, compact toroid (CT) fueling of magnetic fusion plasmas, and jets of capsule shell impurities into DT fusion fuel, which can lead to enhanced impurity mix in inertial fusion implosions. These plasmas may take the form of jets, with open, helical magnetic structure, or plasma “bubbles” with closed magnetic fields (B-fields), such as spheromaks or CT’s. Such structures, both open and closed B-field cases, can transport heat, particles and magnetic flux or magnetic helicity into background plasma regions. For example, the origin of extragalactic magnetic fields may be due, at least in part, to transport by astrophysical jets. The goal of this proposed work was to elucidate the detailed plasma and magnetic field dynamics of high-density plasma jets (open B-field) and bubbles (closed B-field) propagating into lower density background magnetized plasma through controlled laboratory experiments and closely coupled nonlinear MHD modeling. These experiments were conducted in the HelCat (Helicon-Cathode) linear plasma device at the University of New Mexico (UNM). Plasma jets and bubbles were launched via an existing compact coaxial plasma gun, mounted on the HelCat device. This gun produced plasmas tens of cm in scale and lasting tens of microseconds, thereby allowing detailed multipoint, space- and time-resolved measurements to be made routinely. The experiments were directly modeled using the extended magnetohydrodynamic (XMHD) PERSEUS code, developed at Cornell University [23,24]. Both experimental and numerical modeling work are ongoing. The main results to date are reported here. Additional supplemental funding for one year (8/1/2019 – 7/31/2020) supported numerical investigation of photoionization processes important in many low temperature plasmas, including the HelCat device. Initial results of this modeling work is also reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Impact of Cr and Co on 99Tc retention in magnetite: A combined study of ab initio molecular dynamics and experiments

This work explores the effect of co-mingled dopants, Co(II) and Cr(III), on Tc(IV) incorporation and retention in magnetite when heat treated to 625 or 700 °C. Key trends in Tc retention in the high temperature regime were identified using a combination of density-functional-theory based ab initio molecular dynamics (AIMD) simulations, and batch experiments including solid phase characterization techniques, e.g. X-ray absorption spectroscopy. A stabilizing effect on Tc(IV) was observed when the number of Tc and Cr atoms are equal or when the magnetite surface is oversaturated with Tc and Cr inclusions. Here, oversaturation is hypothesized to force Cr from the magnetite surface to form a Cr2O3 phase, which may act as a protective layer that prevents Tc release. With the addition of Co, Tc(IV) is stabilized via redox processes. The presence of Cr in low concentrations interferes with this redox stabilization and Cr is preferentially stabilized as opposed to Tc. As a result, using Co as a stabilizing dopant for Tc may be compromised in the presence of Cr. When the relative concentration of Tc, Cr and Co is the same, or more Co atoms are added to high Cr incorporated systems, the formation of the Cr2O3 phase may be suppressed. Although waste streams with co-mingled Tc and Cr potentially may benefit from a Co dopant, since the formation of a Cr2O3 passivation layer may protect incorporated Tc from being released, the relative concentration of the three elements will be a critical parameter for maximizing effectiveness of this strategy.

Lee, Mal Soon↗
DOE OSTI2022

Tuning Photoexcited Charge Transfer in Imine-Linked Two-Dimensional Covalent Organic Frameworks

The generation of a long-lived charge-separated state in versatile π-conjugated two-dimensional covalent organic frameworks (2D COFs), a process essential to extending their great potentials in advanced semiconducting applications, is yet fully elucidated. Herein, we report a systematic investigation of the photophysical properties of three highly crystalline imine-linked 2D COFs using steady-state and transient absorption spectroscopy accompanied by time-dependent density functional theory (TDDFT) calculations. The different electron affinity between 5,5',5''-(1,3,5-benzenetriyl)-tris(2-pyridinecarboxaldehyde) (BTPA) and three tunable electron-donating/accepting triamine monomers dominated the formation of the excited state, charge transfer direction and lifetime. A prominent charge transfer from electron-rich 4,4',4''-triaminotriphenylamine (TAPA) to BTPA in COFTAPA-BTPA led to the long-lived charge-separated state, which was attributed to a greater degree of delocalization compared to the two other COFs. Furthermore, these results provide fundamental insight into the importance of structure-property correlation for designing advanced photoactive 2D COF materials with the efficient charge transfer and long-lived charge-separated state.

14 SOLAR ENERGY↗
DOE OSTI2025

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗
DOE OSTI2023

Unimolecular Reactions of E -Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules

The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH 2 OH)CHOO], for its reactions with H 2 O and (H 2 O) 2 , and for the reaction of the E-(CH 2 OH)CHOO…H 2 O complex with H 2 O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH 2 OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH 2 OH)CHOO…H 2 O + H 2 O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H 2 O) 2 . In addition, we find that the atmospheric lifetime of E-(CH 2 OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10 –6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH 2 OH)CHOO enhances its reactivity.

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

Nonadiabatic molecular dynamics simulation of C2H22+ in a strong laser field*

We investigate the alignment dependence of the strong laser dissociation dynamics of molecule C 2 H 2 2 + in the frame of real-time and real-space time-dependent density function theory coupled with nonadiabatic quantum molecular dynamics (TDDFT-MD) simulation. This work is based on a recent experiment study “ultrafast electron diffraction imaging of bond breaking in di-ionized acetylene” [Wolter et al , Science 354 , 308–312 (2016)]. Our simulations are in excellent agreement with the experimental data and the analysis confirms that the alignment dependence of the proton dissociation dynamics comes from the electron response of the driving laser pulse. Our results validate the ability of the TDDFT-MD method to reveal the underlying mechanism of experimentally observed and control molecular dissociation dynamics.

Physics↗
DOE OSTI2021

Establishing gold and platinum standards to 1 terapascal using shockless compression

New techniques are advancing the frontier of high-pressure physics beyond 1 terapascal, leading to new discoveries and offering stringent tests for condensed-matter theory and advanced numerical methods. However, the ability to absolutely determine the pressure state remains challenging, and well-calibrated pressure-density reference materials are required. We conducted shockless dynamic compression experiments at the National Ignition Facility and the Z machine to obtain quasi-absolute, high-precision, pressure-density equation-of-state data for gold and platinum. We derived two experimentally constrained pressure standards to terapascal conditions. Establishing accurate experimental determinations of extreme pressure will facilitate better connections between experiments and theory, paving the way toward improving our understanding of material response to these extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Crystal nucleation rates in one-component Yukawa systems

Nucleation in the supercooled Yukawa system is relevant for addressing current challenges in understanding a range of crystallizing systems including white dwarf (WD) stars. We use both brute force and seeded molecular dynamics simulations to study homogeneous nucleation of crystals from supercooled Yukawa liquids. With our improved approach to seeded simulations, we obtain quantitative predictions of the crystal nucleation rate and cluster size distributions as a function of temperature and screening length. These quantitative results show trends towards fast nucleation with short-ranged potentials. They also indicate that for temperatures T > 0.9⁢T m , where T m is the melt temperature, classical homogeneous nucleation is too slow to initiate crystallization but transient clusters of ~100 particles should be common. As a result, we apply these general results to a typical WD model and obtain a delay of ~0.6 Gyr in the onset of crystallization that may be observable.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hamiltonian structure of 2D fluid dynamics with broken parity

Isotropic fluids in two spatial dimensions can break parity symmetry and sustain transverse stresses which do not lead to dissipation. Corresponding transport coefficients include odd viscosity, odd torque, and odd pressure. We consider an isotropic Galilean invariant fluid dynamics in the adiabatic regime with momentum and particle density conservation. We find conditions on transport coefficients that correspond to dissipationless and separately to Hamiltonian fluid dynamics. The restriction on the transport coefficients will help identify what kind of hydrodynamics can be obtained by coarse-graining a microscopic Hamiltonian system. Interestingly, not all parity-breaking transport coefficients lead to energy conservation and, generally, the fluid dynamics is energy conserving but not Hamiltonian. We show how this dynamics can be realized by imposing a nonholonomic constraint on the Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Electrical conductivity of Sn at high pressure and temperature

To date, temperature and conductivity have many outstanding implications in extreme environments but are yet to be fully understood under high pressure and temperature dynamic conditions. Here, we introduce an approach to provide high quality electrical conductivity results under dynamic loading conditions. Emphasis is given to address the skin depth effect's influence in a dynamic loading experiment by using thin films. The thin film samples in this study were at least 100 times thinner than previous samples in dynamic electrical conductivity experiments, increasing the current density to its full potential across the sample's entire cross section. Consideration of the skin depth accounts for at minimum a 4× scaling factor to the final electrical conductivity result that has been neglected in previous dynamic electrical conductivity studies. We also obtained improved signal-to-noise ratio with custom diagnostics optimized for better electrical impedance matching. Furthermore, these considerations were applied to Sn to assess electrical conductivity at elevated pressure and temperature. The high signal-to-noise ratio with reduced skin depth influence results in Sn allowing observation of the conductivity changes related to solid-to-solid and solid-to-liquid phase transitions. Additionally, we calculate the Sn thermal conductivity using the Wiedemann-Franz law for our experiments and compare against thermal transport dependent temperature measurements from previous work.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanism of Solid Electrolyte Interphase Formation on Lithium Anodes Studied by Molecular Dynamics Simulations

We investigated the initial formation of the SEI on lithium-metal anodes in presence of several electrolyte compositions using ReaxFF, in absence of a bias potential. We observed the formation of well-defined regions with different density and oxidation states of Li. Dissolution of Li metal atoms results in an uneven distribution of Li in the liquid phase. The electrolyte decomposition mechanism was studied for various electrolyte mixtures. Overall, the reaction mechanisms observed with ReaxFF are in good agreement with previous reports using DFT. However, we found some discrepancies in specific predictions by the ReaxFF that should be addressed in future work.

25 ENERGY STORAGE↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
DOE OSTI2021

Light–matter interaction of a molecule in a dissipative cavity from first principles

Cavity-mediated light–matter coupling can dramatically alter opto-electronic and physico-chemical properties of a molecule. Ab initio theoretical predictions of these systems need to combine non-perturbative, many-body electronic structure theory-based methods with cavity quantum electrodynamics and theories of open-quantum systems. Here, we generalize quantum-electrodynamical density functional theory to account for dissipative dynamics of the cavity and describe coupled cavity–single molecule interactions in the weak-to-strong-coupling regimes. Specifically, to establish this generalized technique, we study excited-state dynamics and spectral responses of benzene and toluene under weak-to-strong light–matter coupling. By tuning the coupling, we achieve cavity-mediated energy transfer between electronically excited states. This generalized ab initio quantum-electrodynamical density functional theory treatment can be naturally extended to describe cavity-mediated interactions in arbitrary electromagnetic environments, accessing correlated light–matter observables and thereby closing the gap between electronic structure theory, quantum optics, and nanophotonics.

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

Connectivity-Dependent Exciton–Phonon Coupling in Cesium Bismuth Halide Quantum Dots

Metal halide octahedra form the fundamental functional building blocks of metal halide perovskites, dictating their structures, optical properties, electronic structures, and dynamics. Here, in this study, we show that the connectivity of bismuth halide octahedra in Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 quantum dots (QDs) changes with different halide elements. We use first-principles calculations to reveal the key role of the connectivity of bismuth halide octahedra on the wave function symmetry, Huang-Rhys factor, and exciton-phonon interaction strength. Following QD synthesis via a ligand-mediated transport method, the effect of connectivity is verified with transient absorption spectroscopy, where we contrast Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 QD exciton dynamics. In photoexcited Cs 3 Bi 2 I 9 QDs, phonons related to the vibrational motions of face-sharing [BiI 6 ] 3- bioctahedra couple strongly to the electronic state and drive rapid carrier relaxation. Equivalent signals are not observed for photoexcited Cs 3 Bi 2 Br 9 QDs, implying a lack of phonon involvement in band-edge absorption and subsequent exciton relaxation. Our findings suggest that structural engineering can effectively tune the exciton-phonon coupling and therefore influence exciton relaxation and recombination in perovskite nanomaterials.

TDDFT↗
DOE OSTI2024

Pd-Ru pair on Pt surface for promoting hydrogen oxidation and evolution in alkaline media

Hydrogen oxidation reaction in alkaline media is critical for alkaline fuel cells and electrochemical ammonia compressors. The slow hydrogen oxidation reaction in alkaline electrolytes requires large amounts of scarce and expensive platinum catalysts. While transition metal decoration can enhance Pt catalysts’ activity, it often reduces the electrochemical active surface area, limiting the improvement in Pt mass activity. Here, we enhance Pt catalysts’ activity without losing surface-active sites by using a Pd-Ru pair. Utilizing a mildly catalytic thermal pyrolysis approach, Pd-Ru pairs are decorated on Pt, confirmed by extended X-ray absorption fine structure and high-angle annular dark-field scanning transmission electron microscopy. Density functional theory and ab-initio molecular dynamics simulations indicate preferred Pd and Ru dopant adsorption. The Pd-Ru decorated Pt catalyst exhibits a mass-based exchange current density of 1557 ± 85 A g –1 metal for hydrogen oxidation reaction, demonstrating superior performance in an ammonia compressor.

25 ENERGY STORAGE↗