next up previous
Next: The subroutine WAP1LS Up: SAI_SQLS1 Previous: The subroutine W1

The subroutine AWP1LS

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} (37)

While calculating cases where the orbital number l=0, 1, 2, 3 and the shell's occupation number $N > 2$, 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 $a^{(qls)}$ 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:
  1. IK is the array I for the bra function.
  2. BK is the array B for the bra function.
  3. ID is the array I for the ket function.
  4. BD is the array B for the ket function.
  5. K1 is the rank $k_{1}$.
  6. K2 is the rank $k_{2}$.
  7. K3 is the rank $k_{l}$.
  8. BK4 is the rank $k_{s}$.
  9. QM1, QM2 and QM3 are the quasispin projections in (37).
  10. AW is the value of the reduced matrix element (37) which is returned by the subroutine.


next up previous
Next: The subroutine WAP1LS Up: SAI_SQLS1 Previous: The subroutine W1
2001-12-07