SEARCH · Search NASA

Results for “high order 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Counterdiabatic Driving with Performance Guarantees

Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Engineering-Scale Validation of Novel Algae CO 2 Capture and Bioproducts Technology

The purpose of this project was to advance and scale the algae-based technology Helios-NRG has been developing for capture of CO 2 from the flue gas of coal-burning power plants with conversion of the CO 2 to high-value products. In order to encourage the process becoming commercial and be adoptable, it is essential that the cost of CO 2 be close to net zero or slightly revenue positive. Two methods to achieve this are reducing the cost of operations through efficiencies and selling of products derived from conversion. Helios and the project team partners have accomplished the goals of furthering the technology to be more scalable and substantially lower cost with a defined pathway for commercialization.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2026

Dissecting Disorder: Defect-Driven Structural Complexity in Layered Li3InCl6 Solid Electrolyte

Halide solid electrolytes have emerged as promising candidates for solid-state batteries owing to their high oxidative stability and ionic conductivity. Among them, Li3InCl6 (LIC) has attracted significant attention. However, diffraction patterns of LIC synthesized via different methods exhibit distinct differences particularly at low-angle reflectionsindicative of underlying structural disorder. These variations are attributed to deviations from ideal crystallographic order, especially stacking faults, whose impact on structure and ion transport remains poorly understood. Here, we identify and quantify stacking faults in LIC samples prepared under different synthetic conditions. Using X-ray diffraction and time-of-flight neutron diffraction, we construct and refine stacking fault models that accurately reproduce the experimental diffraction features. LIC samples with higher degrees of stacking faults exhibit only negligible differences in ionic conductivities and activation energies. This indicates that stacking faults have a limited impact on altering the Li+ diffusion pathway along the c-axis, likely due to the high concentration of vacancies in the In layers, while Li+ diffusion remains nearly unchanged in the ab-plane. Our results account for the observed differences in diffraction patterns across samples and provide a quantitative assessment of faulting probabilities and stacking sequences. The insights gained from this study are expected to be broadly applicable to other layered halide solid electrolytes and contribute to a deeper understanding of the role of structural disorder in ion transport.

Liu, Jue [ORNL] (ORCID:000000024453910X)↗

DOE OSTI2024

Solution processed high aspect ratio ultra-long vertically well-aligned ZnO nano scintillators for potential X-ray imaging applications

We report the photon (PL), electron (CL) and X-ray (XEL) induced luminescence characteristics of high aspect ratio ultra-long (~ 50 µm) ZnO nanorods (NRs) and discuss the potential for fast X-ray detection based on the consistent and efficient visible emission (~ 580 nm) from ZnO NRs. Nanostructured ZnO scintillators were rearranged to form a vertically well-aligned NR design in order to help light absorption and coupling resulting in luminescent and fast scintillation properties. The design of the nanorod array combines the key advantages of a low-cost growth technique together with environmentally friendly and widely available materials. A low temperature hydrothermal method was adopted to grow ZnO NRs in one cycle growth and their structural, optical and X-ray scintillation properties were investigated. The relatively short (~ 10 µm) ZnO NRs emitting in the near-band-edge region were found to be almost insensitive to X-rays. On the other hand, the higher XEL response of long ZnO NRs, which is a key parameter for evaluation of materials to be used as scintillators for high quality X-ray detection and imaging, along with a decay time response in the order of ns confirmed promising scintillation properties for fast and high-resolution X-ray detector applications.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Quantum fluctuations in dense plasma simulations

Molecular dynamics (MD) simulations are a powerful tool for modeling warm and hot dense matter. Density functional theory (DFT) MD simulations are often preferred in dense plasmas in order to accurately model quantum electronic structure. However, DFT-MD simulations neglect interaction effects due to fluctuations in excited states. In this work, we present an MD approach that uses excited state method pseudoatoms to run dense plasma simulations with many different core-electron configurations at classical MD speeds. We also allow for transitions between different configurations in our simulations and find that these fluctuations are especially important for highly excited states. Our results suggest that finite configuration lifetimes that are comparable to the inverse ion plasma frequency need to be accounted for in order to accurately model ion distributions in dense plasma simulations. We also demonstrate that excited state fluctuations have a direct impact on ion plasma microfields, generate different plasma microfields for different excitation levels, and thereby induce absorption–emission line shape asymmetries even in steady-state plasmas.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Effect of Plutonium on Uranium Oxide Microstructural Fingerprint

The analysis of particulates from environmental sampling is routinely performed for nuclear forensics applications. In order to increase the tools available for nuclear forensics applications, capability development materials (CDMs), or reference particulates, are necessary for developing and benchmarking new analytical methods. Historically, these CDM particulates have been produced with highly controlled and characterized isotopic and size parameters for applications in developing and benchmarking particle sizers and mass spectrometry analytical methods [1- 4]. However, the development of CDMs with highly characterized particle morphology and phase of the particles is vital for aiding in the development and benchmarking of particle analytical methods for those parameters. In many cases, key properties such as crystallographic phase, morphology, and microstructure, may be correlated to processing history of environmental sampling particulates [5]. To that end, this work investigates determining the structural fingerprint of produced Pu-doped uranium oxide CDMs using transmission electron microscopy (TEM). The particulates, one of which shown in Figure 1, were synthesized to Fig. 1. Single particulate of uranium oxide fabricated using the THESEUS technique. develop capabilities for environmental sampling investigations. These particles are monodisperse at diameter of 1 μm and a range of plutonium concentrations from 0- 1,000 ppm. Given their small diameter, TEM analysis was ideal for studying the particulate structural fingerprint in detail. Previous investigations confirmed the external homogeneity of the particulates at different plutonium concentrations, however this analysis focuses on determining the grain structure, phase distribution, and porosity of the particulates.

Mayer, Jack [Univ. of Florida, Gainesville, FL (Un↗

DOE OSTI2024

Investigation of Numerical Methods for Performance Improvement of MOOSE-based System Analysis Codes

The main objective of this study is to investigate the feasibility of implementation of staggered-grid finite volume method (SG-FVM) using the MOOSE framework to support the development of the advanced system analysis code SAM. This study successfully demonstrated the integration of staggered-grid finite volume method in the application level using the MOOSE framework, although not in the framework level which should be investigated in the future. Several important properties of the implemented SG-FVM, e.g., high-order spatial accuracy and monotonicity preserving, have been demonstrated with selected numerical test cases. The superior performance in execution time was also evident based on a selected one- dimensional flow problem in a loop configuration.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

An efficient level set method for tracking many materials

Here, we present an efficient level set method to track an arbitrary number of materials. The algorithm is optimal in the sense that it only needs to store a single unsigned distance-like function and a single integer indicator function, independent of the number of materials or distinct regions being tracked. Furthermore, for smooth velocity fields and smooth interface shape, arbitrarily high order solutions can be demonstrated. For interfaces that are or become kinked, the solution is limited to second-order convergence rates in the L 1 norm and first-order in the L ∞ norm.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Third-order photon correlations extract single-nanocrystal multiexciton properties in solution

Colloidal semiconductor nanocrystals are considered promising materials for high-flux optical applications, including lasing, light-emitting diodes, biological imaging, and quantum optics. In high-flux applications, multiexcitons can significantly contribute to emission, influencing its brightness, spectral purity, and kinetics. As a result, understanding and controlling multiexciton emission in colloidal nanocrystal materials is of the utmost importance. In the past, single-nanocrystal photon correlation methods have been applied to understand biexciton and triexciton efficiencies, lifetimes, and spectra. While powerful, such methods suffer from user selection bias and require stable emission from single nanocrystals. To compensate for this shortcoming, second-order correlation methods were developed to extract sample-averaged biexciton properties from a solution of nanocrystals. Until now, however, the analogous third-order solution photon correlation methods remained unexplored. In this work, we present a pair of third-order photon correlation techniques to obtain the sample-averaged single-nanocrystal triexciton quantum yield and lifetime in a solution-phase experiment. These techniques derive from the relationship between the Poisson probability of nanocrystal photon absorption and the intrinsic probability of nanocrystal photon emission. We validate the theoretical background of these techniques by creating a numerical model to simulate the diffusion and emission of many nanocrystals in solution. Our simulations confirm that the average triexciton quantum yield and triexciton lifetime can be extracted from a solution of nanocrystals. These techniques will enable researchers to gain a better understanding of the fundamental multiexciton properties of colloidal nanocrystals.

Horowitz, Jonah R. [Massachusetts Institute of Tec↗

DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗

DOE OSTI2026

Understanding polaronic transport in complex oxides by combining precise synthesis and first-principles many-body theory

In complex oxides, charge carriers often couple strongly with lattice vibrations to form polarons–entangled electron–phonon quasiparticles whose transport properties remain difficult to characterize. Experimental access to intrinsic polaronic transport requires ultraclean samples, while theoretical description demands methods beyond low-order perturbation theory. Here, we show a predictive theory–experiment workflow to study polaron transport in complex oxides. Focusing on a prototypical polaronic oxide, anatase TiO 2 , we combine growth of high-quality oxygen-vacancy-doped films using hybrid molecular beam epitaxy with a first-principles electron–phonon diagrammatic Monte-Carlo (FEP-DMC) framework recently developed for accurate polaron predictions. Our films exhibit record-high electron mobility for anatase TiO 2 , in excellent agreement with FEP-DMC calculations conducted prior to experiment, which predict a room-temperature mobility of 45 ± 15 cm −2 V −1 s −1 and a mobility-temperature scaling of μ ∝ T −1.9 ± 0.077 . Microscopic analysis using scanning transmission electron microscopy and x-ray photoelectron spectroscopy reveals the role of oxygen vacancies in modulating transport at lower temperatures. FEP-DMC further provides quantitative insight into polaron formation energy, phonon cloud distribution, lattice distortion around the polaron, and the polaronic contribution to mobility. Together, these results provide a deeper microscopic understanding of large-polaron transport in a complex oxide and provide the blueprint to characterize other polaronic materials.

anatase TiO2↗

DOE OSTI2025

Quantum error mitigation for Fourier moment computation

Hamiltonian moments in Fourier space—expectation values of the unitary evolution operator under a Hamiltonian at different times—provide a convenient framework to understand quantum systems. They offer insights into the energy distribution, higher-order dynamics, response functions, correlation information, and physical properties. This paper focuses on the computation of Fourier moments within the context of a nuclear effective field theory on superconducting quantum hardware. The study integrates echo verification and noise renormalization into Hadamard tests using control reversal gates. These techniques, combined with purification and error suppression methods, effectively address quantum hardware decoherence. The analysis, conducted using noise models, reveals a significant reduction in noise strength by two orders of magnitude. Moreover, quantum circuits involving up to 266 gates over five qubits demonstrate high accuracy under these methodologies when run on IBM superconducting quantum devices. Published by the American Physical Society 2025

Kiss, Oriel (ORCID:0000000174613342)↗

DOE OSTI2025

An evaluation of multi-fidelity methods for quantifying uncertainty in projections of ice-sheet mass change

Abstract. This study investigated the computational benefits of using multi-fidelity statistical estimation (MFSE) algorithms to quantify uncertainty in the mass change of Humboldt Glacier, Greenland, between 2007 and 2100 using a single climate change scenario. The goal of this study was to determine whether MFSE can use multiple models of varying cost and accuracy to reduce the computational cost of estimating the mean and variance of the projected mass change of a glacier. The problem size and complexity were chosen to reflect the challenges posed by future continental-scale studies while still facilitating a computationally feasible investigation of MFSE methods. When quantifying uncertainty introduced by a high-dimensional parameterization of the basal friction field, MFSE was able to reduce the mean-squared error in the estimates of the statistics by well over an order of magnitude when compared to a single-fidelity approach that only used the highest-fidelity model. This significant reduction in computational cost was achieved despite the low-fidelity models used being incapable of capturing the local features of the ice-flow fields predicted by the high-fidelity model. The MFSE algorithms were able to effectively leverage the high correlation between each model's predictions of mass change, which all responded similarly to perturbations in the model inputs. Consequently, our results suggest that MFSE could be highly useful for reducing the cost of computing continental-scale probabilistic projections of sea-level rise due to ice-sheet mass change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Conformal BK equation at QCD Wilson-Fisher point

High-energy scattering in pQCD in the Regge limit is described by the evolution of Wilson lines governed by the BK equation. In the leading order, the BK equation is conformally invariant and the eigenfunctions of the linearized BFKL equation are powers. It is a common belief that at d ≠ 4 the BFKL equation is useless since unlike d = 4 case it cannot be solved by usual methods. However, we demonstrate that at critical Wilson-Fisher point of QCD the relevant part of NLO BK restores the conformal invariance so the solutions are again powers. As a check of our approach to high-energy amplitudes at the Wilson-Fisher point, we calculate the anomalous dimensions of twist-2 light-ray operators in the Regge limit j → 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

On optimal control of hybrid dynamical systems using complementarity constraints

Optimal control for switch-based dynamical systems is a challenging problem in the process control literature. In this study, we model these systems as hybrid dynamical systems with finite number of unknown switching points and reformulate them using non-smooth and non-convex complementarity constraints as a mathematical program with complementarity constraints (MPCC). We utilize a moving finite element based strategy to discretize the differential equation system to accurately locate the unknown switching points at the finite element boundary and achieve high-order accuracy at intermediate non-collocation points. We propose a globalization approach to solve the discretized MPCC problem using a mixed NLP/MILP-based strategy to converge to a non-spurious first-order optimal solution. The method is tested on three dynamic optimization examples, including a gas–liquid tank model and an optimal control problem with a sliding mode solution.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

FLAMES─Fast, Low-Storage, Accurate, and Memory-Efficient Adaptive Sampling─Approach to Resolve Spatially Dependent Dynamics of Molecular Liquids

Many critical phenomena in soft matter occur at large length scales, necessitating the resolution of their structure and dynamics at low wavenumbers. However, resolving wavenumber-dependent dynamics computationally via molecular dynamics simulations presents significant challenges, as these phenomena span several orders of magnitude in both time and length scales, resulting in high computational costs and memory demands. Here, this work highlights the computational and memory challenges associated with analyzing molecular trajectories in reciprocal space and demonstrates a method to address them. We introduce FLAMESFast, Low-storage, Accurate, and Memory-Efficient adaptive Sampling, which is a direct method for calculation of structure factors, allowing us to select only the required number of wavevectors for binning. We also use wavenumber-dependent time steps to extract dynamics. Our FLAMES approach effectively mitigates computational and memory/storage bottlenecks. We demonstrate the method using simulations of a model system, liquid octane, at various temperatures. Comparisons with experimental data and real space computation show that the FLAMES technique achieves high accuracy in resolving temperature- and spatially dependent dynamics while being significantly more computationally efficient and requiring less memory and storage than methods based on a uniform wavevector grid and fixed temporal spacing.

Chen, Guang [Argonne National Laboratory (ANL), Ar↗

DOE OSTI2025

Benchmarking third-order cluster perturbation theory for electronically excited states

In this study, we investigate the reliability of cluster perturbation (CP) theory applied to the calculation of electronically excited states through a comprehensive benchmark. In CP theory, perturbative corrections are added to the properties of a parent excitation space, which converge toward the properties of a target excitation space. For the CPS(D-n) model, perturbative corrections through order n are added to the coupled cluster singles (CCS) excitation energies to target the coupled cluster singles and doubles (CCSD) excitation energies. Through a comparative analysis of excitation energy calculations across a diverse set of molecules and wavefunction methods, we present a comprehensive evaluation of the accuracy of the third-order CPS(D) model, CPS(D-3), in calculating excitation energies. Further, our findings demonstrate that CPS(D-3) is a reliable alternative to established methods, particularly CCSD, while systematically overestimating the excitation energies compared to high-level coupled cluster methods such as CC3. These results highlight the strengths and limitations of CPS(D-3), as well as the promising directions for its future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗