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

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Scale effects on core design, fuel costs, and spent fuel volume of pressurized water reactors

The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.

42 ENGINEERING↗
DOE OSTI2025

Window Observables for Benchmarking Parton Distribution Functions

Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel “window observables” that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken-𝑥, and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different window observables that can be defined within a region of 𝑥 where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.

lattice QCD↗
DOE OSTI2025

Window observables for benchmarking parton distribution functions

Global analysis of collider and fixed-target experimental data and calculations from lattice quantum chromodynamics (QCD) are used to gain complementary information on the structure of hadrons. We propose novel ``window observables'' that allow for higher precision cross-validation between the different approaches, a critical step for studies that wish to combine the datasets. Global analyses are limited by the kinematic regions accessible to experiment, particularly in a range of Bjorken-x, and lattice QCD calculations also have limitations requiring extrapolations to obtain the parton distributions. We provide two different ``window observables'' that can be defined within a region of x where extrapolations and interpolations in global analyses remain reliable and where lattice QCD results retain sensitivity and precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Unveiling the strong interaction among hadrons at the LHC

One of the key challenges for nuclear physics today is to understand from first principles the effective interaction between hadrons with different quark content. First successes have been achieved using techniques that solve the dynamics of quarks and gluons on discrete space-time lattices. Experimentally, the dynamics of the strong interaction have been studied by scattering hadrons off each other. Such scattering experiments are difficult or impossible for unstable hadrons and so high-quality measurements exist only for hadrons containing up and down quarks. Here we demonstrate that measuring correlations in the momentum space between hadron pairs produced in ultrarelativistic proton–proton collisions at the CERN Large Hadron Collider (LHC) provides a precise method with which to obtain the missing information on the interaction dynamics between any pair of unstable hadrons. Specifically, we discuss the case of the interaction of baryons containing strange quarks (hyperons). We demonstrate how, using precision measurements of proton–omega baryon correlations, the effect of the strong interaction for this hadron–hadron pair can be studied with precision similar to, and compared with, predictions from lattice calculations. The large number of hyperons identified in proton–proton collisions at the LHC, together with accurate modelling of the small (approximately one femtometre) inter-particle distance and exact predictions for the correlation functions, enables a detailed determination of the short-range part of the nucleon-hyperon interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Structure and good piezoelectric performance in the complex system of Pb[(Zn,Ni)Nb]O 3 –Pb[(In,Yb)Nb]O 3 –Pb(Zr,Hf,Ti)O 3

High-performance piezoelectrics are always demanded for the high-end application. In this study, a complex piezoelectric system of 0.49Pb(Zn 1/2 Ni 1/2 ) 1/3 Nb 2/3 O 3 – x Pb(In 1/2 Yb 1/2 ) 1/2 Nb 1/2 O 3 –(0.51 – x )Pb(Zr 1/2 Hf 1/2 ) 0.1 Ti 0.9 O 3 (0.16 ≤ x ≤ 0.23) was fabricated through the solid-state method. The structure, ferroelectric, piezoelectric, and dielectric properties were investigated. The optimum piezoelectric coefficient d 33 of 761 pC/N, high Curie temperature of 169 °C, dielectric permittivity ( ε r ) of 4557, and electromechanical coupling coefficient ( k p ) of 63% were found at the morphotropic phase boundary composition of x = 0.19, which are superior to other complex piezoelectric materials. In particular, a significant large-signal d 33 * of 913 pm/V and low strain hysteresis (6%) was obtained in the temperature range of 20–170 °C. Temperature-dependent x-ray diffraction (XRD) has demonstrated that good temperature stability is put down to the structure stability. The agreement between the calculated lattice strain from in situ high-energy synchrotron XRD data and the macroscopic measurements suggests that the large lattice strain has a dominant contribution to the high piezoelectric response. The high piezoelectric performance and good temperature stability makes it potential for application.

36 MATERIALS SCIENCE↗
DOE OSTI2026

PWR Core Analysis for Cycle Extension and Uprates with LEU+ Accident Tolerant Fuel and 80 GWd/Tonne Burnup Limit

The U.S. Nuclear Regulatory Commission has recently drafted a rule enabling fuel burnup increase in light water reactors up to 80 GWd/t. In conjunction with use of fuel enrichment up to 10%, and accident tolerant fuel (ATF), this is anticipated to facilitate 24-month cycles in PWRs, along with further power uprates. In this paper, PWR core analysis is performed for 20% increased PWR power output along with cycle extension up to 24 months, in combination with use of chromia-doped fuel and chromium-coated clad, considered to be the most near-term ATF concepts. In combination, these lead to challenging conditions with a core average discharge burnup of up to ~74 GWd/t, challenging even the 80 GWd/t burnup limit. Analysis is performed using the 2-step method with POLARIS (within SCALE) used for the lattice calculations and PARCS for the core calculations. Core designs are first baselined for current operating conditions (LEU, 62 GWd/t discharge burnup limit) and then derived that meet cycle constraints on power distribution and the updated lead pin discharge burnup limit while maintaining at least two batches of fuel in the core. Gadolina loadings in fuel pins of up to 8% are used, with enrichment zoning both within the core and, to a limited extent, within assemblies. Here, doped fuel with coated cladding can utilize the same core designs as the reference UOX cores, exhibiting slightly lower burnup due to higher fuel density, which also offsets the slight reactivity penalty from the doping and coating. For the analysis performed here, doped fuel enabled a core with 24-month cycle and 20% uprate to stay within the 80 GWd/t lead pin discharge burnup limit.

LEU+↗
DOE OSTI2020

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dynamic properties of a polaron coupled to dispersive optical phonons

We study static and dynamic properties of an electron coupled to dispersive quantum optical phonons in the framework of the Holstein model defined on a one-dimensional lattice. Calculations are performed using the Lanczos algorithm based on a highly efficient construction of the variational Hilbert space. Even small phonon dispersion has a profound effect on the low-energy optical response. While the upward phonon dispersion broadens the optical spectra due to single-phonon excitations, the downward dispersion has the opposite effect. With increasing dispersion, a multiphonon excitation (MPE) state becomes the lowest excited state of the system at zero momentum and determines the low-frequency response of the optical conductivity where the threshold for optical absorption moves below the single-phonon frequency. Multiphonon states form a well-defined bandlike feature just above the polaron band as clearly seen in the electron spectral function. Low-energy MPEs should be observable in systems with strong optical phonon dispersion in optical as well as angle-resolved photoemission experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Joint lattice QCD–dispersion theory analysis confirms the quark-mixing top-row unitarity deficit

Recently, the first ever lattice computation of the $γ\textit{W}$-box radiative correction to the rate of the semileptonic pion decay allowed for a reduction of the theory uncertainty of that rate by a factor of ~3. A recent dispersion evaluation of the $γ\textit{W}$-box correction on the neutron also led to a significant reduction of the theory uncertainty, but shifted the value of $V_{ud}$ extracted from the neutron and superallowed nuclear β decay, resulting in a deficit of the Cabibbo-Kobayashi-Maskawa (CKM) unitarity in the top row. A direct lattice computation of the $γ\textit{W}$-box correction for the neutron decay would provide an independent cross-check for this result but is very challenging. Before those challenges are overcome, we propose a hybrid analysis, converting the lattice calculation on the pion to that on the neutron by a combination of dispersion theory and phenomenological input. The new prediction for the universal radiative correction to free and bound neutron β-decay reads $Δ^V_R$ = 0.02477(24), in excellent agreement with the dispersion theory result $Δ^V_R$ = 0.02467(22). Combining with other relevant information, the top-row CKM unitarity deficit persists.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Renormalons and power corrections in pseudo- and quasi-GPDs

High-order behavior of the perturbative expansion for short-distance observables in QCD is intimately related to the contributions of small momenta in the corresponding Feynman diagrams and this correspondence provides one with a useful tool to investigate power-suppressed nonperturbative corrections. We use this technique to study the structure of power corrections to parton quasi- and pseudo-GPDs which are used in lattice calculations of generalized parton distributions. As the main result, we predict the functional dependence of the leading power corrections to quasi(pseudo)-GPDs on x variable for nonzero skewedness parameter ξ . The kinematic point x = ± ξ turns out to be special. We find that the nonperturbative corrections to quasi-GPDs at this point are suppressed by the first power of the hard scale only. These contributions come from soft momenta and have nothing to do with the known UV renormalon in the Wilson line. We also show that power corrections can be strongly suppressed by the normalization procedure. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Stealth dark matter spectrum using Laplacian Heaviside smearing and irreducible representations

We present nonperturbative lattice calculations in the quenched approximation of the low-lying meson and baryon spectrum of the SU(4) gauge theory with fundamental fermion constituents. This theory is one instance of stealth dark matter, a class of strongly coupled theories, where the lowest mass stable baryon is the dark matter candidate. This work constitutes the first milestone in the program to study stealth dark matter self-interactions. Here, we focus on reducing excited state contamination in the single-baryon channel by applying the Laplacian Heaviside method, as well as projecting our baryon operators onto the irreducible representations of the octahedral group. We compare our resulting spectrum to previous work involving Gaussian smeared nonprojected operators and find good agreement with reduced statistical uncertainties. We also present the spectrum of the low-lying odd-parity baryons for the first time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

First-Principles Calculation of Electroweak Box Diagrams from Lattice QCD

We present the first realistic lattice QCD calculation of the $γ W$-box diagrams relevant for beta decays. The nonperturbative low-momentum integral of the $γ W$ loop is calculated using a lattice QCD simulation, complemented by the perturbative QCD result at high momenta. Using the pion semileptonic decay as an example, we demonstrate the feasibility of the method. By using domain wall fermions at the physical pion mass with multiple lattice spacings and volumes, we obtain the axial γ W -box correction to the semileptonic pion decay, $\square$$\frac{VA}{γW}|_π$ = 2.830 (11) stat (26 )syst × 10 –3 , with the total uncertainty controlled at the level of ~1% . This study sheds light on the first-principles computation of the $γ W$-box correction to the neutron decay, which plays a decisive role in the determination of $|V_ud|$ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Upper Limit on the Photoproduction Cross Section of the Spin-Exotic $π$ 1 (1600)

Here, the spin-exotic hybrid meson π 1 (1600) is predicted to have a large decay rate to the ωππ final state. Using 76.6 pb -1 of data collected with the GlueX detector, we measure the cross sections for the reactions γp → ωπ + π - p, γp → ωπ 0 π 0 p, and γp → ωπ - π 0 Δ ++ in the range E γ = 8 - 10 GeV. Using isospin conservation, we set the first upper limits on the photoproduction cross sections of the π$^0_1$(1600) and $π^-_1$(1600). We combine these limits with lattice calculations of decay widths and find that photoproduction of η' π is the most sensitive two-body system to search for the π 1 (1600).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pion and Kaon decays from lattice QCD and |Vud|, |Vus|

We review lattice calculations of the form factors for Kaon semileptonic decays and the leptonic decay constants of Kaons and Pions. We discuss the determinations of the CKM matrix elements $|V_{ud}|$ and $|V_{us}|$ using the form factors and experimental results. We use the CKM matrix elements to test first-row unitarity of the CKM matrix.

Simone, J. N. [Fermilab] (ORCID:0000000185153337)↗