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,387 records · Page 77

Variational Discrete Action Theory

In this work, we propose the variational discrete action theory (VDAT) to study the ground state properties of quantum many-body Hamiltonians. VDAT is a variational theory based on the sequential product density matrix (SPD) ansatz, characterized by an integer $\mathscr{N}$, which monotonically approaches the exact solution with increasing $\mathscr{N}$. To evaluate the SPD, we introduce a discrete action and a corresponding integer time Green’s function. We use VDAT to exactly evaluate the SPD in two canonical models of interacting electrons: the Anderson impurity model and the d = ∞ Hubbard model. For the latter, we evaluate $\mathscr{N}$ = 2 – 4, where $\mathscr{N}$ = 2 recovers the Gutzwiller approximation (GA), and we show that $\mathscr{N}$ = 3, which exactly evaluates the Gutzwiller-Baeriswyl wave function, provides a truly minimal yet precise description of Mott physics with a cost similar to that of the GA. VDAT is a flexible theory for studying quantum Hamiltonians, competing both with state-of-the-art methods and simple, efficient approaches all within a single framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optical and electronic functionality arising from controlled defect formation in nanoscale complex oxide lateral epitaxy

Epitaxial crystallization of complex oxides provides the means to create materials with precisely selected composition, strain, and orientation, thereby controlling their functionalities. Extending this control to nanoscale three-dimensional geometries can be accomplished via a three-dimensional analog of oxide solid-phase epitaxy, lateral epitaxial crystallization. The orientation of crystals within laterally crystallized SrTiO 3 systematically changes from the orientation of the SrTiO 3 substrate. This evolution occurs as a function of lateral crystallization distance, with a rate of approximately 50° μm -1 . The mechanism of the rotation is consistent with a steady-state stress of tens of megapascal over a 100–nanometer scale region near the moving amorphous/crystalline interface arising from the amorphous-crystalline density difference. Second harmonic generation and piezoelectric force microscopy reveal that the laterally crystallized SrTiO 3 is noncentrosymmetric and develops a switchable piezoelectric response at room temperature, illustrating the potential to use lateral crystallization to control the functionality of complex oxides.

36 MATERIALS SCIENCE↗

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗

The role of AdhE mutations in Thermoanaerobacterium saccharolyticum

ABSTRACT Thermoanaerobacterium saccharolyticum is a thermophilic anaerobic bacterium that natively ferments a variety of hemicellulose substrates to organic acids and alcohols. It has recently been engineered to produce ethanol at high yield and titer; however, it uses a unique metabolic pathway for ethanol production that is poorly characterized. One of the distinctive aspects of this pathway is the presence of acetyl-CoA as an intermediate metabolite. In this organism, acetyl-CoA is converted to ethanol by a bifunctional AdhE enzyme. This enzyme has been a frequent target for mutations, and in many cases, the function of these mutations was unknown. Using a combination of genetic modifications, enzyme assays, and computational analysis, we have developed a better understanding of how mutations in AdhE affect ethanol production in the engineered homoethanologen strain. We identify a set of approximately interchangeable AdhE mutations (G544D, T597K, T597I, and T605I), whose function is to disrupt the activity of the alcohol dehydrogenase (ADH) domain of AdhE. This reduces NADH-linked ADH activity, which dramatically increases ethanol tolerance and changes the overall stoichiometry of acetaldehyde to ethanol conversion. Furthermore, our improved understanding of the function of these AdhE mutations calls into question a proposed feature of AdhE enzymes known as substrate channeling—direct transfer of acetaldehyde between the two domains of the AdhE enzyme. This improved the understanding of the role of AdhE mutations in T. saccharolyticum and provides deeper insights into the function of the unique ethanol production pathway in this organism. IMPORTANCE Many anaerobic bacteria maintain redox equilibrium by producing reduced organic compounds such as ethanol. The final two steps of ethanol production are mediated by a bifunctional enzyme, AdhE, and this enzyme is a frequent target of mutations in strains engineered for increased ethanol production. Paradoxically, these mutations increase ethanol production by eliminating the activity of one domain of the AdhE enzyme (the ADH domain). This provides additional support for a redox-imbalance theory of alcohol tolerance, which challenges the prevailing hypothesis that alcohol tolerance is associated with cell membrane effects.

59 BASIC BIOLOGICAL SCIENCES↗

Perspective: Explicitly correlated electronic structure theory for complex systems

The explicitly correlated approach is one of the most important breakthroughs in ab initio electronic structure theory, providing arguably the most compact, accurate, and efficient ansatz for describing the correlated motion of electrons. Since Hylleraas first used an explicitly correlated wave function for the He atom in 1929, numerous attempts have been made to tackle the significant challenges involved in constructing practical explicitly correlated methods that are applicable to larger systems. These include identifying suitable mathematical forms of a correlated wave function and an efficient evaluation of many-electron integrals. R12 theory, which employs the resolution of the identity approximation, emerged in 1985, followed by the introduction of novel correlation factors and wave function ansätze, leading to the establishment of F12 theory in the 2000s. Rapid progress in recent years has significantly extended the application range of explicitly correlated theory, offering the potential of an accurate wave-function treatment of complex systems such as photosystems and semiconductors. In conclusion, this perspective surveys explicitly correlated electronic structure theory, with an emphasis on recent stochastic and deterministic approaches that hold significant promise for applications to large and complex systems including solids.

74 ATOMIC AND MOLECULAR PHYSICS↗

Variational treatment of electron–polyatomic-molecule scattering calculations using adaptive overset grids

In this paper, the complex Kohn variational method for electron–polyatomic-molecule scattering is formulated using an overset-grid representation of the scattering wave function. The overset grid consists of a central grid and multiple dense atom-centered subgrids that allow the simultaneous spherical expansions of the wave function about multiple centers. Scattering boundary conditions are enforced by using a basis formed by the repeated application of the free-particle Green's function and potential G ^ 0 + V ^ on the overset grid in a Born-Arnoldi solution of the working equations. The theory is shown to be equivalent to a specific Padé approximant to the T matrix and has rapid convergence properties, in both the number of numerical basis functions employed and the number of partial waves employed in the spherical expansions. The method is demonstrated in calculations on methane and CF 4 in the static-exchange approximation and compared in detail with calculations performed with the numerical Schwinger variational approach based on single-center expansions. Finally, an efficient procedure for operating with the free-particle Green's function and exchange operators (to which no approximation is made) is also described.

74 ATOMIC AND MOLECULAR PHYSICS↗

Affordable and accurate large-scale hybrid-functional calculations on GPU-accelerated supercomputers

Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding or even out of reach if high quality basis sets are used. We present a highly optimized multiple graphics processing unit implementation of the exact exchange operator which allows one to perform fast hybrid functional density-functional theory (DFT) calculations with systematic basis sets without additional approximations for up to a thousand atoms. With this method hybrid DFT calculations of high quality become accessible on state-of-the-art supercomputers within a time-to-solution that is of the same order of magnitude as traditional semilocal-GGA functionals. The method is implemented in a portable open-source library.

density functional theory↗

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Recovery of the characteristics of plasma turbulence from the radial correlation backscattering diagnostics

Signals of the backscattering radial correlation Doppler diagnostics of plasma density fluctuations in the presence of the cutoff of the probing wave are analyzed theoretically with allowance for the curvature of magnetic surfaces. The scattering of the probing electromagnetic wave is considered in the linear (Born) approximation with respect to the amplitude of fluctuations. Using the Wentzel−Kramers−Brillouin approach, analytical expressions for the scattered signal and the correlation function of two scattered signals corresponding to oblique probing at different frequencies are derived. A criterion is obtained for the tilt angle of the antenna pattern at which the two-point turbulence correlation function can be measured directly. A method is proposed to recover the spectrum of plasma density fluctuations from the data on the radial wavenumbers even if this criterion is violated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The integrated bispectrum in modified gravity theories

Gravity-induced non-Gaussianity can provide important clues to Modified Gravity (MG) Theories. Several recent studies have suggested using the Integrated Bispectrum (IB) as a probe for squeezed configuration of bispectrum. Extending previous studies on the IB, we include redshift-space distortions to study a class of (parametrised) MG theories that include the string-inspired Dvali, Gabadadze and Porrati (DGP) model. Various contributions from redshift-space distortions are derived in a transparent manner, and squeezed contributions from these terms are derived separately. Results are obtained using the Zel'dovich Approximation (ZA). Results are also presented for projected surveys (2D). We use the Press-Schechter (PS) and Sheth-Tormen (ST) mass functions to compute the IB for collapsed objects that can readily be extended to peak-theory based approaches. The cumulant correlators (CCs) generalise the ordinary cumulants and are known to probe collapsed configurations of higher order correlation functions. We generalise the concept of CCs to halos of different masses. We also introduce a generating function based approach to analyse more general non-local biasing models. The Fourier representations of the CCs, the skew-spectrum, or the kurt-spctra are discussed in this context. The results are relevant for the study of the Minkowski Functionals (MF) of collapsed tracers in redshift-space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exfoliation energy, quasiparticle band structure, and excitonic properties of selenium and tellurium atomic chains

Effects that are not captured by the generalized-gradient density-functional theory have a prominent effect on the structural binding and on the electronic and optical properties of reduced-dimensional and weakly bound materials. Here, we report the exfoliation energy of selenium and tellurium atomic chains with nonempirical van der Waals corrections, and their electronic and optical properties with the GW and Bethe-Salpeter formalisms. The exfoliation energy is found to be within 0.547–0.719 eV/u.c. for the selenium atomic chain, and 0.737–0.926 eV/u.c. for the tellurium atomic chain (where u.c. stands for unit cell), depending on the approximation for the van der Waals interaction and the numerical tool chosen. The GW electronic band gap turned out to be 5.22–5.47 (4.44–4.59) eV for the Se (Te) atomic chains, with the lowest bound obtained with the Godby-Needs (GN), and the upper bound to the Hybertsen-Louie (HL) plasmon-pole models (PPMs). The binding energy of the ground-state excitonic state ranges between 2.69 and 2.72 eV for selenium chains within the HL and GN PPM, respectively, and it turned out to be 2.35 eV for tellurium chains with both approximations. In conclusion, the ground-state excitonic wave function is localized within 50 Å along the axis for both types of atomic chains, and its energy lies within the visible spectrum: blue [2.50(GN)–2.78(HL) eV] for selenium and yellow-green [2.09(GN)–2.28(HL) eV] for tellurium, which could be useful for LED applications in the visible spectrum.

74 ATOMIC AND MOLECULAR PHYSICS↗

Ab initio study of the composite phase diagram of Ni–Mn–Ga shape memory alloys

The magnetic and structural properties of a series of nonstoichiometric Ni–Mn–Ga Heusler alloys are theoretically investigated in terms of the density functional theory. Nonstoichiometry is formed in the coherent potential approximation. Concentration dependences of the equilibrium lattice parameter, the bulk modulus, and the total magnetic moment are obtained and projected onto the ternary phase diagram of the alloys. The stable crystalline structures and the magnetic configurations of the austenitic phase are determined.

36 MATERIALS SCIENCE↗

Examining the Effect of Exchange-Correlation Approximations in First-Principles Dynamics Simulation of Interfacial Charge Transfer

We examine the extent to which the exchange-correlation (XC) approximation influences modeling interfacial charge transfer using fewest-switches surface hopping (FSSH) simulations within the single-particle description. A heterogeneous interface between a lithium ion and an extended boron-nitride sheet was considered here, being an extreme case in which wave function localization and energy level alignments are highly sensitive to the XC approximation. The PBE0 hybrid XC approximation yields nonadiabatic couplings (NACs) that are significantly smaller than the values obtained from the PBE-GGA approximation by an order of magnitude for localized electronic states. This difference between the two XC functionals for the calculated NACs was found to derive mainly from the wave function characteristics rather than from the lattice movement although first-principles molecular dynamics trajectories, along which NACs are obtained, differ noticeably between the two XC functionals. Using the NACs and single-particle energy level alignments at different levels of theory, FSSH simulations were performed to model the electron transfer dynamics at the interface. Here, the electron transfer time scale was found to vary as much as, but not more than, 1 order of magnitude. The time scale was found to be quite sensitive to both NACs and energy level alignments. While the order of magnitude consistency for the charge transfer rate is encouraging even for this rather extreme model of heterojunction interface, continued advancement in electronic structure methods is required for quantitatively accurate determination of the transfer rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of the Reactive CO 2 Surface Flux in Electrocatalytic Aqueous Flow Reactors

We study how the mass transfer of a chemical species in shear flow is suppressed by the production of a second species that partially reacts away the first species. The second species is produced at the surface with the first species as a reactant. Our work is directly motivated by electrochemical CO 2 reduction in aqueous flow reactors, where OH– molecules generated by the CO evolution reaction react away CO 2 molecules, ultimately inhibiting the mass transfer of CO 2 to the cathode surface. We derive a simple approximation of the Sherwood number, a dimensionless measure of the mass flux into the surface, as a function of the Péclet number, surface Damköhler number, and bulk Damköhler number, and we compare our approximation to numerical solutions of the governing equations. We find that in the diffusion-limited regime, the Sherwood number is well described by the classical Graetz-Lévêque result reduced by a reaction factor due to the competing bulk reaction; with the stoichiometry relevant to electrochemical CO 2 reduction, this reaction factor is 1/2. While the model problem we solve provides insight into how OH– production affects CO 2 mass transfer, experimental systems often involve more complex chemistry. We thus also show how a common buffered electrolyte, KHCO 3 , affects our results.

30 DIRECT ENERGY CONVERSION↗

Quark Sivers Function at Small-$x$: Leading contribution from the Spin-Dependent Odderon

We present the calculation of the leading contribution to the quark Sivers function at small-Bjorken x as in [1]. This calculation uses the high energy scattering approximation and operator formalism developed in [2,3] to obtain a dominant contribution to the quark Sivers function coming from the spin-dependent odderon, in agreement with the results of [4]. We then calculate this dominant contribution in the diquark model of the proton to obtain a small-x estimate for the Sivers function.

Santiago, M. Gabriel↗

Resonance Self-Shielding Using Deterministic Transport Methods with Pointwise Cross Sections (Invited) - Eugene P. Wigner Reactor Physics Award

The difficulty of predicting resonance absorption has been recognized since the birth of reactor physics. A paper published by Wigner, et al. in 1955 describes pioneering theoretical work done in 1941 to provide an estimate for the resonance absorption in heterogeneous uranium lattices. This work introduced many fundamental concepts, including the narrow resonance (NR) approximation and the Wigner rational expression for the escape probability from a fuel lump, and it established an equivalence relation for resonance absorption in homogeneous and heterogeneous systems. During the 1960's, Bondarenko proposed an approach in which resonance shielding factors are precomputed and parametrized as a function of the background cross section by using the NR approximation and equivalence theory. This provided an efficient technique for treating resonance absorption in reactor lattices. The Bondarenko method and its modern variations (e.g., certain subgroup techniques) have become the predominate method for resonance shielding calculations in reactor physics. However, this approach has several limitations, most of which are well known, such as the treatment of resonance interference effects and its reliance on homogeneous-heterogeneous equivalence theory relationships. Other limitations are more subtle and less generally recognized, such as neglecting the impact of self-shielding the multigroup (MG) elastic removal cross section and the effects of directionally dependent shielding on the total cross section. This paper and the accompanying presentation apply deterministic solutions of the continuous energy (CE) neutron transport equation to examine common approximations made in resonance methods for a uniform lattice. The calculated fine-structure spectra are used to calibrate the impact of the approximations and describe how results from CE flux solution can be incorporated in the Bondarenko self-shielding method to improve the accuracy. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

First-Principle Study of the Electronic and Magnetic Properties of Nd{sub 1−x}La{sub x}FeSi Alloys (x = 0, 0.25, 0.50, 0.75, and 1)

The structural, electronic, and magnetic properties of Nd{sub 1−x}La{sub x}FeSi (x = 0, 0.25, 0.50, 0.75, and 1) alloys are investigated based on the first-principle density functional theory by using the full-potential linear muffin-tin orbital (FP-LMTO) method. The exchange-correlation potential is treated with local density approximation (LDA). In addition, to treat localization of 4f-electrons properly, the LDA + U approach (where U is the Hubbard correlation term) is also employed. The calculated structural, electronic, and magnetic properties of these alloys within LDA and LDA + U approaches are compared. The calculated lattice constants of these alloys as a function of x within the LDA + U approach are in the best agreement with linear Vegard’s rule. Most importantly, we focus on the magnetic properties of these alloys and find that the magnetic interaction between Nd atoms is the Rudermen-Kittel-Kasuya-Yosida (RKKY) interaction type. Iron does not carry magnetic moment in these alloys and has just an induced small magnetic moment from Nd atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

History-Dependent Problems with Applications to Contact Models for Elastic Beams

We prove an existence and uniqueness result for a class of subdifferential inclusions which involve a history-dependent operator. Then we specialize this result in the study of a class of history-dependent hemivariational inequalities. Problems of such kind arise in a large number of mathematical models which describe quasistatic processes of contact. To provide an example we consider an elastic beam in contact with a reactive obstacle. The contact is modeled with a new and nonstandard condition which involves both the subdifferential of a nonconvex and nonsmooth function and a Volterra-type integral term. We derive a variational formulation of the problem which is in the form of a history-dependent hemivariational inequality for the displacement field. Then, we use our abstract result to prove its unique weak solvability. Finally, we consider a numerical approximation of the model, solve effectively the approximate problems and provide numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗