Programs in Physics & Physical Chemistry
|[Licence| Download | New Version Template] adfv_v10_0.tar.gz(532 Kbytes)|
|Manuscript Title: Maple procedures for the coupling of angular momenta. An up-date of the RACAH module.|
|Authors: S. Fritzsche|
|Program title: RACAH|
|Catalogue identifier: ADFV_v10_0|
Distribution format: tar.gz
|Journal reference: Comput. Phys. Commun. 180(2009)2021|
|Programming language: Maple 11 and 12.|
|Computer: All computers with a license for the computer algebra package Maple .|
|Operating system: Suse Linux 10.2+ and ubuntu 8.10.|
|Keywords: Angular momentum theory, bipolar harmonic, Clebsch-Gordan coefficient and expansion, Gaunt coefficient, Racah algebra techniques, recoupling coefficient, Sharp's symbol, spherical harmonic, spinor, vector and tensor spherical harmonic, sum rule evaluation, Wigner n-j symbol, D-function and rotation matrix, Yutsis graph.|
|PACS: 31.10.+z, 32.10-f.|
|Classification: 4.1, 5.|
Does the new version supersede the previous version?: Yes
Nature of problem:
The theories of angular momentum and spherical tensor operators, sometimes known also as Racah's algebra, provide a powerful calculus for studying spin networks and (quantum) many-particle systems. For an efficient use of these theories, however, one requires not only a reliable handling of a large number of algebraic transformations and rules but, more often than not, also a fast access to their standard quantities, such as the Wigner n-j symbols, Clebsch-Gordan coefficients, spherical harmonics of various kinds, the rotation matrices, and many others.
A set of Maple procedures has been developed and maintained during the last decade which supports both, algebraic manipulations as well as fast computations of the standard expressions and symbols from the theory of angular momentum [2, 3]. These procedures are based on a sizeable set of group-theoretical (and often rather sophisticated) relations which has been discussed and proven in the literature; see the monograph by Varshalovich et al.  for a comprehensive compilation. In particular the algebraic manipulation of complex (Racah) expressions may result in considerable simplifications, thus reducing the 'numerical costs', and often help obtain further insight into the behaviour of physical systems.
Reasons for new version:
A revision of the RACAH module became necessary for mainly three reasons: (i) Since the last extension of the RACAH procedures , which was developed within the framework of Maple 8, several updates of Maple were distributed by the vendors (currently Maple 13) and required a number of adaptations to the source code; (ii) the increasing size and program structure of the RACAH module made it advisible to separate the (procedures for the treatment of the) atomic shell model from the manipulation and computation of RACAH expressions. Therefore, the computation of angular coefficients for different coupling schemes, (grand) coefficients of fractional parentage as well as the matrix elements (of various irreducible tensors from the shell model) is to be maintained from now on independently within the JUCYS module; (iii) a number of bugs and inconsistencies have been reported to us and corrected in the present version.
Summary of revisions:
In more detail, the following changes have been made:
The (RACAH) program is based on the concept of Racah expressions [cf. Fig. 1 in Ref. ] which, in principle, may contain any number of Wigner n-j symbols (n ≤ 9), Clebsch Gordan coefficients, spherical harmonics and/or rotation matrices. In practise, of course, the required time and the success of an evaluation procedure depends on the complexity of the expressions and on the storage available, sometimes also on MAPLE's internal garbage treatment. In some cases, it is advisable to attempt first a simplification of the magnetic quantum numbers for a given expression before the summation over further 6-j and 9-j symbols should be taken into account. For all other quantities (that are compiled in Ref. , Tables 1 and 2 and explained in more detail in the Short Reference Manual, Racah-commands-2009.pdf), we currently just facilitate fast numerical computations by exploiting, as far as possible, MAPLE's hardware floating-point model.
The program also supports simplifications on the Wigner rotation matrices. In integrals over the rotation matrices, products of up to three Wigner D-functions or reduced matrices (with the same angular arguments) are recognized; for the integration over a solid angle, however, the domain of integration must be specified explicitly for the Euler angles α and γ in order to force MAPLE to generate a constant of integration.
In the course of the evaluation of Racah expressions, it is, in practice, often difficult to check internally whether all substructures of an expression are defined properly. Therefore, the user must ensure that all angular momenta (if given explicitly) must finally evaluate to integer and half-integer values and that they satisfy proper coupling conditions.
The RACAH program is designed for interactive use and for providing a quick and algebraic evaluation of (complex) expressions from Racah's algebra. In the evaluation, it exploits a large set of sum rules which are known from Racah's algebra and which may include (multiple) summations over dummy indices; see Varshalovich et al.  for a more detailed account of the theory. One strength of the program is that it recognizes automatically the symmetries of the symbols and functions, and that it applies also (some of) the graphical rules due to Yutsis and coworkers . As before, the result of the evaluation process will be provided as Racah expressions, if a further simplification could be achieved, and may hence be used for further derivations and calculations within the given framework.
In dealing with recoupling coefficients, these coefficients can be entered simply as a string of angular momenta (variables), separated by commas, and very similar to how they appear in mathematical texts. This is a crucial advantage of the program, compared with previous developments, for which the angular momenta and coupling coefficients had often to be given in a very detailed format.
A Short Reference Manual to all procedures of the RACAH program is provided by this distribution; it also contains the worksheet Racah-tests-2009-maple12.mw that contains the examples from all previous versions and may help test the installation. This worksheet can serve as a first tutorial to the RACAH procedures.
In the past, the RACAH program has been utilized extensively in a number of applications including angular and polarization studies of heavy ions , angular distributions and correlation functions following photon-induced excitation processes , entanglement studies , in application of point-group symmetries and several others.
The worksheet supplied with the distribution takes about 1 minute to run.
|||Maple is a registered trademark of Waterloo Maple Inc.|
|||S. Fritzsche, Comp. Phys. Commun. 103, 51 (1997).|
|||S. Fritzsche, S. Varga, D. Geschke, and B. Fricke, Comp. Phys. Commun. 111, 167 (1998). T. Ingho , S. Fritzsche, and B. Fricke, Comp. Phys. Commun. 139, 297 (2001). S. Fritzsche, T. Ingho , T. Bastug and M. Tomaselli, Comp. Phys. Commun. 139, 314 (2001).|
|||J. Pagaran, S. Fritzsche and G. Gaigalas; Comput. Phys. Commun. 174, 616 (2006).|
|||D. A. Varshalovich, A. N. Moskalev, V. K. Khersonskii, Quantum Theory of Angular Momentum (World Scientific, Singapore a.o., 1988).|
|||K. Schulten and R. G. Gordon, Comp. Phys. Commun. 11, 269 (1976).|
|||G. Gaigalas, S. Fritzsche, B. Fricke, Comp. Phys. Commun. 135, 219 (2001); G. Gaigalas and S. Fritzsche, Comp. Phys. Commun. 149, 39 (2002); G. Gaigalas, O. Scharf and S. Fritzsche; Comput. Phys. Commun. 166, 141 (2005).|
|||S. Fritzsche, T. Ingho and M. Tomaselli, Comp. Phys. Commun. 153, 424 (2003).|
|||A. P. Yutsis, I. B. Levinson and V. V. Vanagas, The Theory of Angular Momentum (Israel Program for Scientific Translation, Jerusalem, 1962).|
|||S. Fritzsche, P. Indelicato and T. Stöhlker, J. Phys. B38, S707 (2005).|
|||M. Kitajima, M. Okamoto, M. Hoshino et al. J. Phys. B35, 3327 (2002); N. M. Kabachnik, S. Fritzsche, A. N. Grum-Grzhimailo et al. Phys. Reports 451, 155 (2007); S. Fritzsche, A. N. Grum-Grzhimailo, E. V. Gryzlova and N. M. Kabachnik, J. Phys. B41, 165601 (2008).|
|||T. Radtke and S. Fritzsche, Comput. Phys. Commun. 175, 145 (2006).|
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|