Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum simulation”

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.

448 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Trust-Based Detection and Mitigation of Cyber Attacks in Distributed Cooperative Control of Islanded AC Microgrids

In this study, we address the challenge of detecting and mitigating cyber attacks in the distributed cooperative control of islanded AC microgrids, with a particular focus on detecting False Data Injection Attacks (FDIAs), a significant threat to the Smart Grid (SG). The SG integrates traditional power systems with communication networks, creating a complex system with numerous vulnerable links, making it a prime target for cyber attacks. These attacks can lead to the disclosure of private data, control network failures, and even blackouts. Unlike machine learning-based approaches that require extensive datasets and mathematical models dependent on accurate system modeling, our method is free from such dependencies. To enhance the microgrid’s resilience against these threats, we propose a resilient control algorithm by introducing a novel trustworthiness parameter into the traditional cooperative control algorithm. Our method evaluates the trustworthiness of distributed energy resources (DERs) based on their voltage measurements and exchanged information, using Kullback-Leibler (KL) divergence to dynamically adjust control actions. We validated our approach through simulations on both the IEEE-34 bus feeder system with eight DERs and a larger microgrid with twenty-two DERs. The results demonstrated a detection accuracy of around 100%, with millisecond range mitigation time, ensuring rapid system recovery. Additionally, our method improved system stability by up to almost 100% under attack scenarios, showcasing its effectiveness in promptly detecting attacks and maintaining system resilience. These findings highlight the potential of our approach to enhance the security and stability of microgrid systems in the face of cyber threats.

Computer Science↗
DOE OSTI2026

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE↗
DOE OSTI2026

Plasma flow generation and particle acceleration from expanding magnetic bubbles

Impulsive plasma dynamics in the laboratory are often driven by rising electric currents, yet their quantitative plasma response has not been well established. By means of fully kinetic particle-in-cell simulations and laser-driven capacitor-coil experiments, we show that a rising current expels plasma, forming an expanding magnetic bubble and accelerating particles. The expansion front velocity scales with the Alfvén speed determined by the magnetic field at its inner edge and the plasma density at its outer edge. This mechanism establishes impulsive current drive as a fundamental way that generates plasma flows and accelerates particles in laboratory plasmas, with potential relevance to astrophysics.

Zhang, Yang [Princeton University, NJ (United Stat↗
DOE OSTI2025

Detection of visible-wavelength aurora on Mars

Mars hosts various auroral processes despite the planet’s tenuous atmosphere and lack of a global magnetic field. To date, all aurora observations have been at ultraviolet wavelengths from orbit. We describe the discovery of green visible-wavelength aurora, originating from the atomic oxygen line at 557.7 nanometers, detected with the SuperCam and Mastcam-Z instruments on the Mars 2020 Perseverance rover. Near–real-time simulations of a Mars-directed coronal mass ejection (CME) provided sufficient lead-time to schedule an observation with the rover. The emission was observed 3 days after the CME eruption, suggesting that the aurora was induced by particles accelerated by the moving shock front. To our knowledge, detection of aurora from a planetary surface other than Earth has never been reported, nor has visible aurora been observed at Mars. This detection demonstrates that auroral forecasting at Mars is possible, and that during events with higher particle precipitation, or under less dusty atmospheric conditions, aurorae will be visible to future astronauts.

Science & Technology - Other Topics↗
DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗
DOE OSTI2026

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗
DOE OSTI2024

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2023

Investigating permafrost carbon dynamics in Alaska with artificial intelligence

Abstract Positive feedbacks between permafrost degradation and the release of soil carbon into the atmosphere impact land–atmosphere interactions, disrupt the global carbon cycle, and accelerate climate change. The widespread distribution of thawing permafrost is causing a cascade of geophysical and biochemical disturbances with global impacts. Currently, few earth system models account for permafrost carbon feedback (PCF) mechanisms. This research study integrates artificial intelligence (AI) tools and information derived from field-scale surveys across the tundra and boreal landscapes in Alaska. We identify and interpret the permafrost carbon cycling links and feedback sensitivities with GeoCryoAI, a hybridized multimodal deep learning (DL) architecture of stacked convolutionally layered, memory-encoded recurrent neural networks (NN). This framework integratesin-situmeasurements and flux tower observations for teacher forcing and model training. Preliminary experiments to quantify, validate, and forecast permafrost degradation and carbon efflux across Alaska demonstrate the fidelity of this data-driven architecture. More specifically, GeoCryoAI logs the ecological memory and effectively learns covariate dynamics while demonstrating an aptitude to simulate and forecast PCF dynamics—active layer thickness (ALT), carbon dioxide flux (CO 2 ), and methane flux (CH 4 )—with high precision and minimal loss (i.e. ALT RMSE : 1.327 cm [1969–2022]; CO 2 RMSE : 0.697µmolCO 2 m −2 s −1 [2003–2021]; CH 4 RMSE : 0.715 nmolCH 4 m −2 s −1 [2011–2022]). ALT variability is a sensitive harbinger of change, a unique signal characterizing the PCF, and our model is the first characterization of these dynamics across space and time.

Environmental Sciences & Ecology↗
DOE OSTI2026

Update on Radiochemical Assessment of High Burnup Commercially Irradiated Fuel

This work documents an effort to collect burnup measurements on a high burnup rod, designated 6XV, and first cycle accident tolerant fuel (ATF) rod, designated 47I, to enable benchmarking of fuel performance codes and neutronics codes. In addition to measurements, Virtual Environment for Reactor Applications (VERA) full-core-depletion analysis was also performed for the rods that were experimentally analyzed to provide an opportunity for code validation. This effort focuses on collecting data from rods irradiated at Byron Generating Station and shipped to the Oak Ridge National Laboratory (ORNL) hot-cells. This data will also anchor non-destructive examination evaluations of burnup of the various fuel rods undergoing postirradiation examination (PIE) at ORNL. Previous PIE of these fuel rods provides some guidance on the burnup trend across the fuel. Axial gamma spectroscopy scans provide a measure of relative changes in burnup across a fuel pin. Mass spectrometry based burnup measurements performed for this work at specific axial locations in the fuel are fully quantitative. By combining the mass spectrometry data with the gamma scans it is possible to more quantitatively evaluate axial variations in burnup across the entire fuel pin [1]. The combined set of burnup evaluations will be made available to other organizations that have an interest in high burnup radiochemistry data for validation of neutronic simulations and source term evaluation such as the Nuclear Regulatory Commission (NRC).

Harp, Jason [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2026

Activation dynamics of a water-soluble human mu-opioid receptor

The mu-opioid receptor (MOR), a class A G protein-coupled receptor mediates opioid analgesia and remains a central target for pain therapeutics. While crystal structures of MOR exist, they provide limited insight into the receptor’s dynamic conformational landscape underlying function. Here, we engineered a thermostable water-soluble MOR variant (wsMOR) that retains native-like ligand-binding and activation dynamics. This variant enables high-yield production and detailed solution-phase structural studies that are challenging with membrane-embedded MOR, providing a valuable tool for studying receptor activation and aqueous-phase drug screening. Using a combined computational and experimental approach, we performed long-timescale all-atom molecular dynamics simulations together with neutron scattering and single-molecule FRET, revealing a structurally stable receptor with a diverse ensemble of conformations at different temporal resolutions. In the ligand-free state, wsMOR displayed high conformational flexibility, which decreased upon agonist binding, particularly in transmembrane helix 6, a hallmark of G protein-coupled receptor activation. Positive allosteric modulation and G protein binding further stabilized active-like states. These findings highlight wsMOR’s conformational plasticity across picosecond to millisecond timescales and provide a foundation for structure-guided development of next-generation opioid ligands with improved efficacy and safety.

E, Agyemang [University of Tennessee Knoxville]↗
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2024

Tailoring Ion Transport in Li 3‐3y Ho 1+y Cl 6‐x Br x via Transition‐Metal Free Structural Planes and Charge Carrier Distribution

Abstract Localized atomistic disorder in halide‐based solid electrolytes (SEs) can be leveraged to boost Li + mobility. In this study, Li + transport in structurally modified Li 3 HoCl 6 , via Br − introduction and Li + deficiency, is explored. The optimized Li 3‐3 y Ho 1+ y Cl 6‐ x Br x achieves an ionic conductivity of 3.8 mS cm −1 at 25 °C, the highest reported for holmium halide materials. 6,7 Li nuclear magnetic resonance and relaxometry investigations unveil enhanced ion dynamics with bromination, attaining a Li + motional rate neighboring 116 MHz. X‐ray diffraction analyses reveal mixed‐anion‐induced phase transitions with disproportionate octahedral expansions and distortions, creating Ho‐free planes with favorable energetics for Li + migration. Bond valence site energy analysis highlights preferred Li + transport pathways, particularly in structural planes devoid of Ho 3+ blocking effects. Molecular dynamics simulations corroborate enhanced Li + diffusion with Br − introduction into Li 3 HoCl 6 . Li‐Ho electrostatic repulsions in the (001) plane presumably drive Li + diffusion into the Ho‐free (002) layer, enabling rapid intraplanar Li + motion and exchange between the 2d and 4h sites. Li 3‐3 y Ho 1+ y Cl 6‐ x Br x also demonstrates good battery cycling stability. These findings offer valuable insights into the intricate correlations between structure and ion transport and will help guide the design of high‐performance fast ion conductors for all‐solid‐state batteries.

Chemistry↗
DOE OSTI2026

Dynamics of Radiation Damage Buildup in Ultrathin Hexagonal Boron Nitride Films under Ion Bombardment

Two-dimensional hexagonal boron nitride (hBN) is attractive for several emerging applications. Ion bombardment can be used to modify the hBN properties. However, the understanding of radiation damage buildup in hBN remains limited. Here, we investigate the effects of the dose rate and ion mass on radiation damage buildup by studying 40 nm-thick hBN films bombarded at room temperature with 500 keV 4 He, 15 N, 40 Ar, and 129 Xe ions and comparing with results for ion bombardment of polycrystalline hBN ceramics. Raman spectroscopy is used to quantify damage buildup, and transmission electron microscopy is used for microstructural analysis. Experiments are complemented by molecular dynamics simulations of the formation and evolution of point defects. Lighter ions are found to be more efficient at disordering hBN than heavier ions. This observation points to a critical role of intracascade defect processes. In contrast, a negligible dose rate effect observed suggests limited intercascade defect dynamic annealing processes for these irradiation conditions. These findings provide a fundamental basis for hBN defect engineering.

2D materials↗
DOE OSTI2026

Robustness of inertial fusion energy relevant implosions to low-mode asymmetries

Inertial confinement fusion implosions for power production need to be robust. Degradation sources need to be understood and considered for an economic power plant. In this work, the impact of low-mode asymmetries on yield is studied using 2D radiation hydrodynamics simulations for inertial fusion energy relevant implosions. First, a high-gain indirectly driven capsule design is developed. This design is subsequently degraded with low-mode drive asymmetries. The results show that yield is constant with increasing amounts of asymmetry until a cliff is reached where the yield quickly drops below ignition. Robustness arises because the degraded hotspot still has enough time to get to ignition before capsule decompression. This robustness occurs even though the asymmetry couples less energy to the hotspot and causes decompression to occur sooner. The robustness is linked to the generalized Lawson parameter, χ⁠. The value of χ steadily decreases with increasing asymmetry. Once χ is close to one, the hotspot is no longer able to ignite before the capsule decompresses. A simple piston model is used to predict where the robustness cliff is for any design. The trade-off between yield and robustness is also demonstrated.

Physics - Plasma physics↗
DOE OSTI2025

The Role of Asparagine as a Gatekeeper Residue in the Selective Binding of Rare Earth Elements by Lanthanide‐Binding Peptides

Abstract Lanthanide‐binding tag (LBT) peptides selectively complex lanthanide cations (Ln 3+ ) in their binding pockets and are promising for lanthanide separation. However, designing LBTs that selectively target specific Ln 3+ cations remains a challenge due to limited molecular‐level understanding and control of interactions within the lanthanide‐binding pocket. In this study, we reveal that the N5 asparagine residue acts as a gatekeeper in the binding pocket, resulting in a 100‐fold selectivity for smaller Lu 3+ over larger La 3+ cations. Nuclear magnetic resonance spectroscopy and molecular dynamics simulations show that the N5 residue weakly binds to the larger La 3+ cation, permitting H 2 O molecules inside the pocket. For the smaller Lu 3+ cations, the N5 residue forms an inter‐arm hydrogen bond with the E14 glutamic acid residue, locking the Lu 3+ cation in the pocket and preventing H 2 O infiltration. Mutating the N5 asparagine to a D5 aspartic acid prevents such a hydrogen bond, eliminating the gatekeeping mechanism and precipitously reducing selectivity. The resulting binding affinity to Ln 3+ cations is non‐monotonic but generally increases with cation size. These results suggest a molecular design paradigm: the reduced affinity for larger lanthanides is due to open pocket conformations, while the selectivity of smaller Ln 3+ cations over larger ones is due to the gatekeeping hydrogen bond.

Chemistry↗
DOE OSTI2023

Review of Carbon Support Coordination Environments for Single Metal Atom Electrocatalysts (SACS)

This topical review focuses on the distinct role of carbon support coordination environment of single-atom catalysts (SACs) for electrocatalysis. The article begins with an overview of atomic coordination configurations in SACs, including a discussion of the advanced characterization techniques and simulation used for understanding the active sites. A summary of key electrocatalysis applications is then provided. These processes are oxygen reduction reaction (ORR), oxygen evolution reaction (OER), hydrogen evolution reaction (HER), nitrogen reduction reaction (NRR), and carbon dioxide reduction reaction (CO 2 RR). The review then shifts to modulation of the metal atom-carbon coordination environments, focusing on nitrogen and other non-metal coordination through modulation at the first coordination shell and modulation in the second and higher coordination shells. Representative case studies are provided, starting with the classic four-nitrogen-coordinated single metal atom (M$-$N 4 ) based SACs. Bimetallic coordination models including homo-paired and hetero-paired active sites are also discussed, being categorized as emerging approaches. The theme of the discussions is the correlation between synthesis methods for selective doping, the carbon structure–electron configuration changes associated with the doping, the analytical techniques used to ascertain these changes, and the resultant electrocatalysis performance. In conclusion, critical unanswered questions as well as promising underexplored research directions are identified.

36 MATERIALS SCIENCE↗