1 Introduction

Nondiffracting wave-packets are known solutions to Schrödinger equation [1, 2] and Maxwell equations [3, 4], and several experimental studies have been dedicated to them [3,4,5,6,7,8,9,10,11]. These solutions have been observed in experiments and some of their properties have been determined. Some of their features, such as acceleration, are due to the fact that they exhibit non-vanishing Bohm potential [12]. In fact, in Ref. [12] it is established that the usual dispersion relation for massless fields \(k_\mu k^\mu =0\) is valid only when the Bohm potential associated to the solution vanishes, but otherwise the modified right-hand side is exactly equal to the Bohm potential, which may be either positive or negative. It is important to stress that similar results have been found by different authors in the last seventy years [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35], in addition to the ones already cited above.

Nevertheless, studying the implications of the possible existence of nondiffracting gravitational wave propagation seems to have been neglected so far. To understand its importance on gravitational wave behavior, research in this area is needed. This work intends to help filling that gap.

This manuscript is devoted to investigate accelerating and non-accelerating nondiffracting gravitational waves. We show below how accelerating wave packets are solutions of linearized gravity. Additionally, we explicitly obtain exact non-accelerating, nondiffracting gravitational wave solutions that, as far as we know, have not been considered before. This result consists on a new solution for the Ehlers-Kundt waves [36, 37].

2 Accelerating linearized gravitational waves

First, let us show that accelerating gravitational waves are solutions to linearized gravity. As usual, let us consider a small perturbation \(h_{\alpha \beta }\) of a background flat space-time metric \(\eta _{\alpha \beta }=(-1,1,1,1)\), such that the total metric is \(g_{\alpha \beta }=\eta _{\alpha \beta }+h_{\alpha \beta }\). It is straightforward to show that the trace-reversed metric perturbation \({\overline{h}}_{\alpha \beta }=h_{\alpha \beta }-({1}/{2})\eta _{\alpha \beta }\, h_\mu ^\mu \) [37], satisfies the free wave equation

$$\begin{aligned} \partial _\mu \partial ^\mu {\overline{h}}_{\alpha \beta }=0\, , \end{aligned}$$
(1)

under the Lorenz gauge \(\partial _\beta {\overline{h}}_{\alpha }^{\beta }=0\).

In order to exemplify the accelerating solutions for gravitational waves, let us focus in the evolution of the wave polarization \({\overline{h}}_{zz}={\overline{h}}_{zz}(t,x,y)\). Notice that \(\partial _z {\overline{h}}_{zz}=0\), thus satisfying the Lorenz gauge. Now we assume the following form for such polarization

$$\begin{aligned} {\overline{h}}_{zz}(t,x,y)={\mathcal {G}}(\zeta ,y)\exp \left( i k\, \eta \right) \, , \end{aligned}$$
(2)

where k is an arbitrary constant, \(\zeta =x-t\) and \(\eta =x+t\) [38]. Thus, as \({\overline{h}}_{zz}\) evolves according to Eq. (1), we are able to find that \({\mathcal {G}}\) evolves as

$$\begin{aligned} 4ik\frac{\partial {\mathcal {G}}}{\partial \zeta } +\frac{\partial ^2{\mathcal {G}}}{\partial y^2}=0\, . \end{aligned}$$
(3)

This equation is equivalent to a Schrödinger equation for a free particle. It is very well-known that this equation allows solution for accelerating wavepackets [1]. In our case, the accelerating gravitational wavepackets that solve Eq. (3) are

$$\begin{aligned} {\mathcal {G}}(\zeta ,y)={\text{ Ai }}\left( 2 k\, y-k^2\zeta ^2\right) \exp \left( 2ik^2\zeta \, y-\frac{2i}{3}k^3\zeta ^3 \right) \, , \end{aligned}$$
(4)

where \({\text{ Ai }}\) is an Airy function. Thus, the nondiffracting accelerating gravitational wavepacket is the real part of full solution (2). This Airy wavepacket has an acceleration equal to k/2, deflecting its trajectory in a parabolic path in the \(\zeta -y\) plane [3, 39]. Airy waves have been extensively studied in the realm of light propagation in media [40,41,42,43]. Beyond the theoretical importance of general accelerating wavepackets, which have encompassed several decades of research, the most remarkable fact is that they have been observed [9, 11, 44].

Wave solution (4) has infinite energy, which is similar to what happens with plane waves. However, structured nondiffracting gravitational Airy wavepackets with finite energy can be constructed using Gaussian beams [39, 45, 46]. The finite energy gravitational Airy wavepacket is given by

$$\begin{aligned}&{\mathcal {G}}(\zeta ,y)\nonumber \\&\quad ={\text{ Ai }}\left( 2 k\, y-k^2\zeta ^2+ i a k\zeta \right) \, \nonumber \\&\qquad \times \exp \left( a k y-a k^2 \zeta ^2+\frac{i}{4} a^2 k \zeta +2ik^2\zeta \, y-\frac{2i}{3}k^3\zeta ^3 \right) \, , \end{aligned}$$
(5)

in terms of an arbitrary real parameter a that makes possible the normalization of the solution. The expression given in (5) can be straightforwardly proved to satisfy Eq. (3), and when \(a=0\) we recover the wavepacket (4).

We obtain that \(\int dy\ |{\mathcal {G}}|^2=\int dy\ {\mathcal {G}}^*{\mathcal {G}}\), is independent of \(\zeta \), by using Eq. (3). Then we get that this wavepacket is square integrable along \(\zeta =0\) (\(x=t\)) as

$$\begin{aligned} \int _{-\infty }^\infty dy |{\mathcal {G}}(\zeta ,y)|^2 =\int _{-\infty }^\infty dy |{\mathcal {G}}(0,y)|^2=\frac{e^{a^3/12}}{4k\sqrt{a \pi }}\, , \end{aligned}$$
(6)

implying its finite energy.

As it is shown in Ref. [45], the gravitational Airy wavepackets (5) do not accelerate. However, the main maxima lobe of the wave train does show acceleration in the \(\zeta -y\) plane. This is exemplified in Fig. 1, where a contour plot in \(\zeta -y\) plane of function \(|{\mathcal {G}}(\zeta ,y)|=\sqrt{{\mathcal {G}}^*(\zeta ,y){\mathcal {G}}(\zeta ,y)}\) is shown for \(a=1\) and \(a=2\). In red solid lines we display the evolution of the position of the maxima lobe, that moves in a curved trajectory, thus displaying its acceleration in such plane.

Fig. 1
figure 1

Contour plot for \(|{\mathcal {G}}(\zeta ,y)|\), with \({\mathcal {G}}\) given in Eq. (5), in terms of normalized \(k\zeta \) and ky variables, for different values of a. a Plot for \(a=1\). b Plot \(a=2\). For both cases the red solid line shows the dynamical trajectory of the position of local maxima lobes of the wavepacket, which shows acceleration in the \(\zeta -y\) plane. The colors get brighter as the value of \(|{\mathcal {G}}(\zeta ,y)|\) increases

3 Exact nondiffracting gravitational waves

Structured wave-packets are not only solutions to the linearized Einstein equations, but they are also exact solutions to the full Einstein equations. Exact nondiffracting gravitational waves can be obtained from a modification of Ehlers–Kundt solution [36, 37]. Consider the exact metric in cartesian coordinates given by the interval

$$\begin{aligned} ds^2=-dt^2+ L(u)^2dx^2+W(u)^2dy^2+dz^2\, , \end{aligned}$$
(7)

where L and W are arbitrary functions of \(u=z-t\). Einstein equations are satisfied when

$$\begin{aligned} \frac{1}{L}\frac{d^2 L}{du^2}+\frac{1}{W}\frac{d^2 W}{du^2}=0\, . \end{aligned}$$
(8)

These waves, described by metric (7), are pp-waves [37, 47]. This can be seen by applying the coordinate transformation \(x=X/L\), \(y=Y/W\), \(v=V-X^2 d_u L/L-Y^2 d_uW/W\), together with Eq. (8), to put the metric (7) in the Kerr–Schild form \(ds^2=dX^2+dY^2-du\, dV+du^2(X^2-Y^2)d_u^2 L/L\). For these classes of metrics linear and non-linear theories coincide under certain conditions [48].

Here we show that there are certain solutions of Eq. (8) that describe exact nondiffracting and non-accelerating gravitational waves. The simplest non-trivial solution for a structured gravitational wave can be obtained in terms of Airy functions

$$\begin{aligned} L(u)= & {} L_0\, {\text{ Ai }}(u)\, ,\nonumber \\ W(u)= & {} W_0\, {\text{ Ai }}(-u)\, , \end{aligned}$$
(9)

with arbitrary constants \(L_0\) and \(W_0\). These nondiffracting wave-packets have the property that the intensity of the gravitational wave is concentrated in the caustic. However, these exact Airy gravitational waves have infinite energy. This is a feature that they share with other gravitational waves, such plane waves [37]. Nevertheless, what is relevant in this exact solution is its structured configuration.

4 Discussion

Nondiffracting wave-packets have been a fruitful field of study in quantum mechanics [1, 2, 49, 50], electromagnetism and optics [3,4,5,6,7,8,9,10,11]. The solutions (4) and (9) presented along this work are in the same spirit, and complement others solutions for structured gravitational beams, as Bessel or Laguerre–Gauss, found in linearized general relativity [51, 52].

Although simple in their mathematical form, the new linearized accelerating and exact non-accelerating solutions can bring new insights in the form that different gravitational waves propagate in a flat space-time background.

How the nondiffracting features are manifested for gravitational waves in curved space-time is an open question currently under investigation.