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,585 records · Page 88

One-particle many-body Green’s function theory: Algebraic recursive definitions, linked-diagram theorem, irreducible-diagram theorem, and general-order algorithms

A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green’s function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green’s function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green’s function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1 ≤ n ≤ 3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green’s function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green’s function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1 ≤ n ≤ 3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An efficient method to integrate polynomials over polytopes and curved solids

Here in this paper, we present an efficient approach to compute the integral of monomials and polynomials over polyhedra and regions defined by parametric curved boundary surfaces. We use Euler's theorem for homogeneous functions in combination with Stokes's theorem to reduce the integration of a monomial over a three-dimensional solid to its boundary. If the solid is a polytope, through a recursive application of these theorems, the integral is further reduced to just the evaluation of the monomial and its derivatives at the vertices of the polytope. The present approach is simpler than existing techniques that rely on repeated use of the divergence theorem, which require the antiderivative of the monomials and the projection of these functions onto hyperplanes. For convex and nonconvex polytopes, our approach does not introduce any approximation for the integration of monomials. For curved solid regions bounded by surfaces that admit a parameterization, the same approach yields simplified formulas to compute the integral of any homogeneous function, including monomials. For surfaces parameterized by polynomial surfaces (such as Bezier surface triangles and B-spline patches), the method yields machine-precision accuracy for the volumetric integration of monomials with an appropriate quadrature rule. Numerical examples over regions bounded by polynomial surfaces and rational surfaces are presented to establish the accuracy and efficiency of the method.

97 MATHEMATICS AND COMPUTING↗

Unsteady fluid flow in a slightly curved pipe: A comparative study of a matched asymptotic expansions solution with a single analytical solution

The present work is motivated by the fact that blood flow in the aorta and the main arteries is governed by large finite values of the Womersley number α and for such values of α there is not any analytical solution in the literature. The existing numerical solutions, although accurate, give limited information about the factors that affect the flow, whereas an analytical approach has an advantage in that it can provide physical insight to the flow mechanism. Having this in mind, we seek analytical solution to the equations of the fluid flow driven by a sinusoidal pressure gradient in a slightly curved pipe of circular cross section when the Womersley number varies from small finite to infinite values. Initially the equations of motion are expanded in terms of the curvature ratio δ and the resulting linearized equations are solved analytically in two ways. In the first, we match the solution for the main core to that for the Stokes boundary layer. This solution is valid for very large values of α. In the second, we derive a straightforward single solution valid to the entire flow region and for 8 ≤ α < ∞, a range which includes the values of α that refer to the physiological flows. Each solution contains expressions for the axial velocity, the stream function, and the wall stresses and is compared to the analogous forms presented in other studies. The two solutions give identical results to each other regarding the axial flow but differ in the secondary flow and the circumferential wall stress, due to the approximations employed in the matched asymptotic expansion process. The results on the stream function from the second solution are in agreement with analogous results from other numerical solutions. The second solution predicts that the atherosclerotic plaques may develop in any location around the cross section of the aortic wall unlike to the prescribed locations predicted by the first solution. In addition, it gives circumferential wall stresses augmented by approximately 100% with respect to the matched asymptotic expansions, a factor that may contribute jointly with other pathological factors to the faster aging of the arterial system and the possible malfunction of the aorta.

42 ENGINEERING↗

An Accurate, Affordable Density Functional Tight-Binding Model for Excited State Hydrocarbon Polymer Molecular Dynamics

We have developed a density functional tight-binding model for hydrocarbon excited state dynamics by referencing to high-level electronic structure theory and incorporating a many-body repulsive energy. We then validate our model against n-octane geometry optimizations, bond dissociation scans, and vibrational frequencies. Our model is approximately 10 3 times more efficient than hybrid time-dependent density functional theory calculations with comparable accuracy. In conclusion, our efforts enable longer timescale excited state simulations of photochemistry and scattering of incident radiation.

Chemical calculations↗

Layer-dependent topological phase in a two-dimensional quasicrystal and approximant

The electronic and topological properties of materials are derived from the interplay between crystalline symmetry and dimensionality. Simultaneously introducing “forbidden” symmetries via quasiperiodic ordering with low dimensionality into a material system promises the emergence of new physical phenomena. Here, we isolate a two-dimensional (2D) chalcogenide quasicrystal and approximant, and investigate their electronic and topological properties. The 2D layers of the materials with a composition close to Ta 1.6 Te, derived from a layered transition metal dichalcogenide, are isolated with standard exfoliation techniques, and investigated with electron diffraction and atomic resolution scanning transmission electron microscopy. Density functional theory calculations and symmetry analysis of the large unit cell crystalline approximant of the quasicrystal, Ta 21 Te 13 , reveal the presence of symmetry-protected nodal crossings in the quasicrystalline and approximant phases, whose presence is tunable by layer number. Our study provides a platform for the exploration of physics in quasicrystalline, low-dimensional materials and the interconnected nature of topology, dimensionality, and symmetry in electronic systems.

36 MATERIALS SCIENCE↗

Discovery versus precision in nuclear physics: A tale of three scales

At least three length scales are important in gaining a complete understanding of the physics of nuclei. These are the radius of the nucleus, the average inter-nucleon separation distance, and the size of the nucleon. The connections between the different scales are examined by using examples that demonstrate the direct connection between short-distance and high-momentum transfer physics and also that significant high-momentum content of wave functions is inevitable. The nuclear size is connected via the independent-pair approximation to the nucleon-nucleon separation distance, and this distance is connected via the concept of virtuality to the EMC effect. An explanation of the latter is presented in terms of light-front holographic wave functions of QCD. The net result is that the three scales are closely related, so that a narrow focus on any given specific range of scales may prevent an understanding of the fundamental origins of nuclear properties. It is also determined that, under certain suitable conditions, experiments are able to measure the momentum dependence of wave functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗

Self-consistent quasiparticle G W and hybrid functional calculations for Al/InAs/Al heterojunctions: Band offset and spin-orbit coupling effects

The electronic structure of surfaces and interfaces plays a key role in the properties of quantum devices. Here, we study the electronic structure of realistic Al/InAs/Al heterojunctions using a combination of density functional theory with hybrid functionals and state-of-the-art quasiparticle G W (QSGW) calculations. We find a good agreement between QSGW calculations and hybrid functional calculations, which themselves compare favorably well with angle-resolved photoemission spectroscopy experiments. Our paper confirms the need for well-controlled quality of the interfaces to obtain the needed properties of InAs/Al heterojunctions. A detailed analysis of the effects of spin-orbit coupling on the spin splitting of the electronic states shows a linear scaling in k space, related to the two-dimensional nature of some interface states. The good agreement by QSGW and hybrid functional calculations opens the door towards trustable use of an effective approximation to QSGW for studying very large heterojunctions. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Tensor squeezed limits and the Higuchi bound

We point out that tensor consistency relations—i.e. the behavior of primordial correlation functions in the limit a tensor mode has a small momentum—are more universal than scalar consistency relations. They hold in the presence of multiple scalar fields and as long as anisotropies are diluted exponentially fast. When de Sitter isometries are approximately respected during inflation this is guaranteed by the Higuchi bound, which forbids the existence of light particles with spin: de Sitter space can support scalar hair but no curly hair. We discuss two indirect ways to look for the violation of tensor consistency relations in observations, as a signature of models in which inflation is not a strong isotropic attractor, such as solid inflation: (a) graviton exchange contribution to the scalar four-point function; (b) quadrupolar anisotropy of the scalar power spectrum due to super-horizon tensor modes. This anisotropy has a well-defined statistics which can be distinguished from cases in which the background has a privileged direction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Composition Dependence of the Flory–Huggins Interaction Parameters of Block Copolymer Electrolytes and the Isotaksis Point

The thermodynamics of block copolymer/salt mixtures were quantified through the application of Leibler’s random phase approximation to disordered small-angle X-ray scattering profiles. The experimental system is comprised of polystyrene-block-poly(ethylene oxide) (SEO) mixed with lithium bis(trifluoromethanesulfonyl)imide salt (LiTFSI), SEO/LiTFSI. The Flory–Huggins interaction parameter determined from scattering experiments, χ SC , was found to be a function of block copolymer composition, chain length, and temperature for both salt-free and salty systems. In the absence of salt, χ 0,SC is a linear function of (N$f_{EO}$) -1 ; in the presence of salt, a linear approximation is used to describe the effect of salt on χ eff,SC for a given copolymer composition and chain length. The theory of Sanchez was used to determine χ eff from χ eff,SC to predict the boundary between order and disorder as a function of chain length, block copolymer composition, salt concentration, and temperature. At fixed temperature (100 °C), N crit , the chain length of SEO at the order–disorder transition in SEO/LiTFSI mixtures, was predicted as a function of the volume fraction of the salt-containing poly(ethylene oxide)-rich microphase, $f_{EO,salt,}$ and salt concentration. At $f_{EO,salt,}$ > 0.27, the addition of salt stabilizes the ordered phase; at $f_{EO,salt,}$ < 0.27, the addition of salt stabilizes the disordered phase. We propose a simple theoretical model to predict the block copolymer composition at which phase behavior is independent of salt concentration ($f_{EO,salt,}$ = 0.27). Here, we refer to this composition as the “isotaksis point”.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reliability of Open Public Electric Vehicle Direct Current Fast Chargers

The aim was to systematically evaluate the usability of all public electric vehicles (EV) direct current fast chargers (DCFC) in the San Francisco region. To achieve a rapid transition to EVs, a highly reliable and easy to use charging infrastructure is critical to building confidence among consumers. The functionality and usability of all 182 open, public DCFC charging stations with CCS connectors (combined charging system) in the 9 counties of the Bay Area were tested (655 electric vehicle service equipment (EVSE) ports). An EVSE was classified as functional if it charged an EV for 2 minutes. Overall, 73.3% of the 655 EVSEs were functional. The causes of the nonfunctioning EVSEs (23.5%) were blank or unresponsive screens or error messages; payment system failures; charge initiation failures; network failures; or broken connectors. In addition, the cable was too short to reach the EV inlet for 3.2% of the EVSEs. A random sampling of 10% of the EVSEs, approximately 8 days after the first evaluation, found no overall change in functionality. The level of functionality found with field testing conflicts with the 95–98% uptime reported by the EV service providers (EVSPs) who operate the EV charging stations. There is a need for precise and verifiable definitions of uptime, downtime, and excluded time, as applied to public EV chargers. In conclusion, the level of failure of the existing public EV DCFC charge infrastructure highlights the importance of improving the system design and maintenance to improve adoption of EVs.

33 ADVANCED PROPULSION SYSTEMS↗

An Agile Design-to-Simulation Workflow Using a New Conforming Moving Least Squares Method

This report summarizes the accomplishments and challenges of a two year LDRD effort focused on improving design-to-simulation agility. The central bottleneck in most solid mechanics simulations is the process of taking CAD geometry and creating a discretization of suitable quality, i.e., the "meshing" effort. This report revisits meshfree methods and documents some key advancements that allow their use on problems with complex geometries, low quality meshes, nearly incompressible materials or that involve fracture. The resulting capability was demonstrated to be an effective part of an agile simulation process by enabling rapid discretization techniques without increasing the time to obtain a solution of a given accuracy. The first enhancement addressed boundary-related challenges associated with meshfree methods. When using point clouds and Euclidean metrics to construct approximation spaces, the boundary information is lost, which results in low accuracy solutions for non-convex geometries and mate rial interfaces. This also complicates the application of essential boundary conditions. The solution involved the development of conforming window functions which use graph and boundary information to directly incorporate boundaries into the approximation space.

42 ENGINEERING↗

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semi-inclusive deep-inelastic scattering in Wandzura-Wilczek-type approximation

We present the complete cross-section for the production of unpolarized hadrons in semi-inclusive deep-inelastic scattering up to power-suppressed ${\cal O}(1/Q^2)$ terms in the Wandzura--Wilczek-type approximation which consists in systematically assuming that $\bar{q}gq$--terms are much smaller than $\bar{q}q$--correlators. We compute all twist-2 and twist-3 structure functions and the corresponding asymmetries, and discuss the applicability of the Wandzura--Wilczek-type approximations on the basis of available data. We make predictions which can be tested by data from Jefferson Lab, COMPASS, HERMES, and the future Electron-Ion Collider. The results of this paper can be readily used for phenomenology and for event generators, and will help to improve our understanding of the TMD theory beyond leading twist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cluster spin glass correlations and dynamics in Zn 0.5⁢ Mn 0.5⁢ Te

Here, we present a combined magnetometry, muon spin-relaxation (𝜇⁢SR), and neutron-scattering study of the insulating spin glass Zn 0.5 ⁢Mn 0.5 ⁢Te, for which magnetic Mn 2+ and nonmagnetic Zn 2+ ions are randomly distributed on a face-centered cubic lattice. The magnetometry and 𝜇⁢SR results confirm a spin freezing transition around 𝑇 𝑓 ≈ 23 K, with the spin-fluctuation rate decreasing gradually and somewhat inhomogeneously through the sample volume as the temperature decreases toward 𝑇 𝑓 . Characteristic spin-correlation times well above 𝑇 𝑓 are on the order of 10 −10 s, much slower than typically observed in canonical spin glasses but in line with expectations for a cluster spin glass. Using magnetic pair distribution function (mPDF) analysis and reverse Monte Carlo (RMC) modeling of the magnetic diffuse neutron-scattering data, we show that the spin-glass ground state consists of clusters of spins exhibiting short-range-ordered type-III antiferromagnetic correlations with a locally ordered moment of 3.1⁢(1)⁢𝜇 B between nearest-neighbor spins. The type-III correlations decay exponentially as a function of spin separation distance with a correlation length of approximately 5 Å. The diffuse magnetic scattering and corresponding mPDF show no significant changes across 𝑇 𝑓 , indicating that the dynamically fluctuating short-range spin correlations in the paramagnetic state retain the same basic type-III configuration that characterizes the spin-glass state; the only change apparent from the neutron-scattering data is a gradual reduction of the correlation length and locally ordered moment with increasing temperature. Taken together, these results paint a unique and detailed picture of the local magnetic structure and dynamics in Zn 0.5 ⁢Mn 0.5⁢ Te and provide strong evidence that this material is best described as a cluster spin glass. In addition, this work showcases a statistical method for extracting diffuse scattering signals from neutron powder diffraction data, which we developed to facilitate the mPDF and RMC analysis of the neutron data. This method has the potential to be broadly useful for neutron powder diffraction experiments on a variety of materials with short-range atomic or magnetic order.

magnetism↗

The Zeldovich approximation and wide-angle redshift-space distortions

The contribution of line-of-sight (LOS) peculiar velocities to the observed redshift of objects breaks the translational symmetry of the underlying theory, modifying the predicted two-point functions. These 'wide-angle effects' have mostly been studied using linear perturbation theory in the context of the multipoles of the correlation function and power spectrum. Here, we present the first calculation of wide-angle terms in the Zeldovich approximation, which is known to be more accurate than linear theory on scales probed by the next generation of galaxy surveys.We present the exact result for dark matter and perturbatively biased tracers as well as the small angle expansion of the configuration- and Fourier-space two-point functions and the connection to the multifrequency angular power spectrum. We compare different definitions of the LOS direction and discuss how to translate between them. Finally, we show that wide-angle terms can reach tens of percent of the total signal in a measurement at low redshift in some approximations, and that a generic feature of wide-angle effects is to slightly shift the Baryon Acoustic Oscillation scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Effective description of ajar systems with a 𝑈⁡(1) symmetry

We introduce the concept of “ajar systems” as an intermediate case between closed and open systems, where the time scale for charge exchange with the environment is parametrically larger than all other characteristic time scales. The Schwinger-Keldysh effective action for such finite-temperature systems exhibits a symmetry group G 1 × G 2 weakly broken to its diagonal subgroup G diag , which we systematically describe using spurion techniques. For G = U(1), we calculate leading-order corrections to correlation functions in both diffusive and spontaneously broken phases. Unlike systems with approximate symmetries where both real and imaginary parts of dispersion relations receive corrections, ajar systems modify only the imaginary part—preserving gapless modes while adding finite damping.

Effective field theory↗

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↗