Evaluation of the rotation matrices in the basis of real spherical harmonics

Abstract Rotation matrices (or Wigner D functions) are the matrix representations of the rotation operators in the basis of spherical harmonics. They are the key entities in the generation of symmetry-adapted functions by means of projection operators. Although their expression in terms of ordinary (complex) spherical harmonics and Euler rotation angles is well known, an alternative representation using real spherical harmonics is desirable. The aim of this contribution is to obtain a general algorithm to compute the representation matrix of any point-group symmetry operation in the basis of the real spherical harmonics, paying attention to the use of recurrence relationships that allow the treatment of functions with high angular momenta.