Exact and Paraxial Broadband Airy Wave Packets in Free Space and a Temporally Dispersive Medium

1. Introduction

The Airy beam, a remarkable finite-energy solution to the paraxial equation was first formulated analytically by Siviloglou and Christodoulides [1] and subsequently demonstrated experimentally by Siviloglou, Broky, Dogariu and Christodoulides [2]. Their work was motivated by the infinite-energy nonspreading accelerating Airy solution to the Schrödinger equation introduced by Berry and Balazs [3] in the context of quantum mechanics. The Airy beam is slowly diffracting while bending laterally along a parabolic path even though its centroid is constant, it can perform ballistic dynamics akin to those of projectiles moving under the action of gravity, and it is self-healing, that is, it regenerates when part of the “generating” aperture is obstructed; this is due to the reinforcement of the main lobe by the side lobes.

An important question is whether spatiotemporally localized (i.e., pulsed) versions of Airy beams are feasible. For luminal solutions, this question has been answered affirmatively by Saari [4], Valdmann, Piksarv, Valta-Lukner and Saari [5] and Kaganovsky and Heyman [6] within the framework of the paraxial geometry. The time-diffraction technique introduced by Porras [7,8] recently has been motivated by Besieris and Shaarawi [9] in terms of the Lorentz invariance of the equation governing the narrow angular spectrum and narrowband temporal spectrum paraxial approximation and has been used to derive finite-energy spatiotemporally confined subluminal, luminal, and superluminal Airy wave packets. The goal in this article is to provide exact (i.e., non-paraxial) finite-energy broadband spatiotemporally localized Airy solutions (a) to the scalar wave equation in free space; (b) in a dielectric medium moving at its phase velocity; (c) in a lossless second order temporally dispersive medium.

2. Finite-energy (3+1)D spatiotemporally localized Airy splash mode solution of the scalar wave equation in free space

Consider the $\left(3+1\right)D$scalar wave equation in free space, viz.,

$$\left(\frac{{\partial }^2}{\partial X^2}+\frac{{\partial }^2}{\partial Y^2}+\frac{{\partial }^2}{\partial Z^2}-\frac{{\partial }^2}{\partial T^2}\right)\Psi \left(X,Y,Z,T\right)=0,$$

(1)

written in terms of the nondimensional variables $X=x/x_0, Y=y/x_0, Z=z/x_0,$ and $T=ct/x_0.$ The introduction of the characteristic variables ${\Lambda }_{\pm }=Z\mp T$ changes Eq. (1) as follows:

$$\left(\frac{{\partial }^2}{\partial X^2}+\frac{{\partial }^2}{\partial Y^2}+4\frac{{\partial }^2}{\partial {\Lambda }_{+}\partial {\Lambda }_{-}}\right)\Psi \left(X,Y,{\Lambda }_{+},{\Lambda }_{-}\right)=0.$$

(2)

A specific solution of this equation is the infinite energy accelerating Airy wavepacket

$$\Psi \left(X,{\Lambda }_{+},{\Lambda }_{-}\right)=\exp \left[i{\Lambda }_{-}\right]e^{-i\frac{1}{24}{\Lambda }_{+}\left(-6X+\frac{1}{4}{\Lambda }_{+}^2\right)}Ai\left(X-\frac{1}{16}{\Lambda }_{+}^2\right).$$

(3)

This wavefunction moves in a parabolic trajectory along the characteristic variable 

In 1910, Bateman [10,11] discovered a transformation, more general than a conformal change of the metric, which could be used to transform solutions of Maxwell equations into similar ones. In the case of the scalar wave equation, the Bateman transformation in $\left(3+1\right)D$ assumes the form

$${\Psi }_1\left(X,Y,{\Lambda }_{+},{\Lambda }_{-}\right)=\frac{1}{{\Lambda }_{-}}\Psi \left[\frac{X}{{\Lambda }_{-}},\frac{Y}{{\Lambda }_{-}},-\frac{1}{{\Lambda }_{-}},\frac{X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}}{ {\Lambda }_{-}}\right].$$

(4)

The function ${\Psi }_1\left(X,Y,{\Lambda }_{+},{\Lambda }_{-}\right)$ also obeys the scalar wave equation (2).

The Bateman transformation is applied twice to the solution given in Eq. (3). These two sequential operations result in the following new solution to Eq. (2):

$$\begin{array}{l}{\Psi }_2\left(X,Y,{\Lambda }_{+},{\Lambda }_{-}\right)=\frac{1}{X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}}Ai\left[-\frac{{\Lambda }_{+}^2-16X\left(X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}\right)}{16{\left(X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}\right)}^2}\right]\\ \times \exp \left[-i\frac{-{\Lambda }_{+}^3+24X{\Lambda }_{+}\left(X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}\right)+96{\Lambda }_{-}{\left(X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}\right)}^2}{96{\left(X^2+Y^2+{\Lambda }_{+}{\Lambda }_{-}\right)}^3}\right].\end{array}$$

(5)

Next, this expression is complexified by means of the changes ${\Lambda }_{+}\to {\Lambda }_{+}-i{a}_1, {\Lambda }_{-}\to {\Lambda }_{-}+i{a}_2,$ where $a_{1,2}$ are two positive parameters. Consequently, one obtains

$${\psi }_2\left(R,\varphi ,{\Lambda }_{+,}{\Lambda }_{-,}; a_1,a_2\right)=B^{-1}\exp \left(-i\frac{Q}{96B^3}\right)Ai\left[\frac{{\left(a_1+i{\Lambda }_{+}\right)}^2+BR\cos \varphi }{16B^2}\right];$$

(6)

with the exponent $Q$ given as

$$\begin{array}{l}Q=-i{\left(a_1+i\Lambda {}_{+}\right)}^3-i24R\cos \varphi \left(a_1+i\Lambda {}_{+}\right)B+i96\left(a_2-i\Lambda {}_{-}\right)B^2,\\ B=R^2+\left(a_1+i\Lambda {}_{+}\right)\left(a_2-i\Lambda {}_{-}\right)\end{array}$$

(7)

in cylindrical coordinates $\left(X=R\cos \varphi , Y=R\sin \varphi \right).$ This is a finite-energy $\left(3+1\right)D$ spatiotemporally localized luminal wave packet belonging to the class of splash modes studied by Ziolkowski [12]. It will be referred to as the Airy splash mode.

The parameter $A$ in the Bateman conformal transformation is arbitrary. On the other hand, the free positive parameters $a_1$ and $a_2$ entering the solution given in Eq. (7) are critical. As discussed originally by Ziolkowski ([12]; see also [13]), their presence ensures finite energy. Their relative values measure the size of the forward and backward wave components. Only when $a_1>>a_2$ the backward components are minimized, and the solution is almost undistorted. This is further explained in [14], where it is shown that very close replicas of localized waves, such as the one in Eq. (6), can be launched causally from apertures constructed on the basis of the Huygens principle.

Figure 1 shows surface plots of the intensity of Airy splash mode versus ${\Lambda }_{+}$ and $X$ for various values of $T$ the latter defined by the relationship ${\Lambda }_{-}={\Lambda }_{+}+2T.$The parameters $a_1$ and $a_2$ have the values $8\times {10}^{-2}$ and $10,$ respectively. The wave packet is relatively undistorted because $a_1>>a_2.$

The finite-energy wavefunction $U_2^{+}\left(X,Y,{\Lambda }_{+},Z\right)={\Psi }_2\left(X,Y,{\Lambda }_{+},2Z; a_1,a_2\right)$ obeys the paraxial forward pulsed beam equation

$$\left(\frac{{\partial }^2}{\partial X^2}+\frac{{\partial }^2}{\partial Y^2}+\frac{{\partial }^2}{\partial {\Lambda }_{+}\partial Z}\right)U_2^{+}\left(X,Y,{\Lambda }_{+},Z\right)=0.$$

(8)

The transition from Eq. (3) to (6) is effected by means of the modified complexification ${\Lambda }_{+}\to {\Lambda }_{+}-i{a}_1, {\Lambda }_{-}\to 2Z+i{a}_2.$ As a result, one obtains broadband splash mode-type spatiotemporally localized wave solutions to the paraxial equation.

3. Finite-energy accelerating spatiotemporally localized Airy wavepacket solution to the scalar equation in free space

A solution to Eq. (2) is sought in the form

$$\psi \left(X,{\Lambda }_{+},{\Lambda }_{-};\beta \right)=\exp \left(i\beta {\Lambda }_{-}\right)\phi \left(X,{\Lambda }_{+};\beta \right).$$

Then, $\phi \left(X,{\Lambda }_{-};\beta \right)$ is governed by the parabolic equation

$$\left(i4\beta \frac{\partial }{\partial {\Lambda }_{+}}+\frac{{\partial }^2}{\partial X^2}\right)\phi \left(X,{\Lambda }_{+};\beta \right)=0.$$

(9)

A solution to this equation is the “accelerating beam”

$$\begin{array}{l}\phi \left(X,{\Lambda }_{-};\beta \right)=\exp \left[-\frac{1}{12}\left(2a+i{\Lambda }_{+}\right)\left(2a^2-6\sqrt{2\beta }X-i4a{\Lambda }_{+}+{\Lambda }_{+}^2\right)\right]\\ \times Ai\left(\sqrt{2\beta }X+ia{\Lambda }_{+}-\frac{{\Lambda }_{+}^2}{4}\right).\end{array}$$

(10)

With $a$ a small positive parameter and ${\Lambda }_{+}$ replaced by $Z,$ $\phi \left(X,Z\right)$ is essentially the finite-energy monochromatic paraxial accelerating Airy beam solution introduced by Siviloglou and Christodoulides [1]. In contrast, $\psi \left(X,{\Lambda }_{+},{\Lambda }_{-};\beta \right)=\exp \left(i\beta {\Lambda }_{-}\right)\phi \left(X,{\Lambda }_{+};\beta \right)$ is not a finite-energy solution to the scalar wave equation. To achieve a finite-energy spatiotemporal solution an appropriate superposition over the free parameter $\beta$ of the form

$$\psi \left(X,{\Lambda }_{+},{\Lambda }_{-}\right)={\displaystyle \int d\beta }F\left(\beta \right)\exp \left(i\beta {\Lambda }_{-}\right)\phi \left(X,{\Lambda }_{+};\beta \right)$$

(11)

must be undertaken. Such a superposition can be brought to the form

$$\begin{array}{l}\psi \left(X,{\Lambda }_{+},{\Lambda }_{-}\right)=Q{\displaystyle {\int }_{-\infty }^{\infty }d\sigma }\exp \left(-{\sigma }^2\right)Ai\left[A\left(\frac{B}{A}-\sigma \right)\right];\\ Q=\frac{1}{\sqrt{a_1-i{\Lambda }_{-}}}\exp [-\frac{1}{6\sqrt{2}}b\left(2a^2-i4a{\Lambda }_{+}+{\Lambda }_{+}^2\right)]\exp \left[\frac{1}{4}{b}^2\frac{X^2}{a_1-i{\Lambda }_{-}}\right],\\ \mathrm{A}=\frac{\sqrt{a_1-i{\Lambda }_{-}}}{\sqrt{2}X}, b=\sqrt{2}\left(a+i\frac{{\Lambda }_{+}}{2}\right),\\ B=\frac{ia{\Lambda }_{+}-\frac{{\Lambda }_{+}^2}{4}+b\frac{X^2}{\sqrt{2\left(a_1-i{\Lambda }_{-}\right)}}}{A}.\end{array}$$

(12)

where $a_1$ is a positive free parameter. The integral in Eq. (12) is an Airy transform [15]. It can be carried out explicitly, yielding the finite-energy accelerating spatiotemporal Airy wavepacket

$$\psi \left(X,{\Lambda }_{+},{\Lambda }_{-}\right)=Q\exp \left[\frac{1}{4A^3}\left(B+\frac{1}{24A^3}\right)\right]Ai\left(\frac{B}{A}+\frac{1}{16A^4}\right).$$

(13)

Figure 2 shows surface plots of the modulus square of the Airy wavepacket versus ${\Lambda }_{+}$ and $X$ for various values of $T,$ the latter defined by the relationship ${\Lambda }_{-}={\Lambda }_{+}+2T.$ The parameters $a$ and $a_1$ have the values $5\times {10}^{-2}$ and $100,$ respectively.

The finite-energy wavefunction $U_2^{+}\left(X,{\Lambda }_{+},Z\right)\equiv {\Psi }_2\left(X,{\Lambda }_{+},2Z; a_1,a_2\right)$ obeys the paraxial forward pulsed beam equation [cf. Eq. (8)] if the replacement ${\Lambda }_{-}\to 2Z$ is made in Eq. (13). As mentioned previously, one obtains then a broadband spatiotemporally localized accelerating Airy solution to the paraxial equation.

4. Finite-energy accelerating broadband Airy wavepacket solution in the presence of temporal dispersion

Basic Equation

Electromagnetic wave propagation in a linear, homogeneous, transparent, dispersive medium is governed by the scalar equation

$${\nabla }^{2}u(\overrightarrow{r},t)+{\beta }_{op}^2(-i\partial /\partial t)u(\overrightarrow{r},t)=0,$$

(14)

if polarization is neglected. In this expression, $u(\overrightarrow{r},t)$ is a real field and ${\beta }_{op}^2(-i\partial /\partial t)$ a real pseudo-differential operator. A physical interpretation of the latter is provided in the frequency domain.; specifically,

$$F\left\{{\beta }_{op}^2(-i\partial /\partial t) u(\overrightarrow{r},t)\right\}={\beta }^2(\omega ) \tilde{u}(\overrightarrow{r},\omega )$$

(15)

where $F\{\cdot \}$ denotes Fourier transformation and $\tilde{u}(\overrightarrow{r},\omega )$ is the Fourier transform of $u(\overrightarrow{r},t)$ with respect to time. The function $\beta (\omega )$ appearing on the right-hand side of Eq. (15) is a real wave number.

For a physically convenient central radian frequency${\omega }_0$, the real field $u(\overrightarrow{r},t)$ is expressed as follows:

$$u(\overrightarrow{r},t)=\phi (\overrightarrow{r},t) \exp \left[-i{\omega }_0\left(t-z/v_{ph}\right)\right] +cc, z\ge 0.$$

(16)

Here, $\phi (\overrightarrow{r},t)$ is a complex-valued envelope function and $v_{ph}={\omega }_0/\beta ({\omega }_0)$ denotes the phase speed in the medium computed at the central frequency. For pulses as short as a single optical period $T_0=2\pi /{\omega }_0,$ and within the framework of the paraxial approximation, the envelope function obeys the following equation [16,17,18]:

$$\left(1+i\frac{{\beta }_1}{{\beta }_0}\frac{\partial }{\partial \tau }\right)\frac{\partial }{\partial z}\phi \left(\overrightarrow{r},\tau \right)-\left(1+i\frac{{\beta }_1}{{\beta }_0}\frac{\partial }{\partial \tau }\right)iD\left(-i\frac{\partial }{\partial \tau }\right) \phi \left(\overrightarrow{r},\tau \right)-\frac{i}{2\beta {}_0}{\nabla }_t^2\phi \left(\overrightarrow{r},\tau \right)=0.$$

(17)

Here, ${\nabla }_{\perp }^2$ denotes the transverse (with respect to $z$) Laplacian operator and $\tau =t-(z/v_{gr})$ corresponds to a moving reference frame, defined in terms of the group speed $v_{gr}=1/{\beta }_1; {\beta }_1\equiv {d\beta (\omega )/d\omega |}_{\omega ={\omega }_0}$. The operator $D$ is given by the expression

$$D\left(-i\frac{\partial }{\partial \tau }\right)={\displaystyle \sum _{m=2}^{\infty }\frac{{\beta }_m}{m!}}{\left(-i\frac{\partial }{\partial \tau }\right)}^m; {\beta }_m=\frac{d^m}{d{\omega }^m}{\beta (\omega )|}_{\omega ={\omega }_0}.$$

(18)

In the sequel only the first term the series will be retained. This approximation results in the equation

$$\left(1+i\frac{{\beta }_1}{{\beta }_0}\frac{\partial }{\partial \tau }\right)\frac{\partial }{\partial z}\phi \left(\overrightarrow{r},\tau \right)-\left(1+i\frac{{\beta }_1}{{\beta }_0}\frac{\partial }{\partial \tau }\right)\left(-i\frac{{\beta }_2}{2}\frac{{\partial }^2}{\partial {\tau }^2}\right) \phi \left(\overrightarrow{r},\tau \right)-\frac{i}{2\beta {}_0}{\nabla }_t^2\phi \left(\overrightarrow{r},\tau \right)=0.$$

(19)

Here ${\beta }_2\equiv {d^2\beta (\omega )/d^2\omega |}_{\omega ={\omega }_0}$ is the second-order index of dispersion. It is positive for normal dispersion and negative for anomalous dispersion.

Accelerating Airy solution 

A solution to Eq. (19) is sought in the form $\phi \left(\overrightarrow{r},\tau \right)=\tilde{\phi }\left(\overrightarrow{r},\Omega \right)\exp \left(-i \tau \Omega \right).$ Furthermore, with $\tilde{\phi }\left(\overrightarrow{r},\Omega \right)=\tilde{\psi }\left(\overrightarrow{r},\Omega \right)\exp \left[i{\beta }_2{\Omega }^{2}z/2\right],$ one obtains the parabolic equation

$$\left[i{\beta }_0\left(1+\frac{{\beta }_1}{{\beta }_0}\Omega \right)\frac{\partial }{\partial z}+\frac{1}{2}{\nabla }_t^2\right]\tilde{\psi }\left(\overrightarrow{r},\Omega \right)=0.$$

(20)

The azimuthally asymmetric expression

$$\tilde{\psi }\left(\rho ,\varphi ,z;\Omega \right)=\exp \left(im\varphi \right)\frac{{\rho }^m}{{\left(a_1+i z\right)}^{m+1}}\exp \left[-\left({\beta }_0+{\beta }_1\Omega \right)\frac{{\rho }^2}{2\left(a_1+i z\right)}\right],$$

(21)

with $a_1$ a positive parameter, satisfies the paraxial equation (20). A spatiotemporal solution to Eq. (19) can be derived by means of the superposition

$$\begin{array}{l}\phi \left(\rho ,\varphi ,z,\tau \right)=\exp \left(im\varphi \right)\frac{{\rho }^m}{{\left(a_1+i z\right)}^{m+1}}\exp \left[-{\beta }_0\frac{{\rho }^2}{2\left(a_1+i z\right)}\right]\\ \times {\displaystyle {\int }_{-\infty }^{\infty }d\Omega }\exp \left[-i \Omega \left(\tau -i{\beta }_1\frac{{\rho }^2}{2\left(a_1+i z\right)}\right)\right]\exp \left[i\frac{{\beta }_2}{2}{\Omega }^{2}z\right]F(\Omega ).\end{array}$$

(22)

A large class of solutions can be obtained by using different spectra $F(\Omega )$. Choosing the spectrum $F(\Omega )=\exp \left(-a_2{\beta }_2{\Omega }^2/2\right)\exp \left(i {\Omega }^3/3\right); a_2>0,$ results in the solution

$$\begin{array}{l}\phi \left(\rho ,\varphi ,z,\tau \right)=\exp \left(im\varphi \right)\frac{{\rho }^m}{{\left(a_1+i z\right)}^{m+1}}\exp \left[-{\beta }_0\frac{{\rho }^2}{2\left(a_1+i z\right)}\right]\\ \times \exp \left[-i\frac{{\beta }_2}{2}\left(i{a}_2+ z\right)\left(\frac{2}{3}{\left(i{a}_2+ z\right)}^2{\beta }_2^2/4-i{\beta }_1\frac{{\rho }^2}{a_1+i z}+\tau \right)\right]\\ \times Ai\left[-{\beta }_2^2{\left(i{a}_2+ z\right)}^2/4+i{\beta }_1\frac{{\rho }^2}{a_1+i z}-\tau \right].\end{array}$$

(23)

This is a finite-energy accelerating Airy-Gaussian wavepacket. Figure 3 shows the intensity versus $\tau$ and $\rho$ at $z=0$ for $m=0.$ Figure 4 shows surface plots of the modulus of the azimuthally symmetric wavepacket versus $\tau$ and $z$ for three values of $\rho .$

5. Finite-energy accelerating broadband Airy wavepacket solution in a dielectric moving at its phase velocity

An equation arising in the case of a dielectric medium moving at its phase velocity is given by [19]

$$\left(\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}-\frac{2}{v_0}\frac{{\partial }^2}{\partial t\partial z}-\frac{1+{\beta }^2}{v_0^2}\frac{{\partial }^2}{\partial t^2}\right)u\left(x,y,z,t\right)=0.$$

(24)

Here, $v_0$ denotes the phase speed, $\beta =v_0/c<1$ and $u\left(x,y,z,t\right)$ stands for the longitudinal components $E_z$ and $H_z$ in the absence of sources. On the other hand, the equation of acoustic pressure under conditions of uniform flow is given as follows:

$$\left[{\nabla }^2-\frac{1}{u_0^2}{\left(\frac{\partial }{\partial t}+\overrightarrow{u}\cdot \nabla \right)}^2\right]p\left(\overrightarrow{r},t\right)=0.$$

(25)

Here, $u_0$ is the speed of sound in the rest frame of the medium and $\overrightarrow{u}$ is the uniform velocity of the background flow. In the special case where $\overrightarrow{u}=u{\overrightarrow{a}}_z$ and $u=u_0,$ the resulting equation for the acoustic pressure is isomorphic to Eq. (24).

Assuming independence on the transverse variable $y$ in Eq. (24), the ansatz

$$u\left(x,z,t\right)=\varphi (x,z)\exp \left[i\omega \frac{1+{\beta }^2}{2v_0}\left(z-\frac{2v_0}{1+{\beta }^2}t\right)\right]$$

(26)

gives rise to the equation

$$\left(i\frac{\omega }{v_0}\frac{\partial }{\partial z}+\frac{1}{2}\frac{{\partial }^2}{\partial x^2}\right)\varphi \left(x,z\right)=0.$$

(27)

This is an exact parabolic equation, in contradistinction to the paraxial approximation of the Helmholtz equation associated with the temporal Fourier transform of the ordinary scalar wave equation. In addition to the well-known beam solutions of the usual paraxial equation, Eq. (27) has the following “accelerating” one [1]:

$$\begin{array}{l}g\left(x,z\right)=\exp \left[-\frac{1}{12}\left(2a+iz\right)\left(2a^2-i4az+z^2-6x\sqrt{\omega /v_0}\right)\right]\\ \times Ai\left(x\sqrt{\omega /v_0}-\frac{z^2}{4}+iaz\right).\end{array}$$

(28)

Here, $Ai(\cdot )$ denotes the Airy function and the positive parameter $a$ ensures finite energy for the monochromatic solution. The beam follows a parabolic trajectory upon propagation.

Finite-energy broadband pulse solutions can be obtained by using the solution (28) together with the ansatz (26) and undertaking a superposition with respect to the frequency $\omega$ (see, e.g., Ref. [15]). A specific spatiotemporal broadband solution is given as follows:

$$\begin{array}{l}u\left(x,z,t\right)=\frac{1}{\sqrt{a_1-iz}}\exp \left[\frac{b^2{x}^2}{4\left(a_1-iz\right)}\right]\exp \left[\frac{1}{4A^2}\left(Y+\frac{1}{24A^4}\right)\right]\\ \times Ai\left[Y+\frac{1}{16A^4}\right];\\ Y=\frac{b{x}^2}{\sqrt{2}\left(a_1-iz\right)}-\frac{z^2}{4}+iaz;\\ A=\frac{\sqrt{a_1-i\eta }}{2x}; b=\sqrt{2}\left(a+i\frac{z}{2}\right);\\ \eta =\left(1+{\beta }^2\right)\left(z-\frac{2v_0}{1+{\beta }^2}t\right).\end{array}$$

(29)

Figure 5 shows surface plots of the modulus of the wavepacket versus $\tau =t-z\left(1+{\beta }^2\right)/\left(2v_0\right)$ and $x$ for four values of $z.$

6. Concluding Remarks

Space-time paraxial solutions based on the narrow angular spectrum obey the pulsed beam equation mentioned in Section 2 and Section 3. Broadband as well as narrowband spatiotemporal luminal solution to the beam equation were addressed in [4,5,6]. More recently, work on nonluminal spatiotemporally confined Airy wave packets has appeared in the literature [20,21,22]. The time diffraction method [7,8], motivated in terms of the Lorentz invariance of the narrow angular spectrum and narrow temporal spectrum paraxial equation [9], allows the derivation of finite-energy spatiotemporally localized subluminal, luminal and superluminal Airy wave packets.

In this article, four distinct finite-energy broadband spatiotemporally confined Airy-type solutions have been presented. In Section 2, one starts with a solution to the scalar wave equation in the form of an Airy pulse obeying the parabolic equation along one of the characteristic variables of the scalar wave equation, multiplied by a plane wave involving the second characteristic variable. Two sequential applications of the Bateman conformal transformation result in a finite-energy broadband splash mode-type spatiotemporally localized Airy solution to the scalar wave equation in free space. A different exact broadband solution to the scalar wave equation in free space is derived in Section 3. On starts with an infinite energy solution consisting of the product of a plane wave involving one of the characteristic variables of the scalar wave equation and a variant of the Siviloglou-Christodoulides Airy solution obeying the parabolic equation along the second characteristic variable. An integration over a free parameter entering the solution and use of the Airy transform yields a different type of a finite-energy broadband spatiotemporally localized Airy solution to the scalar wave equation in free space.

Different types of Airy solutions in the presence of second-order temporal dispersion have appeared in the literature. The simplest is the analog of the monochromatic Airy beam involving the axial variable $z$ and the transverse variable $x$.The former involves the axial variable $z$ and the “transverse” variable $\tau =t-z/v_g$, where $v_g$denotes the group velocity. Another separable solution is of the form $\psi \left(x,y,z,\tau \right)=\varphi \left(z,\tau \right)u\left(x,y\right).$ The $\left(3+1\right)D$ symplectic solution $\psi \left(\rho ,\varphi ,z,\tau \right)$ of Eq. (19) in cylindrical coordinates given in Eq. (23) of Section 4 is much more complicated. It is a finite-energy paraxial broadband localized Airy solution.

In Section 5, broadband finite-energy spatiotemporally localized Airy solutions are presented to equations arising in two different physical settings: (a) in the case of a dielectric medium moving at its phase velocity; (b) for acoustic pressure under conditions of uniform flow.