Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice calculation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Molecular tail chemistry controls thermal transport in fullerene films

Herein we report on the thermal conductivities of alkyl- and indene-group functionalized fullerene derivative thin films as measured via time domain thermoreflectance and steady state thermoreflectance. The thermal conductivities vary from 0.064 ± 0.007 W m – 1 K – 1 for [6,6]-phenyl C 61 -butyric acid methyl ester (PCBM) to 0.15 ± 0.017 W m – 1 K – 1 for indene- C 60 bisadduct at room temperature and do not exhibit significant temperature dependence from 300 to 375 K. In comparison to the thermal conductivity of PCBM, increasing the length of the alkyl chain, as in the case of [6,6]-phenyl C 61 -butyric acid butyl ester, and [6,6]-phenyl C 61 -butyric acid octyl ester leads to higher thermal conductivities. Likewise, increasing the number of alkyl chains attached to the fullerenes as in the case of bisadduct PCBM leads to a higher thermal conductivity compared to that of PCBM. We present atomistic insights into the role of chemical functionalization on the overall heat transfer in these fullerene derivatives by conducting molecular dynamics simulations and lattice dynamics calculations. The thermal conductivities predicted via our atomistic simulations qualitatively agree with the experimental trends for our fullerene derivatives. Lattice dynamics calculations reveal that one of the main factors dictating the ultralow thermal conductivities in fullerene derivatives is the large reduction in modal diffusivities in the molecular crystals as calculated from the Allen-Feldman model, thus providing an explanation for their largely reduced thermal conductivities as compared to that of C 60 crystals. The low diffusivities result from high degrees of localization of Einstein-like vibrations in fullerene derivatives due to the molecular side chains, providing the ability to dial-in the properties of these low thermal conductivity solids via molecular engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗
DOE OSTI2023

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗
DOE OSTI2020

Weak neutral current axial form factor using $ \left(\overline{\nu}\right)\nu $-nucleon scattering and lattice QCD inputs

We present a determination of the neutral current weak axial charge $G$ $^{Z}_{A}$(0) = –0.654(3) stat (5) sys using the strange quark axial charge $G$ $^{s}_{A}$(0) calculated with lattice QCD. We then perform a phenomenological analysis, where we combine the strange quark electromagnetic form factor from lattice QCD with (anti)neutrino-nucleon scattering differential cross section from MiniBooNE experiments in a momentum transfer region 0.24 ≲ Q 2 ≲ 0.71 GeV 2 to determine the neutral current weak axial form factor $G$ $^{Z}_{A}$(Q 2 ) in the range of 0 ≲ Q 2 ≤ 1 GeV 2 . This yields a phenomenological value of $G$ $^{Z}_{A}$(0) = –0.687(89)stat(40)sys. The value of $G$ $^{Z}_{A}$(0) constrained by the lattice QCD calculation of $G$ $^{s}_{A}$(0), when compared to its phenomenological determination, provides a significant improvement in precision and accuracy and can be used to provide a constraint on the fit to $G$ $^{Z}_{A}$(Q 2 ) for Q 2 > 0. This constrained fit leads to an unambiguous determination of (anti)neutrino-nucleon neutral current elastic scattering differential cross section near Q 2 = 0 and can play an important role in numerically isolating nuclear effects in this region. We show a consistent description of $G$ $^{Z}_{A}$(Q 2 ) obtained from the (anti)neutrino-nucleon scattering cross section data requires a nonzero contribution of the strange quark electromagnetic form factor. We demonstrate the robustness of our analysis by providing a post-diction of the BNL E734 experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Utility of a hybrid approach to the hadronic vacuum polarization contribution to the muon anomalous magnetic moment

An accurate calculation of the leading-order hadronic vacuum polarization (LOHVP) contribution to the anomalous magnetic moment of the muon ( a μ ) is key to determining whether a discrepancy, suggesting new physics, exists between the Standard Model and experimental results. This calculation can be expressed as an integral over Euclidean time of a current-current correlator G ( t ) , where G ( t ) can be calculated using lattice QCD or, with dispersion relations, from experimental data for e + e − → hadrons . The BMW/DMZ collaboration recently presented a hybrid approach in which G ( t ) is calculated using lattice QCD for most of the contributing t range, but using experimental data for the largest t (lowest energy) region. Here we study the advantages of varying the position t = t 1 separating lattice QCD from data-driven contributions. The total LOHVP contribution should be independent of t 1 , providing both a test of the experimental input and the robustness of the hybrid approach. We use this criterion and a correlated fit to show that Fermilab/HPQCD/MILC lattice QCD results from 2019 strongly favor the CMD-3 cross section data for e + e − → π + π − over a combination of earlier experimental results for this channel. Further, the resulting total LOHVP contribution obtained is consistent with the result obtained by BMW/DMZ, and supports the scenario in which there is no significant discrepancy between the experimental value for a μ and that expected in the Standard Model. We then discuss how improved lattice results in this hybrid approach could provide a more accurate total LOHVP across a wider range of t 1 values with an uncertainty that is smaller than that from either lattice QCD or data-driven approaches on their own. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The anomalous magnetic moment of the muon in the Standard Model: an update

We present the current Standard Model (SM) prediction for the muon anomalous magnetic moment, a μ , updating the first White Paper (WP20) [1]. The pure QED and electroweak contributions have been further consolidated, while hadronic contributions continue to be responsible for the bulk of the uncertainty of the SM prediction. Significant progress has been achieved in the hadronic light-by-light scattering contribution using both the data-driven dispersive approach as well as lattice-QCD calculations, leading to a reduction of the uncertainty by almost a factor of two. The most important development since WP20 is the change in the estimate of the leading-order hadronic-vacuum-polarization (LO HVP) contribution. A new measurement of the e + e - → π + π - cross section by CMD-3 has increased the tensions among data-driven dispersive evaluations of the LO HVP contribution to a level that makes it impossible to combine the results in a meaningful way. At the same time, the attainable precision of lattice-QCD calculations has increased substantially and allows for a consolidated lattice-QCD average of the LO HVP contribution with a precision of about 0.9%. Adopting the latter in this update has resulted in a major upward shift of the total SM prediction, which now reads $a^{SM}_{μ}$ = 116 592 033 (62) x $10^{-11}$ (530 ppb). When compared against the current experimental average based on the E821 experiment and runs 1–6 of E989 at Fermilab, one finds $a^{exp}_{μ} -a^{SM}_{μ}= 38 (63)$ x $10^{-11}$, which implies that there is no tension between the SM and experiment at the current level of precision. The final precision of E989 (127 ppb) is the target of future efforts by the Theory Initiative. The resolution of the tensions among data-driven dispersive evaluations of the LO HVP contribution will be a key element in this endeavor.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Skewness-dependent moments of the pion GPD from nonlocal quark-bilinear correlators

We present lattice QCD calculations of the odd Mellin moments of pion valence-quark generalized parton distribution up to fifth order ⟨𝑥 4 ⟩ and for the skewness range [−0.33, 0] using operator product expansion of bilocal quark-bilinear operators. The calculations are performed on an ensemble with lattice spacing 𝑎 = 0.04 fm and valence pion mass 300 MeV, employing boosted pion states with momenta up to 2.428 GeV and momentum transfers reaching 2.748 GeV 2 . We employ ratio-scheme renormalization and next-to-leading logarithmic resummed perturbative matching. At zero skewness, our results are consistent with previous lattice studies. By combining matrix elements at multiple values of skewness and momentum transfer, skewness-dependent moments are obtained through simultaneous polynomiality-constrained fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Methodology for Sensitivity Analysis of Homogenized Cross-Sections to Instantaneous and Historical Lattice Conditions with Application to AP1000® PWR Lattice

In the two-step method for nuclear reactor simulation, lattice physics calculations are performed to compute homogenized cross-sections for a variety of burnups and lattice configurations. A nodal code is then used to perform full-core analysis using the pre-calculated homogenized cross-sections. One source of uncertainty introduced in this method is that the lattice configuration or depletion conditions typically do not match a pre-calculated one from the lattice physics simulations. Therefore, some interpolation model must be used to estimate the homogenized cross-sections in the nodal code. This current study provides a methodology for sensitivity analysis to quantify the impact of state variables on the homogenized cross-sections. This methodology also allows for analyses of the historical effect that the state variables have on homogenized cross-sections. An application of this methodology on a lattice for the Westinghouse AP1000® reactor is presented where coolant density, fuel temperature, soluble boron concentration, and control rod insertion are the state variables of interest. The effects of considering the instantaneous values of the state variables, historical values of the state variables, and burnup-averaged values of the state variables are analyzed. Using these methods, it was found that a linear model that only considers the instantaneous and burnup-averaged values of state variables can fail to capture some variations in the homogenized cross-sections.

Price, Dean (ORCID:0000000309990111)↗
DOE OSTI2025

TMD Handbook

This handbook provides a comprehensive review of transverse-momentum-dependent parton distribution functions and fragmentation functions, commonly referred to as transverse momentum distributions (TMDs). TMDs describe the distribution of partons inside the proton and other hadrons with respect to both their longitudinal and transverse momenta. They provide unique insight into the internal momentum and spin structure of hadrons, and are a key ingredient in the description of many collider physics cross sections. Understanding TMDs requires a combination of theoretical techniques from quantum field theory, nonperturbative calculations using lattice QCD, and phenomenological analysis of experimental data. The handbook covers a wide range of topics, from theoretical foundations to experimental analyses, as well as recent developments and future directions. It is intended to provide an essential reference for researchers and graduate students interested in understanding the structure of hadrons and the dynamics of partons in high energy collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Monolayer-like lattice dynamics in bulk WSe 2

Understanding the microscopic lattice dynamics is essential for regulating the thermal properties in two-dimensional layered materials. In transition metal dichalcogenides, the layered structures result in different but closely related phonon dispersions between monolayer and bulk. Here, by combining inelastic X-ray scattering and first-principles calculations, the lattice dynamics of tungsten diselenide (WSe 2 ) was investigated comprehensively, and a monolayer-like lattice dynamics in the bulk WSe 2 was revealed. We performed the first temperature-dependent phonon dispersions measurement and obtained the mode Grüneisen parameters of bulk WSe 2 , which are found to be in better agreement with the calculations on the monolayer system than those of the bulk. This observation indicates that lattice dynamics in bulk WSe 2 hold the characteristic of monolayers. The result provides valuable insights into the thermal properties of WSe 2 -based devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Site-Specific Structure at Multiple Length Scales in Kagome Quantum Spin Liquid Candidates

Realizing a quantum spin liquid (QSL) ground state in a real material is a leading issue in condensed matter physics research. In this pursuit, it is crucial to fully characterize the structure and influence of defects, as these can significantly affect the fragile QSL physics. In this work, we perform a variety of cutting-edge synchrotron x-ray scattering and spectroscopy techniques, and we advance new methodologies for site-specific diffraction and L-edge Zn absorption spectroscopy. The experimental results along with our first-principles calculations address outstanding questions about the local and long-range structures of the two leading kagome QSL candidates, Zn-substituted barlowite Cu 3 Zn x Cu 1-x (OH) 6 FBr and herbertsmithite Cu 3 Zn(OH) 6 Cl 2 . On all length scales probed, there is no evidence that Zn substitutes onto the kagome layers, thereby preserving the QSL physics of the kagome lattice. Our calculations show that antisite disorder is not energetically favorable and is even less favorable in Zn-barlowite compared to herbertsmithite. Site-specific x-ray diffraction measurements of Zn-barlowite reveal that Cu 2+ and Zn 2+ selectively occupy distinct interlayer sites, in contrast to herbertsmithite. Using the first measured Zn L-edge inelastic x-ray absorption spectra combined with calculations, we discover a systematic correlation between the loss of inversion symmetry from pseudo-octahedral (herbertsmithite) to trigonal prismatic coordination (Zn-barlowite) with the emergence of a new peak. Overall, our measurements suggest that Zn-barlowite has structural advantages over herbertsmithite that make its magnetic properties closer to an ideal QSL candidate: its kagome layers are highly resistant to nonmagnetic defects while the interlayers can accommodate a higher amount of Zn substitution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward a resolution of the NN controversy

Lattice QCD calculations of two-nucleon interactions have been underway for about a decade, but still haven't reached the pion mass regime necessary for matching onto effective field theories and extrapolating to the physical point. Furthermore, results from different methods, including the use of the Lüscher formalism with different types of operators, as well as the HALQCD potential method, do not agree even qualitatively at very heavy pion mass. We investigate the role that different operators employed in the literature may play on the extraction of spectra for use within the Lüscher method. We first explore expectations from Effective Field Theory solved within a finite volume, for which the exact spectrum may be computed given different physical scenarios. We then present preliminary lattice QCD results for two-nucleon spectra calculated using different operators on a common lattice ensemble.

Nicholson, A↗
DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗