Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Methods, software and datasets to verify DVH calculations against analytical values: Twenty years late(r)

Purpose: The authors designed data, methods, and metrics that can serve as a standard, independent of any software package, to evaluate dose-volume histogram (DVH) calculation accuracy and detect limitations. The authors use simple geometrical objects at different orientations combined with dose grids of varying spatial resolution with linear 1D dose gradients; when combined, ground truth DVH curves can be calculated analytically in closed form to serve as the absolute standards. Methods: DICOM RT structure sets containing a small sphere, cylinder, and cone were created programmatically with axial plane spacing varying from 0.2 to 3 mm. Cylinders and cones were modeled in two different orientations with respect to the IEC 1217 Y axis. The contours were designed to stringently but methodically test voxelation methods required for DVH. Synthetic RT dose files were generated with 1D linear dose gradient and with grid resolution varying from 0.4 to 3 mm. Two commercial DVH algorithms—PINNACLE (Philips Radiation Oncology Systems) and PlanIQ (Sun Nuclear Corp.)—were tested against analytical values using custom, noncommercial analysis software. In Test 1, axial contour spacing was constant at 0.2 mm while dose grid resolution varied. In Tests 2 and 3, the dose grid resolution was matched to varying subsampled axial contours with spacing of 1, 2, and 3 mm, and difference analysis and metrics were employed: (1) histograms of the accuracy of various DVH parameters (total volume, D{sub max}, D{sub min}, and doses to % volume: D99, D95, D5, D1, D0.03 cm{sup 3}) and (2) volume errors extracted along the DVH curves were generated and summarized in tabular and graphical forms. Results: In Test 1, PINNACLE produced 52 deviations (15%) while PlanIQ produced 5 (1.5%). In Test 2, PINNACLE and PlanIQ differed from analytical by >3% in 93 (36%) and 18 (7%) times, respectively. Excluding D{sub min} and D{sub max} as least clinically relevant would result in 32 (15%) vs 5 (2%) scored deviations for PINNACLE vs PlanIQ in Test 1, while Test 2 would yield 53 (25%) vs 17 (8%). In Test 3, statistical analyses of volume errors extracted continuously along the curves show PINNACLE to have more errors and higher variability (relative to PlanIQ), primarily due to PINNACLE’s lack of sufficient 3D grid supersampling. Another major driver for PINNACLE errors is an inconsistency in implementation of the “end-capping”; the additional volume resulting from expanding superior and inferior contours halfway to the next slice is included in the total volume calculation, but dose voxels in this expanded volume are excluded from the DVH. PlanIQ had fewer deviations, and most were associated with a rotated cylinder modeled by rectangular axial contours; for coarser axial spacing, the limited number of cross-sectional rectangles hinders the ability to render the true structure volume. Conclusions: The method is applicable to any DVH-calculating software capable of importing DICOM RT structure set and dose objects (the authors’ examples are available for download). It includes a collection of tests that probe the design of the DVH algorithm, measure its accuracy, and identify failure modes. Merits and applicability of each test are discussed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A STATISTICAL STUDY OF THE SOLAR WIND TURBULENCE AT ION KINETIC SCALES USING THE K-FILTERING TECHNIQUE AND CLUSTER DATA

Plasma turbulence at ion kinetic scales in the solar wind is investigated using the multi-point magnetometer data from the Cluster spacecraft. By applying the k-filtering method, we are able to estimate the full three-dimensional power spectral density P(ω{sub sc}, k) at a certain spacecraft frequency ω{sub sc} in wavevector (k) space. By using the wavevector at the maximum power in P(ω{sub sc}, k) at each sampling frequency ω{sub sc} and the Doppler shifted frequency ω{sub pla} in the solar wind frame, the dispersion plot ω{sub pla} = ω{sub pla}(k) is found. Previous studies have been limited to very few intervals and have been hampered by large errors, which motivates a statistical study of 52 intervals of solar wind. We find that the turbulence is predominantly highly oblique to the magnetic field k >> k {sub ∥}, and propagates slowly in the plasma frame with most points having frequencies smaller than the proton gyrofrequency ω{sub pla} < Ω{sub p}. Weak agreement is found that turbulence at the ion kinetic scales consists of kinetic Alfvén waves and coherent structures advected with plasma bulk velocity plus some minor more compressible components. The results suggest that anti-sunward and sunward propagating magnetic fluctuations are of similar nature in both the fast and slow solar wind at ion kinetic scales. The fast wind has significantly more anti-sunward flux than sunward flux and the slow wind appears to be more balanced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Impact of non-Maxwellian electron velocity distribution functions on inferred plasma parameters in collective Thomson scattering

Optical collective Thomson scattering provides precise density and temperature measurements in numerous plasma-physics experiments. The accuracy of such measurements depends on the core assumption that the underlying electron distribution functions in under-dense laser-produced plasmas are Maxwellian. A statistically based, quantitative analysis of the errors in the measured electron density and temperature is presented when synthetic data calculated using a non-Maxwellian electron distribution function is fit assuming a Maxwellian electron distribution. Such analysis can lead to errors of up to 50% in temperature and 30% in density, in the specific case of super-Gaussian distributions characteristic of inverse bremsstrahlung heating. Furthermore including the proper family of non-Maxwellian electron distribution functions, as a fitting parameter, in Thomson-scattering analysis removes the model-dependent errors in the inferred parameters at a minimal cost to the statistical uncertainty.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Workflow for Accelerating Multimodal Data Collection for Electrodeposited Films

Abstract Future machine learning strategies for materials process optimization will likely replace human capital-intensive artisan research with autonomous and/or accelerated approaches. Such automation enables accelerated multimodal characterization that simultaneously minimizes human errors, lowers costs, enhances statistical sampling, and allows scientists to allocate their time to critical thinking instead of repetitive manual tasks. Previous acceleration efforts to synthesize and evaluate materials have often employed elaborate robotic self-driving laboratories or used specialized strategies that are difficult to generalize. Herein we describe an implemented workflow for accelerating the multimodal characterization of a combinatorial set of 915 electroplated Ni and Ni–Fe thin films resulting in a data cube with over 160,000 individual data files. Our acceleration strategies do not require manufacturing-scale resources and are thus amenable to typical materials research facilities in academic, government, or commercial laboratories. The workflow demonstrated the acceleration of six characterization modalities: optical microscopy, laser profilometry, X-ray diffraction, X-ray fluorescence, nanoindentation, and tribological (friction and wear) testing, each with speedup factors ranging from 13–46x. In addition, automated data upload to a repository using FAIR data principles was accelerated by 64x.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Snowmass2021 - Letter of interest cosmology intertwined IV: The age of the universe and its curvature

A precise measurement of the curvature of the Universe is of prime importance for cosmology since it could not only confirm the paradigm of primordial inflation but also help in discriminating between different early-Universe scenarios. Recent observations, while broadly consistent with a spatially flat standard Λ Cold Dark Matter ( Λ CDM) model, show tensions that still allow (and, in some cases, even suggest) a few percent deviations from a flat universe. In particular, the Planck Cosmic Microwave Background power spectra, assuming the nominal likelihood, prefer a closed universe at more than 99% confidence level. While new physics could be at play, this anomaly may be the result of an unresolved systematic error or just a statistical fluctuation. However, since positive curvature allows a larger age of the Universe, an accurate determination of the age of the oldest objects provides a smoking gun in confirming or falsifying the current flat Λ CDM model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Evaluation and comparison of benchmark QSAR models to predict a relevant REACH endpoint: The bioconcentration factor (BCF)

The bioconcentration factor (BCF) is an important bioaccumulation hazard assessment metric in many regulatory contexts. Its assessment is required by the REACH regulation (Registration, Evaluation, Authorization and Restriction of Chemicals) and by CLP (Classification, Labeling and Packaging). We challenged nine well-known and widely used BCF QSAR models against 851 compounds stored in an ad-hoc created database. The goodness of the regression analysis was assessed by considering the determination coefficient (R{sup 2}) and the Root Mean Square Error (RMSE); Cooper's statistics and Matthew's Correlation Coefficient (MCC) were calculated for all the thresholds relevant for regulatory purposes (i.e. 100 L/kg for Chemical Safety Assessment; 500 L/kg for Classification and Labeling; 2000 and 5000 L/kg for Persistent, Bioaccumulative and Toxic (PBT) and very Persistent, very Bioaccumulative (vPvB) assessment) to assess the classification, with particular attention to the models' ability to control the occurrence of false negatives. As a first step, statistical analysis was performed for the predictions of the entire dataset; R{sup 2}>0.70 was obtained using CORAL, T.E.S.T. and EPISuite Arnot–Gobas models. As classifiers, ACD and log P-based equations were the best in terms of sensitivity, ranging from 0.75 to 0.94. External compound predictions were carried out for the models that had their own training sets. CORAL model returned the best performance (R{sup 2}{sub ext}=0.59), followed by the EPISuite Meylan model (R{sup 2}{sub ext}=0.58). The latter gave also the highest sensitivity on external compounds with values from 0.55 to 0.85, depending on the thresholds. Statistics were also compiled for compounds falling into the models Applicability Domain (AD), giving better performances. In this respect, VEGA CAESAR was the best model in terms of regression (R{sup 2}=0.94) and classification (average sensitivity>0.80). This model also showed the best regression (R{sup 2}=0.85) and sensitivity (average>0.70) for new compounds in the AD but not present in the training set. However, no single optimal model exists and, thus, it would be wise a case-by-case assessment. Yet, integrating the wealth of information from multiple models remains the winner approach. - Highlights: • REACH encourages the use of in silico methods in the assessment of chemicals safety. • The performances of nine BCF models were evaluated on a benchmark database of 851 chemicals. • We compared the models on the basis of both regression and classification performance. • Statistics on chemicals out of the training set and/or within the applicability domain were compiled. • The results show that QSAR models are useful as weight-of-evidence in support to other methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

BEYONDPLANCK II. CMB mapmaking through Gibbs sampling

We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Quasi-classical description of molecular dynamics based on Egorov's theorem

Egorov's theorem on the classical propagation of quantum observables is related to prominent quasi-classical descriptions of quantum molecular dynamics as the linearized semiclassical initial value representation, the Wigner phase space method, or the statistical quasiclassical method. The error estimates show that different accuracies are achievable for the computation of expectation values and position densities. Numerical experiments for a Morse model of diatomic iodine and confined Henon–Heiles systems in various dimensions illustrate the theoretical results.

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

Using Poisson-regularized inversion of Bremsstrahlung emission to extract full electron energy distribution functions from x-ray pulse-height detector data

We measured Electron Energy Distribution Functions (EEDFs) from below 200 eV to over 8 keV and spanning five orders-of-magnitude in intensity, produced in a low-power, RF-heated, tandem mirror discharge in the PFRC-II apparatus. The EEDF was obtained from the x-ray energy distribution function (XEDF) using a novel Poisson-regularized spectrum inversion algorithm applied to pulse-height spectra that included both Bremsstrahlung and line emissions. The XEDF was measured using a specially calibrated Amptek Silicon Drift Detector (SDD) pulse-height system with 125 eV FWHM at 5.9 keV. Finally, the algorithm is found to out-perform current leading x-ray inversion algorithms when the error due to counting statistics is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗
DOE OSTI2014

Weighing 'El Gordo' with a precision scale: Hubble space telescope weak-lensing analysis of the merging galaxy cluster ACT-CL J0102–4915 at z = 0.87

We present a Hubble Space Telescope weak-lensing study of the merging galaxy cluster 'El Gordo' (ACT-CL J0102–4915) at z = 0.87 discovered by the Atacama Cosmology Telescope (ACT) collaboration as the strongest Sunyaev-Zel'dovich decrement in its ∼1000 deg{sup 2} survey. Our weak-lensing analysis confirms that ACT-CL J0102–4915 is indeed an extreme system consisting of two massive (≳ 10{sup 15} M {sub ☉} each) subclusters with a projected separation of ∼0.7 h{sub 70}{sup −1} Mpc. This binary mass structure revealed by our lensing study is consistent with the cluster galaxy distribution and the dynamical study carried out with 89 spectroscopic members. We estimate the mass of ACT-CL J0102–4915 by simultaneously fitting two axisymmetric Navarro-Frenk-White (NFW) profiles allowing their centers to vary. We use only a single parameter for the NFW mass profile by enforcing the mass-concentration relation from numerical simulations. Our Markov-Chain-Monte-Carlo analysis shows that the masses of the northwestern (NW) and the southeastern (SE) components are M{sub 200c}=(1.38±0.22)×10{sup 15} h{sub 70}{sup −1} M{sub ⊙} and (0.78±0.20)×10{sup 15} h{sub 70}{sup −1} M{sub ⊙}, respectively, where the quoted errors include only 1σ statistical uncertainties determined by the finite number of source galaxies. These mass estimates are subject to additional uncertainties (20%-30%) due to the possible presence of triaxiality, correlated/uncorrelated large scale structure, and departure of the cluster profile from the NFW model. The lensing-based velocity dispersions are 1133{sub −61}{sup +58} km s{sup −1} and 1064{sub −66}{sup +62} km s{sup −1} for the NW and SE components, respectively, which are consistent with their spectroscopic measurements (1290 ± 134 km s{sup –1} and 1089 ± 200 km s{sup –1}, respectively). The centroids of both components are tightly constrained (∼4'') and close to the optical luminosity centers. The X-ray and mass peaks are spatially offset by ∼8'' (∼62 h{sub 70}{sup −1} kpc), which is significant at the ∼2σ confidence level. The mass peak, however, does not lead the gas peak in the direction expected if we are viewing the cluster soon after first core passage during a high speed merger. Under the assumption that the merger is happening in the plane of the sky, extrapolation of the two NFW halos to a radius r{sub 200a}=2.4 h{sub 70}{sup −1} Mpc yields a combined mass of M{sub 200a}=(3.13±0.56)×10{sup 15} h{sub 70}{sup −1} M{sub ⊙}. This extrapolated total mass is consistent with our two-component-based dynamical analysis and previous X-ray measurements, projecting ACT-CL J0102–4915 to be the most massive cluster at z > 0.6 known to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Correlated materials design: prospects and challenges

The design of correlated materials challenges researchers to combine the maturing, high throughput framework of DFT-based materials design with the rapidly-developing first-principles theory for correlated electron systems. We review the field of correlated materials, distinguishing two broad classes of correlation effects, static and dynamics, and describe methodologies to take them into account. We introduce a material design workflow, and illustrate it via examples in several materials classes, including superconductors, charge ordering materials and systems near an electronically driven metal to insulator transition, highlighting the interplay between theory and experiment with a view towards finding new materials. We review the statistical formulation of the errors of currently available methods to estimate formation energies. We formulate an approach for estimating a lower-bound for the probability of a new compound to form. Correlation effects have to be considered in all the material design steps. Furthermore, these include bridging between structure and property, obtaining the correct structure and predicting material stability. We introduce a post-processing strategy to take them into account.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: growth rate of structure measurement from cosmic voids

ABSTRACT We present a void clustering analysis in configuration-space using the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) DR16 samples. These samples consist of Luminous Red Galaxies (LRGs) combined with the high-redshift tail of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) DR12 CMASS galaxies (called as LRG+CMASS sample), Emission Line Galaxies (ELGs), and quasars (QSOs). We build void catalogues from the three eBOSS DR16 samples using a ZOBOV-based algorithm, providing 2814 voids, 1801 voids, and 4347 voids in the LRG+CMASS, ELG, and QSO samples, respectively, spanning the redshift range 0.6 < z < 2.2. We measure the redshift space distortions around voids using the anisotropic void-galaxy cross-correlation function and we extract the distortion parameter β. We test the methodology on realistic simulations before applying it to the data, and we investigate all our systematic errors on these mocks. We find βLRG(z = 0.74) = 0.415 ± 0.087, βELG(z = 0.85) = 0.665 ± 0.125 and βQSO(z = 1.48) = 0.313 ± 0.134, for the LRG+CMASS, ELG, and QSO sample, respectively. The quoted errors include systematic and statistical contributions. In order to convert our measurements in terms of the growth rate fσ8, we use consensus values of linear bias from the eBOSS DR16 companion papers, resulting in the following constraints: fσ8(z = 0.74) = 0.50 ± 0.11, fσ8(z = 0.85) = 0.52 ± 0.10, and fσ8(z = 1.48) = 0.30 ± 0.13. Our measurements are consistent with other measurements from eBOSS DR16 using conventional clustering techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

GUNDAM: a toolkit for fast spatial correlation functions in galaxy surveys

ABSTRACT We describe the capabilities of a new software package to calculate two-point correlation functions (2PCFs) of large galaxy samples. The code can efficiently estimate 3D/projected/angular 2PCFs with a variety of statistical estimators and bootstrap errors, and is intended to provide a complete framework (including calculation, storage, manipulation, and plotting) to perform this type of spatial analysis with large redshift surveys. Gundam implements a very fast skip list/linked list algorithm that efficiently counts galaxy pairs and avoids the computation of unnecessary distances. It is several orders of magnitude faster than a naive pair counter, and matches or even surpass other advanced algorithms. The implementation is also embarrassingly parallel, making full use of multicore processors or large computational clusters when available. The software is designed to be flexible, user friendly and easily extensible, integrating optimized, well-tested packages already available in the astronomy community. Out of the box, it already provides advanced features such as custom weighting schemes, fibre collision corrections and 2D correlations. Gundam will ultimately provide an efficient toolkit to analyse the large-scale structure ‘buried’ in upcoming extremely large data sets generated by future surveys.

Donoso, E.↗
DOE OSTI2020

Improved limits on Fierz interference using asymmetry measurements from the Ultracold Neutron Asymmetry (UCNA) experiment

The Ultracold Neutron Asymmetry (UCNA) experiment was designed to measure the β-decay asymmetry parameter, A 0 , for free neutron decay. In the experiment, polarized ultracold neutrons are transported into a decay trap, and their β-decay electrons are detected with ≈4π acceptance into two detector packages which provide position and energy reconstruction. The experiment also has sensitivity to b n , the Fierz interference term in the neutron β-decay rate. In this work, we determine bn from the energy dependence of A 0 using the data taken during the UCNA 2011–2013 run. In addition, we present the same type of analysis using the earlier 2010 A dataset. Furthermore, motivated by improved statistics and comparable systematic errors compared to the 2010 data-taking run, we present a new b n measurement using the weighted average of our asymmetry dataset fits, to obtain b n = 0.066 ± 0.041 stat ± 0.024 syst which corresponds to a limit of -0.012 < b n < 0.144 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗
DOE OSTI2026

Parton physics from a heavy-quark operator product expansion: Lattice QCD calculation of the fourth moment of the pion distribution amplitude

The pion light-cone distribution amplitude (LCDA) is an essential nonperturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of 𝑚 𝜋 = 550 MeV. This calculation finds ⟨𝜉 2 ⟩ = 0.202⁢(8)⁢(9) and ⟨𝜉 4 ⟩ = 0.039⁢(28)⁢(11) where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the $\overline{\textrm{MS}}$ scheme at a renormalization scale of 𝜇 = 2 GeV.

Detmold, William [Massachusetts Inst. of Technolog↗