Search NASA⌕ Search

SEARCH · Search NASA

Results for “Polynomials”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Surface analysis insight note: Differentiation methods applicable to noisy data for determination of sp2‐ versus sp3‐hybridization of carbon allotropes and AES signal strengths

The derivatives of the spectra are commonly used for quantification in Auger Electron Spectroscopy (AES) spectra, while the derivative of the KLL C Auger line has proven to be valuable in obtaining a measure of the relative proportions of sp 2 ‐ and sp 3 ‐hybridization using the D‐parameter in both AES and X‐ray Photoelectron Spectroscopy (XPS). Differentiation of X‐ray Photoelectron Spectroscopy (XPS) and Auger Electron Spectroscopy (AES) spectra by numerical means is presented and illustrated for polymeric, such as PEEK and Nylon, as well as for graphitic materials including highly ordered pyrolytic graphite and graphene oxide. The most commonly available Savitzky–Golay method is explained mathematically and developed through the case of constructing a 5‐point quadratic polynomial convolution kernel suitable for differentiating spectra of adequate signal to noise. The concept of differentiation of spectra where signal to noise is less than adequate is also developed. Two alternative strategies to Savitzky–Golay differentiation are presented, which fit curves to data that allow derivatives to be obtained where Savitzky–Golay would otherwise fail. These alternative methods involve constructing a parametric curve that fits data over the entire energy interval of interest. Derivatives of spectra are then obtained by differentiating these parametric curves directly. A comparison of results for different materials for which specific sp 2 ‐ vs sp 3 ‐hybridized carbon proportions are of interest is used to emphasize the importance of characterizing methods used to differentiate spectra and understanding the characteristics of instrumentation used to measure spectra. The case for using Principal Component Analysis noise reduction with C KLL spectra is made for spectra collected from a heterogeneous graphene oxide sample.

Fairley, Neal↗
DOE OSTI2020

Calibration, measurement, and characterization of soil moisture dynamics in a central Amazonian tropical forest

Soil moisture plays a key role in hydrological, biogeochemical, and energy budgets of terrestrial ecosystems. Accurate soil moisture measurements in remote ecosystems such as the Amazon are difficult and limited because of logistical constraints. Time domain reflectometry (TDR) sensors are widely used to monitor soil moisture and require calibration to convert the TDR's dielectric permittivity measurement (K a ) to volumetric water content (θ v ). In this study, our objectives were to develop a field-based calibration of TDR sensors in an old-growth upland forest in the central Amazon, to evaluate the performance of the calibration, and then to apply the calibration to determine the dynamics of soil moisture content within a 14.2 -m-deep vertical soil profile. Depth-specific TDR calibration using local soils in a controlled laboratory setting yielded a novel K a –θ v third-degree polynomial calibration. The sensors were later installed to their specific calibration depth in a 14.2-m pit. The widely used K a –θ v relationship (Topp model) underestimated the site-specific θ v by 22–42%, indicating significant error in the model when applied to these well-structured, clay-rich tropical forest soils. The calibrated wet- and dry-season θ v data showed a variety of depth and temporal variations highlighting the importance of soil textural differentiation, root uptake depths, as well as event to seasonal precipitation effects. Data such as these are greatly needed for improving our understanding of ecohydrological processes within tropical forests and for improving models of these systems in the face of changing environmental conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Veneziano variations: how unique are string amplitudes?

String theory offers an elegant and concrete realization of how to consistently couple states of arbitrarily high spin. But how unique is this construction? In this paper we derive a novel, multi-parameter family of four-point scattering amplitudes exhibiting i) polynomially bounded high-energy behavior and ii) exchange of an infinite tower of high-spin modes, albeit with a finite number of states at each resonance. These amplitudes take an infinite-product form and, depending on parameters, exhibit mass spectra that are either unbounded or bounded, thus corresponding to generalizations of the Veneziano and Coon amplitudes, respectively. For the bounded case, masses converge to an accumulation point, a peculiar feature seen in the Coon amplitude but more recently understood to arise naturally in string theory [1]. Importantly, our amplitudes contain free parameters allowing for the customization of the slope and offset of the spin-dependence in the Regge trajectory. We compute all partial waves for this multi-parameter class of amplitudes and identify unitary regions of parameter space. For the unbounded case, we apply similar methods to derive new deformations of the Veneziano and Virasoro-Shapiro amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Model selection and signal extraction using Gaussian Process regression

We present a novel computational approach for extracting localized signals from smooth background distributions. We focus on datasets that can be naturally presented as binned integer counts, demonstrating our procedure on the CERN open dataset with the Higgs boson signature, from the ATLAS collaboration at the Large Hadron Collider. Our approach is based on Gaussian Process (GP) regression — a powerful and flexible machine learning technique which has allowed us to model the background without specifying its functional form explicitly and separately measure the background and signal contributions in a robust and reproducible manner. Unlike functional fits, our GP-regression-based approach does not need to be constantly updated as more data becomes available. We discuss how to select the GP kernel type, considering trade-offs between kernel complexity and its ability to capture the features of the background distribution. We show that our GP framework can be used to detect the Higgs boson resonance in the data with more statistical significance than a polynomial fit specifically tailored to the dataset. Finally, we use Markov Chain Monte Carlo (MCMC) sampling to confirm the statistical significance of the extracted Higgs signature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Exact quantization and analytic continuation

In this paper we give a streamlined derivation of the exact quantization condition (EQC) on the quantum periods of the Schrödinger problem in one dimension with a general polynomial potential, based on Wronskian relations. We further generalize the EQC to potentials with a regular singularity, describing spherical symmetric quantum mechanical systems in a given angular momentum sector. We show that the thermodynamic Bethe ansatz (TBA) equations that govern the quantum periods undergo nontrivial monodromies as the angular momentum is analytically continued between integer values in the complex plane. The TBA equations together with the EQC are checked numerically against Hamiltonian truncation at real angular momenta and couplings, and are used to explore the analytic continuation of the spectrum on the complex angular momentum plane in examples.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Non-singlet quark helicity $\mathrm{PDF}$s of the nucleon from pseudo-distributions

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional z 2 contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional z 2 contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale μ 2 . The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned z 2 contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Multi-matrix correlators and localization

Abstract We study generating functions of$$ \frac{1}{4} $$ 1 4 -BPS states in$$ \mathcal{N} $$ N = 4 super Yang-Mills at finiteNby attempting to generalize the Harish-Chandra-Itzykson-Zuber integral to multiple commuting matrices. This allows us to compute the overlaps of two or more generating functions; such calculations arise in the computation of two-point correlators in the free-field limit. We discuss the four-matrix HCIZ integral in the U(2) context and lay out a prescription for finding a more general formula forN> 2. We then discuss its connections with the restricted Schur polynomial operator basis. Our results generalize readily to arbitrary numbers of matrices, opening up the opportunity to study more generic BPS operators.

Physics↗
DOE OSTI2023

Orbit averaging coherent states: holographic three-point functions of AdS giant gravitons

We study correlation functions of two AdS giant gravitons in AdS 5 × S 5 and a BPS supergravity mode using holography. In the gauge theory these are described by BPS correlators of Schur polynomials of fully-symmetric representations and a single trace operator. We find full agreement between the semiclassical gravity and gauge theory computations at large N, for both diagonal and off-diagonal structure constants. Our analysis in $\mathcal{N}$ = 4 SYM provides a simpler derivation to the results in the literature, and it can be readily generalized to operators describing bound states of AdS giant gravitons as well as bubbling geometries.

1/N Expansion↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Vector boson pair production at one loop: analytic results for the process $q\bar{q}$$ℓ\bar{ℓ}$$ℓ'\bar{ℓ}'g$

We present compact analytic results for the one-loop amplitude for the process 0 → $q\bar{q}$$ℓ\bar{ℓ}$$ℓ'\bar{ℓ}'g$, relevant for both the production of a pair of Z and W-bosons in association with a jet. We focus on the gauge-invariant contribution mediated by a loop of quarks. We explicitly include all effects of the loop-quark mass m, appropriate for the production of a pair of Z-bosons. In the limit m → 0, our results are also applicable to the production of W-boson pairs, mediated by a loop of massless quarks. Implemented in a numerical code, the results are fast. The calculation uses novel advancements in spinor-helicity simplification techniques, for the first time applied beyond five-point massless kinematics. We make use of primary decompositions from algebraic-geometry, which now involve non-radical ideals, and p-adic numbers from number theory. We show how to infer whether numerator polynomials belong to symbolic powers of non-radical ideals through numerical evaluations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Quantum error correction from complexity in Brownian SYK

We study the robustness of quantum error correction in a one-parameter ensemble of codes generated by the Brownian SYK model, where the parameter quantifies the encoding complexity. The robustness of error correction by a quantum code is upper bounded by the “mutual purity” of a certain entangled state between the code subspace and environment in the isometric extension of the error channel, where the mutual purity of a density matrix ρAB is the difference $\mathcal{F}$ p ($A : B$) ≡ $\mathrm{T}$r $p^{2}_{AB}$ - $\mathrm{T}$r $p^{2}_{A}$ $\mathrm{T}$r $p^{2}_{B}$. We show that when the encoding complexity is small, the mutual purity is O(1) for the erasure of a small number of qubits (i.e., the encoding is fragile). However, this quantity decays exponentially, becoming O(1/N) for O(log N) encoding complexity. Further, at polynomial encoding complexity, the mutual purity saturates to a plateau of O(e -N ). We also find a hierarchy of complexity scales associated to a tower of subleading contributions to the mutual purity that quantitatively, but not qualitatively, adjust our error correction bound as encoding complexity increases. In the AdS/CFT context, our results suggest that any portion of the entanglement wedge of a general boundary subregion A with sufficiently high encoding complexity is robustly protected against low-rank errors acting on A with no prior access to the encoding map. From the bulk point of view, we expect such bulk degrees of freedom to be causally inaccessible from the region A despite being encoded in it.

1/N expansion↗
DOE OSTI2020

q-nonabelianization for line defects

We consider the q-nonabelianization map, which maps links L in a 3-manifold M to combinations of links $\tilde{L}$ in a branched N -fold cover $\tilde{M}$. In quantum field theory terms, q-nonabelianization is the UV-IR map relating two different sorts of defect: in the UV we have the six-dimensional (2, 0) superconformal field theory of type gl(N ) on M × R 2,1 , and we consider surface defects placed on L × {x 4 = x 5 = 0}; in the IR we have the (2, 0) theory of type gl (1) on $\tilde{M}$ × R 2,1 , and put the defects on $\tilde{L}$ × {x 4 = x 5 = 0}. In the case M = R 3 , q-nonabelianization computes the Jones polynomial of a link, or its analogue associated to the group U(N ). In the case M = C × R, when the projection of L to C is a simple non-contractible loop, q-nonabelianization computes the protected spin character for framed BPS states in 4d N = 2 theories of class S. In the case N = 2 and M = C × R, we give a concrete construction of the q-nonabelianization map. The construction uses the data of the WKB foliations associated to a holomorphic covering $\tilde{C}$ →C.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The quantum $\mathrm{UV}$-$\mathrm{IR}$ map for line defects in $$ \mathfrak{gl} $$(3)-type class $\mathcal{S}$ theories

We consider the quantum UV-IR map for line defects in class S theories of $\mathfrak{gl}$(3)-type. This map computes the protected spin character which counts framed BPS states with spin for the bulk-defect system. We give a geometric method of computing this map motivated by the physics of five-dimensional $\mathcal{N}$ = 2 supersymmetric Yang-Mills theory, and compute it explicitly in various examples. As a spin-off we propose a new way of computing a certain specialization of the HOMFLY polynomial for links in $\mathbb{R}^3$, as a sum over BPS webs attached to the link.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Interacting CFTs for all couplings: thermal versus entanglement entropy at large N

In this paper, I calculate the large N limit of marginal O(N) models with non-polynomial potentials in arbitrary odd dimensions d. This results in a new class of interacting pure conformal field theories (CFTs) in d = 3 + 4n for any n ϵ $\mathbb{Z}$ + . Similarly, in d = 3 + 4n I calculate the thermal entropy for all couplings on R 2+4n × S 1 for n = 0, 1, 2, 3. In 2+1 dimensions I find the strong-to-weak coupling ratio of the thermal entropy to be 4/5, matching recent results, and further extend this analysis to higher odd dimensions. Next, I calculated the vacuum entanglement entropy ${s}_{\textrm{EE}}^d$ on S d–2 for all couplings in arbitrary odd d in the large N limit. I find the vacuum entanglement entropy on S d–2 to be not only solvable but also constant for all couplings λ. Thus, in the large N limit, the vacuum entanglement entropy on S d–2 for odd d is constant for all λ, in contrast to the thermal entropy which is shown to also be monotonically decreasing with λ in d = 3 + 4n.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Giant Wilson loops and AdS 2 /dCFT 1

The 1/2-BPS Wilson loop in N = 4 supersymmetric Yang-Mills theory is an important and well-studied example of conformal defect. In particular, much work has been done for the correlation functions of operator insertions on the Wilson loop in the fundamental representation. In this paper, we extend such analyses to Wilson loops in the large-rank symmetric and antisymmetric representations, which correspond to probe D3 and D5 branes with AdS 2 × S 2 and AdS 2 × S 4 worldvolume geometries, ending at the AdS 5 boundary along a one-dimensional contour. We first compute the correlation functions of protected scalar insertions from supersymmetric localization, and obtain a representation in terms of multiple integrals that are similar to the eigenvalue integrals of the random matrix, but with important differences. Using ideas from the Fermi Gas formalism and the Clustering method, we evaluate their large N limit exactly as a function of the ’t Hooft coupling. The results are given by simple integrals of polynomials that resemble the Q-functions of the Quantum Spectral Curve, with integration measures depending on the number of insertions. Next, we study the correlation functions of fluctuations on the probe D3 and D5 branes in AdS. We compute a selection of three- and four-point functions from perturbation theory on the D-branes, and show that they agree with the results of localization when restricted to supersymmetric kinematics. We also explain how the difference of the internal geometries of the D3 and D5 branes manifests itself in the localization computation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗