Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density matrix methods”

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.

560 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
DOE OSTI2023

A Self‐Healing System for Polydicyclopentadiene Thermosets

Self-healing offers promise for addressing structural failures, increasing lifespan, and improving durability in polymeric materials. Implementing self-healing in thermoset polymers faces significant manufacturing challenges, especially due to the elevated temperature requirements of thermoset processing. To introduce self-healing into structural thermosets, the self-healing system must be thermally stable and compatible with the thermoset chemistry. This article demonstrates a self-healing microcapsule-based system stable to frontal polymerization (FP), a rapid and energy-efficient manufacturing process with a self-propagating exothermic reaction (≈200 °C). A thermally latent Grubbs-type complex bearing two N-heterocyclic carbene ligands addresses limitations in conventional G2-based self-healing approaches. Under FP's elevated temperatures, the catalyst remains dormant until activated by a Cu(I) co-reagent, ensuring efficient polymerization of the dicyclopentadiene (DCPD) upon damage to the polyDCPD matrix. The two-part microcapsule system consists of one capsule containing the thermally latent Grubbs-type catalyst dissolved in the solvent, and another capsule containing a Cu(I) coagent blended with liquid DCPD monomer. Using the same chemistry for both matrix fabrication and healing results in strong interfaces as demonstrated by lap-shear tests. In an optimized system, the self-healing system restores the mechanical properties of the tough polyDCPD thermoset. Self-healing efficiencies greater than 90% via tapered double cantilever beam tests are observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantifying the dislocation content of atomically resolved grain boundary line defects using the Nye tensor

The Nye tensor, which quantifies the density of Burgers vector for a given dislocation line direction, can be effectively used to characterize dislocation content in bulk crystals from atomic-resolution transmission electron microscopy images. The Nye tensor can be calculated from these images, in part because the reference state is simply defined by the lattice of the perfect crystal. The application of the Nye tensor to interfacial line defects, for which the natural reference state is the dichromatic pattern of the two grains in their reference orientation, poses additional challenges. In this work, we present a method that employs the Nye tensor to characterize the edge dislocation content of line defects at grain boundaries from atomic-resolution images. This approach enables us to rapidly characterize all edge dislocation content along a grain boundary. Additionally, the Nye tensor provides information about line defect core structure. Finally, we demonstrate this method on two exemplar defects: a twin boundary disconnection and a facet junction in face-centered cubic Au.

Crystallographic defects↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2025

Echo state network for coarsening dynamics of charge density waves

An echo state network (ESN) is a type of reservoir computer that uses a recurrent neural network with a sparsely connected hidden layer. Compared with other recurrent neural networks, one great advantage of ESN is the simplicity of its training process. Yet, despite the seemingly restricted learnable parameters, ESN has been shown to successfully capture the spatial-temporal dynamics of complex patterns. Here we build an ESN to model the coarsening dynamics of charge-density waves (CDWs) in a semiclassical Holstein model, which exhibits a checkerboard electron density modulation at half-filling stabilized by a commensurate lattice distortion. The inputs to the ESN are local CDW order parameters in a finite neighborhood centered around a given site, while the output is the predicted CDW order of the center site at the next time step. Special care is taken in the design of couplings between hidden layer and input nodes to ensure lattice symmetries are properly incorporated into the ESN model. Since the model predictions depend only on CDW configurations of a finite domain, the ESN is scalable and transferrable in the sense that a model trained on dataset from a small system can be directly applied to dynamical simulations on larger lattices. Furthermore, our work opens avenues for efficient dynamical modeling of pattern formations in functional electron materials.

2-dimensional systems↗
DOE OSTI2026

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Real-space visualization of a defect-mediated charge density wave transition

Here, we study the coupled charge density wave (CDW) and insulator-to-metal transitions in the 2D quantum material 1T-TaS 2 . By applying in situ cryogenic 4D scanning transmission electron microscopy with in situ electrical resistance measurements, we directly visualize the CDW transition and establish that the transition is mediated by basal dislocations (stacking solitons). We find that dislocations can both nucleate and pin the transition and locally alter the transition temperature T c by nearly ~75 K. This finding was enabled by the application of unsupervised machine learning to cluster five-dimensional, terabyte scale datasets, which demonstrate a one-to-one correlation between resistance—a global property—and local CDW domain-dislocation dynamics, thereby linking the material microstructure to device properties. This work represents a major step toward defect-engineering of quantum materials, which will become increasingly important as we aim to utilize such materials in real devices.

4D-STEM↗
DOE OSTI2024

Converting Second‐Order Saddle Points to Transition States: New Principles for the Design of 4π Photoswitches

Abstract Molecular photoswitches have demonstrated potential for storing solar energy at the molecular level, with power densities comparable to commercial batteries and hydroelectric energy storage. However, development of efficient photoswitches is hindered by limitations in cyclability and optical properties of existing materials. We here demonstrate that certain limitations in photoswitches based on electrocyclizations stem from the issue of controlling competition between Woodward‐Hoffmann allowed and forbidden pathways. Our approach moves beyond the traditional view of activation barriers and reveals that second‐order saddle points are crucial in dictating the competition between disrotatory and conrotatory pathways. These insights suggest new opportunities to manipulate the competition between these pathways through geometric constraints, fundamentally altering the connectivity of the potential energy surface. Our study also emphasizes the necessity of multi‐reference methods and the need to conduct higher‐dimensional explorations for competing pathways beyond photoswitch design.

Chemistry↗
DOE OSTI2024

High-Density Capacitive Energy Storage in Low-Dielectric-Constant Polymer PMMA/2D Mica Nanofillers Heterostructure Composite

The ubiquitous, rising demand for energy storage devices with ultra-high storage capacity and efficiency has drawn tremendous research interest in developing energy storage devices. Dielectric polymers are one of the most suitable materials used to fabricate electrostatic capacitive energy storage devices with thin-film geometry with high power density. In this work, we studied the dielectric properties, electric polarization, and energy density of PMMA/2D Mica nanocomposite capacitors where stratified 2D nanofillers are interfaced between the multiple layers of PMMA thin films using two heterostructure designs of the capacitors, PMMA/2D Mica/PMMA (PMP) and PMMA/2D Mica/PMMA/2D Mica/PMMA (PMPMP). The incorporation of a 2D Mica nanofiller in the low-dielectric-constant PMMA leads to an enhancement in the dielectric constant, with ∆ε ~ 15% and 53% for PMP and PMPMP heterostructures at room temperature. Additionally, a significant improvement in discharged energy density was measured for the PMPMP capacitor (Ud ~ 38 J/cm3 at 825 MV/m) compared to the pristine PMMA (Ud ~ 9.5 J/cm3 at 522 MV/m) and PMP capacitors (Ud ~ 19 J/cm3 at 740 MV/m). This excellent capacitive and energy storage performance of the PMMA/2D Mica heterostructure nanocomposite may inform the fabrication of thin-film, high-density energy storage capacitor devices for potential applications in various platforms.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Imprints of High-Density Nuclear Symmetry Energy on Crustal Fraction of Neutron Star Moment of Inertia

The density dependence of nuclear symmetry energy E sym (ρ) remains the most uncertain aspect of the equation of state (EOS) of supradense neutron-rich nucleonic matter. Utilizing an isospin-dependent parameterization of the nuclear EOS, we investigate the implications of the observational crustal fraction of the neutron star (NS) moment of inertia ΔI/I for the E sym (ρ). We find that symmetry energy parameters significantly influence the ΔI/I, while the EOS of symmetric nuclear matter has a negligible effect. In particular, an increase in the slope L and skewness J sym of symmetry energy results in a larger ΔI/I, whereas an increase in the curvature K sym leads to a reduction in ΔI/I. Moreover, the ΔI/I is shown to have the potential for setting a lower limit of symmetry energy at densities exceeding 3 ρ 0 , particularly when L is constrained to values less than 60 MeV, thereby enhancing our understanding of supradense NS matter.

equation of state↗
DOE OSTI2026

Characterization of ceramic fuel powder packing fractions to support INFLUX

Triply periodic minimal surface (TPMS)-based structures show marked potential in novel nuclear reactor fuel designs, as their high surface area-to-volume ratio increases the efficiency of heat transfer out of the fuel, enabling safer, more innovative reactor designs. This milestone report addresses the role of dUO 2 powder processing route on the fill behavior of TPMS-based cladding shells to understand and advance the feasibility of manufacturing TPMS-based nuclear fuel forms. dUO 2 powder was processed through either a dry granulation route, varying consolidation pressure, or through milling, varying milling time, milling method and milled size distribution. The lowest tapped bulk densities (TBD), but best powder flowabilities, were obtained when testing unprocessed dUO 2 powder which was prone to self-agglomeration and formed low-density spheroids. The highest TBD and lowest flowabilities were obtained when using powder produced by hammer-milling dUO 2 powder to pass through a 200-mesh sieve, which led to particles with angular morphologies. Powder produced by dry granulation exhibited TBD that varied according to the consolidation pressure used to form the initial pellets and exhibited improved flowabilities when compared to hammer-milled material. Because of the large span of granule sizes formed as well as the irregular shape associated with the granules, a packing fraction of 0.69 was achieved, exceeding the analytical solution for random close packing of mono-sized spheres. TPMS polymer shells were loaded with unprocessed, granulated, and hammer-milled dUO 2 powders, and their qualitative packing behaviors were analyzed using x-ray computed tomography (xCT). TBDs calculated after loading TPMS polymer shells were 10-20% lower when compared to tapped bulk density measurements taken in a glass graduated cylinder, indicating a non-trivial impact on the tapped bulk density of either the TPMS channel size, TPMS channel surface material, powder cohesiveness, or a combination of the two parameters. A metallic zircaloy-4 TPMS shell will be loaded with hammer-milled dUO 2 powder upon receipt of the shell from Oak Ridge National Laboratory (ORNL) and shipped to Idaho National Labs (INL) for subsequent hot isostatic pressing (HIP) densification experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗
DOE OSTI2024

Crack‐Free Single‐Crystalline LiNiO 2 for High Energy Density All‐Solid‐State Batteries

Abstract Single‐crystalline layered oxide (LiNi 1‐ x ‐ y Mn x Co y O 2 ) cathodes have been found to exhibit exceptional electrochemical properties when coupled with various inorganic solid electrolytes (ISEs) in all‐solid‐state batteries (ASSBs). Their advantages stem from the robust morphological integrity with the absence of grain boundaries and the high electrochemical oxidative stability. Here, ASSBs featuring single‐crystalline LiNiO 2 (LNO) with the highest Ni content are reported, offering a high theoretical specific capacity of 275 mAh g ‐1 alongside a high average discharge voltage (3.7 V vs Li + /Li). Through a careful investigation, it is demonstrated that micron‐sized single‐crystalline LNO (µSC‐LNO) composite cathodes with a halide ISE exhibit a high initial discharge capacity of 205 mAh g ‐1 with an outstanding cycle performance over 200 cycles in room‐temperature ASSBs. The significance of engineering parameters is emphasized, such as particle size and specific density, in promoting a homogeneous and fast Li + transport within the composite cathodes. Furthermore, the formation of undesirable interphase between the halide ISE in the cathode and sulfide ISE separator is elucidated, which may be a critical factor impeding long‐term cyclability of ASSBs. This work provides insights into the design of composite cathodes for high‐energy‐density ASSBs.

Chemistry↗
DOE OSTI2026

Density Functional Theory Analysis of Candidate Tritium Permeation Barriers for Vanadium Alloys as a Fusion Structural Material

Vanadium alloys are being considered as a candidate structural material for the breeder blanket component of future fusion reactors due to their good thermal conductivity, low activation, and chemical compatibility with FLiBe. However, vanadium has a high hydrogen solubility, resulting in hydrogen diffusion and trapping in the vanadium structural alloy and lower than expected tritium breeding ratios. To mitigate the permeation of hydrogen into the structural alloy, a tritium permeation barrier will be coated onto the vanadium alloy. Multiple materials have been proposed as candidate tritium permeation barriers, but the extent to which each of these materials impacts hydrogen diffusion and trapping is still unknown. In this project, we perform atomic-scale calculations to narrow down the list of potential tritium permeation barrier candidates and to understand the mechanisms through which hydrogen diffusion is affected in these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiple shape coexistence near {sup 118}Sn: First 0+ 3 lifetime measurement

The intruder bands in Sn isotopes, built on the 2p-2h excitation across the Z = 50 proton shell gap, are well-known examples of shape coexistence near the neutron mid-shell region. Spectroscopic signatures for shape coexistence include enhanced E0 transitions between the 0+ band heads. However, the underlying shape coexistence and mixing has been unclear because lifetime information for the excited 0+ states was incomplete in 118Sn. We thus present here the first measurement of the 0+ 3 lifetime in 118Sn using the fast-timing technique following thermal-neutron capture. The observed enhancement in ρ2(E0;0+ 3 → 0+ 2 ) of 150(30) milliunits provides compelling indications for multiple shape coexistence in 118Sn. Additionally, three distinct shapes in 116,118,120Sn naturally emerged in theoretical calculations based on the quantum-number-projected generator coordinate method employing a relativistic energy density functional.

Wu, F.↗
DOE OSTI2025

Dimensional Evolution Guides Property Control in the A n Cu 4– n TiS 4 Semiconductor Series

Through progressive reduction of the three-dimensional (3D) covalent network of Cu 4 TiS 4 , we isolate seven new members of the A n Cu 4–n TiS 4 family (A = alkali metal; n = 0–4), spanning 3D, 2D, 1D, and 0D structural fragments. The dimensional reduction is rational, as it preserves the edge-sharing connectivity between [CuS 4 ] 7– and [TiS 4 ] 4– tetrahedra across the series. This structural evolution is driven by the stepwise substitution of Cu with alkali metals, guiding the formation of fragments with reduced dimensionality. The effects of “n” and “A” on the crystal structures, stabilities, electronic structures, and optoelectronic properties are profound, demonstrating that the manipulation of alkali metal size and A n Cu 4–n TiS 4 stoichiometry enables predictable variations in structure and properties. For example, the n = 0 and n = 4 end members of the A n Cu 4–n TiS 4 family set the range of achievable band gaps with 2.00 eV for Cu 4 TiS 4 , 2.60 eV for Na 4 TiS 4 , and intermediate values for the n = 1–3 members. Notably, CsCu 3 TiS 4 exhibits exceptional air stability and congruent melting, with density functional theory (DFT) calculating moderate hole and electron effective masses in specific crystallographic directions (mh = 1.24m 0 , me = 0.87m 0 ). Additionally, A 3 CuTiS 4 (A = Na, K, Rb) displays direct band gap behavior and long photoluminescence lifetimes of 2.3–8.6 μs, and K 3 CuTiS 4 has a PLQY of 5.19%. These findings underscore the potential of the A n Cu 4–n TiS 4 family for applications in optoelectronics and demonstrate widely applicable design concepts that unveil rational stoichiometries within a given composition space to generate a series of crystal structures related through an evolving covalent dimensionality that corresponds to a predictable electronic structure and property progression.

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

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY↗
DOE OSTI2025

Generalized geometric speed limits for quantum observables

Leveraging quantum information geometry, we derive generalized quantum speed limits on the rate of change of the expectation values of observables. These bounds subsume and, for Hilbert space dimension ≥3, tighten existing bounds—in some cases by an arbitrarily large multiplicative constant. Our theoretical results are supported by illustrative examples and an experimental demonstration using a superconducting qutrit. We also derive two upper bounds on the generalized quantum Fisher information in terms of the condition number of the density matrix. One of these bounds applies only to coherent dynamics and depends also on the variance of the Hamiltonian. The other bound depends also on the so-called Wigner-Yanase skew information. These bounds generalize well-known bounds on the symmetric logarithmic derivative quantum Fisher information and are tighter than the existing bounds for sufficiently mixed states (e.g., for sufficiently high temperature thermal states).

open quantum systems & decoherence↗