Search NASA⌕ Search

SEARCH · Search NASA

Results for “Recursion”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Universal lowest-twist in CFTs from holography

We probe the conformal block structure of a scalar four-point function in d ≥ 2 conformal field theories by including higher-order derivative terms in a bulk gravitational action. We consider a heavy-light four-point function as the boundary correlator at large central charge. Such a four-point function can be computed, on the gravity side, as a two-point function of the light operator in a black hole geometry created by the heavy operator. We consider analytically solving the corresponding scalar field equation in a near-boundary expansion and find that the multi-stress tensor conformal blocks are insensitive to the horizon boundary condition. The main result of this paper is that the lowest-twist operator product expansion (OPE) coefficients of the multi-stress tensor conformal blocks are universal: they are fixed by the dimension of the light operators and the ratio between the dimension of the heavy operator and the central charge CT. Neither supersymmetry nor unitary is assumed. Higher-twist coefficients, on the other hand, generally are not protected. A recursion relation allows us to efficiently compute universal lowest-twist coefficients. The universality result hints at the potential existence of a higher-dimensional Virasoro-like symmetry near the lightcone. While we largely focus on the planar black hole limit in this paper, we include some preliminary analysis of the spherical black hole case in an appendix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Measurement of the groomed jet mass in PbPb and pp collisions at $ \sqrt{s_{\mathrm{NN}}}=5.02 $ TeV

Here, a measurement of the groomed jet mass in PbPb and pp collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV with the CMS detector at the LHC is presented. Jet grooming is a recursive procedure which sequentially removes soft constituents of a jet until a pair of hard subjets is found. The resulting groomed jets can be used to study modifications to the parton shower evolution in the presence of the hot and dense medium created in heavy ion collisions. Predictions of groomed jet properties from the PYTHIA and HERWIG++ event generators agree with the measurements in pp collisions. When comparing the results from the most central PbPb collisions to pp data, a hint of an increase of jets with large jet mass is observed, which could originate from additional medium-induced radiation at a large angle from the jet axis. However, no modification of the groomed mass of the core of the jet is observed for all PbPb centrality classes. The PbPb results are also compared to predictions from the JEWEL and Q-PYTHIA event generators, which predict a large modification of the groomed mass not observed in the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Feynman rules for scalar conformal blocks

We complete the proof of “Feynman rules” for constructing M-point conformal blocks with external and internal scalars in any topology for arbitrary M in any spacetime dimension by combining the rules for the blocks (based on their Witten diagram interpretation) with the rules for the construction of conformal cross ratios (based on the OPE and “flow diagrams”). The full set of Feynman rules leads to blocks as power series of the hypergeometric type in the conformal cross ratios. We then provide a proof by recursion of the Feynman rules which relies heavily on the first Barnes lemma and the decomposition of the topology of interest in comb structures. Finally, we provide a nine-point example to illustrate the rules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Top tree amplitudes for higher order calculations

We present compact analytic results for tree-level amplitudes containing a $t\bar{t}$ pair accompanied by up to four massless partons $t\bar{t}gg$, $t\bar{t}ggg$, $t\bar{t}gggg$, $t\bar{t}q\bar{q}$, $t\bar{t}q\bar{q}g$, $t\bar{t}q\bar{q}gg$, and $t\bar{t}q\bar{q}q'\bar{q'}$. The results, obtained using BCFW on-shell recursion, are based both on previous published results and on the new calculations performed in this paper. These amplitudes are sufficient to calculate the production of a $t\bar{t}$ pair and zero, one, or two light parton jets, with the option to include the tree-level decays t → bνe + and $\bar{t} → \bar{b}e^-\bar{v}$ efficiently. Our results are part of the NNLO corrections to $t\bar{t}$ production including the decay correlations for on-shell top quarks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Integrating the full four-loop negative geometries and all-loop ladder-type negative geometries in ABJM theory

The decomposition of the four-point ABJM amplituhedron into negative geometries produces compact integrands of logarithmic of amplitudes such that the infrared divergence only comes from the last loop integration, from which we can compute the cusp anomalous dimension of the ABJM theory. In this note, we integrate L – 1 loop momenta of the L-loop negative geometries for all four-loop negative geometries and a special class of all-loop ladder-type negative geometries by a method based on Mellin transformation, and from these finite quantities we extract the corresponding contribution to the cusp anomalous dimension. We find that the infrared divergence of a box-type negative geometry at L = 4 is weaker than other negative geometries, then only tree-type negative geometries contribute to the cusp anomalous dimension at L = 4. For the all-loop ladder-type negative geometries, we prove and conjecture some recursive structures as integral equations in Mellin space and find that they cannot contribute zeta values like ζ 3 , ζ 5 to the cusp anomalous dimension.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Probing universalities in d > 2 CFTs: from black holes to shockwaves

Gravitational shockwaves are insensitive to higher-curvature corrections in the action. Recent work found that the OPE coefficients of lowest-twist multi-stress-tensor operators, computed holographically in a planar black hole background, are insensitive as well. In this paper, we analyze the relation between these two limits. We explicitly evaluate the two-point function on a shockwave background to all orders in a large central charge expansion. In the geodesic limit, we find that the ANEC exponentiates in the multi-stress-tensor sector. To compare with the black hole limit, we obtain a recursion relation for the lowest-twist products of two stress tensors in a spherical black hole background, letting us efficiently compute their OPE coefficients and prove their insensitivity to higher curvature terms. After resumming the lowest-twist stress-tensors and analytically continuing their contributions to the Regge limit, we find a perfect agreement with the shockwave computation. We also discuss the role of double-trace operators, global degenerate states, and multi-stress-tensor conformal blocks. These holographic results suggest the existence of a larger universal structure in higher-dimensional CFTs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

blocks_3d: software for general 3d conformal blocks

We introduce the software blocks_3d for computing four-point conformal blocks of operators with arbitrary Lorentz representations in 3d CFTs. It uses Zamolodchikov-like recursion relations to numerically compute derivatives of blocks around a crossing-symmetric configuration. It is implemented as a heavily optimized, multi-threaded, C++ application. We give performance benchmarks for correlators containing scalars, fermions, and stress tensors. As an example application, we recompute bootstrap bounds on four-point functions of fermions and study whether a previously observed sharp jump can be explained using the “fake primary” effect. We conclude that the fake primary effect cannot fully explain the jump and the possible existence of a “dead-end” CFT near the jump merits further study.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Non-planar BCFW Grassmannian geometries

In this paper, we study non-adjacent BCFW recursion relations and their connection to positive geometry. For an adjacent BCFW shift, the n-point N k MHV tree-level amplitude in N = 4 SYM theory is expressed as a sum over planar on-shell diagrams, corresponding to canonical “dlog” forms on the cells in the positive Grassmannian G+(k, n). Non-adjacent BCFW shifts naturally lead to an expansion of the amplitude in terms of a different set of objects, which do not manifest the cyclic ordering and the hidden Yangian symmetry of the amplitude. We show that these terms can be interpreted as dlog forms on the non-planar Grassmannian geometries, generalizing the cells of the positive Grassmannian G+(k, n) to a larger class of objects which live in G(k, n). We focus mainly on the case of NMHV amplitudes and discuss in detail the Grassmannian geometries. We also propose an alternative way to calculate the associated on-shell functions and dlog forms using an intriguing connection between Grassmannian configurations and the geometry in the kinematical space.

1/N expansion↗
DOE OSTI2022

134 billion intersecting brane models

The landscape of string vacua is very large, but generally expected to be finite in size. Enumerating the number and properties of the vacua is an important task for both the landscape and the swampland, in part to gain a deeper understanding of what is possible and “generic”. We obtain an exact counting of distinct intersecting brane vacua of type IIA string theory on the T 6 /Z 2 × Z 2 orientifold. Care is taken to only count gauge-inequivalent brane configurations. Leveraging the recursive nature by which branes may be added together one-by-one, we use dynamic programming to efficiently count the number of solutions of the tadpole, K-theory and supersymmetry consistency conditions. The distributions of 4D gauge group rank and complex structure moduli for the entire ensemble of intersecting brane vacua are presented. The methods we developed here may be useful in obtaining sharp upper and lower bounds on other corners of the landscape.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Goldstone bosons on celestial sphere and conformal soft theorems

In this paper, we study celestial amplitudes of Goldstone bosons and conformal soft theorems. Motivated by the success of soft bootstrap in momentum space and the important role of the soft limit behavior of tree-level amplitudes, our goal is to extend some of the methods to the celestial sphere. The crucial ingredient of the calculation is the Mellin transformation, which transforms four-dimensional scattering amplitudes to correlation functions of primary operators in the celestial CFT. The soft behavior of the amplitude is then translated to the singularities of the correlator. Only for amplitudes in “UV completed theories” (with sufficiently good high energy behavior) the Mellin integration can be properly performed. In all other cases, the celestial amplitude is only defined in a distributional sense with delta functions. We provide many examples of celestial amplitudes in UV-completed models, including linear sigma models and Z-theory, which is a certain completion of the SU(N) non-linear sigma model. We also comment on the BCFW-like and soft recursion relations for celestial amplitudes and the extension of soft bootstrap ideas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Off-shell Partition Functions in 3d Gravity

We explore three-dimensional gravity with negative cosmological constant via canonical quantization. We focus on chiral gravity which is related to a single copy of PSL(2,R) Chern-Simons theory and is simpler to treat in canonical quantization. Its phase space for an initial value surface Σ is given by the appropriate moduli space of Riemann surfaces. We use geometric quantization to compute partition functions of chiral gravity on three-manifolds of the form Σ×S 1 , where Σ can have asymptotic boundaries. Most of these topologies do not admit a classical solution and are thus not amenable to a direct semiclassical path integral computation. We use an index theorem that expresses the partition function as an integral of characteristic classes over phase space. In the presence of n asymptotic boundaries, we use techniques from equivariant cohomology to localize the integral to a finite-dimensional integral over $\overline{M}$ g,n , which we evaluate in low genus cases. Higher genus partition functions quickly become complicated since they depend in an oscillatory way on Newton's constant. There is a precise sense in which one can isolate the non-oscillatory part which we call the fake partition function. We establish that there is a topological recursion that computes the fake partition functions for arbitrary Riemann surfaces Σ. As a result, there is a scaling limit in which the model reduces to JT gravity and our methods give a novel way to compute JT partition functions via equivariant localization.

Classical and Quantum Gravity↗
DOE OSTI2021

Compact representations of structured BFGS matrices

For general large-scale optimization problems compact representations exist in which recursive quasi-Newton update formulas are represented as compact matrix factorizations. For problems in which the objective function contains additional structure, recent structured quasi-Newton methods exploit available second-derivative information and approximate unavailable second derivatives. Here, this article develops the compact representations of two structured Broyden-Fletcher-Goldfarb-Shanno update formulas. The compact representations enable efficient limited memory and initialization strategies. Two limited memory line search algorithms are described for which extensive numerical results demonstrate the efficacy of the algorithms, including comparisons to IPOPT on large machine learning problems, and to L-BFGS on a real world large scale ptychographic imaging application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Combination of the Laguerre Transform with the Boundary-element Method for the Solution of Integral Equations with Retarded Kernel

We apply the Laguerre transform with respect to time to a time-dependent boundary-value integral equation encountered in the solution of three-dimensional Dirichlet initial-boundary-value problems for the homogeneous wave equation with homogeneous initial conditions by using the retarded potential of single layer. The obtained system of boundary integral equations is reduced to a sequence of Fredholm integral equations of the first kind that differ solely by the recursively dependent right-hand sides. To find their numerical solution, we use the boundary-element method. We establish an asymptotic estimate of the error of numerical solution and present the results of numerical simulations aimed at finding the solutions of retarded-potential integral equations for model examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Reflection and transmission of electromagnetic waves by a multilayered solar cell containing organic materials

Light trapping or absorption is a key concept for design and fabricating solar cells. The propagation of electromagnetic waves through a multilayered solar cell containing organic materials has theoretically investigated. In particular, the reflected, transmitted, and absorbed powers are found and defined. The main parameters of the structure and the required equations for TE polarization are presented in detail. The boundary conditions of the electric and magnetic fields of the incident waves are imposed at each interface and Snell’s law is applied. Expressions for the overall reflection and transmission coefficients are provided by a recursive method. The reflected, transmitted, and absorbed powers are determined using these coefficients. In the numerical results, the mentioned powers are computed and illustrated as a function of the angle of incidence, the wavelength, and the thickness of the active layer when the extinction coefficient of the active layer changes.

14 SOLAR ENERGY↗
DOE OSTI2021

Modeling household online shopping demand in the U.S.: a machine learning approach and comparative investigation between 2009 and 2017

Despite the rapid growth of online shopping and research interest in the relationship between online and in-store shopping, national-level modeling and investigation of the demand for online shopping with a prediction focus remain limited in the literature. Here, this paper differs from prior work and leverages two recent releases of the U.S. National Household Travel Survey (NHTS) data for 2009 and 2017 to develop machine learning (ML) models, specifically gradient boosting machine (GBM), for predicting household-level online shopping purchases. The NHTS data allow for not only conducting nationwide investigation but also at the level of households, which is more appropriate than at the individual level given the connected consumption and shopping needs of members in a household. We follow a systematic procedure for model development including employing Recursive Feature Elimination algorithm to select input variables (features) in order to reduce the risk of model overfitting and increase model explainability. Among several ML models, GBM is found to yield the best prediction accuracy. Extensive post-modeling investigation is conducted in a comparative manner between 2009 and 2017, including quantifying the importance of each input variable in predicting online shopping demand, and characterizing value-dependent relationships between demand and the input variables. In doing so, two latest advances in machine learning techniques, namely Shapley value-based feature importance and Accumulated Local Effects plots, are adopted to overcome inherent drawbacks of the popular techniques in current ML modeling. The modeling and investigation are performed at the national level, with a number of findings obtained. The models developed and insights gained can be used for online shopping-related freight demand generation and may also be considered for evaluating the potential impact of relevant policies on online shopping demand.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion

We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

ZERNIPAX: A fast and accurate Zernike polynomial calculator in Python

Zernike polynomials serve as an orthogonal basis on the unit disc, and have proven to be effective in optics simulations, astrophysics, and more recently in plasma simulations. Unlike Bessel functions, Zernike polynomials are inherently finite and smooth at the disc center (r=0), ensuring continuous differentiability along the axis. This property makes them particularly suitable for simulations, requiring no additional handling at the origin. We developed ZERNIPAX, an open-source Python package capable of utilizing CPU/GPUs, leveraging Google's JAX package and available on GitHub as well as the Python software repository PyPI. Furthermore, our implementation of the recursion relation between Jacobi polynomials significantly improves computation time compared to alternative methods by use of parallel computing while still performing more accurately for high-mode numbers.

Astrophysics↗
DOE OSTI2019

First order formulation of the Yang–Mills theory in a background field

The background gauge renormalization of the first order formulation of the Yang–Mills theory is studied by using the BRST identities. Together with the background symmetry, these identities allow for an iterative proof of renormalizability to all orders in perturbation theory. However, due to the fact that certain improper diagrams which violate the BRST symmetry should be removed, the renormalizability must be deduced indirectly. The recursive method involves rescalings and mixings of the fields, which lead to a renormalized effective action for the background field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗