Search NASA⌕ Search

SEARCH · Search NASA

Results for “Full coupled”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Enabling fast-charging of lithium-ion batteries through printed electrodes

It has been well recognized that introducing secondary porous networks (SPNs) into the electrodes can effectively improve the electrochemical performance of lithium-ion batteries (LIBs), especially under fast-charging operations. However, the process complexity and high cost limit the commercial success of advanced electrodes with SPNs. To address this issue, we developed a facile screen-printing process to produce structured graphite electrodes with SPNs. The experimental results demonstrated that, by tuning the diameter and center-to-center (C2C) distance of emulsion dots on the stencil screen, the pore diameters and C2C pore distances of SPNs in screenprinted electrodes can be precisely controlled in the range of 100 mu m to 1 mm and 100 mu m to 3 mm respectively. In addition, the SPNs with hexagonal and square-shape pore alignments have also been imprinted onto the electrode coatings through adjusting the patterns of screen stencils. Used as anodes, the printed graphite electrodes demonstrated significantly reduced overpotential and voltage fluctuation under fast-charging operations from 2C to 6C. Coupled with LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) cathodes, the full cells with printed graphite anodes exhibited an unprecedently stable performance with almost no capacity decay up to 170 cycles when charged to 80 % SOC at 2C. Observations from electron microscopy showed plated lithium undetectable at the surface of printed graphite electrodes after numerous cycles. The electrochemical analysis on the voltage evolution during the cell rest period indicated the significantly delayed onset of lithium plating in the presence of printed graphite electrodes. In conclusion, all these results suggest that the significantly improved cell performance is associated with the shortened Li-ion diffusion distance, reduced polarization and suppressed Li plating in the printed electrodes with patterned SPNs.

25 ENERGY STORAGE↗
DOE OSTI2022

Atomic structure considerations for the low-temperature opacity of Xenon

Here, we have begun a preliminary investigation into the opacity of Xe at low temperatures (<50 eV). The emissivity and opacity of Xe is a crucial factor in determining the utility of Xe in EUV lithography, with numerous industrial applications. To this end, we have been exploring the accuracy of some approximations used in opacity models for the relevant ion stages of Xe (6 times ionized through 20 times ionized). Due to the complexity of Xe atomic structure, one needs to use full configuration-interaction to properly describe the strong mixing in the various n=4 sub-shells that give rise to the Δn=0 and Δn=1 transitions that dominate the opacity spectrum at low temperatures. Since calculations that include full configuration-interaction for large numbers of configurations quickly become computationally prohibitive, we have explored hybrid calculations, in which full configuration-interaction is retained for the most important transitions, while intermediate coupling is employed for all other transitions. After verifying atomic structure properties, local-thermodynamic-equilibrium (LTE) opacities are generated using the ATOMIC code at selected temperatures and densities and compared to various experimental results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Ab Initio Polariton Transport Dynamics with the Classical Path Approximation

We present an ab initio framework for simulating polariton transport dynamics based on the classical path approximation (CPA). The quantum dynamics of polariton transport involves simulating many electronic degrees of freedom, making a fully ab initio dynamics simulation computationally expensive. We demonstrate that the CPA, which removes the need for excited-state nuclear gradients, is well-suited for polaritonic systems because collective light–matter coupling leads to vanishing excited-state forces. Benchmark comparisons between CPA and full evaluation of the excited-state forces show excellent agreement for polariton transport results in model light–matter systems such as polariton group velocities and mean-squared displacements. Ab initio simulations of polariton transport using CPA reproduce key physical trends that are observed in experiments with BODIPY molecules. Our work establishes the CPA as a highly efficient tool for ab initio investigations of transport and energy flow in hybrid light–matter systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

High-throughput screening of 2D materials identifies p-type monolayer WS2 as potential ultra-high mobility semiconductor

Abstract 2D semiconductors offer a promising pathway to replace silicon in next-generation electronics. Among their many advantages, 2D materials possess atomically-sharp surfaces and enable scaling the channel thickness down to the monolayer limit. However, these materials exhibit comparatively lower charge carrier mobility and higher contact resistance than 3D semiconductors, making it challenging to realize high-performance devices at scale. In this work, we search for high-mobility 2D materials by combining a high-throughput screening strategy with state-of-the-art calculations based on the ab initio Boltzmann transport equation. Our analysis singles out a known transition metal dichalcogenide, monolayer WS 2 , as the most promising 2D semiconductor, with the potential to reach ultra-high room-temperature hole mobilities in excess of 1300 cm 2 /Vs should Ohmic contacts and low defect densities be achieved. Our work also highlights the importance of performing full-blown ab initio transport calculations to achieve predictive accuracy, including spin–orbital couplings, quasiparticle corrections, dipole and quadrupole long-range electron–phonon interactions, as well as scattering by point defects and extended defects.

Chemistry↗
DOE OSTI2020

Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys

We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials’ design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum–niobium–titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.

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

Non-equilibrium rate theory for polariton relaxation dynamics

We derive an analytic expression of the non-equilibrium Fermi’s golden rule (NE-FGR) expression for a Holstein–Tavis–Cumming Hamiltonian, a universal model for many molecules collectively coupled to the optical cavity. These NE-FGR expressions capture the full-time-dependent behavior of the rate constant for transitions from polariton states to dark states. The rate is shown to be reduced to the well-known frequency domain-based equilibrium Fermi’s golden rule (E-FGR) expression in the equilibrium and collective limit and is shown to retain the same scaling with the number of sites in non-equilibrium and non-collective cases. We use these NE-FGR to perform population dynamics with a time-non-local and time-local quantum master equation and obtain accurate population dynamics from the initially occupied upper or lower polariton states. Furthermore, NE-FGR significantly improves the accuracy of the population dynamics when starting from the lower polariton compared to the E-FGR theory, highlighting the importance of the non-Markovian behavior and the short-time transient behavior in the transition rate constant.

Chemical dynamics↗
DOE OSTI2025

Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode

It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025

Laboratory plasma↗
DOE OSTI2023

Electric Vehicle Infrastructure Consequence Assessment

With consumers’ growing interest in electric vehicles, extreme fast charging stations are poised to provide high-power charging to rapidly recharge light-duty passenger vehicles. High-power charging requires high-level communication between vehicle and charger to govern the charging process. The coupling of power and communication increases the potential scale of cyberattacks. Using a full Western Electricity Coordinating Council planning model, load manipulation from high-power charging infrastructure is investigated. Two cases of load manipulation are studied: (i) a discrete, widespread system event and (ii) loads modulated near the Western Interconnect’s resonant frequency. In (i) some generation trips and in (ii) oscillations are observed on the California Oregon Intertie. Neither scenario results in significant adverse effects to the grid.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗
DOE OSTI2026

Auxiliary-state-facilitated phase synchronization phenomena in isolated spin systems

Extending classical synchronization to the quantum domain is of great interest both from the fundamental physics point of view and with a view toward quantum technology applications. This work characterizes phase synchronization of an effective spin-1 system, which is realized by coupling three quantum states with infinite lifetime to auxiliary excited states that have a finite lifetime. Integrating out the excited states, the effective spin-1 model features coherent and incoherent effective couplings. The following are our key findings. (i) Phase synchronization can be controlled by adjusting the phases of the couplings to the excited states. (ii) Unlike in the paradigmatic spin-1 system studied in the literature, where the dissipative couplings describe decay into the limit-cycle state, the effective spin-1 model investigated in this work is governed by a competition between dissipative decay into and out of the limit-cycle state, with the dissipative decay out of the limit-cycle state playing a critical role. (iii) We identify a parameter regime where phase synchronization of the effective spin-1 system is, in the absence of coherent effective couplings, governed entirely by the effective dissipators. The effective spin-1 model is benchmarked through comparisons with master-equation calculations for the full Hilbert space. Physical insights are gained through analytical perturbation theory calculations. In conclusion, our findings, which are expected to hold for a broad class of energy-level and coupling schemes, are examined using hyperfine states of 87 Rb as an example system.

Nonlinear optics↗
DOE OSTI2020

Relativistic two-component projection-based quantum embedding for open-shell systems

In this work, we present a relativistic quantum embedding formalism capable of variationally treating relativistic effects, including scalar-relativity and spin–orbit coupling. We extend density functional theory (DFT)-in-DFT projection-based quantum embedding to a relativistic two-component formalism, where the full spin magnetization vector form is retained throughout the embedding treatment. To benchmark various relativistic embedding schemes, spin–orbit splitting of the nominally t 2g valence manifold of W(CO) 6 , exchange coupling of [(H 3 N) 4 Cr(OH) 2 Cr(NH 3 ) 4 ] 4+ , and the dissociation potential curve of WF 6 are investigated. The relativistic embedding formalism introduced in this work is well suited for efficient modeling of open-shell systems containing late transition metal, lanthanide, and actinide molecular complexes.

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

A Monte-Carlo method for Coulomb collisions in moderately coupled plasmas

Here, we introduce a new approach for coupling cumulative and single-scattering distributions, extending Monte Carlo methods for Coulomb collisions in weakly coupled plasmas into the moderately coupled regime. Using time steps typical of cumulative-scattering methods, we show that this approach not only recovers leading-order corrections proportional to the inverse Coulomb logarithm but also solves the full Boltzmann collision integral. A key feature of the method is its simplicity: it can be integrated into cumulative-scattering methods for weakly coupled plasmas with only ~ 10 additional lines of code.

Boltzmann collision integral↗
DOE OSTI2023

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Explainable machine learning reveals that local structural motifs encode the thermodynamic state across the CuZr metallic glass-forming range

Metallic glasses derive their properties from the statistics of local atomic motifs rather than from long-range order, yet a quantitative, chemistry-specific link between motif populations and the underlying glassy state has remained elusive. In this work we combine large-scale molecular dynamics, Voronoi tessellation, deep neural networks, and SHapley Additive exPlanations (SHAP) to identify which local structural motifs define the glassy state of Cu—Zr metallic glasses. A dataset of 17,180 atomistic configurations spanning ten compositions (Cu 20 Zr 80 –Cu 80 Zr 20 ) and four quench rates (10 9 –10 12 K/s) is used to train a feed-forward neural network that regresses temperature across the 50–2000 K liquid–supercooled–glass range, achieving a mean absolute error of 19.89 K and R 2 = 0.9974, confirming that the local structural state is faithfully encoded in motif-level structure. SHAP analysis then reveals that a tightly coupled near-icosahedral family of motifs (coordination numbers (CN) 11–13, including the full icosahedron 001200 and its single-atom-perturbation sibling 10930) collectively encodes the thermodynamic state of the system across the full glass-forming range. The CN = 11–13 ordered members carry negative SHAP values at high populations, tracking the most deeply-quenched configurations, while 10930 shows the reversed signature consistent with its role as a soft-spot host whose population shrinks as the icosahedral network deepens. The analysis demonstrates that explainable machine learning can isolate the minimal motif vocabulary defining the glassy state and recovers the near-icosahedral building blocks previously identified by data-driven analyses of Cu—Zr. The approach provides a general, chemistry-specific route for characterizing the structural state of disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unraveling the conundrum of electronic leakage in protonic ceramic cells: Operation-specific insights and rational design strategies

Electronic conduction through proton-conducting electrolytes significantly impairs the efficiency of protonic ceramic cells (PCCs). Here, in this study, we explore the electron and ion mixed transport properties of four common protonic ceramics, BaZr 0.8 Y 0.2 O 3-δ (BZY82), BaZr 0.7 Ce 0.2 Y 0.1 O 3-δ (BZCY721), BaZr 04 Ce 0.4 Y 0.1 Yb 0.1 O 3-δ (BZCYYb4411), and BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ (BZCYYb1711). It marks the first instance of investigating these properties under operation-specific scenarios: fuel side of electrolysis cell, air side of electrolysis cell, fuel side of fuel cell, and air side of the fuel cell. BZCYYb1711 exhibits the highest ionic conductivity, but two to three times higher electronic leakage when exposed to oxygen-containing environments than the others. BZY82 exhibits approximately two times higher electronic leakage in a hydrogen-containing environment. BZCY721 demonstrates excellent ion transport numbers (~0.95) across these four operating conditions. BZCYYb4411 behaves quite similarly to BZCY721. The most challenging operating environment for all candidates is the air side of fuel cell mode. This mode leads to a high initial electronic leakage, followed by a significant increase with polarization. The probable cause for this behavior is a H 2 -free, polarization-induced reduction that leads to the formation of V$^•_O$. The electron small polaron associated with V$^•_O$ is released by the electrical field due to the Poole-Frenkel effect. ZnO and NiO sintering aids are found to be detrimental to the ionic conductivity of the electrolytes. In particular, NiO substantially lowers the ion transport number. The correlation of the operation-specific electronic leakage to full cells is discussed. It is suggested that a rational PCC design should synergistically couple BZCYYb1711 at fuel side with BZCY4411 at air side to deliver a well-balanced performance and faradaic efficiency simultaneously, and the high temperature sintering process with a NiO fuel electrode should be shortened or replaced by ultra-fast sintering techniques or using a fuel electrode scaffold-infiltration fabrication strategy.

42 ENGINEERING↗
DOE OSTI2022

Comparison of spherical harmonics method and discrete ordinates method for radiative transfer in a turbulent jet flame

Here, in this study, we systematically compared the accuracy and computational cost of two popular solution methods for the radiative transfer equation (RTE): the spherical harmonics method (P N ) and the discrete ordinates method (DOM). We first investigated convergence characteristics of different orders of P N and DOM in a series of 1D homogeneous configurations with varying optical thicknesses. Both solvers perform better for optically thicker cases. The accuracy of P N methods increases with its order, , but the gain in accuracy reduces with the increase in , i.e., improvement of P 7 over P 5 is less than that of P 3 over P 1 . This decreasing trend becomes more prominent as the optical thickness decreases. On the other hand, DOM’s accuracy increases almost linearly with the increase in the number of ordinates (or polar angles in this study) in all cases. While comparing the directional profile of radiative intensity, both solvers perform better when the radiative intensity is more isotropic. These solvers were then connected with a full spectrum k-distribution (FSK) spectral model and used to perform radiation-coupled simulations of a turbulent jet flame in an axi-symmetric cylindrical domain. Results obtained from P 1 to P 7 approximations for P N , and 2 x 4, 4 x 4, 4 x 8, 8 x 8 finite angles for DOM are compared with that from an optically thin model, and a reference solution from line-by-line (LBL) photon Monte Carlo (PMC) method. The choice of radiation solver shows a noticeable impact on the temperature distribution of the flame. The P N solvers lead to slightly higher radiant fractions and the DOM solvers lead to slightly lower radiant fractions than the PMC benchmark solution. Finally, the computational costs of each of these solvers are also reported and an intermittent evaluation / time blending scheme to improve the computational efficiency of radiation solvers in radiation-coupled simulations are also demonstrated.

42 ENGINEERING↗
DOE OSTI2021

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗