Abstract
Semileptonic electroweak interactions in nucleisuch as β decay, μ capture, charged and neutralcurrent neutrino reactions, and electron scatteringare described by a set of multipole operators carrying definite parity and angular momentum, obtained by projection from the underlying nuclear charge and threecurrent operators. If these nuclear operators are approximated by their onebody forms and expanded in the nucleon velocity through order  over(p, →)  / M, where over(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. Harmonicoscillator singleparticle 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 threemomentum 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 neutralcurrent neutrino scattering in ^{12}C. Program summary: Program title: SevenOperators. Catalogue identifier: AEAY_v1_0. Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEAY_v1_0.html. Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland. Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html. No. of lines in distributed program, including test data, etc.: 2227. No. of bytes in distributed program, including test data, etc.: 19 382. Distribution format: tar.gz. Programming language: Mathematica. Computer: Any computer running Mathematica; tested on Mac OS X PowerPC (32bit) running Mathematica 6.0.0. Operating system: Any running Mathematica. RAM: Memory requirements determined by Mathematica; 512 MB or greater RAM and hard drive space of at least 3.0 GB recommended. Classification: 17.16, 17.19. Nature of problem: Algebraic evaluation of harmonic oscillator nuclear matrix elements for the onebody multipole operators governing semileptonic weak interactions, such as charged or neutralcurrent neutrino scattering off nuclei. Solution method: Mathematica evaluation of associated angular momentum algebra and spherical Bessel function radial integrals. Running time: Depends on the complexity of the onebody density matrix employed, but times of a few seconds are typical.
Original language  English (US) 

Pages (fromto)  345358 
Number of pages  14 
Journal  Computer Physics Communications 
Volume  179 
Issue number  5 
DOIs  
State  Published  Sep 1 2008 
Keywords
 Harmonic oscillator shellmodel matrix elements
 Semileptonic weak and electromagnetic interactions
ASJC Scopus subject areas
 Hardware and Architecture
 Physics and Astronomy(all)
Fingerprint
Dive into the research topics of 'SevenOperators, a Mathematica script for harmonic oscillator nuclear matrix elements arising in semileptonic electroweak interactions'. Together they form a unique fingerprint.Datasets

SevenOperators, a Mathematica script for harmonic oscillator nuclear matrix elements arising in semileptonic electroweak interactions
Haxton, W. (Contributor) & Lunardini, C. (Contributor), Mendeley Data, 2008
DOI: 10.17632/hwzx6vfjmv.1, https://data.mendeley.com/datasets/hwzx6vfjmv
Dataset