Recursive computation of spherical harmonic rotation coefficients of large degree

DOI10.1007/978-3-319-13230-3_5MaRDI QIDQ2950248zbMATH OpenOpenAlexFDO

Authors R. Duraiswami, N. A. Gumerov

Publication date 8 October 2015

Published in Excursions in Harmonic Analysis, Volume 3 (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1403.7698

zbMATH Keywords

\(\mathrm{SO}(3)\)algorithmCourant-Friedrichs-Lewy conditionfast Fourier transformfast multipole methodsfinite difference schemesHelmholtz equationLaplace equationnumerical examplesrecursionsrotationspherical harmonicsstabilityvon Neumann methodWigner \(d\)-matrix

Mathematics Subject Classification ID

Spherical harmonics (33C55) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Finite difference methods for boundary value problems involving PDEs (65N06) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Boundary element methods for boundary value problems involving PDEs (65N38) Numerical methods for discrete and fast Fourier transforms (65T50) Complexity and performance of numerical algorithms (65Y20)

Abstract: Computation of the spherical harmonic rotation coefficients or elements of Wigner's d-matrix is important in a number of quantum mechanics and mathematical physics applications. Particularly, this is important for the Fast Multipole Methods in three dimensions for the Helmholtz, Laplace and related equations, if rotation-based decomposition of translation operators are used. In these and related problems related to representation of functions on a sphere via spherical harmonic expansions computation of the rotation coefficients of large degree $n$ (of the order of thousands and more) may be necessary. Existing algorithms for their computation, based on recursions, are usually unstable, and do not extend to $n$. We develop a new recursion and study its behavior for large degrees, via computational and asymptotic analyses. Stability of this recursion was studied based on a novel application of the Courant-Friedrichs-Lewy condition and the von Neumann method for stability of finite-difference schemes for solution of PDEs. A recursive algorithm of minimal complexity $Oleft(n^{2}ight)$ for degree $n$ and FFT-based algorithms of complexity $Oleft(n^{2}lognight)$ suitable for computation of rotation coefficients of large degrees are proposed, studied numerically, and cross-validated. It is shown that the latter algorithm can be used for $nlesssim10^3$ in double precision, while the former algorithm was tested for large $n$ (up to $10^4$ in our experiments) and demonstrated better performance and accuracy compared to the FFT-based algorithm.

Recommendations

• A fast and stable method for rotating spherical harmonic expansions
• Efficient and accurate rotation of finite spherical harmonics expansions
• Recursions for the Computation of Multipole Translation and Rotation Coefficients for the 3-D Helmholtz Equation
• A fast transform for spherical harmonics
• Fast algorithms for spherical harmonic expansions. II.

Cites work

• A fast algorithm for particle simulations
• A fast and stable method for rotating spherical harmonic expansions
• A fast spherical harmonics transform algorithm
• A wideband fast multipole method for the Helmholtz equation in three dimensions
• Accurate Image Rotation Using Hermite Expansions
• Addition theorems for spherical wave functions
• Diagonal forms of translation operators for the Helmholtz equation in three dimensions
• scientific article; zbMATH DE number 3863589 (Why is no real title available?)
• scientific article; zbMATH DE number 3743203 (Why is no real title available?)
• scientific article; zbMATH DE number 3274385 (Why is no real title available?)
• scientific article; zbMATH DE number 3417938 (Why is no real title available?)
• Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations
• On the Partial Difference Equations of Mathematical Physics
• Recursions for the Computation of Multipole Translation and Rotation Coefficients for the 3-D Helmholtz Equation

Cited in

(4)

• Evaluation of the spherical harmonic coefficients for the external potential of a polyhedral body with linearly varying density
• Recursions for the Computation of Multipole Translation and Rotation Coefficients for the 3-D Helmholtz Equation
• On the role of the double Fourier sphere method in fast algorithms on \(\mathcal{SO} (3)\)
• A fast and stable method for rotating spherical harmonic expansions

This page was built for publication: Recursive computation of spherical harmonic rotation coefficients of large degree

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2950248)