Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix elements”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Three-pion effects in $K^0-\bar{K}^0$ mixing

The rate of mixing between a neutral kaon and an anti-kaon ($K^0$-$\bar{K}^0$) is given, in part, by a long-range matrix element, defined with two insertions of the weak Hamiltonian separated by physical, Minkowski time evolution. For physical quark masses, the kaon mass lies above the two- and three-pion thresholds and, as a result, this long-range matrix element receives contributions from intermediate on-shell $2\pi$ and $3\pi$ states. These contributions cannot easily be captured in a finite Euclidean spacetime, meaning that such matrix elements are not directly accessible via lattice QCD. In this talk, we present a strategy for combining quantities that can be extracted in numerical lattice QCD calculations in order to reproduce the physical, infinite-volume long-range amplitude for $K^0-\bar{K}^0$. The key novelty relative to published work is that we fully include the effects of three-particle states that were previously neglected. The strategy is built on existing formalism for long-range matrix elements with two-particle intermediate states, together with the relativistic-field-theory finite-volume formalism for extracting three-hadron weak decays.

Jackura, Andrew↗
DOE OSTI2025

Toward shell model interactions with credible uncertainties

Background: The nuclear shell model is a powerful framework for predicting nuclear structure observables, but relies on interaction matrix elements fit to experimental data as its inputs. Extending the shell model's applicability, particularly toward dripline nuclei, requires efficient fitting methods and credible uncertainty quantification. Traditional approaches face computational challenges and may underestimate uncertainties. Purpose: We develop and test a framework combining eigenvector continuation and Markov chain Monte Carlo to efficiently fit shell model interaction matrix elements and quantify their uncertainties. Methods: Eigenvector continuation is used to emulate shell model calculations, reducing computational costs. The emulator enables Markov chain Monte Carlo sampling to optimize interaction matrix elements and rigorously assess parametric uncertainties. Here, the framework is benchmarked using the USDB interaction in the 𝑠⁢𝑑 shell. Results: The emulator reproduces the USDB interaction with negligible error, validating its use in shell model fitting applications. However, we find that to obtain credible predictive intervals, the model defect of the shell model itself, rather than experimental or emulator error, must be taken into account in order to obtain credible uncertainties. Conclusions: The proposed framework provides an efficient and rigorous approach for fitting shell model interactions and quantifying uncertainties. Further, the normality assumption used in the past appears sufficient to describe the distribution of interaction matrix elements. However, it is crucial to account for model correlations to avoid underestimating uncertainties.

Nuclear forces↗
DOE OSTI2023

The forward-backward asymmetry and differences of partial moments in inclusive semileptonic $B$ decays

Global fits to moments of kinematic distributions measured in inclusive semileptonic B\rightarrow X_c l \nu_l B → X c l ν l enable the determination of the Cabibbo-Kobayashi-Maskawa matrix element \left|V_{cb}\right| | V c b | together with non-perturbative matrix elements of the heavy quark expansion. In current fits, only two distinct kinematic distributions are employed and, as a consequence, higher moments of these distributions need to be taken into account to extract the relevant non-perturbative matrix elements. The moments of a given distribution are highly correlated and experimental uncertainties increase for higher moments. To address these issues, Turczyk suggested the inclusion of the charged lepton forward-backward asymmetry \mathcal{A}_{FB} 𝒜 F B in global fits, since it provides information on non-perturbative parameters beyond the commonly used moments. It is possible to construct differences of partial moments of kinematic distributions, which can provide additional information on the non-perturbative parameters beyond \mathcal{A}_{FB} 𝒜 F B and are studied in this work for the first time. Further, experimental cuts on the four-momentum transfer square are studied and are shown to preserve the shape of the angular distribution, in contrast to commonly used cuts on the lepton energy. Finally, the impact of final-state radiation and experimental lepton identification requirements on measurements of \mathcal{A}_{FB} 𝒜 F B and differences of partial moments are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Moments of parton distribution functions of any order from lattice QCD

We describe a procedure to determine moments of parton distribution functions of any order in lattice quantum chromodynamics (QCD). The procedure is based on the gradient flow for fermion and gauge fields. The flowed matrix elements of twist-2 operators renormalize multiplicatively, and the matching with the physical matrix elements can be obtained using continuum symmetries and the irreducible representations of Euclidean 4-dimensional rotations. We calculate the matching coefficients at one-loop in perturbation theory for moments of any order in the flavor nonsinglet case. We also give specific examples of operators that could be used in lattice QCD computations. It turns out that it is possible to choose operators with identical Lorentz indices and still have a multiplicative matching. One can thus use twist-2 operators exclusively with temporal indices, thus substantially improving the signal-to-noise ratio in the computation of the hadronic matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Coupled-Cluster Calculations of Neutrinoless Double- β Decay in Ca 48

We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-$\beta$ decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-$\beta$ decay of $^{48}$Ca with a quenching factor deduced from two-body currents in recent ab initio calculation of the Ikeda sum rule in $^{48}$Ca [Gysbers et al., Nat. Phys. 15, 428 (2019)].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Generalization of rectangular element stiffness matrix and thermal load vector associated with a(0) + a(1)x + a(2)y + a(3)xy type interpolation rule

The stiffness-matrix formulation for the rectangular finite element described by Melosh (1963) and Weaver and Johnston (1984) is generalized to orthotropic materials with material axes not necessarily coincident with the x and y axes; i.e., the condition d(13) = d(23) = 0 is removed. Also included are explicit expressions for the element load vector associated with nonuniform temperature increase in the element. Applications to the analysis of thermal stresses in thin Si-crystal ribbons subjected to temperature changes with highly nonuniform lengthwise and transverse gradients (Utku et al., 1986) and to the simulation of the thermoviscoelastic behavior of growing Si ribbons (Utku and Ray, 1986) are indicated.

Ray, Sujit K.↗
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 OSTI2022

E2 rotational invariants of 0$^{+}_{1}$ and 2$^{+}_{1}$ states for 106 Cd: The emergence of collective rotation

The collective structure of 106 Cd is elucidated by multi-step Coulomb excitation of a 3.849 MeV/A beam of 106 Cd on a 1.1 mg/cm 2 208 Pb target using GRETINA-CHICO2 at ATLAS. Fourteen E2 matrix elements were obtained. The nucleus 106 Cd is a prime example of emergent collectivity that possesses a simple structure: it is free of complexity caused by shape coexistence and has a small, but collectively active number of valence nucleons. This work follows in a long and currently active quest to answer the fundamental question of the origin of nuclear collectivity and deformation, notably in the cadmium isotopes. The results are discussed in terms of phenomenological models, the shell model, and Kumar-Cline sums of E2 matrix elements. The < 0$^{+}_{2}$ ||E2||2$^{+}_{1}$ > matrix element is determined for the first time, providing a total, converged measure of the electric quadrupole strength, < Q 2 >, of the first-excited 2$^{+}_{1}$ level relative to the 0$^{+}_{1}$ ground state, which does not show an increase as expected of harmonic and anharmonic vibrations. Strong evidence for triaxial shapes in weakly collective nuclei is indicated; collective vibrations are excluded. This is contrary to the only other cadmium result of this kind in 114 Cd by C. Fahlander et al., Nucl. Phys. A485, 327 (1988), which is complicated by low-lying shape coexistence near midshell.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Quantum Monte Carlo calculation of 𝛿 NS in 10 C using an effective field theory approach

Here, we compute radiative corrections to the superallowed 𝛽 decay of 10 C in an effective field theory approach using nuclear matrix elements obtained from quantum Monte Carlo calculations. These corrections are an important ingredient in the extraction of the Cabibbo-Kobayashi-Masakawa quark mixing matrix element 𝑉 𝑢⁢𝑑 , and the role of this work is to illuminate the uncertainties arising from nuclear structure. Our results provide good agreement with both the traditional extraction of 𝑉 𝑢⁢𝑑 , as well as with a more recent evaluation performed using the no-core shell model and a dispersion formalism. The dominant uncertainty in this approach is the presence of two unknown low-energy constants that enter into the relevant nuclear matrix elements. Future determinations of these low-energy constants—either from QCD or modeling them with two nucleon amplitudes—would improve the precision of the extraction in this formalism.

beta decay↗
DOE OSTI2022

Developments in Performance and Portability for MadGraph5_aMC@NLO

Event generators simulate particle interactions using Monte Carlo techniques, providing the primary connection between experiment and theory in experimental high energy physics. These software packages, which are the first step in the simulation worflow of collider experiments, represent approximately 5 to 20% of the annual WLCG usage for the ATLAS and CMS experiments. With computing architectures becoming more heterogeneous, it is important to ensure that these key software frameworks can be run on future systems, large and small. In this contribution, recent progress on porting and speeding up the Madgraph5_aMC@NLO event generator on hybrid architectures, i.e. CPU with GPU accelerators, is discussed. The main focus of this work has been in the calculation of scattering amplitudes and "matrix elements", which is the computational bottleneck of an event generation application. For physics processes limited to QCD leading order, the code generation toolkit has been expanded to produce matrix element calculations using C++ vector instructions on CPUs and using CUDA for NVidia GPUs, as well as using Alpaka, Kokkos and SYCL for multiple CPU and GPU architectures. Performance is reported in terms of matrix element calculations per time on NVidia, Intel, and AMD devices. The status and outlook for the integration of this work into a production release usable by the LHC experiments, with the same functionalities and very similar user interfaces as the current Fortran version, is also described.

Valassi, Andrea↗
NASA NTRS1974

Evaluation of atomic constants for optical radiation, volume 2

Various atomic constant for 23 elements from helium to mercury were computed and are presented in tables. The data given for each element start with the element name, its atomic number, its ionic state, and the designation and series limit for each parent configuration. This is followed by information on the energy level, parent configuration, and designation for each term available to the program. The matrix elements subtables are ordered by the sequence numbers, which represent the initial and final levels of the transitions. Each subtable gives the following: configuration of the core or parent, designation and energy level for the reference state, effective principal quantum number, energy of the series limit, value of the matrix element for the reference state interacting with itself, and sum of all of the dipole matrix elements listed in the subtable. Dipole and quadrupole interaction data are also given.

Kylstra, C. D.↗
DOE OSTI2022

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

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

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

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Finite-element time evolution operator for the anharmonic oscillator

The finite-element approach to lattice field theory is both highly accurate (relative errors approximately 1/N(exp 2), where N is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly preserved at the lattice sites). In this talk I construct matrix elements for dynamical variables and for the time evolution operator for the anharmonic oscillator, for which the continuum Hamiltonian is H = p(exp 2)/2 + lambda q(exp 4)/4. Construction of such matrix elements does not require solving the implicit equations of motion. Low order approximations turn out to be extremely accurate. For example, the matrix element of the time evolution operator in the harmonic oscillator ground state gives a results for the anharmonic oscillator ground state energy accurate to better than 1 percent, while a two-state approximation reduces the error to less than 0.1 percent.

Milton, Kimball A.↗
DOE OSTI2025

Gluon unpolarized, polarized, and transversity GPDs from lattice QCD: Lorentz-covariant parametrization

We identify the matrix elements necessary to determine the leading-twist gluon generalized parton distributions (GPDs) H g , E g , H ˜ g , E ˜ g , H g T , E g T , H ˜ g T , E ˜ g T in lattice QCD calculations. We present a method to achieve a Lorentz-covariant parameterization of the matrix elements in terms of a linearly independent basis of tensor structures. This parameterization is crucial for projecting lattice QCD matrix elements onto light cone distributions. For the first time, we determine the corresponding components that project onto the linear combinations of invariant amplitudes, which reduce to the different gluon GPDs in the light cone limit and enable their separation in a lattice QCD calculation for spin-0 and spin- 1 2 hadrons. Hence, this work lays the foundation for the numerical determination of the gluon GPDs from first-principle lattice QCD calculations, directly advancing our understanding of the mass and spin structures and mechanical properties of the nucleon, as well as the physics underlying deeply virtual Compton scattering and deeply virtual meson production in a range of experimental processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Toward the first gluon parton distribution from the LaMET

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

Good, William (ORCID:0000000184281003)↗
NASA NTRS1994

Rovibrational intensities of the minor isotopes of the CO X (sup 1)Sigma (sup +) state for v less than or equal to 20 and J less than or equal to 150

Electric dipole transition matrix elements for rovibrational transitions in the X (sup 1)Sigma(sup +) state of the CO minor isotopes (14)C(16)O and (13)C(17)O are calculated for the first time for all the delta v = +1, +2, and +3 transitions for which v less than or equal to 20 and J less than or equal to 150. Improved electric dipole transition matrix elements are also calculated for the minor isotopes (12)C(17)O, (12)C(18)O, (13)C(18)O. We have fitted polynomials to these matrix elements as a function of the parameter m which is defined in terms of the lower state angular momentum quantum number J; the convenient to use polynomial representations are given in tabular form. These results for the minor species of CO complement those previously reported by us for (12)C(16)O and (13)C(16)O.

Goorvitch, D.↗