# Wave packet expression and Fourier transforms

I read in the Wikipedia article on wave packets that a wave packet is defined as: $$\psi(x,t)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{+\infty}g(k)e^{i(kx-\omega t)}\mathrm dk\tag a$$ where $$g(k)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{+\infty}\psi(x,0) e^{-ikx}\mathrm dx\tag b$$ Now I'd like to understand mathematically where these formulas come from.

As an initial example, if you put $t=0$ in $\rm (a)$, you obtain $$\psi(x,0)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{+\infty}g(k)e^{ikx}\mathrm dk\tag c$$ and $$g(k)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{+\infty}\psi(x,0) e^{-ikx}\mathrm dx,\tag d$$ and these last two formulas are clear $-$ it is the Fourier reverse transform theorem.

However, I don't understand where $\rm (a)$ comes from. I thought that maybe $(a)$ was expressed as the Fourier reverse transform of the function of the two variable $\psi(x,t)$, but in this case I obtain $$\psi(x,t)=\frac{1}{{2\pi}}\iint_{R^{2}}{\mathcal{F}[\psi(x,t)](w,k)} e^{i(kx+\omega t)}\:\mathrm dk\:\mathrm d\omega, \tag e$$ which is different from $\rm (a)$. So, finally, my question is: where does $\rm (a)$ comes from?

• Read where? Which book? How a bit different? Difference between which two formulas? Are you asking about sign conventions etc? Aug 6, 2017 at 10:19
• @Qmechanic Read on my book, but also here: en.wikipedia.org/wiki/Wave_packet ; with "bit different" I meant "different" (sorry); I'm not asking about conventions, only the the mathematical method that allows you to write a wave packet in that form. I thought that this method was a simply inverse Fourier transform, but it is not because it doesn't give the same result. Aug 6, 2017 at 19:43
• Mathematically, it's like a generalization of a Fourier expansion: you suppose you can expand your function as an infinite sum of complex exponentials but like a function in general is not periodic you take the period to be infinite and you get the Fourier transform. Aug 7, 2017 at 21:22

When we write down a wavepacket, we're trying to solve the following problem: Given an initial profile for our wavepacket, $\psi(x,0)$, what will it look like in the future? In other words, we want to find $\psi(x,t)$ given $\psi(x,0)$. To solve this problem, we need to have some wave equation that tells us how the wavepacket evolves in time. For different waves in different contexts, we'll have different wave equations. But we'll assume two things about the wave equation:

• It is linear. This means that if $\psi_1(x,t)$ is a solution and $\psi_2(x,t)$ is a solution, then $a_1\psi_1(x,t)+a_2\psi_2(x,t)$ is a solution. Generalizing, it means that if $\psi_k(x,t)$ is a solution, so is $\int dk\ a_k\psi_k(x,t)$.

• If we start with an initial profile $\psi(x,0)=e^{ikx}$, then the solution to our wave equation is $\psi(x,t)=e^{i(kx-\omega_k t)}$, where $\omega_k$ is a constant that may depend on $k$.

You can check for yourself that generic wave equations, such as the equation for a wave on a string, or the Schrodinger equation, satisfy these properties. Now, once we know the wave equation has these properties, we can solve our problem!

First, we realize that by property (2), all functions of the form $e^{i(kx-\omega_k t)}$ are solutions to our wave equation. That means that, by property (1), any linear combination of these functions is also a solution to our wave equation. So we can write down a guess for the solution:

$$\psi(x,t)=\frac{1}{\sqrt{2\pi}}\int dk\ g(k)e^{i(kx-\omega_kt)}$$

So far, $g(k)$ is just some unknown function. Any choice of $g(k)$ will give a solution to the wave equation, by properties (1) and (2). But we also need our solution to match with our initial profile, $\psi(x,0)$. Plugging in $t=0$ thus gives a condition on $g(k)$:

$$\psi(x,0)=\frac{1}{\sqrt{2\pi}}\int dk\ g(k)e^{ikx}$$

By Fourier's theorem, we can then immediately find the $g(k)$.

$$g(k)=\frac{1}{\sqrt{2\pi}}\int dx\ \psi(x,0)e^{-ikx}$$

• Note that the expression for $g(k)$ refers to a specific time, $t=0$. So you shouldn't expect to be able to get $g(k)$ by fourier transforming time, since fourier transforming time eliminates the reference to a specific $t$ point. $g(k)$ is defined at $t=0$, not for general times. Aug 7, 2017 at 18:42
• Why do you write the linear combination between solutions first as a summation and after as an integral? Aug 7, 2017 at 19:13
• @Landau It's equivalent. If you know the sum of two solutions is a solution, it immediately follows that the sum of $N$ solutions is a solution. We want to sum up all the possible solutions (every k), and since $k$ is continuous we use an integral. Aug 7, 2017 at 19:58
• @Landau In general, linear differential equations allow you to integrate their solutions because integration and differentiation commute ($\frac{d}{dx}\int dk = \int dk\ \frac{d}{dx}$), so that you can integrate a family of solutions to get a new solution. Aug 7, 2017 at 20:00
• @Landau If you like, I'm assuming the wave equation satisfies the integral version of linearity. I don't care about the discrete version. Aug 8, 2017 at 0:52

Well, as to where they come from...how about I give you an explanation of what each means in the context of quantum mechanics/physics?

In 1-D, any analytic function $\psi(x)$ on $(-\infty, +\infty)$ can be written as the infinite sum of sinusoidal functions.

Now suppose our function $\psi(x)$ is a $t=0$ "snapshot" of a time dependent function $\psi(x,t)$.

Using this snapshot, your 2nd integral, known as the Fourier transform of $\psi(x,0)$, gives you the relative strengths $g(k)$ of each wave in the infinite sum making up $\psi$. The actual amplitude of each contributing wave is $g(k)dk$, which means it's an infinitesimal amount, but some waves contribute a "bigger" infinitesimal than others. Here, $k$ is real, but $g(k)$ is complex-valued in general.

Now suppose that all of the sinusoidal functions that sum up to the $t=0$ function are actually $t=0$ snapshots of traveling harmonic waves, where each harmonic wave moves at a phase velocity equal to $\omega/k$ when the clock is running. In much of physics, it is usually the case that $\omega$ is a function of $k$. If $\omega(k)=ck$, where $c$ is a constant, then we can see that the phase velocity of every harmonic wave is the same; it's just $c$. However, if $\omega$ is a more complicated function of $k$, then each wave that contributes to the over all $\psi$ function will travel at a different phase velocity. This is an important concept that explains why wave packets representing particles spread out over time.

So what your first integral says is that the time-dependent complex wave function $\psi(x,t)$ is the infinite sum of traveling harmonic waves, each having a complex relative amplitude of $g(k)$.

If the component waves making up your $t=0$ wave function all travel at the same phase velocity, then you will see your wave function travel at the same velocity without changing its shape. However, if the component waves have a range of differing phase velocities, this means that your wave function will change it's shape over time.

If the g(k) function looks like a "bell curve" centered about some average $k_0$, then the real part of your time-dependent wave function will look like an oscillating wave packet that widens and decreases in amplitude over time as it travels.

• Ok, the intuitive meaning of those formulas was quite clear, but I'm looking for the mathematical method that allows you to write a wave packet in that form. Aug 6, 2017 at 15:49
• @Landau I am a bit confused as to what you mean by mathematical method that allows. I wrote an answer which I deleted because it was too basic and did not address your full question. You have the intuitive (physical) picture, so which part, (and I mean this without offence, seriously ), of the simple extension of the Fourier series idea to wave packets is the problem? Have you checked with MathSE, as you have the physical idea, and MathSE deals with formal derivations and proofs, if that is what you are after?
– user163104
Aug 6, 2017 at 16:18
• @Landau are you asking why these are the correct formulas? If so, then the answer has to do with how you change basis states in a Hilbert space from a position basis to a momentum basis or the other way around. Aug 6, 2017 at 20:45