Evolution of Superoscillations in the Dirac Field

Superoscillating functions are band-limited functions that can oscillate faster than their fastest Fourier component. The study of the evolution of superoscillations as initial datum of field equations requires the notion of supershift, which generalizes the concept of superoscillations. The present paper has a dual purpose. The first one is to give an updated and self-contained explanation of the strategy to study the evolution of superoscillations by referring to the quantum-mechanical Schrödinger equation and its variations. The second purpose is to treat the Dirac equation in relativistic quantum theory. The treatment of the evolution of superoscillations for the Dirac equation can be deduced by recent results on the Klein–Gordon equation, but further additional considerations are in order, which are fully described in this paper.


Introduction to Superoscillations and the Supershift Property
Superoscillating functions are band-limited functions that can oscillate faster than their fastest Fourier component. These functions (or sequences) appear in weak values in quantum mechanics (see [3,12,14,26]) and in several other fields of science and technology, as it will be mentioned in the sequel. A natural problem, suggested by Aharonov, is to study the evolution of weak-values-superoscillations as initial datum of Schrödinger equation or as initial condition of other quantum field equations such are Klein-Gordon, see [10], or Dirac equation which is treated in this paper.
In the last decade a systematic study of superoscillating functions has been carried out also from the mathematical point of view. In fact, the rigorous treatment of the evolution of superoscillations, as initial datum of the Schrödinger equation, required some sophisticated mathematical tools, as it has been shown for example in [1, 2, 4-8, 16, 17, 27-29] and partially summarized in [9]. The prototypical function arising in Aharonov's weak values is given by which, using Euler's identities and Newton's binomial, can be written as where a > 1 and If we fix x ∈ ℝ and we let n go to infinity, we obtain that The term superoscillations comes from the fact that the frequencies (1 − 2j∕n) in (2) are in modulus bounded by 1 but the frequency a, in the limit function (4), can be arbitrarily large. The superoscillatory phenomenon appears because thanks to the form the coefficients C j (n, a) the plane waves e i (1−2j∕n)x interfere in such a way that they produce the plane wave e iax , where a can be very large with respect to one.
Superoscillations occur also when n is a fixed number, not only when n goes to infinity, simply the region of superoscillation is a small interval that contains the origin. One can clearly see that if x is very small compared to n it is x n ∼ 0 so in this limit we have We illustrate in the figure the superoscillations for fixed n and a plotting the logarithm of the real part of F n . The logarithm is necessary to re-scale the function in other to better see the superoscillations ( Fig. 1).
We note that there is no conflict with the theory of Fourier series because the sequence F n (x, a) is not a Fourier sequence even though it looks like. In fact, the

3
coefficients C j (n, a) and the wave numbers k j (n) ∶= 1 − 2j∕n depend both on j and on n while for Fourier series the numbers C j (n) and k j (n) depend on j and but not on n. Sequences such as F n (x, a) are called generalized Fourier sequences and can produce superoscillations while Fourier series cannot. In order to understand the meaning of the evolution of superoscillations as initial condition of a given quantum field equation we need the following interpretation of the generalized Fourier sequence F n (x, a) . Consider x as a parameter and let us focus our attention on the points (1 − 2j∕n) . For the sake of clarity we set and with this position we write the sequence F n (x, a) as Now consider x as a fixed parameter and recall we have set e iax = (a, x) . Then the limit F n (x, a) → (a, x) , as n → ∞ means that we compute the function ↦ ( , x) in infinitely many points 1 − 2j∕n that belong to the interval [−1, 1] (band limited frequencies) and we determine the value of ( , x) in the point = a >> 1 . When we consider the evolution of F n (x, a) as initial value the free particle using Schrödinger equation we obtain the solution We define the following function of where t and x are considered parameters, so that studying the evolution of superoscillations is tantamount to studying the limit of the solution for n → ∞ . Since for the particular case of the free particle the solution of Schrödinger equation is written in terms of the exponential function we still speak of superoscillations, but in the general case of nontrivial potentials the functions ( , t, x) are not necessarily the exponential function.
In this case there arises the natural notion of supershift of the solution where we compute the function ↦ ( , t, x) in infinitely many points 1 − 2j∕n that belong to the interval [−1, 1] , and we determine the value of ( , t, x) in the point = a >> 1 . For example for the quantum harmonic oscillator the function ( , t, x) can be computed by elementary functions and it is given by As shown above, superoscillations are a particular case of supershift, which is the property to be proven when we study the evolution of superoscillations as initial datum of given field equation.
A different class of superoscillating functions with respect to the ones considered in this paper can be found in [23]. In antenna theory the phenomenon of superoscillations was discovered by Toraldo di Francia in [36] (see also the paper of Berry [20]). The literature on superoscillations is large, and without claiming completeness we mention the papers [19,21,22,24,25,30,31,33,35]. More information can be found in the introductory papers [11,13,15,32].
Moreover, in the recent paper Roadmap on superoscillations, see [18], some of the most important achievements in the applications of superoscillations are well explained by the leading experts in this field. In fact, the main aim of the paper [18] is to give an exposition of various aspects of superoscillations that are accessible to non experts in this field. There are contributions in: optical superoscillatory focusing and imaging technologies, superoscillatory interference for super-resolution telescope, an antenna array approach to superoscillations, applications of superoscillations in ultrafast optics, superoscillations in magnetic holography and acoustics, superoscillations and information theory. The paper also contains many illustrations to clarify the experiments and it can be freely downloaded from the web. This paper has two main purposes. The first one is to explain and motivate the general strategy to study the evolution of superoscillations in the context of quantum mechanics. The second one is to study superoscillations as initial datum of the Dirac equation. The paper is organized as follows. In Sect. 2 we formulate and state our main results in precise mathematical term for the Dirac equation. In Sect. 3 we discuss the strategy to demonstrate our results, whose full proofs are then given in Sect. 4.

The Main Results for the Dirac Equation
In this section we formulate in precise mathematical terms the problem of the evolution of superoscillating initial datum for the Dirac equation. We summarize the main result on the supershift property for the Dirac equation, but we postpone the proofs to Sect. 4 of the paper. The sequence defined in (2) extends to the sequence of entire function F n (z, a) , for z ∈ ℂ . The tools to study the evolution of superoscillations are convolution operators that naturally act on spaces of entire functions that contains the holomorphic extensions F n (z, a) , so we recall some notions and results on entire functions. Let f be a non-constant entire function of a complex variable z. We define for r ≥ 0 . The non-negative real number defined by is called the order of f. If is finite then f is said to be of finite order and if = ∞ the function f is said to be of infinite order. In case f is of finite order we define the nonnegative real number which is refereed to as the type of f. More to the point, if ∈ (0, ∞) we say that f is of normal type, while we say that f is of minimal type if = 0 and of maximal type if = ∞ . The constant functions are said to be of minimal type of order zero. Let p ≥ 1 , we denote by A p the space of entire functions with either order lower than p or order equal to p and finite type. It consists of functions f for which there exist constants B, C > 0 such that The extension of the functions F n (x, a) defined in (2) are exponentially bounded and more general sequences that generalize the F n (x, a) can be defined as follows: with coefficients c j (n) ∈ ℂ and k j (n) ∈ ℝ . We say that Y n is a superoscillating sequence (or function), if there exists some a ∈ ℝ , with |a| > 1 , such that Foundations of Physics (2020) 50:1356-1375 Superoscillating functions as F n do not only converge uniformly on the compact sets of ℂ but they converge in the space A 1 . More in general we recall that the We now come to the definition of supershift property. Let I ⊆ ℝ be an interval with be a continuous function on I . We set and consider a sequence of points ( k,n ) such that and a sequence of complex numbers (c k (n)) for k = 0, … , n and n = 0, … , +∞ . We define the functions If we have that lim n→∞ n (t, x) = a (t, x) for some a ∈ I and with |a| > 1 , we say that the function → (t, x) , for t and x fixed, admits a supershift in . Finally we state our results for the Dirac equation. Let be the Gamma matrices and let + 0 (x) and − 0 (x) represent the wave functions for the positive and negative energy level. Now consider, for a, b ∈ ℝ and m > 0 , the Cauchy problem where and denote by 1,a (t, x) , 2,b (t, x) , 3,b (t, x) and 4,b (t, x) the components of the solution (t, x) . Suppose that + 0,n (x) and − 0,n (x) are the superoscillating initial data of the Cauchy problem where Then the solution of (11) can be written in terms of the solution of the Cauchy problem (9) as and the supershift property holds, that is The precise mathematical proof of the above result is not trivial and all the details are given in Sect. 4.

Supershift and Techniques Based on Infinite Order Differential Operators
In this section we recall some facts leading to the definition of supershift and we summarize the main mathematical tools used to study superoscillations. From the prototypical example described in the previous section one can extend the definition of superoscillations to the case of several variables. To this aim, we give the definition of generalized Fourier sequence in several variables, which is a variation of the definition introduced in [8].

Definition 3.1 (Generalized Fourier sequence in several variables)
For m ∈ ℕ we assume that (x 1 , … , x m ) ∈ ℝ m . Let ( k,j (n)) , k = 0, … , n for n = 0, … , ∞ , be realvalued sequences for j = 1, … , m . We call generalized Fourier sequence in several variables a sequence of the form where (c k (n)) , k = 0, … , n , n = 0, … , +∞ is a real-valued sequence.  and there exists a compact subset of ℝ m , which will be called a superoscillation set, on which f n converges uniformly to ∏ m j=1 e ix j g j , where |g j | > 1 for j = 1, … , m.
A very important particular case of generalized Fourier sequence in several variables is the following where C k (n, a) are given by (3) and q j ∈ ℕ , for j = 1, … , m.
From an historical point of view the main problem, which was first suggested by Aharonov, was to study the evolution of superoscillations by Schrödinger equation for the force free field with rigorous mathematical tools. This example was the starting point to understand the most natural functional setting in which one can study the evolution of superoscillations. We have already mentioned the solution of the Schrödinger equation with superoscillatory initial datum for the free particle is given by The crucial problem is to show that during the evolution, as n goes to infinity, the superoscillatory behavior of the limit function persists. It turns out that the limit is uniform on compact sets in ℝ 2 . This means that the limit function e iax−ia 2 t has a superoscillatory behavior. The main tools used to compute rigorously such a limit are convolution operators acting on spaces of entire functions. For such operators a crucial fact is to prove a continuity theorem. Specifically, for the free particle the convolution operator associated with this problem is so that the evolution problem can be written as for every x ∈ ℝ and t ∈ ℝ . Clearly, the function in (16) is a generalized Fourier sequence in two variables, and for a > 1 it is a superoscillatory function.
Inspired by the example of the free particle, if in Eq. (15) one replaces the second derivative with respect to the space variable x by the partial derivative p x (x, t) of order p ∈ ℕ , it is possible to prove that we still have a superoscillating function. In this case we have (p is even) which converges, for n → ∞ , and uniformly on compact sets, to (x, t) = e iax−ia p t . More in general, there is a very large class of superoscillatory functions that can be generated stating from convolution equations of the type where the second derivative with respect to time in (15) has been replaced by a suitable convolution operator where the complex coefficients a satisfy give growth conditions.
For the above class of differential equations, or more in general for convolution equations, the solutions remain in the class of generalized Fourier sequences, but when we take into account the Schrödinger equation with a nonconstant potential it turns out that the solution is not a generalized Fourier series in several variable. For instance, when we consider a linear potential, which means that the evolution takes place in a constant electric field, the solution of the Cauchy problem is given by which is not of the form (14). So, by computing the limit for n → ∞ , we obtain Hence, the class of superoscillations proves too narrow. In order to overcome this limitation, the notion of supershift has been introduced.

3
Foundations of Physics (2020) 50:1356-1375 and we consider a sequence of points ( k,n ) such that We define the functions where (c k (n)) is a sequence of complex numbers for k = 0, … , n and n = 0, … , +∞ . If lim n→∞ n (t, x) = a (t, x) for some a ∈ I with |a| > 1 , we say that the function → (t, x) , for t and x fixed, admits a supershift in .
Let us stress that the term supershift comes from the fact that I can be an arbitrarily large (it can also be ℝ ) and so also a can be arbitrarily far away from the interval [−1, 1] where the function is computed, see (20). In the case of the sequence (F n ) in (2), we have k,n = 1 − 2k∕n , k,n (t, x) = e i(1−2k∕n)x and c k (n) are the coefficients C k (n, a) defined in (3).

Schrödinger-Like Equations
By using the notion of supershift, it is now possible to formulate in precise mathematical terms the problem of the evolution of superoscillations for a quantum mechanical system described by the Hamiltonian operator H . We consider the evolution problem and we denote by n (t, x) the solution of the Cauchy problem (21), by a (t, x) we indicate the solution of the Schrödinger equation in (21) where we replace F n by the initial datum e iax . Given the solutions n (t, x) and a (t, x) , owing to the linearity of Schrödinger equation, the function n (t, x) can be written as The persistence of superoscillations means that the supershift property holds. Hence, we have to show that or, equivalently, that the function → (t, x) admits a supershift in the sense Definition 3.3.
When the Green function of the time-dependent Schrödinger equation is known, there is a standard procedure to study the evolution problem (21) and to determine (t, x) ∶= ( , t, x), k,n ∈ [−1, 1] for k = 0, … , n and n = 0, … , +∞.  (23) converges. This fact can be assured by regularizing the integral or by considering the integral defined as a Fresnel integral. Here complex analysis plays an important role in the definition of the integral.
(II) Once that point (I) is enforced, we have to assure that the function is analytic and then we need suitable estimates, which allows us to exchange the series with the integral in such a way that we obtain This is quite delicate fact to handle because the integral is computed, in general, on an unbounded domain Ω.
(III) Finally, suppose formula (24) holds, we set where and by introducing the auxiliary complex variable we represent the function a ↦ ∫ Ω e iay K(x, y, t)dy where U(x, t, ) is the convolution operator defined by acting on e ia and | =0 denotes the restriction to = 0. The above steps make sense when we can prove that the operator U(x, t, ) acts continuously on a space of entire functions that contains the exponential function. For example, if we show that acts continuously, then we have which in turn gives the sought-after supershift property. This follows because the sequence F n ( , a) converges to e ia in A 1 and in particular it converges uniformly on the compact set in ℂ.
The above procedure is very general and it works also for different types of evolution equations recently it has been applied to the Klein-Gordon equation (see [10]). In the next section, by taking advantage of some preliminary results presented in [10], we study the supershift property for the Dirac field.

Proofs of the Main Results for the Dirac Equation
In this section we recall a well known representation formula for the solution of Klein-Gordon equation in ℝ 3 expressed by smooth functions like the Bessel function J 1 . This formula is a consequence of the knowledge of the Green function and it expresses the solution of the Cauchy problem Precisely, we consider y = (y 1 , y 2 , y 3 ) ∈ ℝ 3 with y written in spherical coordinates, i.e., y 1 = sin cos , y 2 = sin sin , y 3 = cos . We will set d 3 y = r 2 drdΩ y , where dΩ y = sin d d , for ∈ [0, ∞) , ∈ [0, ] , ∈ [0, 2 ]. Theorem 4.1 Let f ∈ C 3 (ℝ 3 ) and g ∈ C 2 (ℝ 3 ). Then the solution of the Cauchy problem (26) is given by where J 1 is the Bessel function, the spherical coordinates are given by d 3 y = r 2 drdΩ y , and B t (0) is the ball with radius t centered at the origin. Moreover, belongs to in C 2 (ℝ 4 ).

Conclusions
Superoscillations are an important notion in quantum mechanics, as they appear in Aharonov's theory of weak values. The study of superoscillations requires sophisticated mathematical tools in order describe their evolution in time according to Schrödinger equation or other relativistic field equations. The mathematical methods to investigate the evolution of superoscillations have improved in the last ten years, and in this paper we have offered a summary of these methods by using Schrödinger equation as a model. In particular, we explicitly explained why we need to generalize the notion superoscillation with the notion of supershift. That put us in a position to develop our main technical result, namely the treatment of the evolution of superoscillation via the Dirac equation.
Funding Open access funding provided by Politecnico di Milano within the CRUI-CARE Agreement.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/ licenses/by/4.0/.