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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Theory of optical transitions in π-conjugated macrocycles

We describe a theoretical and computational investigation of the optical properties of π-conjugated macrocycles. Since the low-energy excitations of these systems are Frenkel excitons that couple to high-frequency dispersionless phonons, we employ the quantized Frenkel-Holstein model and solve it via the density matrix renormalization group (DMRG) method. First we consider optical emission from perfectly circular systems. Owing to optical selection rules, such systems radiate via two mechanisms: (i) within the Condon approximation, by thermally induced emission from the optically allowed j = ± 1 states and (ii) beyond the Condon approximation, by emission from the j = 0 state via coupling with a totally non-symmetric phonon (namely, the Herzberg-Teller effect). Using perturbation theory, we derive an expression for the Herzberg-Teller correction and show via DMRG calculations that this expression soon fails as ħ ω/J and the size of the macrocycle increase. Next, we consider the role of broken symmetry caused by torsional disorder. In this case the quantum number j no longer labels eigenstates of angular momentum, but instead labels localized local exciton groundstates (LEGSs) or quasi-extended states (QEESs). As for linear polymers, LEGSs define chromophores, with the higher energy QEESs being extended over numerous LEGSs. Within the Condon approximation (i.e., neglecting the Herzberg-Teller correction) we show that increased disorder increases the emissive optical intensity, because all the LEGSs are optically active. We next consider the combined role of broken symmetry and curvature, by explicitly evaluating the Herzberg-Teller correction in disordered systems via the DMRG method. The Herzberg-Teller correction is most evident in the emission intensity ratio, I{sub 00}/I{sub 01}. In the Condon approximation I{sub 00}/I{sub 01} is a constant function of curvature, whereas in practice it vanishes for closed rings and only approaches a constant in the limit of vanishing curvature. We calculate the optical spectra of a model system, cyclo-poly(para-phenylene ethynylene), for different amounts of torsional disorder within and beyond the Condon approximation. We show how broken symmetry and the Herzberg-Teller effect explain the spectral features. The Herzberg-Teller correction to the 0–1 emission vibronic peak is always significant. Finally, we note the qualitative similarities between the optical properties of conformationally disordered linear polymers and macrocycles in the limit of sufficiently large disorder, because in both cases they are determined by the optical properties of curved chromophores.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2018

First-Principle Calculations of Structural, Elastic, and Electronic Properties of Intermetallic Rare Earth R{sub 2}Ni{sub 2}Pb (R = Ho, Lu, and Sm) Compounds

The structural, elastic, and electronic properties of rare earth intermetallic R{sub 2}Ni{sub 2}Pb (where R = Ho, Lu, and Sm) compounds were investigated with the density functional theory (DFT) calculations. The calculations are performed using the full potential-linearized augmented plane wave (FP-LAPW) method within the framework of local density approximation (LDA). The calculated values of the equilibrium lattice constants were in agreement with the available experimental values. The elastic constants (C{sub ij}) were also calculated to understand the mechanical properties and structural stability of the compounds. Furthermore, the density of states and the charge density distributions of the compounds were calculated to understand the nature of the bonding in the material. The calculated results are in accordance with the available data in the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

From Molecules to Solids: A vdW-DF-C09 Case Study of the Mercury Dihalides

The mercury dihalides show a remarkable diversity in the structural preferences in their minimum energy structure types, spanning molecular to strongly bound ionic solids. A challenge in the development of density functional methods for extended systems is to arrive at strategies that serve equally well such a broad range of bonding modes or structural preferences. The chemical bonding and the stabilities of mercury dihalides and the general utility and reliability of the van der Waals density functional with C09 exchange (vdW-DF-C09) in predicting or describing the energetics and structural preferences in these metal dihalides is examined. Here, we show that, in contrast with the uncorrected generalized gradient approximation of the Perdew-Burke-Erzenhoff (PBE) exchange-correlation functional, qualitative and quantitative patterns in the bonding of the mercury dihalide solids are well reproduced with vdW-DF-C09 for the full series of HgX 2 systems for X = F, Cl, Br, and I. The possible existence of a low-temperature cotunnite polymorph for HgF 2 and PbF 2 is posited.

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

Atomic density functional and diagram of structures in the phase field crystal model

The phase field crystal model provides a continual description of the atomic density over the diffusion time of reactions. We consider a homogeneous structure (liquid) and a perfect periodic crystal, which are constructed from the one-mode approximation of the phase field crystal model. A diagram of 2D structures is constructed from the analytic solutions of the model using atomic density functionals. The diagram predicts equilibrium atomic configurations for transitions from the metastable state and includes the domains of existence of homogeneous, triangular, and striped structures corresponding to a liquid, a body-centered cubic crystal, and a longitudinal cross section of cylindrical tubes. The method developed here is employed for constructing the diagram for the homogeneous liquid phase and the body-centered iron lattice. The expression for the free energy is derived analytically from density functional theory. The specific features of approximating the phase field crystal model are compared with the approximations and conclusions of the weak crystallization and 2D melting theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Perspective: Fundamental aspects of time-dependent density functional theory

In the thirty-two years since the birth of the foundational theorems, time-dependent density functional theory has had a tremendous impact on calculations of electronic spectra and dynamics in chemistry, biology, solid-state physics, and materials science. Alongside the wide-ranging applications, there has been much progress in understanding fundamental aspects of the functionals and the theory itself. This Perspective looks back to some of these developments, reports on some recent progress and current challenges for functionals, and speculates on future directions to improve the accuracy of approximations used in this relatively young theory.

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

Estimation of Sea Spray Aerosol Surface Area Over the Southern Ocean Using Scattering Measurements

This study focuses on methods to estimate dry marine aerosol surface area (SA) from bulk optical measurements. Aerosol SA is used in many models' ice nucleating particle (INP) parameterizations, as well as influencing particle light scattering, hygroscopic growth, and reactivity, but direct observations are scarce in the Southern Ocean (SO). Two campaigns jointly conducted in austral summer 2018 provided co-located measurements of aerosol SA from particle size distributions and lidar to evaluate SA estimation methods in this region. Mie theory calculations based on measured size distributions were used to test a proposed approximation for dry aerosol SA, which relies on estimating effective scattering efficiency (Q) as a function of Ångström exponent (å). For distributions with dry å < 1, Q = 2 was found to be a good approximation within ±50%, but for distributions with dry å > 1, an assumption of Q = 3 as in some prior studies underestimates dry aerosol SA by a factor of 2 or more. We propose a new relationship between dry å and Q, which can be used for –0.2 < å < 2, and suggest å = 0.8 as the cutoff between primary and secondary marine aerosol-dominated distributions. Application of a published methodology to retrieve dry marine aerosol SA from lidar extinction profiles overestimated aerosol SA by a factor of 3–5 during these campaigns. Using Microtops aerosol optical thickness measurements, we derive alternative lidar conversion parameters from our observations, applicable to marine aerosol over the SO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fast Emulation of Expensive Simulations using Approximate Gaussian Processes [Slides]

Nuclear Computational Low-Energy Initiative (NUCLEI) collaboration uses Density Functional Theory (DFT) simulations to predict the structure and binding energies of nuclei over a wide range of proton (Z) and neutron (N) numbers. The DFT simulations utilize a particular parameterization of a Skyrme energy density functional called UNEDF1 which depends on 12 free parameters that must be fit to data (M Kortelainen et al 2014). Fitting involves comparing (e.g.) predicted binding energies of nuclei to experimentally measured values. We use only binding energies as observables, but DFT with UNEDF1 will predict structure (shape) observables as well. In this work, assessing the capability of approximate GP emulators to balance emulator accuracy with computational speed to facilitate improved UNEDF1 calibration. Sparse GPs are straightforward to train and accurate. Calibration is not straightforward with MCMC (using MH or HMC/NUTS). We produced reusable software for continuing and building on this work as well as accessing and using Darwin cluster compute resources

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Ab Initio Study of Electronic and Magnetic Properties in ZnO-Doped and Co-doped by Vanadium and Silver

Based upon the ab initio of spin density functional calculation, the electronic and magnetic properties of Zn{sub 1−x}M{sub x}O (M = Ag or V) and Zn{sub 1−x−y}A{sub x}B{sub y}O (A = V; B = Ag) are investigated, using the Korringa-Kohn-Rostoker (KKR) method coupled with the coherent potential approximation (CPA). The total and partial densities of state are calculated. The effect of concentration values in Zn{sub 1−x}V{sub x}O and Zn{sub 1−x−y}V{sub x}Ag{sub y}O are deduced. Moreover, the magnetic disorder local moment (DLM) and the total energy are obtained for different concentration values of doped and co-doped zinc oxide (ZnO).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2023

Magnetic ground state and perturbations of the distorted kagome Ising metal TmAgGe

Here, we present the magnetic orders and excitations of the distorted kagome intermetallic magnet TmAgGe. Using neutron single crystal diffraction we identify the propagation vectors k = ($\frac{1}{2}$ 0 0) and k = (0 0 0) and determine the magnetic structures of the zero-field and magnetic field-induced phases for H along the $\mathcal{a}$ and [–110] crystal directions. We determine the experimental magnetic field- temperature (H, T)-phase diagram and reproduce it by Monte Carlo simulations of an effective spin exchange Hamiltonian for one distorted kagome layer. Our model includes a strong axial single-ion anisotropy and significantly smaller exchange couplings, which span up to the third-nearest neighbors within the layer. Single crystal inelastic neutron scattering (INS) measurements reveal an almost flat, only weakly dispersive mode around 7 meV that we use alongside bulk magnetization data to deduce the crystal-electric field (CEF) scheme for the Tm 3+ ions. Random phase approximation (RPA) calculations based on the determined CEF wave functions of the two lowest quasidoublets enable an estimation of the interlayer coupling that is compatible with the experimental INS spectra. No evidence for low-energy spin waves associated to the magnetic order was found, which is consistent with the strongly Ising nature of the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Global description of $β$ – decay with the axially deformed Skyrme finite-amplitude method: Extension to odd-mass and odd-odd nuclei

Here we use the finite-amplitude method (FAM), an efficient implementation of the quasiparticle random phase approximation, to compute $β$-decay rates with Skyrme energy-density functionals for 3983 nuclei, essentially all the medium-mass and heavy isotopes on the neutron-rich side of stability. We employ an extension of the FAM that treats odd-mass and odd-odd nuclear ground states in the equal filling approximation. Our rates are in reasonable agreement both with experimental data where available and with rates from other global calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Bayesian and variational Bayesian approaches for flows in heterogeneous random media

In this paper, we study porous media flows in heterogeneous stochastic media. We propose an efficient forward simulation technique that is tailored for variational Bayesian inversion. As a starting point, the proposed forward simulation technique decomposes the solution into the sum of separable functions (with respect to randomness and the space), where each term is calculated based on a variational approach. This is similar to Proper Generalized Decomposition (PGD). Next, we apply a multiscale technique to solve for each term (as in ) and, further, decompose the random function into 1D fields. As a result, our proposed method provides an approximation hierarchy for the solution as we increase the number of terms in the expansion and, also, increase the spatial resolution of each term. We use the hierarchical solution distributions in a variational Bayesian approximation to perform uncertainty quantification in the inverse problem. We conduct a detailed numerical study to explore the performance of the proposed uncertainty quantification technique and show the theoretical posterior concentration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Ponderomotive and weakly relativistic self-focusing of Gaussian laser beam in plasma: Effect of light absorption

This paper presents an influence of light absorption on self-focusing of laser beam propagation in plasma. The differential equation for beam-width parameter is obtained using the Wentzel-Kramers-Brillouin and paraxial approximations through parabolic equation approach. The nonlinearity in dielectric function is assumed to be aroused due to the combined effect of weakly relativistic and ponderomotive regime. To highlight the nature of propagation, behavior of beam-width parameter with dimensionless distance of propagation is presented graphically and discussed. The present work is helpful to understand issues related to the beam propagation in laser plasma interaction experiments where light absorption plays a vital role.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Characterizing Fluence as a Function of Areal Density for 1 and 14 MeV Neutrons Transported Through Air

Computing neutron dosimetry involves transporting neutrons [n] from a source to a sink, determining the neutron heating within the sink and multiplying the neutron heating by a quality factor to compute the dose equivalent. The space and energy dependent neutron fluence is a measure of the amount of neutrons crossing a particular surface area, about a particular energy and is given in units of [n/cm{sup 2}]. Microscopic cross sections are given in units of area and vary by isotope and reaction type as a function of energy. The product of the atom density of a particular atom in a specie and the microscopic cross section of a particular reaction results in the macroscopic cross section of a particular reaction, given in units of per length [1/cm] as a function of energy. Summing the macroscopic cross sections of all isotopes, and all collisions resulting in heating, within a volume, results in a total collision heating macroscopic cross section for that volume. The product of the energy and spatial dependent collision heating macroscopic cross section and the energy and spatial dependent fluence, within an energy group, results in a collision heating reaction density for that particular reaction in that energy group. Weighting each energy group's collision heating reactions, with the corresponding energy dependent particle damage per collision heating event (usually biological damage), and then integrating over energy and volume results in the delivered dose. Therefore understanding dose delivery involves examining the spatial and energy (and time) dependent fluence behavior. Because of the coupled energy and space dependence of the heating weighting and macroscopic cross sections, the fluence is not so easily separable in energy and space; however, the energy integrated fluence can be used as a first order approximation of the relative magnitude of reactions as a function of space. Here we examine the first part in the dosimetry calculation process, transporting neutrons from source to sink, as applied to neutron transport in air, examining characteristics of the fluence. Presented here is analysis regarding the nature of the fluence peak between 10 g/cm{sup 2} and 100 g/cm{sup 2} for a typical air transport calculation. Due to the high concentration of nitrogen and oxygen isotopes of these elements, these isotopes tend to dominate the transport. From the analysis presented here, the fluence peak is the result of scatter back and forth across the peak regime, and (n,2n) has a minimal (<4%) contribution. Further evidence of minimal (n,2n) impact in the 10-100 g/cm{sup 2} span includes the 1 MeV fluence peak being so much larger than the 14 MeV peak; where 1 MeV is not above the (n,2n) threshold. Also presented here is the usefulness of using tally tagging with cell flagging to further dissect the nature of the fluence. These technologies allow us to understand the originating volume that makes a contribution to a scoring regime. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Energy release rate approximation for edge cracks using higher-order topological derivatives

Topological derivatives provide the variation of a functional when an infinitesimal hole is introduced into the domain. In Silva et al. (J Mech Phys Solids 59(5):925–939, 2011), the authors developed a first-order approximation of the energy release rate associated with a small crack at any boundary location and at any orientation using a first-order topological derivative. This approach offers significant computational advantages over other methods because (i) it requires only a single analysis while other methods require an analysis for each crack size-location-direction combination, and (ii) it is performed on the non-cracked domain, removing the need to create very fine meshes in the vicinity of the crack. In the present study, higher-precision approximations of the energy release rate are developed using higher-order topological derivatives which allow the analyst to accurately treat longer cracks and determine the crack lengths for which the first-order approximation is accurate.

42 ENGINEERING↗