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

Resonant Raman in armchair graphene nanoribbons from first-principles

Resonant Raman spectra of armchair graphene nanoribbons (AGNRs) are computed using Density Functional Theory (DFT) and third-order perturbation theory. Results are benchmarked against available experimental data and compared to previously used theoretical approaches based on the Placzek approximation. Comparable agreement with experiments is found for both previously and presently used methods. In addition, a numerical analysis is carried out to provide a justification for the resonant modeling method based on the use of the frequency-dependent dielectric tensor in the Placzek approximation. Finally, this work also provides additional predictions and references for wide AGNRs that might be investigated with Raman scattering experiments in the future.

42 ENGINEERING↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Simple strategy for the simulation of axially symmetric large-area metasurfaces

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.

Martins, Augusto (ORCID:0000000295546481)↗
DOE OSTI2025

Producing High-fidelity Synthetic Population Ensembles at Scale

Used within social simulations, synthetic population ensembles enable uncertainty quantification (UQ) methods for obtaining more robust model inference and prediction. A synthetic population ensemble is a series of plausible virtual reconstructions of an area’s population at the granularity of people and residences, generated stochastically to preserve privacy of the source population survey’s respondents. In this paper, we demonstrate the production of large synthetic population ensembles for the US via Oak Ridge National Laboratory’s UrbanPop framework to support modeling of high spatial resolution energy affordability metrics from nationwide social surveys in collaboration with the fusionACS project. Our initial task involves creating ensembles for 17 US metropolitan areas, each consisting of 41 population instances (a base realization and 40 replicates). To accomplish this task at scale, we configured an integrated system comprised of a research cloud, virtual containerization, GPU-enhanced functionality, and a dual API/CLI to interact with UrbanPop’s maturing Likeness Python ecosystem. We observe a reduction in theoretical execution time while maintaining high-fidelity approximations of residential totals by metropolitan area and the demographic characteristics of neighborhoods. We discuss expansion of our approach to produce synthetic population ensembles for the entire US, particularly plans to establish automated workflows for job orchestration to increase computational efficiency, as well as provide outlook for broadening applications of the ensembles.

Gaboardi, James [ORNL] (ORCID:0000000247766826)↗
DOE OSTI2026

Jordan–Wigner Transformation for the Description of Strong Correlation in Fermionic Systems

Seniority is a useful way of organizing Hilbert space for strongly correlated systems. The exact zero-seniority wave function, doubly occupied configuration interaction (DOCI), provides accurate results (given the right orbitals) for many strongly correlated electronic systems but has a combinatorial computational cost. In many cases, pair coupled cluster doubles provide a polynomial-cost approximation that closely reproduces the energies of DOCI, but it breaks down in some cases and, as shown herein, it does not provide particularly good density matrices. In this article, we demonstrate that by using the Jordan–Wigner transformation to turn the seniority zero problem back into a Fermionic one, we can provide mean-field variational results of DOCI quality for the Hubbard model and a few small molecular dissociation examples, with polynomial cost, both for the energies and for density matrices, all while being protected from collapse. This success is rooted in the proof we provide, showing that the Hartree–Fock wave function on the Jordan–Wigner-transformed Hamiltonian transforms back to variational coupled cluster doubles in the seniority zero representation, but restricted to have determinant rather than permanent amplitude coefficients, without compromising its overall accuracy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Role of effective mass and long-range interactions in the band-gap renormalization of photoexcited semiconductors

Understanding how to control changes in the electronic structure and related dynamical renormalizations by external driving fields is the key for understanding ultrafast spectroscopy and applications in electronics. Here, we focus on the band gap's modulation by external electric fields and uncover the effect of band dispersion on the gap renormalization. We employ the Green's function formalism using the real-time Dyson expansion to account for dynamical correlations induced by photodoping. The many-body formalism captures the dynamics of systems with long-range interactions, carrier mobility, and variable electron and hole effective mass. We also demonstrate that mean-field simulations based on the Hartree-Fock Hamiltonian, which lacks dynamical correlations, yields a qualitatively incorrect picture of band-gap renormalization. We find the trend that increasing effective mass, thus decreasing mobility, leads to as much as a 6% enhancement in band-gap renormalization. Further, the renormalization is strongly dependent on the degree of photodoping. As the screening induced by free electrons and holes effectively reduces any long-range and interband interactions for highly excited systems, we show that there is a specific turnover point with a minimal band gap. Here, we further demonstrate that the optical gap renormalization follows the same trend though its magnitude is altered by the Moss-Burstein effect.

Approximation methods for many-body systems↗
DOE OSTI2025

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2024

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2025

Transition from weak to strong coupling in thermal gauge theories: Lessons from 𝒩 = 4 SYM theory

We investigate the transition between weak and strong coupling in thermal 𝒩 =4 supersymmetric Yang-Mills theory as a function of ’t Hooft coupling 𝜆 for several quantities of phenomenological interest for which next-to-leading order calculations are available in both regimes. We use Padé approximants to interpolate between the weak and strong coupling expansions and determine the location and width of the transition between the asymptotic regimes.

Gauge Theories↗
DOE OSTI2025

Leveraging Artificial Intelligence to Predict Novel Eutectic Alloys

The goal of this project was to train an artificial neural network (ANN) to predict the fractional composition and melting point of eutectic alloys using fundamental atomic properties as inputs. The fundamental properties considered include atomic number, atomic weight, atomic radius, valence electron concentration, electronegativity, and electron affinity. The project involved several phases, starting with data preparation, where phase diagram data was harvested from the ASM International database. Approximately 1300 binary eutectics were collected and cleaned to ensure relevance and accuracy. A regression model was selected for training, utilizing a rectified linear unit as the activation function. Various model configurations were evaluated for predictive accuracy, with validation techniques employed to ensure robustness. The model demonstrated predictive capabilities above random guessing and was able to achieve up to 11% accuracy under certain conditions. An ablative test identified atomic radius and valence electron concentration as critical inputs for model performance. Incorporating the melting point of atomic constituents improved accuracy significantly, although ultimately the model’s predictive capability still fell short of the 80% target. This report details the methodology, results, and implications of the research, contributing to the understanding of employing artificial intelligence to predict the phase transition behavior of eutectic alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The QPO in RE J1034+396 Originates in the Hot Corona

RE J1034+396 is one of the few active galactic nuclei (AGNs) with a significant quasiperiodic oscillation (QPO). The QPO has been observed in over 1 Ms of XMM-Newton observations spanning over a decade. We investigate the power spectral density function (PSD) of seven long (∼90 ks) XMM-Newton observations of the AGN RE J1034+396 in two energy bands. The soft (0.3–0.5 keV) band targets emission from the disk, while the hard (2–7 keV) band isolates the primary X-ray continuum emission from the corona. The QPO is significantly detected in the hard band of five of the seven observations. The best-fitting models indicate that the QPO detection in both bands is entirely attributable to the coronal emission with no additional contribution from the disk. This explains the strong coherence between the hard and soft bands at the QPO frequency. The covariance spectrum is consistent with this picture as the variability at QPO frequencies is attributed solely to fluctuations in the hot corona. The time lag as a function of energy is well described by a ∼2000 s intrinsic soft lag, resulting from the disk responding to emission from the corona, that undergoes phase wrapping at approximately the QPO frequency. By demonstrating that in this system the QPO arises in the corona, we provide new insights into the mechanisms generating QPOs.

Active galactic nuclei↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
DOE OSTI2024

Atomic Layer Deposition of Nickel Using Ni(dmamb)2 and ZnO Adhesion Layer Without Plasma

Abstract In this study, a novel deposition technique that utilizes diethylzinc (C 4 H 10 ZnO) with H 2 O to form a ZnO adhesion layer was proposed. This technique was followed by the deposition of vaporized nickel(II) 1-dimethylamino-2-methyl-2-butoxide (Ni(dmamb) 2 ) and H 2 gas to facilitate the deposit of uniform layers of nickel on the ZnO adhesion layer using atomic layer deposition. Deposition temperatures ranged from 220 to 300 °C. Thickness, composition, and crystallographic structure results were analyzed using spectroscopic ellipsometry, scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and X-ray diffraction (XRD), respectively. An average growth rate of approximately 0.0105 angstroms per cycle at 260 °C was observed via ellipsometry. Uniform deposition of ZnO with less than 1% of Ni was displayed by utilizing the elemental analysis function via SEM, thereby providing high-quality images. XPS revealed ionizations consistent with nickel and ZnO through the kinetic and binding energies of each detected electron. XRD provided supplemental information regarding the validity of ZnO by exhibiting crystalline attributes, revealing the presence of its hexagonal wurtzite structure.

Baker, Kaiya↗
DOE OSTI2024

A miniature jet pump for slurries

Here, this work reports on the design and performance of a jet pump system to pump slurries. The jet pump uses a high-velocity water jet to entrain and pump the slurry. This pumping action requires no valves and does not use moving solid parts. Hence, the system can be readily sealed, and the pumping action creates no observable mechanical damage to slurry particles. We tested carbon black beads suspended in water at concentrations of 0.0, 0.2, and 2.0 wt% at primary jet flow rates of 50 and 100 mL/min. Our measurements show an approximately linear decrease in entrainment ratio (slurry-to-jet flow rate) and absolute slurry flow rate with increasing adverse (load) pressure difference. The entrainment ratio is a strong function of the absolute jet flow rate. Pump performance was insensitive to slurry composition. Micron-resolution particle image velocimetry (micro-PIV) revealed qualitative changes in the pump’s internal flow due to changes in absolute jet flow rate, including evidence of the suppression of a vortex at higher flow rates. These findings establish key operational parameters for the jet pump and highlight its potential in a variety of slurry flow pumping applications. The pump is particularly suited for flow electrode slurries (e.g. as used for energy storage and water treatment) as it maintains particle integrity and can be used with standard tubing, tube fittings, and sealing methods.

42 ENGINEERING↗
DOE OSTI2026

Innovating Distributed Embedded Energy Prize (InDEEP): A Lessons Learned Report

The U.S. Department of Energy's Water Power Technologies Office (WPTO) launched the Innovating Distributed Embedded Energy Prize (InDEEP) in March 2023 to accelerate innovation in Distributed Embedded Energy Conversion Technologies (DEEC-Tec) for ocean wave energy. Administered by the National Laboratory of the Rockies (NLR) with technical support from Sandia National Laboratories (SNL), InDEEP focused on the development of small, distributed, and embeddable energy converters (DEECs) and their integration into scalable DEEC-Tec metamaterials for marine renewable energy applications. Spanning three phases over two years, InDEEP awarded approximately $2.3 million to teams from academia, industry, and startups. Phase I emphasized conceptual design. Phase II moved into the prototyping of individual DEECs. Phase III required integration into functional DEEC-Tec metamaterial prototypes. Across 60 submissions, teams explored a wide range of energy conversion mechanisms - including piezoelectric, variable-capacitance, ionic, and inductive methods. Note, the prize did not include the design nor demonstration of ocean wave energy conversion systems. Rather, the prize only required participants to design and demonstrate individual DEECs and corresponding DEEC-Tec metamaterials. This prize utilized a mix of novel and proven techniques to attract participants from outside marine energy, including an engagement leaderboard, robust recruitment, technical expert mentorship, and a suite of technical trainings. Key insights from the competition emphasized that DEEC-Tec metamaterials must be intentionally designed to produce beneficial emergent behaviors - advantages that go beyond simply combining multiple DEEC units. Top-performing teams showed that thoughtful design of system architecture, coordinated deformation, and systems adaptabilities could unlock meaningful performance gains both at the DEEC system level and DEEC-Tec metamaterial system level. A critical realization was that many DEEC-Tec metamaterials could benefit from being designed to accept lower-frequency energy inputs and shift those into higher-frequencies per each DEEC making up the respective DEEC-Tec metamaterial. Other important takeaways included the need for rigorous and quantitative performance testing, effective integration of power conditioning electronics, and the pivotal role of material science in enabling innovative, adaptive DEEC-Tec-based energy conversion designs. InDEEP also helped establish a growing DEEC-Tec community of practitioners, attracting participants from beyond traditional marine energy sectors. Through a strong support infrastructure, InDEEP fostered early-stage innovation and laid a foundation for future DEEC-Tec-based ocean wave energy conversion solutions - positioning DEEC-Tec as a promising pathway toward scalable, resilient ocean wave energy conversion. Through focused R&D of individual DEECs and their integration into DEEC-Tec metamaterials, alongside a growing, multidisciplinary community catalyzed by InDEEP, there is a strong opportunity to drive a disruptive shift in ocean wave energy conversion design and development. This convergence of novel architectures, emergent behaviors, and collaborative innovation positions DEEC-Tec as a transformative approach, moving the field from rigid, centralized energy conversion-based designs to resilient, modular systems highly adaptable for real-world ocean wave energy conversion applications.

16 TIDAL AND WAVE POWER↗