I first have to say that I only have a limited education in physics, so I hope this is not a nonsensical question.

I'm trying to understand the connection between the wave equation

$$\partial_{tt} u - c^2 \Delta u = f \tag{W}\label{W}$$

and the Helmholtz equation

$$\Delta U + k^2 U= -\frac{1}{c^2} F. \tag{H}\label{H} $$

I think I have quite a good intuition how the wave equation $\eqref{W}$ works:

If we stimulate our medium with some $f$, this "information" is propagated in all directions with a certain velocity $c$.

Then I read that the Helmholtz equation is derived by assuming that

$$u(x,t) = U(x)e^{-i\omega t} \tag{*}\label{*}$$

(and similarly $f(x,t) = F(x)e^{-i\omega t}$) which, when plugged into $\eqref{W}$ results in $\eqref{H}$ for $k = \frac{\omega}{c}$.

But I fail to see how this encodes some kind of propagation as in $\eqref{W}$. This leads to my first question:

In what physical situation does the assumption $\eqref{*}$ make sense?

Then I also read that for $F=0$ a solution of $\eqref{H}$ describes an oscillation mode. By looking at pictures of solutions I saw how this shows e.g. how a membrane or a string vibrates. But it seems this only works for some $\omega$ (in the string example you need nodes at the boundary) But what I do not understand about that is

How is $\omega$ determined or defined?

  • $\begingroup$ In my understanding $\omega$ is a parameter of the wave. The solution(s) $U(x)$ will be dependent on $\omega$. This can be interpreted as dispersion, different frequencies of the wave have different behavior. Dispersion $\endgroup$
    – Korf
    May 19, 2017 at 12:56
  • $\begingroup$ @Korf Note that dispersive media don't fully fit into the wave equation as posed by OP. Operatively, you describe dispersion via the Helmholtz equation but modifying the relationship between $k$ and $\omega$, but if you want to derive that directly from the wave equation, that's a whole lot of work, needing some form of dimensional elimination (like e.g. dispersive 2D shallow-water waves derived from 3D bulk dynamics) or coupling to an external dof (like light coupling to the polarization of a medium via a a causal, linear constitutive equation), none of which is simple. $\endgroup$ May 19, 2017 at 16:45
  • $\begingroup$ (... which is, of course, why physicists often just say, here's my dispersion relation $\omega(k)$, it works and it fits experiment, so just run with it and don't ask me where I got it.) $\endgroup$ May 19, 2017 at 16:46
  • $\begingroup$ @EmilioPisanty Thank you for clearing this. I knew this only from the point of view presented in your last sentence. $\endgroup$
    – Korf
    May 20, 2017 at 23:12

1 Answer 1


The passage from the full time-dependent wave equation $(\mathrm{W})$ to the Helmholtz equation $(\mathrm{H})$ is nothing more, and nothing less, than a Fourier transform. For sufficiently regular functions, both $u$ and $F$ can be written as superpositions of monochromatic fields, i.e. they are both Fourier transforms of the form $$ u(x,t) = \int_{-\infty}^\infty U(x,\omega) e^{-i\omega t} \mathrm d\omega \quad\text{and} \quad f(x,t) = \int_{-\infty}^\infty F(x,\omega) e^{-i\omega t} \mathrm d\omega, $$ where the temporal Fourier coefficients $U(x,\omega)$ and $F(x,\omega)$ depend on the position - or, switching perspectives, they give us functions of $x$ for each $\omega$. Now, all we've done so far is a fancy rewriting of our variables, but there are two crucial aspects of the wave equation that make this useful:

  • it is linear, and
  • the only dependence on time is through $\partial_t^2$, which is a linear operator whose eigenfunctions are precisely the Fourier kernel, i.e. $\partial_t^2 e^{-i\omega t} = -\omega^2 e^{-i\omega t}$.

The linearity allows us to break in the wave equation's linear operators all the way through to the Fourier coefficients, and the eigenvalue relation for $\partial_t$ enables us to switch that partial differentiation to an algebraic factor on that sector, giving us \begin{align} 0 & = -\partial_{t}^2 u(x,t) + c^2 \nabla^2 u(x,t) + f(x,t) \\ & = -\partial_{t}^2 \int_{-\infty}^\infty U(x,\omega) e^{-i\omega t} \mathrm d\omega + c^2 \nabla^2 \int_{-\infty}^\infty U(x,\omega) e^{-i\omega t} \mathrm d\omega + \int_{-\infty}^\infty F(x,\omega) e^{-i\omega t} \mathrm d\omega \\ & = \int_{-\infty}^\infty \left[ -U(x,\omega) \partial_{t}^2 e^{-i\omega t} + c^2 \nabla^2 U(x,\omega) e^{-i\omega t} + F(x,\omega) e^{-i\omega t} \vphantom{\sum}\right]\mathrm d\omega \\ & = \int_{-\infty}^\infty \left[ \omega^2U(x,\omega) + c^2 \nabla^2 U(x,\omega) + F(x,\omega) \vphantom{\sum}\right] e^{-i\omega t} \mathrm d\omega . \tag{1} \end{align} From here, it's easy to see that if $f(x,t)$ is given (so $F(x,\omega)$ is also given), we can find a solution of the original equation by setting $U(x,\omega)$ to be a solution of the Helmholtz equation, $$ c^2 \nabla^2 U(x,\omega) + \omega^2U(x,\omega) = - F(x,\omega) $$ or with the cosmetic change $k=\omega/c$, $$ \nabla^2 U(x,\omega) + k^2U(x,\omega) = - \frac{1}{c^2} F(x,\omega). $$ (In addition, it's easy to show that the Fourier transform in $(1)$ means that this is a necessary condition, but if all you're doing is finding solutions, as opposed to characterizing the general solution, then the sufficiency is enough.)

OK, so that is the formal side. How is this used in the real world? There are three main ways that one uses this.

  1. The clearest is when the wave equation is being forced by a source that is itself monochromatic (or close enough to monochromatic that your situation doesn't care about the difference), or in terms of the Fourier amplitude $F(x,\omega) = F(x) \delta(\omega-\omega_0)$. This is the case, for example, when one considers the electromagnetic emission of an antenna set to a very narrow band of frequencies.

    Physically speaking, the Helmholtz equation $(\mathrm{H})$ does encode propagation, in a very real sense ─ except that you're considering in one single go the coherent superposition of the emission that comes from a source that is always turned on, and oscillating at a constant frequency for all time. In this case, you expect the physical response to be at that same frequency, but the spatial response can be complicated in the presence of reflections, dispersive media, or whatnot; we solve the Helmholtz equation to find that spatial response.

    In this case, $\omega$ is obviously fixed by the external driver.

  2. A separate application is when we solve for resonant modes of the domain in question; these are nonzero solutions to the Helmholtz equation that hold even when the driver $F$ is zero, and they are important e.g. when $F$ is an impulse that's confined in time, like hitting a drum, and the effects are left to resonate in a confined domain which the energy cannot leave easily.

    In this case, $\omega$ is fixed by the domain to one of a discrete set of resonant frequencies that sustain nonzero solutions of $(\mathrm{H})$ even when the driver is zero, and the process of solving the Helmholtz equation includes finding those resonant frequencies. The end goal in this calculation is a set of resonant frequencies $\{\omega_n\}$ with a corresponding set of solutions $\{u_n(x)\}$ which satisfy the homogeneous Helmholtz equation at that frequency and which form a complete basis, in the $L_2$ sense, for functions over the domain in question.

    To reconcile this with the driver, the simplest case is to consider an impulsive driver, i.e. something of the form $f(x,t) = f(x)\delta(t)$, with a flat Fourier transform. In this case, all modes see the impulse, but only the resonant modes are able to respond. In this case, you decompose $f(x)$ as a linear combination of the $u_n(x)$, and this tells you how much each mode gets excited, which determines the temporal evolution after the impulse is gone.

    Does this describe "propagation" in a suitable sense? Well, you're ultimately solving for the propagation of an initial impulsive disturbance, like plucking a string, by finding a clever decomposition of that initial disturbance in terms of modes that evolve cleanly (monochromatically) in time. So, yes.

  3. Finally, there is also the case where you just have some arbitrary driver $f(x,t)$ for the wave equation, and all that you can say about its Fourier transform is that it exists. In this case, $\omega$ isn't 'chosen', as such: instead, it is a continuous parameter of the problem, where you solve a continuous set of separate inhomogeneous Helmholtz equations to get the $U(x,\omega)$, and then you add them all up coherently to get $u(x,t) = \int_{-\infty}^\infty U(x,\omega) e^{-i\omega t} \mathrm d\omega$.

    It's important, however, not to underestimate the important of what you can say about $F$: just by saying "the temporal Fourier transform of $f(x,t)$ exists", you're saying that $f(x,t)$ can be understood as a superposition of monochromatic waves, each of which can be solved independently and which will cause some monochromatic response $U(x,\omega) e^{-i\omega t}$, which can then be added together to give the global response to the driver.

  • $\begingroup$ I cannot thank you enough for this extensive answer, you fully answered all my questions! So thanks a lot! $\endgroup$
    – flawr
    May 19, 2017 at 16:36
  • $\begingroup$ No worries. I will chide you for using $\Delta$ for the laplacian in a physics venue, though :-P. $\endgroup$ May 19, 2017 at 16:37
  • $\begingroup$ I get the impression that physicists are giving math-people a hard time on purpose :D $\endgroup$
    – flawr
    May 19, 2017 at 16:41
  • $\begingroup$ Fantastic answer. As someone who works with electromagnetics from an RF engineering background, let me add another practical aspect. Us electrical engineers basically think in the frequency domain. Material properties are defined as functions of frequency, and we almost always look at observable quantities (e.g. antenna input impedance) as functions of frequency, so it makes sense to consider Maxwell's equations in the frequency domain (i.e. time harmonic). This is kind of a corollary to the first example. $\endgroup$
    – LedHead
    May 19, 2017 at 20:04
  • $\begingroup$ @led23head That's a bit more complicated, since many relevant frequency-dependent characteristics (particularly the permittivity) can't be described by a simple wave equation as in the OP. Instead, your problem is the wave equation for $E$, coupled with a dynamical equation for $P$, and you solve it by going to the Fourier domain first, then solving for the (frequency-dependent) constitutive relation, and only then do you get a wave problem that's purely in $\mathbf E$ - as a Helmholtz equation, not a wave equation. So it's more the third case than the first. $\endgroup$ May 19, 2017 at 20:12

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.