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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Collaboration 51 Correlation Function Analysis Suite (c51_corr_analysis) v0.1.0

This is a data analysis package designed for analyzing correlation functions generated with lattice QCD calculations. The purpose is to make a centralized software suite for use by all members of my collaboration, so that various members can spend less time developing their own analysis codes, and more time extracting the interesting physics from our calculations. There are also 3 independent collaborations that have expressed interest in using this code, so there is some expectation it will be used in the broader international lattice QCD community.

Walker-Loud, Andre↗
DOE OSTI2021

Few-nucleon matrix elements in pionless effective field theory in a finite volume

Here, pionless effective field theory in a finite volume (FVEFT $\notπ$ ) is investigated as a framework for the analysis of multinucleon spectra and matrix elements calculated in lattice QCD (LQCD). By combining FVEFT π with the stochastic variational method, the spectra of nuclei with atomic number A ∈ {2, 3} are matched to existing finite-volume LQCD calculations at heavier-than-physical quark masses corresponding to a pion mass m $\notπ$ = 806 MeV , thereby enabling infinite-volume binding energies to be determined using infinite-volume variational calculations. Based on the variational wave functions that are constructed in this approach, the finite-volume matrix elements of various local operators are computed in FVEFT $\notπ$ and matched to LQCD calculations of the corresponding QCD operators in the same volume, thereby determining the relevant one- and two-body effective field theory counterterms and enabling an extrapolation of the LQCD matrix elements to infinite volume. As examples, the scalar, tensor, and axial matrix elements are considered, in addition to the magnetic moments and the isovector longitudinal momentum fraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

CompF2: Theoretical Calculations and Simulation Topical Group Report

This report summarizes the work of the Computational Frontier topical group on theoretical calculations and simulation for Snowmass 2021. We discuss the challenges, potential solutions, and needs facing six diverse but related topical areas that span the subject of theoretical calculations and simulation in high energy physics (HEP): cosmic calculations, particle accelerator modeling, detector simulation, event generators, perturbative calculations, and lattice QCD (quantum chromodynamics). The challenges arise from the next generations of HEP experiments, which will include more complex instruments, provide larger data volumes, and perform more precise measurements. Calculations and simulations will need to keep up with these increased requirements. The other aspect of the challenge is the evolution of computing landscape away from general-purpose computing on CPUs and toward special-purpose accelerators and coprocessors such as GPUs and FPGAs. These newer devices can provide substantial improvements for certain categories of algorithms, at the expense of more specialized programming and memory and data access patterns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Nucleon Electric Polarizabilities and Nucleon-Pion Scattering at the Physical Pion Mass

We present a lattice QCD calculation of the nucleon electric polarizabilities at the physical pion mass. Our findings reveal the substantial contributions of the N π states to these polarizabilities. Without considering these contributions, the lattice results fall significantly below the experimental values, consistent with previous lattice studies. This observation has motivated us to compute both the parity-negative N π scattering up to a nucleon momentum of ∼ 0.5 GeV in the center-of-mass frame and corresponding N γ * → N π matrix elements using lattice QCD. Our results confirm that incorporating dynamic N π contributions is crucial for a reliable determination of the polarizabilities from lattice QCD. This methodology lays the groundwork for future lattice QCD investigations into various other polarizabilities. Published by the American Physical Society 2024

Wang, Xuan-He↗
DOE OSTI2021

Progress report on computing the disconnected QCD and the QCD plus QED hadronic contributions to the muon's anomalous magnetic moment

We report progress on calculating the contribution to the anomalous magnetic moment of the muon from the disconnected hadronic diagrams with light and strange quarks and the valence QED contribution to the connected diagrams. The lattice QCD calculations use the highly-improved staggered quark (HISQ) formulation. The gauge configurations were generated by the MILC Collaboration with four flavors of HISQ sea quarks with physical sea-quark masses.

Bazavov, Alexei↗
DOE OSTI2020

Understanding dynamic properties of materials using neutron spectroscopy and atomistic simulation

Recent developments in first-principles lattice dynamics and classical force field based molecular dynamics are revolutionising the field of neutron spectroscopy. Herein we present a short review of these methods, their critical role in the supporting of cutting-edge experiments, and how they are improved by matching experimental data. We begin with a brief overview of how lattice dynamics calculations can be compared to inelastic neutron scattering (INS) and molecular dynamics simulations to both INS and quasi-elastic neutron scattering (QENS). We then provide a series of exemplar applications where lattice dynamics and molecular dynamics have been used in conjunction with neutron spectroscopy to bring significant understanding to topical areas of materials science namely: (i) lattice dynamics and INS for the study of hybrid organic-inorganic perovskites (ii) lattice dynamics and INS for the study of flexible porous solids and (iii) molecular dynamics and QENS for probing molecular behaviour in zeolite catalysis. In all three cases, the understanding gained through the synergy of experiment and computation would have been significantly reduced using either in isolation. Finally, we consider the current state of the art, describing outstanding challenges and suggesting future directions in this exciting and fertile area of physical science.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chiral-odd GPDs in the bag model

A study of chiral-odd generalized parton distributions (GPDs) of the nucleon is presented in the bag model demonstrating that in this model all four chiral-odd GPDs are nonzero contrary to other claims in literature. The bag model results for the GPDs H T q ( x , ξ , t ) , E T q ( x , ξ , t ) , H ˜ T q ( x , ξ , t ) agree with other models within a typical quark model accuracy. We present one of the few quark model calculations where polynomiality is satisfied and the sum rule ∫ d x E ˜ T q ( x , ξ , t ) = 0 holds. We confront our results with predictions from the large- N c limit, and with lattice QCD calculations. We conclude that the bag model successfully catches the main features of chiral-odd GPDs. Published by the American Physical Society 2024

Tezgin, Kemal (ORCID:0000000194929512)↗
DOE OSTI2023

Phonon-Induced Localization of Excitons in Molecular Crystals from First Principles

The spatial extent of excitons in molecular systems underpins their photophysics and utility for optoelectronic applications. Phonons are reported to lead to both exciton localization and delocalization. However, a microscopic understanding of phonon-induced (de)localization is lacking, in particular, how localized states form, the role of specific vibrations, and the relative importance of quantum and thermal nuclear fluctuations. Here, we present a first-principles study of these phenomena in solid pentacene, a prototypical molecular crystal, capturing the formation of bound excitons, exciton-phonon coupling to all orders, and phonon anharmonicity, using density functional theory, the ab initio GW-Bethe-Salpeter equation approach, finite-difference, and path integral techniques. In this study, we find that for pentacene zero-point nuclear motion causes uniformly strong localization, with thermal motion providing additional localization only for Wannier-Mott-like excitons. Anharmonic effects drive temperature-dependent localization, and, while such effects prevent the emergence of highly delocalized excitons, we explore the conditions under which these might be realized

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

$B \to π$, $B_{(s)} \to D_{(s)}$ from 2+1+1 Flavor Lattice QCD

We present a lattice-QCD calculation of the hadronic form factors for $B$-meson semileptonic decays computed using the highly improved staggered quark action for both valence and sea quarks on the MILC collaborations 2+1+1-flavor ensembles with lattice spacing ranging from 0.09 fm to 0.03 fm, many with physical pion masses On our finest ensembles, we compute the form factors directly at the physical $b$-quark mass. We discuss the computational setup and analysis strategies for two- and three-point correlation functions. For $B_{(s)} \to D_{(s)}$ we present preliminary results of chiral-continuum fits for the scalar and vector form factors. The goal of this project is a percent-level determination of the scalar and vector form factors to enable high-precision determinations of $|V_{ub}|$ and $|V_{cb}|$. This work fits into a broader program of lattice-QCD studies of weak $B$-meson decays by the Fermilab Lattice and MILC Collaborations.

Cassar, Nicholas [Colorado State U.]↗
DOE OSTI2022

Computational insights into the lattice dynamics of Pu(IV) oxalates

Despite its use in PuO 2 production, the structure of anhydrous Pu(C 2 O 4 ) 2 is still not completely understood. Recently, two candidate structures for Pu(C 2 O 4 ) 2 were proposed via density functional theory (DFT), after which the first experimental optical vibrational spectra were reported. Here, in this study, we calculated the lattice dynamics of the candidate structures using DFT and found that the primary difference between them is the presence of a vibrational mode near 1380 cm -1 in one structure. The frequency and optical activity of this mode agree well with the published experimental results, providing strong support for this calculated structure as that of anhydrous Pu(C 2 O 4 ) 2 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Generalized parton distributions from lattice QCD with asymmetric momentum transfer: Axial-vector case

Recently, we made significant advancements in improving the computational efficiency of lattice QCD calculations for generalized parton distributions (GPDs). This progress was achieved by adopting calculations of matrix elements in asymmetric frames, deviating from the computationally-expensive symmetric frame typically used, and allowing freedom in the choice for the distribution of the momentum transfer between the initial and final states. A crucial aspect of this approach involves the adoption of a Lorentz covariant parametrization for the matrix elements, introducing Lorentz-invariant amplitudes. This approach also allows us to propose an alternative definition of quasi-GPDs, ensuring frame independence and potentially reduce power corrections in matching to light cone GPDs. In our previous work, we presented lattice QCD results for twist-2 unpolarized GPDs ( H and E ) of quarks obtained from calculations performed in asymmetric frames at zero skewness. Building upon this work, we now introduce a novel Lorentz covariant parametrization for the axial-vector matrix elements. We employ this parametrization to compute the axial-vector GPD H ˜ at zero skewness, using an N f = 2 + 1 + 1 ensemble of twisted mass fermions with clover improvement. The light-quark masses employed in our calculations correspond to a pion mass of approximately 260 MeV. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Thermal properties of the metallic delafossite PdCo⁢O 2 : A combined experimental and first-principles study

Metallic delafossite materials (e.g., PdCo⁢O 2 , PtCo⁢O 2 ), have attracted much recent attention due to record-high oxide conductivities, the origins of which remain unclear. Relatively little attention has been paid to their related thermal properties, however. Here, we address this via wide temperature range experimental studies of the crystal structure, thermal expansion, and specific heat of single-crystal PdCo⁢O 2 , combined with density-functional theory (DFT) calculations of the electronic and phononic densities of states, and thus thermal properties. PdCo⁢O 2 is shown to retain the $R\bar{3}m$ space group from 12 to 1000 K, exhibiting a- and c-axis thermal expansion in good quantitative agreement with DFT-based lattice dynamics calculations. The Co-O bond lengths additionally elucidate the stability of the low-spin state of the nominally Co 3+ ions, which is a notable difference between the edge-shared Co-O octahedra in PdCo⁢O 2 and the corner-shared octahedra in Co-based perovskites. Measurements of specific heat from 1.9 to 400 K provide accurate values for the Debye temperature and Sommerfeld coefficient, the phononic part being describable via a combined Debye-Einstein approach (accounting for high-frequency oxygen-related optical phonons), with excess intermediate-temperature specific heat due to a prominent low-energy peak in the phonon density of states. Most significantly, all electronic and phononic contributions to the specific heat are shown to be remarkably closely reproduced by DFT-based calculations, establishing quantitative understanding of key thermal properties of the metallic delafossite PdCo⁢O 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

QCD Predictions for Physical Multimeson Scattering Amplitudes

We use lattice QCD calculations of the finite-volume spectra of systems of two and three mesons to determine, for the first time, three-particle scattering amplitudes with physical quark masses. Our results are for combinations of 𝜋 + and 𝐾 + , at a lattice spacing 𝑎 = 0.063 fm, and in the isospin-symmetric limit. We also obtain accurate results for maximal-isospin two-meson amplitudes, with those for 𝜋 + ⁢𝐾 + and 2⁢𝐾 + being the first determinations at the physical point. Dense lattice spectra are obtained using the stochastic Laplacian-Heaviside method, and the analysis leading to scattering amplitudes is done using the relativistic finite-volume formalism. Results are compared to chiral perturbation theory and to phenomenological fits to experimental data, finding good agreement.

hadron-hadron interactions↗
DOE OSTI2023

Three relativistic neutrons in a finite volume

We generalize the relativistic field-theoretic (RFT) three-particle finite-volume formalism to systems of three identical, massive, spin-1/2 fermions, such as three neutrons. This allows, in principle, for the determination of the three-neutron interaction from the finite-volume spectrum of three-neutron states, which can be obtained from lattice QCD calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Probing for the Trace Estimation of a Permuted Matrix Inverse Corresponding to a Lattice Displacement

We report thatpProbing is a general technique that is used to reduce the variance of the Hutchinson stochastic estimator for the trace of the inverse of a large, sparse matrix A. The variance of the estimator is the sum of the squares of the off-diagonal elements of A -1 . Therefore, this technique computes probing vectors that when used in the estimator annihilate the largest off-diagonal elements. For matrices that display decay of the magnitude of |A$^{-1}_{ij}$| with the graph distance between nodes i and j, this is achieved through graph coloring of increasing powers A k . Equivalently, when a matrix stems from a lattice discretization, it is computationally beneficial to find a distance-k coloring of the lattice. Previously, a hierarchical coloring was proposed so that k can be increased at runtime as needed without discarding previous work. In this work, we study probing for the more general problem of computing the trace of a permutation of A -1 , say PA -1 . The motivation comes from lattice quantum chromodynamics (QCD), where we need to construct “disconnected diagrams” to extract flavor-separated generalized parton functions. In lattice QCD, where the matrix has a four-dimensional toroidal lattice structure, these nonlocal operators correspond to a PA -1 , where P is the permutation relating to some displacement $\vec{p}$ in one or more dimensions. We focus on a single dimension displacement (p), but our methods are general. We show that probing on A k or (PA) k does not annihilate the largest magnitude elements. To resolve this issue, our displacement-based probing works on PA k using a new coloring scheme that works directly on appropriately displaced neighborhoods on the lattice. We prove lower bounds on the number of colors needed and study the effect of this scheme on variance reduction, both theoretically and experimentally on a real-world lattice QCD calculation. We achieve orders of magnitude speedup over the unprobed or the naively probed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Using resonant energy X-ray diffraction to extract chemical order parameters in ternary semiconductors

II–IV–V 2 materials, ternary analogs to III–V materials, are emerging for their potential applications in devices such as LEDs and solar cells. Controlling cation ordering in II–IV–V 2 materials offers the potential to tune properties at nearly fixed compositions and lattice parameters. While tuning properties at a fixed lattice constant through ordering has the potential to be a powerful tool used in device fabrication, cation ordering also creates challenges with characterization and quantification of ordering. In this work, we investigate two different methods to quantify cation ordering in ZnGeP 2 thin films: a stretching parameter calculated from lattice constants , and an order parameter determined from the cation site occupancies ( S ). We use high resolution X-ray diffraction (HRXRD) to determine and resonant energy X-ray diffraction (REXD) to extract S . REXD is critical to distinguish between elements with similar Z -number ( e.g. Zn and Ge). We found that samples with a corresponding to the ordered chalcopyrite structure had only partially ordered S values. The optical absorption onset for these films occurred at lower energy than expected for fully ordered ZnGeP 2 , indicating that S is a more accurate descriptor of cation order than the stretching parameter. Since disorder is complex and can occur on many length scales, metrics for quantifying disorder should be chosen that most accurately reflect the physical properties of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗