Fast Fourier Transforms for Spherical Gauss-Laguerre Basis Functions

Jürgen Prestin,Christian Wülker

DOI: https://doi.org/10.48550/arXiv.1604.05140

2016-11-30

Abstract:Spherical Gauss-Laguerre (SGL) basis functions, i.e., normalized functions of the type $L_{n-l-1}^{(l + 1/2)} (r^2) r^{l} Y_{lm}(\vartheta,\varphi)$, $|m| \leq l < n \in \mathbb{N}$, $L_{n-l-1}^{(l + 1/2)}$ being a generalized Laguerre polynomial, $Y_{lm}$ a spherical harmonic, constitute an orthonormal basis of the space $L^{2}$ on $\mathbb{R}^{3}$ with Gaussian weight $\exp(-r^{2})$. These basis functions are used extensively, e.g., in biomolecular dynamic simulations. However, to the present, there is no reliable algorithm available to compute the Fourier coefficients of a function with respect to the SGL basis functions in a fast way. This paper presents such generalized FFTs. We start out from an SGL sampling theorem that permits an exact computation of the SGL Fourier expansion of bandlimited functions. By a separation-of-variables approach and the employment of a fast spherical Fourier transform, we then unveil a general class of fast SGL Fourier transforms. All of these algorithms have an asymptotic complexity of $\mathcal{O}(B^{4})$, $B$ being the respective bandlimit, while the number of sample points on $\mathbb{R}^{3}$ scales with $B^{3}$. This clearly improves the naive bound of $\mathcal{O}(B^{7})$. At the same time, our approach results in fast inverse transforms with the same asymptotic complexity as the forward transforms. We demonstrate the practical suitability of our algorithms in a numerical experiment. Notably, this is one of the first performances of generalized FFTs on a non-compact domain. We conclude with a discussion, including the layout of a true $\mathcal{O}(B^{3} \log^{2} B)$ fast SGL Fourier transform and inverse, and an outlook on future developments.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How to quickly calculate the Fourier coefficients of a function with respect to the spherical Gauss - Laguerre (SGL) basis functions**. Specifically, the author proposes a class of fast algorithms for efficiently calculating these Fourier coefficients and develops the corresponding fast inverse transform.

### Background and Problem Description

The spherical Gauss - Laguerre (SGL) basis functions are a class of orthonormalized functions in the form:
\[ H_{nlm}(r, \vartheta, \phi) = N_{nl} R_{nl}(r) Y_{lm}(\vartheta, \phi) \]
where:
- \( N_{nl} = \sqrt{\frac{2(n - l - 1)!}{\Gamma(n + 1/2)}} \)
- \( R_{nl}(r) = L^{(l + 1/2)}_{n - l - 1}(r^2) r^l \), \( L^{(l + 1/2)}_{n - l - 1} \) is the generalized Laguerre polynomial
- \( Y_{lm} \) is the spherical harmonic function

These basis functions are widely used in many fields, such as biomolecular dynamics simulations. However, currently, there is a lack of a reliable fast algorithm for calculating the Fourier coefficients of a function with respect to these basis functions. This is exactly the problem that this paper attempts to solve.

### Solution

The author proposes the fast algorithm through the following steps:

1. **SGL Sampling Theorem**: An exact discrete sampling formula for band - limited functions based on SGL basis functions.
2. **Separation of Variables Method**: Treat the radial part and the angular part separately.
3. **Fast Spherical Fourier Transform**: Use the fast spherical Fourier transform (FSFT) to accelerate the calculation of the angular part.
4. **Discrete R - Transform**: Introduce a new discrete R - transform to accelerate the calculation of the radial part.

### Algorithm Complexity

All the proposed fast algorithms have an asymptotic complexity of \( O(B^4) \), where \( B \) is the bandwidth limit, and the number of sample points is proportional to \( B^3 \). This is significantly better than the \( O(B^7) \) complexity of the naive algorithm.

### Application

These fast algorithms have important applications in the simulation of the biomolecular recognition process, such as protein - protein or protein - ligand docking problems.

### Summary

This paper proposes a class of fast Fourier transform algorithms for efficiently calculating the Fourier coefficients of a function with respect to the spherical Gauss - Laguerre basis functions. These algorithms not only improve the computational efficiency but also realize the application of the fast Fourier transform in the non - compact space \( \mathbb{R}^3 \).