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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

57Fe Mössbauer isomer shift of pure iron and iron oxides at high pressure—An experimental and theoretical study

The 57Fe isomer shift (IS) of pure iron has been measured up to 100 GPa using synchrotron Mössbauer spectroscopy in the time domain. Apart from the expected discontinuity due to the α→ε structural and spin transitions, the IS decreases monotonically with increasing pressure. The absolute shifts were reproduced without semi-empirical calibrations by periodic density functional calculations employing extensive localized basis sets with several common density functionals. However, the best numerical agreement is obtained with the B1WC hybrid functional. Extension of the calculations to 350 GPa, a pressure corresponding to the Earth’s inner core, predicted the IS range of 0.00 to -0.85 mm/s, covering the span from Fe(0) to Fe(VI) compounds measured at ambient pressure. The calculations also reproduced the pressure trend from polymorphs of prototypical iron oxide minerals, FeO and Fe2O3. Analysis of the electronic structure shows a strong donation of electrons from oxygen to iron at high pressure. The assignment of formal oxidation to the Fe atom becomes ambiguous under this condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1979

Solar p-mode oscillations as a tracer of radial differential rotation

Photoelectric observations of solar p-modes obtained with improved wavenumber and frequency resolution are presented. The observations are compared with model calculations of the p-modes, and the degree of spatial and temporal coherence of the observed wave pattern is investigated. It is found that the p-mode oscillations pervade the visible surface of the sun with a high degree of coherence in space and time, so that the whole complex pattern of standing waves with its nodes and antinodes can be regarded as a fixed pattern corotating with the solar surface layers. The p-modes are introduced as a tracer of solar rotational flow velocities. The equatorial differential rotation is estimated as a function of effective depth on the basis of the theoretical contribution functions for the p-modes recently derived by Ulrich et al. (1978). The results strongly indicate that the angular speed of rotation is not uniform even in the relatively shallow layer extending about 20,000 km below the photosphere.

Deubner, F.-L.↗
DOE OSTI2020

Crossing symmetry, transcendentality and the Regge behaviour of 1d CFTs

We develop the technology for Polyakov-Mellin (PM) bootstrap in one- dimensional conformal field theories (CFT 1 ). By adding appropriate contact terms, we bootstrap various effective field theories in AdS 2 and analytically compute the CFT data to one loop. The computation can be extended to higher orders in perturbation theory, if we ignore mixing, for any external dimension. We develop PM bootstrap for O(N) theories and derive the necessary contact terms for such theories (which also involves a new higher gradient contact term absent for N = 1). We perform cross-checks which include considering the diagonal limit of the 2d Ising model in terms of the 1d PM blocks. As an independent check of the validity of the results obtained with PM bootstrap, we propose a suitable basis of transcendental functions, which allows to fix the four-point correlators of identical scalar primaries completely, up to a finite number of ambiguities related to the number of contact terms in the PM basis. We perform this analysis both at tree level (with and without exchanges) and at one loop. We also derive expressions for the corresponding CFT data in terms of harmonic sums. Finally, we consider the Regge limit of one-dimensional correlators and derive a precise connection between the latter and the large-twist limit of CFT data. Exploiting this result, we study the crossing equation in the three OPE limits and derive some universal constraints for the large-twist limit of CFT data in Regge-bounded theories with a finite number of exchanges.

1/N expansion↗
DOE OSTI2021

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. Here, we select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.

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

Molecular electric moments calculated by using natural orbital functional theory

The molecular electric dipole, quadrupole, and octupole moments of a selected set of 21 spin-compensated molecules are determined employing the extended version of the Piris natural orbital functional 6 (PNOF6), using the triple-ζ Gaussian basis set with polarization functions developed by Sadlej, at the experimental geometries. The performance of the PNOF6 is established by carrying out a statistical analysis of the mean absolute errors with respect to the experiment. The calculated PNOF6 electric moments agree satisfactorily with the corresponding experimental data and are in good agreement with the values obtained by accurate ab initio methods, namely, the coupled-cluster single and doubles and multi-reference single and double excitation configuration interaction methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

HALE UAS Concept of Operations. Version 3.0

This document is a system level Concept of Operations (CONOPS) from the perspective of future High Altitude Long Endurance (HALE) Unmanned Aircraft Systems (UAS) service providers and National Airspace System (NAS) users. It describes current systems (existing UAS), describes HALE UAS functions and operations to be performed (via sample missions), and offers insight into the user s environment (i.e., the UAS as a system of systems). It is intended to be a source document for NAS UAS operational requirements, and provides a construct for government agencies to use in guiding their regulatory decisions, architecture requirements, and investment strategies. Although it does not describe the technical capabilities of a specific HALE UAS system (which do, and will vary widely), it is intended to aid in requirements capture and to be used as input to the functional requirements and analysis process. The document provides a basis for development of functional requirements and operational guidelines to achieve unrestricted access into the NAS. This document is an FY06 update to the FY05 Access 5 Project-approved Concept of Operations document previously published in the Public Domain on the Access 5 open website. This version is recommended to be approved for public release also. The updates are a reorganization of materials from the previous version with the addition of an updated set of operational requirements, inclusion of sample mission scenarios, and identification of roles and responsibilities of interfaces within flight phases.

Source record↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Giant gravitons, Harish-Chandra integrals, and BPS states in symplectic and orthogonal $\mathcal{N}$ = 4 SYM

We find generating functions for half BPS correlators in $\mathcal{N}$ = 4 SYM theories with gauge groups Sp(2N), SO(2N + 1), and SO(2N) by computing the norms of a class of BPS coherent states. These coherent states are built from operators involving Harish-Chandra integrals. Such operators have an interpretation as localized giant gravitons in the bulk of anti-de-Sitter space. This extends the analysis of [1] to Sp(2N), SO(2N + 1), and SO(2N) gauge theories. We show that we may use ordinary Schur functions as a basis for the sector of states with no cross-caps in these theories. This is consistent with the construction of these theories as orientifold projections of an SU(2N) theory. We make note of some relations between the symmetric functions that appear in the expansion of these coherent states and symplectic Schur functions. We also comment on some connections to Schubert calculus and Gromov-Witten invariants, which suggest that the Harish-Chandra integral may be extended to such problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Comparative analysis of proton and pion scattering on the isotopes {sup 6,8}He within Glauber theory

The differential cross sections for p{sup 6,8}He and p{sup 6,8}He scattering at the RIKEN andGSI energies (0.073 and 0.7 GeV per nucleon) were calculated on the basis of Glauber theory. Pion and proton interaction with a nucleus is determined by the multiple-scattering series (Glauber operator), so that various collision multiplicities and their contributions to the summed cross section can be taken into account. The use of the α-n-n wave function for {sup 6}He and the wave function on the basis of the large-scale shell model (LSSM) for {sup 8}He makes it possible to calculate analytically scattering matrix elements.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Spin orientations of the spin-half Ir{sup 4+} ions in Sr{sub 3}NiIrO{sub 6}, Sr{sub 2}IrO{sub 4}, and Na{sub 2}IrO{sub 3}: Density functional, perturbation theory, and Madelung potential analyses

The spins of the low-spin Ir{sup 4+} (S = 1/2, d{sup 5}) ions at the octahedral sites of the oxides Sr{sub 3}NiIrO{sub 6}, Sr{sub 2}IrO{sub 4}, and Na{sub 2}IrO{sub 3} exhibit preferred orientations with respect to their IrO{sub 6} octahedra. We evaluated the magnetic anisotropies of these S = 1/2 ions on the basis of density functional theory (DFT) calculations including spin-orbit coupling (SOC), and probed their origin by performing perturbation theory analyses with SOC as perturbation within the LS coupling scheme. The observed spin orientations of Sr{sub 3}NiIrO{sub 6} and Sr{sub 2}IrO{sub 4} are correctly predicted by DFT calculations, and are accounted for by the perturbation theory analysis. As for the spin orientation of Na{sub 2}IrO{sub 3}, both experimental studies and DFT calculations have not been unequivocal. Our analysis reveals that the Ir{sup 4+} spin orientation of Na{sub 2}IrO{sub 3} should have nonzero components along the c- and a-axis directions. The spin orientations determined by DFT calculations are sensitive to the accuracy of the crystal structures employed, which is explained by perturbation theory analyses when interactions between adjacent Ir{sup 4+} ions are taken into consideration. There are indications implying that the 5d electrons of Na{sub 2}IrO{sub 3} are less strongly localized compared with those of Sr{sub 3}NiIrO{sub 6} and Sr{sub 2}IrO{sub 4}. This implication was confirmed by showing that the Madelung potentials of the Ir{sup 4+} ions are less negative in Na{sub 2}IrO{sub 3} than in Sr{sub 3}NiIrO{sub 6} and Sr{sub 2}IrO{sub 4}. Most transition-metal S = 1/2 ions do have magnetic anisotropies because the SOC induces interactions among their crystal-field split d-states, and the associated mixing of the states modifies only the orbital parts of the states. This finding cannot be mimicked by a spin Hamiltonian because this model Hamiltonian lacks the orbital degree of freedom, thereby leading to the spin-half syndrome. The spin-orbital entanglement for the 5d spin-half ions Ir{sup 4+} is not as strong as has been assumed.

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

Expression, purification and primary crystallographic study of human androgen receptor in complex with DNA and coactivator motifs

The androgen receptor (AR) is a DNA-binding and hormone-activated transcription factor that plays critical roles in the development and progression of prostate cancer. The transcriptional function of AR is modulated by intermolecular interactions with DNA elements and coactivator proteins, as well as intramolecular interactions between AR domains; thus, the structural information from the full-length AR or a multi-domain fragment is essential for understanding the molecular basis of AR functions. Here we report the expression and purification of full-length AR protein and of a fragment containing its DNA-binding and ligand-binding domains connected by the hinge region in the presence of its natural ligand, dihydrotestosterone. Crystals of ligand-bound full-length AR and of the AR fragment in complex with DNA elements and coactivator motifs have been obtained and diffracted to low resolutions. These results help establish a foundation for pursuing further crystallographic studies of an AR/DNA complex.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1998

Do Bond Functions Help for the Calculation of Accurate Bond Energies?

The bond energies of 8 chemically bound diatomics are computed using several basis sets with and without bond functions (BF). The bond energies obtained using the aug-pVnZ+BF basis sets (with a correction for basis set superposition error, BSSE) tend to be slightly smaller that the results obtained using the aug-pV(n+I)Z basis sets, but slightly larger than the BSSE corrected aug-pV(n+I)Z results. The aug-cc-pVDZ+BF and aug-cc-pVTZ+BF basis sets yield reasonable estimates of bond energies, but, in most cases, these results cannot be considered highly accurate. Extrapolation of the results obtained with basis sets including bond functions appears to be inferior to the results obtained by extrapolation using atom-centered basis sets. Therefore bond functions do not appear to offer a path for obtaining highly accurate results for chemically bound systems at a lower computational cost than atom centered basis sets.

Bauschlicher, Charles W., Jr.↗
NASA NTRS2013

RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers

RUSHMAPS is a new onboard data reduction scheme that gives real-time access to key science parameters (e.g. moments) of a class of heliophysics science and/or solar system exploration investigation that includes plasma particle spectrometers (PPS), but requires moments reporting (density, bulk-velocity, temperature, pressure, etc.) of higher-level quality, and tolerates a lowpass (variable quality) spectral representation of the corresponding particle velocity distributions, such that telemetry use is minimized. The proposed methodology trades access to the full-resolution velocity distribution data, saving on telemetry, for real-time access to both the moments and an adjustable-quality (increasing quality increases volume) spectral representation of distribution functions. Traditional onboard data storage and downlink bandwidth constraints severely limit PPS system functionality and drive cost, which, as a consequence, drives a limited data collection and lower angular energy and time resolution. This prototypical system exploit, using high-performance processing technology at GSFC (Goddard Space Flight Center), uses a SpaceCube and/or Maestro-type platform for processing. These processing platforms are currently being used on the International Space Station as a technology demonstration, and work is currently ongoing in a new onboard computation system for the Earth Science missions, but they have never been implemented in heliospheric science or solar system exploration missions. Preliminary analysis confirms that the targeted processor platforms possess the processing resources required for realtime application of these algorithms to the spectrometer data. SpaceCube platforms demonstrate that the target architecture possesses the sort of compact, low-mass/power, radiation-tolerant characteristics needed for flight. These high-performing hybrid systems embed unprecedented amounts of onboard processing power in the CPU (central processing unit), FPGAs (field programmable gate arrays), and DSP (digital signal processing) elements. The fundamental computational algorithm de constructs 3D velocity distributions in terms of spherical harmonic spectral coefficients (which are analogous to a Fourier sine-cosine decomposition), but uses instead spherical harmonics Legendre polynomial orthogonal functions as a basis for the expansion, portraying each 2D angular distribution at every energy or, geometrically, spherical speed-shell swept by the particle spectrometer. Optionally, these spherical harmonic spectral coefficients may be telemetered to the ground. These will provide a smoothed description of the velocity distribution function whose quality will depend on the number of coefficients determined. Successfully implemented on the GSFC-developed processor, the capability to integrate the proposed methodology with both heritage and anticipated future plasma particle spectrometer designs is demonstrated (with sufficiently detailed design analysis to advance TRL) to show specific science relevancy with future HSD (Heliophysics Science Division) solar-interplanetary, planetary missions, sounding rockets and/or CubeSat missions.

Figueroa-Vinas, Adolfo↗
DOE OSTI2025

Optimizing stochastic algorithms for hadron correlation function computations in lattice QCD using a localized distillation basis

Distillation is a quark-smearing method for the construction of a broad class of hadron operators useful in lattice QCD computations and defined via a projection operator into a vector space of smooth gauge-covariant fields. A new orthonormal basis for this space is constructed which builds in locality. This basis is useful for the construction of stochastic methods to estimate the correlation functions computed in Monte Carlo calculations relevant for hadronic physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

PhenoProfiling: Mapping phenotypic outcomes to molecular determinants of biochemical activity

The aggregate genomes of the trillions of microorganisms within soil, animal hosts, and aquatic systems encode for an extensive functional capacity for myriad biochemical activities. C, N, P, and S metabolism, synthesis of signaling molecules and vitamins, mineralization, and other activities are essential to microbe, community, and plant physiology, and more broadly to climate, water, animal, flora, and human health. The current understanding of the molecular basis for the function of microbial communities stems primarily from comparative metagenomic and metatranscriptomic studies. These same tools are employed to ascertain the impacts to the community resulting from perturbations, such as climate change, emerging pollutants, fires, and seawater infiltration for environmental communities, and dietary changes, xenobiotic exposure, and various disease states for the human gut microbiome. Such studies can identify the potential for a specific function, but they cannot determine that a particular cell is functionally active, nor can they determine the molecular architecture required for function. In short, genes and transcripts alone fail to reveal the complex subcellular arrangement of proteins and molecules that elicit a given phenotype in a microbial cell.

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

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗