![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-6280706e-b7f2-4db3-8b8b-1004de872d77.jpg)

## Time and space, frequency and wavevector: or, what I talk about when I talk about propagation

## Paul Kinsler ∗

Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2AZ, United Kingdom. (Dated: Tuesday 10 th December, 2019)

The existence of waves with negative frequency is a surprising and perhaps controversial claim which has recently been revisted in optics and for water waves. Here I explain a context within which to understand the meaning of the 'negative frequency' conception, and why it appears in some cases.

## Introduction

Claims or analyses involving negative frequency waves occur in a wide range of fields: quantum mechanics, optics &amp; EM [1, 2], and acoustics and water waves [3]. However, some regard the notion of negative frequencies as unphysical. Here I describe four alternative views of the universe with a view to clarifying why negative frequencies might appear, and whether or not we should be worried: the omniscient , the traditional , the directed , and the causal .

The 'traditional' and the 'directed' pictures, involving propagation along time axis and one selected spatial axis respectively, have already been compared for a range of acoustic wave equations [4]. In particular, the convenience of spatial propagation in optics, and the additional utility of the unidirectional approximation has been widely studied (see e.g. [5] and references therein). Here I am most interested in comparing the traditional &amp; directed views, and restrict notation to Cartesian coordinates for simplicity; although extension to alternate systems is of course possible [6]. It is important to note that in practical terms, utilising spatial propagation amounts to making an approximation, since not only the initial conditions but the ongoing propagation require a knowledge of the future which we rarely have.

It is useful to mention two ideas before we start: first, the standard notion of the casual past of a point in time and space, generally called the past light cone; and second, the computational past of that point. This notion of computational past applies primarily to numerical simulations, but can also be useful from more abstract perspective. It consists of both the set of initial conditions and/or computed data which might have affected the calculated state of that chosen point.

In what follows I will denote the profile of the wave field under discussion as F , but rather than festooning it with a plethora of typographical tics (bars, tildes, primes, and the like) to indicate what arguments may or may not have been transformed, I use its arguments alone to define F 's meaning - whether time t or space r =( x y z , , ) , whether it is a spectra with frequency w or wavevector k = ( kx , ky , kz ) , and so on. Further, I will not consider any specifics as to why a specified wave field might develop a certain character or characteristics - pulses, oscillations, spectral features and the like - just what we can say about the wave given the available knowledge.

∗ Electronic address: Dr.Paul.Kinsler@physics.org

In addition to the four sections addressing each of the omniscient, traditional, directed, and causal pictures in turn, along with a brief conclusion, there are appendices which address some issues relevant in particular to the field of nonlinear optics: the handling of time-response models, and the appearance of negative frequency terms.

## I. OMNISCIENT: ALL TIME AND ALL SPACE

The omniscient position is taken when we claim all possible knowledge about how some wave has and will move through its environment. Here we might express our knowledge of some wave-field F by writing it with full space and time arguments: i.e. F t ( , r ) . Most likely, this view is a result of us having an analytic solution to some specific case of the physical model we are interested in.

Given this fully known wave field solution F t ( , r ) , we are at liberty to Fourier transform in either time or space (or both) as we see fit. In the resulting double spectrum F ( w , k ) , we will quite naturally see not only positive and negative frequency components, but also forward and backward wavevector components along each spatial coordinate axes.

Since the spectrum of some function is closely related to the complex conjugate of that spectrum for negative frequencies, we know for the doubly transformed F ( w , k ) that F ( w , k ) = F ( -w , -k ) . Notably, in this case there is no obvious distinction between a positive frequency disturbance evolving in one direction and a negative frequency disturbance evolving in the opposite direction. We therefore need not be troubled by the presence of negative frequencies, since we can reinterpret them as oppositely directed positive ones - we might, for example, refer to the idea that a positron is only an electron travelling backwards in time 1 .

It is also possible to consider a 'limited omniscience' position, where we claim complete knowledge of the wave behaviour in some defined region of space and time. This might be from some (past) experimental results, or be a numerical or analytic solution valid over that limited extent. In such a case, both time and space can be Fourier transformed, and the remarks above about handling and/or interpreting negative frequencies still hold. However, whilst such after-the-fact

1 Of course the handling of most particles is complicated by their mass, whereas photons, being massless, can be their own antiparticle.

FIG. 1: Temporal propagation of waves, where disturbances (or pulses) evolve either forward or backward in space. At any point in the propagation, we know the spatial behaviour of our wave field at all points in space, as indicated by the pale green vertical lines. The pale green shaded triangles indicate the past light cones (i.e. the causal past) of a wave element (black circles) at selected times along the path of the disturbance. In a temporally propagated numerical simulation which has a maximum wave speed, the causal past matches the computational past. A notional interface has been added to the diagram to show how a reflection would behave.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-28e326da-868f-40e2-85e2-8e509eca9e32.jpg)

analysis is typically very useful, caution should be taken when applying the conclusions to an on-going propagation: in such inescapably dynamical situations, the alternatives discussed in the next three sections are more applicable.

## II. TRADITIONAL: TEMPORAL PROPAGATION

The traditional position is taken when we claim all knowledge about how the wave has moved through its environment up until now. This picture is shown on fig. 1, and assumes the normal preferred direction of propagation along the time axis, which is usually from t = -¥ , through t = 0, and towards t =+ ¥ .

Here we might express our knowledge of the wave-field F by writing it at the current time t 0 with full space arguments: i.e. Ft 0 ( r ) ; and where is is implied that we also knew Ft ( r ) at all past times t &lt; t 0 also. The starting state of our propagation matches the early-time boundary conditions - i.e. what would normally be considered to be the initial conditions. An electromagnetic FDTD simulation [7, 8], along with many other finite element approaches, uses this traditional approach; it is also possible to apply factorization methods to get unidirectional versions of temporally propagated wave equations [9]. In this sense the starting state for some temporally propagated simulation - the 'computational initial conditions' - are the same as the traditional physical initial conditions 2 .

2 However, if we were to decide to propagate backwards in time from t =+ ¥

Given some wave field state Ft ( r ) known over all space r at some specified time t = t 0, we are at liberty to Fourier transform in space, but not time. The spatial transform will give us the wavevector or momentum-like properties of the wave field, in a spectrum Ft 0 ( k ) . This will contain both forward and backward wavevector components along each spatial coordinate axis, and we can calculate it not only at the current time t 0, but also all past times t &lt; t 0.

For simplicity, imagine we are following some trajectory through time and space, while counting interesting events that occur locally. Time is always increasing but we might travel any direction through space, although it is easier to think of the case where our position is fixed. Naturally, our count of interesting events will only ever increase, thus any estimate of the time period T between events would be a positive number; as would any frequency estimation f in events-per-second. However, the events themselves may have different spatial characteristics: notably, we may see objects that pass us travelling left, or perhaps travelling right. The differing character of events, accessible to us through our spatial knowledge Ft ( r ) , allows us to impose a sign (or signs) when adjusting our counting total - perhaps +1 for left-going objects, and -1 for right-going ones.

Having discussed how a counting/timing argument that gives us a strictly positive frequency estimate can be converted into a signed count using spatial information, let us revist the temporally propagated spatial spectrum Ft 0 ( k ) of our notional wave field. As already noted, the spatial spectrum of the wave has both positive and negative wavevector components. Starting with the well known velocity relation for waves v = w / k , we can use this to claim that the positive wavevector part of the spectrum represents forward evolving waves (with v &gt; 0), while the negative wavevector part represents backward evolving waves (with v &lt; 0). However, we have to be careful: even a static field profile will have positive and negative wavevector components, since for real-valued Ft ( ) z , Ft ( -k ) = F t ∗ ( k ) . It is the changing complex phase(s) of Ft ( k ) that represents the shift in Ft ( ) z profile from one position to another.

Note that I have so far only discussed frequency estimates . Strictly speaking, in this picture, any system state exists only at an instant in time, so of itself that state has no frequency content at all. Further, even given our knowledge of past states, we cannot calculate a true frequency dependence, because we do not have the entire wave history to hand: we have the past but not the future. The best we might do in getting an up-to-date idea of the spectrum is apply a Laplace transform 3 over the known past information, converting t into a Laplacespectral s , and get a hybrid spectrum F t ′ 0 ( s , k ) . Nevertheless,

to t = -¥ , as is sometimes done when we have preferred final-time boundary conditions, note that those become the computational initial conditions, and that the computational past of any point in the propagation will be its future light cone. The interested reader might try drawing a suitable counterpart to fig. 1 themselves.

3 Although any one-sided transform that seemed appropriate, if applied over past (known) data, would also be fine.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-bb3282e7-b2fd-4421-888c-19119fead8e6.jpg)

it may be possible to characterise the dominant frequencylike properties of the propagation as a function of the know wavevector; and so be able to use a reference frequency W ( k ) [4] as a basis on which to simplify the propagation - perhaps by assuming it to be unidirectional.

However, if we judge that omitting the future behaviour will have a negligible effect, or that it is irrelevant since we only care about an experiment or computation that has already finished, we can simply take data from the appropriate past time interval, and calculate the frequency spectrum of that. This after-the-fact analysis of past data returns us to a version of the omniscient position, albeit one limited in scope, so that the same view of negative frequencies as there can be applied. This is what is usually done (either implicitly or explicitly), which is of course the reason why so many frequency responses &amp; spectra appear in the literature, textbooks, and on datasheets.

So in this traditional picture, we need not be overly troubled by negative frequencies. Either we are being very rigorous, and so deny that a true frequency spectrum exists at all, or we have a spectrum calculated from known historical data, which has negative frequencies reinterpretable as positive ones, just as in the omniscient picture. As a final note, and as a result of the considerations above, the spectroscopists habit of giving spatial (wavevector, or momentumlike) spectra rather than temporal (frequency, or energy-like) ones makes sense from a theoretical perspective, as well as from a practical experimental one.

## III. DIRECTED: SPATIAL PROPAGATION

The directed position is taken when we claim all knowledge about how the wave the wave has moved in a preferred direction along some chosen path through its environment. This picture is shown on fig. 2, where the path is along the z axis from z = -¥ , through z = 0, and towards z =+ ¥ . Here we might express our knowledge of the wave-field F by writing it at the current position (e.g. z = z 0) with full time and transverse space arguments: i.e. Fz 0 ( t , x y , ) ; and where it is implied that we also knew Fz ( t , x y , ) at all prior positions z &lt; z 0 along that path also. This often seems a very natural thing to do, especially when considering unidirectional beam propagation in waveguides or optical fibres. In particular, this position is made most explicit in the PSSD (pseudospectral spatial domain) method for propagating electromagnetic pulses, but also in many others [5, 10-13]. Here, the starting state at z = zi for some spatially propagated simulation - the 'computational initial conditions' - could be written Fzi ( t , x y , ) , and are emphatically not the same as the traditional physical initial conditions at a time t = t i (written as e.g. Fti ( x y z , , ) ). Contrast, for example, the starting states on figs. 1 and 2.

Given some wave field state Fz ( t , x y , ) known at some specified path position z 0, we are at liberty to Fourier transform in x and y space, and over time t , but not along the propgation axis z . The result will be a mixed spectrum Fz ( w , kx , ky ) as calculable for the current position z 0, as well as prior locations z &lt; z 0. The spatial transform of this will give us the trans-

FIG. 2: Spatial propagation of waves, where disturbances (or pulses) evolve either forward or backward in time. At any point in the propagation, we know the full time behaviour of our wave field - both history and future, as indicated by the pale pink horizontal lines.. The light pink shaded triangles indicates the computational past of a wave element at particular points (black circles) along the path of the disturbance. Note that unlike the temporally propagated case shown in fig. 1, the computational past of a spatially propagated system is not the same as the causal past. A notional interface has been added to the diagram to show how reflections behave - i.e. in an unexpected way [4]. This is because a reflection should have been put in the initial conditions, but was not due to an (assumed) lack of knowledge of that future behaviour.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-a059ef45-1921-46ac-9ba8-4ebee8ca53fd.jpg)

verse wavevector spectrum of the wave field, which will contain both positive and negative wavevector components along each transverse coordinate axes. The temporal transform, using the known full past-and-future history of that transverse behaviour, will contain both positive and negative frequency components. This outcome is particularly convenient because (assumed) knowledge of the full time behaviour enables a frequency spectrum to be calculated, and arbitrary temporal material response (dispersion) to be implemented in a numerically efficient way using (see commentary in e.g. [4, 5, 14] and references therein).

For simplicity, and in analogy to the comparable discussion for the traditional picture, imagine we are following some trajectory through space and time, while counting interesting events that occur locally along that path. In this picture, the spatial coordinate along our propagation axis (e.g. z ) is always increasing but our trajectory can move forward or back along the other spatial axes ( x , y ), and even forward of back in time; although it is easier to think of the case where our x , y position is fixed, as is time t . Naturally, our count of interesting events will only ever increase, but rather than measuring second between events (or events per second), it would be in meters between events (or events per meter) - not a temporal period (or frequency), but a spatial interval l or spatial recurrence rate (wavevector) kz = 2 p / l . Consquently, any l (or kz ) estimation we might make would be a positive number. However, the events themselves may have different characteristics: notably, we may see objects that pass us travelling along a transverse spatial axis in either direction. Further, since Fz ( t , x y , )

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-1f77b1a0-617d-48b3-b8e6-5c24cafb0b58.jpg)

contains a full time history, this picture can even encode objects travelling forward or backward in time! The differing character of events, accessible to us through the spatial and temporal knowledge in Fz ( t , x y , ) , allows us to impose a sign (or signs) when adjusting our counting total - perhaps +1 for future-going objects, and -1 for past-going ones.

Having discussed how a counting/distance argument that gives us a strictly positive wavevector estimate can be converted into a signed count using temporal information, let us revist the spatially propagated frequency spectrum Fz 0 ( w ) of our notional wave field. As already noted, the frequency spectrum of the wave has both positive and negative frequency components. Again using the well known velocity relation for waves v = w / k , we can now claim that the positive frequency part of the spectrum represents forward (in time) evolving waves (with v &gt; 0), while the negative wavevector part represents backward (in time) evolving waves (with v &lt; 0). However, we have to be careful: even a static wave will have positive and negative frequency components, since for real-valued F t ( ) , F ( -w ) = F ∗ ( w ) . It is the changing complex phase of F ( w ) that represents the shift in F t ( ) profile from one time to another.

Note that I have so far only discussed estimates of the propagation axis wavevector kz . In this picture, any system state exists only at a specific location (e.g. z = z 0) along its path, and so of itself that state provides no information about a kz at all. Further, even given our knowledge of previous states on the path, we cannot calculate a true kz , because we only know about where we have been ( z ≤ z 0), not where we are yet to go ( z &gt; z 0). The best we might do is apply a Laplace transform over the data from behind us along our path, converting z into a Laplace-spectral qz , and get hybrid spectrum F z ′ 0 ( w , kx , ky , qz ) . Nevertheless, it is usually extremely advantageous to characterise the dominant kz -like properties of the propagation as a function of the known frequency and so be able to use a reference kz ( w ) [4, 5] as a basis on which to simplify the propagation - perhaps by assuming it to be unidirectional.

So in the directed picture, we can directly obtain a true frequency spectrum at any point in our propagation, and that spectrum will have negative components. These negative frequencies have a precise and well-defined meaning, but that meaning results from the convenient (but physically approximate) decision to treat propagation as if it were along a spatial path, rather than forward in time.

## IV. CAUSAL: THE PAST LIGHT CONE

The causal position taken when we claim all knowledge allowed by causal signalling about how the wave has moved through its environment. Here we might express our knowledge of some wave-field F by writing it F ( t , B ) ; where t and B are the proper time interval into the past and B the vector 'rapidity' needed for a signal to travel from the past to our current location. These t and B are constructed in a similar way to Rindler coordinates, and quite naturally respect the light cone. This has been discussed in more detail elsewhere

FIG. 3: Causal propagation, where our knowledge only advances with the edge of a single point's (observer's) lightcone; here I show the point travelling at slightly less than the maximum speed (i.e. of that of light) for clarity. At any point in the propagation, we know only the behaviour in our past lightcone; we show three such past lightcones, one being inscribed with curves showing how the rapidity B varies at a selection of fixed proper time t intervals. The initial conditions are a single point, and the final state the lightcone border. A notional interface has been added to the diagram to show how reflections behave.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-b1aaaea3-1c2a-465d-ae55-dafee5ccf208.jpg)

[15].

Since here our frequency-like quantity relates to the proper time t and not t , and because our wavevector-like quantity relates to B and not r , I leave any more systematic analysis to later work. However, note that t is like time t in the sense that it is one sided - we only know the past; and that B is like r in that spans all (allowed) points in space.

On fig. 3 I indicate in diagram form how a strictly causal simulation might proceed if we consider the knowledge of (an observer) at a single point e.g. one co-moving with a point on the wave profile. The starting state of this causal propagation matches what would we would know the instant we switched on some sensors - i.e. only that of the observer's current location. Here, the starting state for some causally propagated simulation - the 'computational initial conditions' - are a limited (single point) subset of the traditional physical initial conditions, which tend to assume a much greater knowledge of the environment. If desired, these causal initial conditions could be expanded to cover the past light cone of the initial point. This initial knowledge would then expand, because as we propagate, we would integrate our model of the system to add a 'new layer' to the outside of our prior past light cone to get our new (updated) past light cone.

It is debatable whether this rather purist causal picture is of much use in a scientific setting, except perhaps as part of a thought experiment. Nevertheless, it could be invaluable when considering cause and effect or the dynamical responses in specific situations, such as a metamaterial element being driven in the ultrafast and nanoscopic regime. On a more colourful level, for those engaged in spaceship combat at relativistic velocities - or more likely, those science fiction authors writing about such things - it is the only view of the

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-7e1f0b19-4490-4e9b-9555-16dccff10ed1.jpg)

known environment that makes sense.

## V. SUMMARY

The above discussion shows that negative frequencies can have a well grounded and physical basis - as long as we either are sufficiently omniscient and have a complete knowledge of the behaviour we consider relevant, or if we start from the premise that choosing spatial propagation is reasonable. In such cases we cannot object to the appearance of negative frequency components, but it is worth noting that that we are not omniscient, and that spatial propagation - however useful -

- [1] E. Rubino, J. McLenaghan, S. C. Kehr, F. Belgiorno, D. Townsend, S. Rohr, C. E. Kuklewicz, U. Leonhardt, F. König, and D. Faccio, Phys. Rev. Lett. 108 , 253901 (2012), doi:10.1103/PhysRevLett.108.253901.
- [2] M. Conforti, N. Westerberg, F. Baronio, S. Trillo, D. Faccio, Phys. Rev. A 88 , 013829 (2013), doi:10.1103/PhysRevA.88.013829.
- [3] G. Rousseaux, C. Mathis, P. Maïssa, T. G. Philbin, U. Leon-
- hardt, New J. Phys. 10 , 053015 (2008), doi:10.1088/1367-2630/10/5/053015.
- J. Phys. Commun. 10 , 025011 (2018), doi:10.1088/2399-6528/aaa85c,
- [4] P. Kinsler, arXiv:1202.0714.
- [5] P. Kinsler,
- Phys. Rev. A 81 , 013819 (2010), doi:10.1103/PhysRevA.81.013819, arXiv:0810.5689.
- [6] P. Kinsler,
- Directional pulse propagation in beam, rod, pipe, and disk geometries (2012),
- arXiv:1210.6794.
- [7] K. S. Yee,

IEEE Trans. Antennas Propagat.

14

, 302 (1966), doi:10.1109/TAP.1966.1138693

- [8] A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, Comput. Phys. Commun. 181 , 687 (2010), doi:10.1016/j.cpc.2009.11.008.
- [9]
- P. Kinsler, J. Opt. 20 , 025502 (2018),
- doi:10.1088/2040-8986/aaa0fc, arXiv:1501.05569.
- [10] K. J. Blow and D. Wood, IEEE J. Quantum Electronics 25 , 2665 (1989), doi:10.1109/3.40655.
- [11] M. Kolesik, J. V. Moloney, and M. Mlejnek, Phys. Rev. Lett. 89 , 283902 (2002), doi:10.1103/PhysRevLett.89.283902.
- [12] G. Genty, P. Kinsler, B. Kibler, and J. M. Dudley, Opt. Express 15 , 5382 (2007), doi:10.1364/OE.15.005382.

is (in practice) an approximation of reality. This means that the concept of negative frequencies must be treated with some caution in any kind of dynamical situation.

## Acknowledgments

I would like to acknowledge the role of the recent 'The Nonlinear Meeting 2014' in Edinburgh (May 2014), funded by the Max Planck Society, in helping crystallize some of my thoughts on this topic, and motivating me to write them up. I also acknowledge financial support from EPSRC [17].

[13] P. Kinsler, S. B. P. Radnor, and G. H. C. New, Phys. Rev. A 72 , 063807 (2005), doi:10.1103/PhysRevE.75.066603, arXiv:physics/0611215.

[14] J. C. A. Tyrrell, P. Kinsler, G. H. C. New, J. Mod. Opt. 52 , 973 (2005), doi:10.1080/09500340512331334086.

[15] P. Kinsler,

Eur. J. Phys.

32

, 1687 (2011),

doi:10.1088/0143-0807/32/6/022,

arXiv:1106.1792 (a longer and updated version).

[16] P. Kinsler, G. H. C. New, and J. C. A. Tyrrell, (2006),

arXiv:physics/0611213.

[17]

M. W. McCall, P. Kinsler, EPSRC grant number EP/K003305/1.

[18] J. R. Gulley and W. M. Dennis, Phys. Rev. A 81 , 033818 (2010), doi:10.1103/PhysRevA.81.033818.

[19] P. Kinsler and G. H. C. New, Phys. Rev. A 72 , 033804 (2005), doi:10.1103/PhysRevA.72.033804, arXiv:physics/0606111.

[20] R. H. Stolen, J. P. Gordon, W. J. Tomlinson, and H. A. Haus, J. Opt. Soc. Am. B 6 , 1159 (1989), doi:10.1364/JOSAB.6.001159.

[21] P. Kinsler, Phys. Rev. A 82 , 055804 (2010), doi:10.1103/PhysRevA.81.013819, arXiv:1008.2088.

[22]

P. Kinsler, An introduction to spatial dispersion (2019),

arXiv:1904.11957.

## Appendix: Time response models

The time response of a propagation medium is handled rather differently in the temporal and spatial propagation approaches. Consider a set of material properties ni that follow some kind of dynamical (temporal) response models defined by a differential equation. The differential equations for each ni will depends on both the local material states and that of the

local field:

$$\partial _ { t } ^ { k } n _ { i } = r _ { i } ( \{ n _ { j } \}, F ). \quad \quad \ ( 5. 1 ) \quad \partial _ { i }$$

Typical examples might be nonlinear response delays [10], a free carrier density model [18], a Raman response [19, 20], or even just a Drude or Lorentz oscillator [15] giving rise to dispersion. The requirement here for the model to be properly causal is that the order of the time derivative k is greater than any other time derivative parts present in the response function (or operator) ri [21].

- glyph[trianglerightsld] In the traditional time propagated picture, each spatial location needs to not only know its field state F ( r ) , but also its material property state(s) ni ( r ) . In some situations, values of F or ni from earlier times - i.e. values from the computational past of the simulation - may need to be stored for use in subsequent computations. Further, as the propagation proceeds step-by-step in time, both the field F ( r ) and current material properties ni ( r ) need to be updated.
- glyph[trianglerightsld] In the (directed) space propagated picture, each currently held state of the field holds the time history of each point. This means that it is never necessary to store the computational past of the simulation, since the values for earlier times are already incorporated into the current computational state F t ( , x y , ) . We can simply solve equations such as eqn. (5.1) by directly integrating them from the distant past up to the desired time, using only that current state information. And, as already mentioned, if the response model is linear and has a known frequency response, as in the case of typical dispersive behaviours, it can be applied quickly and efficiently as a simple pahse shift applied to the frequency spectrum - a process requiring only two Fourier transforms and a multiplication.

## Appendix: Nonlinear optics and negative frequencies

There are some specific issues relating to nonlinear optics (NLO) and negative frequencies which bear additional examination. Consider the nonlinear Schrödinger equation (NLSE) as derived from Maxwell's equations in the 1+1D regime of ( t , z ) , with some sources of dispersion and a perturbative third order (Kerr) nonlinearity. Generally in NLO it is written in the convenient (i.e. directed) spatially propagated picture, in a unidirectional approximation with fields evolving forward in time only, as the zNLSE [5]

$$\partial _ { z } E ^ { + } ( t ) = \imath K E ^ { + } ( t ) + \sum _ { j } \beta _ { j } \partial _ { t } ^ { j } E ^ { + } ( t ) + \imath \chi _ { z } \left | E ^ { + } ( t ) \right | ^ { 2 } E ^ { + } ( t ), \quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{$$

but note that it is also possible to derive a temporally propagated version, for the case of a unidirectional field evolving forward in space only, as the tNLSE 4

$$( 5. 1 ) \quad \partial _ { t } D _ { + } ( z ) = \imath \Omega D _ { + } ( z ) + \sum _ { j } \gamma _ { j } \partial _ { z } ^ { j } D _ { + } ( z ) + \imath \chi _ { t } | D _ { + } ( z ) | ^ { 2 } D _ { + } ( z ). \\ 0 1. \text{ a}$$

Note that in the 'directed' spatial zNLSE eqn. (5.2) all sources of dispersion - whether due to time-dependent material response or the geometric properties of a confining waveguide - are treated as if they were solely due to a timedependent response (by means of the power series in time derivatives). In contrast, in the 'traditional' temporal zNLSE eqn. (5.3), all sources of dispersion in the tNLSE are treated as if they were solely due to the local geometric properties of that material i.e. by spatial dispersion [22]; typically this is expressed as a power series in either spatial derivatives or wavevector. However, given that most dispersions are rather weak in applications NLSE propagation, either picture can plausibly be assumed to be able to treat dispersion accurately enough for such purposes (although see [9]).

In both these cases we can see that it would be useful to Fourier transform the field to get temporal and spatial spectra ˜ E + ( w ) and ¯ D +( ) k respectively, because then the dispersion model can be directly applied as a simple polynomial phase shift, i.e. either of

$$\beta ( \omega ) = \sum _ { j } \beta _ { j }. \, ( - \imath \omega ) ^ { j } \, \text{ \quad \ \ } \text{(5.4)}$$

$$\gamma ( k ) = \sum _ { j } \gamma _ { j }. \, ( \imath k ) ^ { j }$$

It is often stated that a Kerr nonlinearity causes third harmonic generation, but we now also need to consider what this generation is a 'third harmonic' of. In the usual zNLSE model of eqn. (5.2), an existing harmonic field E + ( ) t glyph[similarequal] E + 0 e -ı w 1 t + c.c. leads to not only some self phase modulation (SPM) but also third harmonic frequency generation at w 3 = 3 w 1. In the tNLSE model of eqn. (5.3), an existing harmonic field D +( ) z glyph[similarequal] D e ık z 0 + 1 + c.c. leads to not only some self phase modulation (SPM) but also third harmonic wavevector generation at k 3 = 3 k 1.

For clarity, let us write this out as carefully as possible, in the case where we spilt the real-valued fields E + or D + into complex conjugate halves, in the style of [16]. With E + ( ) = t E ′ ( ) t e -ı w 1 t + E ′∗ ( ) t e + ı w 1 t , the zNLSE equation can be partitioned into two complex conjugate halves. We can also dispense with the exponential oscillation by setting w 1 = 0 although it is often useful to include such carrier oscillations when the fields are narrow band, having them written explicitly rather than implicitly is not required. The two complex

4 Author's derivation, unpublished

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0128v3-figure-ed4fee27-f1c6-44e0-8ac1-74e767a646ec.jpg)

conjugate zNLSE equations are

$$\partial _ { z } E ^ { \prime } ( t ) e ^ { - i \omega _ { 1 } t } = & \, \imath K E ^ { \prime } ( t ) e ^ { - i \omega _ { 1 } t } + \sum _ { j } \beta _ { j } \, \partial _ { t } ^ { j } E ^ { \prime } ( t ) e ^ { - i \omega _ { 1 } t } \\ & + \imath \chi _ { z } \left [ E ^ { \prime } ( t ) ^ { 3 } e ^ { - 3 i \omega _ { 1 } t } + 3 E ^ { \prime } ( t ) ^ { 2 } \left [ E ^ { \prime } ( t ) \right ] ^ { * } e ^ { - i \omega _ { 1 } t } \right ], \quad \text{ might } \, \imath \, \text{ be to } \, \i \, \text{struct t} \, \i }$$

$$\begin{matrix} \end{matrix}$$

$$\partial _ { z } E ^ { \prime * } ( t ) e ^ { + \imath \omega _ { 1 } t } = \imath K E ^ { \prime * } ( t ) e ^ { + \imath \omega _ { 1 } t } + \sum _ { j } \beta _ { j } \partial _ { t } ^ { j } E ^ { \prime * } ( t ) e ^ { + \imath \omega _ { 1 } t }$$

$$+ \imath \chi _ { z } \left [ E ^ { \prime * } ( t ) ^ { 3 } e ^ { + 3 \imath \omega _ { 1 } t } + 3 E ^ { \prime * } ( t ) ^ { 2 } \left [ E ^ { \prime * } ( t ) \right ] ^ { * } e ^ { + \imath \omega _ { 1 } t } \right ] \underset { \substack { \text{num} \\ ( 5. 7 ) } } { \text{num} }.$$

The sum of these two equations is just eqn. (5.2). Further, since they are exact complex conjugates of one another, we can solve just one version, which automatically gives us a solution to the other, and hence the solution to the real valued field E + ( ) t .

Note in particular that the nonlinearity in this exact mathematical re-expression of the zNLSE equation drives only resonant or third-harmonic frequencies. There is no explicit coupling to negative frequencies, which can be seen most clearly when assuming finite w 1 and a constant E ′ ; then, eqn. (5.6) does not drive negative frequencies of itself. Nevertheless, although the complex E ′ ( w ) might start with content solely in positive frequencies, its frequency bandwidth is not restricted. Consequently after propagating some distance it may well have evolved into a field whose spectrum is very wide, perhaps even extending past the origin. In this state, multifrequency interactions that drive the field at negative frequencies could indeed occur as a result of the nonlinearity - remember that if represented in the frequency domain, the cubic nonlinear term transforms into a double convolution over the entire frequency spectra.

However, we could, if we wanted, make such 'negative frequency' driving terms appear explicitly by repartitioning the whole nonlinear term into different complex conjugate halves, e.g. by writing not eqn. (5.6) but

$$\partial _ { z } E ^ { \prime } ( t ) e ^ { - \imath \omega _ { 1 } t } = \imath K E ^ { \prime } ( t ) e ^ { - \imath \omega _ { 1 } t } + \sum _ { \colon } \beta _ { j } \partial _ { t } ^ { j } E ^ { \prime } ( t ) e ^ { - \imath \omega _ { 1 } t }$$

$$+ \imath \chi _ { z } \left [ E ^ { \prime } ( t ) ^ { 3 } e ^ { - 3 \imath \omega _ { 1 } t } + \frac { 3 } { 2 } E ^ { \prime } ( t ) ^ { 2 } \left [ E ^ { \prime } ( t ) \right ] ^ { * } e ^ { - \imath \imath \omega _ { 1 } t } \quad \text{As} \\ + \frac { 3 } { 2 } \left [ E ^ { \prime } ( t ) \right ] ^ { * 2 } E ^ { \prime } ( t ) e ^ { + \imath \imath \omega _ { 1 } t } \right ], \quad \text{of} \ \imath \\ \quad \text{$\quad \text{$\quad \colon c.}$}$$

$$& \stackrel { | t } { \L } + \sum _ { j } \beta _ { j } \, \partial _ { t } ^ { J } \, E ^ { \prime } ( t ) e ^ { - \iota \omega _ { 1 } t } \\ & e ^ { - 3 \iota \omega _ { 1 } t } + \frac { 3 } { 2 } E ^ { \prime } ( t ) ^ { 2 } \left [ E ^ { \prime } ( t ) \right ] ^ { * } e ^ { - \iota \omega _ { 1 } t } \\ & + \frac { 3 } { 2 } \left [ E ^ { \prime } ( t ) \right ] ^ { * 2 } E ^ { \prime } ( t ) e ^ { + \iota \omega _ { 1 } t } \right ], \quad \, \stackrel { | t } { \L } \\ & \quad ( 5. 8 ) \quad \, \text{$1$}$$

]

. and there is of course also a matching complex congugate counterpart of eqn. (5.7), and the sum of both will be equal to the original zNLSE equation. Whether or not this representation containing explicit driving of negative frequencies might be useful is another matter, but it is certainly possible to (re)express the mathematical model in order to construct them. But, apart from specific numerical difficulties that might occur when solving these propagation equations, the solutions gained from either form should be identical: both are just different ways of representing the same physical model .

Naturally one can apply the same method to the tNLSE as well, setting D +( ) = z D z e ık t ′ ( ) 1 + D ′∗ ( ) z e -ık 1 z , and arriving at two complex conjugate equations

$$\underset { \text{re} - } { \text{ley} } \quad \partial _ { z } D ^ { \prime } ( z ) e ^ { + \imath k _ { 1 } z } = \imath K D ^ { \prime } ( z ) e ^ { + \imath k _ { 1 } z } + \sum _ { z } \gamma _ { j } \partial _ { z } ^ { j } D ^ { \prime } ( z ) e ^ { + \imath k _ { 1 } z }$$

$$\frac { \Gamma \cos ^ { - } } { - \underset { \dots } { \mathfrak { e } } ^ { - } } \partial _ { z } D ^ { \prime * } ( z ) e ^ { - \imath k _ { 1 } z } = \imath K D ^ { \prime * } ( z ) e ^ { - \imath k _ { 1 } z } + \sum _ { \imath } \gamma _ { j } \partial _ { z } ^ { j } D ^ { \prime * } ( z ) e ^ { - \imath k _ { 1 } z }$$

$$+ \imath \chi _ { t } \left [ D ^ { \prime * } ( z ) ^ { 3 } e ^ { - 3 t k _ { 1 } z } + 3 D ^ { \prime * } ( z ) ^ { 2 } \left [ D ^ { \prime * } ( z ) \right ] ^ { * } e ^ { - t k _ { 1 } z } \right ].$$

$$\L \L e ^ { + t k _ { 1 } z } & = \imath K D ^ { \prime } ( z ) e ^ { + t k _ { 1 } z } + \sum _ { j } \gamma _ { j } \partial _ { z } ^ { j } D ^ { \prime } ( z ) e ^ { + t k _ { 1 } z } \\ & \quad + \imath \chi _ { t } \left [ D ^ { \prime } ( z ) ^ { 3 } e ^ { + 3 t k _ { 1 } z } + 3 D ^ { \prime } ( z ) ^ { 2 } \left [ D ^ { \prime } ( z ) \right ] ^ { * } e ^ { + t k _ { 1 } z } \right ], \\ \L \L e ^ { - t k _ { 1 } z } & = \imath K D ^ { \prime * } ( z ) e ^ { - t k _ { 1 } z } + \sum _ { j } \gamma _ { j } \partial _ { z } ^ { j } D ^ { \prime * } ( z ) e ^ { - t k _ { 1 } z } \\ & \quad + \imath \chi _ { t } \left [ D ^ { \prime * } ( z ) ^ { 3 } e ^ { - 3 t k _ { 1 } z } + 3 D ^ { \prime * } ( z ) ^ { 2 } \left [ D ^ { \prime * } ( z ) \right ] ^ { * } e ^ { - t k _ { 1 } z } \right ].$$

As we should expect, these two equations sum to eqn. (5.3). Also, here there again is no explicit coupling to opposite parts of the wavevector spectra D k ′ ( ) in each equation, although (as before) we might repartition the nonlinear term to construct it, if we so desired.
Paul Kinsler
Time and space, frequency and wavevector: or, what I talk about when I talk about propagation
The existence of waves with negative frequency is a surprising and perhaps
controversial claim which has recently been revisted in optics and for water
waves. Here I explain a context within which to understand the meaning of the
"negative frequency" conception, and why it appears in some cases.