An atomic spectral line is characteristic of the element producing the spectrum. The line also depends on the isotope. The program ris3 (Relativistic Isotope Shift) calculates the electron density at the origin and the normal and specific mass shift parameters. Combining these electronic quantities with available nuclear data, isotope-dependent energy level shifts are determined.Keywords: Isotope shift, field shift, mass shift, normal mass shift, specific mass shift, relativistic nuclear recoil, nuclear radius The nuclear motion and volume effects are treated in first order perturbation theory. Taking the zero-order wave function in terms of a configuration state expansion |Ψ P JM J = µ c µ |Φ(γ µ P JM J ) , where P , J and M J are, respectively, the parity and angular quantum numbers, the electron density at the nucleus and the normal and specific mass shift parameters may generally be expressed as µ,ν c µ c ν γ µ P JM J |ϑ|γ ν P JM J where ϑ is the relevant operator. The matrix elements, in turn, can be expressed as sums over radial integrals multiplied by angular coefficients. All the angular coefficients are calculated using routines from grasp2K version 1 1 package [1].
NEW VERSION PROGRAM SUMMARY
Reasons for the new version:This new version takes the nuclear recoil corrections into account within the (αZ) 4 m 2 /M approximation [2] and also allows a storage of the angular coefficients for a series of calculations within a given isoelectronic sequence. Furthermore, the program jj2lsj, a module of grasp2K version 1 1 toolkit that allows a transformation of ASFs from a jj-coupled CSF basis into an LSJ-coupled CSF basis, has been especially adapted to present ris3 results using LSJ labels of the states. This additional tool is called ris3 lsj.
Summary of revisions:This version is compatible with the new angular approach of grasp2K version 1 1 package [1] and can store necessary angular coefficients. According to the formalism of the relativistic nuclear recoil, the "uncorrected" expression of the normal mass shift has been fundamentally modified compared with its expression in [3].
Restrictions:The complexity of the cases that can be handled is entirely determined by grasp2K package [1] used for the generation of the electronic wave functions. Unusual features: Angular data stored on disc and can be reused. LSJ labels are used for the states. Additional comments: ??
2Running time:As an example, we evaluated the isotope shifts parameters and the electron density at the origin using the wave functions of Be-like system. We used MCDHF wave function built on a complete active space (CAS) with n = 8 (296 626 CSFs -62 orbitals) that contains 3 non-interacting blocks of given parity and J values involving 6 different eigenvalues in total. Calculations take around 10 hours on one AMD Opteron 6100 @ 2.3GHz CPU with 8 cores (64GB DDR3 RAM 1.333GHz). If angular files are available the time is reduced to 20 minutes. The storage of the angular data takes 139Mb and 7.2Gb for the one-body and the two-body elements, r...