Search NASAβŒ• Search

Engineering topics

Clark, M. A.

Publications and source records attributed to Clark, M. A..

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Di-nucleons do not form bound states at heavy pion mass

We perform a high-statistics lattice QCD calculation of the low-energy two-nucleon scattering amplitudes. In order to address discrepancies in the literature, the calculation is performed at a heavy pion mass in the limit that the light quark masses are equal to the physical strange quark mass, $m_Ο€= m_K \simeq 714 $ MeV. Using a state-of-the-art momentum space method, we rule out the presence of a bound di-nucleon in both the isospin 0 (deuteron) and 1 (di-neutron) channels, in contrast with many previous results that made use of compact hexaquark creation operators. In order to diagnose the discrepancy, we add such hexaquark interpolating operators to our basis and find that they do not affect the determination of the two-nucleon finite volume spectrum, and thus they do not couple to deeply bound di-nucleons that are missed by the momentum-space operators. Further, we perform a high-statistics calculation of the HAL QCD potential on the same gauge ensembles and find qualitative agreement with our main results. We conclude that two-nucleons do not form bound states at heavy pion masses and that previous identification of deeply bound di-nucleons must have arisen from a misidentification of the spectrum from off-diagonal elements of a correlation function.

FOS: Physical sciences↗
DOE OSTI2022

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Scale setting the MΓΆbius domain wall fermion on gradient-flowed HISQ action using the omega baryon mass and the gradient-flow scales 𝑑 0 and 𝑀 0

We report on a subpercent scale determination using the omega baryon mass and gradient-flow methods. The calculations are performed on 22 ensembles of 𝑁 𝑓 =2 +1 +1 highly improved, rooted staggered sea-quark configurations generated by the MILC and CalLat Collaborations. The valence quark action used is MΓΆbius domain wall fermions solved on these configurations after a gradient-flow smearing is applied with a flowtime of 𝑑 gf = 1 in lattice units. The ensembles span four lattice spacings in the range 0.06 ≲ π‘Ž ≲0.15 fm, six pion masses in the range 130 ≲ π‘š πœ‹ ≲ 400 MeV and multiple lattice volumes. On each ensemble, the gradient-flow scales 𝑑 0 /π‘Ž 2 and 𝑀 0 /π‘Ž and the omega baryon mass π‘Žβ’π‘š Ξ© are computed. The dimensionless product of these quantities is then extrapolated to the continuum and infinite volume limits and interpolated to the physical light, strange and charm quark mass point in the isospin limit, resulting in the determination of $\sqrt{t}_0$ = 0.1422⁒(14) fm and 𝑀 0 = 0.1709⁒(11) fm with all sources of statistical and systematic uncertainty accounted for. The dominant uncertainty in both results is the stochastic uncertainty, though for $\sqrt{t}_0$ there are comparable continuum extrapolation uncertainties. For 𝑀 0 , there is a clear path for a few-per-mille uncertainty just through improved stochastic precision, as recently obtained by the Budapest-Marseille-Wuppertal Collaboration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with π‘š 𝑒⁒𝑑 β‰ˆ π‘š$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of π‘Ž β‰ˆ0.086 fm, lattice volume of 𝐿 = 48β’π‘Ž β‰ˆ 4.1 fm and pion mass π‘š πœ‹ β‰ˆ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

𝐹 𝐾 ⁑/𝐹 πœ‹ from MΓΆbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of 𝐹 𝐾 ⁑/𝐹 πœ‹ using MΓΆbius domain-wall fermions computed on gradient-flowed 𝑁 𝑓 =2 +1 +1 highly improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from 130 ≲ π‘š πœ‹ ≲ 400 MeV , four lattice spacings of π‘Ž ∼ 0.15, 0.12, 0.09 and 0.06 fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the π‘Ž ∼ 0.06 fm ensemble is helpful, but not necessary to achieve a subpercent determination of 𝐹 𝐾 ⁑/𝐹 πœ‹ . We also include an estimate of the strong isospin breaking corrections and arrive at a final result of 𝐹 $\hat{K}$ + ⁑ /𝐹 $\hat{πœ‹}$ + = 1.1942⁒(45) with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of |𝑉 𝑒⁒𝑠 |/|𝑉 𝑒⁒𝑑 | = 0.2311⁒(10).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Performance Portability of a Wilson Dslash Stencil Operator Mini-App Using Kokkos and SYCL

We describe our experiences in creating mini-apps for the Wilson-Dslash stencil operator for Lattice Quantum Chromodynamics using the Kokkos and SYCL programming models. In particular we comment on the performance achieved on a variety of hardware architectures, limitations we have reached in both programming models and how these have been resolved by us, or may be resolved by the developers of these models.

Joo, Balint↗
NASA NTRS1996

Suggestions for crops grown in controlled ecological life-support systems, based on attractive vegetarian diets

Assuming that crops grown in controlled ecological life-support systems (CELSS) should provide a basis for meals that are both nutritious and attractive (to taste and vision), and that CELSS diets on the moon or Mars or in space-craft during long voyages will have to be mostly vegetarian, a workshop was convened at the Johnson Space Center, Houston, Texas, U.S.A. on 19 to 21 January, 1994. Participants consisted of trained nutritionists and others; many of the approximately 18 presenters who discussed possible diets were practicing vegetarians, some for more than two decades. Considering all the presentations, seven conclusions (or points for discussion) could be formulated: nutritious vegetarian diets are relatively easily to formulate, vegetarian diets are healthy, variety is essential in vegetarian diets, some experiences (e.g., Bios-3 and Biosphere 2) are relevant to planning of CELSS diets, physical constraints will limit the choice of crops, a preliminary list of recommended crops can be formulated, and this line of research has some potential practical spinoffs. The list of crops and the reasons for including specific crops might be of interest to professionals in the field of health and nutrition as well as to those who are designing closed ecological systems.

NASA Program CELSS↗
NASA NTRS1988

An analysis of the response of the thermospheric normal modes to temporally varying convection at high latitudes

The response of the high-latitude thermosphere to temporal variations in the forcing due to the plasma convection pattern is investigated using a normal modes model that includes realistic temperature and electron density profiles as well as the effects of viscosity, heat conduction, and Lorentz forcing. The calculations show that the rotational wind component is only slightly modified when the plasma rotation is varied, while the divergent wind component responds more strongly to variations of the rotation rate, with the effect being greatest in the E region. A variation of the rotation rate then results in a greater magnitude and phase variation with the height of the divergent wind component in the E region, indicating that vertically propagating gravity waves are of importance in that height range.

Clark, M. A.β†—