Hyperfine structure parametrisation in Maple

Abstract

In hyperfine structure examinations, routine high resolution spectroscopy methods have to be combined with exact fine structure calculations. The so-called magnetic A and electric B factor of the fine structure levels allow to check for a correct fine structure analysis, to find errors in the level designation, to find new levels and to probe the electron wavefunctions and its mixing coefficients. This is done by parametrisation of these factors into different contributions of the subshell electrons, which are split further into their radial and spin–angular part. Due to the routine with which hyperfine structure measurements are done, a tool for keeping the necessary information together, performing checks online with the experiment and deriving standard quantities is of great help. Maple [Maple is a registered trademark of Waterloo Maple Inc.] is a highly-developed symbolic programming language, often referred to as the pocket calculator of the future. Packages for theoretical atomic calculation exist ( Racah and Jucys) and the language meets all the requirements to keep and present information accessible for the user in a fast and practical way. We slightly extended the Racah package [S. Fritzsche, Comput. Phys. Comm. 103 (1997) 51] and set up an environment for experimental hyperfine structure calculations, the Hfs package. Supplying the fine structure and nuclear data, one is in the position to obtain information about the hyperfine spectrum, the different contributions to the splitting and to perform a least square fit of the radial parameters based on the semiempirical method. Experimentalist as well as theoretical physicist can do a complete hyperfine structure analysis using Maple. Program summary Title of program: H fs Catalogue number: ADXD Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADXD Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Licensing provisions: none Computers for which the program is designed: All computers with a license of the computer algebra package Maple Installations: University of Kassel (Germany) Operating systems under which the program has been tested: Linux 9.0 Program language used: Maple, Release 7, 8 and 9 Memory required to execute with typical data: 5 MB No. of lines in distributed program, including test data, etc.: 34 300 No. of bytes in distributed program, including test data, etc.: 954 196 Distribution format: tar.gz Nature of the physical problem: Atomic state functions of an many configuration many electron atom with several open shells are defined by a number of quantum numbers, by their coupling and selection rules such as the Pauli exclusion principal or parity conservation. The matrix elements of any one-particle operator acting on these wavefunctions can be analytically integrated up to the radial part [G. Gaigalas, O. Scharf, S. Fritzsche, Central European J. Phys. 2 (2004) 720]. The decoupling of the interacting electrons is general, the obtained submatrix element holds all the peculiarities of the operator in question. These so-called submatrix elements are the key to do hyperfine structure calculations. The interaction between the electrons and the atomic nucleus leads to an additional splitting of the fine structure lines, the hyperfine structure. The leading components are the magnetic dipole interaction defining the so-called A factor and the electric quadrupole interaction, defining the so-called B factor. They express the energetic splitting of the spectral lines. Moreover, they are obtained directly by experiments and can be calculated theoretically in an ab initio approach. A semiempirical approach allows the fitting of the radial parts of the wavefunction to the experimentally obtained A and B factors. Method of solution: Extending the existing csf_LS() and asf_LS() to several open shells and implementing a data structure level_LS() for the fine structure level, the atomic environment is defined in Maple. It is used in a general approach to decouple the interacting shells for any one-particle operator. Further submatrix elements for the magnetic dipole and electric quadrupole interaction are implemented, allowing to calculate the A and B factors up to the radial part. Several procedures for standard quantities of the hyperfine structure are defined, too. The calculations are accelerated by using a hyper-geometric approach for three, six and nine symbols. Restrictions onto the complexity of the problem: Only atomic state functions in nonrelativistic LS-coupling with states having l ⩽ 3 are supported. Typical running time: The program replies promptly on most requests. The least square fit depends heavily on the number of levels and can take a few minutes.