Quantum Mechanics and Control Using Fractional Calculus: A Study of the Shutter Problem for Fractional Quantum Fields

The ‘diffraction in space’ and the ‘diffraction in time’ phenomena are considered in regard to a continuously open, and a closed shutter that is opened at an instant in time, respectively. The purpose of this is to provide a background to the principal theme of this article, which is to extend the ‘quantum shutter problem’ for the case when the wave function is determined by the fundamental solution to a partial differential equation with a fractional derivative of space or of time. This involves the development of Green’s function solutions for the space- and time-fractional Schrödinger equation and the time-fractional Klein–Gordon equation (for the semi-relativistic case). In each case, the focus is on the development of primarily one-dimensional solutions, subject to an initial condition which controls the dynamical behaviour of the wave function. Coupled with variations in the fractional order of the fractional derivatives, illustrative example results are provided that are based on presenting space-time maps of the wave function; specifically, the probability density of the wave function. In this context, the paper provides a case study of fractional quantum mechanics and control using fractional calculus.


Introduction
In 1952, Marcos Moshinsky published a seminal article entitled 'Transient Phenomena in Quantum Mechanics: Diffraction in Time' [1]. Using a one-dimensional model, Moshinsky considered the case when a beam of non-relativistic electrons (travelling from left to right in the negative half-space) is incident upon a shutter (placed at x = 0) so that they cannot propagate through to the positive half-space. At a time t = 0, the shutter in opened instantaneously so that the electron beam can propagate through from the negative to the positive half-space. The problem is to evaluate the probability density of the resulting wave-function (determined by the Schrödinger equation) in the positive half-space as a function of both space and time. In other words, if u(x, t) is taken to be the electron wave-field, then what is the (x, t) map, for x > 0 and t > 0 of | u(x, t) | 2 when the shutter is opened at t = 0?
The space-time characteristics of the probability density are governed by the initial condition u 0 (x) ≡ u(x, 0) of the field at t = 0, where u 0 (x) is taken to represent of the electron beam in the negative half-space. In this context, the study undertaken by Moshinsky was an example of a control problem in quantum mechanics. The analysis undertaken by him involved developing a Green's function solution to the homogeneous time-dependent Schrödinger equation [2]. For atomic units (where the Dirac constant and the electron mass are set to 1), this equation is identical to the homogeneous time-dependent diffusion equation for a Diffusivity of 1/2 but with imaginary time. Consequently, the Green's function for the time-dependent diffusion equation (for a Diffusivity of 1/2), which includes the function exp(−x 2 /2t), yields a solution to the Schrödinger equation, which is characterised by the function exp(ix 2 /2t).
For a fixed value of t, this is a 'chirp function' of space (a linear frequency modulated chirp with a 'chirp parameter' of 1/2t). It is therefore analogous to the convolution kernel associated with Fresnel diffraction [3] which is characterised by a quadratic phase factor [4]. However, for a fixed value of x, as t → 0, the function exp(−x 2 /2t) oscillates at an increasingly high frequency, characterised by an inverse phase factor of time. For this reason, the phenomenon has been referred to as 'Diffraction in Time' (or even 'Time Scattering') [5,6]. The principal point here is that, unlike the diffusion equation, the same Green's function solution applied to the Schrödinger equation yields a fundamentally different result, which is compounded in the oscillations of the wave-field in space but also, and, more significantly, in time.
The exact nature of these Fresnel-type oscillations depends on the initial condition, but the fact that such oscillations in time occur at all is due entirely to the characteristics of the Green's function alone. This result illustrates a fundamental physical phenomenon, namely that a stream of electrons described by a plane wave-function with a fixed frequency, will spontaneously undergo rapid oscillations in time as soon as the shutter is opened. While such oscillations decrease in frequency (and amplitude) over time, the result shows that the electron wave-field spontaneously oscillates, due to the instantaneous opening of the shutter alone (with no other physical interactions). Such a transient phenomenon is now recognised as being a ubiquitous phenomenon in quantum dynamics [7,8], and, some fifty years after Moshinsky proposed such a dynamic process on theoretical grounds alone, experimental confirmation of the phenomenon was provided in 1996 using ultra-cold atoms [9].
The experimental verification of the 'Quantum Shutter Problem', as it has come to be known, has led to a number of publications which have looked into the evaluation of the problem in three dimensions [10], general formulations [11], new theories of diffraction in time [12], and an extension of the problem to include relativistic effects for both spin-less and half-spin particle beams [13].
There are numerous applications to which fractional calculus is becoming a central theme [14]. In this respect, its application to quantum mechanics has and continues to be explored. Many relatively recent publications (e.g., [15][16][17]) are paving the way for a future study of fractal geometry in quantum mechanics, for example, using fractional calculus. In this context, the material reported in this paper belongs to the same field of study. It is focusing on the quantum shutter problem and extends it further to include a study of the following: A schematic diagram, which is referred to throughout the paper, is given in Figure 1. Known generically as the 'shutter problem', in this paper, we consider the diffraction-type patterns that occur when a fractional partial differential equation is taken to govern the quantum field that is generated.

Incident Wavefield
(The initial condition at t = 0) Shutter (opened instantaneously at t = 0) What is the wavefield that is created when the shutter is opened for and t > 0 subject to the 'physics' that is taken to occur in this half space ?
Negative half-space Positive half-space What is the wave-field generated in the positive half-space when the shutter is opened instantaneously at t = 0, subject to the 'physics' (quantified by the partial differential equation which defines the physical model) that is taken to occur?

Focus and Context
The principal purpose of this paper is the study of a control problem associated with the scattering and the transient behaviour of a quantum field described by a fractional differential equation. This is an extension to the quantum shutter problem. However, in order to provide a more general and informative context, the paper includes a study of scattering theory for both an Electromagnetic (EM) and a Quantum Mechanical (QM) field. This is done in order to place the quantum shutter problem in a framework that attempts to show the connectivity between the scattering of EM and QM fields and their diffraction from: (i) a shutter that is permanently open and produces a diffraction pattern that is a function of space alone and independent of time; (ii) a shutter is opened at an instant in time and produces a diffraction-type pattern that is a function of both space and time.
By providing a review of the mathematical models associated with points (i) and (ii) above, the paper provides the reader with the background to the focus for this work and the original contributions that are developed. This is specifically concerned with an investigation of point (ii) for a fractional case, which is an example of quantum mechanics and control using Fractional Calculus. In this case, we are concerned with the control of a fractional QM field that is determined by the initial condition(s) that are applied.
In order to develop the themes presented in this paper in a way that is dimensionally compatible with the original QM shutter problem (where a time-dependent solution is developed for a one-dimensional space), the time-independent analysis is presented for the two-dimensional domain. In this way, for the time-independent problem, solutions are sought for r ∈ R 2 and presented graphically in the (x, y)-plane.
With regard to the time-dependent problem, solutions are evolved for r ∈ R 1 and presented graphically (as colour-coded maps) in the (x, t)-plane. There is no attempt to provide any specific quantitative measures of the solutions. Rather, the emphasis is to investigate the field patterns that occur for the different models that are investigated coupled with different example parameter values.
A complementary reason for taking this approach is that a two-dimensional timeindependent scattering theory is relatively less well documented and reported on than is the three-dimensional case. The material provided therefore fills a gap in the literature as well as providing an informative way of comparing solutions in the (x, y)-plane (timeindependent model) and the (x, t)-plane (time-dependent model).
fractional differential equations in relation to the quantum shutter problem, appear not to have been considered in previous works. This is undertaken from a theoretical point of view alone and includes an approximative approach to computing the Green's functions for the fractional differential equations considered coupled with the application of a memory function for defining time-fractional derivatives. In this respect, and, in a broader perspective, the study provided belongs to the field of analytical and numerical solutions to fractional differential equations [20].
In addition to the principal theme of the work, an exact solution to the Schrödinger scattering problem is proposed for investigating multiple scattering problems. This complements the use of the random walk models that are considered to address classical and fractional diffusion processes. Finally, a 'fractional Klein-Gordon equation' is proposed to investigate the shutter problem for a 'semi-relativistic' field.
The overall aim of the paper has been to consider a range of physical effects associated with the interaction of a QM field with a shutter in such a way that a unified theme is evolved in regard to studying the dynamics of standard and non-standard quantum fields. The advantages of the approach considered relate to the use of the Green's functions for solving the fractional partial differential equations considered. These Green functions are based on the application of conditions for studying a field close to the shutter over short time periods. The Green's functions constructed in this case are shown to be compatible with those for the conventional equations when the fractional differential reduces to the integer differential. The disadvantage of this approach is that the solutions are constrained by the approximations and conditions that are applied. However, the approach considered provides Green's function solutions that can relatively easily be investigated numerically (examples of which are presented in the paper) that transcend a numerical evaluation of the Mittag-Leffler function, for example.

Mathematical Preliminaries
For the purpose of completeness and internal referencing, the mathematical preliminaries associated with the material are provided. Some of the results are standard, while others are non-standard; in either case, some selected references are given as to where the results are provided.

The Fourier Transform
For a square integrable function f (r) ∈ L 2 (R n ) : C → C, and, taken over a Schwartz tempered distributional space, we define the Fourier and inverse Fourier transforms in the 'non-unitary form' as respectively, where r ≡| r | and k ≡| k |= 2π/λ is the spatial frequency for wavelength λ. These integral transforms define a Fourier transform pair, which, in this paper, is written using the notation For a function of time t denoted by f (t), the notation becomes and ω is the temporal (angular) frequency, which is related to a wave speed c 0 and the spatial frequency k by In this work c 0 is taken to be the speed of light.

The Convolution Integral
For r ∈ R n , the convolution of two functions f (r) and g(r) is defined as The dimension associated with this integral operator is taken to be inferred from the dimension of the functions to which the convolution operator is applied. In any dimension, the Convolution Theorem is taken to be applicable, i.e., For a function of time, the following notation is used where ⊗ t denotes the (non-causal) convolution integral in time.

The Dirac Delta Function
For r ∈ R n , we define the (n-dimensional) Dirac delta function as where f (r) = δ n (r) ⊗ r f (r)

Specific Fourier Transform Pairs
The following standard Fourier transforms are used in this paper and are available in standard Fourier transform tables and/or at [21]:

The Laplace Transformation
The Laplace and inverse Laplace transforms are defined as [22] Thus, the necessary and sufficient conditions for the existence of the Laplace transform are, respectively, that (i) the function f (t) should be piece-wise continuous in the given closed interval and must be of an exponential order; and (ii) the function f (t) exp(−st) should be absolutely integrable.
The Laplace transform pair is implied using the notation for which the following fundamental results are applicable: where ⊗ t denotes the causal convolution integral (in association with the Laplace transform alone and not the Fourier transform) The following specific Laplace transforms pairs are used in this paper and are available in standard Laplace transform tables and/or at [23]: A proof of the non-standard Relationship (12) is given in [24], for example.

Green's Function
For the classical diffusion operator D∇ 2 − ∂ t , where ∇ 2 is the Laplacian differential operator and D is the diffusivity, the Green's functions are given by [2] G(r, t) = H(t) 1 4πDt For the homogeneous Helmholtz operator ∇ 2 + k 2 , the time-independent and outgoing Green's functions, for r ∈ R n , are given by (functions of r ≡| r |) [2] g(r, k) =

p-and Uniform-Norms
We define the p-norm as with the uniform norm being defined as and standard properties: and f (r) ⊗ r g(r) p ≤ f (r) p g(r) p (18)

Field Notation
Throughout this work, the field, whose solution is sought, is denoted by u(r, t). This function may be the wave-field for a scalar electromagnetic field or a quantum field describing a photon or an electron, respectively, for example. In this case, the measurable quantity is taken to be given by | u(r, t) | 2 (the probability density). However, u(r, t) is also used to represent the density field which is the number of particles (photons or electrons, for example) per unit dimension. The difference in the physical meaning of u(r, t) is inferred through the context in which the physical models are considered, and the governing equations and the solutions thereof that are addressed.

Units
Some of the analysis provided in the latter part of this work involves the development of solutions to equations for which 'Natural units' (for particle and atomic physics) are applied. In this case, the speed of light c 0 , the Dirac constanth, and the electron mass m are all taken to be equal to 1, i.e., c 0 =h = m = 1 In the construction of a fractional derivative, a scaling parameter is required to maintain the dimensionality of the equation which can then be set to 1. Thus, for example, consider the first order spatial derivative d f (x)/dx of a dimensionless function f (x). The dimensions of this derivative are 1/L (the reciprocal of length L), but the fractional derivative d α f (x)/dx α , 0 < α < 1, say, has dimensions of 1/L α . In order to maintain the dimensionless nature of the function f (x) that is being fractionally differentiated, it is necessary to introduce an appropriate scaling factor η (a real constant with the dimensions of L) raised to an appropriate fractional power. Thus, in this particular example, the fractional derivative is written as η α−1 d α f (x)/dx α . However, for η = 1, the introduction of this notation is null and void and can be merely inferred in relation to preserving the dimensionality of the original (non-fractional) differential.

Fractional Partial Differential Equations and the Mittag-Leffler Function
Fractional partial differential equations and their solutions critically depend on the definition(s) that are applied for a fractional differential, and, more specifically, a fractional differo-integral operator. There are a wide range of conventional definitions for such operators, such as the Riemann-Liouville fractional integral and the Caputo fractional derivative [25]. Moreover, there are a broad range of generalised definitions, leading to the issue as to how many fractional derivative there are [26], and, whether there is an upper bound to the number of definition that can, in principle, be constructed, as discussed in [25]. Such definitions often depend upon whether they should incur a boundary and/or initial condition that is consistent with the physical model that is being considered. This is of particular importance in regard to time-fractional differential equations involving a fundamental solution using the Green's function through the application of a Fourier transform of a function over space and a Laplace transformation over time. In regard to the latter case, the Mittag-Leffler function [19] is often used to study the effect of initial conditions in regard to time-fractional control problems. The following sections provide an example of this in regard to solving the fractional diffusion equation.

Fourier-Laplace Transform Solution to the Diffusion Equation
Consider the one-dimensional classical diffusion equation (for unit Diffusivity) where we require a Green's function solution for u(x, t), x ∈ (−∞, ∞), t > 0, subject to the condition that u(x, t) → 0 as | x |→ ∞.
The standard approach to solving this problem is to apply the Fourier and Laplace transforms. Defining the transformed function as and using Relationships (4) and (11), the diffusion equation becomes (in Fourier-Laplace The solution is then given by computing the inverse Fourier and inverse Laplace transforms to this equation, which, in principle, can be done in any order. Taking the inverse Laplace transform first and using Relationship (13), and using Relationship (7) together with the convolution theorem, we obtain [2] u( This result is a convolution in space with the Green's function, i.e., where G(x, t) is given by Equation (14) for n = 1.

Fourier-Laplace Transform Solution to the Time-Fractional Diffusion Equation
Let us now repeat this analysis given in the previous section for the time-fractional diffusion equation subject to the same asymptotic conditions. In this case, following the same transformations, and using Relationship (12), we obtain To obtain an expression for U(k, t), the inverse Laplace transform of the function s α−1 /(s α + k 2 ) is now required; this is where the Mittag-Leffler plays an essential role.

The Mittag-Leffler and Wright Functions
The Mittag-Leffler function is defined as [19] E α,β (z) = ∞ ∑ n=0 z n Γ(αn + β) , α > 0, β > 0 and, importantly, has the Laplace transform relationship [24] E α,β (z) This result immediately allows the inverse Laplace transform of U(k, s), as given by Equation (19), to be written as so that the solution required becomes where and M α (z) is the Wright function [27], The geometry of convergence in regard to the Mittag-Leffler function is studied in [28].

Discussion
The Mittag-Leffler and Wright functions used to develop the Green's function solution to the time-fractional diffusion equation can be applied to solve other partial differential equations including the time-fractional Schrödinger equation [29]. Both these functions are non-standard or 'special functions' which involve definitions based on infinite series. Nevertheless, the numerical evaluation of such functions is available from different sources, e.g., [30,31]. This provides the facility for investigating Equation (20) numerically, based on the application of a discrete convolution operation. On the other hand, it arguable that the inverse Fourier and Laplace transforms are applied numerically to compute u(x, t) from Equation (19) directly, thereby by-passing the need to introduce such non-standard functions and their numerical evaluation. In this context, the approach taken in this paper follows an evaluation of the Green's functions based on various asymptotic's that are specifically designed to obtain Green's function of the type given by Equation (14). This is coupled with a definition of a fractional time derivative that is based on introducing an appropriate memory function m(t), where (for a causal convolution operator ⊗ t )

Review of Diffraction in Space for an Open Shutter-Diffraction by an Aperture
The purpose of this section is to provide the context to the 'Diffraction in Time' phenomenon. This is done by first re-visiting the basic mathematical foundations of diffraction theory, which provides a connectivity between diffraction in space and diffraction in timeessentially, time-independent and time-dependent diffraction, respectively. In this sense, the term 'Diffraction' can refer to a light beam or an electron beam or any other elementary 'particle' beam, at least in the context of the non-relativistic and spin-less case, when particles propagate at a speed much less than light speed. Although the models for diffraction in space are well known, it is revisited in this section so that the reader is made aware of both the differences and basic similarities between the spatial diffraction of an electromagnetic wave and a (non-relativistic) electron wave 'oscillating' at a fixed frequency.
With reference to Figure 1, let us consider the case when the shutter is permanently open, thereby producing an aperture. If the width of the incident wave-field in the negative half space is less than the width of the aperture, then one would expect the beam to travel through the aperture and continue propagating (from left to right) into the positive half space x ∈ [0, ∞). In this context, if the shutter remains closed, then the positive half space will remain free of any wave-field. However, if the beam width is larger than, or at least equal to, the width of the aperture, then part of the wave-field will interact with the edges of the aperture. It is this interaction that produces a diffraction pattern in the positive half-space. The characteristics of this pattern (measured over the coordinate y, say) will depend on the distance in x away from the aperture at x = 0, over which the wave-field is detected (as a function of y). In this case, (x, y) defines the coordinates in the positive (and negative) half-space. It is this issue that is now considered in terms of the basic mathematical formulation of the diffraction patterns that occur in the 'far-field' (the 'Fourier plane') and the 'intermediate-field' (the 'Fresnel zone'), which is undertaken for the two-dimensional case when r ∈ R 2 .

Electromagnetic Wave Equation
A light wave (or any other electromagnetic wave, representing a stream of photons) propagating in a free space at the speed of light c 0 , is governed by the homogeneous wave equation This a scalar model for the electric or magnetic field represented by the wave-field u(r, t). For r ∈ R 2 , and, for unit vectorsx andŷ, and u(r, t) is the Electric or Magnetic field in a two-dimensional space. This model is restrictive, as polarisation effects, due to the fact that electric and magnetic fields are vector fields, are not incorporated. Moreover, the model is two-dimensional. However, the purpose of this section is to contextualise the phenomenon of (spatial) diffraction in terms of the model illustrated in Figure 1. It is for this reason that a two-dimensional scalar wave-field model for diffraction is considered.

Time Harmonic Equation
If we consider the wave-field to be harmonic in time when we can write u(r, t) = U(r, ω) exp(iωt) where ω is a constant (angular) frequency, then we can consider the time-independent (homogenous Helmholtz) equation where λ is the wavelength. Note that this equation applies equally well for the case when there is a frequency spectrum, given that, in this case, where u(r, t) and U(r, ω) are Fourier transform pairs, i.e., For a wave-field U(r, ω) with a characteristic spectrum P(ω) say, and using the convolution theorem, where ⊗ t denotes the (non-causal) convolution integral in time. Thus, for a spectrum consisting of a single frequency ω 0 , if we let P(ω) = 2πδ(ω − ω 0 ), then thereby recovering the time harmonic case for constant frequency ω 0 . In this case, the wavefield is taken to oscillate at the same frequency for all time and only the spatial characteristics of the wave-field U(r, ω) are considered (for constant ω). No initial conditions for the field u(r, t) and/or its gradient, when we need to specify the values of u(r, t) and/or ∂u(r, t)/∂t at t = 0, respectively, are relevant. The model therefore relates to a non-causal system and there are no initial conditions in time that control any time evolution of the 'system'.

Inhomogeneous Helmholtz Equation
With reference to Figure 1, we consider the case when the incident beam interacts with the edges of the aperture. While the frequency is constant, in regard to Equation (23), this interaction can be taken to slightly reduces the propagation speed, which is strictly less than light speed. Thus, in a general context, the light speed is inhomogeneous in space over the region of the aperture. We can model this inhomogeneity by replacing 1/c 2 0 in Equation (22) so that, if V(r) = 0, ∀r, then the space is homogenous. In this case, Equation (23) and it is this equation that we now need to solve for U(r, k) given V(r), which is taken to be of compact support. This result is of course based on considering the conversion of a constant wave-speed to a variable wave-speed and has been achieved in a phenomenological context. However, by decoupling the macroscopic Maxwell's equation it can be shown that, if the aperture is a non-conductive dielectric, then V(r) = r (r) − 1, where r ≥ 1 is the relative permittivity (ignoring polarisation effects) [3]. Similarly, if the aperture is a conductive dielectric, then V(r) = r − 1 − iz 0 σ(r)/k where σ is the conductivity and z 0 is impedance of free space [3]. In either case, the relative permeability is taken to be that of free space and U(r, k) is taken to be the electric field.

Non-Relativistic Quantum Wave Equation
For the non-relativistic case, an electron wave u(r, t), say, propagating in an n-dimensional free space, is described by the Schrödinger equation whereh is the Dirac constant and m is the electron mass. Here, the operator ih∂ t is associated with the energy E =hω of the electron wave-field u(r, t).
If an electron beam is incident upon an open aperture, the edges of the aperture represent a potential energy V(r) of space r, and, given that the total energy is given by E − V, then from Equation (25), we obtain the inhomogeneous equation Thus, for harmonic time variations in the wave-field u(r, t) at a frequency ω, we can again consider the case where u(r, t) = U(r, ω) exp(−iωt) giving (the time-independent inhomogeneous Schrödinger equation) where k 2 = 2mE/h 2 and V(r) := 2mV(r)/h 2 .
Comparing Equation (24) with Equation (26), the basic differences lie in the definitions of k and V(r) and whether or not −k 2 is a coefficient of V(r), respectively. Thus, for constant k, any solution acquired for U(r, k) is the same in either case. In both cases, the time-independent equation describes a non-causal system with no time controlling properties and associated initial conditions. The fundamental solutions to these equations are now considered.

Fundamental Solution: The Green's Function Solution
Since, for constant k, Equations (24) and (26) are effectively the same, we consider the Green's function solution to Equation (26). This is given by the Lippmann-Schwinger equation [2] U where U s (r, k) is the field scattered by the 'scattering function' V(r) given by and ⊗ r denotes the two-dimensional convolution integral for r ∈ R 2 . The function g(r, k) (where r =| r |) is the two-dimensional 'outgoing' free space Green's function, which is the solution to for two-dimensional Dirac delta function δ 2 (r). The 'incident wave-field' U i (r, k) is taken to be the solution to the homogeneous Helmholtz equation Equation (27) is a fundamental solution to Equation (26) since given that For r ∈ R 2 , the Green's function is given by [2] g(r, k) = i 4π 0 is the Hankel function of the first kind and of order zero given by Thus, the asymptotic form for the two-dimensional Green's function is Application of this conditional form is valid when kr >> 1 so that the wavelength of the wave-field is small compared to the distance between the point source of the field δ 2 (r) and the position r at which the field is observed.

Solutions under the Born Approximation
Equation (27) is an integral equation of Fredholm type and has a range of solution method for evaluating U(r, k). The simplest and most common (approximate) solution to this equation is obtained by applying what is commonly referred to as the Born approximation [2]. This is where it is assumed that the convolution integral can be approximated by g(r, k) ⊗ r V(r)U i (r, k). Application of this approximation requires that (for some norm denoted by • ) where the Born scattered field is given by The Born scattering condition given above implies that the scattering is a 'weak effect', i.e., the scattered field is a small perturbation of the incident field. In physical terms, this means that there are no multiple scattering effects taken to be present. Thus, the Born scattered field is a model for single scattering events alone. A quantification of this condition will be considered later.
Multiple scattering events can be taken into account through iteration of Equation (27) (the first iteration being the Born scattered field), which requires that the series converges [2]. This is a formal solution to the scattering problem, a problem that has been studied over many years. In this context, Appendix A provides a complementary solution to Equation (27) for r ∈ R n , coupled with a short study on the type of functions to which V(r) must conform, for the solution(s) to be valid.
Given that U i (r, k) = exp(ik · r) is a solution to Equation (29), the Born scattered field can be written as (for r ∈ R 2 ) It is assumed that the incident field U i (r, k) is a unit plane wave, but the solution to Equation (29) applies to any plane wave with arbitrary amplitude, an amplitude that may also be a function of k.
The evaluation of the scattered field U s (r, k) now depends on an analysis of the function exp(ikr)/ √ r as a kernel of the convolution integral, i.e., Equation (28). This analysis is based on noting that (where s =| s |) The difference between evaluating the scattered field in the Fourier or Fresnel zones is dependent on the incorporation of just the second or the second and third terms of the binomial series given above, as shall now be demonstrated.

Scattering Amplitude in the Fourier Zone
If we assume that s/r << 1, then the term s 2 /2r 2 and higher order terms in the binomial series can be ignored. We can then consider the result In this case, the scattered field is given by is the scattering amplitude. It is immediately clear that the scattering amplitude is, in effect, given by the two-dimensional Fourier transform of V(r).

Scattering Amplitude in the Fresnel Zone
If we assume that the term r 2 /2s 2 is much less than 1, then only higher order terms after this term can be ignored. In this case, we consider a representation of the complex exponential given by This result can now be written in the compound form which allows the scattered field to be written as and α = k/2r is a constant. Note that the term k/2r is not a variable in the convolution integral but a constant coefficient of the variable | r − s |.

Diffraction Patterns
In the context of Figure 1, let us now consider how to construct a model for the diffraction pattern in the Fourier and Fresnel zones.

Diffraction Pattern in the Fourier Zone
The aim is to consider Equation (35) and modify the result in a way that is compatible with the geometry associated with Figure 1, where a one-dimensional screen (upon which the diffraction pattern is taken to be observed as a function of y) is placed a long distance away (in the 'far-field') from the aperture so that x/y << 1. The incident field is taken to be a function of x alone, as it propagates in the x-direction (from left to right, in the negative half-space as illustrated in Figure 1). Thus, given thatm =x, then exp(ikm · r) = exp(ikx). Furthermore, since we are considering the result in the far-field of the positive (twodimensional) half-space (x, y), x > 0 then With these results, Equation (35) becomes (using the notation s =xs x +ŷs y , for Cartesian coordinates), This result now needs to be further quantified with regard to the function V(s x , s y ). In terms of an open aperture relating to Figure 1, for a width W and thickness T, we can consider that where the aperture is taken to be a unit constant for which values other than 1 may of course be applied. In this context, the 'otherwise' condition is simply stating that, apart from the aperture, no radiation can penetrate into the positive half-space as illustrated in Figure 1. Thus, from Equation (37), the scattering amplitude becomes It now becomes clear that the diffraction pattern (as a function of y) for a fixed value of x is determined by the one-dimensional Fourier transform of a constant where s y ∈ [−W/2, W/2], the measurable intensity pattern being given by | A(K) | 2 .
There is another approach that can be used to derive what is essentially the same result and involves a model where the aperture is taken to be infinitely thin. However, in order to eliminate the trivial result A(K) = 0 if T = 0, we resort to the application of a delta function to model an aperture that is infinitely thin (in a hypothetical context). This is accomplished by considering an expression for V(s x , s y ) given by given that (by definition) For this model, the constant T is eliminated, and the scattering amplitude intensity is given by

Diffraction Pattern in the Fresnel Zone
The same basic approach can be used to generate an expression for the diffraction pattern in the Fresnel zone given Equation (36). Specifically, for an infinitely thin aperture and withm =x, then, in Cartesian coordinates, where α = k/2x given that y/x << 1.
The convolution integral is over y alone and the measurable intensity of the Fresnel diffracted field becomes Comparing Equation (39) with Equation (38), for an infinitely thin aperture, we observe that the difference is compounded in the one-dimensional convolution of the (Fresnel) kernel exp(iαy 2 ) with a unit constant of compact support y ∈ [−W/2, W/2] and the Fourier transform of the same constant (of compact support). Figure 2 shows example diffraction patterns based on Equations (38) and (39) computed using a uniformly sampled grid of 1000 elements for y ∈ [−1, 1] with W = 0.05 and k/x = 500. The results have been obtained using the Matlab functions fft (with fftshift) and conv, to compute the Fourier transform and the convolution operations, respectively. In order to enhance the dynamics of the intensity patterns I =| A | 2 , the plots are first normalised and plotted on a logarithmic scale given by log(0.01 + I/ I ∞ ).  (38) and (39) in the (twodimensional) Fresnel and the Fourier zones generated by the scattering of a plane wave propagating through and interacting with a single (infinitely thin) aperture as illustrated-not to scale.

Fourier and Fresnel Zone Equivalence
Equations (38) and (39) are solutions for the scattering cross-section (the intensity of the scattered wave-field) based on Equation (34) under the conditions that s/r << 1 and s 2 /r 2 << 1, respectively. Is there a simple relationship between these two results leading to an equivalence principle? The answer to this question involves using the same equivalence relationship for a Fourier transform that was originally designed to develop a fast Fourier transform algorithm by Leo Bluestein in 1970 [32]. It is based on re-writing the discrete Fourier transform in terms of a discrete convolution sum. However, this approach can also be applied for the case of continuous functions as shall now be addressed.
The key to developing the equivalence principle is to note that a Fourier transform (for r ∈ R 1 ) can be written in the form Thus, using the notation for a convolution integral, we can write Equation (38) in the form Comparing this result with (39), where α = k/2x; then, for the case k/x = 1, we can construct two equations, namely, for the scattering cross-section in the Fresnel zone and the Fourier zone, respectively.
In the context of these two equations, and, specifically for the case when k/x = 1, the difference between the diffraction pattern in the Fresnel zone and the Fourier zone becomes the difference between convolving the chirp function exp(iy 2 /2) with the aperture function V(y) and convolving the same chirp function with the conjugate chirped aperture function V(y) exp(−iy 2 /2).

Quantification of the Born Approximation
Equations (38) and (39) are effectively the result of considering different geometric scenarios, both of them being taken under the Born approximation-Condition (33). It is therefore worth quantifying this approximation within the context of the model considered for the scattering function, namely, an infinitely thin slit.
For this purpose, we consider the condition using a Euclidean norm for r ∈ R 2 , i.e., an integrable function f (r) for which we can define Using inequalities (16), (17) and (18), we can write and Condition (33) therefore becomes Given the expression for g(r, k) provided by Equation (32), and, using polar coordi- Thus, defining the root mean square of V(r) as being the condition for the Born approximation to be valid can be written as This is the condition for the scattering of a non-relativistic electron or ion beam, quantified by Equation (26). For the case of a beam of photons quantified by Equation (24), the difference is that V(r) is replaced by −k 2 V(r). The condition for the Born approximation to be valid is then given by These conditions have been derived for a disc with radius R, which is a measure of the square root of an area. For an infinitely thin slit of the type illustrated in Figure 2, the area is negligible, and, in this context, both the conditions (40) and (41) are valid. Note, however, that, in the latter case, and, in general, we require that the wavelength λ is significantly larger than the physical dimensions of the scatterer, i.e., for a given value of V, we require that λ >> R. In other words, for an infinitely thin aperture (or otherwise), the width of the aperture must be small compared to the wavelength in order to generate a diffraction pattern. In the case of Condition (40), λ << R −3 .

The Electromagnetic Shutter Problem
Having reviewed time-independent scattering for r ∈ R 2 under the Born approximation for an infinitely thin aperture, and obtained expressions for the intensity of the diffraction pattern in the plane, we now consider the time-dependent case for r ∈ R 1 and for a scalar EM wave. In this context, the (x, y) plane considered previously is replaced with an (x, t) map. The purpose of this is to emphasise the fact that, while the characteristics of the time-independent diffraction of EM and QM fields are essentially the same, the time-dependent shutter problem is markedly different for the two fields.
With reference to Figure 1, we consider a unit amplitude plane wave that is incident on a shutter, which is opened instantaneously at t = 0. The width of the wave-field is taken to be less than the width of the aperture so there are no diffraction effects. In this case, for c 0 = 1, Equation (22) becomes , with initial conditions given by To solve this problem, we consider the time-dependent Green's function (for c 0 = 1). In this case, the Green's function is given by the solution to where the notation G(x | y) denotes G(| x − y |). Time dependent problems of this type require application of the Laplace transform. This is because the (one-sided) Laplace transform is a causal transform and therefore provides the facility to implement initial values for the field u(x, t) and/or its derivative at t = 0. Using Relationship (11), Equations (42) and (43) transform to where u(x, s) and where respectively. Pre-Multiplying Equation (44) by G and Equation (45) by u, subtracting the results and then integrating over x and t, we can construct a solution for u(x, s) given by The values of u(x, t) and ∂ x u(x, t) at x = ±∞ are un-specified. Moreover, given that they can be taken to be beyond the spatial domain of interest, they are irrelevant, and can therefore be set to zero. This is consistent with the asymptotic condition that u(x, t) → 0 as | x |→ ∞. Thus, given the initial condition specified, we can consider the solution u(x, s) = u(x, 0) ⊗ x sG(| x |, s) Inverse Laplace transforming, and, noting that we then obtain the time-dependent solution given that G(| x |, 0) = 0. The time-dependent Green's function-the solution to Equation (43)-is given by (for This result can be obtained, by taking the inverse Fourier transform of the timeindependent (out-going) Green's function, which, for r ∈ R 1 , is given by Thus, using Relationship (2), and then applying the convolution theorem, Given that, for x > 0, then, from Equation (46), we have which describes a wave travelling from left to right (using the quantum mechanical convention [33] in the positive half space of Figure 1). Figure 3a shows a space-time map of the real component of u(x, t), x ∈ [0, 1], t ∈ [0, 1] for ω = 20π computed for a regular grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colour map after normalisation-a map of Re[u(x, t)] + 1/2 ∈ [0, 1]. Note that the intensity of the field | u(x, t) | 2 in this case is a constant ∀x > 0, t > 0.
The solution given by Equation (46) is dependent on the initial condition u 0 (x) = 0. For the case when this condition is non-zero, the equivalent solution is

Multiple Scattering and Diffusion
With reference to Figure 1, the analysis presented so far has assumed that the positive half space is homogenous. The scattering effects are due to an open shutter which yields diffraction patters as presented in Sections 5.2.1 and 5.2.2 for the time-independent case. This is based on an application of the Born approximation for an aperture that is infinitely thin (when the Born approximation is valid). The shutter problem discussed in Section 6 is illustrative of the fact that, if a beam of photons is incident on a shutter that is instantaneously opened at t = 0, then they will simply propagate through the aperture without diffracting, given that the width of the beam is taken to be less than the width of the aperture. This is not the case, when a beam of electrons (or other non-relativistic particles) are incident upon a shutter as will be discussed later. However, before discussing the quantum shutter problem, it is informative to consider the effect of multiple scattering in terms of the physics of diffusion and why this comes about. This is because there is a natural connectivity between this effect and the quantum shutter problem, which relates to the form of the Green's function for the diffusion equation. Furthermore, using this approach provides a way of defining a time-fractional derivative via a memory function which is introduced through a derivation of the Generalised Kolmogorov-Feller Equation as addressed in this section.

Multiple Scattering
Consider the case when, with reference to Figure 1, the shutter is opened permanently and a beam of photons passes through the aperture without diffracting from the edges as discussed in Section 6. However, the positive half-space is now taken to be composed of a random distribution of scatterers rather than a homogenous free space. In this case, the incident field will interact with a multitude of scatterers, each of them taken to have a physical size of the order of a wavelength. This multiple scattering phenomenon could describe the interaction of photons with water droplets in a fog or a cloud or a turbulent medium, for example. It could also represent the multiple interactions of free electrons in a plasma, for example.
For multiple scattering interactions, the Born approximation is not valid. This is because the Born approximation assumes that only single scattering events occur and that there are no higher order interactions (double, triple, . . . , scattering). Formally, these higher order interactions, relate, respectively, to the terms associated with the iterative solution of Equation (27)-the 'Born series'-when the scattered is given by: The first term on the RHS represents single scattering events (the Born approximation) and the second term represents double scattering events and so on. There are other approaches to modelling a multiple scattered field includes a solution for evaluating the exact scattered field as given in Appendix A. However, if V(r) describes a random ensemble of scatterers, then one can expect the field that evolves to become akin to a process of diffusion where each photon or electron, for example, is taken to undergo a random walk in the plane for r ∈ R 2 . In addition, this principle is the key to developing a diffusive model for photons (light diffusion) or electron and ion diffusion (which underpins transport theory in a plasma, for example). Modelling multiple scattering as a random walk process requires the problem to be quantified in terms of a time-dependent phenomenon, and, in the following section, the basic field equation-the evolution equation-associated with such a process is addressed.

The Evolution Equation for a Random Walk
For a canonical ensemble of particles undergoing a random walk involving elastic scattering intersections, the evolution of the density field u(x, t) for r ∈ R 1 is given by [34] u where u(x, t) is the number of particles per unit length (a real dependent variable where u(x, t) ≥ 0, ∀x, t > 0) and p(x) is the Probability Density Function (PDF), where Thus, p(x) determines how the particles undergoing a random walk are distributed over x in a short period of time t, thereby providing the initial density function u(x, t). As the particles continue to undergo random walks, the density function will continue to evolve in time, subject to the same statistical properties characterised by the PDF. At any future time t + τ, the density function will be determined by the density function at time t convolved with the PDF. This argument is the basis for Equation (50) which is applicable for any PDF. The type of PDF that is considered defines the statistical characteristics of the random walk that occurs, which in turn determines the evolution of the density field.
In terms of the scattering of light, for example, u(x, t) represents the number of photons per unit length (the photon density) at a time t. The density function is equivalent to | u(x, t) | 2 in Equation (49) for r ∈ R 1 where u(x, t) F ←→ U(x, k), V(x) is replaced by −k 2 V(x) (for the EM case), and where V(x) is a random variable. The incident field at t = 0 is not the wave field u(x, 0) but the intensity | u(x, 0) | 2 = 1 (for a unit amplitude wave). Thus, with reference to Equation (50), the space-time field u(x, t) must be interpreted in terms of the intensity of a multiple scattered wave field-the scattering cross section for strong scattering interactions leading to an evolution in the spatial distribution of photons.
For quantum scattering (the scattering of electrons in a plasma, for example), | u(x, t) | 2 is the probability of a particle existing at a point in space and time, and, in this sense, also represents the electron density, given in Equation (50).

The Evolutionary Equation for τ << 1
Applying a Taylor series in time to the function u(x, t + τ) in Equation (50), then Thus, for τ << 1, the evolution equation has the form This is known as the Kolmogorov-Feller Equation [35].

The Evolution Equation for a Memory Function
For those cases when we can not assume that τ << 1, it is conventional to introduce a 'memory function' m(t) by expressing the Taylor series for the function u(x, t + τ) in Equation (50) as We can then write the evolution equation in the form This equation is known as the Generalised Kolmogorov-Feller Equation [36]. The reason for using the term 'memory function' is that, if m(t) = δ(t), then Equation (52) is recovered, and, since the delta function is localised entirely at t = 0, there is no time memory association in regard to the term ∂ t u(x, t). However, if m(t) = δ(t), then the convolution of m(t) with ∂ t u(x, t) implies that the value of this convolution integral at a time T, say, will depend on the time dependence of m(t), for t < T, i.e., the 'history' of the function over t ∈ [0, T]. This idea provides an approach to defining a fractional time differential in terms of the convolution of an integer time differential with an appropriate memory function as will be introduced later in Section 9.3.

The Diffusion Equation
The process of (classical) diffusion can be taken to be the case when p(x) is a normal distribution which, for standard deviation σ, and variance σ 2 , is given by where P(k) is the characteristic function-the Fourier transform of p(x).
In Fourier space, where u(x, t) , t), and, with application of the convolution theorem, we can write Equation (53) as which yields the equation, using Relationship (4), where D is the 'Diffusivity' given by In this context, the classical diffusion equation is recovered when m(t) = δ(t) In terms of Equation (53), the diffusion equation is a result of considering the case when where it is noted that thereby showing that p(x) conforms to Equation (51). Although we have applied a condition on the variance, it is nevertheless informative to explore the relationship between this measure and the physical nature of a random walk. Within the context of the condition, σ 2 << 1 is used to derive the diffusion equation, and a larger value of σ 2 implies that a particle undergoing a random walk (conforming to a Gaussian distribution) can travel further from the source over a specific period of time. In other words, for a constant value τ, the larger the diffusivity, the greater the range of the random walks so that particles travel further. This implies that the spatial density of scatterers is less so that photons or electrons can diffuse further over the same time interval-the mean free path is greater. This property is compounded in the characteristics of the Green's function for Equation (54), i.e., where the width of the spatial variance of the exponential function at a given time t is determined by the value of D. This is the same for r ∈ R n where the Green's functions have the same basic form, as given by Equation (14).

Green's Function Solution to the Diffusion Equation
We apply exactly the same analysis as that given in Section 6 but for Equation (54) instead of Equation (42). This is subject to the initial condition that u 0 (x) ≡ u(x, 0) = 1, noting that this initial condition is for the intensity of the incident wave field and not the wave field itself, as illustrated in Figure 1. In this case, we obtain the solution where G(x, t) is given by Equation (56).
This result is obtained by using relationship (11), when, for D = 1, Equation (54) and the associated equation for the Green's function, i.e., and respectively. Thus, pre-multiplying these equations by G and u, respectively, subtracting the results and integrating over x ∈ (−∞, ∞), we can construct a solution for u given by (for u(x, t) = 0 and ∂u(x, t)/∂x = 0 at x = ±∞) where G(x, 0) = 0. Inverse Laplace transforming then yields Equation (57) Figure 3b shows a space-time map of u(x, t) for x ∈ [0, 1], t ∈ [0 + , 1] (with D = 1), computed for a regular grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colourmap after normalisation-a map of u(x, t)/ u(x, t) ∞ .
The diffusion of photons or electrons, for example, based on a random walk model in which the 'system' is characterised by a Gaussian PDF, provides a useful background to the (non-relativistic) quantum shutter problem as presented in Section 9. This is essentially because the model used in the analysis of the quantum shutter problem can be said to be based on the diffusion equation for imaginary time, i.e., forh = 1 and m = 1, then with D = 1 2 and replacing t with it, Equations (54) and (25) become one and the same (for r ∈ R 1 ). However, before considering this, we formulate and then consider solutions to the space-fractional diffusion equation. This is presented on the following section.

Fractional Diffusion
Consider the case when we generalise the Characteristic Function to the form where a has the dimensions of [Length] γ and the index γ is the 'Lévy index'. In this case, the distribution can only be evaluated as an asymptote (at least for 1 < γ < 2) given that In the context of the analysis presented in Section 7.3, and, using Relationship (5), the introduction of a Lévy distribution yields the fractional diffusion equation where D γ = a τ which has dimensions of [Length] γ [Time] −1 . The field u γ (x, t) is used to denote that it is the solution to a fractional differential equation with a spatial differential exponent determined by the value of γ.

Green's Function for the Space Fractional Diffusion Equation
We now consider the Green's function for the space fractional diffusion operator for D γ = 1, which is given by the solution of The basis for this is to first of all consider the inter-relationship between the timeindependent and time-dependent Green's function for the wave operator and diffusion operator. This is reviewed in the following sections.

Green's Function for the Wave Operator
We require the solution of the time-independent form of the problem is given by a solution to the equation Expressing g(| x |, ω) in terms of its spatial Fourier transform g(k, ω), i.e., This expression arises from an evaluation of an integral which has two poles at ±ω, the pole for k = ω providing an expression for the out-going Green's function as given.

Green's Function for the Classical Diffusion Operator
In the context of evaluating the Green's function for the wave operator discussed in the previous section, suppose we now consider Equation (62) for γ = 2. In this case, the time-independent Green's function is given by the solution of Applying exactly the same procedure as presented in the previous section, the timeindependent Green's function becomes This result is based on an evaluation of the integral in the equivalent equation for the Green's function, i.e., Hence, inverse Fourier transformation, coupled with Relationship (10) yields thereby recovering the Green's function for the diffusion equation given by Equation (14) for r ∈ R 1 and D = 1.

Green's Functions for the Space Fractional Diffusion Operator
Given the analysis provided in Section 8.3, the time-independent Green's function associated with Equation (62) can be taken to be given by (with s = iω) which is a generalisation of the result for γ = 2 to the case when 1 < γ < 2. To implement the inverse Laplace transform required to obtain G γ (| x |, t), we consider a solution under the condition that | x | 2 << 1 in Laplace space. This condition is analogous to the condition used for developing an expression for the scattering amplitude in the Fresnel zone, as discussed in Section 5.1.2. Thus, we consider a truncated series expression for g γ (| x |, ω), which is based on the series representation of the exponential function and given by Using Relationships (8) and (9), we can then write where it is noted that thereby preserving the negative exponentiation associated with a diffusive process over the range of values for γ considered. The essential 'trick' in developing this expression for G γ (| x |, t) is to note that the delta function δ(t) = 0, ∀t > 0, which eliminates the second term, thereby (after factorising) allowing the remaining terms to be recognised as the first two terms in the series definition of an exponential function. This requires that Furthermore, noting that Γ(1/2) = √ π, Γ(−1/2) = −2 √ π, Γ(1) = 1 and 1/Γ(−1) = 0, we can write, for t > 0, thereby recovering the expression for the Green's function when γ = 2, i.e., Equation (14) for r ∈ R 1 . The result for γ = 1 reflects the Green's function associated with the operator ∂ x − ∂ t , which is given by 1 2 H(x − t). For completeness, Appendix B considers a similar approach for generating the asymptotic fractional diffusive Green's function for r ∈ R 3 ; a result that is used later in Section 10.
Repeating the analysis for Equation (60), we obtain

Green's Function Solution to the Space Fractional Diffusion Equation
With the Green's functions for the fractional diffusion operator evaluated for r ∈ R 1 , and, subject to the initial condition u 0 (x), the fractional diffused field, defined by Equation (60), is given by In reference to Figure 1, the initial condition u 0 (x) = 1, ∀x ≤ 0 represents a constant stream of photons or electrons with a unit density. They then undergo multiple scattering in the positive half space where the multiple scattering is modelled in terms of a space fractional diffusive process; a random walk characterised by a Lévy distribution with a Lévy index between 1 and 2. The field u(x, t) describes the density of photons (the intensity of light) or (non-relativistic) electrons (the number of electrons per unit length) that occurs close to the shutter thereby conforming with the condition that | x | 2 << 1 to evaluate the Green's function that is used.
An example of the visualisation of the properties associated with Equation (65) is given in Figure 4a. The figure shows a space-time map of u γ (x, t) given by Equation (65) for γ = 1.9, where x ∈ [0, 1], t ∈ [0 + , 1] and D γ = 1. The map is computed for a regular grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colour map after normalisation, i.e., a map of u γ (x, t) ∞ , which is visually enhanced using the Matlab function histeq. Figure 4b shows a map with exactly the same parameters but for γ = 1.1. In the former case, and, as expected, the map shows a field that approaches classical diffusion where the density decreases in both space and time. However, for γ → 1, the decay of the field in space becomes increasingly uniform in time. We may consider a demarcation between the values of γ and the type of diffusion that occurs. For γ = 2, the field is classically diffusive, for γ ∈ [1, 2), the process is sub-diffusive and, for γ < 1, the process becomes hyper-and super-diffusive.
The methods considered can of course be extended to two-and three-dimensional models. In this respect, the application of space fractional diffusion includes the prediction of market behaviour under the Fractal Market Hypothesis (e.g., [34,37,38]) and the processing and analysis of images [39][40][41], for example.

The Quantum Shutter Problem
As discussed in the Introduction, the 'quantum shutter' problem leads to the principle of 'diffraction in time' which is a fundamental transient phenomenon in quantum mechanics. In this case, we consider a pencil-line beam of non-relativistic particles described by wave-function u(x, t), which satisfies the time-dependent and one-dimensional Schrödinger's equation. For r ∈ R 1 , and, using natural units m =h = 1, we are interested in a solution to the equation subject to an initial condition u 0 (x) ≡ u(x, t = 0). With references to Figure 1, the initial condition is taken to describe a unit plane wave that is incident upon the shutter before it is instantaneously opened at t = 0 and the particle beam is allowed to pass through. The beam is described by a right-traveling unit amplitude plane wave exp(ikx) which is incident upon a closed shutter placed at x = 0.
The shutter is initially taken to be a perfect absorber so that, in the positive half-space, u(x, t) = 0, x > 0. It is then opened instantaneously at t = 0 after which the particle beam is free to travel into the positive half-space. The problem is to find the transient behaviour of the particle beam once it has been made 'free' to travel in the positive half-space after the shutter has been opened subject to the initial condition If we consider this problem in regard to the propagation of photons, then the wave function u(x, t), is governed by the classical wave equation as discussed in Section 6. Intuitively, one would considers the photons to propagate into the positive half-space after the shutter is opened so that, for t > 0, this half-space is characterised by a linear wave traveling from left to right as verified in the analysis presented in Section 6. However, for a beam on (non-relativistic) electrons, the physical behaviour is entirely different. This is a consequence of the nature of the Green's function, which, given Equation (14), for D = 1/2, r ∈ R 1 , and imaginary time it becomes which is the Green's function for Equation (66), Following the approach to generating a Green's function solution given earlier in the paper, the solution for u(x, t) is given by illustrating that, like the classical diffusion equation, u(x, t) → 0 as t → ∞.

Solution for the Schrödinger Equation
Compared to a beam of photons, the transient behaviour associated with a beam of electrons, for example (subject to the instantaneous opening of a shutter), is determined by a chirp function. This is a direct consequence of Equation (66) being characterised by the (imaginary) time derivative operator −i∂ t compared to Equation (22)-for r ∈ R 1 -which is characterised by a second order (real) time derivative operator ∂ 2 t . This is also the case for the multi-dimensional Schrödinger equation (forh = m = 1) given that, from Equation (14), the Green's functions for this case are given by The solution given by Equation (68) can be studied analytically using a change of variables to obtain an expression for the probability density given by [1]  However, Equation (68) can also be computed numerically via application of a convolution sum and the solutions displayed on a space-time map of the type given in Figure 3. Two examples of this are shown in Figure 5, which provides time-diffraction maps of the probability function | u(x, t) | 2 , where u(x, t) is given by Equation (68). In both cases, the function | u(x, t) | 2 is computed on a uniformly sampled grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colour map after normalisation, i.e., a map of | u(x, t) | 2 ∈ [0, 1]. The maps have been enhanced using the Matlab histogram equalisation function histeq. They show examples of the oscillations in time (and space) that decay over time and decrease in frequency, the effect being dependent on the value of the wavenumber k associated with the initial condition compounded in Equation (67). As the value of k increases, the oscillatory behaviour 'expands' over the (x, t) plane. The examples given in Figure 5 illustrates the oscillatory behaviour and decay of the probability density function for non-relativistic electrons. This is of course very different to the optical case (a beam of photons characterised by the classical wave equation as presented in Section 6) when, by comparison, the intensity function in the positive halfspace is a constant, i.e., Figure 3a would be a map of the value 1. There is a similarity between the expression for u(x, t) given by Equation (68) to the Fresnel zone solution given in Equation (39). This is why the transient phenomenon associated with the quantum shutter problem is dubbed 'diffraction in time' [1]. It is a quantum effect that requires the application of very low velocity particles, observed over small distances, so that the time interval of an experimental observation is very small, when the Fresnel-type oscillations of the density function in time are prevalent as illustrated in Figure 5.

Solution for the Space Fractional Schrödinger Equation
In Section 8, the fractional diffusion equation was constructed from Equation (53) on the basis of a delta memory function and a PDF given by δ(x) + aδ (γ) . Suppose we apply the same (phenomenological) principle for constructing the space fractional Schrödinger equation, but for the imaginary memory function −iδ(t). In this case, because the delta function is imaginary, the field u(x, t) can no longer be taken to be a density function but a complex wave-type function (on a phenomenological basis). In this context, with a = η γh2 /(2m), where η has dimensions of [Length] γ , and τ =h, we recover the fractional Schrödinger equation (forh = m = η = 1) which is an identical equation to Equation (60) for imaginary time it and with D γ = 1/2. Thus, we can use the solution given by Equation (65) for the initial condition given by Equation (67), i.e., but where the Green's function (for | x | 2 << 1) is now given by In this case, it is noted that the Green's function includes an exponential decay which is eliminated as γ → 2 for all values of x and t, conforming to Condition (63), and/or as Two examples of this solution are shown in Figure 6, which provides maps of the probability function | u γ (x, t) | 2 , where u γ (x, t) is given by Equation (70). The function | u γ (x, t) | 2 is computed on a uniformly sample grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colormap after normalisation-a map of | u γ (x, t) | 2 ∈ [0, 1]. The maps have been enhanced using the Matlab histogram equalisation function histeq.
These results are illustrative of the fact that, close to the shutter, and, as time increases (after the shutter has been opened), then, for γ → 2, and, for k = nπ, the probability density forms n fringes which decay over time. However, as γ → 1, the decay over time is eliminated, and, consequently, the density function forms a standing fringe pattern that is constant in time. If such a density function could be generated experimentally, then it could form a periodic grating for the diffraction of an electron beam, for example that is perpendicular to the plane of the grating, for which the diffraction models given in Section 5.2 become relevant.
is given by Equation (70) for k = 4π and γ = 1.9; (b) pace-time map of the function is given by Equation (70) for k = 4π and γ = 1.1.

Time Fractional Solutions
Having discussed the space fractional quantum shutter problem in the previous section, we now consider the time fractional case. This is based on using a memory to represent a fractional time derivative given that ∂ α t = ∂ α−1 t ∂ t , 0 < α < 1. For this range of α, the operator ∂ α−1 t is an anti-fractional derivative-a fractional integrator. The most basic fractional integrator is the Riemann-Liouville integral which can be defined as In the context of this definition, and, with respect to Equation (53) with p(x) = δ(x) + aδ 2 (x), we consider a time fractional diffusion equation given by and it is noted that However, for α = 1, when the classical diffusion equation should be recovered, 1/Γ(0) = 0 and thus a result that is incompatible with the expected limit α → 1. This incompatibility is due to the fact that the transform 1 s α L ←→ 1 Γ(α)t 1−α is strictly only applicable for α < 1. Thus, in order to construct a memory function that is compatible with the expected result when α = 1, we can define it as noting that

Time Fractional Diffusion
Following the approach detailed in Section 7.4, we can now construct a solution for density field u(x, t) given by (for a/τ = 1) Compared to the solution for the space-fractional diffusion equation, the time-fractional case includes a causal convolution in time with the memory function. Furthermore, comparing the fractional indices, it is clear that 1/γ = α/2.

Time Fractional Quantum Shutter
This same principle applies to developing a solution for the time fractional shutter problem. This essentially involves replacing t for it when the solution given by Equation (71)-for the wave function u(x, t)-is given by Two examples of the solution given by Equation (72) are shown in Figure 7, which provides maps of the probability density function | u α (x, t) | 2 . As in previous results, the function | u α (x, t) | 2 is computed on a uniformly sample grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colormap after normalisation giving a map of | u α (x, t) | 2 ∈ [0, 1], the maps having been enhanced using the Matlab histogram equalisation function histeq.
The results illustrated in Figure 7 show that oscillations in time are prevalent when α → 1 as should be expected, when the solution approaches that of the solution presented in Section 9.1, compounded in Figure 5. However, the time oscillations do not reflect that of a Fresnel diffraction pattern and are of a significantly lower frequency for the same initial condition u 0 (x) = exp(ikx), ∀x < 0. On the other hand, as α → 0, the probability density reflects the standing fringe pattern illustrated in Figure 6, for example, for the space fractional case, but with a greater dissipation in both space and time.

The Quantum Shutter Problem for the Relativistic and Semi-Relativistic Case
In the final part of this paper, we consider the quantum shutter problem for the relativistic and then the semi-relativistic case. In the latter case, a phenomenological equation is introduced, which is taken to describe a quantum field which is neither nonrelativistic or fully-relativistic. In this regard, and, in reference to Figure 1, we consider a beam of 'semi-relativistic' particles with velocities of the order of 50-70% light speed. Examples of such particles include the delocalised electrons that exist in Graphene [42]. This is because the charge carriers in Graphene show linear, rather than quadratic, dependence of energy E on momentum p [43,44]. This dependence indicates that they adhere more to the relativistic energy equation E 2 = p 2 + m 2 for a charge carrier with mass m (for h = c 0 = 1) rather than conforming to the non-relativistic energy equation E = p 2 /2m.

Quantum Shutter Problem for the Klein-Gordon Equation
The problem is effectively the same as that considered in Section 9 but instead of solving the problem for the one-dimensional Schrödinger equation, we consider the (onedimensional) Klein-Gordon Equation given by (forh = c 0 = m = 1) subject to initial conditions for u 0 (x) and u (1) 0 (x). The fundamental solution to this equation is the same as that given by Equation (48), but where the Green's function is now the solution of and given by [45] where J 0 is the Bessel function (of the first kind and of order zero). The fact that this solution depends on a convolution with a Bessel function ensures that the field is characterised by oscillations in both space and time. Two examples of this are given in Figure 8.

Quantum Shutter Problem for the Fractional Klein-Gordon Equation
In the context of the solution to the Klein-Gordon equation considered in the previous section, the following 'semi-relativistic' equation is proposed for 0 < α < 1 where ξ = −it [46]: 1 2m (75) is an equation for the semi-relativistic field u(x, t), constructed so that it reduces to the homogeneous Schrödinger equation when α = 0 and to the Klein-Gordon equation when α = 1. Note that Equation (75) is a phenomenology as are the Schrödinger and Klein-Gordon equations. With reference to Figure 1, the aim is to solve Equation (75) for the same initial u(x, 0) = exp(ikx), ∀x ≤ 0 with m = 1, when it can be written in the form Following the approach presented in Section 9.3, by introducing the memory function we can develop a Green's function solution given by The problem is to evaluate the Green's function for this case, a problem which is now addressed.
Fourier-Laplace transformations allow us to study the Green's function in Fourier-Laplace space, when it can be written as where Inverse Fourier transformation, coupled with application of the convolution theorem, yields The purpose of writing the Green's function in this way is so that one can make use of the property that (using the convolution of two functions application available at [47]) Moreover, using the same application available at [47], it can be shown that and similarly that so that, in general, for n convolutions, we can write Thus, through induction, for high frequency fields (which is the case for the semirelativistic case), we can consider an approximate expression for the Green's function given by (for | Ω |>> 1) Consider the case for x > 0. In order to obtain an expression for the time-dependent Green's, we are required to obtain the inverse Laplace transform of G β (x, s). The function i exp(iΩx)/2Ω can be evaluated using the approach taken in Section 9.3. The function exp(iβ 2 x/2Ω) can be evaluated using the same approach as that given in Appendix B. Combining the results, and using the convolution theorem for Laplace transforms, this yields the result However, since the inverse Laplace transform is being evaluated for | x |<< 1, the compatible result allows the second term of the Green's function for x ≤ 0 to be eliminated, leaving the result Two examples of the solution given by Equation (76) are shown in Figure 9, which provides maps of the probability function | u β (x, t) | 2 . The maps are computed using a uniformly sample grid consisting of 10 3 × 10 3 elements and displayed using the Matlab jet colour map after normalisation giving a map of | u β (x, t) | 2 ∈ [0, 1], both maps having been enhanced using the Matlab histogram equalisation function histeq. Figure 7 implies that, for the semi-relativistic model considered, space-time oscillation is not prevalent for α → 0, and the probability density decays relatively slowly. However, as α → 1, the probability density oscillates in space with a period equal to that of the incident wave-field while dissipating in time.
The solution given by Equation (76) implicitly assumes that u (1) 0 (x) = 0. However, by defining the memory function as x may be define thus:

Conclusions: Summary and Further Work
Coupled with a review of the more relevant materials, and, coupled with appropriate references, the principal focus of this paper has been to study the quantum shutter problem for fractional quantum fields. In this context, both the space-fractional and time-fractional Schrödinger equations have been considered. Furthermore, the same shutter problem has been considered for a phenomenological equation, which represents the case when the quantum 'system' is semi-relativistic.
In all cases, the term 'shutter problem' refers to the case when a shutter is opened instantaneously at time t = 0, subject to an initial condition u 0 (x) = exp(ikx), ∀x ≤ 0. The problem is then to evaluate the probability density of the quantum field in the positive half space for all t > 0.

Summary
The study reported in this paper represents an extension of the shutter problem to the field of fractional quantum dynamics. The paper has attempted to explain the models and solutions thereof in terms of a broader connectivity with classical diffraction phenomena and the classical and fractional diffusion equations. This has been done in order to unify the physical models associated with time-independent and then time-dependent phenomena with regard to the interaction of a field with an aperture which is open over all time or opened instantaneously at t = 0, respectively.
In order to contextualise the fractional dynamical models considered, the associated non-fractional models and solutions thereof have been presented. This has been done in order to produce an integrated framework, and to provide reference points to prior art. In both cases, the essential key to developing a solution is to evaluate the Green's functions for the linear differential operator, fractional or otherwise. Once this function has been evaluated, its convolution with some initial condition provides the basic solution. In this context, the paper has explored solutions to the following: Example results have been presented to illustrate the behaviour of the solutions as a function of space and time. This has been done by providing enhanced colour maps of the solutions to illustrate the dynamical behaviour of the probability density. In order to provide an example of the Matlab code used to obtain the result presented, Appendix C provides the .m function used to generate the maps given in Figure 5. This function is just one of many that have been generated to output the field maps presented throughout this paper. It is provided as a basic exemplar for the sake of interested readers who may wish to repeat the results given and extend them further.
In the case of the time-fractional problems considered, the issue of defining an initial condition for a fractional time derivative has been addressed through the application of a memory function. Here, a fractional time derivative is represented in terms of a first derivative in time, convolved with the memory function. The memory function is chosen so that application of the Laplace transform yields an expression which incorporates the required initial condition for the field at t = 0. This provides an explicit indication of the fact that fractional time derivatives are consistent with modelling a memory associative system.

Future Work
The analysis presented in this paper has focused on one-dimensional models, providing (x, t) maps to illustrate the one-dimensional dynamical behaviour of the solutions derived. Thus, an obvious extension is to consider the same approach for multi-dimensional models for both space-and time-fractional cases. In this respect, a further, and, more quantitative, analysis is required, based on the methods presented and their extension to multi-dimensional models, subject to a variety of initial conditions. In other words, the results presented are only illustrative examples of a wide range of results that characterise dynamical behaviour based on variations of the fractional exponents γ and α, for the space-fractional and time-fractional derivatives, respectively.
Another area of interest is to compare the multiple scattering solution provided in Appendix A with the time-fractional diffusion equation for the operator ∇ 2 − ∂ 1+α t , 0 < α < 1. In principle, this operator describes the intermediate case when scattering interactions are neither fully diffusive (α = 0) or fully 'propagative' α = 1. In this context, such an operator can be used to model a scenario where multiple scattering interactions occur for a limited number of scattering sites; where the density of scatterers reduces to 0 as α → 1, yielding a model for propagation in a homogeneous half-space. If all scattering sites are taken to be Coulomb potentials, then the results can then be compared directly with the solution given in Appendix A, subject to ∇ 2 V(r) = 0 as discussed in Corollary A2. In this respect, the development of an exact scattering solution for the space fractional Schrödinger equation may be undertaken using a similar approach.
The numerical results presented have been based on the use of a convolution with the corresponding Green's function, a function that, for both space-and time-fractional models, is based on an approximation in which the field is taken to be observed close to the shutter when | x | 2 << 1 and where t → 0. These conditions have been used so that the Green's function has a similar form to the Green's function for the classical diffusion equation. However, there is another approach to developing a numerical evaluation of the fractional quantum fields that can be considered. This is to use the approach presented in Sections 3 and 3.3, given that the Mittag-Leffler function can be computed numerically [30], coupled with the application of an inverse fast Fourier transform.
For the time fractional Schrödinger equation, written in the form This solution may then be investigated by computing the Mittag-Leffler function numerically, applying an inverse fast Fourier transform to the result, and convolving the output with the initial condition.
For the time fractional Klein-Gordon equation considered in Section 10.2, if we consider the equation (where β = 1 + α) subject to the initial conditions u 0 (x) and u (1) 0 (x), then a solution based on an application of the Mittag-Leffler function can thus be constructed: A numerical evaluation of the fractional quantum field can then be investigated by computing the functions E α,1 and E α,0 numerically [30], applying an inverse fast Fourier transform and convolving the results with the initial condition(s).
Finally, an investigation should be undertaken into the time-fractional Dirac equation for the quantum shutter problem, following the approach taken in [13], for example.

Appendix A. Exact Forward and Inverse Schrödinger Scattering Solutions
The results provided in this appendix are based on an attempt to develop an exact solution to the forward and inverse Schrödinger scattering problem.

Statement of the Problem
For r ∈ R n , the solution to the time-independent Schrödinger equation (∇ 2 + k 2 )U(r, k) = V(r)U(r, k) has a fundamental solution given by the Lippmann-Schwinger equation where u i is the incident wave function which is taken to be a solution to the equation and U s (r, k) is the scattered field given by U s (r, k) = g(r, k) ⊗ r V(r)U(r, k) where g(r, k) is the free-space and out-going Green's function which is the solution of (∇ 2 + k 2 )g(r, k) = δ n (r) On the basis of Equation (A4), the forward scattering problem is 'Given V(r) evaluate U s (r, k)'. The corresponding the inverse scattering problem is 'Given U s (r, k) evaluate V(r)'. Both problems ideally require unconditional exact solutions.
To progress further, we now introduce a crucial condition, which is to let If we then write this result in terms of U(u, k), application of the convolution and product theorems yields the Fourier-space equation 1 (2π) n V(u, k) ⊗ u U i (u, k) = U(u, k) W(u, k) − g(u, k) U(u, k)[ W(u, k)] 2 (A8) where V(u) denotes the n-dimensional Fourier transform of V(r), and U i (u) denotes the n-dimensional Fourier transform of U i (r). Thus, given Equations (A2) and (A4), we have U(r, k) = U i (r, k) + g(r, k) ⊗ r V(r)U(r, k) or, given Equation (A5), and, using the convolution theorem W(u, k) = U(u, k) − U i (u, k) g(u, k) U(u, k) Substituting this expression for W(u, k) into Equation (A8) and noting that U(u, k) − U i (u, k) = U s (u, k), we obtain U i (u, k) U s (u, k) U i (u, k) + U s (u, k) = 1 (2π) n g(u, k)[ V(u, k) ⊗ u U i (u, k)] Finally, taking the inverse Fourier transform, we can write this result in the form U(r, k) = g(r, k) ⊗ r V(r)U i (r, k) where U(r, k) where U B s (r, k) is the scattered field under the condition applied.
This result is equivalent to the Born approximation discussed in Section 5.1, i.e., Equation (34) under condition Equation (33). The exact scattering solution given by Equation (A9) therefore reduces to the Born approximation under the condition that the scattered field is much weaker than the incident field. Corollary A2. Equation (A9) depends explicitly on the condition that Coupled with the condition that (∇ 2 + k 2 )U i (r, k) = 0, it is therefore imperative that we consider the ramifications of this condition on the solution developed-in other words, what are the scattering potentials V(r) for which Equation (A12) is applicable.
The first approach for doing this is to expand Equation (A12) and write it in the form ∇ 2 V(r) + 2∇V(r) · ∇ ln U i (r, k) = 0 It is then apparent that, for any incident plane wave-field U i (r, k) = exp(±ik · r) with a constant amplitude A or amplitude spectrum A(k), we can write ∇ 2 V(r) ± 2ik · ∇V(r) = 0 (A13) which has the general solution (for arbitrary constants C 1 and C 2 ) V(r) = C 1 exp(∓2ik · r) + C 2 Thus, one condition on the scattering potential is that it is a plane wave oscillating at twice the frequency of the incident plane wave. However, if the potential is taken to be real only, then Equation (A13) suggests that the scattering potential is constrained to satisfying the Laplace equation, i.e., ∇ 2 V(r) = 0 (A14) This allows as to consider the Coulomb potential for a point charge (ignoring any scaling constant) However, this condition on the scattering function also allows for the sum of such potentials to be considered, when V(r) = ∞ ∑ n=1 1 | r − r n | , r = r n and where the position vector at which the multiple potentials occur can be regular or random. Thus, the exact scattering solution may, in principle, be applied to a medium consisting of a random distribution of point charges such as a plasma-electrons scattering from electrons. Other potentials V(r) that satisfy Equation (A14) are available depending on the dimension of r ∈ R n . This includes the harmonic functions and their linear combination. For a potential of compact support, the boundary conditions for which a solution to Equation (A14) is required need to be compatible with those that provide a solution for U(r, k) which conforms to Equation (A2). These conditions include U(r, k) = U i (r, k) and ∇U(r, k) = ∇U i (r, k) on the surface of V(r).
Corollary A3. The condition on the scattering potential can be investigated further on the basis of convolving Equation (A12) with the free-space Green's function when we can write g(r, k) ⊗ r ∇ 2 [V(r)U i (r)] + k 2 g(r, k) ⊗ r V(r)U i (r) = g(r, k) ⊗ r ∇ 2 [V(r)U i (r)] + k 2 U(r, k) = 0 where U(r, k) is given by Equation (A9) from which the exact scattered field U s (r, k) can be evaluated. Fourier transformation then yields the u-space equation so that, upon rearrangement and inverse transformation, we can write V(r) = k 2 U * i (r, k) | U i (r, k) | 2 [U(r, k) + k 2 h(r) ⊗ r U(r, k)], h(r) ↔ − 1 u 2 = k 2 U * i (r, k)[U(r, k) + k 2 h(r) ⊗ r U(r, k)], U i (r, k) = exp(±ik · r) which provides an exact inverse scattering solution.
The function h(r) depends on the dimensionality of the solution that is required, and is given by the following: For r ∈ R n , and, r ≡| r |, Remark A1. In terms of using Theorem A1 to compute the exact scattered field, for potentials that conform to those presented in Corollary A2, the exact scattered field is computed through application of Equation (A9) directly. For scattering functions V(r) that do not conform to the conditions presented in Corollary A2, then, from Equation (A15) and Relationship (A10), we can construct an expression for the Fourier transform of the scattered field given by , → 0 We can then consider the asymptotic solution U s (r, k) ∼ − 1 k 2 g(r, k) ⊗ r ∇ 2 [V(r)U i (r, k)], k → ∞