We strongly condemn Russia's aggressive attack on Ukraine for which there is no rational justification. We stand firmly with the Ukrainian people.
Burkhard Schmidt and Petra Zdánská
The method of symmetry adapted wavepackets (SAWP) to solve the time-dependent Schrödinger equation for a highly symmetric potential energy surface is introduced. The angular dependence of a quantum-mechanical wavepacket is expanded in spherical harmonics where the number of close-coupled equations for the corresponding radial functions can be efficiently reduced by symmetry adaption of the rotational basis using the SAWP approach. Various techniques to generate symmetry adapted spherical harmonics (SASHs) for the point groups of highest symmetry (octahedral, icosahedral) are discussed. The standard projection operator technique involves the use of Wigner rotation matrices. Two methods to circumvent numerical instabilities occuring for large azimuthal quantum numbers are suggested. The first is based on a numerical scheme which employs Gaussian integrations yielding exact and stable results. The second is a recursive algorithm to generate higher order SASHs accurately and efficiently from lower order ones. The paper gives a complete set of "seed functions" generated by projection techniques which can be used to obtain SASHs for all irreducible representations of the octahedral and icosahedral point groups recursively.
Comp. Phys. Comm. 127 (2-3), 290-308 (2000)
DOI:10.1016/S0010-4655(99)00524-X
The (time-dependent or time-independent) Schrödinger equation for a highly symmetric potential energy surface can be solved by expanding the angular dependence of a quantum-mechanical wavefunction in symmetry adapted spherical harmonics (SASHs). While the generation of the required SASHs is trivial for simple symmetries, special techniques have been developed for the point groups of highest symmetry. Our work gives complete sets of “seed functions” to obtain SASHs for all irreducible representations of the octahedral and icosahedral point groups recursively. Basically, the technique consists of two steps:
A wavepacket in three-dimensional space is described as a sum of products of time-dependent radial wavefunctions and spherical harmonic functions. This approach yields a set of coupled equations for the radial functions where both the numerical effort and the storage requirement depend quadratically on the number of spherical harmonics used in the expansion scheme. This approach is known in the literature as close coupled wavepackets (CCWP).
In order to reduce this number, symmetry adapted spherical harmonics (SASHs) are used which transform according to an irreducible representation of a point group. Note that the usual technique employing projection operators (based on the great orthogonality theorem of group theory) causes difficulties when applied to cubic point groups and to spherical harmonics for high azimuthal quantum numbers. Instabilities of the conventional scheme based on Wigner rotation matrices [1] encountered for high angular momentum states can be overcome by a numerical scheme based on Gaussian quadratures [2]. An elegant alternative can be found in an iterative scheme [3] where low order SASHs ("seed functions") can be used to recursively generate all higher order SASHs with very little computational effort.
Tables of "seed functions" are available for the octahedral and icosahedral point groups. The tables given in the following are complete in the sense that all SASHs of all irreducible representations of the two point groups under consideration can be recursively generated from the seed functions.
It is a great pleasure to thank F. Neugebaur (Humboldt University, Berlin) and P. Zdanska (Heyrovsky Institute, Prague) for their invaluable help in writing computer codes and generating the above tables of SASHs. Furthermore, Ch. Salzmann's help in preparing the figures is highly appreciated.
C. J. Bradley and A. P. Cracknell
The Mathematical Theory of Symmetry in Solids
Clarendon, Oxford (1972)
G. W. Fernando, M. Weinert, R. E. Watson, and J. W. Davenport
J. Comp. Phys. 112, 282 (1994)
DOI:10.1006/jcph.1994.1100
W. Prandl, P. Schiebel, and K. Wulf
Acta Cryst. A52, 171 (1996)
DOI:10.1107/S0108767395012578