# 3.5: Propagators and Representations

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 wavefunction 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ℏ \dfrac{∂}{\partial t}ψ(x,t) =−\dfrac{ℏ^2}{2m} \frac{∂^2}{∂x^2} ψ(x,t) +V(x)ψ(x,t),\]

Or, as we would now write it

\[iℏ\dfrac{∂}{∂t}|ψ(x,t)⟩=H|ψ(x,t)⟩.\]

Since \( H\) is itself time independent, this is very easy to integrate!

\[|ψ(x,t)⟩=e^{−iH(t−t_0)/ℏ} |ψ(x,t_0)⟩.\]

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\):

\[|ψ(x,t)⟩=U(t−t_0)|ψ(x,t_0) ⟩.\]

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

\[U(t−t_0)=e^{−iH}(t−t_0) \]

so

\[ U^†(t−t_0)=e^{iH^†}(t−t_0)=e^{iH}(t−t_0)=U−1(t−t_0).\]

Since H is Hermitian, \(U\) is unitary. It immediately follows that

\[⟨ψ(x,t)|ψ(x,t)⟩=⟨ψ(x,t_0)|U^† U(t−t_0 )|ψ(x,t_0)⟩=⟨ψ(x,t_0)|ψ(x,t_0)⟩\]

the norm of the ket vector is conserved, or, translating to wavefunction language, a wavefunction 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(t−t_0)/ℏ}\) are the same as the eigenstates of \(H\), and if \(H|ψ_n⟩=E_n |ψ_n⟩\), then

\[e^{−iH(t−t_0)/ℏ} | ψ_n ⟩ =e^{−iEn(t−t_0)/ℏ}|ψ_n⟩.\]

This is of course nothing but the time dependent phase factor for the eigenstates we found before — 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:

\[|ψ(t)⟩=e^{−iH(t−t_0)/ℏ} \sum_{n=1}^∞ |ψ_n⟩⟨ψ_n|ψ(t_0)⟩=\sum_{n=1}^∞e^{−iE_n(t−t_0)/ℏ}|ψ_n⟩⟨ψ_n|ψ(t_0)⟩.\]

Staring at this, we see that it’s just what we had before: at the initial time \(t=t_0\), the wavefunction can be written as a sum over the eigenkets:

\[|ψ(t_0)⟩=∑|ψ_n(t_0) ⟩ ⟨ψ_n (t_0)|ψ(t_0)⟩= \sum c_n|ψ_n(t_0)⟩\]

with

\[c_n=⟨ψ_n| ψ⟩\]

and

\[\sum |c_n|^2=1\]

and the usual generalization for continuum eigenvalues, and the time development is just given by inserting the phases:

\[|ψ(t)⟩=\sum c_n e^{−iE_n(t−t_0) /ℏ} |ψ_n(t_0)⟩.\]

The expectation value of the energy \(E\) in \(|ψ⟩\),

\[⟨E⟩=⟨ψ|H|ψ⟩=\sum |c_n|^2E_n\]

and is (of course) time independent.

The expectation value of the particle position \(x\) is

\[⟨ψ(t)|x|ψ(t)⟩=\sum _{n,m} c^∗_n c_m e^{i(E_n−E_m)(t−t_0)/ℏ}⟨ψ_n(t_0)|x|ψ_m(t_0)⟩\]

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=\dfrac{p^2}{2m}=\dfrac{ℏ^2 k^2}{2m}\]

the energy eigenstates are also momentum eigenstates, we label them \(|k⟩\), so

\[U(t)=e^{−iHt/ℏ}=\in_{−∞}^{∞} e^{−iH t/ℏ dk^2π} |k⟩⟨k|= \int _{−∞}^{∞}e^{−iℏ }k^2t/2mdk^2π|k⟩⟨k|.\]

a particle is at \(x_0\): \(ψ(x,t=0)=δ(x−x_0)=|x_0⟩\): what is the probability amplitude for finding it at \(x\) at a later time \(t\)? (This would be just its wavefunction at the later time.)

\[⟨x|U(t,0)|x_0⟩ = \int_{−∞}^{∞} e^{−iℏk^2t/ 2m} dk^2π⟨x |k⟩⟨k|x0⟩=\int _{−∞}^{∞} e^{−iℏ k^2t/2m} dk2πe^{−ik(x_0−x)}=m^2πℏit−−−−√e^{im(x_0−x)2}{2ℏ}t, \label{eq30}\]

Exercise \(\PageIndex{1}\)

On examining Equation \ref{eq30}, though, it turns out to be nonsense! Noting that the term in the exponent is pure imaginary, \(|ψ(x,t)|^2=\dfrac{m}{2πℏt}\) independent of \(x\)! This particle apparently instantaneously fills all of space, but then its probability dies away as 1/t…

Question: Where did we go wrong?

Answer: Notice first that \(|ψ(x,t)|^2\) is constant throughout space. This means that the normalization, ∫|ψ(x,t)|2dx=∞ ! And, as we’ve seen above, the normalization stays constant in time — the propagator is unitary. Therefore, our initial wavefunction must have had infinite norm. That’s exactly right — we took the initial wavefunction

\][ψ(x,t=0)=δ(x−x_0)=| x_0 ⟩.\]

Think of the δ-function as a limit of a function equal to 1/Δ over an interval of length Δ , with Δ going to zero, and it’s clear the normalization goes to infinity as 1/Δ . This is not a meaningful wavefunction for a particle. Recall that continuum kets like |x0⟩ are normalized by ⟨x|x′⟩=δ(x−x′ ) , they do not represent wavefunctions individually normalizable in the usual sense. The only meaningful wavefunctions are integrals over a range of such kets, such as ∫dxψ(x)|x⟩ . In an integral like this, notice that states |x⟩ within some tiny x -interval of length δx, say, have total weight ψ(x)δx , which goes to zero as δx is made smaller, but by writing ψ(x,t=0)=δ(x−x0)=| x0 ⟩ 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 wavefunction initially localized near a point develops. To find out, we must apply the propagator to a legitimate wavefunction — one that is normalizable to begin with. The simplest “localized particle” wavefunction from a practical point of view is a Gaussian wave packet,

\[ψ(x′,0)=e^{ip_0x′/ℏ}e^{−x′2/2d^2(πd^2)}1/4.\]

(I’ve used d in place of Shankar’s Δ here to try to minimize confusion with Δ x, etc.)

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

\[ψ(x,t)=∫U(x,t;x′,0)eip0x′/ℏ e−x′2 /2d2(πd2)1/4dx′=m2πℏit−−−−√∫eim(x−x′)2/2ℏteip0x′/ℏe−x′2/2d2(πd2)1/4dx′.\]

The integral over x′ is just another Gaussian integral, so we use the same result,

\[∫−∞∞dx′e−ax′2 +bx′=πa−−√eb2/4a .\]

Looking at the expression above, we can see that

\[b=−imℏt⋅(x−p0tm ) , a=12d2 −im2ℏt .\]

This gives

\[ψ(x,t)=π−1/4d(1+iℏtmd2 ) √exp(imx22ℏt)exp(−imℏt(x−p0tm)22(1+iℏtmd2))\]

where the second exponential is the term eb2/4a . As written, the small t limit is not very apparent, but some algebraic rearrangement yields:

\[ψ(x,t)=π−1/4d(1+iℏt/md2) √exp(−(x−p0t/m)22d2(1+iℏt/md2))exp(ip0ℏ(x−p0t/2m)) .\]

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:

\[|ψ(x,t)|2=1π(d2+ℏ2 t2/m2d2)√⋅exp−(x−p0t/m)2(d2+ℏ2t2/m2d2) .\]

This is a Gaussian wave packet, having a width which goes as ℏt/md for large times, where d is the width of the initial packet in x -space — so ℏ/md is the spread in velocities ⟨Δv⟩ within the packet, hence the gradual spreading ⟨Δv⟩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 δ -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≪x. At time t∼mxd/ℏ , the wave packet has spread to x and has |ψ(x,t)|2 of order 1/x at x. Thereafter, it continues to spread at a linear rate in time, so locally |ψ(x,t)|2 must decrease as 1/t to conserve probability. In the δ -function limit d→0 , the wavefunction 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ℏ\dfrac{∂}{∂t} |ψ(x,t)⟩=H|ψ(x,t)⟩\]

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

\[|ψ(x,t)⟩=e^{−iHt/ℏ}|ψ(x,t=0)⟩.\]

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

\[⟨φ(t)|A|ψ(t)⟩=⟨φ(0)|e^{iHt/ℏ} Ae^{−iHt/ℏ} |ψ(0)⟩\]

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 |ψ(0)⟩, but the operators evolved according to:

\[AH(t)=e^{iHt/ℏ} AH (0)e^{−iHt/ℏ}.\]

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

\[iℏ\dfrac{d}{dt} AH(t) =[AH (t),H].\]

The Hamiltonian itself does not change in time — energy is conserved, or, to put it another way, H commutes with e−iHt/ℏ. But for a nontrivial Hamiltonian, say for a particle in one dimension in a potential,

\[H=\dfrac{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\), and \(p\) will be time-independent:

\[[xH(t),pH(t)]=e^{iHt/ℏ} [xH(0),pH(0)]e^{−iHt/ℏ}=e^{iHt/ℏ}iℏe^{−iHt/ℏ}=iℏ.\]

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

Applying the general commutator result [A,BC]=[A,B]C+B[A,C] ,

\[[xH(t), p2H(t) 2m]=iℏpH(t)m\]

so

\[dxH(t)dt =pH (t)m\]

and since

\[[xH(t),pH(t)] =iℏ, pH (t)=−iℏd/dxH(t) ,\]

\[dpH(t)dt = 1iℏ [pH(t),V(xH(t))]=−∇V(xH(t)).\]

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

Exercise \(\PageIndex{1}\)

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:

\[xH(t)=xH(0)cosωt+(pH (0)/mω)sinωtpH(t)=pH(0)cosωt−mωxH(0)sinωt.\]

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(t)=eiHt/ℏa(0)e−iHt/ℏ\]

with

\[H=ℏω(a†(t)a(t)+12/ ).\]

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

\[iℏddta(t)=[a(t),H]=ℏω[a(t),a†(t)a(t)] =ℏω[a(t),a†(t)]a(t)=ℏω a(t)\]

so

\[a(t)=a(0)e−iωt.\]

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

\[e−iHt/ℏ|n⟩=e−in ℏωt/ℏ |n⟩=e−inωt|n⟩.\]

Now the only nonzero matrix elements of the annihilation operator aˆ between energy eigenstates are of the form

\[⟨n−1|a(t)|n⟩=⟨n−1|eiHt/ℏa(0)e−iHt/ℏ |n⟩=eiω(n−1)t⟨n−1|a(0)|n⟩e−iωnt=⟨n−1|a(0)|n⟩e−iωt.\]

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(t) 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=ξ+iπ2√ =12ℏmω√(mωx+ip) , so the complex plane corresponds to the (mωx,p) 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†(t)=a† (0)eiωt .\]

### Contributors

Michael Fowler (Beams Professor, Department of Physics, University of Virginia)