Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation methods”

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 OSTI2020

Capturing multireference excited states by constrained-density-functional theory

The computation of excited electronic states with commonly employed (approximate) methods is challenging, typically yielding states of lower quality than the corresponding ground state for a higher computational cost. In this work, we present a mean-field method that extends the previously proposed excited constrained-density-functional theory (XCDFT) from single Slater determinants to ensemble one-body reduced density matrices for computing low-lying excited states. The method still retains an associated computational complexity comparable to a semilocal density functional theory (DFT) calculation while at the same time it is capable of approaching states with multireference character. Further, we benchmark the quality of this method on well-established test sets, finding good descriptions of the electronic structure of multireference states and maintaining an overall accuracy for the predicted excitation energies comparable to semilocal time-dependent DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Here in this work, we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CEs for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions, and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for nonzero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. II. The self-consistent cumulant expansion

Here in this paper we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattice sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demonstration that the SC-CE corrects many deficits of the standard perturbative CE in the (nondispersive) Holstein model. The natural phonon damping that arises in the more realistic dispersed Holstein model renders the SC-CE stable, allowing for a complete assessment of the method. Here, we find that self-consistency dramatically corrects many of the failures found in the perturbative CE, but also introduces some unphysical features. Finally, we comment on the potential use of SC-CE as a tool for calculating Green's functions in generic many-body problems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonequilibrium DMFT + CPA for correlated disordered systems

We present a solution for the nonequilibrium dynamics of an interacting disordered system. The approach adapts the combination of the equilibrium dynamical mean-field theory and the equilibrium coherent potential approximation methods to the nonequilibrium many-body formalism, using the Kadanoff-Baym-Keldysh complex time contour, for the dynamics of interacting disordered systems away from equilibrium. We use our time domain solution to obtain the equilibrium density of states of the disordered interacting system described by the Anderson-Hubbard model, bypassing the necessity for the cumbersome analytical continuation process. We further apply the nonequilibrium solution to the interaction quench problem for an isolated disordered system. Here, the interaction is abruptly changed from zero (noninteracting system) to another constant (finite) value at which it is subsequently kept. We observe via the time dependence of the potential, kinetic, and total energies the effect of disorder on the relaxation of the system as a function of final interaction strength. Finally, the real-time approach has the potential to shed light on the fundamental role of disorder in the nonequilibrium dynamics of interacting quantum systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Experimental study of the low-lying negative-parity states in 11 Be using the 12 B( d , 3 He) 11 Be reaction

In this work, low-lying negative-parity states in 11 Be having dominant p-wave neutron configurations were studied using the 12 B(d, 3 He) 11 Be proton-removal reaction in inverse kinematics. The 1/2 $\overset{-}{1}$ state at 0.32 MeV, the 3/2 $\overset{-}{1}$ state at 2.56 MeV and one or both of the states including the 5/2 $\overset{-}{1}$ level at 3.89 MeV and the 3/2 $\overset{-}{2}$ level at 3.96 MeV were populated in the present reaction. Spectroscopic factors were determined from the differential cross sections using a distorted wave Born approximation method. The p-wave proton removal strengths were well described by the shell model calculations while the Nilsson model calculation underestimates the spectroscopic factors for the higher excited states. Results from both variational Monte Carlo and no-core shell model calculations were also compared with the experimental observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Numerical and perturbative computations of the fuzzy dark matter model

We investigate nonlinear structure formation in the fuzzy dark matter (FDM) model using both numerical and perturbative techniques. On the numerical side, we examine the virtues and limitations of a Schrödinger-Poisson solver (wave formulation) versus a fluid dynamics solver (Madelung formulation). On the perturbative side, we carry out a computation of the one-loop mass power spectrum, i.e. up to third order in perturbation theory. We find that (1) in many situations, the fluid dynamics solver is capable of producing the expected interference patterns, but it fails in situations where destructive interference causes the density to vanish—a generic occurrence in the nonlinear regime. (2) The Schrödinger-Poisson solver works well in all test cases, but it is demanding in resolution: suppose one is interested in the mass power spectrum on large scales, it is not sufficient to resolve structure on those same scales; one must resolve the relevant de Broglie scale which is often smaller. The fluid formulation does not suffer from this issue. (3) We compare the one-loop mass power spectrum from perturbation theory against the mass power spectrum from the Schrödinger-Poisson solver, and find good agreement in the mildly nonlinear regime. We contrast fluid perturbation theory with wave perturbation theory; the latter has a more limited range of validity. (4) As an application, we compare the Lyman-alpha forest flux power spectrum obtained from the Schrödinger-Poisson solver versus one from an N-body simulation (the latter is often used as an approximate method to make predictions for FDM). At redshift 5, the two, starting from the same initial condition, agree to better than 10% on observationally relevant scales as long as the FDM mass exceeds 2 × 10 –23 eV. Furthermore, we emphasize that the so called quantum pressure is capable of both enhancing and suppressing fluctuations in the nonlinear regime—which dominates depends on the scale and quantity of interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Variational Principle for Mass Transport

A variation principle for mass transport in solids is derived that recasts transport coefficients as minima of local thermodynamic average quantities. The result is independent of diffusion mechanisms and applies to amorphous and crystalline systems. This unifies different computational approaches for diffusion and provides a framework for the creation of new approximation methods with error estimation. It gives a different physical interpretation of the Green function. Lastly, the variational principle quantifies the accuracy of competing approaches for a nontrivial diffusion problem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Continuous Spin Excitations in the Three-Dimensional Frustrated Magnet K 2 Ni 2 ⁢(SO 4 ) 3

Continuous spin excitations are widely recognized as one of the hallmarks of novel spin states in quantum magnets, such as quantum spin liquids (QSLs). Here, we report the observation of such kind of excitations in K 2 Ni 2 (SO 4 ) 3 , which consists of two sets of intersected spin-1 (Ni 2+ ) trillium lattices. Our inelastic neutron scattering measurement on single crystals clearly shows a dominant excitation continuum, which exhibits a distinct temperature-dependent behavior from that of spin waves and is rooted in strong quantum spin fluctuations. Further using the self-consistent-Gaussian-approximation method, we determine that the fourth- and fifth nearest-neighbor exchange interactions are dominant. These two bonds together form a unique three-dimensional network of corner-sharing tetrahedra, which we name as a “hypertrillium” lattice. In conclusion, our results provide direct evidence for the existence of QSL features in K 2 Ni 2 (SO 4 ) 3 and highlight the potential for the hypertrillium lattice to host frustrated quantum magnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electronic Coherences in Molecules: The Projected Nuclear Quantum Momentum as a Hidden Agent

Electronic coherences are key to understanding and controlling photoinduced molecular transformations. Here, we identify a crucial quantum-mechanical feature of electron-nuclear correlation, the projected nuclear quantum momenta, essential to capture the correct coherence behavior. For simulations, we show that, unlike traditional trajectory-based schemes, exact-factorization-based methods approximate these correlation terms and correctly capture electronic coherences in a range of situations, including their spatial dependence, an important aspect that influences subsequent electron dynamics and that is becoming accessible in more experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Negative thermal expansion in the Ruddlesden-Popper calcium titanates

Materials exhibiting negative thermal expansion (NTE) are important for the fabrication and operation of microelectronic devices and optical systems. In this paper, combining density-functional-theory calculations and the self-consistent quasiharmonic approximation method, we investigate the variation of the NTE mechanisms active in layered calcium titanate oxides. We find a NTE trend with layer number for the orthorhombic materials comprising the RP series, but the monoclinic polymorph is an outlier.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Coupling Visualization, Simulation, and Deep Learning for Ensemble Steering of Complex Energy Models

We describe a new framework that allows users to explore and steer ensembles of energy systems simulations by coupling multiple energy models and interactive visualization through a dataflow API. Through the visual interface, users can interactively explore complex parameter spaces populated by hundreds, or thousands, of simulation runs and interactively spawn new simulations to 'fill in' regions of interest in the parameter space. The computational and visualization capabilities reside within a general-purpose dataflow architecture for connecting producers of multidimensional timeseries data, such as energy simulations, with consumers of that data, whether they be visualizations, statistical analyses, or datastores. Fast computation and agile dataflow can enhance the engagement with energy simulations, allowing users to populate the parameter space in real time. However, many energy simulations are far too slow to provide an interactive response. To support interactive feedback, we are creating reduced-form simulations developed through machine learning techniques, which provide statistically sound estimates of the results of the full simulations at a fraction of the computational cost. These reduced-form simulations have response times on the order of seconds, suitable for real-time human-in-the-loop design and analysis. The approximation methods apply to a wide range of computational models, including supply-chain models, electric power grid simulations, and building models. Such reduced-form representations do not replace or re-implement existing simulations, but instead supplement them by enabling rapid scenario design and exploration for large ensembles of simulations. The improved understanding, facilitated by the reduced-form models, dataflow API, and visualization tools, allows researchers to better allocate computational resources to capture informative relationships within the system as well as provide a low-cost method for validating and quality-checking large-scale modeling efforts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Data-Driven Performance Optimization of Gamma Spectrometers With Many Channels

In gamma spectrometers with variable spectroscopic performance across many channels (e.g., many pixels or voxels), a tradeoff exists between including data from successively worse-performing readout channels and increasing efficiency. Brute-force calculation of the optimal set of included channels is exponentially infeasible as the number of channels grows, and approximate methods are required. In this work, we present a data-driven framework for attempting to find near-optimal sets of included detector channels. The framework leverages non-negative matrix factorization (NMF) to learn the behavior of gamma spectra across the detector and clusters similarly-performing detector channels together. Performance comparisons are then made between spectra with channel clusters removed, which is more feasible than brute force. The framework is general and can be applied to arbitrary, user-defined performance metrics depending on the application. We apply this framework to optimizing gamma spectra measured by H3D M400 CdZnTe (CZT) spectrometers, which exhibit variable performance across their crystal volumes. In particular, we show several examples optimizing various performance metrics for uranium and plutonium gamma spectra in non-destructive assay (NDA) for nuclear safeguards, and explore trends in performance versus parameters such as clustering algorithm type. We also compare the NMF + clustering pipeline to several non-machine-learning (ML) algorithms, including several greedy algorithms. Although, we find that the NMF + clustering pipeline tends to find the best-performing set of detector voxels, significantly improving over the unoptimized spectra, but that a greedy accumulation of spectra segmented by detector depth can, in some cases, give similar performance improvements in much less computation time.

Energy resolution↗
DOE OSTI2021

Modeling State Transition and Head-Dependent Efficiency Curve for Pumped Storage Hydro in Look-Ahead Dispatch

As one of the most widely installed utility-scale storage facilities, pumped storage hydro (PSH) plays an essential role in providing flexibility for power systems worldwide. Thus, to accurately quantify the flexibility of PSH units in operational optimization problems is important. The conventional PSH models in the literature rarely consider detailed state transitions due to their hourly-based settings. However, it becomes imperative for operational optimization on short-term intervals, especially with increasing shares of renewable energy in power systems. To this end, this paper presents a novel deterministic PSH model that considers the transition time and trajectory between three states of PSH units in look-ahead dispatch. Moreover, to better characterize the varying efficiency of PSH units with water head and flow rate, this paper proposes to model detailed headdependent efficiency curves (hereafter called input-output curves) in look-ahead dispatch, which encounters heavy computational burdens when short time intervals are applied. In this work, a zig-zag piece-wise linear approximation method is used for input-output curve modeling. This can enable an accurate quantification for variable efficiency and head dependence in a computationally effective manner. Here, numerical results are presented to show performances of the proposed PSH model in both flexibility quantification and computation time.

13 HYDRO ENERGY↗
DOE OSTI2019

Machine Learning-Based Anomaly Detection for Load Forecasting Under Cyberattacks

Accurate load forecasting can create both economic and reliability benefits for power system operators. However, the cyberattack on load forecasting may mislead operators to make unsuitable operational decisions for the electricity delivery. To effectively and accurately detect these cyberattacks, this paper develops a machine learning-based anomaly detection (MLAD) methodology. First, load forecasts provided by neural networks are used to reconstruct the benchmark and scaling data by using the k-means clustering. Second, the cyberattack template is estimated by the naive Bayes classification based on the cumulative distribution function and statistical features of the scaling data. Finally, the dynamic programming is utilized to calculate both the occurrence and parameter of one cyberattack on load forecasting data. Here, a widely used symbolic aggregation approximation method is compared with the developed MLAD method. Numerical simulations on the publicly load data show that the MLAD method can effectively detect cyberattacks for load forecasting data with relatively high accuracy. Also, the robustness of MLAD is verified by thousands of attack scenarios based on Monte Carlo simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Anisotropic Thermal Response of Packed Copper Wire

The apparent thermal conductivity of packed copper wire test specimens was measured parallel and perpendicular to the axis of the wire using laser flash, transient plane source, and transmittance test methods. Approximately 50% wire packing efficiency was produced in the specimens using either 670- or 925-μm-diameter copper wires that both had an insulation coating thickness of 37 μm. The interstices were filled with a conventional varnish material and also contained some remnant porosity. The apparent thermal conductivity perpendicular to the wire axis was about 0.5–1 W/mK, whereas it was over 200 W/mK in the parallel direction. The Kanzaki model and an finite element analysis (FEA) model were found to reasonably predict the apparent thermal conductivity perpendicular to the wires but thermal conductivity percolation from nonideal wire-packing may result in their underestimation of it.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Beyond-classical computation in quantum simulation

Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. However, establishing this capability, especially for impactful and meaningful problems, remains a central challenge. Here, we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schrödinger equation. We demonstrate area-law scaling of entanglement in the model quench dynamics of two-, three-, and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for matrix-product-state approaches. We show that several leading approximate methods based on tensor networks and neural networks cannot achieve the same accuracy as the quantum annealer within a reasonable time frame. Thus, quantum annealers can answer questions of practical importance that may remain out of reach for classical computation.

King, Andrew D. [D-Wave Quantum Inc., Burnaby, BC ↗
DOE OSTI2015

Exactly solvable model of resonance tunneling of an electromagnetic wave in plasma containing short-scale inhomogeneities

Different variants of resonance tunneling of a transverse electromagnetic wave through a plasma layer containing short-scale (subwavelength) inhomogeneities, including evanescence regions to which approximate methods are inapplicable, are analyzed in the framework of an exactly solvable one-dimensional model. Complex plasma density profiles described by a number of free parameters determining the permittivity modulation depth, the characteristic scale lengths of plasma structures, their number, and the thickness of the inhomogeneous plasma layer are considered. It is demonstrated that reflection-free propagation of the wave incident on the layer from vacuum (the effect of wave-barrier transillumination) can be achieved for various sets of such structures, including plasma density profiles containing a stochastic component. Taking into account cubic nonlinearity, it is also possible to obtain an exact solution to the one-dimensional problem on the nonlinear transillumination of nonuniform plasma. In this case, the thicknesses of the evanescence regions decrease appreciably. The problem of resonance tunneling of electromagnetic waves through such barriers is of interest for a number of practical applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗