Abstract. A symbolic programming package NCoperators with applications to atomic physics is introduced. The package runs over Mathematica and it implements NCAlgebra, the noncommutative algebra package. NCoperators features the algebra of irreducible tensor operators, the second quantization representation, the angular momentum theory, and the effective operator approach exploited in many-body perturbation theory, including Wick's theorem. The comprehensiveness is yet another characteristic feature of the present package: The generation of expressions is performed in a way as if it were done by hand. Although the theoretical atomic spectroscopy is a direct target of NCoperators, the package, with minor modifications, if any, is believed to appropriate other areas of theoretical physics as well.
IntroductionTo this day, a number of Mathematica sources featuring the properties of creation and annihilation operator products are observed in literature:Against the variety of areas of theoretical physics the packages supply with, none of them provide an opportunity to manage the Racah algebra [6-8] in a fully compatible fashion; the capability of constructing the irreducible tensor form of complex operators that act on the basis of many-electron open-shell wave functions is advantageous in a special manner. The package NCoperators is just what one needs to solve this kind of tasks [9] along with many other ones, eg. the summation of Clebsch-Gordan products to the 3nj-symbols [10]. Whilst the Racah package [11][12][13], based on Maple, is known to handle several similar procedures such as the calculation of transformation matrices and recoupling coefficients, evaluation of many-particle matrix elements etc., the advantages and distinctive features of NCoperators are that (i) it manages the reduction of the product of more than two SU(2)-irreducible representations therefore giving an opportunity to calculate the matrix elements on the basis of many-open-shell wave functions (eg. six open shells appear even in the second-order Rayleigh-Schrödinger perturbation theory based on effective operator approach [14]); (ii) it is based on the algebraic manipulations of the quantities considered neither the diagrammatic representations such as the Jucys graphs exploited in the angular momentum theory [15], though the diagrammatic visualization, to some extent, is implemented as well; (iii) the output of expressions, including Clebsch-Gordan coefficients and 3nj-symbols, obtained within Mathematica interface fits the standard text-based form; (iv) Racah refers to Varshalovich et al [16], while NCoperators implements the sum rules given by Jucys et al [15; 17]; in many