Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Parity-violating electric-dipole transitions in helium

The paper examines parity-violating electric-dipole transitions in He in order to gain insight into the reliability of approximate calculations which are carried out for transitions in many-electron atoms. The contributions of the nearest-lying states are computed with a variety of wave functions, including very simple product wave functions, Hartree-Fock functions and Hylleraas-type wave functions with up to 84 parameters. It is found that values of the matrix elements of the parity-violating interaction can differ considerably from the values obtained from the good wave functions, even when these simple wave functions give accurate values for the matrix elements in question

Hiller, J.↗
DOE OSTI2026

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
DOE OSTI2020

Nucleon pair shell model in $\textit{M}$ scheme

The nucleon pair shell model (NPSM) is cast into the so-called $\textit{M}$ scheme for the cases with isospin symmetry and without isospin symmetry. The odd system and even system are treated on the same footing. The uncoupled commutators for nucleon pairs, which are suitable for the $\textit{M}$ scheme, are given. Explicit formula of matrix elements in $\textit{M}$ scheme for overlap, one-body operators, and two-body operators are obtained. In this work, it is found that the CPU time used in calculating the matrix elements in $\textit{M}$ scheme is much shorter than that in the $\textit{J}$ scheme of NPSM.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1992

Derivation of phase statistics from the Mueller matrix

Experimental observations show a strong dependence of phase differences of scattering matrix elements on the physical parameters of random media. Here, the statistical behavior of the phase differences is studied for distributed targets. The pdfs of the phase differences are derived from the Mueller matrix of the target. In deriving the density functions, it is assumed that the real and imaginary parts of the copolarized and cross-polarized terms of the scattering matrix are jointly Gaussian and their covariance matrices are found in terms of the Mueller matrix elements. The functional forms of the copolarized and cross-polarized density functions are similar and are obtained independently. It is shown that the density function of the phase difference is completely determined in terms of only two parameters.

Sarabandi, K.↗
DOE OSTI2025

Three-dimensional imaging of pion using lattice QCD: generalized parton distributions

In this work, we report a lattice calculation of x-dependent valence pion generalized parton distributions (GPDs) at zero skewness with multiple values of the momentum transfer −t. The calculations are based on an N f = 2 + 1 gauge ensemble of highly improved staggered quarks with Wilson-Clover valence fermion. The lattice spacing is 0.04 fm, and the pion valence mass is tuned to be 300 MeV. We determine the Lorentz-invariant amplitudes of the quasi-GPD matrix elements for both symmetric and asymmetric momenta transfers with similar values and show the equivalence of both frames. Then, focusing on the asymmetric frame, we utilize a hybrid scheme to renormalize the quasi-GPD matrix elements obtained from the lattice calculations. After the Fourier transforms, the quasi-GPDs are then matched to the light-cone GPDs within the framework of large momentum effective theory with improved matching, including the next-to-next-to-leading order perturbative corrections, and leading renormalon and renormalization group resummations. We also present the 3-dimensional image of the pion in impact-parameter space through the Fourier transform of the momentum transfer −t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Radiative processes on a quantum computer

Radiative processes, where a photon/neutrino is emitted as a result of a collision or decay of a particle, play a central role in atomic, nuclear, and particle physics. Their rate is determined by certain off-diagonal matrix elements between different initial and final states. We propose a method to compute them using quantum computers. It relies on a single extra qubit that, in a certain sense, represents the photon/neutrino. The generic formula relating this matrix element to the amplitude and frequency of oscillations of the extra qubit is derived for the near-resonance case. Furthermore, we demonstrate the feasibility of the method by using it in actual quantum computations and simulations of simple systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Finite element solution for energy conservation using a highly stable explicit integration algorithm

Theoretical derivation of a finite element solution algorithm for the transient energy conservation equation in multidimensional, stationary multi-media continua with irregular solution domain closure is considered. The complete finite element matrix forms for arbitrarily irregular discretizations are established, using natural coordinate function representations. The algorithm is embodied into a user-oriented computer program (COMOC) which obtains transient temperature distributions at the node points of the finite element discretization using a highly stable explicit integration procedure with automatic error control features. The finite element algorithm is shown to posses convergence with discretization for a transient sample problem. The condensed form for the specific heat element matrix is shown to be preferable to the consistent form. Computed results for diverse problems illustrate the versatility of COMOC, and easily prepared output subroutines are shown to allow quick engineering assessment of solution behavior.

Baker, A. J.↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Beyond leading twist: 𝜌 meson decay constants and distribution amplitudes in a self-consistent light-front quark model

In this study, we present a comprehensive analysis of decay constants and chiral-even and chiral-odd distribution amplitudes (DAs) up to twist 4 for the 𝜌 meson in the standard light-front quark model (LFQM) based on the Bakamjian-Thomas (BT) construction. For the 𝜌 meson, which possesses both longitudinal (ℎ=0) and transverse (ℎ =±1) polarizations, two types of decay constants, 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$, arises accordingly. We demonstrate that these decay constants can be self-consistently extracted from both local (𝑧𝜇 =0) and nonlocal (𝑧𝜇 ≠0) matrix elements ⟨0⁢|$\overline{𝑞}$(𝑧)⁢Γ⁢𝑞⁡(−𝑧)|⁢𝜌⁡(𝑃,ℎ)⟩, with Γ=(𝛾 𝜇 ,𝜎 𝜇⁢𝜈 ,𝛾 𝜇 ⁢𝛾 5 ,𝟏), in a manner independent of current components, polarizations, and reference frames. In particular, we emphasize the role of nonlocal matrix elements involving axial-vector and scalar currents, where mixing between 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$ occurs. We show that this mixing is consistently resolved through the BT construction, ensuring the proper extraction of these decay constants.

Electroweak interaction↗
DOE OSTI2024

Consistency of the pion form factor and unpolarized transverse momentum dependent parton distributions beyond leading twist in the light-front quark model

We investigate the interplay among the pion’s form factor, transverse momentum dependent distributions (TMDs), and parton distribution functions (PDFs) extending our light-front quark model (LFQM) computation based on the Bakamjian-Thomas construction for the two-point function [1,2] to the three-point and four-point functions. Ensuring the four-momentum conservation at the meson-quark vertex from the Bakamjian-Thomas construction, the meson mass is taken consistently as the corresponding invariant meson mass both in the matrix element and the Lorentz factor in our LFQM computation. We achieve the current-component independence in the physical observables such as the pion form factor and delve into the derivation of unpolarized TMDs and PDFs associated with the forward matrix element. We address the challenges posed by twist-4 TMDs and exhibit the fulfillment of the sum rule. Effectively, our LFQM successfully handles the light-front zero modes and offers insights for broader three-point and four-point functions and related observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Timelike meson form factors beyond the elastic region from lattice QCD

We present a calculation of the vector-isovector timelike form factors of the pion and the kaon using lattice quantum chromodynamics. We calculate two-point correlation functions with m π ∼ 280 MeV , extracting both the finite-volume spectrum and matrix elements for these states created from the vacuum by a vector current. After determining the coupled-channel π π , K K ¯ scattering amplitudes, we perform the necessary correction for the significant finite-volume effects present in the current matrix elements, leading to the timelike form factors. We find these to be dominated by the presence of the ρ resonance, and we extract its decay constant by an analytic continuation of the amplitudes to the resonance pole. In addition, the spacelike pion form factor is determined on the same lattice configurations, and a dispersive parametrization is used to simultaneously describe the spacelike and elastic timelike regions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of 48 Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. In this work, we use our scheme to compute the matrix element that governs the neutrinoless double beta decay of 48 Ca to 48 Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76 Ge, 130 Te, and 136 Xe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2023

Spectroscopic quadrupole moments in 124 $\mathrm{Xe}$

The Xe isotopic chain with four valence protons above the Z=50 shell closure is an ideal laboratory for the study of the evolution of nuclear deformation. At the N=82 shell closure, 136 Xe presents all characteristics of a doubly closed shell nucleus with a spherical shape. In the very neutron-deficient isotopes close to N=50, the α-decay chain of Xe was investigated to probe the radioactive decay properties near the drip-line and the magicity of 100 Sn. Additionally, the Xe isotopes present higher order symmetries in the nuclear deformation such as the octupole degree of freedom near N=60 and N=90 or O(6) symmetry in stable isotopes. The relevance of the O(6) symmetry has been investigated by measuring the spectroscopic quadrupole moment of the first excited states in 124 Xe. In the O(6) symmetry limit, the spectroscopic quadrupole moment of collective states is expected to be null. A stable 124 Xe beam with energies of 4.03A MeV and 4.11 A MeV was used to bombard a nat W target at the GANIL facility. Excited states were populated via the safe Coulomb excitation reaction. The collision of the heavy ions with a large Z at low energy make this reaction sensitive to the diagonal E2 matrix element of the excited states. The recoils were detected in the VAMOS++ magnetic spectrometer and the γ rays in the AGATA tracking array. The least squares fitting code gosia was used for the analysis to extract both E2 and M1 transitional and E2 diagonal matrix elements. The rotational ground state band was populated up to the 8$^{+}_{1}$ state as well as the 2$^{+}_{2}$ and 4$^{+}_{2}$ states. Using high precision spectroscopic data to constrain the gosia fit, the spectroscopic quadrupole moments of the 2$^{+}_{1}$, 4$^{+}_{1}$, and 6$^{+}_{1}$ states were determined for the first time. The spectroscopic quadrupole moments were found to be negative, large, and constant in the ground state band underlining the prolate axially deformed ground state band of 124 Xe. Here, the present experimental data confirm that the O(6) symmetry is substantially broken in 124 Xe .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,2 1 A and 1 3 A states of NH 3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

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

Exploration of mass splitting and muon/tau mixing parameters for an eV-scale sterile neutrino with IceCube

We present the first three-parameter fit to a 3+1 sterile neutrino model using 7.634 years of data from the IceCube Neutrino Observatory on $ν$ μ + $\overline{ν}$ μ charged-current interactions in the energy range 500–9976 GeV. Our analysis is sensitive to the mass-squared splitting between the heaviest and lightest mass state (Δm$^{2}_{41}$), the mixing matrix element connecting muon flavor to the fourth mass state (|U μ4 | 2 ), and the element connecting tau flavor to the fourth mass state (|U $\tau$4 | 2 ). Predicted propagation effects in matter enhance the signature through a resonance as atmospheric neutrinos from the Northern Hemisphere traverse the Earth to the IceCube detector at the South Pole. The remaining sterile neutrino matrix elements are left fixed, with |U e4 | 2 = 0 and δ 14 = 0, as they have a negligible effect, and δ 24 = π is set to give the most conservative limits. The result is consistent with the no-sterile neutrino hypothesis with a probability of 4.3%. Profiling the likelihood of each parameter yields the 90% confidence levels: 2.4eV 2 < Δm$^{2}_{41}$) < 9.6eV 2 , 0.0081< |U μ4 | 2 < 0.10, and |U $\tau$4 | 2 < 0.035, which narrows the allowed parameter-space for |U $\tau$4 | 2 . However, the primary result of this analysis is the first map of the 3+1 parameter space exploring the interdependence of Δm$^{2}_{41}$, |U μ4 | 2 , and |U $\tau$4 | 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Non-singlet quark helicity $\mathrm{PDF}$s of the nucleon from pseudo-distributions

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional z 2 contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional z 2 contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale μ 2 . The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned z 2 contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics

We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.

First-principles calculations↗