Search NASA⌕ Search

SEARCH · Search NASA

Results for “Real Functions”

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

Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core

In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Improved Methods for Modeling Functional Transition Metal Compounds in Complex Environments: Ground States, Excited States, and Spectroscopies

Rapid advances in energy applications require new theory and computational models to provide guidance for interpretation of experimental results and mechanistic understanding. New theory development is necessary to treat systems of increasing complexity, size, and relevance to real applications. Functional transition metal compounds, including molecules, clusters, nanoparticles, surfaces, and solids, provide particular promise for magnetic, optical, and catalytic applications. However, such systems can be exceptionally challenging to model. To make progress in understanding and designing transition metal compounds for energy applications, theory must be able to simulate such systems in complex environments, as well as simulate the spectra of such complex systems to provide direct connections with experiment. Leveraging the independent expertise of the team’s members, this project aimed to make inroads to the theoretical and computational challenges associated with studying transition metal compounds, their reaction chemistry, photophysics and photochemistry, and response to spectroscopic interrogation.

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

Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations

In this article, we report a scheme to analyze and visualize the energy density fluctuations during the real-time time-dependent density functional theory (RT-TDDFT) simulations. Using Ag 4 –N 2 complexes as examples, here we show that the grid-based Kohn–Sham energy density can be computed at each time step using a procedure from Nakai and coworkers. Then the instantaneous energy of each molecular fragment (such as Ag 4 and N 2 ) can be obtained by partitioning the Kohn–Sham energy densities using Becke or fragment-based Hirshfeld (FBH) scheme. A strong orientation-dependence is observed for the energy flow between the Ag 4 cluster and a nearby N 2 molecule in the RT-TDDFT simulations. Future applications of such an energy density analysis in electron dynamics simulations are discussed.

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

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method

In this paper, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Electron dynamics in extended systems within real-time time-dependent density-functional theory

Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density-functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. Here we extensively explore the shortcomings for describing electron–electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward.

2D materials↗
DOE OSTI2022

Instantaneous tracking of earthquake growth with elastogravity signals

Rapid and reliable estimation of large earthquake magnitude (above 8) is key to mitigating the risks associated with strong shaking and tsunamis. Standard early warning systems based on seismic waves fail to rapidly estimate the size of such large earthquakes. Geodesy-based approaches provide better estimations, but are also subject to large uncertainties and latency associated with the slowness of seismic waves. Recently discovered speed-of-light prompt elastogravity signals (PEGS) have raised hopes that these limitations may be overcome, but have not been tested for operational early warning. Here we show that PEGS can be used in real time to track earthquake growth instantaneously after the event reaches a certain magnitude. We develop a deep learning model that leverages the information carried by PEGS recorded by regional broadband seismometers in Japan before the arrival of seismic waves. After training on a database of synthetic waveforms augmented with empirical noise, we show that the algorithm can instantaneously track an earthquake source time function on real data. Our model unlocks ‘true real-time’ access to the rupture evolution of large earthquakes using a portion of seismograms that is routinely treated as noise, and can be immediately transformative for tsunami early warning.

58 GEOSCIENCES↗
DOE OSTI2022

Quasiparticle excitations in a one-dimensional interacting topological insulator: Application for dopant-based quantum simulation

Here, we study the effects of electron-electron interactions on the charge excitation spectrum of the spinful Su-Schrieffer-Heeger (SSH) model, a prototype of a one-dimensional bulk obstructed topological insulator. In view of recent progress in the fabrication of dopant-based quantum simulators we focus on experimentally detectable signatures of interacting topology in finite lattices. To this end we use Lanczos-based exact diagonalization to calculate the single-particle spectral function in real space which generalizes the local density of states to interacting systems. Its spatial and spectral resolution allows for the direct investigation and identification of edge states. By studying the noninteracting limit, we demonstrate that the topological in-gap states on the boundary are robust against both finite-size effects as well as random bond and onsite disorder which suggests the feasibility of simulating the SSH model in engineered dopant arrays in silicon. While edge excitations become zero-energy spinlike for any finite interaction strength, our analysis of the spectral function shows that the single-particle charge excitations are gapped out on the boundary. Despite the loss of topological protection, we find that these edge excitations are quasiparticlelike as long as they remain within the bulk gap. Above a critical interaction strength of U c ≈5t these quasiparticles on the boundary lose their coherence which is explained by the merging of edge and bulk states. This is in contrast to the many-body edge excitations which survive the limit of strong coupling, as established in the literature. Our findings show that for moderate repulsive interactions the nontrivial phase of the interacting SSH model can be detected through remnant signatures of topological single-particle states using single-particle local measurement techniques such as scanning tunneling spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

The Next Stage for Smart Cities: Metrics to Spur an Evolution from Technology for Efficiency to Technology for People's Needs

The Smart Cities movement is about a decade old, energized by the USDOT Smart Cities grant that sparked cities to dream big. Now, 10 years later, and multitudes of conferences and papers later, the biggest question regarding Smart Cities is whether any progress has been made? Smart Cities forums are typically led by the technology sector, either telecommunications, internet access, or most recently artificial intelligence (AI) concerns. Each of these areas have grown and found new applications in the context of Smart Cities, however, these technology applications are just tools toward ends. Do digital communications and assistance enable or foster improved access to real, tactile, non-virtual goods and services such as transportation, food, housing, trash removal, health care, education, social interaction? These and other topics are not 'big tech', but are each critical to quality of life outcomes and are day-to-day necessities, comprising many of the touch points through which citizens interact with municipal agencies. Regarding Smart Cities, relationships have not been fully refined as to how best to merge the emerging generation of technology and related trends with more mundane but essential requirements of tactile services. How should municipal managers and operators bring the range of new technology tools (e.g. on-demand transit, automated vehicles, digital connectivity) to bear to improve functions serving real- physical- non-virtual needs? How do we make cities more efficient in providing access to services and opportunities for improved quality of life? This paper examines the current Smart City movement, emerging trends, and the need to put metrics with respect to the delivery, access, and provision of real-life needs. Using a few case studies of past, present, and projected future scenarios, this paper attempts to bring Smart Cities into the real cities domain, perhaps pushing the moniker away from 'Smart' to 'Efficient' - something that can be measured, counted, and assessed. Through the consideration of the next stage of the evolution of cities from 'Smart' to 'Efficient', the combined efforts and experience of multiple researchers at the NREL and their network of collaborators highlights the need for updates metrics that can be used to start guiding efforts toward productive initiatives to improve quality of life.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

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

Stochastic evaluation of fourth-order many-body perturbation energies

A scalable, stochastic algorithm evaluating the fourth-order many-body perturbation (MP4) correction to energy is proposed. Three hundred Goldstone diagrams representing the MP4 correction are computer generated and then converted into algebraic formulas expressed in terms of Green’s functions in real space and imaginary time. They are evaluated by the direct (i.e., non-Markov, non-Metropolis) Monte Carlo (MC) integration accelerated by the redundant-walker and control-variate algorithms. The resulting MC-MP4 method is efficiently parallelized and is shown to display O(n 5.3 ) size-dependence of cost, which is nearly two ranks lower than the O(n 7 ) dependence of the deterministic MP4 algorithm. Furthermore, it evaluates the MP4/aug-cc-pVDZ energy for benzene, naphthalene, phenanthrene, and corannulene with the statistical uncertainty of 10 mE h (1.1% of the total basis-set correlation energy), 38 mE h (2.6%), 110 mE h (5.5%), and 280 mE h (9.0%), respectively, after about 10 9 MC steps.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Precision redshift-space galaxy power spectra using Zel'dovich control variates

Numerical simulations in cosmology require trade-offs between volume, resolution and run-time that limit the volume of the Universe that can be simulated, leading to sample variance in predictions of ensemble-average quantities such as the power spectrum or correlation function(s). Sample variance is particularly acute at large scales, which is also where analytic techniques can be highly reliable. This provides an opportunity to combine analytic and numerical techniques in a principled way to improve the dynamic range and reliability of predictions for clustering statistics. In this paper we extend the technique of Zel'dovich control variates, previously demonstrated for 2-point functions in real space, to reduce the sample variance in measurements of 2-point statistics of biased tracers in redshift space. We demonstrate that with this technique, we can reduce the sample variance of these statistics down to their shot-noise limit out to k ~ 0.2 h Mpc -1 . This allows a better matching with perturbative models and improved predictions for the clustering of e.g. quasars, galaxies and neutral Hydrogen measured in spectroscopic redshift surveys at very modest computational expense. We discuss the implementation of ZCV, give some examples and provide forecasts for the efficacy of the method under various conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dark energy survey year 3 results: cosmological constraints from the analysis of cosmic shear in harmonic space

ABSTRACT We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-Cℓ method and complement the analysis of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, making this analysis an important cross-check. Using the same fiducial Lambda cold dark matter model as in the DES Y3 real-space analysis, we find ${S_8 \equiv \sigma _8 \sqrt{\Omega _{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to S8 = 0.784 ± 0.026 when including shear ratios. This result is within expected statistical fluctuations from the real-space constraint, and in agreement with DES Y3 analyses of non-Gaussian statistics, but favours a slightly higher value of S8, which reduces the tension with the Planck 2018 constraints from 2.3σ in the real space analysis to 1.5σ here. We explore less conservative intrinsic alignments models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}\, {h}\, {\rm Mpc}^{-1}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, found to be about 20 per cent lower than predicted by Planck 2018, as reflected by the lower S8 value.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Microscopic scattering approach to in-gap states: Cr adatoms on superconducting 𝛽−Bi 2 ⁢Pd

We develop a microscopic scattering formalism to describe Yu-Shiba-Rusinov (YSR) states due to a single Cr adatom on the Bi-terminated surface of 𝛽−Bi 2 ⁢Pd, by combining ab initio Wannier functions with a real-space Green's function approach in the Bogoliubov–de Gennes formalism. Our framework reproduces key scanning tunneling spectroscopy features, including a single particle-hole asymmetric YSR peak and isotropic 𝑑⁢𝐼/𝑑⁢𝑉 maps around the impurity. Decomposing the YSR states reveals contributions from four nearly degenerate 𝐶 4⁢𝑣 representations, with energy broadening masking their individual signatures. Spin-orbit coupling induces partial spin polarization, while the spatial asymmetry between particle and hole components arises from Cr 𝑑−Bi 𝑝 hybridization. These results highlight the importance of realistic band structures and microscopic modeling for interpreting STM data and provide a foundation for studying impurity chains hosting topological excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Automated label flows for excited states of correlation functions in lattice gauge theory

Extracting excited states from lattice gauge theory correlation functions can be achieved through $χ^2$ minimization fits or algebraic approaches such as the variational method and Prony's method. Performing any kind of error analysis often leads to overlapping confidence regions of model parameters, even when the spectrum is not particularly dense. To correctly estimate errors, one must beware of mislabeling the states. In this work, we provide an algorithm that we call automated label flows which consistently and systematically identifies a deterministic labeling of states. This is a black-box approach in the sense that it gives a sensible set of labels without user guidance. Further, as an example, we pair one black box method with another, analyzing a lattice correlation function from real data using automated label flows in the context of Prony's method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Investigation of Defect Physics for Efficient, Durable and Ubiquitous Perovskite Solar Modules (Final Technical Report)

Organometal halide perovskite solar cells have experienced eye-catching improvements in its recent few years. It serves as one of the most promising candidates to replace the currently widely used silicon-based solar modules. To implement its final step to the real application, functional longevity becomes the dernier continent to conquer. As a polycrystalline material, defect plays critical role in the efficiency and stability of the perovskite solar cells. Thus, the investigation of the defect physics of the perovskite layer is indispensable in this research field. However, hard evidence and a consensus are still lacking in terms of the specific nature of the defects and their effects on performance and hysteresis, and perhaps even more importantly, there is absence of fundamental understanding of the correlations between the defects and long-term operational stability of the device. A more fundamental understanding of the nature of defects in perovskite materials is of paramount importance to progress their efficiency and durability. In this work we propose in-depth studies of correlations of defects with performance and stability of perovskite solar cells. Our project aimed 1) to investigate the defects physics in perovskite solar cells, and 2) to develop a comprehensive understanding and physical model of defects and its influence on performance and stability of perovskite solar cells. With the support from program manager, Peter Lobaccaro, and the Solar Energy Technologies Office of U.S. Department of Energy, the project ends with impact achievements. Our research results have systematically provided strategies to analyze the influences of constructive molecular configurations to the charged defects in the perovskite lattices and developed in-depth understanding of chemical additive approach to improve the perovskite solar cell performance and stability. As history has shown us, control over defect properties of semiconductor materials is the key to achieving high performance and low cost devices. Therefore, the potential impact of unlocking the understanding and manipulation of defects in perovskites is great, enabling this technology to realize SETO goals. The research project is highly productive with 18 published papers in three years in top-level journals such as Science, Nature, Nature Materials, Nature Communications, Journal of American Chemistry Society, Joule, Advanced Materials, and Nano Letters. These works have drawn great attention nationwide with notable total citations over 700 times from 2020 to 2022.

14 SOLAR ENERGY↗
DOE OSTI2025

Fermilab PIP II Machine Protection System Digitized Data Noise Elimination Scheme and Its FPGA Implementation

In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can con-taminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is re-quired. Several real-time functional blocks were simu-lated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid base-line by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.

Wu, Jinyuan [Fermilab] (ORCID:0000000344329521)↗