## Abstract

Research on time-varying media has recently enjoyed renewed interest, especially in photonics. Despite the large amount of research done in this field in the past few years, the attention has been focused on electromagnetic waves solely, while a comprehensive framework describing how wave phenomena in general are influenced by time-varying media has not been fully developed yet. To this aim, we study the implications of time-varying wave mechanics, and show how the standard wave equation is modified if the speed of a wave is not constant in time. In particular, waves that experience longitudinal acceleration are shown to have clear relativistic properties when a constant reference speed exists. Moreover, the accelerating wave equation admits only solutions propagating forward in time, which are continuous across material interfaces. We then consider the special case of electromagnetic waves, finding that the Abraham–Minkowski controversy is caused by relativistic effects, and the momentum of light is in fact conserved between different media. Furthermore, we show that the accelerating waves conserve energy when the wave is moving along a geodesic and demonstrate two example solutions. We conclude with some remarks on the role of the accelerating wave equation in the context of the arrow of time.

Published by Optica Publishing Group under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.

## 1. INTRODUCTION

Time-varying media are unconventional types of matter, where the properties of the medium change uniformly in space, as a function of time. These types of materials have been studied at least since the 1950s in the context of electromagnetism and plasma science [1–3], as well as in the case of quantum physics [4,5]. Time-varying media have recently enjoyed renewed interest in wave mechanics, in particular for the case of electromagnetic radiation [6]. This is due to the emergence of highly nonlinear materials, such as epsilon-near-zero (ENZ) media, which may feature large and fast modulation to the dielectric function under specific conditions. Time-varying media have been shown to host a range of interesting effects, such as temporal scattering [1], time crystals [7–9], photonic time crystals [10,11], and antireflection temporal coatings [12], to name a few. What has been neglected so far is that when the properties of matter change, the wave propagating inside it experiences longitudinal acceleration. This is true no matter what type of wave we are discussing.

The speed of light may be the most important constant of nature. It governs the physics of the small and large, sets boundaries in the theory of relativity [13], and is directly linked to a large number of observables [14–16]. In the framework of modern physics, it is normal practice to take the speed of light as a strictly constant value [17]. However, we encounter situations where it varies every day: whenever light interacts with matter, it appears to slow down. Obviously, this is not a modification to the vacuum speed of light but it is an important property nonetheless.

However, this presents a problem: although electromagnetic waves clearly experience longitudinal acceleration, such accelerating waves are *not* solutions of the standard wave equation, although they are solutions to Maxwell’s equations. Motivated by this observation, it is clear that the same problem persists in all wave mechanics, whether electromagnetic or mechanical. In the present work, we examine the implications of a varying speed in wave propagation, and derive a general wave equation for this case. The considerations here apply to propagation through free space as well as any media. We start by deriving a wave equation for arbitrary accelerating waves, and present a general solution to it. The solutions have a preferred direction of time, therefore establishing the microscopic arrow of time. We then continue by concentrating on the special case of electromagnetic waves, and show that an electromagnetic wave that has a varying speed will experience relativistic effects. Moreover, when the wave is moving along a geodesic, its energy and momentum are conserved. This has important implications on one of the longest standing debates in photonics. We then consider a curved spacetime geometry, which causes the wave to gain energy upon propagation.

## 2. ACCELERATING WAVE EQUATION

For smooth functions, differentiation is commutative, i.e., the result of iterated derivation does not depend on the order in which the single derivatives are applied. Therefore, multiplication with first order differential operators on a set of smooth functions makes perfect sense in a quite general way [18]. For the sake of simplicity, in what follows we limit ourselves to consider differential operators on twice differentiable functions, in which case higher order derivatives may not be defined, but the multiplication with differential operators is still valid.

#### A. Derivation and a General Solution

Let us consider ${1 + 1}$-dimensional motion for simplicity, with the understanding that the following calculations can be generalized to three spatial coordinates as shown in Appendix A via differential geometry. Let us say that we have a function $f$ that exists on a spacetime grid $(x,t)$ such that $f \equiv f(x,t)$. With the use of the triple product rule, we can write [19]

Since we are considering motion, the third partial derivative can be interpreted as the inverse speed of the function $f$ on the grid, that is,

which we may plug into Eq. (1) without loss of generality.Furthermore, if we assume that the function $g = \partial x/\partial f$ is a smooth enough function on a given domain $D = \{g \in {C^1}:g(f) \ne 0\}$, then it is well behaved, and in particular, its inverse, ${g^{- 1}} = \partial f/\partial x$, exists and is well behaved as well. We may then multiply both sides of Eq. (1) with $c{g^{- 1}}$, and collect everything on the l.h.s. to obtain

where we have left the vertical evaluation lines implicit for brevity. Furthermore, if we multiply both sides with $\partial /\partial t - c\partial /\partial x$, and assume that $c$ is constant, we arrive at the familiar wave equationThe reader should keep in mind that so far, the derivation presented above does not offer any information on the nature of the field $f$, and such information has to be inferred from somewhere else (e.g., from Maxwell’s equations). Remarkably, the only assumption that is required to arrive at the standard wave equation is that the speed of the wave is constant everywhere and at all times.

Let us now consider the case in which we allow the speed of the wave to vary with time, i.e., we set $c \equiv c(t)$. In this case, since the speed of the wave is not constant anymore, when we multiply Eq. (3) by $(\partial /\partial t - c\partial /\partial x)$ on the left, we will get an extra term, proportional to the temporal derivative of the velocity. In this case, the above treatment then yields a wave equation for an accelerating wave:

Let us now analyze Eq. (5). If we take the spatial Fourier transform of the wave, defined as $\tilde f(k,t) = \int_{- \infty}^\infty f(x,t)\exp (ikx){\rm d}x$, where $k = \omega /{c_0}$ is the wavenumber, $\omega$ is the angular frequency, and ${c_0}$ is some constant reference speed, we attain an equation of the form

where ${\Omega ^2}(t) = {k^2}c{(t)^2} - ik\dot c(t)$. This equation is analogous to a harmonic oscillator with a time dependent frequency, and analyzing it usually requires the use of invariants or other sophisticated mathematical tools [21–24]. However, in the present case where $\Omega (t)$ depends on the speed and acceleration of the wave, there is a simple plane wave solution of the form where the minus sign is not a solution. Note that the integral is indefinite, which introduces a constant phase shift to the solution. However, this is not detectable, and can therefore be ignored. Alternatively, one can also write the integral as definite and introduce upper and lower bounds as necessary.This solution was obtained with a trial method. Simply by inserting a plane wave with laboratory time into Eq. (6), one can see from the remaining terms that the true solution must contain an integral over the speed of the wave, thus motivating a trial with Eq. (7). Inserting this into Eq. (6) reveals that the sign must be positive, or otherwise there are terms that do not cancel out. As a side note, it should be possible to find good approximate solutions to Eq. (5) with standard perturbation theory.

Intriguingly, we can now draw direct parallels to the theory of relativity by choosing the reference speed, ${c_0}$, as the vacuum speed of light, and writing the speed variations as a modulation to this speed such that $c(t) = {c_0}n(t)$, where $n(t) \in [0, 1]$. With this choice, we immediately identify the integral in Eq. (7) to be analogous to the *proper time* found in general relativity [17], namely, ${t^\prime} = \int \sqrt {{g_{\textit{tt}}}} {\rm d}t$, where ${g_{\textit{tt}}}$ is the time component of the metric ${g_{\mu \nu}}$ characterizing the curved spacetime in which the time experienced by the wave is calculated. Therefore, $n(t) \equiv \sqrt {{g_{\textit{tt}}}}$, and we can identify the time $t$ as the *coordinate time* measured by a stationary observer. The ratio between these two times is what defines *time dilation*.

Going back to the spatial domain, we can write the general solution as a superposition of waves moving towards the positive and negative directions along the space axis:

where the amplitudes $A$ and $B$ depend on the particular solution. Let us briefly consider the spatial coordinates in the different frames with a constant $n(t) = n$ for simplicity. In general, a constant $n$ defines a geodesic in a Minkowski-type spacetime with a mostly positive signature [13]. In the stationary frame, we have $dx/dt = c(t)$ such that the leading edge of the wave is at $x = \int c(t){\rm d}t = {c_0}{t^\prime}$, where ${t^\prime} = tn$. If we multiply both sides by $k$, we find that $kx = \omega {t^\prime}$. For a stationary observer who describes the dynamics with respect to the laboratory time, the relation $k{x^\prime} = \omega t$ must hold, which leads to ${x^\prime} = x/n$. In other words, the wave will undergo a*length contraction*in the stationary frame. The above reasoning shows that not only matter, but also waves experience relativistic effects. In particular, a wave that is moving at a significant fraction of the speed of light will experience special relativistic effects, whereas a longitudinally accelerating wave experiences general relativistic effects, reminiscent of gravitation. It needs to be stressed that these considerations apply to all waves, whether they are electromagnetic or mechanical in nature.

#### B. Arrow of Time

Let us briefly review our understanding of time. In everyday experience, time always seems to have a preferred direction. That is, time always moves towards the future and away from the past. However, there is no concrete reason why this should be. As Sir Arthur Eddington wrote [25]: “Let us draw an arrow arbitrarily. If as we follow the arrow we find more and more of the random element in the state of the world, then the arrow is pointed towards the future; if the random element decreases, the arrow points towards the past. That is the only distinction known to physics.” This was how Eddington coined the term “arrow of time,” and tied it explicitly to the second law of thermodynamics. Up to now, this has been the best explanation for the existence of a preferred direction of time. Nonetheless, it is only an observation of the behavior of macroscopic objects, and not an actual law of nature.

To see why this is, consider a collection of repulsive particles in a highly ordered configuration and inside a bounded region, which then evolves towards *either* positive *or* negative time. This situation is depicted schematically in Fig. 1. In both cases, the total entropy will increase. However, if time were to be reversed at any point during the evolution (dashed lines in Fig. 1), entropy would decrease until the particles go back to their initial highly ordered configuration, and then start rising again. Therefore, increasing entropy is a sign of a macroscopic arrow of time, and it does not tell us anything about the direction of microscopic time (Ref. [25], p. 79), or why time cannot be reversed (Ref. [25], p. 93–94).

Here, we shed light on these questions. An observation of great importance is that when a wave experiences longitudinal acceleration, Eq. (5) allows only solutions that propagate **towards positive time**, as seen from Eq. (7). Moreover, reflection along the temporal axis (i.e., true time reversal) is allowed only when $\dot c(t) = 0$, but any reflection process must induce an acceleration. Thus, a wave must experience a positive time, which cannot be reversed. Note that this is distinct from time reflection discussed, e.g., in Refs. [26–31] in the context of electromagnetic time dependent media. This ambiguity in terminology requires some consideration.

In the context of electromagnetic waves, time dependent media have been studied at least since the 1950s [1,6], but the term time refraction has been introduced only relatively recently [26,27]. Multiple studies on spatiotemporal effects can be found in existing literature, where the term time reflection or time refraction is often used explicitly [26–29]. However, this usually does not mean a reflection or refraction along the time dimension, but rather, a spatial reflection $x \to - x$, due to a time dependent change in the material properties. For example, in Refs. [27,28], the case of an abrupt and uniform change of the physical properties of a bulk material were investigated, and the properties of time reflection and refraction were discussed in detail.

This is the most common use of the term “time reflection,” and it is entirely valid in the form it has been introduced in the literature. In the present study, we discuss reflection as a mathematical process *along* the temporal dimension with respect to the origin of time, i.e., $t \to - t$, while leaving the spatial coordinate unchanged. This constitutes a *true time reversal* for an arbitrary wave, which is an entirely different concept. Indeed, if the time of a wave could be exactly reversed, it would propagate in a manner that allows it to reverse its direction and even any diffraction effects.

Further, it needs to be noted that studies on time reversal have been successfully carried out with optical beams [30]. This has been realized with Pendry’s four-wave mixing method [31], where a transition of ${-}2\omega$ is induced on a probe beam of frequency $\omega$ to cause a shift to ${-}\omega$, i.e., a reflection along the frequency axis $\omega \to - \omega$. Since a plane wave has a time dependence of the form $\exp (i\omega t)$, this has a similar effect as reversing time. Our theoretical results are not contradicting this in any way, since even in this type of experiment, time still flows forward.

The time asymmetry of the solutions we report is in fact the first theoretical result on the *microscopic* arrow of time. In other words, we have found a plausible explanation as to why time has a preferred direction, at least in terms of wave mechanics. Moreover, we show that there is a plausible theoretical reason why time cannot be reversed, and thus entropy must increase in all cases, thus constituting an important addendum to the macroscopic arrow of time. Since the accelerating wave equation is universal, this appears to be widely applicable.

## 3. ACCELERATING ELECTROMAGNETIC WAVES

To gain further insight into the accelerating wave equation and its solutions, we need to choose the type of wave to be studied. Electromagnetic fields constitute an important special case of longitudinally accelerating waves, since they travel at the vacuum speed of light by default, and we can only decelerate them. This leads to some unintuitive properties, which will be discussed at length. In the case of electromagnetic radiation, $c(t)$ is the speed of light in the above equations, and $f$ is the electric (or magnetic) field. For electromagnetic waves, the modulation to the vacuum speed is more naturally written as $c(t) = {c_0}/n(t)$, so that $n(t) \in [1,\infty]$, and it can be interpreted as a refractive index. In principle, lower values of $n$ are also possible, but such cases require dispersion for them to be physical.

Now, we can define an *intrinsic* time for the electric field $f(x,t)$, as ${t^\prime} = \int n{(t)^{- 1}}{\rm d}t$, which is a time experienced by the electromagnetic wave. The notion of time in the frame of an electromagnetic wave may seem faulty, since photons move along null geodesics. Still, intrinsic time is well defined, since electromagnetic waves accumulate a phase on propagation, which depends on time. Notice how in general a constant $n(t) = n$ defines a geodesic in a Minkowski-type spacetime with a mostly positive signature [13], but is characterized by an effective speed of light $c = {c_0}/n$. The usual Minkowski spacetime with $c = {c_0}$ is then obtained for the particular value of $n = 1$, i.e., when the wave is propagating in vacuum. In other words, the acceleration $\dot c(t)$ effectively induces a nonstationary metric ${g_{\mu \nu}}$ in the moving reference frame of the electric field. In such a frame, the dynamics of the field $f(x,t)$ are that of a field evolving on a curved background. However, an external observer will measure a time $t$ and describe the dynamics of $f(x,t)$ in a flat spacetime, as usual.

If we again consider the leading edge of the wave at $x = {c_0}{t^\prime}$, where the intrinsic time is ${t^\prime} = t/n$, we find that the electromagnetic wave experiences *time contraction* instead of time dilation. This follows from the fact that we are decelerating the wave from its usual speed, as opposed to the acceleration of mechanical waves considered above. Moving on in the same manner as in the mechanical case, we multiply both sides by $k$ to get $kx = \omega {t^\prime}$, and for a stationary observer who describes the dynamics with respect to the laboratory time, we get ${x^\prime} = xn$ so that the wave undergoes *length dilation*. Again, this is entirely due to the fact that we are causing the wave to slow down, and a positive acceleration would cause the usual length contraction. However, this effect is usually absorbed into the wavenumber, and thus we can define ${k^\prime} = kn$. This is exactly the relation that is employed in conventional optics and photonics.

Note that when the speed of the wave from an outside observer’s point of view is equal to the vacuum speed, the intrinsic time corresponds to the laboratory time, and the acceleration term in Eq. (5) is null, so that the above reasoning reduces to the usual wave equation and a plane wave solution.

#### A. Energy and Momentum Conservation

Next, we briefly consider the energy and momentum carried by a longitudinally accelerating wave, as seen by a stationary observer. Conservation laws are most elegantly found via Noether’s theorem. However, in this work, in order to simplify the treatment, we do not employ an action formalism, and we therefore examine the conservation of energy in a more roundabout way. We can find the energy flux with the Poynting vector ${\textbf S}$, and the average energy density simply by dividing $S = |{\textbf S}|$ with the speed of the wave. Furthermore, the overall energy $E$ can be attained by integrating the average energy density over the whole volume occupied by the wave, which in this case would be of the form

where $w$ is a constant that comes from integrating along the transverse dimensions, and the coordinates $({x_0},{x_1})$ define the length of the volume. The two spatial coordinates can be expressed in terms of intrinsic time as ${x_0} = {c_0}t_0^\prime $, and ${x_1} = {c_0}t_1^\prime $. Therefore, performing the integral yields $E = wSn(t)\Delta {t^\prime}$, where $\Delta {t^\prime} = t_1^\prime - t_0^\prime $ can be chosen to be a constant. If $n(t)$ defines a geodesic, the above relation yields $E = wS\Delta t$, for any constant value of $n$, and thus, energy is conserved. However, this is generally not the case in a curved spacetime, which is a well-known result in general relativity [32]. We will examine one such case later on, in Section 3.C.Since it is well known that the momentum of an electromagnetic wave is found by dividing its energy with its phase velocity, conservation of energy immediately implies conservation of momentum. We can see this also from our plane wave solution: no matter what type of longitudinal acceleration the wave experiences, the wavenumber $k$ remains constant in the moving reference frame of the wave. Since the wavenumber is directly related to the momentum the wave carries, we can conclude that the momentum of a longitudinally accelerating wave is a conserved quantity. More precisely, it is a Lorentz covariant scalar, and thus remains the same under Lorentz transformations.

#### B. Abraham–Minkowski Controversy

Momentum conservation has a direct consequence on one of the longest standing fundamental questions in photonics: what is the momentum of a photon that enters a dielectric medium? Historically, there have been two contending views, and as such it has been named the Abraham–Minkowski controversy [33–37].

Let us start by applying our plane wave solution to an electromagnetic wave crossing a dielectric interface. Here, the speed of light abruptly changes from some initial value, ${c_0}$, to something else at time $t = 0$, and we can write

where $m$ is a free parameter. This is essentially a smooth approximation of the Heaviside step function, when $m$ is sufficiently large. Here, ${n_1}$ is a proportionality factor by which the speed of light changes after $t = 0$, i.e., $c(t \gt 0) = {c_0}/{n_1}$, and we take ${n_1}$ to be equal to the refractive index of the material. Integrating with respect to $t $, we obtain the intrinsic time asThis ensures that before $t = 0$, the time of the wave will flow as in free space (${t^\prime} = t$), whereas after entering the medium, intrinsic time becomes ${t^\prime} = t/{n_1}$. A depiction of the above scenario is shown in Fig. 2, for $\omega = 20$, ${n_1} = 3$, ${c_0} = 1$, and amplitude coefficients are calculated with Fresnel equations. The reflected wave does not appear in the temporal reference frame for clarity.

As can be seen from Fig. 2, our approach with a smooth, continuous electric field crossing a dielectric boundary is equivalent to the standard formulation of a piecewise continuous electric field defined via electromagnetic boundary conditions (at least in this simple case). In fact, it would appear that we can reformulate many of the problems in optics and photonics in terms of accelerating waves and do away with dielectric boundaries. Moreover, the interpretation of the wave properties becomes more straightforward within the present framework; if we insist on the intrinsic time description, both energy and momentum are conserved since $\omega$ and $k$ are both constant. If we instead say that time should always be taken in the outside observer’s reference frame, then the variation can be transferred to the frequency, such that ${\omega ^\prime} = \omega /{n_1}$. If we insist on this explanation, then we are forced to choose a new wavelength inside the medium to conserve energy. However, if we do that, then momentum is no longer conserved, leading to the famous Abraham–Minkowski controversy.

Here, we argue that the whole controversy stems from relativistic effects, and is in fact, not a problem at all when inspected in the correct frame of reference. Indeed, one can make compelling arguments for both ${n_1}\lambda$ and $\lambda /{n_1}$ cases that are equally correct. However, neither reaches the root of the problem, which is that the apparent change in the wavelength (and thus momentum) is a relativistic effect, and the momentum of the wave is in fact unchanged. In other words, one can find a *proper wavelength*, which is invariant in this case. In usual interactions such as crossing the boundary of an ordinary linear medium, only the time experienced by the wave changes when we are in the moving reference frame of the wave. This in turn means that the relativistic momentum of the wave remains unchanged. Our findings are supported by the heuristic arguments recently laid out in Ref. [38]. In contrast to earlier resolutions to the controversy, the present formalism does not require a division between material and optical momenta, greatly simplifying the physical picture.

#### C. Time-Varying Media

As a second example of the present study, we look at disordered photonic time crystals [10,11], which is a prime example of a time-varying medium. We examine this special case, since the physical situation is governed by the accelerating wave equation. A disordered photonic time crystal is defined as a medium where the dielectric function changes abruptly and randomly as a function of time. So far, the photonic time crystal has been studied with the use of discrete time modulations to the dielectric constant. In Ref. [10], the authors demonstrated via numerical simulations that when the modulation to the dielectric function is fast and strong enough, then a pulse propagating inside a disordered photonic time crystal features an exponentially decreasing speed, while the overall energy of the pulse increases exponentially. This was further corroborated in Ref. [11], and this is also the case we shall consider next.

As is well known, the microscopic explanation for the slowing down of an electromagnetic wave’s phase velocity inside a material is due to the disturbance in the charges of each atom. As the electromagnetic field oscillates inside a material, the charges will be accelerated back and forth at the same frequency. This acceleration causes the charges to radiate their own electromagnetic wave at the same frequency, but usually with a phase delay, as the charges may move out of phase with the driving wave. The total electromagnetic wave propagating inside the medium is the superposition of all the waves inside the medium, and since some are out of phase with the driving wave, it will cause the total wave to slow down. This is directly quantified with the refractive index.

Similarly, in a photonic time crystal each time the wave experiences a change in the dielectric function, it will be reflected and refracted. If the modulation of the dielectric function is fast enough, the waves will overlap. The total wave is then again the superposition of all waves present in the material. Continuing to temporally modulate the dielectric function quickly enough will cause an ever increasing number of waves (that are generally out of phase) to overlap inside the material. This will then cause the total wave to slow down, which is accompanied by an increase of amplitude. The effect of a disordered photonic time crystal can therefore be thought to result from the increase of an effective refractive index.

In our study, we do not consider the underlying mechanism that causes the wave to slow down exponentially. Instead, we simply suppose that there is some material where this occurs, and then show how the solutions to the accelerating wave equation behave in this situation. Our analytical solutions would allow us to model the case of discrete modulation as well, which would just be a superposition of a large number of waves similar to the one presented in Fig. 2. While applicable to the discrete case, our results can be thought of as the continuous limit where the length of each modulation segment, $L$, is taken to be vanishingly small, i.e., $L \to 0$. Indeed, our approach allows for continuous solutions to problems where discontinuous (or piecewise continuous) methods have been used before, such as when light is incident at the interface of two media.

Let us then take $c(t)$ to be of the exponentially decaying form such that $c(t) = {c_0}\exp (- t/\tau)$, where $\tau$ is the decay half-life. For this particular case, we can write a plane wave solution in the form

Further employing Eq. (9), we find that $E \propto \exp (t/\tau)$. In other words, the curved spacetime will cause the wave to violate energy conservation, as expected. However, it must be noted that the violation of energy conservation in a time-varying medium is a local effect, which is maintained by pumping the medium somehow. Energy is still conserved in the global sense.

Since the accelerating wave equation is linear, we can write a superposition integral in the form

As can be seen from Fig. 3, the pulse slows down, increases in energy, and becomes shorter. Apart from the shortening of the pulse, these results correspond to earlier works. Therefore, it would seem that disordered photonic time crystals are potential analogue systems for studying wave dynamics in a curved spacetime. This is particularly interesting from an experimental point of view, since it has been proposed that such crystals can be created in a laboratory environment, by illuminating an ENZ medium with ultra-short pulses [10].

Although the present work is entirely theoretical, it is of great interest to consider how such experiments could be realistically carried out. Obviously, the first example in the main text is entirely trivial, as it corresponds to a simple dielectric interface. However, the second example is not as straightforward.

As discussed in the supplementary material of Ref. [10], it is possible to greatly simplify this problem by considering radio waves. In this regime, the frequency as well as the large temporal period make it relatively simple to temporally modulate a signal such that these effects are seen also experimentally. Such experiments have been carried out by modulating the capacitance of transmission lines. Since radio waves and light are both electromagnetic waves, it is to be expected that results with such a simplified system are directly applicable to optical frequencies. However, the experimental implementation is necessarily more complicated.

Again, as outlined in the supplementary material of Ref. [10], possible experiments in the optical regime are based on nonlinear epsilon-near-zero (ENZ) materials. ENZ materials are materials with a real part of the permittivity that is close to zero. Due to their strong nonlinear response, the properties of such materials can be greatly modified by illuminating them with ultrashort optical pulses at frequencies sufficiently higher than their plasma frequency. In such a case, an ultrashort pump pulse is absorbed, exciting a large number of electrons to the conduction band, thus modifying the free charge density. At the ENZ wavelength, this causes changes to the refractive index that can be on the order of unity. Since the response time of these materials is on the order of a few femtoseconds, it is then only a matter of structuring the pump pulses in a suitable manner, such that they induce suitable fluctuations to the refractive index that are necessary for experimentally demonstrating the disordered photonic time crystals.

Although, the imaginary part of the dielectric function may cause considerable losses, such that the ENZ material has to be optically thin for observable effects to appear. Moreover, the ultrashort pulses may cause considerable optical damage to the ENZ material, thus limiting the applicability of this approach. Nonetheless, nonlinear ENZ materials are probably the most promising candidates for experimentally realizing disordered photonic time crystals. Other possible avenues are temporally modulated multilayer structures and metasurfaces [39,40].

## 4. DISCUSSION AND CONCLUSIONS

To summarize, in this work, we have derived a wave equation for longitudinally accelerating waves. The novel accelerating wave equation is applicable to mechanical, as well as electromagnetic waves, and could be of interest to the wider physics community. We then concentrated on electromagnetic waves and, in particular, situations where the speed of light is modulated away from its vacuum value. In the cases we consider, this modulation corresponds to either a simple change in refractive index, or a disordered photonic time crystal where the speed of light decreases exponentially. We do not consider instances where the speed of light is greater than the vacuum value. Within the analysis, we find several important features of accelerating waves.

First, the solutions to the accelerating wave equation are continuous everywhere. This is in stark contrast to the usual formulation of photonics, where material interfaces have to be treated via continuity conditions. Second, such waves experience relativistic effects, with a clear connection to the Abraham–Minkowski controversy: if there is a constant reference speed ${c_0}$, the energy and momentum of the wave are conserved when the wave is moving along a geodesic (even across material interfaces). Third, the exponential increase in the energy of a pulse propagating in a disordered photonic time crystal is due to a curved spacetime that the pulse experiences, and may offer experimental tests of general relativity.

We consider these to be literal relativistic effects in the sense that they are observer dependent, and one can measure them (at least in principle). In fact, many of them have been measured countless times. The general relativistic effects we discuss can also be thought of as a mathematical analogy. This is because general relativity is usually considered to be a theory of gravity, and we do not consider any effects caused by massive objects. This analogy is readily understood if one thinks of the refractive index $n(t)$ as a mass acquired by a photon when entering a medium. A time dependent refractive index induces an effective curvature on the spacetime in the reference frame attached to the wave. This is fully analogous to a mass that curves spacetime, which generates a gravitational field.

Last, the accelerating wave equation allows only solutions that propagate towards *positive* time. This is a very peculiar asymmetry, which has been encountered before only in the context of the weak force [41]. This may have important implications for the arrow of time: since the accelerating wave equation is general, all light should obey it. This may even include force carrier photons, which are not described by any known wave equation. If this is indeed the case, then forces that are mediated via virtual photons *cannot* be time reversed, and they have a preferred direction of time. It needs to be emphasized that this is a novel result, distinct from the thermodynamic arrow of time.

There are experimental tests that can be carried out to confirm some of the findings presented here, most notably in disordered photonic time crystals. By engineering the curved spacetime properly, it may be possible to generate short and intense pulses. Moreover, it would be of interest to see how light behaves in the present formalism as it approaches the event horizon of a black hole. Finally, we will examine the properties of the ${3 + 1}$-dimensional accelerating wave equation in detail in a separate work.

## APPENDIX A: ALTERNATIVE DERIVATION OF Eq. (5)

The aim of this appendix is to provide a more general and elegant derivation of Eq. (5) that complements and enriches the derivation given in the main text based on the triple product rule [Eq. (1)]. Although there are multiple different mathematical methods available to derive the accelerating wave equation, based on either a Green’s function [42] or integral [43] approach, the derivation we give here is instead based on differential geometry and as such is reference frame and coordinate independent. This—as will be clear at the end of this appendix—allows for an easy generalization of our results to an arbitrary number of spatial dimensions.

**Exterior Derivative, Hodge Dual, and the Wave Equation.** Before presenting our alternative derivation, it is instructive to recall some concepts of differential geometry, in particular those of exterior derivative, exterior product, and Hodge dual. We follow the notation convention from Nakahara [18], which we also refer the interested reader to for the details about the rigorous definitions and derivations of these quantities. Here, we only give an operational definition.

In this paper, we are interested in differential forms over a two-dimensional pseudo-Riemannian manifold, i.e., the two-dimensional Minkowski space ${\mathbb{M}^2}$. To fix the notation, and write the exterior derivative in coordinate form, we choose a local coordinate frame on ${\mathbb{M}^2}$, i.e., $\{x,\tau \}$, such that the metric tensor will be given as ${d^2}s = {g_{\mu \nu}}d{x^\mu}d{x^\nu} = d{x^2} - d{\tau ^2}$. Let us consider a $0 -$ form, i.e., a function, $f \in {\Omega ^0}({\mathbb{M}^2})$ and represent it in the local coordinate frame as $f(x,\tau)$. The exterior derivative $df$ of the 0–form $f$ is a 1–form, and can be then thought of as the generalization of the concept of the differential of a function to the case of differential forms, and its explicit expression in the local coordinate frame is given as

If the exterior derivative would instead act on a 1–form, the result would be a 2–form. For example, if we consider the 1–form $g(x,\tau)dx$, the exterior derivative $dg$ would be given by

We now want to use Eq. (A4) to derive the accelerating wave equation.

## 1. Derivation of Eq. (5)

Following Eq. (A4), the exterior derivative of the 0-form $f(x,\tau)$ is given by

The action of the Hodge dual on $df$ simply exchanges the $1$–forms $dx$ and $d\tau$, giving

Now, the second exterior derivative acts on the $1$–form $df$, thus generating the $2$–form

The equation above reduces to the usual wave equation in $1 + 1$–dimensions if we identify $\tau = ct$, with $c = {\rm constant}$, since

If, on the other hand, one assumes that the speed of the wave is not constant, i.e., $c = c(t)$, we get instead

Inverting this, we can rewrite the $\dot c$-dependent term above in terms of a spatial derivative, thus obtaining

## 2. Generalization to Higher Dimensions

Unlike the triple product method presented in the main text, the method above has the advantage that it can be readily generalized to an arbitrary number of spatial dimensions. If we now assume $f \in {\Omega ^0}({\mathbb{M}^{n + 1}})$, and endow the $(n + 1) -$ dimensional Minkowski space ${\mathbb{M}^{n + 1}}$ with the metric ${d^2}s = - d{\tau ^2} + \def\LDeqbreak{} dx_1^2 + dx_2^2 + \cdots + dx_n^2$ (where, again, $\tau = ct$) we can, following the steps highlighted above, derive the accelerated wave equation in higher dimensions as follows:

## Funding

Academy of Finland (333938, 320166, 320165).

## Disclosures

The authors declare no conflicts of interest.

## Data availability

No data were generated or analyzed in the presented research.

## REFERENCES

**1. **F. R. Morgenthaler, “Velocity modulation of electromagnetic waves,” IRE Trans. Microw. Theory Tech. **6**, 167–172 (1958). [CrossRef]

**2. **C. Elachi, “Waves in active and passive periodic structures: a review,” Proc. IEEE **64**, 1666–1698 (1976). [CrossRef]

**3. **D. K. Kalluri, V. R. Goteti, and A. M. Sessler, “WKB solution for wave propagation in a time-varying magnetoplasma medium: longitudinal propagation,” IEEE Trans. Plasma Sci. **21**, 70–76 (1993). [CrossRef]

**4. **J. H. Shirley, “Solution of the Schrödinger equation with a Hamiltonian periodic in time,” Phys. Rev. **138**, B979–B987 (1965). [CrossRef]

**5. **C. K. Law, “Effective Hamiltonian for the radiation in a cavity with a moving mirror and a time-varying dielectric medium,” Phys. Rev. A **49**, 433–437 (1994). [CrossRef]

**6. **E. Galiffi, R. Tirole, S. Yin, H. Li, S. Vezzoli, P. A. Huidobro, M. G. Silveirinha, R. Sapienza, A. Alù, and J. B. Pendry, “Photonics of time-varying media,” Adv. Photon. **4**, 014002 (2022). [CrossRef]

**7. **K. Sacha, “Modeling spontaneous breaking of time-translation symmetry,” Phys. Rev. A **91**, 033617 (2015). [CrossRef]

**8. **D. V. Else, B. Bauer, and C. Nayak, “Floquet time crystals,” Phys. Rev. Lett. **117**, 090402 (2016). [CrossRef]

**9. **N. Y. Yao, A. C. Potter, I.-D. Potirniche, and A. Vishwanath, “Discrete time crystals: rigidity, criticality, and realizations,” Phys. Rev. Lett. **118**, 030401 (2017). [CrossRef]

**10. **Y. Sharabi, E. Lustig, and M. Segev, “Disordered photonic time crystals,” Phys. Rev. Lett. **126**, 163902 (2021). [CrossRef]

**11. **R. Carminati, H. Chen, R. Pierrat, and B. Shapiro, “Universal statistics of waves in a random time-varying medium,” Phys. Rev. Lett. **127**, 094101 (2021). [CrossRef]

**12. **V. Pacheco-Peña and N. Engheta, “Antireflection temporal coatings,” Optica **7**, 323–331 (2020). [CrossRef]

**13. **A. O. Barut, *Electrodynamics and Classical Theory of Fields and Particles* (Dover, 1981).

**14. **K. Huang, *Fundamental Forces of Nature—The Story of Gauge Fields* (World Scientific, 2007).

**15. **J. D. Jackson, *Classical Electrodynamics*, 3rd ed. (Wiley, 1998).

**16. **K. S. Mendelson, “The Story of c,” Am. J. Phys. **74**, 995 (2006). [CrossRef]

**17. **L. D. Landau and E. M. Lifshitz, *The Classical Theory of Fields* (Butterworth-Heinemann, 1975).

**18. **M. Nakahara, *Geometry, Topology and Physics*, 2nd ed. (Taylor & Francis, 2003).

**19. **A. H. Carter, *Classical and Statistical Thermodynamics* (Prentice Hall, 2001).

**20. **D. A. Angus, “The one-way wave equation: a full-waveform tool for modeling seismic body wave phenomena,” Surv. Geophys. **35**, 359–393 (2014). [CrossRef]

**21. **H. R. Lewis, “Classical and quantum systems with time-dependent harmonic-oscillator-type hamiltonians,” Phys. Rev. Lett. **18**, 510 (1967). [CrossRef]

**22. **D. C. Khandekar and S. V. Lawande, “Exact propagator for a time–dependent harmonic oscillator with and without a singular perturbation,” J. Math. Phys. **16**, 384–388 (1975). [CrossRef]

**23. **C. J. Eliezer and A. Gray, “A note on the time-dependent harmonic oscillator,” SIAM J. Appl. Math. **30**, 463–468 (1976). [CrossRef]

**24. **C. W. Robson, Y. Tamashevich, T. T. Rantala, and M. Ornigotti, “Path integrals: from quantum mechanics to photonics,” APL Photon. **6**, 071103 (2021). [CrossRef]

**25. **A. Eddington, *The Nature of the Physical World* (Cambridge University, 1929).

**26. **J. T. Mendonça and P. K. Shukla, “Time refraction and time reflection: two basic concepts,” Phys. Scr. **65**, 160 (2002). [CrossRef]

**27. **Y. Xiao, D. N. Maywar, and G. P. Agrawal, “Reflection and transmission of electromagnetic waves at a temporal boundary,” Opt. Lett. **39**, 574–577 (2014). [CrossRef]

**28. **B. W. Plansinis, W. R. Donaldson, and G. P. Agrawal, “What is the temporal analog of reflection and refraction of optical beams?” Phys. Rev. Lett. **115**, 183901 (2015). [CrossRef]

**29. **G. A. Menendez and B. Maes, “Time reflection and refraction of graphene plasmons at a temporal discontinuity,” Opt. Lett. **42**, 5006–5009 (2017). [CrossRef]

**30. **S. Vezzoli, V. Bruno, C. DeVault, T. Roger, V. M. Shalaev, A. Boltasseva, M. Ferrera, M. Clerici, A. Dubietis, and D. Faccio, “Optical time reversal from time-dependent epsilon-near-zero media,” Phys. Rev. Lett. **120**, 043902 (2018). [CrossRef]

**31. **J. Pendry, “Time reversal and negative refraction,” Science **322**, 71 (2008). [CrossRef]

**32. **C. W. Misner, K. Thorne, J. A. Wheeler, and D. I. Kaiser, *Gravitation* (Princeton University, 2017).

**33. **H. Minkowski and H. Die, “Grundgleichungen für die elektromagnetischen Vorgänge in bewegten Körpern,” Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse **53**, 111 (1908).

**34. **M. Abraham and M. Zur, “Elektrodynamik bewegter Körper,” Rendiconti del Circolo Matematico di Palermo **28**, 1 (1909). [CrossRef]

**35. **S. Barnett, “Resolution of the Abraham–Minkowski dilemma,” Phys. Rev. Lett. **104**, 070401 (2010). [CrossRef]

**36. **D. F. Nelson, “Momentum, pseudomomentum, and wave momentum: toward resolving the Minkowski–Abraham controversy,” Phys. Rev. A **44**, 3985 (1991). [CrossRef]

**37. **K. T. McDonald, “Bibliography of the Abraham-Minkowski debate,” 2017, https://physics.princeton.edu//mcdonald/examples/ambib.pdf.

**38. **G. Feng and J. Huang, “A heuristic resolution of the Abraham–Minkowski controversy,” Eur. Phys. J. Plus **136**, 520 (2021). [CrossRef]

**39. **M. S. Masyukov and A. N. Grebenchukov, “Temporally modulated metamaterial based on a multilayer graphene structure,” Phys. Rev. B **104**, 165308 (2021). [CrossRef]

**40. **X. Wang, M. S. Mirmoosa, V. S. Asadchy, C. Rockstuhl, S. Fan, and S. A. Tretyakov, “Metasurface-based realization of photonic time crystals,” Sci. Adv. **9**, eadg7541 (2023). [CrossRef]

**41. **M. E. Peskin and D. V. Schroeder, *An Introduction to Quantum Field Theory* (CRC Press, 2018).

**42. **F. W. Byron and R. W. Fueller, *Mathematics of Classical and Quantum Physics* (Dover, 1992).

**43. **A. Nerukh and T. Benson, *Non-Stationary Electromagnetics—An Integral Equations Approach*, 2nd ed. (Jenny Stanford Publishing, 2018).

**44. **K. F. Warnick and P. St. J. Russel, “Differential forms and electromagnetic theory,” Prog. Electromagn. Res. **148**, 83 (2014). [CrossRef]