Next: The subroutine WAP1LS
Up: SAI_SQLS1
Previous: The subroutine W1
The routine determines the value of the matrix elements:
![\begin{displaymath}
( l^{N} QLS \vert\vert [ a^{(qls)}_{m_{q1}} \times [ a^...
...(k_{1}k_{2})} ]^{(k_{l}k_{s})} \vert\vert l^{N'} Q'L'S' ).
\end{displaymath}](img389.png) |
(37) |
While calculating cases where the orbital number l=0, 1, 2, 3 and the
shell's occupation number
, the program relies on the expression (31)
from the paper P1 [3]. In that case,
the subroutine IZAS1 checks that the subshell has a state
with the specified characteristics. The subroutine ITLS2 finds the first and
the last numbers of the state from the running intermediate sum in array
MT. RUMT finds the shell's total angular momentum
LS and quasispin Q for
each intermediate state. The routine IXJTIK checks all triads.
The subroutine C0T5S finds the Clebsch-Gordan coefficient giving the
dependence on the shell's occupation number and SLS finds the reduced matrix
element of
tensor operator
(see in Section 3.2). The second part of
the expression is calculated by the routine W1
(see in Section 3.2.2). The routine SIXJ finds 6j- symbol.
In other cases, the program calculates according to the
expression (40) from paper P2 [4].
The subroutine has the arguments:
- IK is the array I for the bra function.
- BK is the array B for the bra function.
- ID is the array I for the ket function.
- BD is the array B for the ket function.
- K1 is the rank
.
- K2 is the rank
.
- K3 is the rank
.
- BK4 is the rank
.
- QM1, QM2 and QM3 are the quasispin projections in (37).
- AW is the value of the reduced matrix element (37) which is
returned by the subroutine.
Next: The subroutine WAP1LS
Up: SAI_SQLS1
Previous: The subroutine W1
2001-12-07