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 55 records · Page 3

Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models

A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗

Associative polymers with controlled sticker placement: How reversible bond distribution and density govern polymer dynamics

Associative polymers with precisely arranged stickers offer opportunities to program material properties with molecular precision. Yet, it remains unclear how the placement and fraction of stickers dictate structure, dynamics, and macroscopic properties. By developing a model unentangled polymer system with hydrogen-bonding stickers, we show that randomly distributed stickers neither form clusters nor change flow properties, whereas stickers placed at chain ends drive nanocluster formation even at low concentrations. Adding more end stickers produces a rubbery plateau spanning eight decades in frequency with two distinct relaxation timescales, in contrast to the single plateau predicted by the classic sticky Rouse model. These results demonstrate that sticker distribution dictates whether associative polymers undergo nanocluster formation or microphase separation, while substantial alterations in dynamics and viscoelasticity require both sticker aggregation and thermomechanical stability of associated domains. Our findings resolve a longstanding debate on associative polymer dynamics and provide molecular design rules for programmable soft materials.

36 MATERIALS SCIENCE↗

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. Here, in this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Hamiltonians↗

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula

Here, we study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Incorporating point defect generation due to jog formation into the vector density-based continuum dislocation dynamics approach

During plastic deformation of crystalline materials, point defects such as vacancies and interstitials are generated by jogs on moving dislocations. A detailed model for jog formation and transport during plastic deformation was developed within the vector density-based continuum dislocation dynamics framework. Here as a part of this model, point defect generation associated with jog transport was formulated in terms of the volume change due to the non-conservative motion of jogs. Balance equations for the vacancies and interstitials including their rate of generation due to jog transport were also formulated. A two-way coupling between point defects and dislocation dynamics was then completed by including the stress contributed by the eigen-strain of point defects. A jog drag stress was further introduced into the mobility law of dislocations to account for the energy dissipation during point defects generation. A number of test problems and a fully coupled simulation of dislocation dynamics and point defects generation and diffusion were performed. The results show that there is an asymmetry of vacancy and interstitial generation due to the different formation energies of the two types of defects. The results also show that a higher hardening rate and a higher dislocation density are obtained when the point defect generation mechanism is coupled to dislocation dynamics.

36 MATERIALS SCIENCE↗

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synergistic Density Functional Theory and Molecular Dynamics Approach to Elucidate PNIPAM–Water Interaction Mechanisms

This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.

08 HYDROGEN↗

Charge correlations and magnetoelastic coupling in intercalated transition metal dichalcogenides

The large van der Waals gap in transition metal dichalcogenides (TMDs) offers an avenue to tune the ground state of 2D materials through the intercalation of magnetic atoms. Here, we investigate the charge correlations in Fe 1/3 ⁢TaS 2 , Co 1/3 ⁢TaS 2 , and Fe 0.35 ⁢NbS 2 by combining angle-resolved photoemission spectroscopy (ARPES), x-ray scattering, magnetometry, and density functional theory (DFT). We find that, while short-range charge fluctuations develop in Ta-based compounds, Fe 0.35 ⁢NbS 2 exhibits long-range charge order which is strongly coupled with magnetic order and tunable by external magnetic field. Our electronic structure analysis reveals that intercalation reconstructs the Fermi surface via charge transfer and band renormalization, yet does not generate the nesting conditions compatible with the observed ordering vectors. Complementary phonon calculations further exclude a conventional electron-phonon origin of charge order. Together, these results establish magnetoelastic coupling as the dominant mechanism behind charge ordering in Fe 0.35 ⁢NbS 2 and highlight the contrasting role of Nb and Ta hosts in stabilizing correlated ground states in intercalated TMDs.

36 MATERIALS SCIENCE↗

DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various density functional theory packages

We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗

Local observation of linear-T superfluid density and anomalous vortex dynamics in URu2Si2

The heavy fermion superconductor URu2Si2 is a candidate for chiral, time-reversal symmetry-breaking superconductivity with a nodal gap structure. Here, we microscopically visualized superconductivity and spatially inhomogeneous ferromagnetism in URu2Si2. We observed linear-T superfluid density, consistent with d-wave pairing symmetries including chiral d wave, but did not observe the spontaneous magnetization expected for chiral d wave. Local vortex pinning potentials had either four- or twofold rotational symmetries with various orientations at different locations. Taken together, these data support a nodal gap structure in URu2Si2 and suggest that chirality either is not present or does not lead to detectable spontaneous magnetization.

36 MATERIALS SCIENCE↗

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Dynamics of Ions on Phased Radio-frequency Carpets in High Pressure Gases and Application for Barium Tagging in Xenon Gas Time Projection Chambers

Radio-frequency (RF) carpets with ultra-fine pitches are examined for ion transport in gases at atmospheric pressures and above. We develop new analytic and computational methods for modeling RF ion transport at densities where dynamics are strongly influenced by buffer gas collisions. An analytic description of levitating and sweeping forces from phased arrays is obtained, then thermodynamic and kinetic principles are used to calculate ion loss rates in the presence of collisions. This methodology is validated against detailed microscopic SIMION simulations. We then explore a parameter space of special interest for neutrinoless double beta decay experiments: transport of barium ions in xenon at pressures from 1 to 10 bar. Our computations account for molecular ion formation and pressure dependent mobility as well as finite temperature effects. We discuss the challenges associated with achieving suitable operating conditions, which lie beyond the capabilities of existing devices, using presently available or near-future manufacturing techniques.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Magnetism and finite-temperature effects in UZr 2 : A density functional theory analysis

The structure and the thermophysical properties of d-UZr 2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Finally, our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.

36 MATERIALS SCIENCE↗

Defect Production and Microstructural Feature Impact for Radiation Damage in Additively Manufactured 316 Stainless Steel

This milestone presents multi-scale modeling research results for additively manufactured 316 stainless steel. A combination of phase field, cluster dynamics, molecular dynamics, and density functional theory with machine learning is used, allowing for predictions of radiation-driven microstructural evolution in additively manufactured 316 stainless steel over a range of temperatures, damage rates, neutron spectra, and microstructures, and supporting the development of combined ion and neutron irradiation for material qualification. Informed by ion irradiation and neutron irradiation results across the Advanced Materials and Manufacturing Technologies program, we investigate the unique aspects of radiation-driven microstructure evolution in additively manufactured 316 stainless steel. In particular, we focus on understanding the impact of carbon concentration (varying, for example, between the 316L and 316H standards) on void formation; radiation-induced segregation at dislocation cells and grain boundaries; and the evolution of dislocation loops and network dislocation populations. We find that the unique characteristics of the additively manufactured microstructures must be accounted for in understanding the evolution of dislocation populations under thermal and irradiation conditions, such as the variation in sink strengths arising due to the variation in dislocation density. We also find that the radiation-induced segregation of Cr and Ni to grain boundaries and cell walls differs due to the differences in their defect sink biases. We also find that increasing the Ni content can slow vacancy diffusion, which may provide a mechanism for the observed reduction in transient swelling rate for austenitic Fe-Cr-Ni alloys with increasing Ni content. In addition, ion irradiations have shown that increasing carbon content in 316 SS results in a larger population of smaller voids, suggesting reduced vacancy diffusion. Our results show that the carbon content of additively manufactured 316 SS has a significant impact on the migration rate of defect clusters. The presence of carbon atoms results in carbon-vacancy trapping, significantly reducing the diffusion rate of vacancies. Carbon atoms may also be trapped near the surface of a void, which may reduce void growth by trapping vacancies that diffuse toward the void.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗