A systematic method for calculating reduced matrix elements
Matrix elements systematic calculation method in terms of inelastic transition densities to evaluate nuclear properties
SEARCH · Engineering Papers
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.
Matrix elements systematic calculation method in terms of inelastic transition densities to evaluate nuclear properties
The elements of the scattering or Mueller Matrix for three polydisperse systems of irregular, randomly-oriented particles have been measured in absolute terms as a function of scattering angle for several visible wavelengths. The samples consisted of commercially available silicon dioxide particles that fit three distinct lognormal size distributions. The measured matrix elements were compared with the matrix elements calculated for spheres that had the same refractive index and fitted the same distribution of projected areas. Correlations and discrepancies between the two sets of matrix elements will be discussed.
Processes that violate baryon number, most notably proton decay and $n\bar{n}$ transitions, are promising probes of physics beyond the Standard Model (BSM) needed to understand the lack of antimatter in the Universe. To interpret current and forthcoming experimental limits, theory input from nuclear matrix elements to UV complete models enters. Thus, an interplay of experiment, effective field theory, lattice QCD, and BSM model building is required to develop strategies to accurately extract information from current and future data and maximize the impact and sensitivity of next-generation experiments. Here, we briefly summarize the main results and discussions from the workshop ‘INT-25-91W: Baryon Number Violation: From Nuclear Matrix Elements to BSM Physics,’ held at the Institute for Nuclear Theory, University of Washington, Seattle, WA, 13–17 January 2025.
We present a new measurement of the Cabibbo-Kobayashi-Maskawa matrix element
A graphical representation of angular momentum was used to evaluate relativistic matrix elements between antisymmetrized states of many particle configurations having any number of open shells. The antisymmetrized matrix element was expanded as a sum of semisymmetrized matrix elements. The diagram representing a semisymmetrized matrix element was composed of four diagram blocks; the bra block, the ket block, the spectator block, and the interaction block. The first three blocks indicate the couplings of the two interacting configurations while the last depends on the interaction and is the replaceable component. Interaction blocks for relativistic operators and commonly used potentials were summarized in ready to use forms. A simple step by step procedure was prescribed generally for calculating antisymmetrized matrix elements of one and two particle operators.
Recent work [M. L. Wagman, Lanczos, the transfer matrix, and the signal-to-noise problem, .] found that an analysis formalism based on the Lanczos algorithm allows energy levels to be extracted from Euclidean correlation functions with faster ground-state convergence than effective masses, convergent estimators for multiple states from a single correlator, and two-sided error bounds. After filtering out spurious eigenvalues and using outlier-robust estimators within a nested bootstrap framework, Lanczos estimators behave more like multistate fit results than effective masses—but without involving statistical fitting. We extend this formalism to the determination of matrix elements from three-point correlation functions and provide a physical picture of “spurious-state filtering” involving restriction to a Hermitian subspace. We demonstrate similar advantages for matrix elements as for spectroscopy through example applications to noiseless mock-data and (bare) forward matrix elements of the strange scalar current between both ground and excited states with the quantum numbers of the nucleon.
The calculation of off-diagonal matrix elements has various applications in fields such as nuclear physics and quantum chemistry. In this paper, we present a noisy intermediate scale quantum algorithm for estimating the diagonal and off-diagonal matrix elements of a generic observable in the energy eigenbasis of a given Hamiltonian without explicitly preparing its eigenstates. By means of numerical simulations we show that this approach finds many of the matrix elements for the one and two qubits cases. Specifically, while in the first case, one can initialize the ansatz parameters over a broad interval, in the latter the optimization landscape can significantly slow down the speed of convergence and one should therefore be careful to restrict the initialization to a smaller range of parameters.
Neutrinoless double beta decay nuclear matrix element (M0ν) for 136Xe was recently analyzed using a statistical approach (Phys. Rev. C 107, 045501 (2023)). In the analysis, three initial shell model effective Hamiltonians were randomly altered, and their results for 23 measured observables were used to infer credibility for the M0ν nuclear matrix element (NME) based on a Bayesian Model Averaging approach. In that analysis, a reasonable Gamow-Teller quenching factor of 0.7 was assumed for each starting effective Hamiltonian. Given that the result of the statistical analysis was sensible to this choice, we are here improving that analysis by assuming that the Gamow-Teller quenching factor is also randomly chosen within reasonabe limits for all three starting Hamiltonians. The outcomes are slightly higher expectation values and uncertainties for the M0ν NME.
Recursion relations for Coulomb matrix elements, with psi and psi prime wave functions describing bound state
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.
We examined the ratios of neutron and proton transition matrix elements, 𝑀 𝑛 /𝑀 𝑝 , for the 0$^{+}_{𝑔𝑠}$ → 2$^{+}_{1}$ transitions in 48 even-even stable nuclei with 𝑁 > 20 for which electromagnetic matrix elements were compiled by Pritychenko et al. and for which high-quality inelastic proton-scattering data were available. Several deformed rare-earth nuclei have (𝑀 𝑛 /𝑀 𝑝 )/(𝑁/𝑍) values significantly below 1.0, which is not consistent with a simple liquid-drop picture. However, this phenomenon can be explained using a schematic picture in which 𝑀 𝑝 reaches a maximum at proton midshell (𝑍 = 66) and 𝑀𝑛 reaches its maximum at neutron midshell (𝑁 = 104). Several midmass vibrational nuclei have 𝑀 𝑛 /𝑀 𝑝 values significantly below 𝑁/𝑍, which is not consistent with the expectation that 𝑀 𝑛 /𝑀 𝑝 = 𝑁/𝑍 in such nuclei. As a result, a shell-model investigation of these observations might yield insights about this behavior.
We calculate basis-space converged neutrinoless ββ-decay nuclear matrix elements for the lightest candidates: 48 Ca, 76 Ge, and 82 Se. Starting from initial two- and three-nucleon forces, we apply the ab initio in-medium similarity renormalization group to construct valence-space Hamiltonians and consistently transformed ββ-decay operators. Here, we find that the tensor component is non-negligible in 76 Ge and 82 Se, and the resulting nuclear matrix elements are overall 25%–45% smaller than those obtained from the phenomenological shell model. While a final matrix element with uncertainties still requires substantial developments, this work nevertheless opens a path toward a true first-principles calculation of neutrinoless ββ decay in all nuclei relevant for ongoing large-scale searches.
The cross sections of certain inelastic atomic collision processes can be determined from the matrix elements of the collision operator, d/dR, where R is the separation of the colliding atoms. At one extreme, the matrix element may pass through zero near a pseudocrossing of potential energy curves, while at the other extreme, it may pass through a maximum. The resulting cross sections are entirely different in magnitude and in energy dependence. An attempt is made to predict the qualitative behavior of the collision matrix elements with variations in R from an analysis of the Born-Oppenheimer adiabatic Hamiltonian. The modifications caused by a second pseudocrossing with a third adiabatic state are studied.
We report on the measurement of spin-density matrix elements for the K∗(892) → K+π0 photoproduction process with a recoiling Λ hyperon. The measurement used high-statistics GlueX data from photoproduction off a proton target at photon energies between 8.2 and 8.8 GeV, in a kinematic regime dominated by t-channel exchange processes. GlueX is a fixed-target experiment with a linearly polarized photon beam and a nearly 4π hermetic detector, allowing extraction of the full spin observ- ables from both the K+π0 and Λ → π−p systems. The GlueX data show clear evidence of the K∗(892) resonance in the K+π0 mass spectrum, along with additional structures at higher masses that may correspond to K∗0 or K∗2 resonances. The precise measurement of the spin-density matrix elements of the K∗(892) could serve as a standard candle for con- tinuing studies of higher-mass excited K∗ states. In addition, t-channel K∗ photoproduction differs from previous measurements of non-strange vector meson photoproduction, as it is free from Pomeron exchange and is virtually unexplored in this energy regime. Furthermore, the unprece- dented data collected by GlueX allow for investigations of correlations between the Λ polarization and the spin-density matrix elements of K∗ production, paving the way for future partial-wave analyses with full spin information involving the recoil hyperon.
In a recent paper, “Nonlocal nucleon matrix elements in the rest frame” [1], it was observed that the next-to-leading order calculations of the renormalization factor can describe, to a few percent accuracy, the logarithm of the lattice quantum chromodynamics (QCD) rest frame matrix elements with separations up to distances of 0.6 fm on multiple lattice spacings. We argue that perturbative QCD breaks down at such a distance scale after resumming the associated large logarithms, while the Ansatz used in the analysis there is not justified in perturbation theory. Besides, we explain the observation in Ref. [1] and demonstrate that the Ansatz fails to describe the data for 𝑧 > 0.3 fm, showing an opposite trend. Finally, although Ref. [1] proposes multiplying the Ansatz by a Gaussian correction model, which is shown to reduce the discrepancy with the data, this does not legitimize the use of perturbative QCD at such distance scales.
Demonstration that the vanishing of certain coupling matrix elements at level crossings follow from angular momentum commutation relations. A magnetic dipole transition having delta M = plus or minus 1, induced near a crossing of the levels in a nonzero magnetic field, is found to have a dipole matrix element comparable to or smaller than the quotient of the level separation and the field. This result also applies in the analogous electric field electric dipole case.
Neutrinoless double-beta decay (0vbb) is a hypothetical nuclear decay that is only possible if the neutrino is a Majorana fermion. Experimental searches for this process with ever-increasing sensitivity have placed strong constraints on the 0vbb half-lives of relevant isotopes. Relating these experimental half-lives to the underlying particle physics -- the effective Majorana mass of the neutrino -- requires understanding of the nuclear matrix elements for the transition. These matrix elements can be computed within a nuclear effective field theory framework, but input from lattice QCD is necessary to constrain low-energy constants relevant for the decay. This talk will discuss calculations of these matrix elements using lattice QCD and the implications for determination of nuclear EFT parameters.
Variational principles are considered for the approximate evaluation of the diagonal matrix elements of an arbitrary known linear Hermitian operator. A method is derived that is immediately applicable to the variational determination of both the off-diagonal and diagonal matrix elements of normal and modified Green's functions.