Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis function”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Efficient Ensemble-Based Stochastic Gradient Methods for Optimization Under Geological Uncertainty

Ensemble-based stochastic gradient methods, such as the ensemble optimization (EnOpt) method, the simplex gradient (SG) method, and the stochastic simplex approximate gradient (StoSAG) method, approximate the gradient of an objective function using an ensemble of perturbed control vectors. These methods are increasingly used in solving reservoir optimization problems because they are not only easy to parallelize and couple with any simulator but also computationally more efficient than the conventional finite-difference method for gradient calculations. In this work, we show that EnOpt may fail to achieve sufficient improvement of the objective function when the differences between the objective function values of perturbed control variables and their ensemble mean are large. On the basis of the comparison of EnOpt and SG, we propose a hybrid gradient of EnOpt and SG to save on the computational cost of SG. We also suggest practical ways to reduce the computational cost of EnOpt and StoSAG by approximating the objective function values of unperturbed control variables using the values of perturbed ones. We first demonstrate the performance of our improved ensemble schemes using a benchmark problem. Results show that the proposed gradients saved about 30–50% of the computational cost of the same optimization by using EnOpt, SG, and StoSAG. As a real application, we consider pressure management in carbon storage reservoirs, for which brine extraction wells need to be optimally placed to reduce reservoir pressure buildup while maximizing the net present value. Results show that our improved schemes reduce the computational cost significantly.

58 GEOSCIENCES↗
DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Waste Chemistry Envelope for the Determination of Pitting Corrosion Susceptibility - 19554

Liquid radioactive wastes at the Hanford and Savannah River sites are stored in large underground carbon steel tanks. The high-level wastes are processed in several of the tanks and then transferred by piping to other site facilities for further processing before they are stabilized in a vitrified or grout waste form. Based on waste removal and processing schedules, many of the tanks, will be required to be in service for times exceeding the initial intended life (i.e., greater than 60 years). Until the waste is removed from storage, transferred, and processed, the materials and structures of the tanks must maintain a confinement function by providing a barrier to the environment and by maintaining acceptable structural stability during design basis events, which include loadings from both normal service and abnormal (e.g., earthquake) conditions. Administrative programs are in place to mitigate potential corrosion mechanisms and thereby maintain the structural and leak integrity functions of these waste tanks throughout their intended service life. Waste chemistry and temperature control are the means by which corrosion degradation of the waste tanks is minimized. Chemistry controls are based on in-situ coupon tests, waste tank experience, and laboratory testing. These requirements comprise the technical standards for corrosion control. The standards are utilized then as guidelines for operational procedures. The present waste chemistry controls for the Hanford DSTs were established in the 1980's in response to tank failures caused by stress corrosion cracking (SCC) at the DOE Savannah River Site (SRS). Portions of the SRS chemistry limits, which were developed at the Savannah River National Laboratory (SRNL) were coupled with corrosion testing results from Pacific Northwest National Laboratories (PNNL) and adopted for corrosion control of the double-shell tank (DST) waste at the Hanford Site. Since that time, the temperatures of the waste in the DSTs have decreased to less than 50 deg. C in most cases as the radioactive material has decayed. As a consequence, new chemistry controls were formulated for SCC in 2010. This testing however, did not address that potential for pitting in the supernatant and at the liquid-air interface. Indeed, initial tests indicated that the inhibitor requirements for SCC would be insufficient for the prevention of pitting. In addition to lower waste temperatures, future retrieval and process changes may impact the DST waste chemistry. Initial projections suggest that the waste chemistry of the DSTs may shift to a broader range of pH and higher aggressive anion concentrations (e.g., chloride, sulfate, etc.) than the wastes that are presently stored. Additionally, processes within the vitrification facility may result in streams that are returned to the tank farm for evaporation. These streams may also contain high concentrations of aggressive species. SRNL investigated the propensity for pitting of the carbon steel waste tanks at current and projected DST waste chemistries. A statistically designed series of tests provided a pitting factor that relates the ratio of the inhibitor species to the aggressive species and the probability of observing pitting. An acceptable pitting factor was determined and was utilized for the development of a new chemistry envelope. The new chemistry envelope was compared to the envelope that was established for SCC and was determined to be bounding. Validation of the new chemistry envelope was performed via model simulations and comparison with historical experimental results and approaches to establishing chemistry envelopes. The new chemistry envelope will eventually be part of an overall assessment of the condition of the tanks and any changes made to the overall structural integrity program for the Hanford DSTs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Coupled-cluster based basis sets for valence correlation calculations

Novel basis sets are generated that target the description of valence correlation in atoms H through Ar. The new contraction coefficients are obtained according to the Atomic Natural Orbital (ANO) procedure from CCSD(T) (coupled-cluster singles and doubles with perturbative triples correction) density matrices starting from the primitive functions of Dunning et al. [J. Chem. Phys. 90, 1007 (1989); ibid. 98, 1358 (1993); ibid. 100, 2975 (1993)] (correlation consistent polarized valence X-tuple zeta, cc-pVXZ). The exponents of the primitive Gaussian functions are subject to uniform scaling in order to ensure satisfaction of the virial theorem for the corresponding atoms. These new sets, named ANO-VT-XZ (Atomic Natural Orbital Virial Theorem X-tuple Zeta), have the same number of contracted functions as their cc-pVXZ counterparts in each subshell. The performance of these basis sets is assessed by the evaluation of the contraction errors in four distinct computations: correlation energies in atoms, probing the density in different regions of space via 〈r{sup n}〉 (−3 ≤ n ≤ 3) in atoms, correlation energies in diatomic molecules, and the quality of fitting potential energy curves as measured by spectroscopic constants. All energy calculations with ANO-VT-QZ have contraction errors within “chemical accuracy” of 1 kcal/mol, which is not true for cc-pVQZ, suggesting some improvement compared to the correlation consistent series of Dunning and co-workers.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Polarizable Frozen Density Embedding with External Orthogonalization

In this work, we report a polarizable subsystem density functional theory to describe electronic properties of molecules embedded on a metal cluster. Interaction between the molecule and metal cluster is described using frozen density embedding (FDE). Substituting the nonadditive kinetic potential (NAKP) by approximate functionals is circumvented by enforcing external orthogonality (EO) through a projection operator. The computationally expensive freeze/thaw (FT) cycles are bypassed by including a polarization term in the embedding operator. Furthermore, the combination of polarization and EO permits supermolecular basis set calculations, which was not possible for strongly interacting systems with existing kinetic energy functionals. To test the method, we described the ground state density of pyridine, water, and benzene on a silver cluster. Performing FT on top of EO results in exact density embedding for this category of systems and is thus used for benchmarking the method. We find that the density is reproduced to within 0.15e, and the dipole and quadrupole moments are within 18% of the reference points for subsystem separations ranging from bonding to noninteracting distances. Additionally, our formalism allows the flexibility of incorporating different density functionals to the molecular and the metallic subsystems reducing the overall computational cost.

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

Predicting Surface Energies and Particle Morphologies of Boehmite (y-AlOOH) from Density Functional Theory

Particle morphology is an important property affecting the reactivity of solid phases. The aluminum oxyhydroxide boehmite (?-AlOOH) is one example for which particle morphology has relevance to geochemistry, environmental management, and catalysis. Morphologies can be predicted from surface and attachment energies obtained from electronic structure calculations but the dependence of these predictions on the level of theory and complexity of the physics involved has been overlooked. Density functional theory calculations were therefore performed for a wide range of exchange-correlation (XC) functionals, including a hybrid functional, to calculate the energetics of the four main low-index surfaces of boehmite both in vacuum and with a water monolayer adsorbed. The approximation used for describing core electrons, the size of the plane-wave basis set, and the need for dispersion corrections were also investigated. The results highlight the critical need for converged plane-wave basis sets to obtain accurate surface energies as well as the strong influence the XC functional can have on surface energies and water adsorption energies. The predicted morphologies were compared to particle morphologies obtained from a variety of synthesis routes both carried out in this work and published in the literature. The comparison showed that, at basic pH, growth morphologies dominated at low temperatures and short time periods whereas equilibrium morphologies were achieved as the temperature and/or aging time increased, which could be explained by the slow stacking growth mode of boehmite along the [010] direction. Finally, calorimetry experiments are often used to derive surface energies, but excess water on hydrophilic surfaces can add significant uncertainty. An approach for correcting the surface energetics obtained by calorimetry using electronic structure calculations is presented.

Prange, Micah P.↗
DOE OSTI2012

THE FAINT END OF THE LUMINOSITY FUNCTION AND LOW SURFACE BRIGHTNESS GALAXIES

Smithsonian Hectospec Lensing Survey (SHELS) is a dense redshift survey covering a 4 deg{sup 2} region to a limiting R = 20.6. In the construction of the galaxy catalog and in the acquisition of spectroscopic targets, we paid careful attention to the survey completeness for lower surface brightness dwarf galaxies. Thus, although the survey covers a small area, it is a robust basis for computation of the slope of the faint end of the galaxy luminosity function to a limiting M{sub R} = -13.3 + 5log h. We calculate the faint-end slope in the R band for the subset of SHELS galaxies with redshifts in the range 0.02 {<=}z < 0.1, SHELS{sub 0.1}. This sample contains 532 galaxies with R < 20.6 and with a median surface brightness within the half-light radius of SB{sub 50,R} = 21.82 mag arcsec{sup -2}. We used this sample to make one of the few direct measurements of the dependence of the faint end of the galaxy luminosity function on surface brightness. For the sample as a whole the faint-end slope, {alpha} = -1.31 {+-} 0.04, is consistent with both the Blanton et al. analysis of the Sloan Digital Sky Survey and the Liu et al. analysis of the COSMOS field. This consistency is impressive given the very different approaches of these three surveys. A magnitude-limited sample of 135 galaxies with optical spectroscopic redshifts with mean half-light surface brightness, SB{sub 50,R} {>=} 22.5 mag arcsec{sup -2} is unique to SHELS{sub 0.1}. The faint-end slope is {alpha}{sub 22.5} = -1.52 {+-} 0.16. SHELS{sub 0.1} shows that lower surface brightness objects dominate the faint-end slope of the luminosity function in the field, underscoring the importance of surface brightness limits in evaluating measurements of the faint-end slope and its evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Algebraic function operator expectation value based quantum eigenstate determination: A case of twisted or bent Hamiltonian, or, a spatially univariate quantum system on a curved space

Recent studies on quantum evolutionary problems in Demiralp’s group have arrived at a stage where the construction of an expectation value formula for a given algebraic function operator depending on only position operator becomes possible. It has also been shown that this formula turns into an algebraic recursion amongst some finite number of consecutive elements in a set of expectation values of an appropriately chosen basis set over the natural number powers of the position operator as long as the function under consideration and the system Hamiltonian are both autonomous. This recursion corresponds to a denumerable infinite number of algebraic equations whose solutions can or can not be obtained analytically. This idea is not completely original. There are many recursive relations amongst the expectation values of the natural number powers of position operator. However, those recursions may not be always efficient to get the system energy values and especially the eigenstate wavefunctions. The present approach is somehow improved and generalized form of those expansions. We focus on this issue for a specific system where the Hamiltonian is defined on the coordinate of a curved space instead of the Cartesian one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A new time dependent density functional algorithm for large systems and plasmons in metal clusters

A new algorithm to solve the Time Dependent Density Functional Theory (TDDFT) equations in the space of the density fitting auxiliary basis set has been developed and implemented. The method extracts the spectrum from the imaginary part of the polarizability at any given photon energy, avoiding the bottleneck of Davidson diagonalization. The original idea which made the present scheme very efficient consists in the simplification of the double sum over occupied-virtual pairs in the definition of the dielectric susceptibility, allowing an easy calculation of such matrix as a linear combination of constant matrices with photon energy dependent coefficients. The method has been applied to very different systems in nature and size (from H{sub 2} to [Au{sub 147}]{sup −}). In all cases, the maximum deviations found for the excitation energies with respect to the Amsterdam density functional code are below 0.2 eV. The new algorithm has the merit not only to calculate the spectrum at whichever photon energy but also to allow a deep analysis of the results, in terms of transition contribution maps, Jacob plasmon scaling factor, and induced density analysis, which have been all implemented.

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

Performance of a nonempirical density functional on molecules and hydrogen-bonded complexes

Recently, Tao and Mo derived a meta-generalized gradient approximation functional based on a model exchange-correlation hole. In this work, the performance of this functional is assessed on standard test sets, using the 6-311++G( 3df,3pd ) basis set. These test sets include 223 G3/99 enthalpies of formation, 99 atomization energies, 76 barrier heights, 58 electron affinities, 8 proton affinities, 96 bond lengths, 82 harmonic vibrational frequencies, 10 hydrogen-bonded molecular complexes, and 22 atomic excitation energies. Here our calculations show that the Tao-Mo functional can achieve high accuracy for most properties considered, relative to the local spin-density approximation, Perdew-Burke-Ernzerhof, and Tao-Perdew-Staroverov-Scuseria functionals. In particular, it yields the best accuracy for proton affinities, harmonic vibrational frequencies, hydrogen-bond dissociation energies and bond lengths, and atomic excitation energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Mesoscale Origin of the Enhanced Cycling-Stability of the Si-Conductive Polymer Anode for Li-ion Batteries

Electrode used in lithium-ion battery is invariably a composite of multifunctional components. The performance of the electrode is controlled by the interactive function of all components at mesoscale. Fundamental understanding of mesoscale phenomenon sets the basis for innovative designing of new materials. Here we report the achievement and origin of a significant performance enhancement of electrode for lithium ion batteries based on Si nanoparticles wrapped with conductive polymer. This new material is in marked contrast with conventional material, which exhibit fast capacity fade. In-situ TEM unveils that the enhanced cycling stability of the conductive polymer-Si composite is associated with mesoscale concordant function of Si nanoparticles and the conductive polymer. Reversible accommodation of the volume changes of Si by the conductive polymer allows good electrical contact between all the particles during the cycling process. In contrast, the failure of the conventional Si-electrode is probed to be the inadequate electrical contact.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

How Accurate Is Density Functional Theory at Predicting Dipole Moments? An Assessment Using a New Database of 200 Benchmark Values

We present that dipole moments are a simple, global measure of the accuracy of the electron density of a polar molecule. Dipole moments also affect the interactions of a molecule with other molecules as well as electric fields. To directly assess the accuracy of modern density functionals for calculating dipole moments, we have developed a database of 200 benchmark dipole moments, using coupled cluster theory through triple excitations, extrapolated to the complete basis set limit. This new database is used to assess the performance of 88 popular or recently developed density functionals. The results suggest that double hybrid functionals perform the best, yielding dipole moments within about 3.6-4.5% regularized RMS error versus the reference values - which is not very different from the 4% regularized RMS error produced by coupled cluster singles and doubles. Many hybrid functionals also perform quite well, generating regularized RMS errors in the 5-6% range. Lastly, some functionals, however, exhibit large outliers, and local functionals in general perform less well than hybrids or double hybrids.

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

On the averaging of cardiac diffusion tensor MRI data: the effect of distance function selection

Diffusion tensor magnetic resonance imaging (DT-MRI) allows a unique insight into the microstructure of highly-directional tissues. The selection of the most proper distance function for the space of diffusion tensors is crucial in enhancing the clinical application of this imaging modality. Both linear and nonlinear metrics have been proposed in the literature over the years. The debate on the most appropriate DT-MRI distance function is still ongoing. Here, in this paper, we presented a framework to compare the Euclidean, affine-invariant Riemannian and log-Euclidean metrics using actual high-resolution DT-MRI rat heart data. We employed temporal averaging at the diffusion tensor level of three consecutive and identically-acquired DT-MRI datasets from each of five rat hearts as a means to rectify the background noise-induced loss of myocyte directional regularity. This procedure is applied here for the first time in the context of tensor distance function selection. When compared with previous studies that used a different concrete application to juxtapose the various DT-MRI distance functions, this work is unique in that it combined the following: (i) metrics were judged by quantitative—rather than qualitative—criteria, (ii) the comparison tools were non-biased, (iii) a longitudinal comparison operation was used on a same-voxel basis. The statistical analyses of the comparison showed that the three DT-MRI distance functions tend to provide equivalent results. Hence, we came to the conclusion that the tensor manifold for cardiac DT-MRI studies is a curved space of almost zero curvature. The signal to noise ratio dependence of the operations was investigated through simulations. Finally, the 'swelling effect' occurrence following Euclidean averaging was found to be too unimportant to be worth consideration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

A simple finite element method for the Stokes equations

The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Effects of La Doping and Zn or O Vacancy on the Magnetic Property of ZnO

Research on the magnetic origin and mechanism of La-doped ZnO system has been controversial. To solve this problem, we studied the effects of La doping and Zn or O vacancy on the magnetic property of ZnO using geometry optimization and energy calculation on the basis of first-principle generalized gradient approximation method of density functional theory. The electronic structures and magnetic properties of undoped ZnO, La-doped ZnO system, and systems with coexisting La doping and Zn or O vacancies were calculated. Results showed that the systems of La-doped ZnO and the coexistence of La doping and O vacancy are nonmagnetic. In addition, the system of La replacing Zn and one Zn vacancy exhibits long-range orderly ferromagnetism, and the Curie temperature of the doping system can achieve room temperature. The magnetism source of the systems with La doping and Zn vacancy coexisting in ZnO demonstrates strong hybrid coupling electron exchange effects existing among Zn-4s, Zn-3p, O-2p, and La-5s orbits that are nearest to the Zn vacancy. The next nearest distance between doping and Zn vacancy leads to the lowest formation energy and highest stability under the same doping condition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

molgw 1: Many-body perturbation theory software for atoms, molecules, and clusters

Here, we summarize the MOLGW code that implements density-functional theory and many-body perturbation theory in a Gaussian basis set. The code is dedicated to the calculation of the many-body self-energy within the GW approximation and the solution of the Bethe–Salpeter equation. These two types of calculations allow the user to evaluate physical quantities that can be compared to spectroscopic experiments. Quasiparticle energies, obtained through the calculation of the GW self-energy, can be compared to photoemission or transport experiments, and neutral excitation energies and oscillator strengths, obtained via solution of the Bethe–Salpeter equation, are measurable by optical absorption. The implementation choices outlined here have aimed at the accuracy and robustness of calculated quantities with respect to measurements. Furthermore, the algorithms implemented in MOLGW allow users to consider molecules or clusters containing up to 100 atoms with rather accurate basis sets, and to choose whether or not to apply the resolution-of-the-identity approximation. Finally, we demonstrate the parallelization efficacy of the MOLGW code over several hundreds of processors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗