A Mathematica script for harmonic oscillator nuclear matrix elements arising in semileptonic electroweak interactions

Abstract

Semi-leptonic electroweak interactions in nuclei - such as βdecay, μcapture, charged- and neutral-current neutrino reactions, and electron scattering - are described by a set of multipole operators carrying definite parity and angular momentum, obtained by projection from the underlying nuclear charge and three-current operators. If these nuclear operators are approximated by their one-body forms and expanded in the nucleon velocity through order |\vec{p}|/M, where \vec{p} and M are the nucleon momentum and mass, a set of seven multipole operators is obtained. Nuclear structure calculations are often performed in a basis of Slater determinants formed from harmonic oscillator orbitals, a choice that allows translational invariance to be preserved. Harmonic-oscillator single-particle matrix elements of the multipole operators can be evaluated analytically and expressed in terms of finite polynomials in q^2, where q is the magnitude of the three-momentum transfer. While results for such matrix elements are available in tabular form, with certain restriction on quantum numbers, the task of determining the analytic form of a response function can still be quite tedious, requiring the folding of the tabulated matrix elements with the nuclear density matrix, and subsequent algebra to evaluate products of operators. Here we provide a Mathematica script for generating these matrix elements, which will allow users to carry out all such calculations by symbolic manipulation. This will eliminate the errors that may accompany hand calculations and speed the calculation of electroweak nuclear cross sections and rates. We illustrate the use of the new script by calculating the cross sections for charged- and neutral-current neutrino scattering in ^{12}C.

15 pages, 2 tables, Mathematica notebook included in the form of figures. Mathematica package and documentation available at http://www.int.washington.edu/users/lunardi/7o.htm. Replaced version has improved graphics; text unchanged

Cited by in corpus (12)

• The Effective Field Theory of Dark Matter Direct Detection
• Model-independent WIMP Scattering Responses and Event Rates: A Mathematica Package for Experimental Analysis
• Coherent elastic neutrino-nucleus scattering: EFT analysis and nuclear responses
• Analysis strategies for general spin-independent WIMP-nucleus scattering
• Nuclear structure aspects of spin-independent WIMP scattering off xenon
• An $\text{ab initio}$ strategy for taming nuclear-structure dependence of $V_{ud}$ extractions: the ${}^{10}\mathrm{C}\to {}^{10}\mathrm{B}$ superallowed transition
• Nuclear-level Effective Theory of $\mu \to e$ Conversion: Formalism and Applications
• Inelastic nuclear scattering from neutrinos and dark matter
• Multipole decomposition of tensor interactions of fermionic probes with composite particles and BSM signatures in nuclear reactions
• A formalism to assess the accuracy of nuclear-structure weak interaction effects in precision $\beta$-decay studies
• Translationally invariant matrix elements of general one-body operators
• Tensor interaction in coherent elastic neutrino-nucleus scattering