Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aabi_v2_0.gz(35 Kbytes)
Manuscript Title: A Fortran 77 version of 'a function subprogram in order to calculate the matrix elements of rotation operators'.
Authors: F. Brut
Program title: ROTATION MATRIX ELEMENTS DD
Catalogue identifier: AABI_v2_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 39(1986)297
Programming language: Fortran.
Computer: IBM 3081/K.
Operating system: VM/CMS.
RAM: 10K words
Word size: 32
Keywords: General purpose, Algebras, Quantum mechanics, Nuclear physics, Particle physics, Atomic physics, Molecular physics, Angular momentum projection techniques.
Classification: 4.2.

Nature of problem:
The reduced matrix elements of rotation operators are calculated. The sum is done for all the values for which the argument of the factorials are greater than or equal to zero. Each matrix element is written explicitly as a literal polynomial on the variables cos (beta/2) and sin (beta/2).

Solution method:
The reduced matrix elements djmm'(Beta) are homogeneous polynomials of degree 2j on the variables cos(Beta/2) and sin(Beta/2). With the phase convention of Wigner, namely:
djmm'(Beta) = (-1)m-m' djm'm(Beta) = dj-m',-m(Beta)
only matrix elements which satisfy simultaneously the following conditions
                 -j <= m <= 0                                            
                 |m'| <= |m|  
are considered. The others can be deduced from the above relations. Each single term of the polynomial includes three factors: a signed numerical constant, the appropriate power of cos(Beta/2) and the appropriate power of sin(Beta/2). Each of these factors is stored, as data for the first of them, in a variable name and the polynomial is then written explicitly.

Unusual features:
In the present version, the program assumes that the quantum number j can take integer and half-integer values from 0 up to and including 15/2. The range of the values of the spin j can be changed easily. In fact, this function subprogram is the output resulting from a FORTRAN program. Thus, it is somewhat easy to make reasonable changes on the boundaries of the spin j, if required.

Running time:
For instance, 15000 calls, for various values of Euler angle Beta and for all the values of j between 0 and 15/2, with 0 <= |m| or |m'| <= j, took less than .8 s CPU time for execution on a IBM 3081/K.