Integral Involving Bessel Functions Arising in Propagation Phenomena

: An new integral identity involving the product of two modiﬁed Bessel functions K 0 and K 1 , previously unreported, is presented and proved.


Introduction
Bessel functions are ubiquitous in mathematical physics and engineering. Yet they continue to surprise us with amazing new properties that, to the best of our knowledge, seem not to have been recorded anywhere, at least not in classic literature on the subject [1][2][3]. In this paper we state and prove one such property. We also speculate on a possible application of this new property to an eventual theory of quantum gravity.
Concerning our use of Bessel functions, we follow the notations and conventions of ref. [2].

Statement and Proof of Theorem
Theorem 1. The following identity holds: where x, y ∈ R with y = 0.
Proof. Let us consider the contour integral in the complex ξ-plane along a contour C to be specified presently. For a fixed x, the change of variables To begin with, we assume y > 0. Let the contour C extend along the straight segment from −x + iy to x + iy, while closing on the upper half-plane along a semicircle centred at the point (0, y) with radius R. By Cauchy's theorem where the sum extends over all the residues enclosed by C.
We will need the asymptotics of K 0 (z) and K 1 (z) on the complex plane. As z → 0 we have while for z → ∞ Away from the negative real axis in the complex plane, the functions K 0 and K 1 are analytic. In Equation (7), K 0 (|z|) and K 1 ( y 2 + (z + x) 2 ) both vanish as R → ∞. Then the sum over residues in Equation (5) reduces to the single residue located at In Equation (6), in a neighbourhood of z 0 the function on the right-hand side of (5) behaves as Altogether which proves the statement made in Equation (1) when y > 0. An analogous reasoning holds when y < 0, the contour C now being closed by a semicircle in the lower half plane and the residue being located at z 0 = −x − iy, which introduces an extra minus sign. This completes the proof.

Quantum Gravity?
The authors' own interest in the subject arose out of contact [4][5][6][7] with problems in the theory of heat conduction [8] and quantum field theory [9]. Without aspirations to completeness (the subject is too vast to summarise) let us briefly mention one such occurrence of Bessel functions.
In the quantum field theory of a scalar field with mass m in Euclidean spacetime R d [9] one is interested in the Feynman propagator G d . The latter is a Green function for the Klein-Gordon operator on R d satisfying a specific set of boundary conditions. With x denoting a vector in R d and ||x|| = √ x 2 its Euclidean norm, it turns out that We see that the propagator G d depends on the vector argument x ∈ R d only through the norm ||x||. The latter is a real number that we may denote by x. Further setting m = 1 for convenience, we simplify Equation (12) to a function of a single real variable x := ||x||: The Bessel functions K 1−d/2 of lowest integer order arise when d = 2 and d = 4, respectively. When d = 2 we have that the Feynman propagator G 2 (x) is proportional to the function K 0 (x): When d = 4, the identity The new identity (1) is thus reminiscent of the self-reproducing property of integral kernels (such as the propagators G 2 and G 4 ), yet it is definitely different, because it involves propagation in different dimensions. In fact we speculate that it might be related to propagation phenomena in an (eventual) theory of quantum gravity. In order to see how this might come about, let us recall the following facts [10,11].
The propagators (14) and (15) have been obtained under the assumption that spacetime remains a continuum all the way down to the shortest classical length, namely zero length. Now there are good theoretical reasons to believe that the continuous nature of spacetime breaks down when one reaches lengths of the order of the Planck scale L P = √h G/c 3 . Hereh is Planck's quantum of action, G is Newton's gravitational constant, and c is the speed of light. The numerical value of L P is roughly 10 −33 cm, well beyond what current technology can probe. All this notwithstanding, it has been convincingly argued [10,11] that an effective, or low-energy, manifestation of quantum gravity effects on the propagator functions G d (x) can be captured by the simple replacement within the argument of G d . In other words, a true propagator that encodes overall quantum-gravity effects will no longer be given by G d (x) = G d ( √ x 2 ), but instead by G d ( x 2 + L 2 P ). Thus even after setting x = 0, spacetime points will continue to be separated by a nonzero length L P . Among other desirable properties, this has the welcome effect that the singularities present in the propagators (14) and (15) as x → 0 are smoothed out. For the sake of simplicity let us assume choosing units such that L P = 1. Then the quantum-gravity corrected propagator functions (14) and (15) become The above corrected propagators duly reduce to their counterparts (14) and (15) for large x, while no longer exhibiting a singular behaviour at the origin. This correction is denoted by the superindex QG, standing for quantum gravity. As a rule one generally expects quantum effects to mollify (at least some of) the singular behaviour that classical theories may exhibit.
Altogether, substituting Equations (14) and (17) into Equation (1) and setting y = 1 = L P in the latter, the identity we have just proved reads in terms of scalar propagation functions in Euclidean space. Although reminiscent of the self-reproducing property of Feynman propagators, it is intriguing for several reasons. First, it involves different dimensions (2 and 4). Second, it links propagators in the presence of quantum-gravity effects to propagators in the absence of quantum-gravity effects. Last but not least, it requires that y be nonvanishing, so it suggests setting y = L P -a key feature of quantum-gravity theories is the existence of a quantum of length. One possible reading of Equation (18) is the following: starting from the 2-dimensional, quantum-gravity free propagator G 2 , it suffices to compute its convolution with the 4-dimensional, quantum-gravity corrected propagator G (QG) 4 , in order to obtain the quantum-gravity corrected propagator G (QG) 2 in two dimensions. We hope to report on this issue in the future [12].
Author Contributions: All authors contributed equally.