Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-752: Three-Dimensional Carcinogenic Maps Induced by Photons and Protons

Purpose: Evaluate the excess relative risk (ERR) induced by photons and protons in each voxel of the lung, and display it as a three-dimensional map, known as the ERRM (i.e. excess relative risk map) along with the dose distribution map. In addition, we also study the effect of variations in the linear energy transfer (LET) distribution on ERRM for a given proton plan. Methods: The excess relative risk due to radiation is estimated using the initiation-inactivation-proliferation formalism. This framework accounts for three biological phenomenon: mutation induction, cell kill and proliferation. Cell kill and mutation induction are taken as a function of LET using experimental data. LET distributions are calculated using a Monte Carlo algorithm. ERR is then estimated for each voxel in the organ, and displayed as a three dimensional carcinogenic map. Results: The differences in the ERR’s between photons and protons is seen from the three-dimensional ERR map. In addition, we also varied the LET of a proton plan and observed the differences in the corresponding ERR maps demonstrating variations in the ERR maps depend on features of a proton plan. Additionally, our results suggest that any two proton plans that have the same integral dose does not necessarily imply identical ERR maps, and these changes are due to the variations in the LET distribution map. Conclusion: Clinically, it is important to have a three dimensional display of biological end points. This study is an effort to introduce 3D ERR maps into the treatment planning workflow for certain sites such as pediatric head and neck tumors.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2010

NICMOS Spectroscopy of HD 189733b

Spectral features corresponding to methane and water opacity were reported based on transmission spectroscopy of HD 189733b with Hubble/NICMOS. Recently, these data, and a similar data set for XO-1b, have been reexamined in Gibson et al. 2010, who claim they cannot reliably reproduce prior results. We examine the methods used by the Gibson team and identify two specific issues that could act to increase the formal uncertainties and to create instability in the minimization process. This would also be consistent with the GPA10 finding that they could not identify a way to select among the several instrument models they constructed. In the case of XO-1b, the Gibson team significantly changed the way in which the instrument model is defined (both with respect to the three approaches they used for HD 189733b, and the approach used by previous authors); this change, which omits the effect of the spectrum position on the detector, makes direct intercomparison of results difficult. In the experience of our group, the position of the spectrum on the detector is an important element of the instrument model because of the significant residual structure in the NICMOS spectral flat field. The approach of changing instrument models significantly complicates understanding the data reduction process and interpreting the results. Our team favors establishing a consistent method of handling NICMOS instrument systematic errors and applying it uniformly to data sets.

infrared↗
NASA NTRS2018

Qualifying Bulk Metallic Glass Gear Materials for Spacecraft Applications

NASA is evaluating bulk metallic glass (BMG) gears for extreme environment (cryogenic) applications; e.g., Europa Lander. The main purpose of this report is to recommend a quality assurance (QA) protocol for the production of consistent and reliable gear castings. Currently, there are two, separate manufacturers involved; Materion Corporation produces re-melt stock and Visser Precision Cast, Inc. produces BMG castings. Division of the existing alloy specifications document into material supplier-specific documents is proposed. Until a composition-based specification is established, processability- and performance-based methodologies may be an option during manufacturing development. It is recommended that each of the documents has a distinct focus; the processing behavior of crystalline feedstock, and the mechanical behavior of amorphous castings. The report addresses the “gray area” in QA testing between these materials suppliers by outlining an equitable division of responsibilities. The use of a number of dedicated reference dies to produce witness materials for physical and mechanical property evaluation is strongly recommended. Allocation of QA tests between screen testing of gear castings, and formal testing of witness materials, including cryogenic properties, is also proposed. The most effective suite of tests for both physical and mechanical property determination are down-selected from available national and international standards. The contents of this report create a solid basis for selection of QA tests, thereby setting the stage for negotiations between the materials suppliers and the customer.This technical paper:a) provides recommendations to advance materials specifications document JPL D-56223,Rev. A from “draft” to “preliminary” status;b) proposes division of the current document into two supplier-specific documents;c) separates testing methods by “ingots/castings” and “witness materials”;d) classifies test methodology based on “processability” and “performance”;e) assembles candidate ASTM specifications for “screen testing” and “formal testing”;f) identifies standard tests compatible with evaluation of BMG flexspline gears;g) suggests division of responsibilities for processability and performance testing; andh) addresses the QA “gray area” between production of ingots and castings.The report recommends that:a) processability- and performance-based QA methods might be used in lieu ofcompositional specifications during manufacturing development only;b) reference dies be designed for “acceptance” and “qualification” castings;c) witness materials provide physical and mechanical property data;d) formal testing be conducted at ambient and cryogenic temperatures; ande) processing/properties be correlated with compositional variations.

Hales, Stephen J.↗
DOE OSTI2025

Granger causal inference for climate change attribution

Abstract Climate change detection and attribution (D&A) is concerned with determining the extent to which anthropogenic activities have influenced specific aspects of the global climate system. D&A fits within the broader field of causal inference, the collection of statistical methods that identify cause and effect relationships. There are a wide variety of methods for making attribution statements, each of which require different types of input data and focus on different types of weather and climate events and each of which are conditional to varying extents. Some methods are based on Pearl causality (direct experimental interference) while others leverage Granger (predictive) causality, and the causal framing provides important context for how the resulting attribution conclusion should be interpreted. However, while Granger-causal attribution analyses have become more common, there is no clear statement of their strengths and weaknesses relative to Pearl-causal attribution and no clear consensus on where and when Granger-causal perspectives are appropriate. In this prospective paper, we provide a formal definition for Granger-based approaches to trend and event attribution and a clear comparison with more traditional methods for assessing the human influence on extreme weather and climate events. Broadly speaking, Granger-causal attribution statements can be constructed quickly from observations and do not require computationally-intesive dynamical experiments. These analyses also enable rapid attribution, which is useful in the aftermath of a severe weather event, and provide multiple lines of evidence for anthropogenic climate change when paired with Pearl-causal attribution. Confidence in attribution statements is increased when different methodologies arrive at similar conclusions. Moving forward, we encourage the D&A community to embrace hybrid approaches to climate change attribution that leverage the strengths of both Granger and Pearl causality.

Risser, Mark D. (ORCID:0000000319561783)↗
DOE OSTI2014

Analysis of the methods for the derivation of binary kinetic equations in the theory of fluorescence concentration quenching

In the framework of unified many-particle approach the familiar problem of fluorescence concentration quenching in the presence of pumping (light pulse) of arbitrary intensity is considered. This process is a vivid and the simplest example of multistage bulk reaction including bimolecular irreversible quenching reaction and reversible monomolecular transformation as elementary stages. General relation between the kinetics of multistage bulk reaction and that of the elementary stage of quenching has been established. This allows one to derive general kinetic equations (of two types) for the multistage reaction in question on the basis of general kinetic equations (differential and integro-differential) of elementary stage of quenching. Relying on the same unified many-particle approach we have developed binary approximations with the use of two (frequently employed in the literature) many-particle methods (such as simple superposition approximation and the method of extracting pair channels in three-particle correlation evolution) to the derivation of non-Markovian binary kinetic equations. The possibility of reducing the obtained binary equations to the Markovian equations of formal chemical kinetics has been considered. As an example the exact solution of the problem (for the specific case) is examined, and the applicability of two many particle methods of derivation of binary equations is analyzed.

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

Considering nonlocality in the optical potentials within eikonal models

Background: For its simplicity, the eikonal method is the tool of choice to analyze nuclear reactions at high energies (E > 100 MeV/nucleon), including knockout reactions. However, so far, the effective interactions used in this method are assumed to be fully local. Purpose: Given the recent studies on nonlocal optical potentials, in this work we assess whether nonlocality in the optical potentials is expected to impact reactions at high energies and then explore different avenues for extending the eikonal method to include nonlocal interactions. Method: We compare angular distributions obtained for nonlocal interactions (using the exact R-matrix approach for elastic scattering and the adiabatic distorted wave approximation for transfer) with those obtained using their local-equivalent interactions. Results: Our results show that transfer observables are significantly impacted by nonlocality in the high-energy regime. Because knockout reactions are dominated by stripping (transfer to inelastic channels), nonlocality is expected to have a large effect on knockout observables too. Three approaches are explored for extending the eikonal method to nonlocal interactions, including an iterative method and a perturbation theory. Conclusions: None of the derived extensions of the eikonal model provide a good description of elastic scattering. Here, this paper suggests that nonlocality removes the formal simplicity associated with the eikonal model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Modeling techniques for quantum cascade lasers

Quantum cascade lasers are unipolar semiconductor lasers covering a wide range of the infrared and terahertz spectrum. Lasing action is achieved by using optical intersubband transitions between quantized states in specifically designed multiple-quantum-well heterostructures. A systematic improvement of quantum cascade lasers with respect to operating temperature, efficiency, and spectral range requires detailed modeling of the underlying physical processes in these structures. Moreover, the quantum cascade laser constitutes a versatile model device for the development and improvement of simulation techniques in nano- and optoelectronics. This review provides a comprehensive survey and discussion of the modeling techniques used for the simulation of quantum cascade lasers. The main focus is on the modeling of carrier transport in the nanostructured gain medium, while the simulation of the optical cavity is covered at a more basic level. Specifically, the transfer matrix and finite difference methods for solving the one-dimensional Schrödinger equation and Schrödinger-Poisson system are discussed, providing the quantized states in the multiple-quantum-well active region. The modeling of the optical cavity is covered with a focus on basic waveguide resonator structures. Furthermore, various carrier transport simulation methods are discussed, ranging from basic empirical approaches to advanced self-consistent techniques. The methods include empirical rate equation and related Maxwell-Bloch equation approaches, self-consistent rate equation and ensemble Monte Carlo methods, as well as quantum transport approaches, in particular the density matrix and non-equilibrium Green's function formalism. The derived scattering rates and self-energies are generally valid for n-type devices based on one-dimensional quantum confinement, such as quantum well structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

SU-G-JeP2-02: A Unifying Multi-Atlas Approach to Electron Density Mapping Using Multi-Parametric MRI for Radiation Treatment Planning

Purpose: MRI has a number of advantages over CT as a primary modality for radiation treatment planning (RTP). However, one key bottleneck problem still remains, which is the lack of electron density information in MRI. In the work, a reliable method to map electron density is developed by leveraging the differential contrast of multi-parametric MRI. Methods: We propose a probabilistic Bayesian approach for electron density mapping based on T1 and T2-weighted MRI, using multiple patients as atlases. For each voxel, we compute two conditional probabilities: (1) electron density given its image intensity on T1 and T2-weighted MR images, and (2) electron density given its geometric location in a reference anatomy. The two sources of information (image intensity and spatial location) are combined into a unifying posterior probability density function using the Bayesian formalism. The mean value of the posterior probability density function provides the estimated electron density. Results: We evaluated the method on 10 head and neck patients and performed leave-one-out cross validation (9 patients as atlases and remaining 1 as test). The proposed method significantly reduced the errors in electron density estimation, with a mean absolute HU error of 138, compared with 193 for the T1-weighted intensity approach and 261 without density correction. For bone detection (HU>200), the proposed method had an accuracy of 84% and a sensitivity of 73% at specificity of 90% (AUC = 87%). In comparison, the AUC for bone detection is 73% and 50% using the intensity approach and without density correction, respectively. Conclusion: The proposed unifying method provides accurate electron density estimation and bone detection based on multi-parametric MRI of the head with highly heterogeneous anatomy. This could allow for accurate dose calculation and reference image generation for patient setup in MRI-based radiation treatment planning.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2004

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗
DOE OSTI2026

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗
NASA NTRS1997

Algorithms for Performance, Dependability, and Performability Evaluation using Stochastic Activity Networks

Modeling tools and technologies are important for aerospace development. At the University of Illinois, we have worked on advancing the state of the art in modeling by Markov reward models in two important areas: reducing the memory necessary to numerically solve systems represented as stochastic activity networks and other stochastic Petri net extensions while still obtaining solutions in a reasonable amount of time, and finding numerically stable and memory-efficient methods to solve for the reward accumulated during a finite mission time. A long standing problem when modeling with high level formalisms such as stochastic activity networks is the so-called state space explosion, where the number of states increases exponentially with size of the high level model. Thus, the corresponding Markov model becomes prohibitively large and solution is constrained by the the size of primary memory. To reduce the memory necessary to numerically solve complex systems, we propose new methods that can tolerate such large state spaces that do not require any special structure in the model (as many other techniques do). First, we develop methods that generate row and columns of the state transition-rate-matrix on-the-fly, eliminating the need to explicitly store the matrix at all. Next, we introduce a new iterative solution method, called modified adaptive Gauss-Seidel, that exhibits locality in its use of data from the state transition-rate-matrix, permitting us to cache portions of the matrix and hence reduce the solution time. Finally, we develop a new memory and computationally efficient technique for Gauss-Seidel based solvers that avoids the need for generating rows of A in order to solve Ax = b. This is a significant performance improvement for on-the-fly methods as well as other recent solution techniques based on Kronecker operators. Taken together, these new results show that one can solve very large models without any special structure.

Deavours, Daniel D.↗
DOE OSTI2015

Superposition of Fragment Excitations for Excited States of Large Clusters with Application to Helium Clusters

We develop a local excited-state method, based on the configuration interaction singles (CIS) wave function, for large atomic and molecular clusters. This method exploits the properties of absolutely localized molecular orbitals (ALMOs), which strictly limits the total number of excitations, and results in formal scaling with the third power of the system size for computing the full spectrum of ALMO-CIS excited states. Here, the derivation of the equations and design of the algorithm are discussed in detail, with particular emphasis on the computational scaling. Clusters containing ~500 atoms were used in evaluating the scaling, which agrees with the theoretical predictions, and the accuracy of the method is evaluated with respect to standard CIS. Finally, a pioneering application to the size dependence of the helium cluster spectrum is also presented for clusters of 25-231 atoms, the largest of which results in the computation of 2310 excited states per sampled cluster geometry.

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

Effects of overlapping sources on cosmic shear estimation: Statistical sensitivity and pixel-noise bias

The next generation of dark-energy imaging surveys — so called “Stage-IV” surveys, such as that of the Rubin Observatory Legacy Survey of Space and Time (LSST) — will cross a threshold in the number density of detected sources on the sky that requires qualitatively different image analysis and measurement techniques compared to the current generation of Stage-III surveys. In Stage-IV surveys, a significant amount of the cosmologically useful information is due to sources whose images overlap with those of other sources on the sky. Here, we focus on the weak gravitational lensing probe, for which we expect the largest impact since the cosmic shear signal is primarily encoded in the estimated shapes of observed galaxies and thus directly impacted by overlaps. We introduce a framework based on the Fisher formalism to analyze the effect of the overlapping sources (“blending”) on the estimation of cosmic shear. This method gives concrete predictions for the minimum loss of information due to noise and blending for any choice of “deblending” scheme and shape-measurement algorithm. Our studies account for undetected sources but do not address their full effects and biases they may introduce. We use simulated images and predict this impact of blending for three surveys: the Dark Energy Survey (DES), the Hyper-Suprime Cam Subaru Strategic Program (HSC-SSP), and the Rubin LSST. Our methodology successfully estimates the statistical sensitivity to weak lensing for DES and HSC early results. For LSST, we present the expected loss in statistical sensitivity for the ten-year survey due to blending. We find that for approximately 62% of galaxies that are likely to be detected in full-depth LSST images, at least 1% of the flux in their pixels is from overlapping sources. We also find that the statistical correlations between measures of overlapping galaxies and, to a much lesser extent (0.2%) the higher shot noise level due to their presence, decrease the effective number density of galaxies, N eff , by ~ 18%. We calculate an upper limit on N eff of 39.4 galaxies per arcmin 2 in r band. We study the impact of stars on as a function of stellar density and illustrate the diminishing returns of extending the survey into lower Galactic latitudes. We extend the simulation-based Fisher formalism to predict the expected increase in pixel-noise bias due to blending for maximum-likelihood (ML) shape estimators. We find that noise bias depends sensitively on the particular shape estimator and measure of ensemble-average shape that is used, and properties of the galaxy that include redshift-dependent quantities such as size and luminosity. The source code for these studies is available online.[The documented software developed for the catalog-level studies are available in the open-source LSST DESC github repository https://github.com/LSSTDESC/WeakLensingDeblending. The software for analyzing one or two galaxies with user-defined parameters is in the open-source github repository https://github.com/ismael-mendoza/ShapeMeasurementFisherFormalism.]

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Quantum conductance of zigzag graphene oxide nanoribbons

The electronic properties of zigzag graphene oxide nanoribbons (ZGOR) are presented. The results show interesting behaviors which are considerably different from the properties of the perfect graphene nanoribbons (GNRs). The theoretical methods include a Huckel-tight binding approach, a Green's function methodology, and the Landauer formalism. The presence of oxygen on the edge results in band bending, a noticeable change in density of states and thus the conductance. Consequently, the occupation in the valence bands increase for the next neighboring carbon atom in the unit cell. Conductance drops in both the conduction and valence band regions are due to the reduction of allowed k modes resulting from band bending. The asymmetry of the energy band structure of the ZGOR is due to the energy differences of the atoms. The inclusion of a foreign atom's orbital energies changes the dispersion relation of the eigenvalues in energy space. These novel characteristics are important and valuable in the study of quantum transport of GNRs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Real and complex valued geometrical optics inverse ray-tracing for inline field calculations

In this work, a 3-D ray based model for computing laser fields in dissipative and amplifying media is presented. The eikonal equation is solved using inverse ray-tracing on a dedicated nonstructured 3-D mesh. Inverse ray-tracing opens the possibility of using Complex Geometrical Optics (CGO), for which we propose a propagation formalism in a finite element mesh. Divergent fields at caustics are corrected using an etalon integral method for fold-type caustics. This method is successfully applied in dissipative media by modifying the ray-ordering and root selection rules, thereby allowing one to reconstruct the field in the entire caustic region. In addition, we demonstrate how caustics in the CGO framework can disappear entirely for sufficiently dissipative media, making the complex ray approach valid in the entire medium. CGO is shown to offer a more precise modeling of laser refraction and absorption in a dissipative medium when compared to Geometrical Optics (GO). In the framework of Inertial Confinement Fusion (ICF), this occurs mostly at intermediate temperatures or at high temperatures close to the critical density. Additionally, GO is invalid at low temperatures if an approximated expression of the permittivity is used. The inverse ray-tracing algorithm for GO and CGO is implemented in the IFRIIT code, in the framework of a dielectric permittivity described in 3-D using a piecewise linear approximation in tetrahedrons. Fields computed using GO and CGO are compared to results from the electromagnetic wave solver LPSE. Excellent agreement is obtained in 1-D linear and nonlinear permittivity profiles. Good agreement is also obtained for ICF-like Gaussian density profiles in 2-D. Finally, we demonstrate how the model reproduces Gaussian beam diffraction using CGO. The IFRIIT code will be interfaced inline to 3-D radiative hydrodynamic codes to describe the nonlinear laser plasma interaction in ICF and high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The separate universe approach to soft limits

We develop a formalism for calculating soft limits of n -point inflationary correlation functions using separate universe techniques. Our method naturally allows for multiple fields and leads to an elegant diagrammatic approach. As an application we focus on the trispectrum produced by inflation with multiple light fields, giving explicit formulae for all possible single- and double-soft limits. We also investigate consistency relations and present an infinite tower of inequalities between soft correlation functions which generalise the Suyama-Yamaguchi inequality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

QCD Predictions for Physical Multimeson Scattering Amplitudes

We use lattice QCD calculations of the finite-volume spectra of systems of two and three mesons to determine, for the first time, three-particle scattering amplitudes with physical quark masses. Our results are for combinations of 𝜋 + and 𝐾 + , at a lattice spacing 𝑎 = 0.063 fm, and in the isospin-symmetric limit. We also obtain accurate results for maximal-isospin two-meson amplitudes, with those for 𝜋 + ⁢𝐾 + and 2⁢𝐾 + being the first determinations at the physical point. Dense lattice spectra are obtained using the stochastic Laplacian-Heaviside method, and the analysis leading to scattering amplitudes is done using the relativistic finite-volume formalism. Results are compared to chiral perturbation theory and to phenomenological fits to experimental data, finding good agreement.

hadron-hadron interactions↗