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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Lattice QCD Benchmark of Proton Helicity and Flavor-Dependent Unpolarized Transverse Momentum-Dependent Parton Distribution Functions at Physical Quark Masses

We present the first lattice QCD calculations of the isovector helicity transverse momentum-dependent parton distribution function (TMDPDF) and the flavor-dependent unpolarized TMDPDFs for up and down quarks in the proton. Our computations utilize domain-wall fermion discretization with physical quark masses. Employing Coulomb-gauge-fixed quark correlation functions within the large-momentum effective theory framework, we access nonperturbative transverse quark separations 𝑏 𝑇 up to approximately 1 fm, corresponding to transverse momenta as low as 200 MeV. Based on the quasi-TMD factorization theorem, we construct renormalization-group–invariant ratios that are equal to the corresponding light-cone TMDPDF ratios. At moderate 𝑥, our results reveal that the isovector helicity and unpolarized TMDPDFs exhibit nearly identical transverse structure up to a normalization factor, and the unpolarized distributions display only mild flavor dependence. These findings not only support key trends observed in recent global analyses but also provide robust, nonperturbative constraints that can distinguish between different parameterizations. This Letter establishes a first-principles benchmark for TMDPDFs, offering valuable input for ongoing and future experimental efforts to map the proton’s three-dimensional structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

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. 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, 𝑚 𝜋 = 𝑚 𝐾 ≃ 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. 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. Furthermore, 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 di-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.

Physics - Physics of elementary particles and fiel↗
DOE OSTI2024

Separation of infrared and bulk in thermal QCD

A new thermal regime of QCD, featuring decoupled scale-invariant infrared glue, has been proposed to exist both in pure-glue (N f =0) and “real-world” (N f =2+1 at physical quark masses) QCD. In this IR phase, elementary degrees of freedom flood the infrared, forming a distinct component independent from the bulk. This behavior necessitates non-analyticities in the theory. In pure-glue QCD, such non-analyticities have been shown to arise via Anderson-like mobility edges in Dirac spectra (λ IR = 0, ±λ A ≠ 0), as manifested in the dimension function d IR (λ). Here, we present the first evidence, based on lattice QCD calculation at a=0.105 fm, that this mechanism is also at work in real-world QCD, thus supporting the existence of the proposed IR regime in nature. An important aspect of our results is that, while at T = 234 MeV we find a dimensional jump between zero modes and lowest near-zero modes very close to unity (d IR = 3 to d IR ≃ 2), similar to the IR phase of pure-glue QCD, at T = 187 MeV we observe a continuous λ-dependence. This suggests that thermal states just above the chiral crossover are non-analytically (in T) connected to thermal state at T = 234 MeV, supporting the key original proposition that the transition into the IR regime occurs at a temperature strictly above the chiral crossover.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Signs of nonmonotonic finite-volume corrections to 𝑔 𝐴

We study finite-volume (FV) corrections to determinations of 𝑔 𝐴 via lattice quantum chromodynamics (QCD) using analytic results and numerical analysis. We observe that 𝑆⁢𝑈⁡(2) heavy Baryon chiral perturbation theory does not provide an unambiguous prediction for the sign of the FV correction, which is not surprising when one also considers large-𝑁 𝑐 constraints on the axial couplings. We further show that nonmonotonic FV corrections are naturally allowed when one considers either including explicit Δ-resonance degrees of freedom or one works to higher orders in the chiral expansion. We investigate the potential impact of these FV corrections with a precision study of 𝑔 𝐴 using models of FV corrections that are monotonic and nonmonotonic. Using lattice QCD data that is approximately at the 1% level of precision, we do not see significant evidence of nonmonotonic corrections. Looking forward to the next phase of lattice QCD calculations, we estimate that calculations that are between the 0.1% and 1% level of precision may be sensitive to these FV artifacts. Finally, we present an update of the CalLat prediction of 𝑔 𝐴 in the isospin limit with subpercent precision, 𝑔$^{QCD}_{𝐴}$ = 1.2674⁢(96).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Electronic coherence evolution at the nearly commensurate-incommensurate CDW boundary of 1⁢T−TaS 2

Transition-metal dichalcogenides host a variety of charge-density-wave phases that couple lattice, charge, and correlation effects. In 1⁢𝑇−TaS 2 , the commensurate and nearly commensurate states are well characterized, yet the transition near 350 K into the incommensurate phase has lacked direct momentum-resolved insight. In this work, we use temperature-dependent angle-resolved photoemission spectroscopy to track the electronic structure across this transition. We observe a suppression of quasiparticle spectral weight at the Brillouin-zone center, coincident with the transport anomaly, but without clear evidence of a full band-gap opening. The transition appears to involve momentum-dependent redistribution of spectral weight, consistent with a loss of coherence that reshapes the Fermi surface while leaving conduction dispersions largely intact. These results suggest that the nearly commensurate–incommensurate transition may not align with a conventional metal-insulator transition picture, but rather as an electronic reconstruction driven by loss of coherence. Our work provides new microscopic insight into the resistivity anomaly near room temperature and may guide design principles for collective electronic switching in transition-metal dichalcogenides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hadronic vacuum polarization for the muon g - 2 from lattice QCD: Complete short and intermediate windows

We present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment a μ in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to a μ , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain a μ l l , SD ( conn ) = 48.139 ( 11 ) stat ( 91 ) syst [ 92 ] total × 10 - 10 and a μ l l , W ( conn ) = 206.90 ( 14 ) stat ( 61 ) syst [ 63 ] total × 10 - 10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are a μ SD = 69.05 ( 1 ) stat ( 21 ) syst [ 21 ] total × 10 - 10 and a μ W = 236.45 ( 17 ) stat ( 83 ) syst [ 85 ] total × 10 - 10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hadronic vacuum polarization for the muon 𝑔 − 2 from lattice QCD: Complete short and intermediate windows

We present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment 𝑎 𝜇 in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to 𝑎 𝜇 , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain 𝑎$^{𝑙⁢𝑙,SD}_{𝜇⁡}$(conn) = 48.139⁢(11) stat ⁢(91) syst ⁢[92] total × 10 −10 and 𝑎$^{𝑙⁢𝑙,W}_{𝜇}$⁡(conn) = 206.90⁢(14) stat ⁢(61) syst ⁢[63] total × 10 −10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are 𝑎$^{SD}_{𝜇}$ = 69.05⁢(1) stat⁢ (21) syst ⁢[21] total × 10 −10 and 𝑎$^{W}_{𝜇}$ = 236.45⁢(17) stat⁢ (83) syst ⁢[85] total × 10 −10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Semiconductor Thermal and Electrical Properties Decoupled by Localized Phonon Resonances

Abstract Thermoelectric materials convert heat into electricity through thermally driven charge transport in solids or vice versa for cooling. To compete with conventional energy‐conversion technologies, a thermoelectric material must possess the properties of both an electrical conductor and a thermal insulator. However, these properties are normally mutually exclusive because of the interconnection between scattering mechanisms for charge carriers and phonons. Recent theoretical investigations on sub‐device scales have revealed that nanopillars attached to a membrane exhibit a multitude of local phonon resonances, spanning the full spectrum, that couple with the heat‐carrying phonons in the membrane and cause a reduction in the in‐plane thermal conductivity, with no expected change in the electrical properties because the nanopillars are outside the pathway of voltage generation and charge transport. Here this effect is demonstrated experimentally for the first time by investigating device‐scale suspended silicon membranes with GaN nanopillars grown on the surface. The nanopillars cause up to 21% reduction in the thermal conductivity while the power factor remains unaffected, thus demonstrating an unprecedented decoupling in the semiconductor's thermoelectric properties. The measured thermal conductivity behavior for coalesced nanopillars and corresponding lattice‐dynamics calculations provide evidence that the reductions are mechanistically tied to the phonon resonances. This finding paves the way for high‐efficiency solid‐state energy recovery and cooling.

Spann, Bryan T.↗
DOE OSTI2020

Dialing in Direct Air Capture of CO 2 by Crystal Engineering of Bis-iminoguanidines

Direct air capture (DAC) technologies that extract carbon dioxide from the atmosphere via chemical processes have the potential to restore the atmospheric CO 2 concentration to an optimal level. This study elucidates structure-property relationships in DAC by crystallization of bis-iminoguanidine (BIG) carbonate salts. Here, their crystal structures are analyzed by X-ray and neutron diffraction to accurately measure key structural parameters including molecular conformations, hydrogen bonding, and π-stacking. Experimental measurements of key properties, such as aqueous solubilities and regeneration energies and temperatures, are complemented by first-principles calculations of lattice and hydration free energies, as well as free energies of reactions with CO 2 , and BIG regenerations. Minor structural modifications in the molecular structure of the BIGs are found to result in major changes in the crystal structures and the aqueous solubilities within the series, leading to enhanced DAC.

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

Electroweak three-body decays in the presence of two- and three-body bound states

Recently, formalism has been derived for studying electroweak transition amplitudes for three-body systems both in infinite and finite volumes. The formalism provides exact relations that the infinite-volume amplitudes must satisfy, as well as a relationship between physical amplitudes and finite-volume matrix elements, which can be constrained from lattice QCD calculations. This formalism poses additional challenges when compared with the analogous well-studied two-body equivalent one, including the necessary step of solving integral equations of singular functions. In this work, we provide some non-trivial analytical and numerical tests on the aforementioned formalism. In particular, we consider a case where the three-particle system can have three-body bound states as well as bound states in the two-body subsystem. For kinematics below the three-body threshold, we demonstrate that the scattering amplitudes satisfy unitarity. We also check that for these kinematics the finite-volume matrix elements are accurately described by the formalism for two-body systems up to exponentially suppressed corrections. Finally, we verify that in the case of the three-body bound state, the finite-volume matrix element is equal to the infinite-volume coupling of the bound state, up to exponentially suppressed errors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Influence of isoelectronic mass modulation on phonon decay and transport in layered pnictogen-carbon systems

This study systematically investigates the structural, electronic, and thermal transport properties of layered pnictogen–carbon (Pn 2 C 2 ; Pn = P, As, Sb, Bi) monolayers, revealing a profound influence of isoelectronic mass modulation on phonon dynamics and thermal conductivity. Through first-principles calculations and lattice dynamics analysis, we demonstrate that increasing the atomic mass of pnictogens leads to elongated Pn-C bonds, weakened interatomic interactions, and enhanced phonon-phonon scattering, resulting in a dramatic reduction in lattice thermal conductivity (κ). Specifically, κ decreases exponentially from 65.6 W/m-K in P 2 C 2 to an ultralow 0.37 W/m-K in Bi 2 C 2 , driven by suppressed acoustic-optical phonon gaps Δ$^{Γ}_{A-0}$, increased anharmonicity, and reduced phonon lifetimes and group velocities. The transition from semiconducting to metallic behavior in heavier pnictogen-based systems further enhances phonon-electron scattering, contributing to thermal conductivity suppression. Structural analysis highlights the strengthening of out-of-plane C—C bonds and the contraction of C—Pn—C bond angles, which disrupt in-plane phonon transport. Further, these findings establish a design framework for engineering ultralow-κ materials through chemical substitution and structural tuning, offering significant potential for thermoelectric and thermal management applications. This work provides fundamental insights into phonon decay mechanisms and thermal transport in 2D materials, paving the way for advanced phononic and energy-efficient devices.

2D materials↗
DOE OSTI2022

Large negative thermal expansion in GdFe(CN) 6 driven by unusual low-frequency modes

Understanding the negative thermal expansion (NTE) mechanism is of great importance. In this work, we consider the new NTE compound GdFe(CN) 6 (αv = -34.2×10 -6 K -1 ) as a case study to investigate the NTE mechanism from the perspective of the lattice vibrational dynamics. The atomic mean-square displacements suggest that the NTE of GdFe(CN) 6 comes from the strong tension effect induced by the transverse vibrations of the atomic –Fe–Ctriple bondN–Gd– linkages, with the largest contribution given by N atoms. Lattice dynamics calculations show that three low-frequency optical modes at about 50 cm -1 show the largest negative Grüneisen parameters thus providing the largest contribution to the NTE. In conclusion, the existence of these unusual low-frequency vibrational modes can be ascribed to the presence of GdN 6 trigonal prisms in the framework structure of GdFe(CN) 6 .

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

First principles prediction of the Al-Li phase diagram including configurational and vibrational entropic contributions

The whole Al-Li phase diagram is predicted from first principles calculations and statistical mechanics including the effect of configurational and vibrational entropy. The formation enthalpy of different configurations at different temperatures was accurately predicted by means of cluster expansions that were fitted from first principles calculations. The vibrational entropic contribution of each configuration was determined from the bond length vs. bond stiffness relationships for each type of bond and the Gibbs free energy of the different phases was obtained as a function of temperature from Monte Carlo simulations. The predicted phase diagram was in excellent agreement with the currently accepted experimental one in terms of the stable (AlLi, Al 2 Li 3 , AlLi 2 , Al 4 Li 9 ) and metastable (Al 3 Li) phases, of the phase boundaries between them and of the maximum stability temperature of line compounds. In addition, it provided accurate information about the gap between Al 3 Li and AlLi solvus lines. Finally, the influence of the vibrational entropy on the correct prediction of the phase diagram is discussed. Overall, the methodology shows that accurate phase diagrams of alloys of technological interest can be predicted from first principles calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING↗
DOE OSTI2022

Transverse momentum structure of proton within the basis light-front quantization framework

We obtain the leading-twist valence quark transverse-momentum-dependent parton distribution functions (TMD PDFs) for the proton within the basis light-front quantization (BLFQ) framework. Our results are consistent with lattice QCD calculations and our previous results for the collinear limit. We also obtain consistency with the Soffer-type bounds. Within our approach, we find that six T-even TMDs in the leading twist are all independent of each other, and previously found model-dependent relations do not hold. This is a promising sign that our results are representative of future, more extensive treatments of QCD. Furthermore, we obtain a non-trivial x-dependence of the <(p ⊥ ) 2 > and some consistency with the Gaussian ansatz but only in the small <(p ⊥ ) 2 > region. Those features suggest our results may be a useful alternative in future experimental extractions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

New determination of the branching ratio of the structure dependent radiative K + → e + ν e γ

The branching ratio of the structure dependent (SD) radiative K + → e + ν e γ decay relative to that of the K + → e + ν e (γ) decay including the internal bremsstrahlung (IB) process (K e2 (γ) ) has been measured in the J-PARC E36 experiment using plastic scintillator/lead sandwich detectors, in contrast to the previous E36 measurement, which used a CsI(Tl) calorimeter. In the analysis, the effect of IB was also taken into account in the SD radiative decay as $K$$^{SD}_{e2γ(γ)}$. By combining the new data with the previous E36 result after revision for the IB correction for $K$$^{SD}_{e2γ(γ)}$, a new value Br($K$$^{SD}_{e2γ(γ)}$)/Br(Ke 2(γ) ) = 1.20 ± 0.07 has e2γ(γ) e2γ(γ) been determined. This is consistent with a recent lattice QCD calculation, but larger than the expectation of Chiral Perturbation Theory (ChPT) at order O(p 4 ) and the previous KLOE value. Using the method to relate form factor and branching ratio described in the KLOE paper, the present result is also consistent with the form factor prediction based on a gauged nonlocal chiral quark model, but larger than that from ChPT at order O(p 6 ).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Charm degrees of freedom in hot matter from lattice QCD

We study the nature of charm degrees of freedom in hot strong interaction matter by performing lattice QCD calculations of the second and fourth-order cumulants of charm fluctuations, and their correlations with net baryon number, electric charge and strangeness fluctuations. We show that below the chiral crossover temperature thermodynamics of charm can be very well understood in terms of charmed hadrons. Above the chiral transition charm quarks show up as new degrees of freedom contributing to the partial charm pressure. However, up to temperatures as high as 175 MeV charmed hadron-like excitations provide a significant contribution to the partial charm pressure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Transverse single-spin asymmetries in $γ$ SIDIS as a direct probe of quark-gluon-quark longitudinal momentum structure

Transverse single-spin asymmetries in the semi-inclusive deep-inelastic production of isolated photons ($\gamma$SIDIS), $A_{UT}^{\gamma {\rm SIDIS}}$, provide an unprecedented opportunity to extract the quark-gluon-quark correlators $F_{FT}(x,x')$ and $G_{FT}(x,x')$ point-by-point in their full support $x,x'$. We utilize realistic models for these functions, based on input from the Sivers transverse momentum dependent parton distribution function and imposing constraints from the $d_2$ matrix element calculated in lattice QCD, in order to provide numerical estimates for $A_{UT}^{\gamma {\rm SIDIS}}$ at the Electron-Ion Collider (EIC). We thoroughly explore the EIC phase space in order to isolate in which regions the asymmetry can be sizable, finding it can be as much as $10\%$ or larger for certain kinematics. Given that $F_{FT}(x,x')$ and $G_{FT}(x,x')$ are basically unknown, $A_{UT}^{\gamma {\rm SIDIS}}$ will be an important future measurement to learn about multi-parton correlations in the nucleon.

Electron-Ion Collider↗