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

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

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

Exchange–correlation entropy from the generalized thermal adiabatic connection

Warm dense matter is a highly energetic phase characterized by strong correlations, thermal effects, and quantum mechanical electrons. Thermal density functional theory is commonly used in simulations of this challenging phase, driving the development of temperature-dependent approximations to the exchange–correlation free energy. Approaches using the adiabatic connection formula are well known at zero temperature and have been recently leveraged at non-zero temperatures as well. In this work, a generalized thermal adiabatic connection (GTAC) formula is proposed, introducing a fictitious temperature parameter. This allows extraction of the exchange–correlation entropy S XC using simulated interaction strength scaling. Here, this procedure uses a Hellmann–Feynman approach to express the exchange–correlation entropy in terms of a temperature- and interaction strength-dependent exchange–correlation potential energy. In addition, analysis of S XC as a function of interaction strength suggests new forms for approximations, and GTAC itself offers a new framework for exploring both the exact and approximate interplay of temperature, density, and interaction strength across a wide range of conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

CGC for ultra-peripheral Pb+Pb collisions at the Large Hadron Collider: a more realistic calculation

We provide the first calculation of two-gluon production at mid-rapidity in ultra-peripheral collisions in the Color Glass Condensate framework. To estimate systematic uncertainty associated with poor understanding of the wave function of the nearly real photon, we consider two diametrically different models: the dilute quark-antiquark dipole approximation and a vector meson, in which color charge density is approximated by McLerran-Venugopalan model. In the experimentally relevant range, the target nucleus can be faithfully approximated by a highly saturated state. This simplification enables us to perform efficient numerical simulations and extract the two-gluon correlation functions and the associated azimuthal harmonics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Elastic and resonance structures of the nucleon from the hadronic tensor in lattice QCD: Implications for neutrino-nucleon scattering and hadron physics

We compute the Euclidean hadronic tensor from charge density operators and extract elastic and resonance structures by employing exponential fits to the four-point function correlator, as well as a Bayesian reconstruction inverse algorithm to obtain the corresponding spectral density for qualitative comparison. We present the determination of the nucleon’s Sachs electric form factor using the hadronic tensor formalism and verify that it is consistent with that from the conventional three-point function calculation. Beyond the elastic peak, we observe a structure located approximately 0.5–0.7 GeV above the nucleon mass in the Bayesian reconstruction. This structure is interpreted as a mixture of the Roper resonance [𝑁⁡(1440)], and states with both positive and negative parities in this mass region, as well as multihadron states. Assuming the observed structure is dominated by 𝐽 𝑃 = 1/2 ± states, we extract the transition electric form factor 𝐺$^*_𝐸$⁡(𝑄 2 ) and the corresponding longitudinal helicity amplitude 𝑆 1/2 ⁡(𝑄 2 ), and compare them with those determined from the CLAS experimental data of nucleon-to-Roper transition. Although fitting to the four-point correlation function or using the inverse algorithm does not resolve individual resonances, it nevertheless enables the determination of total inclusive lepton–nucleon scattering cross sections in appropriate energy bins. This lattice QCD calculation presents the first major step toward studying the inclusive 𝑁 → 𝑋 contributions within the hadronic tensor formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Investigation of the Thermodynamic and Electronic Properties of Double Perovskite Ca{sub 2}CoNbO{sub 6}

In the present work, a self-consistent ab initio calculation using the full-potential linearized augmented plane wave (FP-LAPW) method within the framework of the spin-polarized density functional theory (DFT) was used to study the structural, electronic, and thermodynamic properties of Ca{sub 2}CoNbO{sub 6} double perovskite compound. The generalized gradient approximation (GGA) described by Perdew–Burke–Ernzerhof (PBE) and GGA + U were used. The results obtained for the electronic properties show a ferrimagnetic and half-metallic behavior of the compound. The novelty of our work is the study of the thermodynamic properties of Ca{sub 2}CoNbO{sub 6} double perovskite such as heat capacity and Debye temperature which showed an important effect of pressure compared to the temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Pd(II) and Rh(I) Catalytic Precursors for Arene Alkenylation: Comparative Evaluation of Reactivity and Mechanism Based on Experimental and Computational Studies

We combine experimental and computational investigations to compare and understand catalytic arene alkenylation using the Pd(II) and Rh(I) precursors Pd(OAc) 2 and [(η 2 -C 2 H 4 ) 2 Rh(µ-OAc)] 2 with arene, olefin and Cu(II) carboxylate at elevated temperature (> 120 °C). Under specific conditions, previous computational and experimental efforts have identified heterotrimetallic cyclic PdCu 2 (η 2 -C 2 H 4 ) 3 (μ-OPiv) 6 and [(η 2 -C 2 H 4 ) 2 Rh(µ-OPiv) 2 ] 2 (µ-Cu) (OPiv = pivalate) species as likely active catalysts for these processes. Further studies of catalyst speciation suggest a complicated equilibrium between Cu(II)-containing complexes containing one Rh or Pd atom with complexes containing two Rh or Pd atoms. At 120 °C, Rh catalysis produces styrene > 20-fold more rapidly than Pd. Also, at 120 °C, Rh is ~98% selective for styrene formation while Pd is ~82% selective. Furthermore, our studies indicate that Pd catalysis has a higher predilection toward olefin functionalization to form undesired vinyl ester, while Rh catalysis is more selective for arene/olefin coupling. However, at elevated temperatures, Pd converts vinyl ester and arene to vinyl arene, which is proposed to occur through low valent Pd(0) clusters that are formed in situ. Regardless of arene functionality, the regioselectivity for alkenylation of mono-substituted arenes with the Rh catalyst gives an approximate 2:1 meta:para ratio with minimal ortho C–H activation. In contrast, Pd selectivity is significantly influenced by arene electronics with electron-rich arenes giving an approximate 1:2:2 ortho:meta:para ratio while the electron deficient (α,α,α)-trifluorotoluene gives a 3:1 meta:para ratio with minimal ortho functionalization. Kinetic intermolecular arene ethenylation competition experiments find that Rh reacts most rapidly with benzene, and the rate of mono-substituted arene alkenylation does not correlate with arene electronics. In contrast, with Pd catalysis, electron-rich arenes react more rapidly than benzene while electron-deficient arenes react less rapidly than benzene. These experimental findings, in combination with computational results, are consistent with the arene C–H activation step for Pd catalysis involving significant η 1 -arenium character due to Pd-mediated electrophilic aromatic substitution character. In contrast, the mechanism for Rh catalysis is not sensitive to arene substituent electronics, which we propose indicates less electrophilic aromatic substitution character for the Rh-mediated arene C–H activation.

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

Paraxial, Thin-Lens Analysis of Fixed-Tune, Non-Scaling FFAs with Two Magnets per Cell

Fixed-field accelerators (FFAs) have closed orbits that change as a function of beam momentum. It is sometimes useful to avoid various resonances by keeping the tunes constant as these orbits change. A well-known example is the scaling FFA where the entire beam orbit and optics are geometrically scaled as a function of momentum. However, this is stricter than necessary: cells with three or more lenses can have fixed tunes even when the focussing strengths in the lenses change. Recently, Dejan Trbojevic has found a pair of nonlinear magnets that produce fixed tunes and fixed beta functions while not being a scaling FFA. This improves on a more approximate solution in. Notably a scaling FFA would require one magnet to be entirely reverse-bending, whereas the Trbojevic solution does not do this and instead resembles an intermediate point between the traditional scaling field profiles and the nonscaling profiles centred around a momentum with equal and positive fields. This suggests there are at least three levels of stringency that one can apply to a fixed-tune FFA design: 1. Fixed cell tunes (in both planes) as a function of momentum; 2. Fixed optics (beta functions) as a function of momentum; 3. Similarity of all orbits via a scaling symmetry law ↔ traditional scaling FFA. This note studies the interesting case #2 above (#3 being fully characterised by the orbit at a single energy) in the simplest possible example: a cell of two thin lenses in the small angle (paraxial) approximation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Sequential Learning of Active Subspaces

This study shows that in recent years, active subspace methods (ASMs) have become a popular means of performing subspace sensitivity analysis on black-box functions. Naively applied, however, ASMs require gradient evaluations of the target function. In the event of noisy, expensive, or stochastic simulators, evaluating gradients via finite differencing may be infeasible. In such cases, often a surrogate model is employed, on which finite differencing is performed. When the surrogate model is a Gaussian process (GP), we show that the ASM estimator is available in closed form, rendering the finite-difference approximation unnecessary. We use our closed-form solution to develop acquisition functions focused on sequential learning tailored to sensitivity analysis on top of ASMs. We also show that the traditional ASM estimator may be viewed as a method of moments estimator for a certain class of GPs. We demonstrate how uncertainty on GP hyperparameters may be propagated to uncertainty on the sensitivity analysis, allowing model-based confidence intervals on the active subspace. Our methodological developments are illustrated on several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Investigating the minimum post-burst energy required to function an exploding bridgewire detonator

The minimum energy required after bridge-burst to function a commercially available RP-1 exploding bridgewire detonator with a gold bridgewire was investigated using two capacitor discharge units (CDUs) making use of state-of-the-art on/off high-voltage high-current switches. Both switches are capable of shutting off the current flow in ≈100ns and the conduction time was varied systematically to limit the duration of electrical energy delivered from the CDU until a detonation go/no-go threshold was reached. It was found that the minimal extra electrical energy required for a soft-fire (a statistically marginal go) was approximately half that required for a hard-fire (a definite go). Furthermore, the measured energy required to function the RP-1 detonator from the 11.5 kV CDU was approximately 20 mJ greater than for the 600 V CDU, and it is postulated that this was due to arcing across the detonator terminals and an associated current shunt that was not efficient at depositing energy conducive to detonator functioning. Measurements of the bridgewire temperatures bursting in air showed that the maximum temperature reached was consistent (≈18kK) when sufficient energy was delivered that the detonator functioned (either soft- or hard-fire), but that the temperature for a sub-threshold bridgewire burst was significantly cooler at approximately 8 kK. Furthermore, it was demonstrated that the intensity of the light emitted at burst scaled with deposited energy. A related finding is that the light emission from the bridgewire peaks some tens of nanoseconds after both peak voltage and peak electrical power.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Two-point sum-rules in three-dimensional Yang-Mills theory

We compute the stress-tensor two-point function in three-dimensional Yang-Mills theory to three-loops in perturbation theory. Using its calculable shape at high momenta, we test the notion that its Borel transform is saturated at low energies by the lowest glueball state(s). This assumption provides relatively stable estimates for the mass of the lightest glueball that we compare with lattice simulations. We also provide estimates for the coupling of the lightest glueball to the stress tensor. Along the way, we comment on the extent that such estimates are non-rigorous. Lastly, we discuss the possibility of applying the sum-rule analysis to two-point functions of higher-spin operators and obtain a crude approximation for the glueball couplings to these operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Theoretical Model Study of Interplay of Coulomb Interaction and Electron-Phonon Interaction in the Thermal Properties of Monolayer Graphene

We propose here a tight-binding (TB) model Hamiltonian for monolayer graphene-on-substrate describing the nearest-neighbor-hopping, on-site Coulomb interaction on the sub-lattices and the electron-phonon interaction under the high-frequency limit of phonon vibration. Applying Lang-Firsov canonical transformation, the electron and phonon systems are decoupled in the atomic Hamiltonian, such that the effective Coulomb interaction and effective nearest-neighbor-hopping integral respectively appear as U-tilde =U−2t{sub 1}λ and t-tilde {sub 1}=t{sub 1}e{sup ((−t{sub 1}λ)/(ω{sub 0}))}, where U, t{sub 1}, λ and ω{sub 0} are respectively Coulomb energy, nearest-neighbor-hopping integral, electron-phonon (e-ph) coupling and phonon frequency. The effective Coulomb interaction in the Hamiltonian is considered within mean-field approximation. The Hamiltonian is solved by Zubarev’s Green’s function technique. The temperature-dependent electronic entropy and specific heat are calculated from the free energy of graphene system and are computed numerically. The temperature-dependent electronic specific heat exhibits a charge gap peak at room temperature arising due to the effect of Coulomb interaction and electron-phonon interaction. The evolution of these peaks in specific heat is investigated by varying the model parameters of the system.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Nonlinear dynamics from linear quantum evolutions

Linear dynamics restricted to invariant submanifolds generally gives rise to nonlinear dynamics. Submanifolds in the quantum framework may emerge for several reasons: one could be interested in specific properties possessed by a given family of states, either as a consequence of experimental constraints or inside an approximation scheme. In this work we investigate such issues in connection with a one parameter group ϕ{sub t} of transformations on a Hilbert space, H, defining the unitary evolutions of a chosen quantum system. Two procedures will be presented: the first one consists in the restriction of the vector field associated with the Schrödinger equation to a submanifold invariant under the flow ϕ{sub t}. The second one makes use of the Lagrangian formalism and can be extended also to non-invariant submanifolds, even if in such a case the resulting dynamics is only an approximation of the flow ϕ{sub t}. Such a result, therefore, should be conceived as a generalization of the variational method already employed for stationary problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Parameterization of the Stoner-Wohlfarth model of magnetic hysteresis

The Stoner-Wohlfarth (SW) is the most used model of magnetic hysteresis, but its computation is time-consuming. We approximate piecewise this model by easy-to-compute analytic functions. Overall, our parametrization is suitable for fast quantitative evaluations and fitting experimental data, which we exemplify. Applicability of this parametrization is as broad as that of the SW model, which can be extended to include more parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

Stabilizing lattice oxygen redox in layered sodium transition metal oxide through spin singlet state

Abstract Reversible lattice oxygen redox reactions offer the potential to enhance energy density and lower battery cathode costs. However, their widespread adoption faces obstacles like substantial voltage hysteresis and poor stability. The current research addresses these challenges by achieving a non-hysteresis, long-term stable oxygen redox reaction in the P3-type Na 2/3 Cu 1/3 Mn 2/3 O 2 . Here we show this is accomplished by forming spin singlet states during charge and discharge. Detailed analysis, including in-situ X-ray diffraction, shows highly reversible structural changes during cycling. In addition, local CuO 6 Jahn-Teller distortions persist throughout, with dynamic Cu-O bond length variations. In-situ hard X-ray absorption and ex-situ soft X-ray absorption study, along with density function theory calculations, reveal two distinct charge compensation mechanisms at approximately 3.66 V and 3.99 V plateaus. Notably, we observe a Zhang-Rice-like singlet state during 3.99 V charging, offering an alternative charge compensation mechanism to stabilize the active oxygen redox reaction.

25 ENERGY STORAGE↗
DOE OSTI2024

Structural diversity and clustering of bacterial flagellar outer domains

Supercoiled flagellar filaments function as mechanical propellers within the bacterial flagellum complex, playing a crucial role in motility. Flagellin, the building block of the filament, features a conserved inner D0/D1 core domain across different bacterial species. In contrast, approximately half of the flagellins possess additional, highly divergent outer domain(s), suggesting varied functional potential. In this study, we report atomic structures of flagellar filaments from three distinct bacterial species: Cupriavidus gilardii , Stenotrophomonas maltophilia , and Geovibrio thiophilus . Our findings reveal that the flagella from the facultative anaerobic G. thiophilus possesses a significantly more negatively charged surface, potentially enabling adhesion to positively charged minerals. Furthermore, we analyze all AlphaFold predicted structures for annotated bacterial flagellins, categorizing the flagellin outer domains into 682 structural clusters. This classification provides insights into the prevalence and experimental verification of these outer domains. Remarkably, two of the flagellar structures reported herein belong to a distinct cluster, indicating additional opportunities on the study of the functional diversity of flagellar outer domains. Our findings underscore the complexity of bacterial flagellins and open up possibilities for future studies into their varied roles beyond motility.

Science & Technology - Other Topics↗
DOE OSTI2019

Ferroelectric barium titanate derivatives containing Mo and Mg for transparent photovoltaic applications

Using first-principles methods, we investigate the electronic properties of the [Ba(Mo 1/2 ,Mg 1/2 )O 3 ] x -[BaTiO 3 ] 1-x solid solution derived from barium titanate as a potential candidate to be used in photovoltaic devices. Focusing on the bandgap and its origin, we study the effect of different possible Mo and Mg contents, arrangements, and phases of [Ba(Mo 1/2 ,Mg 1/2 )O 3 ]x-[BaTiO 3 ] 1-x . We find that [Ba(Mo 1/2 ,Mg 1/2 )O 3 ] 0.25 -[BaTiO 3 ] 0.75 is a viable candidate for use in transparent photovoltaics due to its energy bandgap of 2.6 eV in the rhombohedral phase. In all cases, [Ba(Mo 1/2 ,Mg 1/2 )O 3 ]x-[BaTiO 3 ] 1-x materials exhibit spontaneous polarization that allows the exploitation of the bulk photovoltaic effect and in principle may allow high power conversion efficiency exceeding the Shockley-Queisser limit for these materials.

14 SOLAR ENERGY↗
DOE OSTI2020

ABINIT: Overview and focus on selected capabilities

ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗