U111: Wigner 3-j, 6-j, 9-j Symbols; Clebsch-Gordan, Racah W-, Jahn U-Coefficients

Author(s): K.S. Kölbig Library: MATHLIB
Submitter: Submitted: 15.10.1994
Language: Fortran Revised:

Function subprograms RWIG3J, DWIG3J; RWIG6J, DWIG6J; RWIG9J, DWIG9J; RCLEBG, DCLEBG; RRACAW, DRACAW and RJAHNU, DJAHNU calculate the Wigner 3-j, 6-j and 9-j symbols, the Clebsch-Gordan coefficients, the Racah W-coefficients and the Jahn U-coefficients, respectively.

On CDC and Cray computers, the double-precision versions DWIG3J etc. are not available.

Structure:

FUNCTION subprograms
User Entry Names:
RWIG3J, RWIG6J, RWIG9J, RCLEBG, RRACAW, RJAHNU
DWIG3J, DWIG6J, DWIG9J, DCLEBG, DRACAW, DJAHNU

Usage:

In any arithmetic expression, for tex2html_wrap_inline164 (type REAL), or tex2html_wrap_inline166 (type DOUBLE PRECISION),
tWIG3J(A,B,C,X,Y,Z) has the value of tex2html_wrap_inline168 ;
tWIG6J(A,B,C,X,Y,Z) has the value of tex2html_wrap_inline170 ;
tWIG9J(A,B,C,P,Q,R,X,Y,Z) has the value of tex2html_wrap_inline172 ;
tCLEBG(A,B,C,X,Y,Z) has the value of tex2html_wrap_inline174 ;
tRACAW(A,B,C,D,E,F) has the value of tex2html_wrap_inline176 ;
tJAHNU(A,B,C,D,E,F) has the value of tex2html_wrap_inline178 .

All the arguments must have integral or half-integral values (see Notes). They have the same type as the function name. For definitions and notations see References.
The following relations hold (see Refs. 1 and 3):
Clebsch-Gordan coefficient (in terms of the Wigner 3-j symbol):

eqnarray93

Racah W-coefficient (in terms of the Wigner 6-j symbol):

eqnarray100

Jahn U-coefficient (in terms of the Wigner 6-j symbol and the Racah W-coefficient):

eqnarray106

Method:

The Wigner 3-j symbol and the Clebsch-Gordan coefficient are calculated from formulas (5.1) and (5.10) of Ref. 1, respectively. The Wigner 6-j symbol, the Racah W- and the Jahn U-coefficient are calculated from formulas (5.23) and (5.24) of Ref. 1. In both cases, the factorials are replaced by their logarithms during the calculation. The Wigner 9-j symbol is calculated from formula (5.37) of Ref. 1 in terms of Wigner 6-j symbols.

Notes:

A Wigner-3j symbol tex2html_wrap_inline206 is considered to be zero unless simultaneously
(i) tex2html_wrap_inline208 and tex2html_wrap_inline210 have both either integral or half-integral values (each i),
(ii) tex2html_wrap_inline214 (each i),
(iii) tex2html_wrap_inline218 ,
(iv) tex2html_wrap_inline220 is an integer,
(v) tex2html_wrap_inline222 is an integer and \ tex2html_wrap_inline224 .

The conditions (v) are often denoted by tex2html_wrap_inline226 and are called the triangle relations.

For a Clebsch-Gordan coefficient tex2html_wrap_inline228 , condition (iii) reads tex2html_wrap_inline230 and condition (iv) disappears.

A Wigner-6j symbol tex2html_wrap_inline234 is considered to be zero unless simultaneously
(i) all tex2html_wrap_inline236 and tex2html_wrap_inline238 have non-negative integral or half-integral values,
(ii) the four triangle relations tex2html_wrap_inline240 hold.

A Wigner-9j symbol tex2html_wrap_inline244 is considered to be zero unless simultaneously
(i) all tex2html_wrap_inline246 have non-negative integral or half-integral values,
(ii) the arguments in each row and in each column satisfy the triangle relations.

Restrictions:

The sum of arguments in any triangle relation must not exceed 100. No test is made.

References:

  1. R.D. Cowan, The theory of atomic structure and spectra, (Univ. of California Press, Berkeley CA 1981).
  2. A.F. Nikiforov, V.B. Uvarov and Yu.L. Levitan, Tables of Racah coefficients (Pergamon Press, Oxford 1965).
  3. M. Rotenberg, R. Bivins, N. Metropolis and J.K. Wooten, Jr., The 3-j and 6-j symbols (Crosby Lockwood, London 1959).
  4. D.A. Varshalovich, A.N. Moskalev and V.K. Khersonskii, Quantum theory of angular momentum (World Scientific, Singapore 1988).
tex2html_wrap_inline252

Michel Goossens Wed Jun 5 08:08:43 METDST 1996