Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

04. Spatial Degrees of Freedom

Introduction

In these notes we develop the theory of wave functions in configuration space, building it up from the ket formalism and the postulates of quantum mechanics. We consider a system with one particle moving in space. In this set of notes we will not deal with the dynamics of the particle, so we will not need to say what the Hamiltonian is (that is taken care of in 05. Time Evolution in Quantum Mechanics). We assume the particle has no spin or other internal degrees of freedom, or that if such degrees of freedom exist they can be ignored. This allows us to establish an isomorphism between the Hilbert space of kets and the Hilbert space of wave functions on configuration space. We then introduce translation operators and use them and the classical correspondence to motivate the definition of momentum in quantum mechanics. We then explore the two representations, position and momentum, and the relationship between them. Finally, we discuss minimum uncertainty wave packets.

The Position Representation; Wave Functions

For simplicity we begin with the one-dimensional case, in which the physical system is a particle moving in one dimension. According to the postulates of quantum mechanics, this system is associated with a Hilbert space E\ES. We assume the position xx of a particle in one dimension can be measured, and that the results of the measurement are continuous. Thus, we are dealing with the case of the continuous spectrum. We denote the operator corresponding to measuring xx by x^\xhat, so the eigenvalue-eigenket problem is

x^∣x⟩=x∣x⟩.\xhat \ket{x} = x \ket{x}.

Here x^\xhat is the operator, xx is the eigenvalue, and ∣x⟩\ket{x} is the eigenket with eigenvalue xx. Since xx belongs to the continuous spectrum, the eigenkets ∣x⟩\ket{x} are not normalizable and lie outside Hilbert space.

We also assume that x^\xhat by itself forms a complete set of commuting observables, which, as we shall see, means that the wave function is a scalar (it has only one component). In practice, this means that we are dealing with a spinless particle (a particle of spin 0, for example, the π\pi-meson, the nucleus of ordinary helium, or the hydrogen atom in its ground state). It may also mean that we have a particle with spin, such as the electron, but the spin degrees of freedom are not important for the physics we are interested in. This point was discussed in Resolving Degeneracies; Complete Sets of Commuting Observables.

Since x^\xhat is a complete set of commuting observables, the eigenkets ∣x⟩\ket{x} are nondegenerate and form a basis in the Hilbert space in the continuum sense. These eigenkets are only defined by (1) to within a normalization and a phase. We fix the normalization by requiring

⟨x1∣x2⟩=δ(x1−x2),\braket{x_1}{x_2} = \delta(x_1-x_2),

the usual orthonormality condition in the continuum. The phase conventions for the eigenkets ∣x⟩\ket{x} are a question we shall return to later. Under these assumptions, the resolution of the identity takes the form

1=∫−∞+∞dx ∣x⟩⟨x∣.1 = \int_{-\infty}^{+\infty} dx\, \ketbra{x}{x}.

Now let ∣ψ⟩\ket{\psi} represent a normalized (pure) state of the system, let I=[x0,x1]I=[x_0,x_1] be an interval on the xx-axis, and let PIP_I be the corresponding projection operator,

PI=∫x0x1dx ∣x⟩⟨x∣.P_I= \int_{x_0}^{x_1} dx \, \ketbra{x}{x}.

Then, according to the postulates of quantum mechanics (Postulates of Quantum Mechanics (Incomplete Form)), we have

textProb(x0≤x≤x1)=⟨ψ∣PI∣ψ⟩=∫x0x1dx ⟨ψ∣x⟩⟨x∣ψ⟩.\\textProb(x_0 \le x \le x_1) = \matrixelement{\psi}{P_I}{\psi} = \int_{x_0}^{x_1} dx \, \braket{\psi}{x} \braket{x}{\psi}.

We now introduce the definition,

ψ(x)=⟨x∣ψ⟩,\boxed{\psi(x) = \braket{x}{\psi},}

of the wave function ψ(x)\psi(x), whereupon the probability (5) can be written,

textProb(x0≤x≤x1)=∫x0x1dx ∣ψ(x)∣2.\\textProb(x_0 \le x \le x_1) = \int_{x_0}^{x_1} dx \, |\psi(x)|^2.

Since the interval [x0,x1][x_0,x_1] is arbitrary, we see that ∣ψ(x)∣2|\psi(x)|^2 must be interpreted as the probability density of finding the particle on the xx-axis. We have derived wave functions from the ket formalism plus the postulates of quantum mechanics, as promised in 01. The Mathematical Formalism of Quantum Mechanics.

(6) is important because it shows how to go from kets to wave functions. It takes the place of sloppy notation such as (7). That is, in addition to the Hilbert space E\ES of kets, we have another Hilbert space, this one consisting of normalizable wave functions on configuration space, ψ(x)\psi(x). (6) expresses the mapping between these spaces, which identifies kets with wave functions. To go the other direction, we multiply ∣ψ⟩\ket{\psi} on the left by the resolution of the identity (3),

∣ψ⟩=∫dx ∣x⟩⟨x∣ψ⟩,\ket{\psi} = \int dx \, \ket{x} \braket{x}{\psi},

which, by using (6), becomes

∣ψ⟩=∫dx ∣x⟩ψ(x).\boxed{\ket{\psi} = \int dx \, \ket{x} \psi(x).}

This is the inverse of (6), allowing one to go from a wave function ψ(x)\psi(x) to the corresponding ket ∣ψ⟩\ket{\psi}. We see that ψ(x)\psi(x) is the set of expansion coefficients of the ket ∣ψ⟩\ket{\psi} in the position eigenbasis.

Another consequence of this formalism follows easily. Suppose we have two kets ∣ψ⟩\ket{\psi} and ∣ϕ⟩\ket{\phi}, and we wish to compute the scalar product ⟨ψ∣ϕ⟩\braket{\psi}{\phi} in wave function language. We simply insert the resolution of the identity (3) between the bra and the ket to obtain,

⟨ψ∣ϕ⟩=∫dx ⟨ψ∣x⟩⟨x∣ϕ⟩=∫dx ψ†(x)ϕ(x),\braket{\psi}{\phi} = \int dx \, \braket{\psi}{x}\braket{x}{\phi} = \int dx \, \psi^\cc(x)\phi(x),

a familiar formula that translates the scalar product on the Hilbert space of kets into the language of wave functions.

This formalism is easily generalized to three dimensions. Let x\xvec be a position vector in three-dimensional space, with components (x,y,z)(x,y,z) or (x1,x2,x3)(x_1,x_2,x_3). The measurement of the three components of position corresponds to a vector of operators, which we denote by x^{\hat{\mathbf{x}}} (with a hat), to distinguish the operators from their eigenvalues (the results of the position measurement). In quantum mechanics, when we speak of a “vector operator” we usually mean a vector of operators, in this case three operators x^=(x^,y^,z^)=(x^1,x^2,x^3){\hat{\mathbf{x}}}= (\xhat,\yhat,\zhat) = (\xhat_1, \xhat_2, \xhat_3). We assume that the components of x^{\hat{\mathbf{x}}} commute with one another,

[x^i,x^j]=0,[\xhat_i,\xhat_j]=0,

which as explained in 02. The Postulates of Quantum Mechanics\ can be tested experimentally. For example, if we filter a beam in xx and then yy by means of small slits, or do it in the reverse order, we find that the statistical results of arbitrary measurements on the filtered beam are the same in both cases.

As above we assume that we can ignore any internal degrees of freedom of the particle, so that the three operators x^{\hat{\mathbf{x}}} by themselves form a complete set of commuting observables. Then the simultaneous eigenkets of x^{\hat{\mathbf{x}}} are nondegenerate and form a basis by themselves. The eigenket-eigenvalue problem is

x^i∣x⟩=xi∣x⟩,\xhat_i\ket{\xvec} = x_i\ket{\xvec},

for i=1,2,3i=1,2,3, a set of three simultaneous equations satisfied by an single eigenket ∣x⟩=∣x,y,z⟩\ket{\xvec}=\ket{x,y,z}, labeled by the three eigenvalues of the three commuting operators. In obvious generalizations of the one-dimensional case, we normalize these kets according to

⟨x1∣x2⟩=δ3(x1−x2)=δ(x1−x2)δ(y1−y2)δ(z1−z2),\braket{\xvec_1}{\xvec_2} = \delta^3(\xvec_1-\xvec_2) = \delta(x_1-x_2) \delta(y_1-y_2) \delta(z_1-z_2),

we have a resolution of the identity,

1=∫d3x ∣x⟩⟨x∣,1 = \int d^3\xvec \, \ketbra{\xvec}{\xvec},

and the transformation between wave functions and kets is given by

ψ(x)=⟨x∣ψ⟩,\psi(\xvec) = \braket{\xvec}{\psi},

and

∣ψ⟩=∫d3x ∣x⟩ψ(x).\ket{\psi} = \int d^3\xvec \, \ket{\xvec} \psi(\xvec).

Also, the probability that the particle will be found in a region RR of three-dimensional space is

∫Rd3x ∣ψ(x)∣2,\int_R d^3\xvec \, |\psi(\xvec)|^2,

so ∣ψ(x)∣2|\psi(\xvec)|^2 must be interpreted as the probability density.

Equations (6) or (15) define the wave function ψ(x)\psi(x) or ψ(x)\psi(\xvec), respectively, in quantum mechanics. A question that arises sometimes when solving the Schr"odinger equation is whether ψ\psi is single-valued. But by our definitions, ψ(x)\psi(x) or ψ(x)\psi(\xvec) is automatically single-valued, since it is just the expansion coefficient of the state ∣ψ⟩\ket{\psi} with respect to the basis ∣x⟩\ket{x} or ∣x⟩\ket{\xvec}.

Equations~(6) and (9) (in one dimension) or (15) and (16) (in three dimensions) make explicit the transformations back and forth between the Hilbert space of ket vectors ∣ψ⟩\ket{\psi} and the Hilbert space of wave functions ψ(x)\psi(x) or ψ(x)\psi(\xvec). As discussed in Dirac Ket Notation, we regard these two Hilbert spaces as isomorphic but distinct, and we say that the space of wave functions ψ(x)\psi(x) or ψ(x)\psi(\xvec) forms the configuration representation of the space of kets ∣ψ⟩\ket{\psi}. There are other representations, and while the configuration representation is often used in practice one should note that in the postulates of quantum mechanics there is no privileged role assigned to it. Any calculation that can be carried out in the configuration representation can be carried out in any other representation with the same results from a physical standpoint.

In effect, a representation is merely a choice of basis, allowing one to work with the expansion coefficients of a state vector ∣ψ⟩\ket{\psi} with respect to the chosen basis instead of the abstract ket vectors themselves. As noted below (9), the configuration space wave function ψ(x)\psi(x) or ψ(x)\psi(\xvec) is the expansion coefficients of the state ket ∣ψ⟩\ket{\psi} with respect to the basis of position eigenkets {∣x⟩}\{ \ket{x} \} or {∣x⟩}\{ \ket{\xvec} \}. In quantum mechanics we frequently choose as a basis the simultaneous eigenkets of a complete set of commuting observables; thus, the basis is labeled by the observables in question, and we speak, for example, of the position representation, momentum representation, etc. Nothing prevents us from using other bases, however (which have no correspondence with any simple set of commuting observables).

Just as the wave function ψ(x)\psi(x) represents the ket vector ∣ψ⟩\ket{\psi} in the configuration representation, so also is there a representation of various operators that act on ket vectors. Consider, for example, the operator x^\xhat (working in one dimension for simplicity). We know that the ket ∣ψ⟩\ket{\psi} corresponds to the wave function ψ(x)=⟨x∣ψ⟩\psi(x)=\braket{x}{\psi}. What wave function does the ket x^∣ψ⟩\xhat\ket{\psi} correspond to? Let ∣ϕ⟩=x^∣ψ⟩\ket{\phi} = \xhat\ket{\psi}, so that

ϕ(x)=⟨x∣ϕ⟩=⟨x∣x^∣ψ⟩=x⟨x∣ψ⟩=xψ(x),\phi(x) = \braket{x}{\phi} = \matrixelement{x}{\xhat}{\psi} = x \braket{x}{\psi} = x\psi(x),

where we have used ⟨x∣x^=x⟨x∣\bra{x}\xhat = x \bra{x}, the Hermitian conjugate of (1). We see that the effect of multiplying a ket vector ∣ψ⟩\ket{\psi} by the operator x^\xhat is to multiply the corresponding configuration wave function ψ(x)\psi(x) by xx. That is, we can write

(x^ψ)(x)=xψ(x).(\xhat \psi)(x)=x\psi(x).

As we say, the operator x^\xhat is represented by multiplication by xx in the configuration representation. It is represented by other operations in other representations. Similarly, in three dimensions we have

(x^iψ)(x)=xiψ(x).(\xhat_i \psi)(\xvec) = x_i\psi(\xvec).

You may worry about the physical meaning of the unnormalizable eigenkets ∣x⟩\ket{x} or ∣x⟩\ket{\xvec}. In reality we never measure the position of a particle exactly, instead the best we can do is to localize it in some small region of space. The eigenkets ∣x⟩\ket{x} or ∣x⟩\ket{\xvec} are an idealization of this process, in which the size of the region is allowed to approach zero. This limit also leads to singular mathematics, since the eigenkets ∣x⟩\ket{x} or ∣x⟩\ket{\xvec} have infinite norm and do not belong to Hilbert space.

Measuring the position of a particle is really a nonrelativistic concept, because as we localize a particle to smaller and smaller regions, by the uncertainty principle the momentum increases, ultimately taking on relativistic values. Then processes such as the creation of particle-antiparticle pairs come into play, and we are really dealing with a multi-particle situation, which is properly handled by the methods of quantum field theory. We will find later in the course that the position operator for a particle in relativistic quantum mechanics is one that is fraught with difficulties. The position operator is really a nonrelativistic concept.

Translation Operators

The subject of position measurements leads naturally into the question of translational invariance. The energy and other physical properties of an isolated system do not change if it is moved somewhere else, a fact that is more precisely stated by saying that the Lagrangian or Hamiltonian of such systems is invariant under spatial translations. This symmetry in turn leads to the conservation of momentum, which first appeared in physics as Newton’s third law.

If a system is not isolated then its energy does depend on its position, for example, the potential energy of an electron in an electric field changes if the electron is moved. But the electric field is created by other charges somewhere else, and if those charges are included in the “system,” then the energy does not change when the entire system is moved.

In summary, spatial translations are a fundamental symmetry of nature which is closely related to the law of conservation of momentum. We now develop translation operators in quantum mechanics, working initially in one dimension for simplicity.

Let aa be a displacement. We imagine a displacement operation as one that acts on a physical system, moving all particles from their initial positions (say, xx) to their new positions (x+ax+a), both measured relative to an inertial frame. This is sometimes called the active point of view, because our operations take a given physical system and transform it into a new system, in this case in a different location.

In quantum mechanics, we can define a translation operator which carries out this operation on a physical system. The translation operator T(a)T(a) is a linear operator acting on the Hilbert space of a physical system, parameterized by the displacement aa. We define the translation operator in one dimension by

T(a)∣x⟩=∣x+a⟩.T(a)\ket{x} = \ket{x+a}.

This definition makes sense, because physically ∣x⟩\ket{x} is the state of the system after a measurement has placed the particle in a small region around position xx, and similarly for ∣x+a⟩\ket{x+a}. (21) is actually a definition of the operators T(a)T(a) because the kets ∣x⟩\ket{x} form a basis. If we specify the action of a linear operator on a set of basis vectors, then by linear superposition its action on an arbitrary vector becomes known.

(21) gives the action of the translation operator on the position eigenkets. Let us also work out its action on wave functions. Let ∣ψ⟩\ket{\psi} be a state with wave function ψ(x)=⟨ψ∣x⟩\psi(x) =\braket{\psi}{x}, and let ∣ϕ⟩=T(a)∣ψ⟩\ket{\phi}=T(a)\ket{\psi} be the translated state with wave function ϕ(x)=⟨x∣ϕ⟩\phi(x) =\braket{x}{\phi}. What is the relationship between the old wave function ψ(x)\psi(x) and the new one ϕ(x)\phi(x)? We answer this by writing

ϕ(x)=⟨x∣ϕ⟩=⟨x∣T(a)∣ψ⟩=∫dx′ ⟨x∣T(a)∣x′⟩⟨x′∣ψ⟩=∫dx′ ⟨x∣x′+a⟩⟨x′∣ψ⟩\begin{aligned} \phi(x) &= \braket{x}{\phi} =\matrixelement{x}{T(a)}{\psi} = \int dx' \, \matrixelement{x}{T(a)}{x'} \braket{x'}{\psi} =\int dx' \, \braket{x}{x'+a} \braket{x'}{\psi} \end{aligned}
=∫dx′ δ(x−x′−a)ψ(x′)=ψ(x−a),\begin{aligned} &= \int dx' \, \delta(x-x'-a) \psi(x') = \psi(x-a), \end{aligned}

where we have inserted a resolution of the identity, used (21) and (2), and carried out the integral. We write the result as

(T(a)ψ)(x)=ψ(x−a),\bigl(T(a)\psi\bigr)(x) = \psi(x-a),

where we have replaced ϕ\phi by T(a)ψT(a)\psi. This is a companion to (21), which gives the action of translation operators on the basis kets; this gives their action on wave functions.

There are several remarks concerning (24). First, notice the minus sign in this equation, compared to the plus sign in (21). The minus sign is necessary to get a wave function that has been moved forward under the displacement operation, as illustrated in Figure 1. To remember the signs it helps to write (23) in the form, (T(a)ψ)(x+a)=ψ(x)\bigl(T(a)\psi\bigr) (x+a) =\psi(x) and to say, “the value of the new wave function at the new point equals the value of the old wave function at the old point.”

The action of the translation operator T(a) on a wave function \psi(x).

Figure 1:The action of the translation operator T(a)T(a) on a wave function ψ(x)\psi(x).

Another remark is that that (24) uses the translation operator T(a)T(a) in a different sense from its original definition, because it is acting on a configuration space wave function ψ\psi instead of a ket ∣ψ⟩\ket{\psi}. As we say, (24) gives the representation of the operator T(a)T(a) on configuration space wave functions.

Finally, we remark that many books would write (24) without the parentheses, that is, as

T(a)ψ(x)=ψ(x−a).T(a)\psi(x) = \psi(x-a).

The problem with this notation is that it is not clear what T(a)T(a) acts on. ψ\psi is a function, and ψ(x)\psi(x) is the value of that function at a point xx, that is, it is a number. Does T(a)T(a) act on the function or the value of the function? Obviously, it acts on the function, which is what the extra parentheses in (24) make explicit.

The translation operator is easily generalized to three dimensions, where the displacement a\avec is a vector, and the translation operator is defined by

T(a)∣x⟩=∣x+a⟩.T(\avec)\ket{\xvec} = \ket{\xvec+\avec}.

The other formulas and results of this are easily generalized to the three-dimensional case.

Properties of the Translation Operators

The translation operators in one dimension satisfy the following properties:

T(0)=1,\begin{aligned} T(0) &= 1, \end{aligned}
T(a)T(b)=T(a+b)=T(b)T(a),\begin{aligned} T(a)T(b) &= T(a+b) = T(b)T(a), \end{aligned}
T(a)−1=T(−a),\begin{aligned} T(a)^{-1} &= T(-a), \end{aligned}
T(a)−1=T(a)†.\begin{aligned} T(a)^{-1} &= T(a)^\dagger. \end{aligned}

Property (27) follows immediately from the definition (21). Property (28) is the composition law that is also obvious from the definition (21); it can be proved explicitly by writing

T(a)T(b)∣x⟩=T(a)∣x+b⟩=∣x+b+a⟩=T(a+b)∣x⟩.T(a)T(b)\ket{x} = T(a)\ket{x+b} = \ket{x+b+a} = T(a+b)\ket{x}.

Note that the parameter x+b+ax+b+a of the ket ∣x+b+a⟩\ket{x+b+a} is just a label of the ket, so these labels obey the usual (commutative) rules for addition. This implies that the translation operators are commutative, as indicated in (28). The third property (29) is proved by setting b=−ab=-a in (28). The translation operators are invertible. The last property (30) states that the translation operators are unitary. We prove this by letting ψ(x)\psi(x) be an arbitrary wave function, and ∣ϕ⟩=T(a)∣ψ⟩\ket{\phi} = T(a)\ket{\psi}, so that ϕ(x)=ψ(x−a)\phi(x)=\psi(x-a). Then

∫dx ∣ϕ(x)∣2=∫dx ∣ψ(x−a)∣2=∫dx ∣ψ(x)∣2,\int dx \, |\phi(x)|^2 = \int dx\, |\psi(x-a)|^2= \int dx \, |\psi(x)|^2,

as we prove by substituting x′=x−ax'=x-a in the middle integral. Thus, T(a)T(a) preserves the norm of arbitrary states. But a linear operator that does this is necessarily unitary (see Prob. Paragraph(c)).

Altogether, properties (27)--(30) qualify the set of translation operators as a group of unitary operators. Unitary operators appear frequently in symmetry operations because these are the only linear operators that preserve probabilities (hence, the results of physical measurements). In this case, the group is Abelian, which means that the translation operators commute with one another (property (28)).

The Generator of Translations

A simple idea that arises in the case of continuous symmetries is that a given, finite symmetry operation can be built up as a composition of smaller symmetry operations. For example, a displacement of one meter is the composition of a thousand displacements of one millimeter. In the limit we can imagine a finite symmetry operation as being built up out of an infinite number of infinitesimal symmetry operations. For this reason, special attention is attached to infinitesimal symmetry operations. It turns out that infinitesimal versions of unitary symmetry operators in quantum mechanics are always expressible in terms of certain Hermitian operators, which are called the generators of the symmetry. We will now see how this works out in the case of translations.

Thinking of a small displacement aa, we expand the translation operator T(a)T(a) in a Taylor series in powers of aa. This series begins with

T(a)=T(0)+adTda(0)+….T(a) = T(0) + a \dod{T}{a}(0) + \ldots.

The first (zeroth order) term is T(0)=1T(0)=1, and as for the second (first order) term, we notice that dT/da∣a=0dT/da|_{a=0} does not depend on aa so it is some operator with no parameters. We put this into a more convenient form by defining

k^=idTda(0),\khat = i\dod{T}{a}(0),

where the hat emphasizes that k^\khat is an operator (in contrast to aa, which is a number).

In these notes we frequently use a hat to distinguish an operator from an ordinary number, or from the classical counterpart of the operator. We normally omit the hat when there is no danger of confusion (for example, TT above is an operator, but we put no hat on it). But we also use the hat for unit vectors, and certain other purposes. We will explain these as they arise, so the meaning will be clear in any particular context.

Through first order, we can write the expansion of the translation operator as

T(a)=1−iak^+….T(a) = 1 -ia\khat+ \ldots.

We have split off a factor of ii in the definition (34) so that k^\khat will be Hermitian. This follows when we use (35) to write out the series for T(a)†T(a)^\hc and T(−a)=T(a)−1T(-a)=T(a)^{-1}:

T(a)†=1+iak^†+…,T(−a)=1+iak^+…,\begin{aligned} T(a)^\hc &= 1 + ia\khat^\dagger + \ldots, \\ T(-a) &= 1 + ia \khat + \ldots, \\ \end{aligned}

which by (29) and (30) must be equal. But this implies

k^=k^†.\khat = \khat^\dagger.

The Hermitian operator k^\khat appears in the first correction term in (35), which, when aa is small, is an infinitesimal translation. For this reason k^\khat is regarded as the generator of translations in one dimension. That is, if we knew how to act on states with the operator k^\khat, we could bring about an infinitesimal translation; and finite translations are compositions of a large number of small translations. As yet we have no physical interpretation for k^\khat, but since it is Hermitian it must correspond to the measurement of some physical quantity. We will see in a moment what that is.

Exponential Form for Translation Operators

We have only written out the first two terms in the Taylor series (35) for T(a)T(a), but the entire series can be summed and put into a neat form. To do this we will obtain a differential equation for T(a)T(a) and then solve it. See Prob. Paragraph(a), which uses similar techniques. We first write out the definition of the derivative of T(a)T(a) as a limit,

dT(a)da=lim⁡ϵ→0T(a+ϵ)−T(a)ϵ.\dod{T(a)}{a} = \lim_{\epsilon\to0} {T(a+\epsilon) - T(a) \over \epsilon}.

By (28), the first term in the numerator can be written as T(ϵ)T(a)T(\epsilon)T(a), which allows the operator T(a)T(a) to be factored out of the entire numerator to the right:

dT(a)da=(lim⁡ϵ→0T(ϵ)−1ϵ)T(a).\dod{T(a)}{a} = \Bigl( \lim_{\epsilon\to0} {T(\epsilon)-1 \over \epsilon} \Bigr) T(a).

But the remaining limit is just the derivative dT(a)/dadT(a)/da evaluated at a=0a=0, which by (34) is −ik^-i\khat. Altogether, we obtain

dT(a)da=−ik^ T(a).\dod{T(a)}{a} = -i\khat \,T(a).

This is a differential equation that we must solve subject to the initial conditions T(0)=1T(0)=1. The solution is immediate,

T(a)=exp⁡(−iak^).T(a)= \exp(-ia\khat).

Now we can easily extend the power series (35) to higher order,

T(a)=1−iak^−12a2k^2+….T(a)=1 -ia\khat -{1\over2} a^2\khat^2+ \ldots.

We see that the generator of translations k^\khat, which first appeared in infinitesimal translation operators, also appears in the exponential form for finite translation operators.

There is another way to obtain (41) that is closer to the idea of building up finite transformations from infinitesimal ones. Suppose aa is a displacement that is not small. We break it up into a large number NN of small displacements, writing ϵ=a/N\epsilon=a/N. Then property (28) implies

T(a)=T(ϵ)N.T(a) = T(\epsilon)^N.

But if ϵ\epsilon is small, we can approximate T(ϵ)T(\epsilon) by the first two terms of the Taylor series, as in (35), without assuming at this point that we know what the higher order terms are. Thus we have

T(a)≈(1−iϵk^)N.T(a) \approx (1-i\epsilon\khat)^N.

This approximation should get better as ϵ\epsilon gets smaller, that is, as NN increases. This makes it plausible that we should have the limit,

T(a)=lim⁡N→∞(1−iak^N)N.T(a) = \lim_{N\to\infty} \Bigl(1-{ia\khat\over N}\Bigr)^N.

This is similar to the limit of elementary calculus,

lim⁡N→∞(1+xN)N=ex.\lim_{N\to\infty}\Bigl(1+{x\over N}\Bigr)^N = e^x.

Applying this to (45) we obtain the exponential expression (41).

Action of k^\khat on Kets and Wave Functions

It is easy to work out the action of k^\khat on the basis kets ∣x⟩\ket{x}. We simply write

k^∣x⟩=idT(a)da∣x⟩∣a=0=idda∣x+a⟩∣a=0=iddx∣x⟩.\khat \ket{x} = i\Bigl.\dod{T(a)}{a} \ket{x}\Bigr|_{a=0} = i\Bigl.\dod{}{a}\ket{x+a}\Bigr|_{a=0} = i\dod{}{x} \ket{x}.

The derivative of the ket ∣x⟩\ket{x} may look strange. In wave function language it is the derivative of a δ\delta-function. One can use expressions with derivatives of kets as in (47), formally integrating by parts as if they were ordinary functions, and obtain correct answers.

We obtain an equivalent result that looks more familiar by working out the action of k^\khat on a wave function. The procedure is the same:

(k^ψ)(x)=i(dT(a)daψ)(x)∣a=0=iddaψ(x−a)∣a=0=−idψ(x)dx.(\khat\psi)(x) = i\Bigl. \Bigl(\dod{T(a)}{a}\psi\Bigr)(x) \Bigr|_{a=0} = i \Bigl. \dod{}{a} \psi(x-a)\Bigr|_{a=0} = -i\dod{\psi(x)}{x}.

We start to see the appearance of the usual momentum operator on wave functions.

Translations in Three Dimensions, and Commutation Relations

In a similar manner, for translation operators in three dimensions we define

k^i=i∂T(a)∂ai∣a=0,\khat_i = i \Bigl.\frac{\partial T(\avec)}{\partial a_i}\Bigr|_{\avec=0},

where i=1,2,3i=1,2,3. This gives us a Hermitian vector operator (that is, a vector of Hermitian operators) k^{\hat{\mathbf{k}}}, with components k^i\khat_i, i=1,2,3i=1,2,3. Now the exponential form of the translation operator is

T(a)=exp⁡(−ia⋅k^)=1−ia⋅k^−12(a⋅k^)2+….T(\avec) = \exp(-i\avec\cdot{\hat{\mathbf{k}}}) = 1 -i\avec\cdot{\hat{\mathbf{k}}} -{1\over 2}(\avec\cdot{\hat{\mathbf{k}}})^2 +\ldots.

This series can be used to obtain the commutation relations of the operators k^=(k^1,k^2,k^3){\hat{\mathbf{k}}}=(\khat_1,\khat_2,\khat_3). We expand the product T(a)T(b)T(\avec)T(\bvec) in power series, carrying everything to second order,

T(a)T(b)=[1−i(a⋅k^)−12(a⋅k^)2+…][1−i(b⋅k^)−12(b⋅k^)2+…]\begin{aligned} T(\avec)T(\bvec) & = \Bigl[1-i(\avec\cdot{\hat{\mathbf{k}}}) -{1\over2}(\avec\cdot{\hat{\mathbf{k}}})^2 + \ldots \Bigr] \Bigl[1-i(\bvec\cdot{\hat{\mathbf{k}}}) -{1\over2}(\bvec\cdot{\hat{\mathbf{k}}})^2 + \ldots \Bigr] \end{aligned}
=1−i(a+b)⋅k^−12(a⋅k^)2−(a⋅k^)(b⋅k^)−12(b⋅k^)2+….\begin{aligned} &= 1-i(\avec+\bvec)\cdot{\hat{\mathbf{k}}} -{1\over2}(\avec\cdot{\hat{\mathbf{k}}})^2 -(\avec\cdot{\hat{\mathbf{k}}})(\bvec\cdot{\hat{\mathbf{k}}}) -{1\over2}(\bvec\cdot{\hat{\mathbf{k}}})^2 + \ldots. \end{aligned}

But since T(a)T(b)=T(b)T(a)T(\avec)T(\bvec) =T(\bvec) T(\avec), the answer must be the same if we swap a\avec and b\bvec. Most of the series is obviously symmetric in a\avec and b\bvec, but the one term (a⋅k^)(b⋅k^)(\avec\cdot{\hat{\mathbf{k}}}) (\bvec\cdot{\hat{\mathbf{k}}}) is not. When we subtract the swapped series from the original series, we obtain a vanishing commutator,

[a⋅k^,b⋅k^]=0.[\avec\cdot{\hat{\mathbf{k}}}, \bvec\cdot{\hat{\mathbf{k}}}]=0.

Since this is true for all choices of a\avec and b\bvec, we can choose each of these vectors to be one of the unit vectors along the coordinate axes (nine choices in all), and we find

[k^i,k^j]=0.[\khat_i, \khat_j]=0.

We see that k^{\hat{\mathbf{k}}} is a vector of commuting, Hermitian operators.

The relations (47) and (48) are also easily generalized to three dimensions. These give

k^i∣x⟩=i∂∂xi∣x⟩,\khat_i \ket{\xvec} = i\pop{}{x_i} \ket{\xvec},

and

(k^iψ)(x)=−i∂ψ(x)∂xi,\bigl(\khat_i\psi\bigr)(\xvec) = -i \frac{\partial \psi(\xvec)}{\partial x_i},

or, in vector notation,

(k^ψ)(x)=−i∇ψ(x).\bigl({\hat{\mathbf{k}}}\psi\bigr)(\xvec) = -i\del\psi(\xvec).

(56) provides another way to derive the commutation relations (54), that is, the derivatives ∂/∂xi\partial/\partial x_i and ∂/∂xj\partial/\partial x_j commute with one another.

It is also easy to work out the commutation relations among the operators x^i\xhat_i and k^j\khat_j. We just apply them to a wave function in opposite orders, obtaining

x^ik^jψ=xi(−i∂∂xj)ψ=−ixi∂ψ∂xj,k^jx^iψ=(−i∂∂xj)xiψ=−i δij ψ−ixi∂ψ∂xj.\begin{aligned} \xhat_i\khat_j\psi &= x_i \Bigl(-i\pop{}{x_j}\Bigr)\psi = -ix_i\frac{\partial \psi}{\partial x_j}, \\ \khat_j\xhat_i\psi &= \Bigl(-i\pop{}{x_j}\Bigr)x_i\psi = -i\,\delta_{ij}\,\psi -i x_i\frac{\partial \psi}{\partial x_j}. \end{aligned}

Subtracting these, we find

[x^i,k^j]=i δij.[\xhat_i, \khat_j] = i\, \delta_{ij}.

Combined with (11) and (54), this gives a complete set of commutation relations among the components of the operators x^{\hat{\mathbf{x}}} and k^{\hat{\mathbf{k}}}.

The Momentum in Classical Mechanics

We shall make a few remarks about the momentum in classical mechanics, in preparation for our discussion of the momentum in quantum mechanics.

As explained in The Momentum, there is more than one definition of the momentum of particle in classical mechanics. One is the kinetic momentum,

p=mx˙,\pvec = m \dot\xvec,

and the other is the canonical momentum,

p=∂L∂x˙,\pvec= \frac{\partial L}{\partial \dot\xvec},

where LL is the Lagrangian. In the case of a particle of charge qq moving in static electric and magnetic fields E=−∇Φ\Evec=-\del \Phi and B=∇×A\Bvec = \del \cross \Avec, the Lagrangian is

L(x,x˙)=m2∣x˙∣2−qΦ(x)+qcx˙⋅A(x),L(\xvec,\dot\xvec) = {m \over 2} | \dot \xvec|^2 -q\Phi(\xvec) + {q \over c} \dot\xvec \cdot \Avec(\xvec),

so the canonical momentum is given by (35), reproduced here,

p=mx˙+qcA(x).\pvec= m \dot\xvec + {q\over c} \Avec(\xvec).

In the presence of a magnetic field (more precisely, in the presence of a vector potential) the definitions of the classical momentum (60) and (61) are not the same. We have not defined the momentum operator in quantum mechanics yet, but before we do, we must ask which of these two classical momenta will it correspond to, the kinetic or the canonical?

We have been developing the properties of the operator k^\khat in quantum mechanics, which is the generator of displacements. It turns out that there is a similar role played by the momentum in classical mechanics, but it is the canonical momentum, not the kinetic, that does this. Let us write f(x,p)f(x,p) be a function of the position and momentum for a particle moving in one dimension. We can think of ff as a classical observable. Then we can define a classical translation operator that acts on classical observables according to

(Tcl(a)f)(x,p)=f(x−a,p).\bigl( T_cl(a) f\bigr)(x,p) =f(x-a,p).

This operator moves the function f(x,p)f(x,p) forward in xx while leaving pp alone. If now aa is a small displacement, then the result can be expanded to first order,

(Tcl(a)f)(x,p)=f(x,p)−a∂f∂x.\bigl( T_cl(a) f\bigr)(x,p) = f(x,p) - a\frac{\partial f}{\partial x}.

But the correction term can be written in terms of a Poisson bracket of ff with the momentum,

−a∂f∂x=a{p,f}.-a \frac{\partial f}{\partial x}= a\{ p, f \}.

Finite displacements can be built up by composing infinitesimal displacements, and the result can be expressed in terms of an exponential series of iterated Poisson brackets. We introduce the following notation for the Poisson bracket,

LXA={X,A},L_X A = \{X,A\},

where XX and AA are arbitrary classical observables. Then a finite displacement can be expressed as

Tcl(a)=exp⁡(aLp),T_cl(a) = \exp(a L_p),

that is,

f(x−a,p)=(Tcl(a)f)(x,p)=(exp⁡(aLp)f)(x,p)=f(x,p)+a{p,f}+a22!{p,{p,f}}+…=f(x,p)−a∂f∂x(x,p)+a22!∂2f∂x2(x,p)+…\begin{aligned} f(x-a,p) &= \bigl(T_cl(a)f\bigr)(x,p) = \bigl(\exp(aL_p) f\bigr)(x,p) \\ &=f(x,p) + a \{ p,f\} + {a^2\over 2!} \{p,\{p,f\}\} + \ldots \\ &= f(x,p) - a\frac{\partial f}{\partial x}(x,p) + {a^2\over 2!} {\partial^2 f \over \partial x^2}(x,p) +\ldots \\ \end{aligned}

which may be compared to (41) and (42). The result is the Taylor series expansion of f(x,p)f(x,p) in the variable xx. It may also be compared to the exponential series developed in Prob. Paragraph(a).

These results are easily generalized to three dimensions. In three dimensions, the distinction between the kinetic and canonical momentum becomes important. It is the canonical momentum, not the kinetic, that generates translations in classical mechanics (in cases where it makes a difference). That is, the Poisson brackets are computed with respect to the canonical momentum, not the kinetic. See (100).

The role that momentum plays in classical mechanics as the generator of displacements is closely connected with the conservation of momentum in systems with translational symmetry. This conservation law is a special case of Noether’s theorem, which presents a general relationship between continuous symmetries and conserved quantities.

The Momentum in Quantum Mechanics

Let us work in three dimensions. According to the postulates of quantum mechanics, measurements of the momentum p\pvec of a particle must correspond to some operator (actually, a vector of three operators), call it p^{\hat{\mathbf{p}}}, where the hat distinguishes the classical or cc-number momentum p\pvec from the operator. At this point we know what the Hilbert space is for particles moving in three dimensions; it is a ket space E\ES that isomorphic to the space of wave functions ψ(x)\psi(\xvec). But what does the momentum operator p^{\hat{\mathbf{p}}} do to states in this space? That is, what is the momentum operator explicitly?

The operator k^{\hat{\mathbf{k}}} is the generator of displacements in quantum mechanics, since by (35) it provides the small correction necessary when an infinitesimal displacement is carried out. Also, as noted, k^{\hat{\mathbf{k}}} corresponds to the measurement of some physical quantity, and both k^{\hat{\mathbf{k}}} and p\pvec (the latter being the generator of translations in classical mechanics) are vectors. So the suspicion arises that k^{\hat{\mathbf{k}}} is related to the momentum operator p^{\hat{\mathbf{p}}} in quantum mechanics. The vector operators k^{\hat{\mathbf{k}}} and p^{\hat{\mathbf{p}}} cannot be equal, however, because they have different units: k^{\hat{\mathbf{k}}} has units of inverse length, and p^{\hat{\mathbf{p}}} has units of momentum. To make the units come out right we guess that there is a proportionality factor between p^{\hat{\mathbf{p}}} and k^{\hat{\mathbf{k}}} with units of action, a constant that we call ℏ\hbar:

p^=ℏk^.{\hat{\mathbf{p}}} = \hbar {\hat{\mathbf{k}}}.

We take this as the definition of momentum in quantum mechanics.

The postulates of quantum mechanics do not tell us what quantum operator corresponds to a given classical quantity, so we must not think that the definition of momentum p^{\hat{\mathbf{p}}} is engraved in stone. The best we can do is to find a vector of operators in quantum mechanics that has properties similar to the classical momentum and that goes over to the classical momentum in the classical limit. In fact, the similarity with the classical momentum (as the generator of translations) has motivated the definition (70), and, as for the classical limit and the other expected properties of momentum, they will be laid bare as we proceed.

The Constant ℏ\hbar

Notice that ℏ\hbar did not appear in the postulates of quantum mechanics; it appears for the first time in our development of the formalism in (70). The value of ℏ\hbar in ordinary units must be determined experimentally. The discussion so far has been based on the analysis of a single particle moving in three-dimensional space, so one might question whether different particles have different values of ℏ\hbar. We would not get the usual classical limit if this were so, so it seems theoretically doubtful that there are different values of ℏ\hbar, but it is a question that can be tested experimentally (all the ℏ\hbar’s turn out to be the same to within experimental accuracy). Given that there is only one ℏ\hbar, there is the option of choosing “natural units” in which ℏ=1\hbar=1.

The relation (70) implies the usual de~Broglie relation connecting momentum and wavelength,

λ=2πℏp,\lambda={2\pi\hbar \over p},

as we shall show momentarily, so the value of ℏ\hbar can be determined by measuring the de~Broglie wavelength of a particle of a known momentum. The most elegant experiments along these lines have involved Bragg diffraction of a neutron beam from a single, large (about 10cm) silicon crystal with no dislocations or grain boundaries. The manufacture of such crystals has become feasible with modern semiconductor technology. The predictions of the relation (70) are fully confirmed. Actually, if (70) were not correct, there are very few of the theoretical predictions of quantum mechanics that would agree with experiment.

The silicon neutron interferometer just mentioned was used in the 1970’s to measure the phase shift of neutron wave functions as the neutrons fall in the earth’s gravitational field. This was the first demonstration that gravitational potentials enter the Schr"odinger equation in the same manner as other potentials. The experiment and results are discussed in Sakurai, Modern Quantum Mechanics and in Commins, Quantum Mechanics: An Experimentalist’s Approach. Commins has a particulary nice explanation of the neutron interferometer.

The fact that ℏ\hbar is a universal constant was known to Planck and others in the 1890’s, even before quantum mechanics existed. Planck knew that a constant with dimensions of action must occur in the expression for the spectrum of black body radiation, even before the correct expression was known. He also realized that h=2πℏh=2\pi\hbar, cc and GG (Newton’s constant of gravitation) could be combined to create a natural system of units for distance, time and mass. The Planck unit of distance is

LPlanck=ℏGc3=1.6×10−33cm.L_Planck = \sqrt{{\hbar G \over c^3}} = 1.6 \times 10^{-33}cm.

It is the length scale at which the effects of both quantum mechanics and gravity are expected to be important. Physics at such scales is currently a matter of much speculation, and is likely to remain so for some time, in view of the near impossibility of direct experimental tests. Nevertheless the importance of the length scale itself has been known for a long time.

The Momentum in Quantum Mechanics: The History

The pseudo-axiomatic approach we are taking here to the development of the momentum operator in quantum mechanics is quite different from the historical one. In a series of papers between 1905 and 1916, Einstein developed the idea that light of frequency ω\omega and wave vector k\kvec is associated with particles of energy E=ℏωE=\hbar\omega and momentum p=ℏk\pvec=\hbar\kvec. We now call these particles “photons,” although that name was not coined until later. According to Maxwell’s equations the frequency and wave number of a light wave are related by ω=ck\omega=ck, so Einstein’s formulas give the energy-momentum relation for the photon as E=cpE=cp, exactly what the special theory of relativity predicts for a massless particle traveling at the speed of light. Einstein’s arguments relied heavily on statistical mechanics. Most physicists of the time regarded the idea of particles of light as completely crazy, since it had been established a hundred years earlier that light was a wave. But in 1922 Compton showed experimentally that collisions between x-rays and electrons obey the relativistic rules of conservation of energy and momentum, assuming Einstein’s expressions for the energy and momentum of the photon. With this physicists began to take the wave-particle duality more seriously. In 1924, de~Broglie argued that if light could have both wave and particle aspects, then so might ordinary particles also have wave aspects. This was in part an appeal to the unity of physics, and a suggestion offered by the equivalence of energy and mass, as shown by special relativity. De~Broglie suggested that the formulas E=ℏωE=\hbar\omega and p=ℏk\pvec=\hbar\kvec should apply also to massive particles, giving the frequency and wave vector of a wave associated with the particles. Somewhat later Debye remarked that if there was a wave, there should be a wave equation, and Schr"odinger, on hearing of this, set out to find it. He first tried a relativistic wave equation without success, only later trying the nonrelativistic version which gave an explanation of the hydrogen spectra and other physical facts. The latter equation is the one usually associated with his name.

A good recent book that covers this history and much else is Einstein and the Quantum by A. Douglas Stone.

Properties of the Momentum

We will provisionally accept the definition (70) of the momentum operator in quantum mechanics, and proceed to examine the consequences. First, the translation operators can now be written in the form,

T(a)=exp⁡[−iℏ(a⋅p^)]=1−iℏ(a⋅p^)+….T(\avec) = \exp\Bigl[-{i\over\hbar}(\avec\cdot{\hat{\mathbf{p}}}) \Bigr] = 1 - {i\over\hbar}(\avec\cdot{\hat{\mathbf{p}}}) + \ldots.

The various commutation relations we have evaluated may now be expressed in terms of p^{\hat{\mathbf{p}}} instead of k^{\hat{\mathbf{k}}} and gathered together in one place,

[x^i,x^j]=0,\begin{aligned} [\xhat_i, \xhat_j] &=0, \end{aligned}
[x^i,p^j]=iℏ δij,\begin{aligned} [\xhat_i, \phat_j] &= i\hbar \, \delta_{ij}, \end{aligned}
[p^i,p^j]=0.\begin{aligned} [\phat_i,\phat_j] &=0. \end{aligned}

These are the Heisenberg-Born commutation relations. They may be compared to the classical canonical Poisson bracket relations, (108).

Also, the properties (55) and (57) of the operators k^{\hat{\mathbf{k}}} are easily translated into the properties of the momentum. We have

p^∣x⟩=iℏ∇∣x⟩,{\hat{\mathbf{p}}} \ket{\xvec} = i\hbar \del \ket{\xvec},

and

(p^ψ)(x)=−iℏ∇ψ(x).({\hat{\mathbf{p}}}\psi)(\xvec) = -i\hbar \del\psi(\xvec).

(77) gives the action of the momentum operator in ket language, whereas (78) gives it in the position representation. It has other forms in other representations. Another form of these equations is

⟨x∣p^∣ψ⟩=−iℏ∇⟨x∣ψ⟩.\matrixelement{\xvec}{{\hat{\mathbf{p}}}}{\psi} =-i\hbar\del \braket{\xvec}{\psi}.

This is really a rephrasing of (78), and can be obtained by taking the Hermitian conjugate of (77) and then multiplying by the ket ∣ψ⟩\ket{\psi}.

The Momentum Representation

Because p^{\hat{\mathbf{p}}} is a vector of commuting, Hermitian operators, these operators possess a simultaneous eigenbasis. We write the eigenvalue-eigenket problem for momentum in the form,

p^∣p⟩=p∣p⟩.{\hat{\mathbf{p}}} \ket{\pvec} = \pvec \ket{\pvec}.

We solve this in the configuration representation, that is, we multiply (80) on the left by the bra ⟨x∣\bra{\xvec}, define the wave function of the momentum eigenket ∣p⟩\ket{\pvec} by

ψp(x)=⟨x∣p⟩,\psi_\pvec(\xvec) = \braket{\xvec}{\pvec},

and use (78) to obtain

−iℏ∇ψp(x)=pψp(x).-i\hbar \del \psi_\pvec(\xvec) = \pvec \psi_\pvec(\xvec).

Thus, ψp(x)\psi_\pvec(\xvec) is the eigenfunction of the differential operator −iℏ∇-i\hbar\del (really three commuting differential operators) with eigenvalue p\pvec. The solution exists for all values of p\pvec and is unique to within a constant,

ψp(x)=Aeip⋅x/ℏ,\psi_\pvec(\xvec) = A e^{i\pvec\cdot\xvec/\hbar},

where AA is a normalization and a phase. Thus, momentum has a continuous spectrum, so we normalize the eigenstates according to

⟨p1∣p2⟩=δ3(p1−p2).\braket{\pvec_1}{\pvec_2} = \delta^3(\pvec_1-\pvec_2).

Using the integral

∫d3x ei(p2−p1)⋅x/ℏ=(2πℏ)3δ(p1−p2),\int d^3\xvec \, e^{i(\pvec_2-\pvec_1)\cdot \xvec/\hbar} = (2\pi\hbar)^3 \delta(\pvec_1 - \pvec_2),

the normalization constant AA in (83) is determined to within a phase, which we fix by demanding that AA be real and positive. The result is

⟨x∣p⟩=1(2πℏ)3/2 eip⋅x/ℏ.\boxed{\braket{\xvec}{\pvec} = {1 \over (2\pi\hbar)^{3/2}} \, e^{i\pvec\cdot\xvec/\hbar}.}

The eigenfunctions of momentum are plane waves, whose wave length is given by the de~Broglie relation (71).

Since the momentum eigenstates are nondegenerate, the momentum operators p^{\hat{\mathbf{p}}} form a complete set of commuting observables, and we may speak of the momentum representation. The resolution of the identity in this representation is

1=∫d3p ∣p⟩⟨p∣.1 = \int d^3\pvec \, \ketbra{\pvec}{\pvec}.

We will denote the wave function of the state ∣ψ⟩\ket{\psi} in the momentum representation by

ϕ(p)=⟨p∣ψ⟩\phi(\pvec) = \braket{\pvec}{\psi}

(compare (15) for the configuration representation), so that by multiplying (87) onto the ket ∣ψ⟩\ket{\psi} we obtain

∣ψ⟩=∫d3p ∣p⟩ϕ(p),\ket{\psi} = \int d^3\pvec \, \ket{\pvec} \phi(\pvec),

similar to (16) in the configuration representation. The wave function ϕ(p)\phi(\pvec) is the expansion coefficients of the state ∣ψ⟩\ket{\psi} with respect to the momentum eigenbasis {∣p⟩}\{ \ket{\pvec} \}.

By using (86) and resolutions of the identity it is easy to switch back and forth between wave functions ψ(x)\psi(\xvec) and ϕ(p)\phi(\pvec), representing the same quantum state in two different representations. For example, we have

ψ(x)=⟨x∣ψ⟩=∫d3p ⟨x∣p⟩⟨p∣ψ⟩=∫d3p(2πℏ)3/2 eip⋅x/ℏ ϕ(p),\psi(\xvec) = \braket{\xvec}{\psi} = \int d^3\pvec \, \braket{\xvec}{\pvec} \braket{\pvec}{\psi} = \int {d^3\pvec \over (2\pi\hbar)^{3/2}} \, e^{i\pvec\cdot\xvec/\hbar} \, \phi(\pvec),

and its inverse,

ϕ(p)=∫d3x(2πℏ)3/2 e−ip⋅x/ℏ ψ(x).\phi(\pvec) = \int {d^3\xvec \over (2\pi\hbar)^{3/2}} \, e^{-i\pvec\cdot\xvec/\hbar} \, \psi(\xvec).

The wave functions ψ(x)\psi(\xvec) and ϕ(p)\phi(\pvec) are Fourier transforms of one another, modulo the insertion of factors of ℏ\hbar to account for the physical units.

In the momentum representation, the operator p^{\hat{\mathbf{p}}} is represented simply by multiplication by p\pvec, the vector of cc-numbers. That is, if ϕ(p)\phi(\pvec) is given in terms of the state ∣ψ⟩\ket{\psi} by (88), then the momentum space wave function of the state p^∣ψ⟩{\hat{\mathbf{p}}}\ket{\psi} is

⟨p∣p^∣ψ⟩=p⟨p∣ψ⟩=p ϕ(p),\matrixelement{\pvec}{{\hat{\mathbf{p}}}}{\psi} = \pvec \braket{\pvec}{\psi} = \pvec \,\phi(\pvec),

where we have allowed p^{\hat{\mathbf{p}}} to act to the left on the bra ⟨p∣\bra{\pvec}, bringing out the eigenvalue p\pvec. We can write this as

(p^ϕ)(p)=pϕ(p).({\hat{\mathbf{p}}}\phi)(\pvec) = \pvec\phi(\pvec).

Similarly, the operator x^{\hat{\mathbf{x}}} in the momentum representation is given by

(x^ϕ)(p)=iℏ∂ϕ(p)∂p.({\hat{\mathbf{x}}} \phi)(\pvec) = i\hbar \frac{\partial \phi(\pvec)}{\partial \pvec}.

Compare this to (78), and notice the difference in sign. This follows since if ϕ(p)\phi(\pvec) is related to ∣ψ⟩\ket{\psi} by (88), then what we mean by the left-hand side of (94) is

⟨p∣x^∣ψ⟩=∫d3x ⟨p∣x⟩⟨x∣x^∣ψ⟩=∫d3x(2πℏ)3/2 e−ip⋅x/ℏ xψ(x)\begin{aligned} \matrixelement{\pvec}{{\hat{\mathbf{x}}}}{\psi} &= \int d^3\xvec \, \braket{\pvec}{\xvec} \matrixelement{\xvec}{{\hat{\mathbf{x}}}}{\psi} = \int {d^3\xvec \over (2\pi\hbar)^{3/2}} \, e^{-i\pvec\cdot\xvec/\hbar} \, \xvec \psi(\xvec) \end{aligned}
=iℏ∂∂p∫d3x(2πℏ)3/2 e−ip⋅x/ℏ ψ(x)=iℏ∂ϕ(p)∂p.\begin{aligned} &= i\hbar \pop{}{\pvec} \int {d^3\xvec \over (2\pi\hbar)^{3/2}} \, e^{-i\pvec\cdot\xvec/\hbar} \, \psi(\xvec) = i\hbar \frac{\partial \phi(\pvec)}{\partial \pvec}. \end{aligned}

Finally, we summarize the representations of the operators x^i\xhat_i and p^i\phat_i in the configuration and momentum representations in Table 1.

Table 1:Representations of the operators x^i\xhat_i and p^i\phat_i when acting on configuration and momentum space wave functions. \par

x^i\xhat_i

p^i\phat_i

Configuration

mult by xi\text{mult by } x_i

-iℏ∂∂xii\hbar\pop{}{x_i}

Momentum

iℏ∂∂pii\hbar\pop{}{p_i}

mult by pi\text{mult by } p_i

Multiparticle Wave Functions

In a system of NN particles, the positions of the individual particles are independent observables that commute with one another, so a complete set (ignoring spin for now) consists of the operators (x^1,…,x^N)({\hat{\mathbf{x}}}_1, \ldots, {\hat{\mathbf{x}}}_N), with eigenvalues (x1,…,xN)(\xvec_1, \ldots, \xvec_N). In this case the wave function is defined by

ψ(x1,…,xN)=⟨x1,…,xN∣ψ⟩.\psi(\xvec_1, \ldots, \xvec_N) = \braket{\xvec_1, \ldots, \xvec_N}{\psi}.

The wave function is defined over configuration space, that is, the space in which a single point specifies the positions of all the particles. This is the same configuration space as in classical mechanics, which is discussed in Constrained Systems, Configuration Space, and Generalized Coordinates. Configuration space coincides with physical space only in the case of a single particle. Similarly, one can define a multiparticle, momentum space wave function ϕ(p1,…,pN)\phi(\pvec_1, \ldots,\pvec_N).

If one is careful about the application of the postulates of quantum mechanics, one will see certain subtleties in the derivation of (97) in the case of identical particles. We will examine this question more carefully in 29. Identical Particles.

The Sign of ii

The following is a remark concerning the the definition of k^{\hat{\mathbf{k}}} in (34), which led to the definition of momentum in (70). We split off a factor of ii in (34) to make k^{\hat{\mathbf{k}}} Hermitian, but the same would have been achieved if we had split off −i-i (thereby changing the definition of k^{\hat{\mathbf{k}}} by a sign). This would lead to the opposite sign in the definition of p^{\hat{\mathbf{p}}}, and changes in signs in many of the subsequent formulas. Would this change lead to any physical consequences or contradictions with experiment?

The answer is no, but it would change most of the familiar formulas in quantum mechanics, by replacing ii by −i-i. For example, a plane wave with wave vector k\kvec would become e−ik⋅xe^{-i\kvec\cdot\xvec} instead of the usual eik⋅xe^{i\kvec\cdot\xvec}. It is a matter of convention to choose the sign of ii in quantum mechanics, and our choice has been made in (34) and (70). Once this choice has been made, however, then the sign of the Pauli matrix σy\sigma_y is determined, so that spin angular momentum has the same commutation relations as orbital angular momentum. This was a question addressed in Prob. Paragraph(d).

The Position-Momentum Uncertainty Relation

The generalized uncertainty relations (53), applied to the case of position xx and momentum pp of a particle in one-dimension, for which [x,p]=iℏ[x,p]=i\hbar, give

Δx Δp≥ℏ2.\Delta x\,\Delta p \ge {\hbar\over 2}.

This can be used in a number of applications to obtain quick estimates of orders of magnitude of physical quantities, including three-dimensional problems where for simplicity we ignore the vector nature of x\xvec and p\pvec. We give several examples.

Consider a beam of particles of well defined momentum p=p0 z^\pvec=p_0\, {\hat{\mathbf{z}}} directed at a screen in the plane z=0z=0. The screen has a hole of dimension LL (the hole need not be circular, LL is just an order of magnitude estimate of the size of the hole). The hole confines the particles in the xx- and yy-directions to a distance of order LL, so by (98) the particles acquire a transverse momentum (in the xx- or yy-directions) roughly of magnitude

p⊥=ℏL,p_\perp = {\hbar\over L},

where we drop constants of order unity since this is only an estimate. Thus downstream from the hole the beam should spread at an angle given roughly by

θ=p⊥p0=ℏLp0=λ2πL,\theta = {p_\perp \over p_0} = {\hbar \over Lp_0} = {\lambda\over 2\pi L},

where λ\lambda is the wave length of the incident waves. This calculation assumes the angle θ\theta is small, which is equivalent to λ≪L\lambda \ll L, that is, the hole should be at least a few wavelengths in diameter. This very simple estimate of the spreading angle of the beam is confirmed by more complicated calculations of diffraction theory in the far-field (Fraunhofer) region.

If we eliminate all momenta pp in favor of wave numbers kk, by p=ℏkp=\hbar k, then all the ℏ\hbar’s disappear and we have results that can be applied also to classical waves (light waves, sound waves, etc).

For another example, let aa be an estimate of the size of a hydrogen atom. The electron is confined to a distance of order aa and so must have a momentum of order p=ℏ/ap=\hbar/a. In a circular orbit this implies a centrifugal force of

mv2a=p2ma=ℏ2ma3.{mv^2\over a} = {p^2\over ma} = {\hbar^2\over ma^3}.

Equating this to the electrostatic force on the electron from the proton, e2/a2e^2/a^2, and solving for aa, we obtain

a=ℏ2me2.a = {\hbar^2\over me^2}.

This is the usual expression for the Bohr radius. Given an estimate for the size of an atom and the known density of ordinary matter, one also obtains an estimate for the Avogadro number and many other physical quantities.

Minimum Uncertainty Wave Packets

Let us loosely define a wave packet in one dimension as a wave function ψ(x)\psi(x) whose dispersions Δx\Delta x and Δp\Delta p, defined as in Statistical Aspects of Measurement by

(Δx)2=⟨x^2⟩−⟨x^⟩2,(Δp)2=⟨p^2⟩−⟨p^⟩2,\begin{aligned} (\Delta x)^2 &= \xpecval{\xhat^2} - \xpecval{\xhat}^2, \\ (\Delta p)^2 &= \xpecval{\phat^2} - \xpecval{\phat}^2, \\ \end{aligned}

are small, that is, close to the minimum values allowed by the inequality,

ΔxΔp≥ℏ2.\Delta x \Delta p \ge {\hbar\over 2}.

See (54). Similarly, let us define a minimum uncertainty wave packet as one for which the product ΔxΔp\Delta x \Delta p takes on its minimum value of ℏ/2\hbar/2.

Let ψ(x)\psi(x) be a normalized, minimum uncertainty wave packet with

⟨ψ∣x^∣ψ⟩=⟨x^⟩=0,⟨ψ∣p^∣ψ⟩=⟨p^⟩=0,\begin{aligned} \matrixelement{\psi}{\xhat}{\psi} &= \xpecval{\xhat} =0, \\ \matrixelement{\psi}{\phat}{\psi} &= \xpecval{\phat} =0, \\ \end{aligned}

and let Δx=L\Delta x = L, so that Δp=ℏ/2L\Delta p = \hbar/2L. Now let

∣ζ⟩=(x^+2iL2ℏp^)∣ψ⟩.\ket{\zeta} = \Bigl(\xhat + {2iL^2\over\hbar}\phat\Bigr) \ket{\psi}.

The operator appearing on the right seems to have been pulled out of the air, but it is a version of an annihilation operator of the general form x^+ip^\xhat+i\phat, with coefficients adjusted to make the following argument come out right. We will see such annihilation operators again when we discuss the harmonic oscillator. Now squaring (106), we obtain

⟨ζ∣ζ⟩=⟨x^2⟩+2iL2ℏ⟨x^p^−p^x^⟩+4L4ℏ2⟨p^2⟩=L2−2L2+L2=0.\braket{\zeta}{\zeta} = \xpecval{\xhat^2} +{2iL^2 \over \hbar} \xpecval{\xhat\phat-\phat\xhat} + {4L^4\over\hbar^2} \xpecval{\phat^2} = L^2 - 2L^2 + L^2 = 0.

This implies ∣ζ⟩=0\ket{\zeta}=0 (see (28), so (106) in wave function language becomes

2L2dψdx+xψ=0,2L^2 \dod{\psi}{x} + x\psi=0,

which has the normalized solution,

ψ(x)=1L2π e−x2/4L2.\psi(x) = {1\over{\sqrt{L\sqrt{2\pi}}}}\, e^{-x^2/4L^2}.

The minimum uncertainty wave packet is a Gaussian. If we allow ⟨x^⟩=a\xpecval{\xhat}=a and ⟨p^⟩=b\xpecval{\phat}=b to take on arbitrary values aa, bb as shown, then we find

ψ(x)=1L2πexp⁡[−(x−a)24L2+ib(x−a)ℏ+iγ],\psi(x) = {1\over{\sqrt{L\sqrt{2\pi}}}} \exp\Bigl[ -{(x-a)^2 \over 4L^2} + i{b(x-a) \over \hbar} + i\gamma\Bigr],

where γ\gamma is a phase. The wave packet is still a Gaussian, but it has been shifted in position and momentum.

Problems

Problem 1.

Let ∣ψ⟩\ket{\psi} be the state of a spinless particle in three dimensions, and let ϕ(p)=⟨p∣ψ⟩\phi(\pvec) = \braket{\pvec}{\psi} be its momentum space wave function. Find the momentum space wave function of the state T(a)∣ψ⟩T(\avec)\ket{\psi}, that is, find the action of the translation operator T(a)T(\avec) in the momentum representation.

Problem 2.

Consider a particle moving in one dimension. In classical mechanics, the state of the particle is given by its position xx and momentum pp, and other observables such as energy are functions of these. The classical position and momentum are ordinary numbers, and obey the usual (commutative) rules of algebra, for example, xp=pxxp=px. In quantum mechanics, we have the operators x^\xhat and p^\phat which are analogous to the classical xx and pp, and other observables such as the Hamiltonian are functions of x^\xhat and p^\phat. On the other hand, x^\xhat and p^\phat do not obey a commutative algebra, in view of the relation x^p^=p^x^+iℏ\xhat\phat = \phat\xhat + i\hbar.

A question that has exercised various people is the relation between classical observables (that is, functions of }xx and pp) and quantum observables (that is, functions of x^\xhat and p^\phat). One problem is to quantize a classical observable, that is, given a classical function A(x,p)A(x,p), what is the corresponding quantum observable? Dirac suggested that xx and pp should be associated with x^\xhat and p^\phat, but this gives rise to an ambiguity in cases like the classical function xpxp. Should it correspond to x^p^\xhat\phat or to p^x^\phat\xhat, or maybe to their average (which at least would be Hermitian)?

Another answer to this question was given by Weyl. We introduce the Weyl transform in this problem because it is good exercise in the position and momentum representations.

We denote operators with a hat, as in x^\hat x or A^\hat A, and we denote eigenvalues or classical quantities without a hat, as in xx or A(x,p)A(x,p). We work in one dimension, and think of a wave function ψ(x)\psi(x) or ψ(x,t)\psi(x,t).

If A^\hat A is an operator, we define the {\sl Weyl transform} of A^\hat A, denoted A(x,p)A(x,p), by

A(x,p)=∫−∞+∞ds e−ips/ℏ ⟨x+s/2∣A^∣x−s/2⟩.A(x,p) = \int_{-\infty}^{+\infty} ds \, e^{-ips/\hbar} \, \matrixelement{x+s/2}{\hat A}{x-s/2}.

Here the notation ∣x−s/2⟩\ket{x-s/2}, for example, means the eigenket of x^\hat x with eigenvalue x−s/2x-s/2. It is useful to think of A(x,p)A(x,p) as a function defined on the classical (x,p)(x,p) phase space which is in some sense the classical observable corresponding to the quantum operator A^\hat A.

(a) Show that if A(x,p)A(x,p) is the Weyl transform of operator A^\hat A, then A(x,p)†A(x,p)^\cc is the Weyl transform of A^†\hat A^\hc. In particular, this shows that the Weyl transform of a Hermitian operator is a real function on phase space.

(b) Show that if operators A^\hat A and B^\hat B have Weyl transforms A(x,p)A(x,p) and B(x,p)B(x,p), respectively, then

tr⁡(A^†B^)=∫dx dp2πℏ A(x,p)†B(x,p).\tr(\hat A^\hc \hat B) = \int {dx\,dp \over 2\pi\hbar} \, A(x,p)^\cc B(x,p).

Notice how the right hand side looks like the “scalar product” of two classical observables on phase space.

(c) Find the Weyl transforms of the following operators: 1 (the identity operator); x^\hat x; p^\hat p; x^p^\hat x \hat p; p^x^\hat p \hat x; p^2/2m+V(x^)\hat p^2/2m + V(\hat x).

(d) The classical probability density ρ(x,p)\rho(x,p) in phase space (see Classical Statistical Mechanics) describes a system whose dynamical state can only be described statistically. Since ρ(x,p)\rho(x,p) is a probability density, it is nonnegative, ρ(x,p)≥0\rho(x,p)\ge0, and it is normalized,

∫ρ(x,p) dx dp=1.\int \rho(x,p)\,dx\,dp=1.

The probability density in xx alone (call it F(x)F(x)) or in pp alone (call it G(p)G(p)) is obtained by integrating ρ\rho over the other variable (pp or xx),

∫ρ(x,p) dp=F(x),∫ρ(x,p) dx=G(p),\int \rho(x,p)\,dp = F(x), \qquad \int \rho(x,p)\,dx = G(p),

and each of these is normalized,

∫F(x) dx=1,∫G(p) dp=1.\int F(x)\,dx=1, \qquad \int G(p)\,dp=1.

The functions F(x)F(x) or G(p)G(p) are called marginal distributions in statistics.

It turns out that the Weyl transform W(x,p)W(x,p) of the density operator ρ^\hat\rho in quantum mechanics is very similar in its properties to the classical probability density ρ(x,p)\rho(x,p). The function W(x,p)W(x,p) is called the Wigner function. Because ρ^\hat\rho is Hermitian, the Wigner function is real.

For simplicity let us work with a pure state, so ρ^=∣ψ⟩⟨ψ∣\hat\rho=\ketbra{\psi}{\psi}. Express the marginal distribution

∫W(x,p) dp\int W(x,p)\,dp

in terms of the wave function ψ(x)\psi(x). Now let ϕ(p)\phi(p) be the momentum space wave function, and express the marginal distribution

∫W(x,p) dx\int W(x,p)\,dx

in terms of ϕ(p)\phi(p). What is the normalization integral,

∫W(x,p) dx dp?\int W(x,p)\,dx\,dp?

These results suggest that W(x,p)W(x,p) is a distribution function of particles in phase space whose statistics reproduces the statistics inherent in quantum measurement, apart from the normalization. Unlike a classical distribution function ρ(x,p)\rho(x,p), however, W(x,p)W(x,p) can take on negative values. These “negative probabilities” have no meaning in any statistical sense, but they only arise, in a certain sense, when we attempt to measure xx and pp simultaneously to a precision greater than that allowed by the uncertainty principle. The Wigner function can be used to express quantum statistical mechanics in a manner that is surprisingly similar to classical statistical mechanics, but one which is still fully quantum mechanical.

Problem 3.

Given that the range of the nuclear forces is approximately 10-13~cm, estimate the velocity of the proton or neutron in a deuteron (a bound state of a proton and a neutron). Compare this v/cv/c to the v/cv/c of the electron in the ground state of hydrogen. It is believed that nuclear forces are independent of the charge state of the nucleon; thus, two protons should feel the same nuclear force as two neutrons. Calculate the nuclear force between two protons at a distance of 10-13~cm and compare to the electrostatic force.

Problem 4.

Let f(x)f(\xvec) be an arbitrary function of x\xvec and g(p)g(\pvec) and arbitrary function of p\pvec. Evaluate the commutators, [pi,f(x)][p_i,f(\xvec)] and [xi,g(p)][x_i,g(\pvec)].

Problem 5.

The Large Hadron Collider is designed to reach an energy of 7~TeV (7×10127\times 10^{12}~eV). What is the smallest distance scale that can be probed at these energies? How does it compare to the Planck length, at which the effects of quantum gravity are expected to be important?