Search NASA⌕ Search

SEARCH · Search NASA

Results for “function approximation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beam Signal Recovery Using Wiener Deconvolution

This project explores the use of Wiener deconvolution to recover an original signal that has been distorted by a known transfer function and by noise. A simulated Gaussian pulse was used as the test signal, and a transfer function was applied in the frequency domain to model system distortion. Controlled noise was then introduced to approximate real-world signal degradation. A Wiener filter was implemented to reverse the effects of the transfer function while minimizing the influence of noise. The recovered signal was compared to the original pulse to evaluate the effectiveness of the filter. Results demonstrate that Wiener deconvolution offers a stable and effective approach to signal recovery, balancing complete transfer function inversion and noise suppression.

Campbell, Savanna [El Camino Coll.]↗
DOE OSTI2024

Electronic Properties of Ultra‐Wide Bandgap B x Al 1− x N Computed from First‐Principles Simulations

Abstract Ultra‐wide bandgap (UWBG) materials such as AlN and BN hold great promise for future power electronics due to their exceptional properties. They exhibit large bandgaps, high breakdown fields, high thermal conductivity, and high mechanical strengths. AlN and BN have been extensively researched, however, their alloys, B x Al 1− x N, are much less studied despite their ability to offer tunable properties by adjusting x . In this article, the electronic properties of 17 recently predicted ground states of B x Al 1− x N in the x = 0 − 1 range are predicted using first‐principles density functional theory and many‐body perturbation theory within GW approximation. All the B x Al 1− x N structures are found to be UWBG materials and have bandgaps that vary linearly from that of wurtzite‐phase ( w ) AlN (6.19 eV) to that of w ‐BN (7.47 eV). The bandstructures of B x Al 1− x N show that a direct‐to‐indirect bandgap crossover occurs near x = 0.25. Furthermore, it is found that B x Al 1− x N alloys have much larger dielectric constants than the constituent bulk materials (AlN = 9.3 ɛ 0 or BN = 7.3 ɛ 0 ), with values reaching as high as 12.1 ɛ 0 . These alloys are found to exhibit large dielectric breakdown fields in the range 9–35 MV cm −1 with a linear dependence on x . This work provides the much needed advancement in the understanding of the properties of B x Al 1− x N to aid their application in next‐generation devices.

Milne, Cody L.↗
DOE OSTI2025

Sm 2 Ru 3 Sn 5 : A Noncentrosymmetric Cubic Member of the Ln 2 M 3 X 5 Family

An optimized synthetic method is presented for Sm 2 Ru 3 Sn 5 and investigate its physical properties and electronic structure. Sm 2 Ru 3 Sn 5 is prepared by arc-melting stoichiometric ratios of the elements and is confirmed by single crystal and powder X-ray diffraction. An antiferromagnetic transition is observed at T N = 3.8 K. A modified Curie-Weiss fit to the data in the range 50–150 K yields a Curie-Weiss temperature: θ CW = −36.6 K and an effective magnetic moment: μ eff = 0.83 μ B , in agreement with a Sm 3+ oxidation state. Field-dependent magnetization up to H = 7 T at 2 K shows a maximum response of 0.06 μ B , which is significantly lower than the expected Sm 3+ saturation moment (0.71 μ B ). Resistivity measurements indicate metallic behavior, and analysis of the magnetic entropy from the heat capacity reveals a doublet ground state due to crystal electric field splitting. The electronic structure and density of states are calculated with density function theory and further supported by the local density approximation with dynamical mean-field theory. Finally, the experimental and computational results highlight localized Sm 3+ moments and suggest a possible interplay between Ruddelman–Kitel–Kasuya–Yosida and Kondo interactions, positioning Sm 2 Ru 3 Sn 5 as a promising material for studying topology and complex physical phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Orientation-dependent structural properties during growth and growth mechanism of CoO films

The orientation-dependent local structural properties of CoO films on sapphire substrates during growth were investigated through linearly-polarized extended X-ray absorption fine structure (EXAFS) measurements. Specifically, CoO(1 1 1) and (100) crystals with a rock-salt structure (Fm3m) were epitaxially grown on alpha-Al2O3(0001) and (1012) substrates, respectively, at 700 degrees C using a radio-frequency sputtering system. The local structural properties of CoO films in the in-plane and out-of-plane orientations were quantitatively determined using linearly-polarized EXAFS at the Co K-edge during growth. The EXAFS analysis revealed that during the initial stages of growth, the local structural properties exhibit significant differences compared to thick films, with short atomic distances and large (small) Debye-Waller factors observed in the out-of-plane (in-plane) orientations. The local structural strain mostly diminished when approximately 20 CoO layers accumulated on the substrate. Density functional theory (DFT) calculations further supported these findings by confirming that cobalt atoms initially form stable bonds with the sapphire surface, leading to the simultaneous growth of oxygen and cobalt layers in a coordinated manner through layer-by-layer growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗
DOE OSTI2024

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

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

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗
DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Wiener Filter Deconvolution for Analog Signals

This project explores the use of Wiener deconvolution to recover an original signal that has been distorted by a known transfer function and by noise. A simulated Gaussian pulse was used as the test signal, and a transfer function was applied in the frequency domain to model the system distortion. Controlled noise was then introduced to approximate real-world signal degradation. A Wiener filter was implemented to reverse the effects of the transfer function while minimizing the influence of noise. The recovered signal was compared with the original pulse to evaluate the effectiveness of the filter. The results demonstrate that Wiener deconvolution offers a stable and effective approach to signal recovery, balancing complete transfer function inversion and noise suppression.

Campbell, Savanna [Fermilab]↗
DOE OSTI2025

The Langdon effect in laser plasmas: Absorption and conduction

A plasma heated by inverse bremsstrahlung absorption of laser light develops a non-Maxwellian electron distribution function, called the Langdon effect [A. B. Langdon, Phys. Rev. Lett. 44, 575 (1980)]. These non-Maxwellian distributions are sufficiently long-lived to impact the absorption processes itself as well as the transport of heat by electrons. The theory of the Langdon effect in a homogeneous plasma is reviewed to clarify some aspects of Langdon's derivation as well as to confirm that the widely used super-Gaussian approximation works fairly well to describe the shape of the distribution function and reduction of the absorption rate. The Langdon effect on thermal conduction in an inhomogeneous plasma is developed by considering perturbations in a homogeneous absorbing plasma, which develops a heat flux due to both temperature and density gradients. A practical theory of the heat flux is developed by fitting the results of Vlasov–Fokker–Planck simulations, which avoids several approximations that compromised the usefulness of past theoretical predictions, most critically, the effect of electron–electron collisions on the fluxes. The present fits parameterize the coefficients of the temperature gradient (thermal conductivity) and the density gradient for a plasma of any ionization state and for any laser intensity where the theory of the Langdon effect remains locally valid. It is expected that this generalized theory of heat flow in an absorbing plasma will improve the predictive capability of radiation-hydrodynamics simulations of laser-produced plasmas, especially those formed in inertial confinement fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗
DOE OSTI2025

An Accurate, Affordable Density Functional Tight-Binding Model for Excited State Hydrocarbon Polymer Molecular Dynamics

We have developed a density functional tight-binding model for hydrocarbon excited state dynamics by referencing to high-level electronic structure theory and incorporating a many-body repulsive energy. We then validate our model against n-octane geometry optimizations, bond dissociation scans, and vibrational frequencies. Our model is approximately 10 3 times more efficient than hybrid time-dependent density functional theory calculations with comparable accuracy. In conclusion, our efforts enable longer timescale excited state simulations of photochemistry and scattering of incident radiation.

Chemical calculations↗
DOE OSTI2026

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗
DOE OSTI2024

Self-consistent quasiparticle G W and hybrid functional calculations for Al/InAs/Al heterojunctions: Band offset and spin-orbit coupling effects

The electronic structure of surfaces and interfaces plays a key role in the properties of quantum devices. Here, we study the electronic structure of realistic Al/InAs/Al heterojunctions using a combination of density functional theory with hybrid functionals and state-of-the-art quasiparticle G W (QSGW) calculations. We find a good agreement between QSGW calculations and hybrid functional calculations, which themselves compare favorably well with angle-resolved photoemission spectroscopy experiments. Our paper confirms the need for well-controlled quality of the interfaces to obtain the needed properties of InAs/Al heterojunctions. A detailed analysis of the effects of spin-orbit coupling on the spin splitting of the electronic states shows a linear scaling in k space, related to the two-dimensional nature of some interface states. The good agreement by QSGW and hybrid functional calculations opens the door towards trustable use of an effective approximation to QSGW for studying very large heterojunctions. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cluster spin glass correlations and dynamics in Zn 0.5⁢ Mn 0.5⁢ Te

Here, we present a combined magnetometry, muon spin-relaxation (𝜇⁢SR), and neutron-scattering study of the insulating spin glass Zn 0.5 ⁢Mn 0.5 ⁢Te, for which magnetic Mn 2+ and nonmagnetic Zn 2+ ions are randomly distributed on a face-centered cubic lattice. The magnetometry and 𝜇⁢SR results confirm a spin freezing transition around 𝑇 𝑓 ≈ 23 K, with the spin-fluctuation rate decreasing gradually and somewhat inhomogeneously through the sample volume as the temperature decreases toward 𝑇 𝑓 . Characteristic spin-correlation times well above 𝑇 𝑓 are on the order of 10 −10 s, much slower than typically observed in canonical spin glasses but in line with expectations for a cluster spin glass. Using magnetic pair distribution function (mPDF) analysis and reverse Monte Carlo (RMC) modeling of the magnetic diffuse neutron-scattering data, we show that the spin-glass ground state consists of clusters of spins exhibiting short-range-ordered type-III antiferromagnetic correlations with a locally ordered moment of 3.1⁢(1)⁢𝜇 B between nearest-neighbor spins. The type-III correlations decay exponentially as a function of spin separation distance with a correlation length of approximately 5 Å. The diffuse magnetic scattering and corresponding mPDF show no significant changes across 𝑇 𝑓 , indicating that the dynamically fluctuating short-range spin correlations in the paramagnetic state retain the same basic type-III configuration that characterizes the spin-glass state; the only change apparent from the neutron-scattering data is a gradual reduction of the correlation length and locally ordered moment with increasing temperature. Taken together, these results paint a unique and detailed picture of the local magnetic structure and dynamics in Zn 0.5 ⁢Mn 0.5⁢ Te and provide strong evidence that this material is best described as a cluster spin glass. In addition, this work showcases a statistical method for extracting diffuse scattering signals from neutron powder diffraction data, which we developed to facilitate the mPDF and RMC analysis of the neutron data. This method has the potential to be broadly useful for neutron powder diffraction experiments on a variety of materials with short-range atomic or magnetic order.

magnetism↗
DOE OSTI2025

Effective description of ajar systems with a 𝑈⁡(1) symmetry

We introduce the concept of “ajar systems” as an intermediate case between closed and open systems, where the time scale for charge exchange with the environment is parametrically larger than all other characteristic time scales. The Schwinger-Keldysh effective action for such finite-temperature systems exhibits a symmetry group G 1 × G 2 weakly broken to its diagonal subgroup G diag , which we systematically describe using spurion techniques. For G = U(1), we calculate leading-order corrections to correlation functions in both diffusive and spontaneously broken phases. Unlike systems with approximate symmetries where both real and imaginary parts of dispersion relations receive corrections, ajar systems modify only the imaginary part—preserving gapless modes while adding finite damping.

Effective field theory↗