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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

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↗

NASA NTRS1976

Effect of polarization on superfluidity in low density neutron matter

The singlet-state quasi-particle interaction in neutron matter is examined on the basis of results of a detailed evaluation of the Landau Fermi-liquid parameters for pure neutron effects, including polarization effects. This means that the interaction induced by exchange of density and spin-density excitations is taken into account. It is shown that polarization actually works to suppress the pairing matrix elements, owing to the spin dependence of the quasi-particle interaction and, ultimately, the balance of attraction, repulsion, and spin dependence in the fundamental two-neutron interaction. Since the isotropic energy gap and the condensation energy in low-density neutron-star matter are extremely sensitive functions of the pairing matrix elements, they will also be suppressed by the polarizability of the neutron medium.

Clark, J. W.↗

NASA NTRS2004

Examination of Surface Roughness on Light Scattering by Long Ice Columns by Use of a Two-Dimensional Finite-Difference Time-Domain Algorithm

Natural particles such as ice crystals in cirrus clouds generally are not pristine but have additional micro-roughness on their surfaces. A two-dimensional finite-difference time-domain (FDTD) program with a perfectly matched layer absorbing boundary condition is developed to calculate the effect of surface roughness on light scattering by long ice columns. When we use a spatial cell size of 1/120 incident wavelength for ice circular cylinders with size parameters of 6 and 24 at wavelengths of 0.55 and 10.8 mum, respectively, the errors in the FDTD results in the extinction, scattering, and absorption efficiencies are smaller than similar to 0.5%. The errors in the FDTD results in the asymmetry factor are smaller than similar to 0.05%. The errors in the FDTD results in the phase-matrix elements are smaller than similar to 5%. By adding a pseudorandom change as great as 10% of the radius of a cylinder, we calculate the scattering properties of randomly oriented rough-surfaced ice columns. We conclude that, although the effect of small surface roughness on light scattering is negligible, the scattering phase-matrix elements change significantly for particles with large surface roughness. The roughness on the particle surface can make the conventional phase function smooth. The most significant effect of the surface roughness is the decay of polarization of the scattered light.

Sun, W.↗

DOE OSTI2025

Systematic Study of the Self-Renormalized Nucleon Gluon PDF in Large-Momentum Effective Theory

We present a systematic study of the nucleon gluon parton distribution function (PDF) using the self-renormalized large-momentum effective theory (LaMET) approach in lattice QCD. This work extends previous gluon-PDF extractions by performing a detailed analysis of key systematic effects, including gauge-link smearing, lattice spacing, pion mass, and nucleon boost momentum. The self-renormalization framework mitigates ultraviolet divergences associated with Wilson-line self-energy and renormalon contributions by combining lattice matrix elements with perturbative short-distance information, thereby preserving the correct infrared structure. Calculations are performed on $N_f=2+1+1$ HISQ ensembles generated by the MILC Collaboration at three lattice spacings and two pion masses, with boosted nucleon states reaching momenta up to 2.2~GeV. We determine renormalization factors from zero-momentum matrix elements and apply hybrid renormalization to suppress discretization artifacts. After extrapolating large-separation behavior and performing Fourier transforms, we reconstruct quasi-PDFs and match them to lightcone PDFs using next-to-leading order Wilson coefficients. Our results demonstrate that smearing and lattice-spacing effects are under control, and pion-mass and lattice-spacing dependence is mild relative to the current $O(10^6)$ statistics; however, momentum dependence remains a significant source of uncertainty. Future work including even larger boost momenta will be essential to reduce systematics in lattice determinations of the gluon PDF and to advance toward precision QCD phenomenology at the LHC and the future Electron-Ion Collider.

FOS: Physical sciences↗

NASA NTRS1977

Variational calculations of vibrational properties of ozone

A variational method is used to obtain vibrational-rotational properties for ozone from an experimental quartic force field. Band positions, average structures, matrix elements for calculating infrared intensities, and effective rotational constants are reported for (O-16)3 and its O-18 isotopic species. Also, the degree to which the vibrational energies and properties are converged is investigated as a function of the basis set parameters and basis set size, and of the method of obtaining the Hamiltonian matrix elements. A convenient procedure for assigning the vibrational states is developed for variational vibrational wavefunction expansions.

Carney, G. D.↗

NASA NTRS1980

Atmospheric particulate analysis using angular light scattering

Using the light scattering matrix elements measured by a polar nephelometer, a procedure for estimating the characteristics of atmospheric particulates was developed. A theoretical library data set of scattering matrices derived from Mie theory was tabulated for a range of values of the size parameter and refractive index typical of atmospheric particles. Integration over the size parameter yielded the scattering matrix elements for a variety of hypothesized particulate size distributions. A least squares curve fitting technique was used to find a best fit from the library data for the experimental measurements. This was used as a first guess for a nonlinear iterative inversion of the size distributions. A real index of 1.50 and an imaginary index of -0.005 are representative of the smoothed inversion results for the near ground level atmospheric aerosol in Tucson.

Hansen, M. Z.↗

NASA NTRS1980

Atmospheric particulate analysis using angular light scattering

Using light scattering matrix elements measured by a polar nephelometer, a procedure for estimating the characteristics of atmospheric particulates was developed. A theoretical library data set of scattering matrices derived from Mie theory was tabulated for a range of values of the size parameter and refractive index typical of atmospheric particles. Integration over the size parameter yielded the scattering matrix elements for a variety of hypothesized particulate size distributions. A least squares curve fitting technique was used to find a best fit for the experimental measurements. This was used as a first guess for a nonlinear iterative inversion of the size distribution. A real index of 1.50 and an imaginary index of -0.005 are representative of the smoothed inversion results for the near ground level atmospheric aerosol in Tucson.

Hansen, M. Z.↗

DOE OSTI2024

Generalized parton distributions from the pseudodistribution approach on the lattice

Generalized parton distributions (GPDs) are key quantities for the description of a hadron’s three-dimensional structure. They are the current focus of all areas of hadronic physics—phenomenological, experimental and theoretical, including lattice QCD. Synergies between these areas are desirable and essential to achieve precise quantification and understanding of the structure of, particularly, nucleons, as the basic ingredients of matter. In this paper, we investigate, for the first time, the numerical implementation of the pseudodistribution approach for the extraction of zero-skewness GPDs for unpolarized quarks. Pseudodistributions are Euclidean parton correlators computable in lattice QCD that can be perturbatively matched to the light-cone parton distributions of interest. Although they are closely related to the quasidistributions and come from the same lattice-extracted matrix elements, they are, however, subject to different systematic effects. We use the data previously utilized for quasi-GPDs and extend it with other momentum transfers and nucleon boosts, in particular a higher one ( P 3 = 1.67 GeV ) with eightfold larger statistics than the largest one used for quasidistributions ( P 3 = 1.25 GeV ). We renormalize the matrix elements with a ratio scheme and match the resulting Ioffe time distributions to the light cone in coordinate space. The matched distributions are then used to reconstruct the x dependence with a fitting . We investigate some systematic effects related to this procedure, and we also compare the results with the ones obtained in the framework of quasi-GPDs. Our final results involve the invariant four-momentum transfer squared ( − t ) dependence of the flavor nonsinglet ( u − d ) H and E GPDs. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1994

Calculation of (12)C(16)O and (13)C(16)O X(1)Sigma(+) rovibrational intensities for v less than or equal to 20 and J less than or equal to 150

Improved electric dipole transition matrix elements for rovibrational transitions in the ground state X(1)Sigma(+) of (12)C(6)O and (13)C(16)O are calculated 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. 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. These convenient to use polynomial representations are given in Tables 1-4 for (12)C(16)O and in Tables 5-8 for (13)C(16)O. We observe that there is intensity enhancement due to vibration-rotation interaction for the P-branch transitions at the expense of the R-branch transitions for delta v = +1. This enhancement can be as large as 40% at the highest J. For the delta v = +2 and +3 transitions, the R-branch transitions are enhanced by as much as a factor of 2.75 and 10 at the highest J, respectively. The P-branch transitions exhibit only minor decreases. Comparisons with previous calculations show good agreement for the delta v = +1 transitions. The comparison for delta v = +2 and +3 transitions show differences as large as a factor of 5.

Goorvitch, D.↗

NASA NTRS1997

The Influence of Correlated Crustal Signals in Modelling the Main Geomagnetic Field

Algorithms used in geomagnetic main-field modelling have for the most part treated the noise in the field measurements as if it were white. A major component of the noise consists of the field due to magnetization in the crust and it has been realized for some time that such signals are highly correlated at satellite altitude. Hence approximation by white noise, while of undoubted utility, is of unknown validity. In this paper we study two plausible statistical models for the crustal magnetization, in which the magnetization is a realization of a stationary, isotropic, random process. At a typical satellite altitude the associated fields exhibit significant correlation over ranges as great as 15 deg. or more, which introduces off-diagonal elements into the covariance matrix, elements that have usually been neglected in modelling procedures. Dealing with a full covariance matrix for a large data set would present a formidable computational challenge, but fortunately most of the entries in the covariance matrix are so small that they can be replaced by zeros. The resultant matrix comprises only about 3 per cent non-zero entries and thus we can take advantage of efficient sparse matrix techniques to solve the numerical system. We construct several main-field models based on vertical-component data from a selected 5 deg. by 5 deg. data set derived from the Magsat mission. Models with and without off-diagonal terms are compared.

Rygaard-Hjalsted, C.↗

DOE OSTI2026

Developments in Performance and Portability of BlockGen

For more than a decade Monte Carlo event generators with the current matrix element algorithms have been used for generating hard scattering events on CPU platforms, with excellent flexibility and good efficiency. While the HL-LHC is approaching and precision requirements are becoming more demanding, many studies have been made to solve the bottleneck in the current Monte Carlo event generator tool chains. The novel BlockGen family of fast matrix element algorithms shown in this report, is one of the new developments that are more suitable for GPU acceleration. We report the development experience of porting BlockGen using Kokkos. Moreover, we discuss the performance of the Kokkos version in comparison with the dedicated GPU version in CUDA.

Bothmann, E. [Gottingen U.]↗

DOE OSTI2020

Long-range quadrupole electron-phonon interaction from first principles

Lattice vibrations in materials induce perturbations on the electron dynamics in the form of long-range (dipole and quadrupole) and short-range (octopole and higher) potentials. The dipole Fröhlich term can be included in current first-principles electron-phonon (e-ph) calculations and is present only in polar materials. The quadrupole e-ph interaction is present in both polar and nonpolar materials, but currently it cannot be computed from first principles. In this work, we show an approach to compute the quadrupole e-ph interaction and include it in ab initio calculations of e-ph matrix elements. The accuracy of the approach is demonstrated by comparing with direct density functional perturbation theory calculations. We apply our method to silicon as a case of a nonpolar semiconductor and tetragonal PbTiO 3 as a case of a polar piezoelectric material. In both materials we find that the quadrupole term strongly impacts the e-ph matrix elements. Analysis of e-ph interactions for different phonon modes reveals that the quadrupole term mainly affects optical modes in silicon and acoustic modes in PbTiO 3 , although the quadrupole term is needed for all modes to achieve quantitative accuracy. The effect of the quadrupole e-ph interaction on electron scattering processes and transport is shown to be important. Our approach enables accurate studies of e-ph interactions in broad classes of nonpolar, polar, and piezoelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗