Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Filtered Rayleigh-Ritz is all you need

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Field Theory of the Fermi Function

The Fermi function F ( Z , E ) accounts for QED corrections to beta decays that are enhanced at either small electron velocity β or large nuclear charge Z . For precision applications, the Fermi function must be combined with other radiative corrections and with scale- and scheme-dependent hadronic matrix elements. We formulate the Fermi function as a field theory object and present a new factorization formula for QED radiative corrections to beta decays. We provide new results for the anomalous dimension of the corresponding effective operator complete through three loops, and resum perturbative logarithms and π enhancements with renormalization-group methods. Our results are important for tests of fundamental physics with precision beta decay and related processes. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2018

Scattering and Radiative Properties of Morphologically Complex Carbonaceous Aerosols: A Systematic Modeling Study

This paper provides a thorough modeling-based overview of the scattering and radiative properties of a wide variety of morphologically complex carbonaceous aerosols. Using the numerically-exact superposition T-matrix method, we examine the absorption enhancement, absorption Angstroem exponent (AAE), backscattering linear depolarization ratio (LDR), and scattering matrix elements of black-carbon aerosols with 11 different model morphologies ranging from bare soot to completely embedded soot-sulfate and soot-brown carbon mixtures. Our size-averaged results show that fluffy soot particles absorb more light than compact bare-soot clusters. For the same amount of absorbing material, the absorption cross section of internally mixed soot can be more than twice that of bare soot. Absorption increases as soot accumulates more coating material and can become saturated. The absorption enhancement is affected by particle size, morphology, wavelength, and the amount of coating. We refute the conventional belief that all carbonaceous aerosols have AAEs close to 1.0. Although LDRs caused by bare soot and certain carbonaceous particles are rather weak, LDRs generated by other soot-containing aerosols can reproduce strong depolarization measured by Burton et al. for aged smoke. We demonstrate that multi-wavelength LDR measurements can be used to identify the presence of morphologically complex carbonaceous particles, although additional observations can be needed for full characterization. Our results show that optical constants of the host/coating material can significantly influence the scattering and absorption properties of soot-containing aerosols to the extent of changing the sign of linear polarization. We conclude that for an accurate estimate of black-carbon radiative forcing, one must take into account the complex morphologies of carbonaceous aerosols in remote sensing studies as well as in atmospheric radiation computations.

carbonaceous aerosols; scattering matrix; polariza↗
DOE OSTI2024

Measurement of single top-quark production in association with a W boson in pp collisions at $\sqrt{s}$ =13 TeV with the ATLAS detector

The inclusive cross section for the production of a single top quark in association with a Formula Presented boson is measured using Formula Presented of proton-proton collision data collected with the ATLAS detector at $\sqrt{s}$ =13 TeV. Events containing two charged leptons and at least one jet identified as originating from a b-quark are selected. A multivariate discriminant is constructed to separate the tW signal from the $t\bar{t}$ background. The cross section is extracted using a profile likelihood fit to the signal and control regions and it is measured to be σ tW = 75$^{+15}_{–14}$, in good agreement with the Standard Model prediction. The measured cross section is used to extract a value for the left-handed form factor at the Wtb vertex times the Cabibbo-Kobayashi-Maskawa matrix element |$f$ LV $V$ tb | of 0.97±0.10.

Cabibbo–Kobayashi–Maskawa matrix↗
DOE OSTI2023

Measurement of differential distributions of B → D*$\ell\bar{v}_{\ell}$ and implications on | V cb |

We present a measurement of the differential shapes of exclusive B → D*$\ell\bar{v}_{\ell}$ (B = B - , $\bar{B}$ 0 and $\ell$ = e, μ) decays with hadronic tag-side reconstruction for the full 711 fb -1 Belle dataset. We extract the Caprini-Lellouch-Neubert (CLN) and Boyd-Grinstein-Lebed (BGL) form factor parameters and use an external input for the absolute branching fractions to determine the Cabibbo-Kobayashi-Maskawa matrix element and find |V cb | CLN = (40.2 ± 0.9) x 10 -3 and |V cb | BGL = (40.7 ± 1.0) x 10 -3 with the zero-recoil lattice QCD point $\mathcal{F}$(1) = 0.906 ± 0.013. We also perform a study of the impact of beyond zero-recoil lattice QCD calculations on the |V cb | determinations. Additionally, we present the lepton-flavor universality ratio R eμ = $\mathcal{B}$(B → D*e$\bar{v}$ e )/$\mathcal{B}$(B → D*μ$\bar{v}$ μ ) = 0.993 ± 0.023 ± 0.023, the electron and muon forwardbackward asymmetry and their difference ΔA FB = 0.028 ± 0.028 ± 0.008, and the electron and muon D* longitudinal polarization fraction and their difference ΔF$^{D*}_L$ = 0.030 ± 0.025 ± 0.007. The uncertainties quoted correspond to the statistical and systematic uncertainties, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Lattice Symmetry-Guided Charge Transport in 2D Supramolecular Polymers Promotes Triplet Formation

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = m b /M h « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $αα^2_sL^k$ in the three-loop decay amplitude, where $L=ln(-M^2_h/m^2_b)$ and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $~αα^n_sL^{2n+1}$

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The light-ray OPE and conformal colliders

We derive a nonperturbative, convergent operator product expansion (OPE) for null-integrated operators on the same null plane in a CFT. The objects appearing in the expansion are light-ray operators, whose matrix elements can be computed by the generalized Lorentzian inversion formula. For example, a product of average null energy (ANEC) operators has an expansion in the light-ray operators that appear in the stress-tensor OPE. An important application is to collider event shapes. The light-ray OPE gives a nonperturbative expansion for event shapes in special functions that we call celestial blocks. As an example, we apply the celestial block expansion to energy-energy correlators in N = 4 Super Yang-Mills theory. Using known OPE data, we find perfect agreement with previous results both at weak and strong coupling, and make new predictions at weak coupling through 4 loops (NNNLO).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Measurement of hadronic event shapes in high-p T multijet final states at $ \sqrt{s} $ = 13 TeV with the ATLAS detector

A measurement of event-shape variables in proton-proton collisions at large momentum transfer is presented using data collected at $ \sqrt{s} $ = 13 TeV with the ATLAS detector at the Large Hadron Collider. Six event-shape variables calculated using hadronic jets are studied in inclusive multijet events using data corresponding to an integrated luminosity of 139 fb -1 . Measurements are performed in bins of jet multiplicity and in different ranges of the scalar sum of the transverse momenta of the two leading jets, reaching scales beyond 2 TeV. These measurements are compared with predictions from Monte Carlo event generators containing leading-order or next-to-leading order matrix elements matched to parton showers simulated to leading-logarithm accuracy. At low jet multiplicities, shape discrepancies between the measurements and the Monte Carlo predictions are observed. At high jet multiplicities, the shapes are better described but discrepancies in the normalisation are observed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Measurement of double-parton scattering in inclusive production of four jets with low transverse momentum in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A measurement of inclusive four-jet production in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. The transverse momenta of jets within |η| < 4.7 are required to exceed 35, 30, 25, and 20 GeV for the first-, second-, third-, and fourth-leading jet, respectively. Differential cross sections are measured as functions of the jet transverse momentum, jet pseudorapidity, and several other observables that describe the angular correlations between the jets. The measured distributions show sensitivity to different aspects of the underlying event, parton shower modeling, and matrix element calculations. In particular, the interplay between angular correlations caused by parton shower and double-parton scattering contributions is shown to be important. The double-parton scattering contribution is extracted by means of a template fit to the data, using distributions for single-parton scattering obtained from Monte Carlo event generators and a double-parton scattering distribution constructed from inclusive single-jet events in data. The effective double-parton scattering cross section is calculated and discussed in view of previous measurements and of its dependence on the models used to describe the single- parton scattering background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Even-point multi-loop unitarity and its applications: exponentiation, anomalies and evanescence

Abstract We identify novel structure in newly computed multi-loop amplitudes and quantum actions for even-point effective field theories, including both the nonlinear sigma model (NLSM) and double-copy gauge theories such as Born-Infeld and its supersymmetric generalizations. We exploit special properties of all even-point theories towards establishing an efficient unitarity based amplitude construction. In doing so, we find evidence that the leading IR divergence of NLSM amplitudes exponentiates when the target space is$$ \mathbbm{CP} $$ CP 1 ≅ SU(2)/U(1). We then systematically compute the two-loop anomalous behavior of Born-Infeld, and find that the counterterms needed to restore U(1) invariant behavior at loop-level can be constructed via a symmetric-structure double-copy. We also demonstrate that the divergent part of the one-minus (−+++) two-loop anomaly vanishes upon introducing an evanescent operator. In addition to these purely photonic counterterms, we verify through explicit calculation that the anomalous matrix elements that violate U(1) duality invariance can be alternatively cancelled by summing over internal$$ \mathcal{N} $$ N = 4 DBIVA superfields. Finally we find that$$ \mathcal{N} $$ N = 4 Dirac-Born-Infeld-Volkov-Akulov (DBIVA) amplitudes admit double-copy construction through two-loop order by reproducing our unitarity based result with a double copy between color-dual$$ \mathcal{N} $$ N = 4 super-Yang-Mills and our two-loop NLSM amplitudes. This result supports the possibility of color-dual representations for NLSM beyond one-loop. We conclude with an overview of howD-dimensional four-photon counterterms can be constructed in generality with the symmetric-structure double-copy, and outline a convenient way of counting evanescent operators using Hilbert series as generating functions.

Physics↗
DOE OSTI2024

A formalism for extracting track functions from jet measurements

The continued success of the jet substructure program will require widespread use of tracking information to enable increasingly precise measurements of a broader class of observables. The recent reformulation of jet substructure in terms of energy correlators has simplified the incorporation of universal non-perturbative matrix elements, so called “track functions”, in jet substructure calculations. These advances make it timely to understand how these universal non-perturbative functions can be extracted from hadron collider data, which is complicated by the use jet algorithms. In this paper we introduce a new class of jet functions, which we call (semi-inclusive) track jet functions, which describe measurements of the track energy fraction in identified jets. These track jet functions can be matched onto the universal track functions, with perturbatively calculable matching coefficients that incorporate the jet algorithm dependence. We perform this matching, and present phenomenological results for the charged energy fraction in jets at the LHC and EIC/HERA at collinear next-to-leading logarithmic accuracy. We show that higher moments of the charged energy fraction directly exhibit non-linear Lorentzian renormalization group flows, allowing the study of these flows with collider data. Our factorization theorem enables the extraction of universal track functions from jet measurements, opening the door to their use for a precision jet substructure program.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Measurement of $t\overline{t}$ production in association with additional b -jets in the eμ final state in proton–proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

This paper presents measurements of top-antitop quark pair ($t\overline{t}$) production in association with additional b-jets. The analysis utilises 140 fb -1 of proton–proton collision data collected with the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of 13 TeV. Fiducial cross-sections are extracted in a final state featuring one electron and one muon, with at least three or four b-jets. Results are presented at the particle level for both integrated cross-sections and normalised differential cross-sections, as functions of global event properties, jet kinematics, and b-jet pair properties. Observable quantities characterising b-jets originating from the top quark decay and additional b-jets are also measured at the particle level, after correcting for detector effects. The measured integrated fiducial cross sections are consistent with $t\overline{t}$$b\overline{b}$ predictions from various next-to-leading-order matrix element calculations matched to a parton shower within the uncertainties of the predictions. State of-the-art theoretical predictions are compared with the differential measurements; none of them simultaneously describes all observables. Differences between any two predictions are smaller than the measurement uncertainties for most observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Moments of axial-vector GPD from lattice QCD: quark helicity, orbital angular momentum, and spin-orbit correlation

In this work, we present a lattice QCD calculation of the Mellin moments of the twist-2 axial-vector generalized parton distribution (GPD), $\overset{\sim }{H}\left(x,\xi, t\right)$ , at zero skewness, ξ, with multiple values of the momentum transfer, t. Our analysis employs the short-distance factorization framework on ratio-scheme renormalized quasi-GPD matrix elements. The calculations are based on an N f = 2 + 1 + 1 twisted mass fermions ensemble with clover improvement, a lattice spacing of a = 0.093 fm, and a pion mass of m π = 260 MeV. We consider both the iso-vector and iso-scalar cases, utilizing next-to-leading-order perturbative matching while omitting the disconnected contributions and gluon mixing in the iso-scalar case. For the first time, we determine the Mellin moments of $\overset{\sim }{H}$ up to the fifth order. From these moments, we discuss the quark helicity and orbital angular momentum contributions to the nucleon spin, as well as the spin-orbit correlations of the quarks. Additionally, we perform a Fourier transform over the momentum transfer, which allows us to explore the spin structure in the impact-parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

An analytic result for the 0 → ggHHH amplitude

We present a fully analytic calculation of the leading-order one-loop amplitude for triple Higgs production via gluon fusion, gg → HHH, retaining full dependence on the mass of the heavy quark circulating in the loop. This amplitude provides a direct probe of the triple and quartic Higgs self-couplings, the measurement of which is a central goal of current and future colliders. The amplitude can be presented in compact form thanks to the use of analytic reconstruction techniques, based on finite-field and p-adic evaluations, multivariate partial fraction decompositions, and primary decompositions to identify common numerator factors. Although full analytic results are given in the text and in the supplementary material, the main thrust of this paper is to further test and illustrate these analytic reconstruction techniques in a concrete physical example. Our results provide a compact and efficient representation of the matrix element for this process, enabling evaluations that are more than an order of magnitude faster than existing numerical alternatives. Full analytic control of the leading-order, loop-induced amplitude is an important step towards handling more complex 2-loop or real-radiation corrections to this and related processes.

Campbell, John M. [Fermilab] (ORCID:00000002136193↗
DOE OSTI2022

Bootstrapping 2d ϕ 4 theory with Hamiltonian truncation data

We combine the methods of Hamiltonian Truncation and the recently proposed generalisation of the S-matrix bootstrap that includes local operators to determine the two-particle scattering amplitude and the two-particle form factor of the stress tensor at s > 0 in the 2d Φ 4 theory. We use the form factor of the stress tensor at s ≤ 0 and its spectral density computed using Lightcone Conformal Truncation (LCT), and inject them into the generalized S-matrix bootstrap set-up. The obtained results for the scattering amplitude and the form factor are fully reliable only in the elastic regime. We independently construct the “pure” S-matrix bootstrap bounds (bootstrap without including matrix elements of local operators), and find that the sinh-Gordon model and its analytic continuation the “staircase model” saturate these bounds. Surprisingly, the Φ 4 two-particle scattering amplitude also very nearly saturates these bounds, and moreover is extremely close to that of the sinh-Gordon/staircase model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Renormalization of the Standard Model Effective Field Theory from geometry

S-matrix elements are invariant under field redefinitions of the Lagrangian. They are determined by geometric quantities such as the curvature of the field-space manifold of scalar and gauge fields. We present a formalism where scalar and gauge fields are treated together, with a metric on the combined space of both types of fields. Scalar and gauge scattering amplitudes are given by the Riemann curvature Rijkl of this combined space, with indices i, j, k, l chosen to be scalar or gauge indices depending on the type of external particle. One-loop divergences can also be computed in terms of geometric invariants of the combined space, which greatly simplifies the computation of renormalization group equations. We apply our formalism to the Standard Model Effective Field Theory (SMEFT), and compute the renormalization group equations for even-parity bosonic operators to mass dimension eight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗