SEARCH · Search NASA

Results for “lattice quantum chromodynamics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

New state of matter between the hadronic phase and the quark-gluon plasma?

Lattice-quantum chromodynamics (QCD) simulations and theoretical arguments hint at the existence of an intermediate phase of strongly interacting matter between a confined hadron gas and a deconfined quark-gluon plasma (QGP). We qualitatively and semiquantitatively explore and differentiate the phase structures in the temperature window from the QCD pseudocritical temperature 𝑇 c ≃ 160 MeV to the pure gluonic deconfinement temperature 𝑇 d ≃ 285 MeV. We propose a three-regime picture using a hadron resonance gas description augmented with the exponential spectrum of strings, corresponding to highly excited mesons and glueballs, based on the analysis of a large number 𝑁 c of colors. We estimate the entropy density from our model to confirm that the lattice-QCD data are bracketed with three regimes, i.e., a hadron gas, a QGP, and a new phase for 𝑇 c ≲ 𝑇 ≲ 𝑇 d . In this new phase, which we name a spaghetti of quarks with glueballs (SQGBs), thermal degrees of freedom of quarks are liberated, yet gluons remain confined in glueballs. Since the Hagedorn temperature 𝑇 H ∼ 285 MeV is universal in the meson and the glueball sectors, in the infinite-𝑁 c limit, the phase diagram in the plane of the baryon chemical potential and the temperature is reduced to one with the confined and deconfined phases and quarkyonic matter at high density. At large but finite 𝑁 c , an SQGB window may open between these phases. We point out that the SQGB has interesting similarities with quarkyonic matter and that this matter in the large-𝑁 c limit is confined as measured by the interaction between heavy quarks, but behaves in other respects like a quasifree gas of quarks. As a result of the extrapolation to 𝑁 c = 3, we present a revised phase diagram with the SQGB phase bounded by thermal crossovers. Finally, we give a quantitative analysis of chiral-symmetry restoration in the SQGB phase.

Color confinement↗

DOE OSTI2024

Calculating the two-photon exchange contribution to K L → μ + μ − decay

We present a theoretical framework within which both the real and imaginary parts of the complex, two-photon exchange amplitude contributing to K L → μ + μ − decay can be calculated using lattice quantum chromodynamics. The real part of this two-photon amplitude is of approximately the same size as that coming from a second-order weak strangeness-changing neutral-current process. Thus a test of the standard model prediction for this second-order weak process depends on an accurate result of this two-photon amplitude. A limiting factor of our proposed method comes from low-energy three-particle π π γ states. The contribution from these states will be significantly distorted by the finite volume of our calculation—a distortion for which there is no available correction. However, a simple estimate of the contribution of these three-particle states suggests their contribution to be at most a few percent allowing their neglect in a lattice calculation with a 10% target accuracy. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2023

Generalized parton distributions through universal moment parameterization: non-zero skewness case

We present the first global analysis of generalized parton distributions (GPDs) combing lattice quantum chromodynamics (QCD) calculations and experiment measurements including global parton distribution functions (PDFs), form factors (FFs) and deeply virtual Compton scattering (DVCS) measurements. Following the previous work where we parameterize GPDs in terms of their moments, we extend the framework to allow for the global analysis at non-zero skewness. Together with the constraints at zero skewness, we fit GPDs to global DVCS measurements from both the recent JLab and the earlier Hadron-Electron Ring Accelerator (HERA) experiments with two active quark flavors and leading order QCD evolution. With certain choices of empirical constraints, both sea and valence quark distributions are extracted with the combined inputs, and we present the quark distributions in the proton correspondingly. We also discuss how to extend the framework to accommodate more off-forward constraints beyond the small ξ expansion, especially the lattice calculated GPDs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Towards a real-time computation of timelike hadronic vacuum polarization and light-by-light scattering: Schwinger Model tests

Hadronic vacuum polarization (HVP) and light-by-light scattering (HLBL) are crucial for evaluating the Standard Model predictions concerning the muon’s anomalous magnetic moment. However, direct first-principle lattice gauge theory-based calculations of these observables in the timelike region remain challenging. Discrepancies persist between lattice quantum chromodynamics (QCD) calculations in the spacelike region and dispersive approaches relying on experimental data parametrization from the timelike region. Here, we introduce a methodology employing 1+1-dimensional quantum electrodynamics (QED), i.e. the Schwinger Model, to investigate the HVP and HLBL. To that end, we use both tensor network techniques, specifically matrix product states, and classical emulators of digital quantum computers. Demonstrating feasibility in a simplified model, our approach sets the stage for future endeavors leveraging digital quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Two-neutrino double-𝛽 decay in pionless effective field theory from a Euclidean finite-volume correlation function

Two-neutrino double-β decay of certain nuclear isotopes is one of the rarest Standard Model processes observed in nature. Its neutrinoless counterpart is an exotic lepton-number nonconserving process that is widely searched to determine if the neutrinos are Majorana fermions. In order to connect the rate of these processes to the Standard Model and beyond the Standard Model interactions, it is essential that the corresponding nuclear matrix elements are constrained reliably from theory. Lattice quantum chromodynamics (LQCD) and low-energy effective field theories (EFTs) are expected to play an essential role in constraining the matrix element of the two-nucleon subprocess, which could in turn provide the input into ab initio nuclear-structure calculations in larger isotopes. Focusing on the two-neutrino process $nn\rightarrow pp(ee\bar{v}_{e}\bar{v}_e)$, the amplitude is constructed in this work in pionless EFT at next-to-leading order, demonstrating the emergence of a renormalization-scale independent amplitude and the absence of any new low-energy constant at this order beyond those present in the single-weak process. Most importantly, it is shown how a LQCD four-point correlation function in Euclidean and finite-volume spacetime can be used to constrain the Minkowski infinite-volume amplitude in the EFT. The same formalism is provided for the related single-weak process, which is an input to the double-β decay formalism. The LQCD-EFT matching procedure outlined for the double-weak amplitude paves the road toward constraining the two-nucleon matrix element entering the neutrinoless double-β decay amplitude with a light Majorana neutrino.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Status of Lattice QCD Determination of Nucleon Form Factors and Their Relevance for the Few-GeV Neutrino Program

Calculations of neutrino–nucleus cross sections begin with the neutrino–nucleon interaction, making the latter critically important to flagship neutrino oscillation experiments despite limited measurements with poor statistics. Alternatively, lattice quantum chromodynamics (LQCD) can be used to determine these interactions from the Standard Model with quantifiable theoretical uncertainties. Recent LQCD results of g A are in excellent agreement with data, and results for the (quasi-)elastic nucleon form factors with full uncertainty budgets are expected within a few years. We review the status of the field and LQCD results for the nucleon axial form factor, F A ( Q 2 ), a major source of uncertainty in modeling sub-GeV neutrino–nucleon interactions. Results from different LQCD calculations are consistent but collectively disagree with existing models, with potential implications for current and future neutrino oscillation experiments. We describe a road map to solidify confidence in the LQCD results and discuss future calculations of more complicated processes, which are important to few-GeV neutrino oscillation experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Kinematically enhanced interpolating operators for boosted hadrons

We propose to use interpolating operators for lattice quantum chromodynamics calculations of highly boosted pions and nucleons with kinematically enhanced ground-state overlap factors at large momentum. Because this kinematic enhancement applies to the signal but not the variance of the correlation function, these interpolating operators can achieve better signal-to-noise ratios at large momentum. We perform proof-of-principle calculations with boosted pions and nucleons using close-to-physical and larger quark masses to explore the utility of our proposal. Results for effective energies and matrix elements, as well as Lanczos ground-state energy estimators, are consistent with theoretical expectations for signal-to-noise improvement at large momenta.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Empirical investigation of nuclear correlation function distributions in lattice QCD

Two-point correlation functions of systems with baryon number 𝐵 ∈ {1,2,3,4} are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the 𝑂⁡(𝑁) model. The theoretical 𝑂⁡(𝑁) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number 𝐵 QCD correlation function at large temporal separation is well reproduced by the 𝑂⁡(𝑁 ∼ 2/𝐵) model distribution.

Lattice field theory↗

DOE OSTI2023

Measurement of the branching fraction of ${\textrm{D}}_{\textrm{s}}^{+}$ → τ + ν τ via τ + → μ+ν μ $\overline{\nu}$τ

Utilizing 7 . 33 fb - 1 of e + e - collision data taken at the center-of-mass energies of 4.128, 4.157, 4.178, 4.189, 4.199, 4.209, 4.219, and 4.226 GeV with the BESIII detector, the branching fraction of the leptonic decay ${D}_s^{+}$ → τ + ν τ via τ + → μ + ν μ $\overline{\nu}$ τ is measured to be ${\mathcal{B}}_{D_s^{+}\to {\tau}^{+}{\nu}_{\tau }}$ = (5 . 37 ± 0 . 17 stat ± 0 . 15 syst )%. Combining this branching fraction with the world averages of the measurements of the masses of τ + and ${D}_s^{+}$ as well as the lifetime of ${D}_s^{+}$, we extract the product of the decay constant of ${D}_s^{+}$ and the c → s Cabibbo-Kobayashi-Maskawa matrix element to be ${f}_{D_s^{+}}$ |V cs | = (246 . 7 ± 3 . 9 stat ± 3 . 6 syst ) MeV. Taking |V cs | from a global fit in the standard model we obtain ${f}_{D_s^{+}}$ = (253 . 4 ± 4 . 0 stat ± 3 . 7 syst ) MeV. Conversely, taking ${f}_{D_s^{+}}$ from lattice quantum chromodynamics calculations, we obtain |V cs | = 0 . 987 ± 0 . 016 stat ± 0 . 014 syst .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Closing in on critical net-baryon fluctuations at LHC energies: Cumulants up to third order in Pb–Pb collisions

Fluctuation measurements are important sources of information on the mechanism of particle production at LHC energies. This article reports the first experimental results on third-order cumulants of the net-proton distributions in Pb–Pb collisions at a center-of-mass energy $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector. The results on the second-order cumulants of net-proton distributions at $\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV are also discussed in view of effects due to the global and local baryon number conservation. The results demonstrate the presence of long-range rapidity correlations between protons and antiprotons. Such correlations originate from the early phase of the collision. The experimental results are compared with HIJING and EPOS model calculations, and the dependence of the fluctuation measurements on the phase-space coverage is examined in the context of lattice quantum chromodynamics (LQCD) and hadron resonance gas (HRG) model estimations. The measured third-order cumulants are consistent with zero within experimental uncertainties of about 4% and are described well by LQCD and HRG predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Measurement of the axial vector form factor from antineutrino–proton scattering

Scattering of high energy particles from nucleons probes their structure, as was done in the experiments that established the non-zero size of the proton using electron beams. The use of charged leptons as scattering probes enables measuring the distribution of electric charges, which is encoded in the vector form factors of the nucleon. Scattering weakly interacting neutrinos gives the opportunity to measure both vector and axial vector form factors of the nucleon, providing an additional, complementary probe of their structure. The nucleon transition axial form factor, F A , can be measured from neutrino scattering from free nucleons, ν μ n → μ – p and ν¯ μ p → μ + n , as a function of the negative four-momentum transfer squared (Q 2 ). Up to now, FA(Q 2 ) has been extracted from the bound nucleons in neutrino–deuterium scattering, which requires uncertain nuclear corrections. Here we report the first high-statistics measurement, to our knowledge, of the ν¯ μ p → μ + n cross-section from the hydrogen atom, using the plastic scintillator target of the MINERvA experiment, extracting F A from free proton targets and measuring the nucleon axial charge radius, r A , to be 0.73 ± 0.17 fm. The antineutrino–hydrogen scattering presented here can access the axial form factor without the need for nuclear theory corrections, and enables direct comparisons with the increasingly precise lattice quantum chromodynamics computations. Finally, the tools developed for this analysis and the result presented are substantial advancements in our capabilities to understand the nucleon structure in the weak sector, and also help the current and future neutrino oscillation experiments to better constrain neutrino interaction models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

A regression algorithm for accelerated lattice QCD that exploits sparse inference on the D-Wave quantum annealer

We propose a regression algorithm that utilizes a learned dictionary optimized for sparse inference on a D-Wave quantum annealer. In this regression algorithm, we concatenate the independent and dependent variables as a combined vector, and encode the high-order correlations between them into a dictionary optimized for sparse reconstruction. On a test dataset, the dependent variable is initialized to its average value and then a sparse reconstruction of the combined vector is obtained in which the dependent variable is typically shifted closer to its true value, as in a standard inpainting or denoising task. Here, a quantum annealer, which can presumably exploit a fully entangled initial state to better explore the complex energy landscape, is used to solve the highly non-convex sparse coding optimization problem. The regression algorithm is demonstrated for a lattice quantum chromodynamics simulation data using a D-Wave 2000Q quantum annealer and good prediction performance is achieved. The regression test is performed using six different values for the number of fully connected logical qubits, between 20 and 64. The scaling results indicate that a larger number of qubits gives better prediction accuracy.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Lossy compression of statistical data using quantum annealer

Abstract We present a new lossy compression algorithm for statistical floating-point data through a representation learning with binary variables. The algorithm finds a set of basis vectors and their binary coefficients that precisely reconstruct the original data. The optimization for the basis vectors is performed classically, while binary coefficients are retrieved through both simulated and quantum annealing for comparison. A bias correction procedure is also presented to estimate and eliminate the error and bias introduced from the inexact reconstruction of the lossy compression for statistical data analyses. The compression algorithm is demonstrated on two different datasets of lattice quantum chromodynamics simulations. The results obtained using simulated annealing show 3–3.5 times better compression performance than the algorithm based on neural-network autoencoder. Calculations using quantum annealing also show promising results, but performance is limited by the integrated control error of the quantum processing unit, which yields large uncertainties in the biases and coupling parameters. Hardware comparison is further studied between the previous generation D-Wave 2000Q and the current D-Wave Advantage system. Our study shows that the Advantage system is more likely to obtain low-energy solutions for the problems than the 2000Q.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

From lattice QCD to in-medium heavy-quark interactions via deep learning

Bottomonium states are key probes for experimental studies of the quark-gluon plasma (QGP) created in high-energy nuclear collisions. Theoretical models of bottomonium productions in high-energy nuclear collisions rely on the in-medium interactions between the bottom and antibottom quarks, which can be characterized by real (V R (T, r)) and imaginary (V I (T, r)) potentials, as functions of temperature and spatial separation. Recently, the masses and thermal widths of up to 3S and 2P bottomonium states in QGP were calculated using lattice quantum chromodynamics (LQCD). Starting from these LQCD results and through a novel application of deep neural network (DNN), here, we obtain model-independent results for V R (T, r) and V I (T, r). The temperature dependence of V R (T, r) was found to be very mild between T ≈ 0 - 330 MeV. Meanwhile, VI(T, r) shows rapid increase with T and r, which is much larger than the perturbation theory based expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Finalizing Transition to the New Data Center at BNL

Computational science, data management and analysis have been key factors in the success of Brookhaven National Laboratory's scientific programs at the Relativistic Heavy Ion Collider (RHIC), the National Synchrotron Light Source (NSLS-II), the Center for Functional Nanomaterials (CFN), and in biological, atmospheric, and energy systems science, Lattice Quantum Chromodynamics (LQCD) and Materials Science, as well as our participation in international research collaborations, such as the ATLAS Experiment at Europe's Large Hadron Collider (LHC) at CERN (Switzerland) and the Belle II Experiment at KEK (Japan). The construction of a new data center is an acknowledgement of the increasing demand for computing and storage services at BNL in the near term and enable the Lab to address the needs of the future experiments at the High-Luminosity LHC at CERN and the Electron-Ion Collider (EIC) at BNL in the long term. The Computing Facility Revitalization (CFR) project is aimed at repurposing the former National Synchrotron Light Source (NSLS-I) building as the new data center for BNL. The construction of the new data center was finished in 2021Q3, and it was delivered for production in early FY2022 for all collaborations supported by the Scientific Data and Computing Center (SDCC), including STAR, PHENIX and sPHENIX experiments at RHIC collider at BNL, the Belle II Experiment at KEK (Japan), and the Computational Science Initiative at BNL (CSI). This paper highlights the key mechanical, electrical, and networking components of the new data center in its final configuration as used in production since 2021Q4 and gives an overview for the extension of the central network systems into the new data center and the migration of a significant portion of IT load and services from the old data center to the new data center carried out in 20212023, with expected completion of the main phase of the gradual IT equipment replacement and migration from the old data center into the new one set to the end of FY2023 (Sep 30, 2023).

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Toward the first gluon parton distribution from the LaMET

Abstract We present progress towards the first unpolarized gluon quasi-parton distribution function (PDF) from lattice quantum chromodynamics using high-statistics measurements for hadrons at two valence pion masses M π ≈ 310 and 690 MeV computed on an a ≈ 0.12 fm ensemble with 2 + 1 + 1-flavors of highly improved staggered quark generated by the MILC collaboration. In this study, we consider two gluon operators for which the hybrid-ratio renormalization matching kernels have been recently derived and a third operator that has been used in prior pseudo-PDF studies of the gluon PDFs. We compare the matrix elements for each operator for both the nucleon and pion, at both pion masses, and using two gauge-smearing techniques. Focusing on the more phenomenologically studied nucleon gluon PDF, we compare the ratio and hybrid-ratio renormalized matrix elements at both pion masses and both smearings to those reconstructed from the nucleon gluon PDF from the CT18 global analysis. We identify the best choice of operator to study the gluon PDF and present the first gluon quasi-PDF under some caveats. Additionally, we explore the recent idea of Coulomb gauge fixing to improve signal at large Wilson-line displacement and find it could be a major help in improving the signal in the gluon matrix elements. This work helps identify the best operator for studying the gluon quasi-PDF, shows higher hadron boost momentum is needed to implement hybrid-ratio renormalization reliably, and suggests the need to study more diverse set of operators with their corresponding perturbative calculations for hybrid-ratio renormalization to further gluon quasi-PDF study.

Good, William (ORCID:0000000184281003)↗

DOE OSTI2025

Sensitivity of transverse momentum correlations to early-stage and thermal fluctuations

Transverse momentum correlations were recently measured by the ALICE collaboration at the LHC [Acharya et al (2020 Phys. Lett. B 804 135375)]. A long-range structure in terms of relative pseudorapidity of particle pairs is observed. This may imply some signal of the initial state owing to the sheer spread of the correlation. However, the fluctuations inside a thermally equilibrated medium have to be taken into account, serving as motivation for this paper. Using lattice Quantum Chromodynamics constraints, we predicted the development and spread of balancing correlations caused by energy-momentum conservation. At the same time, we propagated the Gaussian-shaped initial correlation using hydrodynamics to estimate its effects. Our findings suggest that the resulting correlation, known as ‘the ridge,’ is sensitive both to flucutations seeded in the pre-equilibrium stage and to those seeded in the equilibrated medium. This can provide important insight into the early stages of the collision.

ALICE↗