Fourier Transform for Locally Integrable Functions with Rotational and Dilation Symmetry

: The Fourier transform for slowly increasing functions is deﬁned by the Parseval equation for tempered distributions. This deﬁnition was supplemented by a novel method of performing practical calculations by computing the Fourier transform for a suitably tempered function and then by integration by parts. The application of this method is illustrated both for the toy case, in which the function is integrable, so its Fourier transform can also be computed using the standard formula, and for the case of Coulomb-like potentials, which are only locally integrable functions. All of them have spherical symmetry, and two of them additionally have dilation symmetry. The proposed novel method does not violate these symmetries at any stage of the calculation.


Introduction
The Fourier transform (FT), defined for an integrable function, i.e., ψ ∈ L 1 (R 3 ) as [1] is a very useful concept that simplifies solutions of many problems in physical sciences, engineering, etc. Unfortunately, for several important physical quantities, the above definition cannot be applied, e.g., for the Coulomb potential V c (x) = 1 r , with r = |x| = √ x · x, which plays an important role in the description of electrostatic [2] and gravitational interactions, but evidently V c / ∈ L 1 (R 3 ). Not surprisingly, the standard definition of FT, when integrated over spherical angles in the position space R 3 , leads to an ill-defined integral where R + = {x ∈ R : x ≥ 0} and k = |k|. To eliminate this deficiency, a damping factor e −αr , with α > 0, is commonly added; for example, see [3][4][5], which converts the long-range Coulomb potential into the short-range Yukawa potential V α (x) Then, it is argued that the FT of the Coulomb potential V c can be defined by the limit α 0 of F {V α }, which is based on the analogy The advantage of such a procedure is that the spherically symmetric Coulomb potential is replaced by a spherically symmetric Yukawa potential, which means that rotational sym-metry remains unbroken. However, this does not apply to the dilation symmetry associated with the transformation x → x = λx, where the Coulomb potential is dilational covariant whereas the Yukawa potential explicitly breaks this symmetry. It seems obvious that any damping factor will break the dilation symmetry, so you have to invent a completely different computational procedure to preserve the dilation symmetry. This trouble of calculating the FT for the Coulomb potential can be avoided by saying that this FT is determined uniquely by the Poisson equation [6], but we did not pursue such reasoning further.
In this work, we present a novel method for calculating FT, which, when applied directly to the Coulomb potential, will reproduce the above result but without recourse to the Yukawa potential, so that both rotational and dilational symmetries remain intact at each stage of the procedure. Here, the dilational symmetry represents a physical symmetry that should be easily controlled at each stage of the computation. For more complicated cases, recovering such symmetry in the final result can lead to quite serious problems such as noncommutativity of limits, etc.). It will use the distributions with the test functions of the Schwartz class S (R 3 ), which due to the Parseval equation for the tempered distributions will lead to FT in the sense of distributions. For checking the consistency of this method, one may take a toy case W(x) = 1 x 2 ∈ L 1 (R 3 ), which gives the same FT, as the standard formula (1). Next, we consider the Coulomb-like potentials, which are slow-growing functions and thus have no standard FT, according to (2). First, we considered the softened Coulomb potential with a decaying exponential term V p (x) = 1 r e −c/r , with c > 0, [7], which vanishes for r → 0. Then, we studied the soft-core Coulomb potential U a (x) = 1 √ r 2 +a 2 , with an arbitrary parameter a ∈ R + [8], which is a special case of a generic soft-core Coulomb potential U a,q (x) = 1 (r q +a q ) 1/q , for q ≥ 1 [9][10][11], which is used for modeling potentials in atomic and molecular physics. At last, we considered the pure Coulomb potential V c , which has both rotational and dilation symmetry, unlike the previous two cases, which have only rotational symmetry. Our novel method is shown to be equivalent to the standard one, with the exponential damping factor, when the distributional derivatives are the standard derivatives-this is presented in the Appendix A. Finally, we mention how our method can be compared to the one used in [12], where dilation symmetry is used as a key step in the computation.

Calculation of Distributional FT-Methods
We begin by recalling some important properties of the Fourier transform and tempered decompositions that can be found in the vast literature, where the fundamental monograph by Gelfand and Shilov (see [12] pp. 71-74 and pp. 190-200) should be mentioned, but we can also select some items relevant to physics: [13][14][15]. The standard Fourier transform defined by (1) for integrable functions in L 1 (R 3 ) is a continuous and bounded function on R 3 , and it vanishes at infinity. For φ, ψ ∈ L 1 (R 3 ) one finds ψF {φ} ∈ L 1 (R 3 ) and φF {ψ} ∈ L 1 (R 3 ), so the Parseval theorem follows from Fubbini's theorem, where we used a compact notation for x, y ∈ R 3 . Evidently, this Parseval theorem is valid also for the Schwartz test function φ ∈ S (R 3 ), because S is a dense subspace of L p with 1 ≤ p < ∞. Moreover in this case both sides of (6) are the regular distributions S (R 3 ). This is a consequence of the fact that the Fourier transform is a bijection map F : S (R 3 ) → S (R 3 ) from the space S (R 3 ) unto itself, and hence F {φ} ∈ S (R 3 ). The situation changes for φ ∈ S (R 3 ) and L 1 loc F {ψ} cannot be given by the standard definition (1), so we cannot calculate the left-hand side in (6). In contrast, the right-hand side in (6) is a tempered distribution S (R 3 ) provided ψ is a slow-growing function, i.e., it has polynomial growth as |x| → ∞, which gives ψF {φ} ∈ S (R 3 ). This allows us to reinterpret (6) as a definition of distributional FT in the sense of distributions S (R 3 ).

Definition 1.
For a slow-growing locally integrable function ψ its distributional Fourier transform F D {ψ} is given by the relation for tempered distributions S (R 3 ) This equation will be called the Parseval equation for the distributional FT or the Parseval equation for short. Evidently if ψ ∈ L 1 (R 3 ), then its distributional Fourier transform F D {ψ} coincides with the standard one F {ψ} given by (1).
Before we move on to practical computations with the help of (8), we can slightly change the notation for the Schwartz test function by introducing where the bar denotes the complex conjugate. Thus, we may re-express (7) and (8) as with explicit notation for vectors x, k ∈ R 3 . Furthermore, in the following discussion, without loss of generality, we can drop the subscript in φ F = φ. The calculation of F D {ψ} is based on the analysis of the tempered distribution ψ[φ], and the standard procedure starts with a suitable modification of ψ[φ], where a damping factor is inserted under the integral sign. Such a calculation can be illustrated using the Coulomb potential ψ = V c / ∈ L 1 (R 3 ) as an example, where relations (9) and (10) look as However, for the tempered distribution we can introduce a damping factor e −αr , α > 0 and calculate using the complex conjugation of (1) for φ(x), where we can switch the order of integrations by Fubini's theorem, because V α (x) = V c (x)e −αr ∈ L 1 (R 3 ) and take F {V α } from (3). Then, for the last convergent integral, we can insert the limit α 0 under the sign of the integral, thus reaching the left-hand side of the Parseval equation in the form (10) and yielding the distributional FT for the Coulomb potential Accordingly, the ad hoc calculations presented in the introduction are closely related to the standard calculation of the distributional FT. Therefore, the previous drawbacks mentioned in the introduction also apply to this method.
We propose a different computational procedure, (a comparison of these two computational methods is presented in Appendix A), that relies on specific properties of Schwartz where n ∈ N and M ∈ R + were chosen arbitrarily [14] (also see pp. 16-17 in [12]). From the definition (1), one checks the formula for FT: k is the Laplace operator in the space of k ∈ R 3 . All this leads to the tempered distribution where we can take n ∈ N 0 and M ∈ R + arbitrarily, provided that Then, one can integrate by parts with respect to k, which gives the final equivalent form, which can be compared with the Parseval equation for the distributional FT (10) where and it has a polynomial growth for |x| → ∞ and ψ/ (r 2 + M 2 ) n ∈ L 1 (R 3 ) with appropriate choice of n and M, then its distributional FT is given in the sense of distributions S (R 3 ) by In case of rotational symmetry one has ψ(x) = ψ(r), so the integral over the unit sphere embedded in R 3 can be easily performed done leading to the result, which depends only on the radial variable k = |k| This in turn allows us to move on to integrals with respect to k ∈ R 3 , which can be expressed as the radial integral where the radial FT F rad {φ} is defined as the integral over the unit sphere embedded in R 3 where dω k is the hypersurface element on the unit sphere, and Ω 3 = 4π is the hypersurface area of the unit sphere in R 3 (see [12] p. 71). (Evidently F rad {φ}(k) is the mean value of F rad {φ}(k) on the sphere of radius k (see [12] p. 71), and one may express (18) by means of the mean value of φ(x) on the sphere of radius k which agrees with the FT for the spherically symmetric functions.) In practical calculations, one must first make the definite integral with respect to r, which gives a well-defined function of k. Then, integrating by parts with respect to k one finally obtains thus indicating that the distributional FT F D {ψ}(k) is also rotationally symmetric.
Although the tempered distribution ψ[φ] is defined in (13) as an expression independent of M, its final form explicitly includes the parameter M. It is therefore desirable to see how this presence effectively cancels itself, and this can be done by computing the derivative of the last line in (13) with respect to the parameter M where in (16) one may push the derivative under the integral sign in F ψ/(r 2 + M 2 ) n that is given by a convergent integral. These two tempered distributions S (R 3 ) are equal to each other, which can be easily checked by integrating by parts in the second case, where ∇ 2 k can be inserted under the integral sign F ψ/(r 2 + M 2 ) n+1 , which is given by the convergent integral. This brings us to the expected result ∂ ∂M 2 ψ[φ] = 0, but we discovered the crucial role played by integration by parts for checking the M-independence of ψ [φ].
If one imposes the dilation transformation on the function ψ(r) → ψ λ (r) = ψ(λr), then the radial Parseval equation (17) for the distribution ψ λ [φ] becomes where we changed the integral variables: r = ξλ −1 , k = pλ and the parameter M = M 1 λ. Thus, in the sense of distributions, we have the distributional FT for the dilation transformed function F D {ψ λ }(k) = λ −3 F D {ψ}(kλ −1 ). However, if the function ψ has dilation symmetry, then its distributional FT also has dilation symmetry This shows that our novel method of computing FT for functions in R 3 does not conflict with rotational and dilational symmetry.
Example 3. One may check the above tools for a toy function W(x) = 1 r 2 , which is both rotationally symmetric and dilation covariant W λ (x) = W(λx) = λ −2 W(x) and W ∈ L 1 (R 3 ), so its FT can be also computed from the standard definition (1) Then, we used the Parseval equation (17) with n = 1, for the functional W[φ], which gives the distribution, where first we performed the integral with respect to r, and then we integrated by parts with respect to k , so it agrees with the standard FT, as expected. Moreover, the dilation symmetry relation for FT is in agreement with the Formula (23).

Results for Coulomb-like Potentials
Next, we calculated FT for Coulomb-type potentials in turn: the smooth potential V P (x) = 1 r e −b 2 /r , the soft-core potential U a (x) = 1 √ r 2 +a 2 for a ∈ R + , and the plain Coulomb potential V c (x) = 1 r . All of them are locally integrable in R 3 ; thus, one can define the respective tempered distributions, for test function φ ∈ S (R 3 ), , we may take the simplest form of the Parseval equation (17) with n = 1 and M = 0, because the exponential factor smooths the vicinity of r = 0; then, we computed the integral with respect to r using formula (3.957.1) in [16] where K n (z) = K −n (z) is the modified Bessel function of the first kind of the n-th order. The integration by parts with respect to k introduces no boundary terms and one needs to take the second order derivative for the modified Bessel function Therefore, by analogy with (20), we obtain differential equation with the distributional FT This FT agrees with the corresponding result in [7], provided we use (iz) = (z) and take into account the slightly different definition of FT therein. We also note that in [7] the damping factor e −κr with κ > 0 was added before computing the FT, and the limit of κ 0 was taken for the result obtained. Then, we calculated the distributional FT for the soft-core Coulomb potential U a (x) = 1 √ r 2 +a 2 , where we can take (17) with n = 1 and M = 0, which gives the distribution that allows to compute the integral with respect to r by using formula (3.754 The integration by parts with respect to k introduces no boundary terms; therefore, we obtained with the distributional FT For this soft-core Coulomb potential we may make alternative calculations that begin with separation of the leading term (for r a) as U a (x) = V c (x) + U a sub (x), where the subtracted soft-core potential U a sub (x) follows from the relation This leads to the subtracted functional where we integrated by parts twice with respect to k, inserting the differentiation under the sign of the integral. Then, the integrals with respect to r and µ can be calculated successively, using formula (3.754.3) in [16] a 0 R + µ r sin(kr) giving the final form of the subtracted functional From the formula U a (x) = V c (x) + U a sub (x), one directly obtains the formula for the functionals V c [φ] = U a [φ] − U a sub [φ] and next the formula for FTs In this way, we obtained the distributional FT for the Coulomb potential, which agrees with the expression obtained earlier by different methods. This clearly shows that a FT that is symmetric under dilation can appear as a difference of two FTs that are not symmetric under dilation. At last, for checking consistency, we calculated directly the distributional FT for the Coulomb potential, so we took (17) with n = 1 and M > 0, and, as before, first we needed to calculate the integral with respect to r. Using formula (8.217.1) in [16], we obtained where Ei(z) is the exponential integral function and then we may integrate by parts with respect to k (Above we have used the differential properties of F ± , that are discussed in Appendix B). Thus, the distributional FT is F D {V c }(k) = 4π k 2 , which agrees with the previous result (35). The distributional FT for the pure Coulomb potential V c can be used to prove the Poisson equation for V c in the sense of distribution S (R 3 ). To demonstrate this, we calculates a new distribution V c [∇ 2 φ] for test function φ ∈ S (R 3 ), which allows for performing distributional derivatives, applying the Parseval equation for tempered distributions and the definition of Dirac delta singular function, successively Thus, we may integrate by parts, which leads to distributional derivatives,

Discussion and Further Research
The proposed novel method of computing FT for slowly increasing functions can be applied to functions that are locally integrable and can be tempered with a polynomial of a finite order. The most troublesome part of this procedure is computation of a convergent improper integral, which may additionally depend on some arbitrary auxiliary parameterin this study, we used the formulas for integrals from the table of integrals [16].
Due to the rotational symmetry of the analyzed functions, we obtained radial integrals and radial test functions. The dilation symmetry manifested itself in a less evident way but was still not broken at any stage of the calculation.
For the softened Coulomb potential V p (x) and the soft-core Coulomb potential U a (x), due to a smooth dependence on r → 0, the Parseval equation may be taken with M = 0, which simplifies our computation and leads to differentiable auxiliary functions. Then, the distributional FT for each potential becomes the second-order derivative of the respective auxiliary function. Both the auxiliary functions and the FTs depend on the modified Bessel functions. However, the plain Coulomb potential V c (x) is singular for r → 0, so we needed to take M = 0, and further steps of computation are more difficult. We obtained the auxiliary function that depends on the exponential integral function of Ei(z). However, it appeared only in an intermediate step of the computation, while in the final result this special function disappeared.
Furthermore, we showed the usefulness of calculating the FT for a subtracted function, U a sub (x), that is already an integrable function, at the cost of introducing an additional proper integral with respect to an auxiliary parameter, in our case µ.
In general, it can be said that the increased rigor of the method comes at a price, as the calculations become somewhat more complicated. The applicability of our new method, according to the lemma in Section 2, is limited to locally integrable functions on R 3 that grow slowly as |x| → ∞. However, these functions can also contain additional variables that are not integrated when computing FT. In physics, such an additional variable might be a time coordinate. For example, in quantum field theory, our method can be used for calculation of the Wightman functions [17], for a free scalar field. The resulting expression has the Lorentz symmetry and contains distributional derivatives-this will be published soon in a separate article. This new result deviates from previous approaches, where one either uses mathematically inaccurate tricks, as in [18], or introduces damping factors that break an important physical symmetry (here, the Lorentz symmetry) [19] or uses Lorentz symmetry to analyze divergent integrals [20]. The main feature of the new expressions for Wightman functions is the appearance of distributional derivatives, which has not been used before. We plan to apply our novel approach to FT calculations in the context of light front quantization [21], where existing previous results for Wightman functions [22,23] should be checked and/or improved by analysis from a different point of view. Our computational method can be applied to other mathematical problems arising in quantum mechanics that require the computation of FT for functions that are not in L 1 (R n ); for example, regarding the time-evolution operator in the free Schrödinger equation, see [24].
We also need to comment on the FT calculated for the spherically symmetric generalized function r λ for λ = −n, −n − 2, . . . in R n in [12]. A key part of this calculation is the dilation covariance (homogeneity), so it cannot be used for the case of Coulomb potentials considered in this study. Therefore, we plan to apply our novel procedure to compute the FT for the above generalized function in a separate study, which should provide a basis for comparing the two methods.
which can be used for doing partial integration in the linear functionals where primes denote derivatives with respect to k.