modern atomic spectroscopy studies the structure and properties of practically any atom of the periodic table as well as of ions of any ionization degree .
particular attention is paid to their energy spectra . for the investigations of many - electron atoms and ions ,
it is of great importance to combine experimental and theoretical methods .
nowadays the possibilities of theoretical spectroscopy are much enlarged thanks to the wide use of powerful computers .
theoretical methods utilized must be fairly universal and must ensure reasonably accurate values of physical quantities studied .
many - electron atom usually is considered as many - body problem and is described by the wave function constructed from the wave functions of one electron , moving in the central nuclear charge field and in the screening field of the remaining electrons
. then the wave function of this electron may be represented as a product of radial and spin - angular parts .
the radial part is usually found by solving various modifications of the hartree - fock equations and can be represented in a numerical or analytical forms ( froese fischer  @xcite ) whereas the angular part is expressed in terms of spherical functions .
then the wave function of the whole atom can be constructed in some standard way ( cowan  @xcite , jucys and savukynas  @xcite , nikitin and rudzikas  @xcite ) starting with these one - electron functions and may be used further on for the calculations of any matrix elements representing physical quantities .    for obtaining the values of atomic quantities
it is necessary to solve so - called eigenvalue problem    @xmath1    where @xmath2 is the wave function of the system under investigation and @xmath3 is its hamiltonian . in order to obtain accurate values of atomic quantities
it is necessary to account for relativistic and correlation effects .
it turned out that for very large variety of atoms and their ionization degrees the relativistic effects may be taken into account fairly accurately as breit - pauli corrections ( nikitin and rudzikas  @xcite , rudzikas  @xcite ) .
it is convenient to present the hamiltonian as consisting of two parts in breit - pauli approach , namely ,    @xmath4    where @xmath5 is the non - relativistic hamiltonian and @xmath6 stands for the relativistic contibution .
non - relativistic hamiltonian is the sum of the kinetic energy of electrons @xmath7 , their potential energy @xmath8 , and the electrostatic electron interaction @xmath9    @xmath10    the first two terms , @xmath7 and @xmath8 are the one - particle operators , whereas the third term , @xmath9 , is two - particle operator .
the @xmath6 may also be subdivided into non - fine structure and fine structure contributions    @xmath11    the non - fine structure contributions    @xmath12    shift non - relativistic energy levels without any splitting of them . the mass - velocity term
@xmath13 describes the variation of mass with velocity .
the one- and two - body darwin ( contact ) terms @xmath14 and @xmath15 are the corrections of the one - electron dirac equation due to the retardation of the electromagnetic field produced by an electron ( contact interaction ) .
the spin - spin - contact term @xmath16 accounts for the interaction of the spin magnetic moments of two electrons occupying the same space .
the orbit - orbit interaction @xmath17 accounts for the interaction of two orbital moments .
the fine - structure contributions    @xmath18    split the non - relativistic energy levels ( terms ) into a series of closely - spaced levels .
the most important of these is the nuclear spin - orbit interaction @xmath19 ( spin - own - orbit ) representing the interaction of the spin and angular magnetic moments of an electron in the field of the nucleus .
the spin - other - orbit @xmath20 and spin - spin @xmath21 contributions in a rough sense may be viewed as corrections to the nuclear spin - orbit interaction due to the presence of the other electrons in the system .
unfortunately , practical calculations show that all realistic atomic hamiltonians do not lead straightforwardly to eigenvalue problem ( [ eq : sch ] ) .
actually we have to calculate all non - zero matrix elements of the hamiltonian considered including those non - diagonal with respect to electronic configurations , then to form energy matrix , to diagonalize it , obtaining in this way the values of the energy levels as well as the eigenfunctions ( the wave functions in the intermediate coupling scheme ) .
the latter may be used then to calculate electronic transitions as well as the other properties and processes .
such a necessity raises special requirements for the theory .
the total matrix element of each term of the energy operator in the case of complex electronic configuration will consist of matrix elements , describing the interaction inside each shell of equivalent electrons as well as between these shells . going beyond the single - configuration approximation we have to be able to take into account in the same way non - diagonal , with respect to configurations , matrix elements .    to find
the expressions for the matrix elements of all terms of the hamiltonian considered for complex electronic configurations , having several open shells , is a task very far from the trivial one .
a considerable part of the effort must be devoted to coping with integrations over spin - angular variables , occurring in the matrix elements of the operators under consideration .
this paper presents the general methodology , leading to optimal expressions for operators and matrix elements .
a number of methodologies to calculate the angular parts of matrix elements exists in literature .
many existing codes for integrating the spin - angular parts of matrix elements ( glass  @xcite , glass and hibbert  @xcite , grant  @xcite , burke _
@xcite ) are based on the computational scheme proposed by fano  @xcite .
this methodology is based on having the total wavefunction of an atom built from the antisymmetrized wavefunctions of separate shells , and this antisymmetrization is done via coefficients of fractional parentage .
the shells are coupled one to another via their angular momenta .
so , the finding of matrix elements amounts to finding the recoupling matrices and the coefficients of fractional parentage .
suppose that we have a bra function with @xmath22 shells in @xmath0 coupling : @xmath23{c } \psi _ u^{bra}\left ( lsm_lm_s\right ) \equiv \\ \equiv ( n_1l_1n_2l_2 ... n_ul_u\alpha _ 1l_1s_1\alpha _ 2l_2s_2 ... \alpha _ ul_us_u{\cal a}lsm_lm_s| \end{array}\ ] ] and a ket function : @xmath24{c } \psi _ u^{ket}\left ( l^{\prime } s^{\prime } m^{\prime}_lm^{\prime}_s\right ) \equiv \\
n_ul_u\alpha _ 1^{\prime } l_1^{\prime } s_1^{\prime } \alpha _ 2^{\prime } l_2^{\prime } s_2^{\prime } ... \alpha _
u^{\prime } l_u^{\prime } s_u^{\prime } { \cal a}^{\prime } l^{\prime } s^{\prime}m^{\prime}_lm^{\prime}_s ) , \end{array}\ ] ] where @xmath25 stands for all intermediate quantum numbers , depending on the order of coupling of momenta @xmath26 . using the wigner - eckart theorem in @xmath0 space we shift from the matrix element of any two - particle operator @xmath27 between functions ( [ eq : ma - a ] ) and ( [ eq : ma - b ] ) to the submatrix element ( reduced matrix )    @xmath28 of this operator .
when the two - particle operator acts upon four distinct shells , then the finding general expressions of matrix elements , according to methodology by fano  @xcite , is based upon the formula    @xmath29{c } ( \psi _ u^{bra}\left ( ls\right ) || \widehat{g}^{\left ( \kappa _ 1\kappa _ 2k,\sigma_1\sigma _ 2k\right ) }     \sim \displaystyle { \sum_{n_{i } \lambda _ { i } , n_{j } \lambda _ { j } , n_{i^{\prime } } \lambda _ { i^{\prime } } , n_{j^{\prime } } \lambda _ { j^{\prime } } } } \left ( -1\right ) ^\delta \left [ n_{i } \left ( n_{j}-\delta(i , j ) \right ) n_{i^{\prime}}^{\prime } \left ( n_{j^{\prime}}^{\prime}-\delta(i^{\prime},j^{\prime } ) \right ) \right]^{1/2 } \times \\ \times \left ( l_{i}^{n_{i}}\;\alpha_{i } l_{i}s_{i}||l_{i}^{n_{i}-1}\;\left ( \alpha_{i}^{\prime } l^{\prime}_{i } s^{\prime } _ { i } \right ) , l_{i } \right ) \left ( l_{j}^{n_{j}}\;\alpha_{j } l_{j}s_{j}||l_{j}^{n_{j}-1}\;\left ( \alpha_{j}^{\prime } l^{\prime}_{j } s^{\prime } _ { j } \right ) , l_{j } \right ) \times \\ \times \left ( l_{i^{\prime } } ^{n_{i^{\prime } } } \;\alpha_{i^{\prime } } l_{i^{\prime } } s_{i^{\prime } } ||l_{i^{\prime } } ^{n_{i^{\prime } } -1}\;\left ( \alpha_{i^{\prime } } ^{\prime } l^{\prime}_{i^{\prime } } s^{\prime } _ { i^{\prime } } \right ) , l_{i^{\prime } } \right ) \left ( l_{j^{\prime } } ^{n_{j^{\prime } } } \;\alpha_{j^{\prime } } l_{j^{\prime } } s_{j^{\prime } } ||l_{j^{\prime } } ^{n_{j^{\prime } } -1}\;\left ( \alpha_{j^{\prime } } ^{\prime } l^{\prime}_{j^{\prime } } s^{\prime } _ { j^{\prime } } \right ) , l_{j^{\prime } } \right ) \times \\ \times \left\ { r_{d } \left ( \lambda _ i,\lambda _ j,\lambda _ i^{\prime } , \lambda _ j^{\prime } , \kappa_1 , \kappa_2 , \kappa , \sigma_1 , \sigma_2 , \sigma , \lambda ^{bra } , \lambda ^{ket } \right ) \right .
\\ \times \left ( 1 + \delta(i , j ) \delta(i^{\prime } , j^{\prime } ) \right)^{-1 } \left ( n_i\lambda _ in_j\lambda _ j|| g^{\left ( \kappa _ 1\kappa _ 2k,\sigma_1\sigma _ 2k\right ) }    - r_{e }   \left ( \lambda _ i,\lambda _ j,\lambda _ i^{\prime } , \lambda _ j^{\prime } , \kappa_1 , \kappa_2 , \kappa , \sigma_1 , \sigma_2 , \sigma , \lambda ^{bra } , \lambda ^{ket } \right ) \\ \times   \left .
\left ( 1 + \delta(i , j ) ) ( 1 - \delta(i^{\prime } , j^{\prime } ) \right ) \left ( n_i\lambda _ in_j\lambda _
j|| g^{\left ( \kappa _ 1\kappa _ 2k,\sigma_1\sigma _ 2k\right ) }    \right\ } , \end{array}\ ] ]    where @xmath30 is a phase factor , ( see for example in @xcite ) ,    @xmath31 is the array for the bra function shells terms , and similarly for @xmath32 .
the coefficient @xmath33 is a fractional parentage coefficient , and coefficients    @xmath34 and    @xmath35 are the recoupling matrices in @xmath36 and @xmath37 - spaces of direct and exchange terms , respectively . for more detailes on recoupling matrices
see in grant  @xcite , burke _
et al _  @xcite .
the summation in expression ( [ eq : ma - c ] ) implies that the quantum numbers @xmath38 , @xmath39 of all participating shells are included .
there are four such pairs of @xmath38 , @xmath39 in the sum .
in essence , the fano calculation scheme consists of evaluating recoupling matrices .
although such an approach uses classical racah algebra @xcite on the level of coefficients of fractional parentage , it may be necessary to carry out multiple summations over intermediate terms .
due to these summations and the complexity of the recoupling matrix itself , the associated computer codes become rather time consuming .
jucys and vizbarait .
e  @xcite proposed to use the two - electron coefficients of fractioanl parentage instead of ordinary ones , in matrix elements calculations , but even that did not solve the abovementioned problems . a solution to this problem
was found by burke _
they tabulated separate standard parts of recoupling matrices along with coefficients of fractional parentage at the beginning of a calculation and used them further on to calculate the needed coefficients .
computer codes by glass  @xcite , glass and hibbert  @xcite , grant  @xcite , burke _
@xcite utilize the program njsym ( burke  @xcite ) or njgraf ( bar - shalom and klapisch  @xcite ) for the calculation of recoupling matrices .
both are rather time consuming when calculating matrix elements of complex operators or electronic configurations with many open shells .    in order to simplify the calculations ,
cowan  @xcite suggested that matrix elements be grouped into  classes
( see cowan  @xcite figure 13 - 5 ) .
unfortunately , this approach was not generalized to all two - electron operators . perhaps for this reason
cowan s approach is not very popular although the program itself , based on this approach , is widely used .
many approaches for the calculation of spin - angular coefficients ( glass  @xcite , glass and hibbert  @xcite , grant  @xcite , burke _ et al _
@xcite ) are based on the usage of racah algebra only on the level of coefficients of fractional parentage .
a few authors ( jucys and savukynas  @xcite , cowan  @xcite ) utilize the unit tensors , simplifying the calculations in this way , because use can be made of the tables of unit tensors and selection rules can be used prior to computation to check whether the spin - angular coefficients are zero or not .
moreover , the recoupling matrices themselves have a simpler form .
unfortunately , these ideas were applied only to diagonal matrix elements with respect to configurations , though cowan  @xcite suggested the usage of unit tensors for non - diagonal ones as well .
all the above mentioned approaches were applied in the coordinate representation . the second quantization formalism ( judd  @xcite , rudzikas and kaniauskas  @xcite and rudzikas  @xcite ) has a number of advantages compared to coordinate representation .
first of all , it is much easier to find algebraic expressions for complex operators and their matrix elements , when relying on second quantization formalism .
it has contributed significantly to the successful development of perturbation theory ( see lindgren and morrison  @xcite , merkelis _ et al _  @xcite , and orthogonal operators ( uylings  @xcite ) , where three - particle operators already occur .
uylings  @xcite suggested a fairly simple approach for dealing with separate cases of three - particle operators .
moreover , in the second quantization approach the quasispin formalism was efficiently developed by innes  @xcite , pakauskas _ et al _
@xcite , rudzikas and kaniauskas  @xcite , fano and rau  @xcite .
the main advantage of this approach is that applying the quasispin method for calculating the matrix elements of any operator , we can use all advantages of the new version of racah algebra ( see rudzikas  @xcite ) for integration of spin - angular part of any one- and two - particle operator .
for example , the reduced coefficients of fractional parentage are independent of the occupation number of the shell ( see rudzikas and kaniauskas  @xcite , gaigalas _ et al _
all this enabled merkelis and gaigalas  @xcite to work out a general perturbation theory approach for complex cases of several open shells . in the paper by merkelis  @xcite a detailed review of a version of graphical methodology
is presented that allows one to represent the operators graphically and to find the matrix elements of these operators using diagrammatic technique .
the majority of methods and computer codes of finding angular coefficients discussed above were faced with a number of problems , the main of these being :    * the high demand of cpu time for calculating the angular parts of matrix elements even on the modern computer .
therefore the high accuracy of characteristics of atomic quantities somtimes is even unattainable .
* the methods are applicable in practice only for the comparatively simple systems , because in the programs based on classical racah algebra the treatment of recoupling matrices is rather complicated , especialy when finding the matrix elements of breit - pauli operators between complex configurations . * the configurations with open @xmath40-shells
must often be included in theoretical calculations .
this causes problems in a number of methodologies , because the complete account of @xmath40- shells implies using a large number of coefficients of fractional parentage .
gaigalas and rudzikas  @xcite , gaigalas _
@xcite , gaigalas _
@xcite , gaigalas _ et al _  @xcite and gaigalas and rudzikas  @xcite suggested an efficient and general approach for finding the spin - angular parts of matrix elements of atomic interactions , relying on the combination of the second - quantization approach in the coupled tensorial form , the generalized graphical technique and angular momentum theory in orbital , spin and quasispin spaces as well as on the symetry properties of the quantities considered .
this approach is free of previous shortcomings .
it is the main goal of this work to present all this methodology consistently and in a unifying manner , paying the special attention to its main ideas
. also , we aim at a detailed discussion of obtaining the efficient tensorial expressions of a two - particle operator , as well as the analytical expressions for recoupling matrices .
in addition to that , in this work we aim at comparing this methodology to other calculation of angular coefficients schemes .    from what is said above
we see that the treatment of angular parts of matrix elements is a many - sided question .
the classical scheme by fano was developed in various aspects .
the methods were developed that used racah algebra @xcite on a higher level , or used the new version of the racah algebra ( angular momentum theory in orbital , spin and quasispin spaces ) ( see rudzikas  @xcite ) .
the ways were searched for , to obtain the expressions that had simpler recoupling matrices . in order to discuss all that in more detail , to compare the existing methods of obtaining the angular parts and to mark the advantages of one method or another , the finding angular parts must be looked at from different angles .
thus , the first thing to discuss is the expression for any physical operator ( see sections 2,3,4 ) , because already the form of it determines the level of application of tensor algebra in the calculation of matrix elements of this operator .
also , it is very important to make clear to what extent the racah algebra is exploited in one or another methodology of angular parts treatment .
this helps to mark the advantages of one method against another ( see section 5 ) .
in addition , it is of importance to compare the ways of calculating the recoupling matrix in various methods ( see section 6 ) .
it is well known in the literature that a scalar two - particle operator may be presented the following tensorial form ( see jucys and savukynas  @xcite , glass  @xcite ) :    @xmath41_{p - p}^{(k~k ) } , \ ] ]    where @xmath42 is the radial part of operator , @xmath43 is a tensor acting upon the orbital and spin variables of the @xmath44-th function , @xmath45 , @xmath46 are the ranks of operator acting in orbital space , and @xmath47 , @xmath48 are the ranks of operator acting in spin space .
all the above mentioned approaches were usually applied in the coordinate representation .
now we will investigate into the second quantization formalism , which is broadly applied in atomic physics as well .
a two - particle operator in second quantization method is written as follows :    @xmath49    where @xmath50 , @xmath51 is the two - electron matrix element of operator @xmath52 , and @xmath53 is the electron creation and @xmath54 electron annihilation operators .
meanwhile two tensorial forms are well known in second quantization . in the first form
the operators of second quantization follow in the normal order :    @xmath55{c } g_i^{\left ( \kappa_1 \kappa_2 k , \sigma_1 \sigma_2 k \right)}= \displaystyle { \sum_{n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } } } \widehat{g}_i(ij , i^{\prime } j^{\prime } ) = \\
= -\frac 12 \displaystyle { \sum_{n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } } } \displaystyle { \sum_{\kappa _ { 12}\kappa _ { 12}^{^{\prime } } \sigma _ { 12}\sigma _ { 12}^{^{\prime } } } } \displaystyle { \sum_p}\left ( -1\right ) ^{k - p}\left [ \kappa _ { 12},\kappa _ { 12}^{\prime } , \sigma _ { 12},\sigma _ { 12}^{\prime } \right ] ^{1/2}\times \\ \times \left ( n_i\lambda _ in_j\lambda _ j||g^{\left ( \kappa _ 1\kappa _ 2k,\sigma _ 1\sigma _ 2k\right ) } ||n_i^{\prime } \lambda _ i^{\prime } n_j^{\prime } \lambda _ j^{\prime } \right ) \times \\ \times \left\ { \begin{array}{ccc } l_i^{\prime } & l_j^{\prime } & \kappa _ { 12}^{\prime } \\
\kappa _ 1 & \kappa _ 2 & k \\ l_i & l_j & \kappa _ { 12 } \end{array } \right\ } \left\ { \begin{array}{ccc } s & s & \sigma _ { 12}^{\prime }
\\ \sigma _ 1 & \sigma _ 2 & k \\ s & s & \sigma _ { 12 } \end{array } \right\ } \times \\ \times \left [ \left [ a^{\left ( \lambda _ i\right ) } \times a^{\left ( \lambda _ j\right ) } \right ] ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times \left [ \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } \times \stackrel{% \sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \right ] _ { p - p}^{\left ( kk\right ) } , \end{array}\ ] ]    where @xmath56 = \left ( 2a+1\right ) \left ( 2b+1\right ) $ ] , @xmath57 , @xmath58 is the two - electron submatrix ( reduced matrix ) element of operator @xmath59 and tensor @xmath60 is defined as ( see for example rudzikas  @xcite )    @xmath61    the product of tensors @xmath62 ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times \left [ \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } \times \stackrel{% \sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \right ] _ { p - p}^{\left ( kk\right ) } $ ] denotes tensorial part of operator @xmath63 .    in another form
the second quantization operators are coupled by pairs consisting of electron creation and annihilation operators . in
coupled tensorial form :    @xmath64{c } g_{ii}^{\left ( \kappa_1 \kappa_2 \kappa , \sigma_1 \sigma_2 \sigma \right)}= \displaystyle { \sum_{n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } } } \widehat{g}_{ii}(ij , i^{\prime } j^{\prime } ) = \\
= \frac 12 \displaystyle { \sum_{n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } } } \displaystyle { \sum_p}\left ( -1\right ) ^{k - p}\left ( n_i\lambda _ in_j\lambda _ j||g^{\left ( \kappa _ 1\kappa _ 2k,\sigma _ 1\sigma _ 2k\right ) } ||n_i^{\prime }
\lambda _ i^{\prime } n_j^{\prime } \lambda _ j^{\prime } \right ) \times \\ \times \left\ { \left [ \kappa _ 1,\kappa _ 2,\sigma _ 1,\sigma _ 2\right ] ^{-1/2}\times \left [ \left [ a^{\left ( \lambda _ i\right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } \right ] ^{\left ( \kappa _ 1\sigma _ 1\right ) } \times \left [ a^{\left ( \lambda _ j\right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } \right ] _ { p - p}^{\left ( kk\right ) } - \right .
-\left ( -1\right ) ^{l_i+l_j^{\prime } } \left\ { \begin{array}{ccc } \kappa _ 1 & \kappa _ 2 & k
\\ l_j^{\prime } & l_i & l_j \end{array } \right\ } \left\ { \begin{array}{ccc } \sigma _ 1 & \sigma _ 2 & k \\ s & s & s \end{array } \right\ } \left [ a^{\left ( \lambda _ i\right ) } \times \stackrel{\sim } { a}% ^{\left ( \lambda _ j^{\prime } \right ) } \right ] _ { p - p}^{\left ( kk\right ) } \delta \left ( n_jl_j , n_i^{\prime } l_i^{\prime } \right ) \right\}. \end{array}\ ] ]    the expression ( [ eq : op - c ] ) consists of only one tensorial product whereas ( [ eq : op - e ] ) has two , but the summation in the first formula is also over intermediate ranks @xmath65 , @xmath66 , @xmath67 and @xmath68 , complicating in this way the calculations
. the advantages or disadvantages of these alternative forms of arbitrary two - electron operator may be revealed in practical applications .    in these forms
the product of second quantization operators denotes tensorial part of operator @xmath52 .
for instance , the tensorial structure of electrostatic ( coulomb ) electron interaction operator is @xmath69 ( jucys and savukynas  @xcite ) , and only the two - electron submatrix elements @xmath70 of these operators are different . in the case of electrostatic interaction :    @xmath71{c } \left ( n_i\lambda _ in_j\lambda _ j||g_{coulomb}^{\left ( kk0,000\right ) } ||n_i^{\prime } \lambda _ i^{\prime } n_j^{\prime } \lambda _ j^{\prime } \right ) = \\
= 2\left [ k\right ] ^{1/2}\left ( l_i||c^{\left ( k\right ) } ||l_i^{\prime } \right ) \left ( l_j||c^{\left ( k\right ) } ||l_j^{\prime } \right ) r_k\left ( n_il_in_i^{\prime } l_i^{\prime } , n_jl_jn_j^{\prime } l_j^{\prime } \right ) .
\end{array}\ ] ]    from ( [ eq : op - f ] ) , by ( [ eq : op - c ] ) and ( [ eq : op - e ] ) , we finally obtain the following two secondly quantized expressions for coulomb operator ( see merkelis _
@xcite ) :    @xmath72{c } \widehat{v}_i=-\frac 12 \displaystyle { \sum_{n_il_in_jl_jn_i^{\prime } l_i^{\prime } n_j^{\prime } l_j^{\prime } } } \displaystyle { \sum_{\kappa _ { 12}\sigma _ { 12}k}}\left ( -1\right ) ^{l_j+l_i^{\prime } + k+\kappa _ { 12}}\left [ \kappa _ { 12},\sigma _ { 12}\right ] ^{1/2}\left ( l_i||c^{\left ( k\right ) } ||l_i^{\prime } \right ) \times \\ \times \left ( l_j||c^{\left ( k\right ) } ||l_j^{\prime } \right ) r_k\left ( n_il_in_i^{\prime } l_i^{\prime } , n_jl_jn_j^{\prime } l_j^{\prime } \right ) \left\ { \begin{array}{ccc } l_i & l_i^{\prime } & k
\\ l_j^{\prime } & l_i & \kappa _ { 12 } \end{array } \right\ } \times \\ \times \left [ \left [ a^{\left ( \lambda _ i\right ) } \times a^{\left ( \lambda _ j\right ) } \right ] ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times \left [ \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } \times \stackrel{% \sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \right ] ^{\left ( 00\right ) } , \end{array}\ ] ]    @xmath73{c } \widehat{v}_{ii}= \displaystyle { \sum_{n_il_in_jl_jn_i^{\prime } l_i^{\prime } n_j^{\prime } l_j^{\prime } } } \displaystyle { \sum_k}\left ( l_i||c^{\left ( k\right ) } ||l_i^{\prime } \right ) \left ( l_j||c^{\left ( k\right ) } ||l_j^{\prime } \right ) r_k\left ( n_il_in_i^{\prime } l_i^{\prime } , n_jl_jn_j^{\prime } l_j^{\prime } \right ) \times \\ \times \left\ { \left [ k\right ] ^{-1/2}\left [ \left [ a^{\left ( \lambda _ i\right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } \right ] ^{\left ( k0\right ) } \times \left [ a^{\left ( \lambda _ j\right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( k0\right ) } \right ] ^{\left ( 00\right ) } + \right . \\ \left .
+ \left ( 2\left [ l_i\right ] \right ) ^{-1/2}\left [ a^{\left ( \lambda _ i\right ) } \times \stackrel{\sim } { a}% ^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( 00\right ) } \delta \left ( n_jl_j , n_i^{\prime } l_i^{\prime } \right ) \right\ } , \end{array}\ ] ]    the tensorial expressions for orbit - orbit and other physical operators in second quantization form may be obtained in the same manner .
it is worth mentioning that the expressions ( [ eq : op - g ] ) and ( [ eq : op - h ] ) embrace , already in an operator form , both the diagonal interaction terms , relative to configurations , and the non - diagonal ones .
non - diagonal terms define the interaction between all the possible electron distributions over the configurations considered , differing by quantum numbers of not more than two electrons .    the merits of representing operators in one form or another ( [ eq : op - g ] ) or ( [ eq : op - h ] ) are mostly determined by the technique used to find their matrix elements and quantities in terms of which they are expressed .    in the paper gaigalas and rudzikas  @xcite it was shown that the tensorial forms ( [ eq : op - a ] ) , ( [ eq : op - c ] ) , ( [ eq : op - e ] ) of two - particle operator do not take full advantage of tensor algebra .
the most characteristic examples are when configurations considered have many open shells , or when the non - diagonal matrix elements are seeked .    in the paper gaigalas
_ et al _
@xcite the following optimal expression of two - particle operator is proposed , which allows one to make the most of the advantages of racah algebra ( see racah @xcite ) .
@xmath74{c } \widehat{g}^{\left ( \kappa_1 \kappa_2 k , \sigma_1 \sigma_2 k \right)}\sim \displaystyle { \sum_{\alpha } } \displaystyle { \sum_{\kappa _ { 12},\sigma _ { 12},\kappa _ { 12}^{\prime } , \sigma _ { 12}^{\prime } } } \theta \left ( \xi \right ) \left\ { a_{p ,- p}^{\left ( kk\right ) } \left ( n_\alpha \lambda _ \alpha , \xi \right ) \delta \left ( u,1\right ) \right . \\
+ \displaystyle { \sum_{\beta } } \left [ b^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \left ( n_\alpha \lambda _ \alpha , \xi \right ) \times c^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \left ( n_\beta \lambda _ \beta , \xi \right ) \right ] _ { p ,- p}^{\left ( kk\right ) } \delta \left ( u,2\right ) \\ + \displaystyle { \sum_{\beta \gamma } } \left [ \left [ d^{\left ( l_\alpha s\right ) } \times d^{\left ( l_\beta s\right ) } \right ] ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times e^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \left ( n_\gamma \lambda _ \gamma , \xi \right ) \right ] _ { p ,- p}^{\left ( kk\right ) } \delta \left ( u,3\right ) \\ \left
. + \displaystyle { \sum_{\beta \gamma \delta } } \left [ \left [ d^{\left ( l_\alpha s\right ) } \times d^{\left ( l_\beta s\right ) } \right ] ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times \left [ d^{\left ( l_\gamma s\right ) } \times d^{\left ( l_\delta s\right ) } \right ] ^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \right ] _ { p ,- p}^{\left ( kk\right ) } \delta \left ( u,4\right ) \right\ } .
\end{array}\ ] ]    whereas in traditional expressions , e. g. ( [ eq : op - c ] ) , the summation runs over the principle and the orbital quantum numbers of open shells without detailing these , in the expression written above the first term represents the case of a two - particle operator acting upon the same shell @xmath75 , the second term corresponds to operator @xmath76 acting upon two different shells @xmath75 , @xmath77 .
when operator @xmath76 acts upon three shells the third term in ( [ eq : op - i ] ) must be considered and when it acts upon four - the fourth one .
we define in this expression the shells @xmath75 , @xmath77 , @xmath78 , @xmath79 to be different .
the tensorial part of a two - particle operator is expressed in terms of operators of the type @xmath80 , @xmath81 , @xmath82 , @xmath83 , @xmath84 .
they correspond to one of the forms :    @xmath85    @xmath86 ^{\left ( \kappa _ 1\sigma _ 1\right ) } , \ ] ]    @xmath87 ^{\left ( \kappa _ 1\sigma _ 1\right ) } \right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } , \ ] ]    @xmath88 ^{\left ( \kappa _ 1\sigma _ 1\right ) } \times a_{m_{q3}}^{\left ( q\lambda \right ) } \right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } , \ ] ]    @xmath89 ^{\left ( \kappa _ 1\sigma _ 1\right ) } \times \left [ a_{m_{q3}}^{\left ( q\lambda \right ) } \times a_{m_{q4}}^{\left ( q\lambda \right ) } \right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } \right ] ^{\left ( kk\right ) } , \ ] ]    for example , if we take a two - particle operator acting upon two shells , then we see from expression ( [ eq : op - i ] ) that the angular part of two - particle operator is expressed via operators @xmath90 and @xmath91 . in a case when the operator @xmath76 acts in such a manner that two operators of second quantization act upon one shell and two act upon another , the @xmath90 and @xmath91 are expressed as ( [ eq : op - k ] ) .
but in a case when three operators of second quantization act upon one shell and one acts upon another , then @xmath90 and @xmath91 are expressed either as ( [ eq : op - j ] ) and ( [ eq : op - l ] ) or ( [ eq : op - j ] ) and ( [ eq : op - m ] ) .    in writing down the expressions ( [ eq : op - j ] ) - ( [ eq : op - n ] ) the quasispin formalism was used , where @xmath92 and @xmath93 are components of the tensor @xmath94 , having in additional quasispin space the rank @xmath95 and projections @xmath96 , i.e.    @xmath97    and    @xmath98    in the expression ( [ eq : op - i ] ) @xmath22 is the overall number of shells acted upon by a given tensorial product of creation / annihilation operators .
parameter @xmath99 implies the whole array of parameters ( and sometimes an internal summation over some of these is implied , as well ) that connect the amplitudes @xmath100 of tensorial products of creation / annihilation operators in the expression ( [ eq : op - i ] ) to these tensorial products ( see gaigalas _ et al _
also , attention must be paid to the fact that the ranks @xmath101 , @xmath102 , @xmath103 , @xmath104 , @xmath105 and @xmath106 are also included into the parameter @xmath99 .
the amplitudes @xmath107 are all proportional to the submatrix element of a two - particle operator @xmath108 ,    @xmath109    to obtain the expression of a concrete physical operator , analogous to expression ( [ eq : op - i ] ) , the tensorial structure of the operator and the two - particle matrix elements ( [ eq : op - r ] ) must be known .
we shall investigate this now .
the _ electrostatic ( coulomb ) electron interaction _
operator @xmath110 itself contains the tensorial structure    @xmath111{c } h^{coulomb } \equiv \displaystyle { \sum_{k } } h_{coulomb}^{(kk0,000 ) } \end{array}\ ] ]    and its submatrix element is    @xmath112{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{coulomb}^{(kk0,000 ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = \\ 2 [ k]^{1/2 } \left ( l_i\left\| c^{\left ( k \right ) } \right\| l_{i^{\prime } } \right ) \left ( l_j\left\| c^{\left ( k \right ) } \right\| l_{j^{\prime } } \right ) r_{k}\left ( n_il_i n_{i^{\prime } } l_{i^{\prime } } , n_jl_j n_{j^{\prime } } l_{j^{\prime } } \right ) .
\end{array}\ ] ]    the _ spin - spin _
operator @xmath113 itself contains tensorial structure of two different types , summed over @xmath114 ( gaigalas and rudzikas  @xcite ) ,    @xmath115{c } h^{ss } \equiv \displaystyle { \sum_{k } } \left [ h_{ss}^{(k+1 k-1 2,112 ) } + h_{ss}^{(k-1 k+1 2,112 ) } \right ] .
\end{array}\ ] ]    their submatrix elements are ( jucys and savukynas  @xcite )    @xmath116{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{ss}^{\left ( k+1k-12,112\right ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = \frac 3 { \sqrt{5}}\sqrt{\left ( 2k+3\right ) ^{\left ( 5\right ) } } \times   \\
\times \left ( l_i\left\| c^{\left ( k+1\right ) } \right\| l_{i^{\prime } } \right ) \left ( l_j\left\| c^{\left ( k-1\right ) } \right\| l_{j^{\prime } } \right ) n^{k-1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \end{array}\ ] ]    @xmath117{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{ss}^{\left ( k-1k+12,112\right ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = \frac 3 { \sqrt{5}}\sqrt{\left ( 2k+3\right ) ^{\left ( 5\right ) } } \times   \\
\times \left ( l_i\left\| c^{\left ( k-1\right ) } \right\| l_{i^{\prime } } \right ) \left ( l_j\left\| c^{\left ( k+1\right ) } \right\| l_{j^{\prime } } \right ) n^{k-1}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) , \end{array}\ ] ]    where we use a shorthand notation    @xmath118 and radial integral ( [ eq : ss - b ] ) , ( [ eq : ss - c ] ) is defined as in glass and hibbert  @xcite :    @xmath119{c } n^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \\ = \frac{\alpha ^2}4\int_0^\infty \int_0^\infty p_i\left ( r_1\right ) p_j\left ( r_2\right ) \frac{r_2^k}{r_1^{k+3}}\epsilon ( r_1-r_2)p_{i^{\prime } } \left ( r_1\right ) p_{j^{\prime } } \left ( r_2\right ) dr_1dr_{2 , } \end{array}\ ] ]    where @xmath120 is a heaviside step - function ,    @xmath121    the _ spin - other - orbit _
operator @xmath122 itself contains tensorial structure of six different types , summed over @xmath114 ( see gaigalas _ et al _
@xcite ) : @xmath123{c } h^{sso } \equiv \displaystyle { \sum_{k } } \left [ h_{sso}^{(k-1 k 1,101 ) } + h_{sso}^{(k-1 k 1,011)}+ \right . \\
h_{sso}^{(k k 1,101 ) } + h_{sso}^{(k k 1,011)}+ h_{sso}^{(k+1 k 1,101 ) } + h_{sso}^{(k+1 k 1,011 ) } \right ] .
\end{array}\ ] ] their submatrix elements are :    @xmath124{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{soo}^{\left ( k-1k1,\sigma _ 1\sigma _ 21\right ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = 2\cdot 2^{\sigma _ 2}\left\ { \left ( 2k-1\right ) \left ( 2k+1\right ) \right . \\ \times \left . \left ( l_i+l_{i^{\prime } } -k+1\right ) \left ( k - l_i+l_{i^{\prime } } \right ) \left ( k+l_i - l_{i^{\prime } } \right ) \left ( k+l_i+l_{i^{\prime } } + 1\right ) \right\ } ^{1/2 }
\\ \times \left ( k\right ) ^{-1/2}\left ( l_i\left\| c^{\left ( k\right ) } \right\| l_{i^{\prime } } \right ) \left ( l_j\left\| c^{\left ( k\right ) } \right\| l_{j^{\prime } } \right ) n^{k-2}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) , \end{array}\ ] ]    @xmath125{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{soo}^{\left ( kk1,\sigma _ 1\sigma _ 21\right ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = -2\cdot 2^{\sigma _ 2}\left ( 2k+1\right ) ^{1/2}\left ( l_i\left\| c^{\left ( k\right ) } \right\| l_{i^{\prime } } \right ) \\
\times \left ( l_j\left\| c^{\left ( k\right ) } \right\| l_{j^{\prime } } \right ) \left\ { \left ( k\left ( k+1\right ) \right ) ^{-1/2}\left ( l_i\left ( l_i+1\right ) -k\left ( k+1\right ) -l_{i^{\prime } } \left ( l_{i^{\prime } } + 1\right ) \right ) \right .
\\ \times \left\ { \left ( k+1\right ) n^{k-2}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) -kn^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \right\ } \\ \left . -2\left
( k\left ( k+1\right ) \right ) ^{1/2}v^{k-1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \right\ } , \end{array}\ ] ]    @xmath126{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{soo}^{\left ( k+1k1,\sigma _ 1\sigma _ 21\right ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = 2\cdot 2^{\sigma _ 2}\left\ { \left ( 2k+1\right ) \left ( 2k+3\right ) \right . \\ \times \left . \left ( l_i+l_{i^{\prime } } -k\right ) \left ( k - l_i+l_{i^{\prime } } + 1\right ) \left ( k+l_i - l_{i^{\prime } } + 1\right ) \left ( k+l_i+l_{i^{\prime } } + 2\right ) \right\ } ^{1/2 } \\ \times \left ( k+1\right ) ^{-1/2}\left ( l_i\left\| c^{\left ( k\right ) } \right\| l_{i^{\prime } } \right ) \left ( l_j\left\| c^{\left ( k\right ) } \right\| l_{j^{\prime } } \right ) n^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) .
\end{array}\ ] ]    the radial integrals in ( [ eq : soo - b ] ) - ( [ eq : soo - d ] ) are ( see glass and hibbert  @xcite ) :    @xmath127{c } v^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right )
\\ = \frac{\alpha ^2}4\int_0^\infty \int_0^\infty p_i\left ( r_1\right ) p_j\left ( r_2\right ) \frac{r_{<}^{k-1}}{r_{>}^{k+2}}r_2\frac \partial { \partial r_1}p_{i^{\prime } } \left ( r_1\right ) p_{j^{\prime } } \left ( r_2\right ) dr_1dr_2 . \end{array}\ ] ]    the tensorial form of _ orbit - orbit _ operator is ( see eissner _
@xcite )    @xmath128    the sum of submatrix elements of three terms @xmath129 , @xmath130 and @xmath131 is equal to ( see badnell  @xcite ) :    @xmath132{c } \left ( n_i\lambda _ in_j\lambda _ j\left\| h_{oo1}^{(kk0,000 ) } + h_{oo2}^{(kk0,000 ) } + h_{oo4}^{(kk0,000 ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) =    \\
= \left\ { -2k\left ( k+1\right ) \left ( t^{k+1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) -t^{k-1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \right ) \right . \\ + \left ( l_i\left ( l_i+1\right ) -k\left ( k+1\right ) -l_{i^{\prime } } \left ( l_{i^{\prime } } + 1\right ) \right ) \\ \times \left ( u^{k+1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) -u^{k-1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \right )   \\ + \left ( l_j\left ( l_j+1\right ) -k\left ( k+1\right ) -l_{j^{\prime } } \left ( l_{j^{\prime } } + 1\right ) \right ) \\ \times \left ( u^{k+1}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) -u^{k-1}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) \right )
\\ + \frac 12\left ( l_i\left ( l_i+1\right ) -k\left ( k+1\right ) -l_{i^{\prime } } \left ( l_{i^{\prime } } + 1\right ) \right ) \left ( l_j\left ( l_j+1\right ) -k\left ( k+1\right ) -l_{j^{\prime } } \left ( l_{j^{\prime } } + 1\right ) \right )
\\ \times \left [ \frac{k-2}{k\left ( 2k-1\right ) } \left ( n^{k-2}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) + n^{k-2}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) \right ) \right .    \\
-\frac{k+3}{\left ( k+1\right ) \left ( 2k+3\right ) } \left ( n^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) + n^k\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) \right ) \right ] \right\ }
\\ \times \left ( 1-\delta \left ( k,0\right ) \right ) , \end{array}\ ] ]    where @xmath133{c } t^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) = \frac{\alpha ^2}{4\left ( 2k+1\right ) } \\
\times \int_0^\infty \int_0^\infty p_i\left ( r_1\right ) p_j\left ( r_2\right ) \frac{r_<^k}{r_>^{k+1 } } \left(\frac{\partial}{\partial r_1}+\frac{1}{r_1}\right ) p_{i^{\prime } } \left ( r_1\right ) \left(\frac{\partial}{\partial r_2}+\frac{1}{r_2}\right ) p_{j^{\prime } } \left ( r_2\right ) dr_1dr_{2 } , \end{array}\ ] ]    @xmath134{c } u^k\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) = \frac{\alpha ^2}{4\left ( 2k+1\right ) } \\
\times \int_0^\infty \int_0^\infty p_i\left ( r_1\right ) p_j\left ( r_2\right)\left ( ( k-1)\frac{r_2^k}{r_1^{k+2}}\epsilon ( r_1-r_2 ) -(k+2)\frac{r_1^{k-1}}{r_2^{k+1}}\epsilon ( r_2-r_1)\right ) \\ \times p_{i^{\prime } } \left ( r_1\right ) \left(\frac{\partial}{\partial r_2}+\frac{1}{r_2}\right ) p_{j^{\prime } } \left ( r_2\right ) dr_1dr_{2}. \end{array}\ ] ]    the submatrix element of remaining term @xmath135 is :    @xmath136{c } \left ( n_i\lambda _ in_j\lambda _ j\left\|h_{oo3}^{(kk0,000 ) } \right\| n_{i^{\prime } } \lambda _ { i^{\prime } } n_{j^{\prime } } \lambda _ { j^{\prime } } \right ) = 2 \sqrt{2k+1}\frac 1{k(k+1 ) }   \\
\times \left ( \left ( l_i+l_{i^{\prime } } + k+2\right ) \left ( l_i+l_{i^{\prime } } -k\right ) \left ( l_i - l_{i^{\prime } } + k+1\right ) \left ( l_{i^{\prime } } -l_i+k+1\right ) \right .   \\
\times \left . \left ( l_j+l_{j^{\prime } } + k+2\right ) \times \left ( l_j+l_{j^{\prime } } -k\right ) \left ( l_j - l_{j^{\prime } } + k+1\right ) \left ( l_{j^{\prime } } -l_j+k+1\right ) \right ) ^{1/2 }   \\
\times \left ( l_i||c^{\left ( k+1\right ) } ||l_{i^{\prime } } \right ) \left ( l_j||c^{\left ( k+1\right ) } ||l_{j^{\prime } } \right )   \\
\times \left ( n^{k-1}\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) + n^{k-1}\left ( n_jl_jn_il_i , n_{j^{\prime } } l_{j^{\prime } } n_{i^{\prime } } l_{i^{\prime } } \right ) \right ) .
\end{array}\ ] ]    the rest of two - particle breit - pauli operators that we did not investigate so far are the two - body darwin and spin - spin - contact terms .
they do not bring any additional difficulties into the investigation of hamiltonian , but for the sake of completeness of presentation we will discuss them briefly    the _ two - body darwin _
operator @xmath15 ( see for more detail nikitin and rudzikas  @xcite ) , as well as the _ spin - spin - contact _ operator @xmath16 ( see shalit and talmi  @xcite and feneuille  @xcite ) , both have the following tensorial structure :    @xmath137{c } h \equiv \displaystyle { \sum_{k } } h^{(kk0,000)}. \end{array}\ ] ]    these two terms are included into calculation by adding to the radial integral @xmath138 a term    @xmath139    where    @xmath140{c } x\left ( n_il_in_jl_j , n_{i^{\prime } } l_{i^{\prime } } n_{j^{\prime } } l_{j^{\prime } } \right ) \\
= \frac{\alpha ^2}4\int_0^\infty \int_0^\infty p_i\left ( r_1\right ) p_j\left ( r_2\right ) \frac{1}{r_1^{2 } } \delta \left ( r_1-r_2 \right ) p_{i^{\prime } } \left ( r_1\right ) p_{j^{\prime } } \left ( r_2\right ) dr_1dr_{2}. \end{array}\ ] ]    the expression ( [ eq : op - i ] ) has a series of terms , and thus at a first glance seems to be difficult to apply . for this purpose in the next sections we shall discuss in more detail :    * obtaining the expression by wick s theorem . *
the compact written form of all terms , using the extended graphical technique . * obtaining the values of recoupling matrix and of the standard quantities .
we shall also compare the existing methodologies of finding angular parts , showing the advantages and shortcomings of one methodology or another .
this methodology is fully applicable also for the one - particle operator ( see gaigalas and rudzikas  @xcite , gaigalas _ et al _
as this operator does not cause big problems in the atomic physics , we will not stop for its details .
wick s theorem in the second quantization formalism is formulated as follows ( see wick  @xcite ; bogoliubov and shirkov  @xcite ) : if @xmath141 is a product of creation and annihilation operators , then    @xmath142    where @xmath143 represents the normal form of @xmath141 and @xmath144 represents the sum of the normal - ordered terms obtained by making all possible single , double , ... contractions within @xmath141 .
based on the identification in bogoliubov and shirkov  @xcite , the operator is presented in normal form when all of the operators of annihilation included in it are transferred to the right of the creation operators .    usually , wick s theorem is applied when treating complex operators that are represented by a large number of second quantization operators in a non - normal product form . in atomic physics
such operators are used in perturbation theory ( see fetter and waleka  @xcite ; lindgren and morrison  @xcite ; merkelis _ et al _  @xcite ) and in the orthogonal operator method ( see uylings  @xcite ) .
most often the wick s theorem is applied to the products of second quantization operators that are not tensorially coupled ( see lindgren and morrison  @xcite ) . while applying the perturbation theory in an extended model space , two different groups of second quantization operators
are defined ( see for details in lindgren and morrison  @xcite ) .
the second quantization operators acting upon core states belong to one group , whereas the operators acting upon open and excited shells belong to another one .
these two groups are very different in applying wick s theorem to them .    in the first group ,
the creation operators are re - named to annihilation operators and are called the hole absorption operators , while the annihilation operators are re - named to creation operators and are called the hole creation operators .
the creation operators of the second group are called the particle creation operators , while annihilation ones are called the particle absorption operators .
such a division of second quantization operators into two groups is called the particle - hole formalism .
merkelis _ et al _  @xcite have proposed to use the so - called graphical analogue of wick s theorem in perturbation theory ( see gaigalas _ et al _  @xcite , gaigalas  @xcite ) .
it is an efficient tool for obtaining the normal products of second quantization operators in a coupled tensorial form . before applying this theorem ,
particular second quantization terms are in a normal product in coupled form . in addition , this theorem is applied in the particle - hole formalism , too ( see merkelis _
@xcite , merkelis _ et al _
@xcite ) .    in all the cases mentioned above , the wick s theorem is applied for the most general case of operators , i.e. the shells that are acted upon are not detailed .
but , in the case of the extended model space , the group that the second quantization operators belong to , depending on the electronic structure of atom or ion under investigation , is defined .
wick s theorem is not applied in investigations of ordinary physical operators .
@xcite proposed a new version of the wick s theorem application , where the optimal tensorial expression of any two - particle operator is easily obtained .
the specificity of wick s theorem application in this case lies in applying it only when the shells that are acted upon by the secondary quantization operators are known , i.e. it is applied for each term @xmath145 separately .
such an interpretation of wick s theorem bears similarity with the particle - hole formalism . the only difference is that in this case the second quantization operators are differentiated formally not on the basis of structure of atom under investigation , but on the basis of shells acted upon .
this is done in the following way .
the second quantization operators acting upon a shell with a lowest index are attributed to the first group .
those acting upon a shell with a next - lowest index are attributed to the second group , etc . in the most general case
we have four distinct groups .    for example
, suppose we have a two - particle operator @xmath146 ,    @xmath147    ( where @xmath148 , @xmath149 , @xmath150 ) , the matrix element of which between the functions    @xmath151 and @xmath152 we must obtain .
then the operators acting upon the second shell are attributed to the first group , the ones acting upon the third shell - to the second , and the ones acting upon the fifth - to the third group .
assuming that all the operators from first group are the creation ones , and the rest are annihilation operators , we apply the wick s theorem .
thus we obtain that the operators from the first group are one beside another , and all are positioned after the operators from the first group .
thus we obtain :    @xmath153    after that , we apply wick s theorem assuming that the operators from the first and the second groups are creation ones , and the rest of them are annihilation operators .
thus we obtain that the operators from the second group are one beside another , and all are positioned after the first group operators :    @xmath154    if in the product that we investigate there are operators of second quantization acting upon four distinct shells , then we apply wick s theorem once again , assuming that operators from the first , second and third groups are creation ones , and from the fourth group - annihilation ones .
in this case the wick s theorem is applied to the second quantization operators in uncoupled form .    from ( [ eq : op - b ] ) we see that in second quantization a two - particle operator is written as a sum , where parameters @xmath44 , @xmath155 , @xmath156 , @xmath157 run over all possible arrays of quantum numbers .
so , the greater the number of open shells in bra and ket functions , the greater the number of terms @xmath145 in the expression of two - particle operator .
it is obvious that all these terms must be systematized in order to obtain in general case the most efficient tensorial expression of a two - particle operator , in the way described above .
.distributions of shells , upon which the second quantization operators are acting , that appear in the submatrix elements of any two - particle operator , when bra and ket functions have @xmath22 open shells ( see gaigalas _
@xcite ) [ cols=">,^,^,^,^,^ " , ]     in the work by gaigalas _
et al _  @xcite it is chosen an optimal number of distributions , which is necessary to obtain the matrix elements of any two - particle operator , when the bra and ket functions consist of arbitrary number of shells .
this is presented in table 1 .
we point out that for distributions 2 - 5 and 19 - 42 the shells sequence numbers @xmath158 , @xmath159 , @xmath160 , @xmath161 ( in bra and ket functions of a submatrix element ) satisfy the condition @xmath162 , while for distributions 6 - 18 no conditions upon @xmath158 , @xmath159 , @xmath160 , @xmath161 are imposed ( this permits to reduce the number of distributions ) . for distributions 19 - 42
this condition is imposed only for obtaining simple analytical expressions for the recoupling matrices @xmath163 .
this will be discussed in more detail in the next section .
so , in the way that is described earlier , the wick s theorem is applied , assuming that the second quantization operators acting upon shells @xmath164 , @xmath165 , @xmath160 and @xmath161 belong to different groups .    in the next section
we discuss the way to obtain irreducible tensorial form of these distributions .
in addition , the arguments will be given in evidence of superiority of the obtained tensorial expressions against other expressions known in the literature .
the methodology presented in this section demonstrates the way to obtain optimal arrangement of the second quantization operators , for any two - particle operator .
it can be applied without restrictions for obtaining the optimal tensorial form of two - particle terms of orthogonal operators and of perturbation theory operators , too .
the graphical technique of angular momentum is widely used in the atomic physics : see yutsis _
@xcite , jucys and bandzaitis  @xcite , brink and satcher  @xcite , el - baz  @xcite .
it is applied efficiently both in the coordinate representation ( see for example jucys and bandzaitis  @xcite ) , and in the second quantization formalism ( gaigalas _ et al _
the use of it allows one to obtain the analytical expressions for the recoupling matrices conveniently ( see for example kaniauskas and rudzikas @xcite ) , to investigate the tensorial products of operators ( see for example jucys _
@xcite ) , and to seek for the matrix elements of operators ( see for example huang and starace  @xcite ) . gaigalas and merkelis @xcite have proposed a graphical way to obtain the values of matrix elements when the operator is a many - particle ( one- , two- , three- , etc . )
one and has irreducible tensorial form .
for example , when a two - particle operator is considered , it has the form ( [ eq : op - c ] ) or ( [ eq : op - e ] ) . the matrix elements in this methodology are expressed not only in the terms of coefficients of fractional parentage or reduced coefficients of fractional parentage , but in terms of standard quantities @xmath166 and @xmath167 , too .
@xcite have proposed to calculate the matrix elements by using the tensorial expressions for such two - particle operator , that take full advantage of racah algebra . in this case
the tensorial form of operator depends on the shells that the operator acts upon ( distributions 1 - 42 from table  1 ) .
this is the difference of this methodology from other .
it is most convenient to obtain the tensorial expressions for 42 distributions graphically , using the generalized graphical technique by gaigalas _
et al _  @xcite .
in such a case not only the similarities between different distributions are easily seen , but also the compact graphical representation of the obtained expressions is possible .
we will stop for details on this in the present section .
a two - particle operator may be represented graphically by a feynman - goldstone diagram @xmath168 from figure  [ fig : dd - a ] ( lindgren and morrison  @xcite ) . as it is shown in the paper bolotin _
@xcite , the feynman - goldstone diagrams are topologically equivalent to the angular momentum graphs . due to that , an irreducible tensorial form for every feynman - goldstone diagram may be obtained ( see merkelis _
the graph @xmath169 is the angular momentum graph corresponding to the diagram @xmath168 .
so the two - particle operator will be written down as follows :    @xmath170{c } \widehat{g } ( ij , i^{\prime } j^{\prime } ) = d_{1}= - \frac 12 \displaystyle { \sum_{m_{\lambda _ i}m_{\lambda _ j } m_{\lambda ^{\prime } _
i}m_{\lambda ^{\prime } _ j } } } \displaystyle { \sum_p } \left [ \kappa _ { 1},\kappa _ { 2},\sigma _ { 1},\sigma _ { 2}\right ] ^{-1/2 } \times \\ \times \left ( n_i\lambda _ in_j\lambda _ j||g^{\left ( \kappa _ 1\kappa _ 2k,\sigma _ 1\sigma _ 2k\right ) } ||n_i^{\prime } \lambda _ i^{\prime } n_j^{\prime } \lambda _ j^{\prime }
\right ) ~d_{2}~~ a^{\left ( \lambda _ i\right ) } _ { m_{\lambda _ i } } a^{\left ( \lambda _ j\right ) } _ { m_{\lambda _ j } } \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } _ { m_{\lambda ^{\prime } _ i } } \stackrel{\sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } _ { m_{\lambda ^{\prime } _ j } } , \end{array}\ ] ]    ( 148,180 ) ( 10,140)(0,1)10 ( 10,150)(0,1)2 ( 5,148)(0,0)[t]@xmath171 ( 20,140)(0,1)10 ( 20,150)(0,1)2 ( 25,148)(0,0)[t]@xmath172 ( 11,140)(2,0)5(2,1)[t ] ( 15,143)(0,0)@xmath173 ( 10,130)(0,1)10 ( 10,130)(0,1)3 ( 5,133)(0,0)@xmath174 ( 10,128)(0,1)3 ( 20,130)(0,1)10 ( 20,130)(0,1)3 ( 25,133)(0,0)@xmath175 ( 20,128)(0,1)3 ( 15,120)(0,0)@xmath176 ( 45,140)(0,1)12 ( 45,140 ) ( 40,148)(0,0)[t]@xmath177 ( 75,140)(0,12)12 ( 75,152)(0,-1)3 ( 75,140 ) ( 81,148)(0,0)[t]@xmath178 ( 45,140.6)(6,0)6 ( 45,140.3)(6,0)6 ( 45,140)(20,0)30 ( 45,139.7)(6,0)6 ( 45,139.4)(6,0)6 ( 52,143)(0,0)@xmath179 ( 69,140.6)(6,0)6 ( 69,140.3)(6,0)6 ( 69,139.7)(6,0)6 ( 69,139.4)(6,0)6 ( 68,143)(0,0)@xmath180 ( 60,140 ) ( 60,145)(0,0)@xmath181 ( 60.6,135)(0,1)5 ( 60.3,135)(0,1)5 ( 60,135)(0,1)5 ( 59.7,135)(0,1)5 ( 59.4,135)(0,1)5 ( 60,132)(0,0)@xmath182 ( 45,130)(0,1)10 ( 39,134)(0,0)@xmath183 ( 40,140)(0,0)@xmath184 ( 45,130)(0,-1)3 ( 75,130)(0,1)10 ( 81,134)(0,0)@xmath185 ( 80,140)(0,0)@xmath181 ( 75,127)(0,1)3 ( 60,120)(0,0)@xmath186 ( 100,136)(0,1)12 ( 100,147)(5,5)[b ] ( 94,144)(0,0)[t]@xmath171 ( 100,120)(0,0)@xmath187 ( 120,126)(0,1)10 ( 114,130)(0,0)@xmath172 ( 120,123)(0,1)3 ( 120,125)(5,5)[t ] ( 120,120)(0,0)@xmath188 ( 10,90)(0,1)10 ( 10,100)(0,1)2 ( 5,98)(0,0)[t]@xmath171 ( 20,90)(0,1)10 ( 20,100)(0,1)2 ( 25,98)(0,0)[t]@xmath172 ( 10,90)(1,0)10 ( 15,93)(0,0)@xmath173 ( 10,80)(0,1)10 ( 10,80)(0,1)3 ( 5,83)(0,0)@xmath174 ( 10,78)(0,1)3 ( 20,80)(0,1)10 ( 20,80)(0,1)3 ( 25,83)(0,0)@xmath175 ( 20,78)(0,1)3 ( 15,70)(0,0)@xmath189 ( 80,90)(0,1)12 ( 80,101)(5,5)[b ] ( 80,90 ) ( 75,98)(0,0)[t]@xmath171 ( 110,90)(0,12)12 ( 110,102)(0,-1)3 ( 110,101)(5,5)[b ] ( 110,90 ) ( 115,98)(0,0)[t]@xmath174 ( 80,90.6)(6,0)6 ( 80,90.3)(6,0)6 ( 80,90)(20,0)30 ( 80,89.7)(6,0)6 ( 80,89.4)(6,0)6 ( 87,93)(0,0)@xmath190 ( 104,90.6)(6,0)6 ( 104,90.3)(6,0)6 ( 104,89.7)(6,0)6 ( 104,89.4)(6,0)6 ( 103,93)(0,0)@xmath191 ( 95,90 ) ( 95,95)(0,0)@xmath181 ( 95.6,85)(0,1)5 ( 95.3,85)(0,1)5 ( 95,85)(0,1)5 ( 94.7,85)(0,1)5 ( 94.4,85)(0,1)5 ( 95,82)(0,0)@xmath182 ( 80,80)(0,1)10 ( 74,84)(0,0)@xmath172 ( 75,90)(0,0)@xmath184 ( 80,80)(0,-1)3 ( 80,79)(5,5)[t ] ( 110,80)(0,1)10 ( 116,84)(0,0)@xmath175 ( 115,90)(0,0)@xmath181 ( 110,77)(0,1)3 ( 110,79)(5,5)[t ] ( 95,70)(0,0)@xmath192 ( 10,40)(0,1)10 ( 10,50)(0,1)2 ( 5,48)(0,0)[t]@xmath171 ( 10,40)(1,0)10 ( 15,37)(0,0)@xmath173 ( 16,48)(0,0)@xmath193 ( 16,45)(-1,0)2 ( 20,30)(0,1)10 ( 15,40)(10,10)[t ] ( 20,30)(0,1)3 ( 26,33)(0,0)@xmath175 ( 20,28)(0,1)3 ( 15,20)(0,0)@xmath194 ( 45,40)(0,1)12 ( 45,51)(5,5)[b ] ( 45,40 ) ( 40,48)(0,0)[t]@xmath171 ( 75,40)(0,1)12 ( 75,51)(5,5)[b ] ( 75,40 ) ( 80,48)(0,0)[t]@xmath172 ( 45,40.6)(6,0)6 ( 45,40.3)(6,0)6 ( 45,40)(20,0)30 ( 45,39.7)(6,0)6 ( 45,39.4)(6,0)6 ( 52,43)(0,0)@xmath179 ( 69,40.6)(6,0)6 ( 69,40.3)(6,0)6 ( 69,39.7)(6,0)6 ( 69,39.4)(6,0)6 ( 68,43)(0,0)@xmath195 ( 60,40 ) ( 60,45)(0,0)@xmath181 ( 60.6,35)(0,1)5 ( 60.3,35)(0,1)5 ( 60,35)(0,1)5 ( 59.7,35)(0,1)5 ( 59.4,35)(0,1)5 ( 60,32)(0,0)@xmath182 ( 45,30)(0,1)10 ( 39,34)(0,0)@xmath174 ( 40,40)(0,0)@xmath184 ( 45,27)(0,1)3 ( 45,29)(5,5)[t ] ( 75,30)(0,1)10 ( 81,34)(0,0)@xmath175 ( 80,40)(0,0)@xmath181 ( 75,27)(0,1)3 ( 75,29)(5,5)[t ] ( 60,20)(0,0)@xmath196 ( 110,42 ) ( 110,37)(0,0)@xmath184 ( 110.6,42)(0,1)5 ( 110.3,42)(0,1)5 ( 110,42)(0,1)5 ( 109.7,42)(0,1)5 ( 109.4,42)(0,1)5 ( 115,45)(0,0)@xmath182 ( 100,32)(1,1)10 ( 100,27)(0,1)5 ( 94,34)(0,0)@xmath171 ( 100,30)(0,-1)3 ( 100,29)(5,5)[t ] ( 120,32)(-1,1)10 ( 120,27)(0,1)5 ( 126,34)(0,0)@xmath175 ( 120,27)(0,1)3 ( 120,29)(5,5)[t ] ( 110,20)(0,0)@xmath197    it must be noted that in expression ( [ eq : pas - a ] ) the projection @xmath198 of @xmath199 , as well as that of the momentum line @xmath200 in graph @xmath169 are the same .
this is also to be said about the remaining operators of second quantization and the three open lines af graph @xmath169 .
as it has been mentiond in section 2 , this operator has two tensorial forms , ( [ eq : op - c ] ) and ( [ eq : op - d ] ) .
these may be represented graphically , since the creation operator @xmath201 , as well as operator @xmath202 respectivelly , are graphically denoted by diagrams @xmath203 and @xmath204 .
the first form ( [ eq : op - c ] ) of two - particle operator @xmath205 is represented as :    @xmath206{c } g_i(ij , i^{\prime } j^{\prime } ) = d_1=-\frac 12 \displaystyle { \sum_{\kappa _ { 12}\kappa _ { 12}^{^{\prime } } \sigma _ { 12}\sigma _ { 12}^{^{\prime } } } } \displaystyle { \sum_p}\left ( -1\right ) ^{k - p}\left [ \kappa _ 1,\kappa _ 2,\sigma _ 1,\sigma _ 2\right ] ^{1/2}\times \\ \times \left ( n_i\lambda _ in_j\lambda _ j||g^{\left ( \kappa _ 1\kappa _ 2k,\sigma _ 1\sigma _ 2k\right ) } ||n_i^{\prime } \lambda _ i^{\prime } n_j^{\prime } \lambda _ j^{\prime } \right ) \left\ { \begin{array}{ccc } l_i^{\prime } & l_j^{\prime } & \kappa _ { 12}^{\prime } \\ \kappa _ 1 & \kappa _ 2 & k \\ l_i & l_j & \kappa _ { 12 } \end{array } \right\ } \times \\ \times \left\ { \begin{array}{ccc } s & s & \sigma _ { 12}^{\prime } \\
\sigma _ 1 & \sigma _ 2 & k \\ s & s & \sigma _ { 12 } \end{array } \right\ } d_6 \end{array}\ ] ]    whereas the second ( [ eq : op - d ] ) :    @xmath207{c } g_{ii}(ij , i^{\prime } j^{\prime } ) = d_5+d_7= \\
= \frac 12 \displaystyle { \sum_p}\left ( -1\right ) ^{k - p}\left ( n_i\lambda _ in_j\lambda _ j||g^{\left ( \kappa _ 1\kappa _ 2k,\sigma _ 1\sigma _ 2k\right ) } ||n_i^{\prime } \lambda _ i^{\prime } n_j^{\prime } \lambda _ j^{\prime } \right ) \times \\ \times \{\left [ \kappa _ 1,\kappa _ 2,\sigma _ 1,\sigma _ 2\right ] ^{-1/2}d_8-\left ( -1\right ) ^{l_i+l_j^{\prime } } \left\ { \begin{array}{ccc } \kappa _ 1 & \kappa _ 2 & k \\
l_j^{\prime } & l_i & l_j \end{array } \right\ } \times \\ \times \left\ { \begin{array}{ccc } \sigma _ 1 & \sigma _ 2 & k \\ s & s & s \end{array } \right\ } \delta \left ( n_jl_j , n_i^{\prime } l_i^{\prime } \right ) d_9\}. \end{array}\ ] ]    we emphasize here that the winding line of interaction in the feynman - goldstone diagram corresponds to the operators of second quantization in the normal order ( figure  [ fig : dd - a ] , @xmath208 ) .
whereas the dotted interaction line indicates that the second quantization operators are ordered in pairs of creation - annihilation operators . in the latter case
first comes the pair on the left side of a feynman - goldstone diagram ( figure  [ fig : dd - a ] , @xmath209 ) .
such a notation of two kinds for an interaction line is meaningful only for two - particle ( or more ) operators , since for any one - particle operator both the winding and dotted lines correspond to the same order of creation and annihilation operators .    from expressions ( [ eq : pas - b ] ) , ( [ eq : pas - c ] ) we see that the two - particle operator in the first form is represented by one feynman - goldstone diagram @xmath208 , whereas in the second - by two diagrams @xmath209 and @xmath210 .
the diagrams , corresponding to tensorial product , have the following algebraic expressions :    @xmath211 ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times \left [ \stackrel{\sim } { a}^{\left ( \lambda _ i^{\prime } \right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \right ] _ { p - p}^{\left ( kk\right ) } , \ ] ]    @xmath212 ^{\left ( \kappa _ 1\sigma _ 1\right ) } \times \left [ a^{\left ( \lambda _ j\right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ j^{\prime } \right ) } \right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } \right ] _ { p - p}^{\left ( kk\right ) } \ ] ]    @xmath213 _ { p - p}^{\left ( kk\right ) } \ ] ]    the positions of the second quantization operators in the diagrams @xmath214 , @xmath215 and @xmath216 define their order in tensorial products : the first place in tensorial product occupies the upper right second quantization operator , the second - lower right , after them the upper left and lower left operators follow .
the angular momenta diagram defines their coupling scheme into tensorial product . for more detail
see gaigalas and rudzikas  @xcite .
as it has been mentioned earlier , these two forms do not always take full advantage of the racah algebra ( see gaigalas and rudzikas  @xcite ) . the expression ( [ eq : op - i ] ) has no such shortcomings .
now we will demonstrate the way to obtain graphically a tensorial expressions for particular distributions 1 - 42 from table  1 .
( 148,120 ) ( 10,105)(0,1)10 ( 10,115)(0,1)2 ( 5,113)(0,0)[t]@xmath217 ( 20,105)(0,1)10 ( 20,115)(0,1)2 ( 25,113)(0,0)[t]@xmath218 ( 11,105)(2,0)5(2,1)[t ] ( 15,110)(0,0)@xmath173 ( 10,95)(0,1)10 ( 10,95)(0,1)3 ( 5,98)(0,0)@xmath219 ( 10,93)(0,1)3 ( 20,95)(0,1)10 ( 20,95)(0,1)3 ( 25,98)(0,0)@xmath217 ( 20,93)(0,1)3 ( 15,85)(0,0)@xmath220 ( 70,105)(0,1)12 ( 70,105 ) ( 65,113)(0,0)[t]@xmath221 ( 100,105)(0,12)12 ( 100,117)(0,-1)3 ( 100,105 ) ( 106,113)(0,0)[t]@xmath222 ( 70,105.6)(6,0)6 ( 70,105.3)(6,0)6 ( 70,105)(20,0)30 ( 70,104.7)(6,0)6 ( 70,104.4)(6,0)6 ( 77,108)(0,0)@xmath179 ( 94,105.6)(6,0)6 ( 94,105.3)(6,0)6 ( 94,104.7)(6,0)6 ( 69,104.4)(6,0)6 ( 93,108)(0,0)@xmath180 ( 85,105 ) ( 85,110)(0,0)@xmath181 ( 85.6,100)(0,1)5 ( 85.3,100)(0,1)5 ( 85,100)(0,1)5 ( 84.7,100)(0,1)5 ( 84.4,100)(0,1)5 ( 85,97)(0,0)@xmath182 ( 70,95)(0,1)10 ( 65,99)(0,0)@xmath223 ( 65,105)(0,0)@xmath184 ( 70,95)(0,-1)3 ( 100,95)(0,1)10 ( 106,99)(0,0)@xmath224 ( 105,105)(0,0)@xmath181 ( 100,92)(0,1)3 ( 85,85)(0,0)@xmath225 ( 10,50)(0,1)12 ( 10,61)(5,5)[b ] ( 10,50 ) ( 5,58)(0,0)[t]@xmath226 ( 40,50)(0,12)12 ( 40,62)(0,-1)3 ( 40,61)(5,5)[b ] ( 40,50 ) ( 46,58)(0,0)[t]@xmath227 ( 10,50.6)(6,0)6 ( 10,50.3)(6,0)6 ( 10,50)(20,0)30 ( 10,49.7)(6,0)6 ( 10,49.4)(6,0)6 ( 17,53)(0,0)@xmath190 ( 34,50.6)(6,0)6 ( 34,50.3)(6,0)6 ( 34,49.7)(6,0)6 ( 34,49.4)(6,0)6 ( 33,53)(0,0)@xmath191 ( 25,50 ) ( 25,55)(0,0)@xmath181 ( 25.6,45)(0,1)5 ( 25.3,45)(0,1)5 ( 25,45)(0,1)5 ( 24.7,45)(0,1)5 ( 24.4,45)(0,1)5 ( 25,42)(0,0)@xmath182 ( 10,40)(0,1)10 ( 4,44)(0,0)@xmath226 ( 5,50)(0,0)@xmath184 ( 10,40)(0,-1)3 ( 10,39)(5,5)[t ] ( 40,40)(0,1)10 ( 46,44)(0,0)@xmath228 ( 45,50)(0,0)@xmath181 ( 40,37)(0,1)3 ( 40,39)(5,5)[t ] ( 25,20)(0,0)@xmath229 ( 70,63 ) ( 100,63 ) ( 70,63)(1,-1)30 ( 100,63)(-1,-1)30 ( 70,63.6)(6,0)6 ( 70,63.3)(6,0)6 ( 70,63)(20,0)30 ( 70,62.7)(6,0)6 ( 70,62.4)(6,0)6 ( 77,66)(0,0)@xmath179 ( 94,63.6)(6,0)6 ( 94,63.3)(6,0)6 ( 94,62.7)(6,0)6 ( 94,62.4)(6,0)6 ( 93,66)(0,0)@xmath180 ( 85,63 ) ( 85,68)(0,0)@xmath184 ( 85.6,63)(0,-1)30 ( 85.3,63)(0,-1)30 ( 85,63)(0,-1)30 ( 84.7,63)(0,-1)30 ( 84.4,63)(0,-1)30 ( 81,59)(0,0)@xmath182 ( 77,55)(0,0)[t]@xmath230 ( 65,48)(0,0)[t]@xmath221 ( 70,68)(0,0)@xmath181 ( 70,63)(0,-1)30 ( 95,55)(0,0)[t]@xmath221 ( 106,48)(0,0)[t]@xmath222 ( 100,68)(0,0)@xmath181 ( 100,63)(0,-1)30 ( 85,33 ) ( 85,28)(0,0)@xmath184 ( 70,28)(0,0)@xmath181 ( 70,33 ) ( 70,33.6)(6,0)6 ( 70,33.3)(6,0)6 ( 70,33)(20,0)30 ( 70,32.7)(6,0)6 ( 70,32.4)(6,0)6 ( 77,30)(0,0)@xmath190 ( 100,28)(0,0)@xmath184 ( 100,33 ) ( 94,33.6)(6,0)6 ( 94,33.3)(6,0)6 ( 94,32.7)(6,0)6 ( 94,32.4)(6,0)6 ( 93,30)(0,0)@xmath191 ( 85,20)(0,0)@xmath231    we take the distribution @xmath160 , @xmath159 , @xmath158 , @xmath160 ( 13 form table 1 ) as an example for investigation
. then the feynman - goldstone diagram of operator @xmath232 is @xmath233 , the angular momentum graph is @xmath234 , and the second quantization operators are in the following order : @xmath235 applying the wick s theorem as described in section 3 , and assuming that the operators acting upon shell @xmath160 belong to the first group , the once acting upon @xmath159 belong to the second , and the ones acting upon @xmath158 belong to the third group , we obtain the following order of operators : @xmath236 using the generalized graphical technique of angular momentum by gaigalas _
@xcite , we couple the operators of second quantization into tensorial product @xmath237 ( see figure  [ fig : dd - c ] ) :    @xmath238 ^{\left ( \kappa _ { 12}\sigma _ { 12}\right ) } \times \left [ \stackrel{\sim } { a}^{\left ( \lambda _ { \beta } \right ) } \times \stackrel{\sim } { a}^{\left ( \lambda _ { \alpha}\right ) } \right ] ^{\left ( \kappa _ { 12}^{\prime } \sigma _ { 12}^{\prime } \right ) } \right ] _ { p - p}^{\left ( kk\right ) } .\ ] ]    in the course of obtaining @xmath237 graphically , a recoupling matrix @xmath239 appears , whose analytical expression is readily obtained by using the graphical technique of jucys and bandzaitis  @xcite .
all the needed expressions are obtained in the same way .
( 148,105 ) ( 20,90)(0,1)10 ( 20,100)(0,1)2 ( 15,98)(0,0)[t]@xmath240 ( 30,90)(0,1)10 ( 30,100)(0,1)2 ( 35,98)(0,0)[t]@xmath241 ( 20,90)(1,0)10 ( 25,93)(0,0)@xmath173 ( 20,80)(0,1)10 ( 20,80)(0,1)3 ( 15,83)(0,0)@xmath242 ( 20,78)(0,1)3 ( 30,80)(0,1)10 ( 30,80)(0,1)3 ( 35,83)(0,0)@xmath243 ( 30,78)(0,1)3 ( 25,70)(0,0)@xmath244 ( 80,90)(0,1)12 ( 80,101)(5,5)[b ] ( 80,90 ) ( 74,98)(0,0)[t]@xmath240 ( 110,90)(0,1)12 ( 110,101)(5,5)[b ] ( 110,90 ) ( 115,98)(0,0)[t]@xmath241 ( 80,90.6)(6,0)6 ( 80,90.3)(6,0)6 ( 80,90)(20,0)30 ( 80,89.7)(6,0)6 ( 80,89.4)(6,0)6 ( 87,93)(0,0)@xmath179 ( 104,90.6)(6,0)6 ( 104,90.3)(6,0)6 ( 104,89.7)(6,0)6 ( 104,89.4)(6,0)6 ( 103,93)(0,0)@xmath195 ( 95,90 ) ( 95,95)(0,0)@xmath181 ( 95.6,85)(0,1)5 ( 95.3,85)(0,1)5 ( 95,85)(0,1)5 ( 94.7,85)(0,1)5 ( 94.4,85)(0,1)5 ( 95,82)(0,0)@xmath182 ( 80,80)(0,1)10 ( 74,84)(0,0)@xmath242 ( 75,90)(0,0)@xmath184 ( 80,77)(0,1)3 ( 80,79)(5,5)[t ] ( 110,80)(0,1)10 ( 116,84)(0,0)@xmath243 ( 115,90)(0,0)@xmath181 ( 110,77)(0,1)3 ( 110,79)(5,5)[t ] ( 95,70)(0,0)@xmath245 ( 10,40)(0,1)12 ( 10,51)(5,5)[b ] ( 10,40 ) ( 4,48)(0,0)[t]@xmath241 ( 40,40)(0,1)12 ( 40,51)(5,5)[b ] ( 40,40 ) ( 45,48)(0,0)[t]@xmath240 ( 10,40.6)(6,0)6 ( 10,40.3)(6,0)6 ( 10,40)(20,0)30 ( 10,39.7)(6,0)6 ( 10,39.4)(6,0)6 ( 17,43)(0,0)@xmath179 ( 34,40.6)(6,0)6 ( 34,40.3)(6,0)6 ( 34,39.7)(6,0)6 ( 34,39.4)(6,0)6 ( 33,43)(0,0)@xmath195 ( 25,40 ) ( 25,45)(0,0)@xmath184 ( 28,37)(0,0)@xmath246 ( 25.6,35)(0,1)5 ( 25.3,35)(0,1)5 ( 25,35)(0,1)5 ( 24.7,35)(0,1)5 ( 24.4,35)(0,1)5 ( 25,32)(0,0)@xmath182 ( 10,30)(0,1)10 ( 4,34)(0,0)@xmath243 ( 5,40)(0,0)@xmath184 ( 10,27)(0,1)3 ( 10,29)(5,5)[t ] ( 40,30)(0,1)10 ( 46,34)(0,0)@xmath242 ( 45,40)(0,0)@xmath181 ( 40,27)(0,1)3 ( 40,29)(5,5)[t ] ( 25,20)(0,0)@xmath247 ( 80,40)(0,1)12 ( 80,51)(5,5)[b ] ( 80,40 ) ( 74,48)(0,0)[t]@xmath241 ( 110,40)(0,1)12 ( 110,51)(5,5)[b ] ( 110,40 ) ( 115,48)(0,0)[t]@xmath240 ( 80,40.6)(6,0)6 ( 80,40.3)(6,0)6 ( 80,40)(20,0)30 ( 80,39.7)(6,0)6 ( 80,39.4)(6,0)6 ( 87,43)(0,0)@xmath179 ( 104,40.6)(6,0)6 ( 104,40.3)(6,0)6 ( 104,39.7)(6,0)6 ( 104,39.4)(6,0)6 ( 103,43)(0,0)@xmath195 ( 95,40 ) ( 95,45)(0,0)@xmath181 ( 98,37)(0,0)@xmath246 ( 95.6,35)(0,1)5 ( 95.3,35)(0,1)5 ( 95,35)(0,1)5 ( 94.7,35)(0,1)5 ( 94.4,35)(0,1)5 ( 95,32)(0,0)@xmath182 ( 80,30)(0,1)10 ( 74,34)(0,0)@xmath243 ( 75,40)(0,0)@xmath184 ( 80,27)(0,1)3 ( 80,29)(5,5)[t ] ( 110,30)(0,1)10 ( 116,34)(0,0)@xmath242 ( 115,40)(0,0)@xmath181 ( 110,27)(0,1)3 ( 110,29)(5,5)[t ] ( 95,20)(0,0)@xmath248    in obtaining these expressions , as well as in representing them graphically , it is very convenient to use the rule of changing the sign of a node , existing in the graphical technique of angular momentum ( see jucys and bandzaitis  @xcite )
. therefore now we will treat an example of using such a rule .
suppose , we have the following correspondence between diagrams ( figure  [ fig : dd - c ] ) :    @xmath249    in which the second quantization operators are in the order @xmath250
@xmath251 @xmath252 @xmath253 .
our goal is to obtain the diagram corresponding to the order @xmath254
@xmath253 @xmath255 @xmath251 .
bearing in mind that the second quantization operators anticommute with each other and they all act on different electronic shells and we are not changing the order of their coupling into tensorial product , we arrive at    @xmath256    let us also discuss another situation : we have defined the disposition of the operators and we want to couple them into certain tensorial product .
suppose that we want to represent graphically the following tensorial product :    @xmath257 ^{\left ( \kappa _ 1\sigma _ 1\right ) } \times \left [ a^{(\lambda _ 3)}\times \tilde a^{(\lambda
_ 4)}\right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } \right ] ^{\left ( \kappa \sigma \right ) } .\ ] ]    for this purpose we have to rearrange the generalized clebsch - gordan coefficient , which is defining the scheme of coupling of the operators into the tensorial product .
it is easy to notice that this coefficient will properly define the tensorial product , if we change the sign of the vertex
@xmath258  in diagram @xmath259 :    @xmath260    the procedures described are fairly simple , however , they are sufficient for the majority of cases .
the more complete description of this generalized graphical approach may be found in gaigalas _
@xcite , gaigalas  @xcite , gaigalas and merkelis  @xcite .
all the analytical expressions for distributions 1 - 42 from table 1 are presented in the paper gaigalas _ et al _
they are written down using the generalized graphical methodology of angular momentum , and the vortex sign change rule , which was discussed in this section . as a consequence of that
, the analytical expressions for 42 terms may be written down via 6 different expressions .
this , undoubtedly , facilitates a lot the implementation of methodology proposed in gaigalas _
et al _  @xcite .
in this section we will discuss several ways to obtain matrix elements of a two - particle operator . as it was mentioned earlier , up to now the fano calculation scheme  @xcite is the most popular one .
its general expression when a two - particle operator acts upon different shells is presented in ( [ eq : ma - c ] ) .
the general expression for a matrix element in other cases is similar .
for example , when the operator acts only upon one shell , we have    @xmath261{c } ( \psi _ u^{bra}\left ( ls\right ) || \widehat{g}^{\left ( \kappa _ 1\kappa _ 2k,\sigma_1\sigma _ 2k\right ) }    \sim \displaystyle { \sum_{n_{i } \lambda _ { i } } } \left ( -1\right ) ^\delta n_{i } \left ( n_{i}-1 \right ) \times \\ \times \displaystyle { \sum _ { \ { t \ } } } \left ( l_{i}^{n_{i}}\;\alpha_{i } l_{i}s_{i}||l_{i}^{n_{i}-1}\;\left ( \alpha_{i}^{\prime \prime } l^{\prime \prime}_{i } s^{\prime \prime}_{i } \right ) , l_{i } \right ) \left ( l_{i}^{n_{i}-1}\;\alpha_{i}^{\prime \prime } l_{i}^{\prime \prime}s_{i}^{\prime \prime }    \alpha_{i}^{\prime \prime \prime } l_{i}^{\prime \prime \prime } s^{\prime \prime \prime}_{i } \right ) , l_{i } \right ) \times   \\
\times \left ( l_{i}^{n_{i}^{\prime } } \ ; \alpha_{i}^{\prime \prime \prime } l_{i}^{\prime \prime \prime}s_{i}^{\prime \prime \prime }    \alpha_{i}^{\prime v } l^{\prime v}_{i } s^{\prime v}_{i } \right ) , l_{i } \right ) \left ( l_{i}^{n_{i}^{\prime } -1}\;\alpha_{i}^{\prime v } l_{i}^{\prime v}s_{i}^{\prime v}||l_{i}^{n_{i}^{\prime } -2}\;\left ( \alpha_{i}^{\prime } l^{\prime}_{i } s^{\prime } _ { i } \right ) , l_{i } \right ) \times   \\
\times r_{d } \left ( \lambda _ i,\lambda _ i,\lambda _ i,\lambda _ i , \kappa_1 , \kappa_2 , \kappa , \sigma_1 , \sigma_2 , \sigma , \lambda ^{bra } , \lambda ^{ket } \right )   \times \\ \times \left ( n_i\lambda _ in_i\lambda _ i|| g^{\left ( \kappa _ 1\kappa _ 2k,\sigma_1\sigma _ 2k\right ) }    \end{array}\ ] ]    as we see here , in contrast to ( [ eq : ma - c ] ) , the summation is performed only over one array @xmath38 , @xmath39 of quantum numbers , because the operator acts only upon one shell .
but here a summation over @xmath262 occurs , though , which indicates the summation over arrays of intermediate terms @xmath263 , @xmath264 , @xmath265 .    remembering the relationship between a coefficient of fractional parentage and a reduced matrix element of a second quantization operator ( see pakauskas  _ et al _ @xcite , rudzikas and kaniauskas  @xcite ) : @xmath266{c } \left ( l^{n}\;\alpha ls||a^{(ls)}|| l^{n-1}\ ; \alpha^{\prime } l^{\prime } s^{\prime } \right ) = \\ = ( -1)^{n+(n+1)\varphi ( n ) } \sqrt{n[l , s ] } \left ( l^{n}\;\alpha ls||l^{n-1}\;\left ( \alpha^{\prime } l^{\prime } s^{\prime } \right ) , l \right ) \end{array}\ ] ] we see that the racah algebra in expressions ( [ eq : ma - c ] ) , ( [ eq : ma - d ] ) is used only on the level of coefficients of fractional parentage . in separate cases ,
e.g. when the two - particle operator acts upon one or two shells , it is possible to use expressions which exploit the racah algebra at a higher level , i.e. to take more advantage of the tensor algebra elements ( see judd  @xcite , jucys and savukynas  @xcite ) .
for example , let us investigate the case when a matrix element is calculated for bra and ket functions having one shell only . the tensorial forms ( [ eq : op - c ] ) and ( [ eq : op - e ] ) are of value here .
taking the second one of these , we have @xmath267{c } ( nl^{n}\alpha ls|\widehat{g}_{ii}|nl^{n}\alpha ^{\prime}l^{\prime}s^{\prime})= \\
= \frac 12 \left ( n\lambda n\lambda    n\lambda \right ) \left\ { \left [ \kappa _ 1,\kappa _ 2,\sigma_1,\sigma _ 2\right ] ^{-1/2 } \right .
\times \\ \times ( nl^{n}\alpha ls| \left [ \left [ a^{\left ( \lambda \right ) } \times \stackrel{\sim } { a}^{\left ( \lambda \right ) } \right ] ^{\left ( \kappa _ 1\sigma _ 1\right ) } \times \left [ a^{\left ( \lambda \right ) } \times \stackrel{\sim } { a}^{\left ( \lambda \right ) } \right ] ^{\left ( \kappa _ 2\sigma _ 2\right ) } \right ] ^{\left ( kk\right ) }    l^{\prime } s^{\prime } ) - \\ -\left\ { \begin{array}{ccc } \kappa _ 1 & \kappa _ 2 & k \\ l & l & l \end{array } \right\ } \left\ { \begin{array}{ccc } \sigma _ 1 & \sigma _ 2 & k \\ s & s & s \end{array } \right\ } \times \\ \left . \times ( nl^{n}\alpha ls| \left [ a^{\left ( \lambda \right ) } \times \stackrel{\sim } { a}^{\left ( \lambda \right ) } \right ] ^{\left ( kk\right ) }   |nl^{n}\alpha l^{\prime } s^{\prime } ) \right\}. \end{array}\ ] ]    using the relationships between the tensorial product of creation and annihilation operators and the tensorial quantities @xmath166 and @xmath167 ( see rudzikas and kaniauskas @xcite ) , the expression ( [ eq : ma - f ] ) for matrix elements can be writen down in terms of @xmath166 and @xmath167 . in comparing ( [ eq : ma - d ] ) to ( [ eq : ma - f ] ) we see that the summation over intermediate terms @xmath263 , @xmath264 , @xmath265 is already performed in expression ( [ eq : ma - f ] ) .
so , in this case the racah algebra is exploited at the level of standard quantities @xmath166 ir @xmath167 .
this simplifies calculations a lot :    * for zero matrix elements are easily tracked down from triangular conditions even before the actual calculation is performed . in case ( [ eq : ma - f ] ) only the triangular conditions @xmath268 and @xmath269 are present , but their number may be greater in other cases .
( in the above , the notation @xmath268 means the triangular condition @xmath270 . ) * the tables of standard quantities ( see nielson and koster @xcite , karazija _ et al _
@xcite , cowan  @xcite ) may be used .
* the recoupling matrix is simpler in this case , and it has an analytical expression .
so , the expressions exploiting the racah algebra at the level of @xmath166 and @xmath167 are much more advantageous than ( [ eq : ma - c ] ) .
such expressions are obtained for all physical operators .
for example , the expressions for spin - other - orbit operator are presented in papers horie @xcite , karazija _ et al _
@xcite and vizbarait .
e _ et al _
@xcite , the ones for spin - spin operator - in papers horie @xcite and karazija _ et al _  @xcite , and the ones for orbit - orbit operator in the monograph jucys and savukynas @xcite .
the shortcoming of the expressions of this type is that the racah algebra is exploited to its full extent in separate cases only .
this is discussed in detail in paper by gaigalas _
_  @xcite .
gaigalas _ et al _
@xcite have proposed a methodology which allows one to take all the advantages of the racah algebra in the most general case . according to the approach by gaigalas _
et al _ @xcite , a general expression of submatrix element for any two - particle operator between functions with @xmath22 open shells can be written down as follows :    @xmath271{c } ( \psi _ u^{bra}\left ( ls\right ) || \widehat{g}^{\left ( \kappa _ 1\kappa _ 2k,\sigma_1\sigma _ 2k\right ) }    } \right ) ) = \\
= \displaystyle { \sum_{n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } } } ( \psi _ u^{bra}\left ( ls\right ) ||\widehat{g } \left ( n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } \right )    } \right ) ) = \\ = \displaystyle { \sum_{n_il_i , n_jl_j , n_i^{\prime } l_i^{\prime } , n_j^{\prime } l_j^{\prime } } } \displaystyle { \sum_{\kappa _ { 12},\sigma _ { 12},\kappa _ { 12}^{\prime } , \sigma _ { 12}^{\prime } } } \sum \left ( -1\right ) ^\delta \theta ^{\prime } \left ( n_i\lambda _ i , n_j\lambda _ j , n_i^{\prime } \lambda _ i^{\prime } , n_j^{\prime } \lambda _ j^{\prime } , \xi \right ) \times \\ \times t\left ( n_i\lambda _ i , n_j\lambda _ j , n_i^{\prime } \lambda _ i^{\prime } , n_j^{\prime } \lambda _ j^{\prime } , \lambda ^{bra},\lambda ^{ket},\xi , \gamma \right ) \times \\ \times r\left ( \lambda _ i,\lambda _ j,\lambda _ i^{\prime } , \lambda _ j^{\prime } , \lambda ^{bra},\lambda ^{ket},\gamma \right ) . \end{array}\ ] ]    in calculating the spin - angular part of a submatrix element using ( [ eq : ma - g ] ) , one has to compute the following quantities ( for more detail see gaigalas  @xcite ) :    1 .   the recoupling matrix @xmath272 .
this recoupling matrix accounts for the change in going from matrix element + @xmath273 , which has @xmath22 open shells in the bra and ket functions , to the submatrix element + @xmath274 , which has only the shells being acted upon by the two - particle operator in its bra and ket functions . 2 .   the submatrix element @xmath275 , which denotes the submatrix elements of operators ( [ eq : op - j ] ) - ( [ eq : op - n ] ) . here
@xmath276 refers to the array of coupling parameters connecting the recoupling matrix @xmath272 to the submatrix element .
phase factor @xmath277 ( for more detail see gaigalas  @xcite ) .
4 .   @xmath278 , which is proportional to the radial part and corresponds to one of @xmath279, ... ,@xmath280 .
it consists of a submatrix element @xmath281 , and in some cases of simple factors and 3@xmath282-coefficients ( for more detail see gaigalas  @xcite ) .
in the next section we shall discuss finding of the recoupling matrix .
now we shall analyse the submatrix element @xmath275 . as the angular part of expression ( [ eq : ma - g ] ) contains the tensors ( [ eq : op - j ] ) - ( [ eq : op - n ] ) , so we will discuss the derivation of submatrix elements of these operators , and present the expressions for these quantities .
it is worth noting that these tensorial quantities all act upon the _ same _ shell .
so , all the advantages of tensor algebra and the quasispin formalism may be exploited efficiently .
we obtain the submatrix elements of operator ( [ eq : op - j ] ) by straightforwardly using the wigner - eckart theorem in quasispin space ( see rudzikas @xcite ) : @xmath283{c } \left ( l^n\;\alpha qls||a_{m_q}^{\left ( qls\right ) } ||l^{n^{\prime } } \;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } \right ) = -\left [ q\right ] ^{-1/2}\left [ \begin{array}{ccc } q^{\prime } & 1/2 & q \\
m_q^{\prime } & m_q & m_q \end{array } \right ] \times \\ \times \left ( l\;\alpha qls|||a^{\left ( qls\right ) } |||l\;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } \right ) , \end{array}\ ] ] where the last multiplier in ( [ eq : ma - h ] ) is the so - called completely reduced ( reduced in the quasispin , orbital and spin spaces ) matrix element .
the value of the submatrix element of operator ( [ eq : op - k ] ) is obtained by    @xmath284{c } \left ( nl^n\;\alpha qls||\left [ a_{m_{q1}}^{\left ( q\lambda \right ) } \times a_{m_{q2}}^{\left ( q\lambda \right ) } \right ] ^{\left ( k_1k_2\right ) } ||nl^{n^{\prime } } \;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } \right ) = \\
= \displaystyle { \sum_{\epsilon , m_\epsilon } } \left [ q\right ] ^{-1/2}\left [ \begin{array}{ccc } q & q & \epsilon \\ m_{q1 } & m_{q2 } & m_\epsilon \end{array } \right ] \left [ \begin{array}{ccc } q^{\prime } & \epsilon & q \\ m_q^{\prime } & m_\epsilon & m_q \end{array } \right ] \times \\ \times \left ( nl\;\alpha qls|||w^{\left ( \epsilon k_1k_2\right ) } |||nl\;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } \right ) .
\end{array}\ ] ]    on the right - hand side of equations ( [ eq : ma - h ] ) and ( [ eq : ma - i ] ) only the clebsch - gordan coefficient @xmath285 $ ] depends on the number @xmath286 of equivalent electrons .
@xmath287 denotes reduced in quasispin space submatrix element ( completely reduced matrix element ) of the triple tensor @xmath288 ^{\left ( \epsilon k_1k_2\right ) } $ ] .
it is related to the completely reduced coefficients ( subcoefficients ) of fractional parentage in a following way :    @xmath289{c } \left ( nl\;\alpha qls|||w^{\left ( \epsilon k_1k_2\right ) } |||nl\;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } \right ) = \\ = \left ( -1\right ) ^{q+l+s+q^{\prime } + l^{\prime } + s^{\prime } + \epsilon + k_1+k_2}\left [ \epsilon , k_1,k_2\right ] ^{1/2}\times \\ \times \displaystyle { \sum_{\alpha ^{\prime \prime } q^{\prime \prime } l^{\prime \prime } s^{\prime \prime } } } \left ( l\;\alpha qls|||a^{\left ( qls\right ) } |||l\;\alpha ^{\prime \prime } q^{\prime \prime } l^{\prime \prime } s^{\prime \prime } \right ) \times \\ \times \left ( l\;\alpha ^{\prime \prime } q^{\prime \prime } l^{\prime \prime } s^{\prime \prime } |||a^{\left ( qls\right ) } |||l\;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } \right ) \times \\ \times \left\ { \begin{array}{ccc } q & q & \epsilon \\ q^{\prime } & q & q^{\prime \prime } \end{array } \right\ }
\left\ { \begin{array}{ccc } l & l & k_1 \\ l^{\prime } & l & l^{\prime \prime } \end{array } \right\ } \left\ { \begin{array}{ccc } s & s & k_2 \\ s^{\prime } & s & s^{\prime \prime } \end{array } \right\ } .
\end{array}\ ] ]    in the other three cases ( [ eq : op - l ] ) , ( [ eq : op - m ] ) , ( [ eq : op - n ] ) we obtain the submatrix elements of these operators by using ( 2.28 ) of jucys and savukynas @xcite : @xmath290{c } ( nl^n\;\alpha qls||\left [ f^{\left ( \kappa _ 1\sigma _ 1\right ) } \left ( n\lambda \right ) \times g^{(\kappa _ 2\sigma _ 2)}\left ( n\lambda \right ) \right ] ^{\left ( kk\right ) } ||nl^{n^{\prime } } \;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } ) = \\ = \left ( -1\right ) ^{l+s+l^{\prime } + s^{\prime } + 2k}\left [ k\right ] \times \\ \times \displaystyle { \sum_{\alpha ^{\prime \prime } q^{\prime \prime } l^{\prime \prime } s^{\prime \prime } } } ( nl^n\;\alpha qls||f^{\left ( \kappa _ 1\sigma _ 1\right ) } \left ( n\lambda \right ) ||nl^{n^{\prime \prime } } \;\alpha ^{\prime \prime } q^{\prime \prime } l^{\prime \prime } s^{\prime \prime } ) \times \\ \times ( nl^{n^{\prime \prime } } \;\alpha ^{\prime \prime } q^{\prime \prime } l^{\prime \prime } s^{\prime \prime } ||g^{(\kappa _ 2\sigma _ 2)}\left ( n\lambda \right ) ||nl^{n^{\prime } } \;\alpha ^{\prime } q^{\prime } l^{\prime } s^{\prime } ) \times \\ \times \left\ { \begin{array}{ccc } \kappa _ 1 & \kappa _ 2 & k \\ l^{\prime } & l & l^{\prime \prime } \end{array } \right\ } \left\ { \begin{array}{ccc } \sigma _ 1 & \sigma _ 2 & k \\ s^{\prime } & s & s^{\prime \prime } \end{array } \right\ } , \end{array}\ ] ] where @xmath291 , @xmath292 are one of ( [ eq : op - j ] ) or ( [ eq : op - k ] ) and the submatrix elements correspondingly are defined by ( [ eq : ma - h ] ) , ( [ eq : ma - i ] ) and ( [ eq : ma - j ] ) .
@xmath293 is defined by the second quantization operators occurring in @xmath294 and @xmath295 .    as is seen , by using this approach gaigalas _
et al _ @xcite , the calculation of the angular parts of matrix elements between functions with @xmath22 open shells is reduced to requiring the reduced coefficients of fractional parentage or the tensors ( for example @xmath296 ) , which are independent of the occupation number of the shell and are acting on one shell of equivalent electrons .
the main advantage of this approach is that the standard data tables in such a case will be much smaller in comparison with tables of the usual coefficients @xmath297 @xmath298 ( see jucys and savukynas @xcite ) and , therefore , many summations will be less time - consuming . also one can see that in such an approach the submatrix elements of standard tensors and subcoefficients of fractional parentage actually can be treated in a uniform way as they all are the completely reduced matrix elements of the second quantization operators . hence , all methodology of calculation of matrix elements will be much more universal in comparison with the traditional one ( see cowan  @xcite , jucys and savukynas @xcite , wybourne @xcite ) .
while seeking the matrix elements of one- or two- particle operators , it is necessary to obtain the values of a recoupling matrix    @xmath272 , if we use the methodology by gaigalas _
et al _  @xcite ( see expresion ( [ eq : ma - g ] ) ) , or the recoupling matrices    @xmath34 and    @xmath35 , if we use the methodology by fano  @xcite ( see expresion ( [ eq : ma - c ] ) ) . in the case of several open shells
the expressions for matrix elements of every physical operator are published , where the recoupling matrices are in the form of simple factors .
usually all these are written in the coordinate representation .
they can be found in karazija _
@xcite , jucys and savukynas  @xcite , karazija  @xcite , rudzikas  @xcite .
meanwhile , for the more complex configurations , i.e. the ones having many shells , the recoupling matrices are much more complicated . beside that , the complexity of a two - particle operator adds to this .
when the tensorial structure of an operator is complex , the recoupling matrix is rather complex , too , e.g. the spin - other - orbit operator ( see gaigalas  @xcite ) . while attempting to calculate the angular part of matrix elements in all the mentioned cases , a general methodology for calculating the recoupling matrices
is necessary .
it has to be efficient , too , because the speed of calculation of angular parts of matrix elements depends on that .
the majority of methodologies to obtain angular parts are based on the fano  @xcite calculation scheme ( see for example glass  @xcite , glass and hibbert  @xcite , grant  @xcite ) . in finding the matrix elements using this , one of the tasks is to obtain the recoupling matrices for direct and exchange terms .
let us treat a matrix element of a two - particle operator acting upon four distinct shells    @xmath299    for such an operator , the recoupling matrix in @xmath300 - space , using graphical technique of jucys and bandzaitis  @xcite , is represented in figure  [ fig : a ] .
as we speak of calculations in @xmath0 - coupling , the analogous recoupling matrix in @xmath301 - space has to be calculated . in the figure  [ fig : a ] the tree containing * b * nodes represents the bra function , and the one with * k * nodes represents the ket function .
bra , as well as ket function , both contain @xmath22 open shells in this case .
now let us treat a case when three distinct shells are acted upon by the operator .
for example , the matrix element of operator @xmath302 is sought .
then in the recoupling matrix the node iii will be on @xmath303 line ( see figure  [ fig : a ] ) .
this means that there will be two nodes on the resulting orbital moment @xmath303 line of shell @xmath44 .
these will be nodes i and iii .
the resulting orbital moment @xmath304 of the shell @xmath44 of bra function is attributed to the line connecting nodes @xmath305 and i. the resulting orbital moment @xmath306 of the shell @xmath44 of ket function is attributed to the line between nodes iii and @xmath307 .
finally , characteristics @xmath308 of all possible intermediate momenta for that particular shell , which are summed over , are attributed to the line between nodes i and iii .
so , in treating the operator acting upon three distinct shells , an intermediate sum in the recoupling matrix appears .
similarly , two intermediate sums occur when two distinct shells are acted upon by the operator . and
when it acts upon one shell only , three intermediate sums are present .
the lines @xmath114 , @xmath309 and @xmath310 in the recoupling matrix represent the structure of two- particle operator in @xmath36 - space .
see tutlys  @xcite , grant  @xcite for details .
this is also valid for the exchange term , only the line @xmath311 must be connected to the node * d * , and line @xmath312 - to the node * c*.    the first program to calculate the recoupling matrices of this type , njsym , was written by burke  @xcite .
it performs the calculations in two stages : 1 ) the recoupling matrix is expressed as a sum of products of the @xmath313 - coefficients ; 2 ) this expression is used in calculation .
tutlys  @xcite wrote a program to calculate angular parts of matrix elements , angula , which expressed the recoupling matrix in terms of clebsch - gordan coefficients before the actual calculations .
while finding the recoupling matrix by the clebsch - gordan coefficient summation , this program eliminates trivial coefficients from the expression .
bar - shalom and klapisch  @xcite developed a new program njgraf .
this program calculates the recoupling matrix in several stages . on the basis of graphical methodology by yutsis , levinson and vanagas  @xcite , the recoupling matrix is analysed graphically and
an optimal expression is found .
afterwards , the value of recoupling matrix itself is calculated .
an analogous program recoup was written by lima  @xcite , and a program newgraph was written by fack _
et al _  @xcite .
all these ( njgraf , recoup and newgraph ) are based on the same principle .
an optimal analytical expression for the recoupling matrix is obtained by graphical method , and then the calculations are carried out according to it .
but the optimal expressions they find are different quite often , and are not really optimal .
as it was mentioned above , the methodology of angular calculation based on the fano calculation scheme has a shortcoming that the intermediate sums appear in complex recoupling matrices . due to these summations and the complexity of the recoupling matrix itself ,
the associated computer codes become rather time consuming . a solution to this problem
was found by burke _
they tabulated separate standard parts of recoupling matrices along with coefficients of fractional parentage at the beginning of a calculation and then used them later to calculate the coefficients needed .
computer codes by glass  @xcite , glass and hibbert  @xcite , burke _
@xcite , fischer  @xcite , fischer  @xcite and dyall _ et al _
@xcite utilize the program njsym ( burke  @xcite ) or njgraf ( bar - shalom and klapisch  @xcite ) for the calculation of recoupling matrices .
both are rather time consuming when calculating matrix elements of complex operators or electronic configurations when calculating matrix elements of complex operators or electronic configurations with many open shells . in order to simplify the calculations , cowan  @xcite suggested grouping matrix elements into classes ( see cowan  @xcite , figure 13 - 15 ) .
unfortunately this approach was not generalized to all two - electron operators .
perhaps this is the reason why cowan s approach is not widely used although the program itself , based on this approach is widely used .
gaigalas _ et al _
@xcite proposed a methodology where the analytical expressions for recoupling matrices are obtained for the most general case . in this methodology , analogically as in cowan  @xcite , the matrix elements are attributed to four different groups .
the operators acting upon only one shell belong to the first group ( distribution 1 from table 1 ) , the ones acting upon two - to the second ( distributions 2 - 10 from table 1 ) , upon three - to the third ( distributions 11 - 18 from table 1 ) , and upon four - to the fourth group ( distributions 19 - 42 from table 1 ) respectively .
each group has a different recoupling matrix , @xmath314 - @xmath315 .
they all are shown in figures  [ fig : b ] - [ fig : e ] .     representing the recoupling matrix @xmath272 when a two - particle operator acts upon a single shell.,height=510 ]     representing the recoupling matrix @xmath272 when a two - particle operator acts upon a two shell.,height=510 ]     representing the recoupling matrix @xmath272 when a two - particle operator acts upon a three shell.,height=510 ]     representing the recoupling matrix @xmath272 when a two - particle operator acts upon a four shells.,height=510 ]      let us treat the recoupling matrix @xmath316 first ( see figure  [ fig : b ] ) .
it represents the case when a physical two - particle operator acts upon a single shell . here
, similarly as in figure  [ fig : a ] , the graphical technique of jucys and bandzaitis  @xcite is used .
the bra and ket functions are represented as in figure  [ fig : a ] . as the two - particle operator acts upon only one shell in this case
, this operator may be treated as a single - particle one from the recoupling matrix point of view . using a graphical rule allowing one to cut two lines and connect the loose ends
, we disconnect the nodes @xmath317 and @xmath318 from the general recoupling matrix . the part separated from recoupling
corresponds to delta - functions only , @xmath319 , @xmath320 .
in addition , using graphical technique it is possible to cut all the nodes until @xmath305 and @xmath321 out of the general recoupling matrix .
all the cut nodes contribute the same delta functions which may be written as @xmath322 , @xmath323 .    in the next stage , it remains to treat the part of recoupling matrix from the nodes @xmath305 and @xmath324 up to @xmath325 and @xmath326 . beside that , the lines @xmath327 and @xmath328 are connected in this remaining diagram . before obtaining the analytical expression for the remaining recoupling matrix we introduce a notation @xmath329 .
it will help us to describe the procedure of obtaining analytical expression for the remaining recoupling matrix in a simpler way .
so , the diagram in figure  [ fig : f ] is denoted as @xmath329 .
it equals to    .,height=529 ]    @xmath330 ^{1/2}\left\ { \begin{array}{ccc } j_{1 } & j_{2 } & j_{12 } \\ j_{3 } & j      & j_{32 } \end{array } \right\ } .\ ] ]    now , using the graphical rule of cutting three lines , we cut out the nodes @xmath331 , @xmath332 and @xmath333 from the recoupling matrix .
thus we obtain a diagram @xmath334 .
this coefficient in the paper by gaigalas _
@xcite is denoted as @xmath335 , and equals to @xmath336 ^{1/2}\left\ { \begin{array}{ccc } k & l_a^{\prime } & l_a \\ j & t & t^{\prime } \end{array } \right\ } .\ ] ]    the diagram is cut from the nodes @xmath337 and @xmath338 up to @xmath339 and @xmath340 in the same way .
it is expressed in terms of diagrams    @xmath341{c } c_2\left ( k , k_{\min } , k_{\max } \right ) = \displaystyle { \stackrel{k_{\max } } { \prod_{i = k_{\min } } } } a_{1}(+ - + - ; k l'_{12 ... i } l_{12 ... i } ; l_{i } l_{12 ... i-1 } l'_{12 ... i-1 } ) .
\end{array}\ ] ]    it is easy to notice that @xmath342 corresponds to the coefficient @xmath342 defined in gaigalas _
@xcite .      before going into investigation of the recoupling matrix pictured in figure  [ fig : c ]
, we denote the diagram @xmath343 ( see figure  [ fig : g ] ) as    @xmath344 .
its analytical expression is    .,height=453 ]    @xmath345 ^{1/2}\left\ { \begin{array}{ccc } j_{1 }   & j_{2 }   & j_{12 } \\ j_{3 }   & j_{4 }   & j_{34 } \\
j_{13 } & j_{24 } & j \end{array } \right\ } .\ ] ]    in the @xmath346 diagram , first the signs of nodes 1 - 6 , and then the moments are presented .
while investigating the recoupling matrix ( see figure  [ fig : c ] ) it is easy to notice that there is an additional graphic element . after cutting this out , we get a diagram the further investigation of which is very similar to the investigation of @xmath314 . that additional matrix element is    @xmath347 .
it is obtained by cutting the nodes @xmath337 , @xmath348 , @xmath349 , @xmath338 out of the diagram @xmath350 . after doing this ,
the diagram under investigation is split into two diagrams closely resembling the recoupling matrix @xmath316 .
as they are cut in the very same way as @xmath316 is , we will not stop at recoupling matrix @xmath350 for details .      in the methodology presented in gaigalas _
@xcite , the recoupling matrix is represented by a diagram @xmath351 ( see figure  [ fig : d ] ) when the two - particle operator is acting upon three shells .
similarly as for the recoupling matrix @xmath350 , the cutting out of nodes @xmath337 , @xmath348 , @xmath349 and @xmath352 leads to a diagram @xmath353 , and the cutting out of nodes @xmath354 , @xmath355 , @xmath356 ir @xmath357 leads to a diagram    @xmath358 . in that case
the diagram splits up into three simple diagrams , the further investigation of which is analogous to that of @xmath316 recoupling matrix .
the most complex recoupling matrix @xmath359 ( see figure  [ fig : e ] ) comes into view when all the second quantization operators act upon different shells .
the part of this recoupling matrix from nodes @xmath317 and @xmath318 to nodes @xmath360 and @xmath361 , and also from nodes @xmath362 and @xmath363 to nodes @xmath325 and @xmath326 is very similar to the recoupling matrices described above , therefore their investigation is the same as that of the diagrams described earlier
. we will investigate in more detail the part of diagram from nodes @xmath360 and @xmath361 to nodes @xmath362 and @xmath363 .
that is pictured in figure  [ fig : e - b ] . after cutting the diagram across four lines @xmath364 , @xmath365 , @xmath366 and @xmath367 and connecting these according the graphical technique of jucys and bandzaitis  @xcite with a generalized clebsch - gordan coefficient , then connecting the lines @xmath368 , @xmath366 and @xmath369 with an ordinary clebsch - gordan coefficient
, we obtain a diagram which is expressed via two @xmath313-coefficients . the exact analytical expression of this diagram is presented in the paper by gaigalas _
@xcite ( see expresion ( 29 ) ) .
.,height=510 ]    after cutting the nodes @xmath362 , @xmath356 , @xmath370 , @xmath371 and @xmath363 out of the remaining part and investigating the obtained diagram according to the jucys and bandzaitis  @xcite technique , we obtain a diagram which corresponds to a @xmath372- coefficient .
the analytical expression of that is presented in paper by gaigalas _
@xcite , expr .
( 31 ) .
if a shell @xmath373 is beside shell @xmath374 , i.e. there are no additional nodes between nodes @xmath360 , @xmath361 and @xmath362 , @xmath363 , then the additional graphical elements vanish .
so , we obtained the full analytical expression of recoupling matrix @xmath375 .
if between the abovementioned shells a single additional shell is present , then we obtain an additional diagram which is expressed as a product of two @xmath313-coefficients ( see ( 30 ) from gaigalas _ et al _
if the number of shells is greater , then the obatined diagram in the most general case is written via ( 32 ) from gaigalas _ et al _  @xcite .
the methodology that is based on the second quantization in coupled tensorial form , on the angular momentum theory in three spaces ( orbital , spin and quasispin ) , on the wick s theorem and on the generalized graphical technique of angular momentum , give the possibility to efficiently calculate the matrix elements of energy operators in general case ( gaigalas and rudzikas  @xcite , gaigalas _
@xcite , this paper ) .
the main ideas of this methodology are :    1 .
a number of theoretical methods is known in atomic physics that facilitate a lot the treatment of angular part of matrix elements .
these are the theory of angular momentum , its graphical representation , the quasispin , the second quantization and its coupled tensorial form . but while treating the matrix elements of physical operators in general , the methods mentiond above are applied only some part , or very inefficiently .
+ an idea is presented and carried out in the work gaigalas and rudzikas  @xcite , gaigalas _ et al _  @xcite and this paper , of unifying all these methods in order to optimize the way to treat matrix elements of physical operators , to investigate efficiently even the most complex cases of atoms and ions .
the wick s theorem is widely known in the theory of atom .
up to now it was applied to the most general products of the operators of second quantization , i.e. when the particular arrays of quantum numbers for each operator of second quantization were not yet defined . in the methodology
described there is proposed to apply the wick s theorem for the products of second quantization operators where they have the values of quantum numbers already defined .
this allows one , using the methodology of second quantization in coupled tensorial form , to abtain immediately the optimal tensorial expressions for any operator .
then , in treting the matrix elements of physical operators , the advantages of a new modification of the racah algebra are exploited to their full extent .
3 .   while analysing the physical operators presented as products of tensors @xmath376 , @xmath377 , @xmath378^{(k_q k_l k_s)}$ ] , @xmath379^{(k_q k_l k_s)}$ ] , + @xmath380^{(k_q k_l k_s)}$ ] , it is possible to obtain convenient analytical expressions for recoupling matrices that must be taken into account in finding the matrix elements ( non - diagonal with respect to configuration included ) of any physical operator between complex configurations ( with any number of open shells ) .
an idea is proposed and carried out , concerning the most efficient way to apply the tensor algebra and the quasispin formalism in the most general case , for the diagonal and the non - diagonal matrix elements as well , when the bra- and ket- functions have any number of open shells .
the combination of all these improvements allows to efficiently account for correlation efects practically for any atom or ion of periodical table .
the author is grateful to professor z. rudzikas for encouraging and valuable remarks .                                  ,
`` on the method of calculation of matrix elements of the operators of atomic quantities in the case of complex configurations '' , proceedings of the academy of sciences of lithuanian ssr , b series v. 4 ( 27 ) , p. 45 - 57
( 1961 ) ( in russian ) .                , `` irreducible tensorial form of the effective hamiltonian of an atom and the diagrammatic representation in the first two orders of the stationary perturbation theory '' , liet .
fiz . rink .
( english translation sov .
phys . coll . ) , v. 25 , no .
6 , p. 14 - 31
( 1985 ) .    , `` energies of @xmath286 equivalent electrons expressed in terms of two - electron energies and independent three - electron parameters : a new complete set of orthogonal operators : i. theory '' , j. phys .
b : at . mol .
phys . , v. 17 , p. 2375 - 2392
( 1984 ) .
, `` reduced in quasi - spin space coefficients of fractional parentage and matrix elements of tensorial operators '' , liet .
fiz . rink . , v. 18 , no .
293 ( 1978 ) .
( english translation sov .
phys . coll . * 18*(3 ) , 1 ( 1978 ) from allerton press , new york ) .
, `` expressions for reduced matrix elements of non - relativistic effective hamiltonian of atom in the first orders of stationary perturbation theory for configurations with two unifilled shells @xmath383 '' , in : _ spectroscopy of autoionized states of atoms and ions _
( moscow , scientific council of spectroscopy , 1985 ) , p. 20 - 42 .                                  , `` diagramatic approach to irreducible tensorial operators and its application to perturbation theory '' , in : _ spectroscopy of autoionized states of atoms and ions _
( moscow , scientific council of spectroscopy , 1985 ) , p. 43 - 61 .    ,
`` application of the graphical methods of angular momentum theory to the investigation of the expansion of stationary perturbation theory '' , izv .
uchebn . zaved .
50 , no . 7 , p. 1403 - 1410
( 1986 ) .                          , _
tables for the calculation of matrix elements of atomic quantities _ , ( moscow , 1967 ) ; english translation by e.k .
wilip , anl - trans-563 ( national technical information service , springfield , va . , 1968 ) .
, `` a program to calculate the matrix elements in multi - configuration approximation '' , in : _ sbornik programm po matematicheskomy obespecheniju atomnyh raschotov_. vipusk 4 .
( vilnius , 1980 ) ( in russian ) .