Search NASA⌕ Search

SEARCH · Search NASA

Results for “functional 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Approximating the particle distribution in rotating and tandem mirror traps

Steady-state distribution functions can be used to calculate stability conditions for modes, radiation energy losses and particle loss rates. Heuristic analytic approximations to these distributions can capture key behaviors of the true distributions such as the relative speeds of different transport processes while possessing computational advantages over their numerical counterparts. In this paper, we motivate and present a closed-form ana- lytic model for a distribution of particles in a centrifugal or tandem mirror. We find that our model outperforms other known models in approximating numerical steady- state simulations outside of a narrow range of low confining potentials. We demonstrate the model’s suitability in the high confining potential regime for applications such as loss-cone stability thresholds, fusion yields and available energy.

Li, G.X.↗
DOE OSTI2025

Approximating the particle distribution in rotating and tandem mirror traps

Steady-state distribution functions can be used to calculate stability conditions for modes, radiation energy losses and particle loss rates. Heuristic analytic approximations to these distributions can capture key behaviors of the true distributions such as the relative speeds of different transport processes while possessing computational advantages over their numerical counterparts. In this paper, we motivate and present a closed-form analytic model for a distribution of particles in a centrifugal or tandem mirror. We find that our model outperforms other known models in approximating numerical steady-state simulations outside of a narrow range of low confining potentials. We demonstrate the model’s suitability in the high confining potential regime for applications such as loss-cone stability thresholds, fusion yields and available energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Weak decays in superheavy nuclei

Superheavy nuclei represent the extreme atoms and nuclides known at the limit of mass and charge. The observed superheavy nuclei are all proton-rich; they decay primarily by emitting 𝛼 particles and by fission with a possible small electron capture (EC) branch. Here, due to the huge atomic numbers and associated relativistic effects, EC decays of superheavy systems are expected to differ from what is known in lighter nuclei. In this letter, using the quantified relativistic nuclear density functional theory and the quasiparticle random-phase approximation with the interaction optimized to experimental EC/𝛽 ± -decay half-lives, and Gamow-Teller resonance energies, we study the EC/𝛽 ± -decays in 𝑍=101–118 nuclei. Both allowed (1 + ) and first-forbidden (0 − ,1 − and 2 − ) transitions are considered. We show that the first-forbidden 1 − transitions dominate the decay rates in almost all studied nuclei. For proton-rich nuclei, EC dominates over 𝛽 + decay. Based on calculations with two relativistic energy density functionals, we identify 45 candidate nuclei in which a competition between weak decays and 𝛼 decay and spontaneous fission is expected.

A ≥ 220↗
DOE OSTI2025

The Effect of Core-Hole Shape on Attosecond Valence Electron Dynamics

Rapid X-ray ionization of a core electron is known to induce the attosecond motion of valence electrons; however, the effect of core-hole shape on the triggered dynamics remains relatively unknown. In this work, the sub-four fs response of prototypical functionalized/heterocyclic/polycyclic molecules was simulated using real-time time-dependent density functional theory (RT-TDDFT), a sudden approximation core-hole, and phenomenologically Auger–Meitner (AM) decay. These molecules included fluorobenzene, chlorobenzene, bromobenzene, phenol, thiophenol, pyridine, phosphorus, and azulene. It is observed that the valence electron dynamics are essentially independent of the core-hole created, provided that it is ionized from an inner-shell orbital and not an inner-valence orbital. This has broad implications for free-electron laser studies of X-ray pumped attosecond processes since the flexibility in edge allows for a wide range of experimental modalities, core-holes with longer AM lifetimes, and molecular targets.

aromatic compounds↗
DOE OSTI2025

A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters

We introduce a novel method for bounding high-order multi-dimensional polynomials in finite element approximations. The method involves precomputing optimal piecewise-linear bounding boxes for polynomial basis functions, which can then be used to locally bound any combination of these basis functions. This approach can be applied to any element/basis type at any approximation order, can provide local (i.e., subcell) extremum bounds to a desired level of accuracy, and can be evaluated efficiently on-the-fly in simulations. Furthermore, we show that this approach generally yields more accurate bounds in comparison to traditional methods based on convex hull properties (e.g., Bernstein polynomials). Furthermore, the efficacy of this technique is shown in applications such as mesh validity checks and optimization for high-order curved meshes, where positivity of the element Jacobian determinant can be ensured throughout the entire element, and continuously bounds-preserving limiters for hyperbolic systems, which can enforce maximum principle bounds across the entire solution polynomial.

Bounding box↗
DOE OSTI2026

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2024

Regioselective On-Surface Synthesis of [3]Triangulene Graphene Nanoribbons

The integration of low-energy states into bottom-up engineered graphene nanoribbons (GNRs) is a robust strategy for realizing materials with tailored electronic band structure for nanoelectronics. Low-energy zero-modes (ZMs) can be introduced into nanographenes (NGs) by creating an imbalance between the two sublattices of graphene. This phenomenon is exemplified by the family of [n]triangulenes (n ϵ $\mathbb{N}$). Here, we demonstrate the synthesis of [3]triangulene-GNRs, a regioregular one-dimensional (1D) chain of [3]triangulenes linked by five-membered rings. Hybridization between ZMs on adjacent [3]triangulenes leads to the emergence of a narrow band gap, E g,exp ~ 0.7 eV, and topological end states that are experimentally verified using scanning tunneling spectroscopy. Tight-binding and first-principles density functional theory calculations within the local density approximation corroborate our experimental observations. Our synthetic design takes advantage of a selective on-surface head-to-tail coupling of monomer building blocks enabling the regioselective synthesis of [3]triangulene-GNRs. Detailed ab initio theory provides insights into the mechanism of on-surface radical polymerization, revealing the pivotal role of Au-C bond formation/breakage in driving selectivity.

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

Machine-learning guided search for phonon-mediated superconductivity in boron and carbon compounds

We present a workflow that iteratively combines ab-initio calculations with a machine-learning (ML) guided search for superconducting compounds with both dynamical stability and instability from imaginary phonon modes, the latter of which have been largely overlooked in previous studies. Electron-phonon coupling (EPC) properties and critical temperature (T c ) of 417 boron, carbon, and borocarbide compounds have been calculated with density functional perturbation theory (DFPT) and isotropic Eliashberg approximation. Our study addresses T c convergence of Brillouin zone sampling with an ansatz test, stabilizing imaginary phonon modes for significant EPC contributions, and comparing the performance of two ML models, especially when including compounds of dynamical instability. We predict a few promising superconducting compounds with formation energy just above the ground state convex hull, such as Ca 5 B 3 N 6 (35 K), TaNbC 2 (28.4 K), Nb 3 B 3 C (16.4 K), Y 2 B 3 C 2 (4.0 K), Pd 3 CaB (7.0 K), MoRuB 2 (15.6 K), RuVB 2 (15.0 K), RuSc 3 C 4 (6.6 K) among others.

Nepal, Niraj K. [Ames Laboratory (AMES), Ames, IA ↗
DOE OSTI2026

Electronic transport, thermal transport, thermal expansion, and magnetization in the strongly correlated metal LaNi⁢O 3

Perovskite structured LaNiO 3 is a strongly correlated metal with intriguing thermal and magnetic properties. The volume dependence of calculated and measured physical properties can add additional critical information to develop a more in-depth understanding of this strongly correlated phenomenon. Taking advantage of recent single crystal LaNiO 3 growth using the floating-zone method, we have measured the thermal expansion, the magnetostriction, and the pressure dependence of the magnetic susceptibility, which then allows derivation of the Grüneisen parameters γ e = $\frac{dlnN(E_F)}{d lnV}$, γ χ = $\frac{dlnχ}{d lnV}$, as well as of electric and thermal transport properties. We simulate the volume dependence of structural and magnetic properties using Density Functional Theory calculations at the Generalized Gradient Approximation level. A large discrepancy between experimental values and calculated ones suggests that strong correlations are likely to be dynamic in nature. This study also provides a side-by-side comparison of measurements in single crystal and polycrystalline samples of LaNiO 3 to elucidate intrinsic materials properties. A broad hump at high temperatures in the temperature dependence of magnetization found in the single crystal sample of LaNiO 3 has been rationalized by a model that includes the influence of electron correlations on the Landau diamagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Real-space local self-motion of protonated and deuterated water

Here, we report on the self-part of the Van Hove correlation function, the correlation function describing the dynamics of a single molecule, of water and deuterated water. The correlation function is determined by transforming inelastic scattering spectra of neutrons or x rays over a wide range of momentum transfer Q and energy transfer E to space R and time t. The short-range diffusivity is estimated from the Van Hove correlation function in the framework of the Gaussian approximation. The diffusivity has been found to be different from the long-range macroscopic diffusivity, providing information about local atomic dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗
DOE OSTI2024

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗
DOE OSTI2026

Large-scale calculations of 𝛽-decay rates and implications for 𝑟-process nucleosynthesis

Nuclear 𝛽 decay is a key element of the astrophysical rapid neutron capture process (𝑟 process). In this work, we present state-of-the-art global 𝛽-decay calculations based on the quantified relativistic nuclear energy density functional theory and the deformed proton-neutron quasiparticle random-phase approximation. Our analysis considers contributions from allowed and first-forbidden transitions. We used two point-coupling functionals with carefully calibrated time-odd terms and isoscalar pairing strength. The new calculations display consistent results for both employed functionals, especially near the neutron drip line, suggesting slower 𝛽 decays past the 𝑁=126 neutron shell closure than in commonly used 𝛽-decay models. The new rates, along with the existing rates based on the recent nonrelativistic global calculations, are found to slow down the synthesis of heavy elements in the 𝑟 process and significantly reduce the contribution of neutron-induced fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗