# Time-Dependent Solutions: Propagators and Representations

*Michael Fowler, UVa *

### Introduction

We’ve spent most of the course so far concentrating on the eigenstates of the Hamiltonian, states whose time-dependence is merely a changing phase. We did mention much earlier a superposition of two different energy states in an infinite well, resulting in a wave function sloshing backwards and forwards. It’s now time to cast the analysis of time dependent states into the language of bras, kets and operators. We’ll take a time-independent Hamiltonian $H,$ with a complete set of orthonormalized eigenstates, and as usual

$$i\hslash \frac{\partial \psi \left(x,t\right)}{\partial t}=-\frac{{\hslash}^{2}}{2m}\frac{{\partial}^{2}\psi \left(x,t\right)}{\partial {x}^{2}}+V\left(x\right)\psi \left(x,t\right),$$

Or, as we would now write it

$$i\hslash \frac{\partial}{\partial t}|\psi \left(x,t\right)\rangle =H|\psi \left(x,t\right)\rangle .$$

Since $H$ is itself time independent, this is very easy to integrate!

$|\psi \left(x,t\right)\rangle ={e}^{-iH\left(t-{t}_{0}\right)/\hslash}|\psi \left(x,{t}_{0}\right)\rangle .$

The exponential operator that generates the
time-dependence is called the *propagator*,
because it describes how the wave propagates from its initial configuration,
and is usually denoted by $U:$

$|\psi \left(x,t\right)\rangle =U\left(t-{t}_{0}\right)|\psi \left(x,{t}_{0}\right)\rangle .$

It’s appropriate to call the propagator $U$, because it’s a unitary operator:

$U\left(t-{t}_{0}\right)={e}^{-iH\left(t-{t}_{0}\right)},\text{so}\text{\hspace{0.17em}}{U}^{\u2020}\left(t-{t}_{0}\right)={e}^{i{H}^{\u2020}\left(t-{t}_{0}\right)}={e}^{iH\left(t-{t}_{0}\right)}={U}^{-1}\left(t-{t}_{0}\right).$

Since $H$ is Hermitian, $U$ is unitary. It immediately follows that

$\begin{array}{c}\langle \psi \left(x,t\right)|\psi \left(x,t\right)\rangle =\langle \psi \left(x,{t}_{0}\right)|{U}^{\u2020}U\left(t-{t}_{0}\right)|\psi \left(x,{t}_{0}\right)\rangle \\ =\langle \psi \left(x,{t}_{0}\right)|\psi \left(x,{t}_{0}\right)\rangle \end{array}$

the norm of the ket vector is conserved, or, translating to wave function language, a wave function correctly normalized to give a total probability of one stays that way. (This can also be proved from the Schrödinger equation, of course, but this is quicker.)

This is all very succinct, but unfortunately the
exponential of a second-order differential operator doesn’t sound too easy to
work with. Recall, though, that any
function of a Hermitian operator has the same set of eigenstates as the
original operator. This means that the
eigenstates of ${e}^{-iH\left(t-{t}_{0}\right)/\hslash}$ are the *same*
as the eigenstates of $H$,
and if $H|{\psi}_{n}\rangle ={E}_{n}|{\psi}_{n}\rangle ,$ then

${e}^{-iH\left(t-{t}_{0}\right)/\hslash}|{\psi}_{n}\rangle ={e}^{-i{E}_{n}\left(t-{t}_{0}\right)/\hslash}|{\psi}_{n}\rangle .$

This is of course nothing but the time dependent phase factor for the eigenstates we found before$\u2014$and, as before, to find the time dependence of any general state we must express it as a superposition of these eigenkets, each having its own time dependence. But how do we do that in the operator language? Easy: we simply insert an identity operator, the one constructed from the complete set of eigenkets, thus:

$\begin{array}{c}|\psi \left(t\right)\rangle ={e}^{-iH\left(t-{t}_{0}\right)/\hslash}{\displaystyle \sum _{n=1}^{\infty}|{\psi}_{n}\rangle \langle {\psi}_{n}}|\psi \left({t}_{0}\right)\rangle \\ ={\displaystyle \sum _{n=1}^{\infty}{e}^{-i{E}_{n}\left(t-{t}_{0}\right)/\hslash}|{\psi}_{n}\rangle \langle {\psi}_{n}}|\psi \left({t}_{0}\right)\rangle .\end{array}$

Staring at this, we see that it’s just what we had before: at the initial time $t={t}_{0},$ the wave function can be written as a sum over the eigenkets:

$\begin{array}{c}|\psi \left({t}_{0}\right)\rangle ={\displaystyle \sum |{\psi}_{n}\left({t}_{0}\right)\rangle}\langle {\psi}_{n}\left({t}_{0}\right)|\psi \left({t}_{0}\right)\rangle \\ ={\displaystyle \sum {c}_{n}|{\psi}_{n}\left({t}_{0}\right)\rangle}\end{array}$

with ${c}_{n}=\langle {\psi}_{n}|\psi \rangle $, $\sum {\left|{c}_{n}\right|}^{2}=1$, and the usual generalization for continuum eigenvalues, and the time development is just given by inserting the phases:

$|\psi \left(t\right)\rangle ={\displaystyle \sum {c}_{n}{e}^{-i{E}_{n}\left(t-{t}_{0}\right)/\hslash}|{\psi}_{n}\left({t}_{0}\right)\rangle}.$

The expectation value of the energy $E$ in $|\psi \rangle $,

$\langle E\rangle =\langle \psi |H|\psi \rangle ={\displaystyle \sum {\left|{c}_{n}\right|}^{2}{E}_{n}}$

and is (of course) time independent.

The expectation value of the particle position *x*
is

$\langle \psi \left(t\right)|x|\psi \left(t\right)\rangle ={\displaystyle \sum _{n,m}{c}_{n}^{*}{c}_{m}{e}^{i({E}_{n}-{E}_{m})\left(t-{t}_{0}\right)/\hslash}\langle {\psi}_{n}\left({t}_{0}\right)|x|{\psi}_{m}\left({t}_{0}\right)\rangle}$_{}

and is *not* in general time-independent. (It *is* real, of course, on adding the $n,m$ term to the $m,n$ term.)

This analysis is only valid for a *time-independent*
Hamiltonian. The important extension to a system in a time-dependent external
field, such as an atom in a light beam, will be given later in the course.

### The Free Particle Propagator

To gain some insight into what the propagator $U$ looks like, we’ll first analyze the case of a particle in one dimension with no potential at all. We’ll also take ${t}_{0}=0$ to make the equations less cumbersome.

For a free particle in one dimension $E={p}^{2}/2m={\hslash}^{2}{k}^{2}/2m$ the energy eigenstates are also momentum eigenstates, we label them $|k\rangle $, so

$$\begin{array}{c}U\left(t\right)={e}^{-iHt/\hslash}={\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{e}^{-iH\text{\hspace{0.05em}}t/\hslash}\frac{dk}{2\pi}|k\rangle \langle k|}\\ ={\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{e}^{-i\hslash \text{\hspace{0.05em}}{k}^{2}t/2m}\frac{dk}{2\pi}|k\rangle \langle k|}.\end{array}$$

Let’s consider (following Shankar and others) what seems the simplest example.

Suppose that at $t={t}_{0}=0$, a particle is at ${x}_{0}$: $\psi (x,t=0)=\delta (x-{x}_{0})=|{x}_{0}\rangle $: what is the probability amplitude for finding it at $x$ at a later time $t$? (This would be just its wave function at the later time.)

$$\begin{array}{c}\langle x|U\left(t,0\right)|{x}_{0}\rangle \text{\hspace{0.17em}}={\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{e}^{-i\hslash {k}^{2}t/2m}\frac{dk}{2\pi}\langle x|k\rangle \langle k|{x}_{0}\rangle}\\ ={\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{e}^{-i\hslash {k}^{2}t/2m}\frac{dk}{2\pi}{e}^{-ik\left({x}_{0}-x\right)}}\\ =\sqrt{\frac{m}{2\pi \hslash it}}{e}^{im{\left({x}_{0}-x\right)}^{2}/2\hslash t},\end{array}$$

On examining the above expression, though, it turns out to be nonsense! Noting that the term in the exponent is pure imaginary, ${\left|\psi \left(x,t\right)\right|}^{2}=m/2\pi \hslash t$ independent of $x$! This particle apparently instantaneously fills all of space, but then its probability dies away as $1/t\dots $

*Question*: Where did we go wrong?

*Answer*:
Notice first that ${\left|\psi \left(x,t\right)\right|}^{2}$ is constant throughout space. This means that the normalization, $\int {\left|\psi \left(x,t\right)\right|}^{2}}dx=\infty $! And, as we’ve seen above, the normalization
stays constant in time$\u2014$the
propagator is unitary. Therefore, our
initial wave function must have had infinite norm. That’s exactly right$\u2014$we took
the initial wave function $\psi (x,t=0)=\delta (x-{x}_{0})=|{x}_{0}\rangle $.

Think of the $\delta $ -function as a limit of a function equal to $1/\Delta $ over an interval of length $\Delta $,
with $\Delta $ going to zero, and it’s clear the
normalization goes to infinity as $1/\Delta $. This is not a meaningful wave function for a
particle. Recall that continuum kets
like $|{x}_{0}\rangle $ are normalized by $\langle x|{x}^{\prime}\rangle =\delta \left(x-{x}^{\prime}\right)$,
they do *not *represent wave functions
individually normalizable in the usual sense.
The only meaningful wave functions are* integrals over a range* of such kets, such as $\int dx\psi \left(x\right)}|x\rangle $. In an integral like this, notice that states $|x\rangle $ within some tiny $x$ -interval of length $\delta x,$ say, have total weight $\psi \left(x\right)\delta x$,
which goes to zero as $\delta x$ is made smaller, but by writing $\psi (x,t=0)=\delta (x-{x}_{0})=|{x}_{0}\rangle $ we took a single such state and gave it a
finite weight. This we can’t do.

Of course, we do want to know how a wave function
initially localized *near* a point
develops. To find out, we must apply the
propagator to a legitimate wave function$\u2014$one that
is normalizable to begin with. The simplest “localized particle” wave function
from a practical point of view is a Gaussian wave packet,

$$\psi \left({x}^{\prime},0\right)={e}^{i{p}_{0}{x}^{\prime}/\hslash}\frac{{e}^{-{{x}^{\prime}}^{2}/2{d}^{2}}}{{\left(\pi {d}^{2}\right)}^{1/4}}.$$

(I’ve used $d$ in place of Shankar’s $\Delta $ here to try to minimize confusion with $\Delta \text{\hspace{0.05em}}x,$ etc.)

The wave function at a later time is then given by the operation of the propagator on this initial wave function:

$$\begin{array}{c}\psi (x,t)={\displaystyle \int U(x,t;{x}^{\prime},0){e}^{i{p}_{0}{x}^{\prime}/\hslash}\frac{{e}^{-{{x}^{\prime}}^{2}/2{d}^{2}}}{{\left(\pi {d}^{2}\right)}^{1/4}}d{x}^{\prime}}\\ =\sqrt{\frac{m}{2\pi \hslash it}}{\displaystyle \int {e}^{im{(x-{x}^{\prime})}^{2}/2\hslash t}}{e}^{i{p}_{0}{x}^{\prime}/\hslash}\frac{{e}^{-{{x}^{\prime}}^{2}/2{d}^{2}}}{{\left(\pi {d}^{2}\right)}^{1/4}}d{x}^{\prime}.\end{array}$$

The integral over ${x}^{\prime}$ is just another Gaussian integral, so we use the same result,

$\underset{-\infty}{\overset{\infty}{\int}}d{x}^{\prime}{e}^{-a{{x}^{\prime}}^{2}+b{x}^{\prime}}=\sqrt{\frac{\pi}{a}}{e}^{{b}^{2}/4a}$.

Looking at the expression above, we can see that

$$b=-\frac{im}{\hslash t}\cdot \left(x-\frac{{p}_{0}t}{m}\right),$$ $$a=\frac{1}{2{d}^{2}}-\frac{im}{2\hslash t}.$$

This gives

$$\psi \left(x,t\right)=\frac{{\pi}^{-1/4}}{\sqrt{d\left(1+\frac{i\hslash t}{m{d}^{2}}\right)}}\mathrm{exp}\left(\frac{im{x}^{2}}{2\hslash t}\right)\mathrm{exp}\left(-\frac{\frac{im}{\hslash t}{\left(x-\frac{{p}_{0}t}{m}\right)}^{2}}{2\left(1+\frac{i\hslash t}{m{d}^{2}}\right)}\right)$$

where the second exponential is the term ${e}^{{b}^{2}/4a}$. As written, the small $t$ limit is not very apparent, but some algebraic rearrangement yields:

$$\psi \left(x,t\right)=\frac{{\pi}^{-1/4}}{\sqrt{d\left(1+i\hslash t/m{d}^{2}\right)}}\mathrm{exp}\left(-\frac{{\left(x-{p}_{0}t/m\right)}^{2}}{2{d}^{2}\left(1+i\hslash t/m{d}^{2}\right)}\right)\mathrm{exp}\left(\frac{i{p}_{0}}{\hslash}\left(x-{p}_{0}t/2m\right)\right).$$

It is clear that this expression goes to the initial wave packet as $t$ goes to zero. Although the phase has contributions from all three terms here, the main phase oscillation is in the third term, and one can see the phase velocity is one-half the group velocity, as discussed earlier.

The resulting *probability
density*:

$${\left|\psi \left(x,t\right)\right|}^{2}=\frac{1}{\sqrt{\pi \left({d}^{2}+{\hslash}^{2}{t}^{2}/{m}^{2}{d}^{2}\right)}}\cdot {\mathrm{exp}}^{-\frac{{\left(x-{p}_{0}t/m\right)}^{2}}{\left({d}^{2}+{\hslash}^{2}{t}^{2}/{m}^{2}{d}^{2}\right)}}.$$

This is a Gaussian wave packet, having a width
which goes as $\hslash t/md$ for large times, where $d$ is the width of the initial packet in $x$-space$\u2014$so $\hslash /md$ is the spread in *velocities *$\langle \Delta v\rangle $ within the packet, hence the gradual
spreading $\langle \Delta v\rangle t$ in $x$-space.

It’s amusing to look at the limit of this as the width $d$ of the initial Gaussian packet goes to zero, and see how that relates to our $\delta $-function result. Suppose we are at distance $x$ from the origin, and there is initially a Gaussian wave packet centered at the origin, width $d\ll x.$ At time $t\sim mxd/\hslash $, the wave packet has spread to $x$ and has ${\left|\psi \left(x,t\right)\right|}^{2}$ of order $1/x$ at $x.$ Thereafter, it continues to spread at a linear rate in time, so locally ${\left|\psi \left(x,t\right)\right|}^{2}$ must decrease as $1/t$ to conserve probability. In the $\delta $-function limit $d\to 0$, the wave function instantly spreads through a huge volume, but then goes as $1/t$ as it spreads into an even huger volume. Or something.

### Schrödinger and Heisenberg Representations

Assuming a Hamiltonian with no explicit time dependence, the time-dependent Schrödinger equation has the form

$i\hslash \frac{\partial}{\partial t}|\psi \left(x,t\right)\rangle =H|\psi \left(x,t\right)\rangle $

and as discussed above, the formal solution can be expressed as:

$|\psi \left(x,t\right)\rangle ={e}^{-iHt/\hslash}|\psi \left(x,t=0\right)\rangle .$

Now, any measurement on a system amounts to measuring a matrix element of an operator between two states (or, more generally, a function of such matrix elements).

In other words, the *physically significant time
dependent quantities* are of the form

$\langle \phi \left(t\right)|A|\psi \left(t\right)\rangle =\langle \phi \left(0\right)|{e}^{iHt/\hslash}A{e}^{-iHt/\hslash}|\psi \left(0\right)\rangle $

where $A$ is an operator, which we are assuming has no explicit time dependence.

So in this Schrödinger picture, the time dependence of the measured value of an operator like $x$ or $p$ comes about because we measure the matrix element of an unchanging operator between bras and kets that are changing in time.

Heisenberg took a different approach: he assumed
that the ket describing a quantum system did *not* change in time, it remained at $|\psi \left(0\right)\rangle ,$ but the *operators* evolved according to:

${A}_{H}(t)={e}^{iHt/\hslash}{A}_{H}(0){e}^{-iHt/\hslash}.$

Clearly, this leads to the same physics as before. The equation of motion of the operator is:

$i\hslash \frac{d{A}_{H}(t)}{dt}=\left[{A}_{H}(t),H\right].$

The Hamiltonian itself does not change in time$\u2014$energy is conserved, or, to put it another way, $H$ commutes with ${e}^{-iHt/\hslash}.$ But for a nontrivial Hamiltonian, say for a particle in one dimension in a potential,

$H={p}^{2}/2m+V(x)$

the separate components *will* have
time-dependence, parallel to the classical case: the *kinetic* energy of a swinging pendulum varies with time. (For a particle in a potential in an energy
eigenstate the expectation value of the kinetic energy is constant, but this is
not the case for any other state, that is, for a superposition of different
eigenstates.) Nevertheless, the *commutator* of $x,p$ will be time-independent:

$\begin{array}{c}\left[{x}_{H}(t),{p}_{H}(t)\right]={e}^{iHt/\hslash}\left[{x}_{H}(0),{p}_{H}(0)\right]{e}^{-iHt/\hslash}\\ ={e}^{iHt/\hslash}i\hslash {e}^{-iHt/\hslash}\\ =i\hslash .\end{array}$

(The Heisenberg operators are identical to the Schrödinger operators at $t=0.$ )

Applying the general commutator result $\left[A,BC\right]=\left[A,B\right]C+B\left[A,C\right]$,

$$\left[{x}_{H}(t),\text{\hspace{0.17em}}\frac{{p}^{2}{}_{H}(t)}{2m}\right]=\frac{i\hslash {p}_{H}(t)}{m}$$

so

$$\frac{d{x}_{H}(t)}{dt}=\frac{{p}_{H}(t)}{m}$$

and since $\left[{x}_{H}(t),{p}_{H}(t)\right]=i\hslash ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{p}_{H}(t)=-i\hslash d/d{x}_{H}(t)$,

$$\frac{d{p}_{H}(t)}{dt}=\text{\hspace{0.17em}}\frac{1}{i\hslash}\left[{p}_{H}(t),V({x}_{H}(t))\right]=-\nabla V({x}_{H}(t)).$$

This result could also be derived by writing $V\left(x\right)$ as an expansion in powers of $x,$ then taking the commutator with $p.$

*Exercise*: check this.

Notice from the above equations that the operators
in the Heisenberg Representation obey the *classical
*laws of motion! *Ehrenfest’s
Theorem*, that the *expectation values* of operators in a quantum state
follow the classical laws of motion, follows immediately, by taking the
expectation value of both sides of the operator equation of motion in a quantum
state.

### Simple Harmonic Oscillator in the Heisenberg Representation

For the simple harmonic oscillator, the equations are easily integrated to give:

$\begin{array}{l}{x}_{H}(t)={x}_{H}(0)\mathrm{cos}\omega t+({p}_{H}(0)/m\omega )\mathrm{sin}\omega t\\ {p}_{H}(t)={p}_{H}(0)\mathrm{cos}\omega t-m\omega {x}_{H}(0)\mathrm{sin}\omega t.\end{array}$

We have put in the ${}_{H}$ subscript to emphasize that these are operators. It is usually clear from the context that the Heisenberg representation is being used, and this subscript may be safely omitted.

The time-dependence of the *annihilation operator *$a$ is:

$a\left(t\right)={e}^{iHt/\hslash}a\left(0\right){e}^{-iHt/\hslash}$

with

$H=\hslash \omega \left({a}^{\u2020}\left(t\right)a\left(t\right)+{\scriptscriptstyle \raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}\right).$

Note again that although $H$ is itself time-independent, it is necessary to include the time-dependence of individual operators within $H.$

$\begin{array}{c}i\hslash \frac{d}{dt}a\left(t\right)=\left[a\left(t\right),H\right]\\ =\hslash \omega \left[a\left(t\right),{a}^{\u2020}\left(t\right)a\left(t\right)\right]\\ =\hslash \omega \left[a\left(t\right),{a}^{\u2020}\left(t\right)\right]a\left(t\right)\\ =\hslash \omega \text{\hspace{0.05em}}\text{\hspace{0.05em}}a\left(t\right)\end{array}$

so

$a\left(t\right)=a\left(0\right){e}^{-i\omega t}.$

Actually, we could have seen this as follows: if $|n\rangle $ are the energy eigenstates of the simple harmonic oscillator,

${e}^{-iHt/\hslash}|n\rangle ={e}^{-in\text{\hspace{0.05em}}\hslash \omega t/\hslash}|n\rangle ={e}^{-in\omega t}|n\rangle .$

Now the *only*
nonzero matrix elements of the annihilation operator $\widehat{a}$ between energy eigenstates are of the form

$\begin{array}{c}\langle n-1|a\left(t\right)|n\rangle =\langle n-1|{e}^{iHt/\hslash}a\left(0\right){e}^{-iHt/\hslash}|n\rangle \\ ={e}^{i\omega \left(n-1\right)t}\langle n-1|a\left(0\right)|n\rangle {e}^{-i\omega nt}\\ =\langle n-1|a\left(0\right)|n\rangle {e}^{-i\omega t}.\end{array}$

Since this time-dependence is true of *all* energy matrix elements (trivially so
for most of them, since they’re identically zero), and the eigenstates of the
Hamiltonian span the space, it is true as an *operator* equation.

Evidently, the expectation value of the operator $a\left(t\right)$ in any state goes clockwise in a circle centered at the origin in the complex plane. That this is indeed the classical motion of the simple harmonic oscillator is confirmed by recalling the definition $a=\frac{\xi +i\pi}{\sqrt{2}}=\frac{1}{\sqrt{2\hslash m\omega}}\left(m\omega x+ip\right)$, so the complex plane corresponds to the $\left(m\omega x,p\right)$ phase space discussed near the beginning of the lecture on the Simple Harmonic Oscillator. We’ll discuss this in much more detail in the next lecture, on Coherent States.

The time-dependence of the creation operator is just the adjoint equation: ${a}^{\u2020}\left(t\right)={a}^{\u2020}\left(0\right){e}^{i\omega t}.$