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.

21. Parity

Introduction

We have now completed our study of proper rotations in quantum mechanics, one of the important space-time symmetries. In these notes we shall examine parity, another space-time symmetry. In a latter set of notes we shall take up the third, which is time reversal. All these symmetries are special cases of Lorentz transformations, which are united in a relativistic treatment.

Spatial Inversion

Parity begins life as an operation on three-dimensional, physical space, represented in an inertial frame by the matrix

P=(−1000−1000−1)=−1.\begin{aligned} \Pmat=\begin{pmatrix}-1 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & -1 \\\end{pmatrix} = -\Imat. \end{aligned}

See also (16) and the discussion in Properties of Rotation Matrices; the Group O(3)O(3). This operation, which we call spatial inversion, inverts vectors through the origin, mapping x\xvec into −x-\xvec. The most important properties of P\Pmat for our purposes are

P2=1,\Pmat^2 = \Imat,

and

PRP−1=R,\Pmat \Rmat \Pmat^{-1} = \Rmat,

for all proper rotations R∈SO(3)\Rmat \in SO(3). These properties are trivial in view of (1); in particular, P\Pmat, being a multiple of the identity, commutes with all matrices, so PRP−1=RPP−1=R\Pmat \Rmat \Pmat^{-1} = \Rmat \Pmat \Pmat^{-1} = \Rmat.

The spatial inversion P\Pmat is an example of an improper rotation, that is, it belongs to O(3)O(3) but not SO(3)SO(3). In fact, every proper rotation is mapped into an improper one by multiplying by P\Pmat, and vice versa. See Fig.~\figr\classrot.5 to visualize the spaces of proper and improper rotations.

Parity in Quantum Mechanics

In this we introduce the operator π\pi, called the parity operator, which corresponds to the spatial inversion operation P\Pmat.

Recall that in the case of proper rotations, we sought and found operators U(R)U(\Rmat), acting on the Hilbert space of various quantum mechanical systems, that implement the effect of classical, proper rotations R\Rmat. The operators U(R)U(\Rmat) were required to satisfy certain postulates, including one that they be unitary, and another that they reproduce the multiplication law of the classical, proper rotations. See Postulates for Rotation Operators. We represent this mapping between classical rotations and rotation operators by

R↦U(R),\Rmat \mapsto U(\Rmat),

and we say that it specifies a unitary representation of the classical rotation group SO(3)SO(3).

Similarly, in the case of spatial inversion, we seek an operator π\pi acting on the Hilbert space of various quantum mechanical systems that is associated with the spatial inversion operation P\Pmat,

P↦π,\Pmat \mapsto \pi,

which is required to satisfy certain postulates. These are that π\pi be unitary,

π†π=1,\pi^\hc \pi =1,

and that it should satisfy

π2=1,\pi^2 = 1,

and

πU(R)π†=U(R)\pi U(\Rmat) \pi^\hc = U(\Rmat)

for all R∈SO(3)\Rmat \in SO(3). The latter two postulates are representations of the multiplication laws (2) and (3) at the level of spatial transformations. The first two postulates imply that π\pi is both unitary and Hermitian, so that

π=π†=π−1,\pi=\pi^\hc=\pi^{-1},

while the third postulate (8) is equivalent to

[π,J]=0,[\pi,\Jvec] = 0,

where J\Jvec is the angular momentum (the generator of the rotation operators U(R)U(\Rmat)). In short, π\pi is a scalar operator [see (8) and (10)]. These postulates do not uniquely determine the parity operator for all quantum mechanical systems, but they narrow the possibilities considerably so that it is easy to find a suitable definition of π\pi in specific cases.

Taken together with the properties of U(R)U(\Rmat), the postulates for π\pi amount to requiring that the set of operators U(R)U(\Rmat) and πU(R)\pi U(\Rmat) for R∈SO(3)\Rmat \in SO(3) form a unitary representation of the full rotation group O(3)O(3). That is, generalizing the map (4) to include improper rotations, we can say

π=U(P).\pi = U(\Pmat).

Parity for a Spinless Particle in Three Dimensions

Let us begin with the case of a spinless particle moving in three dimensions, for which the basis states can be taken to be the position eigenkets ∣x⟩\ket{\xvec} and the Hilbert space is isomorphic to the space of configuration space wave functions ψ(x)\psi(\xvec). In this case we defined the rotation operators by

U(R)∣x⟩=∣Rx⟩,U(\Rmat) \ket{\xvec} = \ket{\Rmat \xvec},

which we found satisfactory. See (3) and the discussion in Rotation Operators for Spatial Degrees of Freedom. In a similar spirit, we now guess that the definition of parity for such systems should be

π∣x⟩=∣Px⟩=∣−x⟩.\pi\ket{\xvec} = \ket{\Pmat\xvec} = \ket{-\xvec}.

This transformation law on basis kets implies that wave functions transform according to

ψ(x)text⟶πψ(−x).\psi(\xvec) \mathrel{\mathop{\\text{\rightarrowfill}}^{\pi}} \psi(-\xvec).

We must check that this definition satisfies the postulates (6)--(8).

First, by applying π\pi twice, we have two changes of sign,

π2∣x⟩=∣x⟩,\pi^2 \ket{\xvec} =\ket{\xvec},

so the postulate (7) is satisfied. Next, since

∫d3x ∣ψ(x)∣2=∫d3x ∣ψ(−x)∣2\int d^3\xvec \, |\psi(\xvec)|^2 = \int d^3\xvec \, |\psi(-\xvec)|^2

for all wave functions ψ\psi, the operator π\pi is unitary (it preserves the normalization of states), and the postulate (6) is satisfied. See Problem Paragraph(c). Finally, since

πU(R)∣x⟩=π∣Rx⟩=∣−Rx⟩=U(R)π∣x⟩\pi U(\Rmat) \ket{\xvec} = \pi \ket{\Rmat \xvec} = \ket{-\Rmat\xvec} = U(\Rmat)\pi\ket{\xvec}

we have πU(R)=U(R)π\pi U(\Rmat) = U(\Rmat)\pi, which is equivalent to the postulate (8). All postulates are satisfied, so we shall take (13) as the definition of π\pi for a spinless particle in three-dimensional space.

As for the last postulate (8) and its equivalent form (10), since the angular momentum for a spinless particle in three-dimensional space is the orbital angular momentum L=x×p\Lvec=\xvec \cross \pvec, we have

πLπ†=L.\pi \Lvec \pi^\hc = \Lvec.

This follows from the postulates, but can also be verified directly.

It is also of interest to compute the conjugation relations of other vector operators with respect to π\pi. We start with the position operator x\xvec. For the sake of the following demonstration, we will write x^{\hat{\mathbf{x}}} for the operator, and x\xvec for associated vectors of cc-numbers. In the configuration representation, the operator x^{\hat{\mathbf{x}}} means multiplication of the wave function by the position of the particle, that is,

(x^ψ)(x)=xψ(x).({\hat{\mathbf{x}}} \psi)(\xvec) = \xvec \psi(\xvec).

Now computing the effect of πx^π†\pi{\hat{\mathbf{x}}}\pi^\hc on an arbitrary wave function, we have

ψ(x)text⟶π†ψ(−x)text⟶x^xψ(−x)text⟶π−xψ(x)=−(x^ψ)(x),\psi(\xvec) \mathrel{\mathop{\\text{\rightarrowfill}}^{\pi^\hc}} \psi(-\xvec) \mathrel{\mathop{\\text{\rightarrowfill}}^{{\hat{\mathbf{x}}}}} \xvec\psi(-\xvec) \mathrel{\mathop{\\text{\rightarrowfill}}^{\pi}} -\xvec\psi(\xvec)=-({\hat{\mathbf{x}}}\psi)(\xvec),

where we use π†=π\pi^\hc = \pi.

Thus we have

πxπ†=−x,\pi \xvec \pi^\hc = -\xvec,

where we have dropped the hats so that now x\xvec represents the operator again. We see that conjugating the position operator x\xvec by π\pi converts it into −x-\xvec. This is reasonable, in view of the action of P\Pmat on points of space. Similarly, we find

πpπ†=−p.\pi\pvec\pi^\hc = -\pvec.

This is fairly obvious, since p=−iℏ∇\pvec=-i\hbar\del in the position representation, and ∇=∂/∂x\del=\partial/\partial \xvec should change sign when x\xvec does. Taking the cross product of (21) and (22), we obtain (18).

Parity for Spin Systems

Next let us consider just the spin degrees of freedom for a particle with spin, ignoring the spatial degrees of freedom for the time being. The Hilbert space is

Espin=span⁡{∣sm⟩,m=−s,…,+s},\ES_spin = \mathspan\{\ket{sm}, m=-s,\ldots,+s\},

with the standard angular momentum basis shown. We would like to find a reasonable definition of π\pi that satisfies the postulates (6)--(8). This means finding the action of π\pi on the basis states ∣sm⟩\ket{sm}. Since the angular momentum in this case is the spin vector S\Svec, the form (10) of the third postulate is

[π,S]=0.[\pi,\Svec]=0.

The states ∣sm⟩\ket{sm} are nondegenerate eigenstates of SzS_z on the Hilbert space (23), and (24) implies [π,Sz]=0[\pi,S_z]=0. Therefore by Theorem Theorem, ∣sm⟩\ket{sm} must also be an eigenstate of parity,

π∣sm⟩=ηm∣sm⟩,\pi\ket{sm} = \eta_m \ket{sm},

where as indicated the eigenvalue ηm\eta_m might depend on mm, as far as we know at this point. Theorem Theorem does not tell us the value of this eigenvalue, but it must be real since π\pi is Hermitian. Now by applying raising and lowering operators and using (24) again, we have

S±(π∣sm⟩)=πS±∣sm⟩=(s∓m)(s±m+1)ℏ π∣s,m±1⟩\begin{aligned} S_\pm (\pi\ket{sm}) &= \pi S_\pm\ket{sm} =\sqrt{(s\mp m)(s\pm m+1)} \hbar \, \pi\ket{s,m\pm1} \end{aligned}
=ηm(s∓m)(s±m+1)ℏ ∣s,m±1⟩,\begin{aligned} &=\eta_m \sqrt{(s\mp m)(s\pm m+1)} \hbar \, \ket{s,m\pm1}, \end{aligned}

or, on canceling the square roots,

π∣s,m±1⟩=ηm∣s,m±1⟩.\pi \ket{s,m\pm1} = \eta_m \ket{s,m\pm1}.

By comparing this to (25), we see that ηm=ηm±1\eta_m = \eta_{m\pm1}, or, by induction, that η\eta is independent of mm for all 2s+12s+1 values it can take on. This fact could have been seen from another standpoint, for since π\pi is a scalar operator, by the Wigner-Eckart theorem for scalars its eigenvalues must be independent of mm. See Secs. Energy Eigenstates in Isolated Systems-The 57Fe\ironfiftyseven Nucleus and the M"ossbauer Effect. Thus we shall write simply

π∣sm⟩=η∣sm⟩.\pi\ket{sm} = \eta \ket{sm}.

In fact, since the states ∣sm⟩\ket{sm} form a basis for Espin\ES_spin, all kets in this space are eigenstates of π\pi with eigenvalue η\eta, that is, π\pi is η\eta times the identity operator.

Now using the fact that π2=1\pi^2=1, we find that η2=1\eta^2=1, or η=±1\eta=\pm1, and we have

π∣sm⟩=±∣sm⟩.\pi\ket{sm} = \pm \ket{sm}.

The postulates required of π\pi reduce its definition to a choice of a sign. This sign is the same for all vectors in the Hilbert space Espin\ES_spin of (23), so it is really a characteristic of the particle, that is, a constant like the spin quantum number ss.

Can this sign be determined by any further arguments or by appeal to experimental results? Within nonrelativistic quantum mechanics, the answer is no, the final choice of sign is without physical consequences and can be made arbitrarily. The simplest choice is η=1\eta=1, and for as long as we remain with the nonrelativistic theory, we shall assume that the parity operator satisfies

π∣sm⟩=∣sm⟩.\pi\ket{sm} = \ket{sm}.

That is, the parity operator has no effect on the spin state of the particle.

In relativistic quantum mechanics, however, it becomes possible to create and destroy particles. In some cases, an initial state of a definite parity decays or otherwise evolves into a state in which a particle has been created or destroyed, and in order to account for conservation of parity it is necessary to attribute a certain parity to a particle because of the kind of particle it is. This is called intrinsic parity, and it is equivalent to making one of the two choice η=±1\eta=\pm1 for each type of particle in the process. Some details of this assignment depend on the type of particle, fermion or boson. (It turns out that for fermions, the intrinsic parity of a particle-antiparticle pair is odd; it is arbitrary whether the particle is interpreted as having even intrinsic parity and the antiparticle odd, or vice versa. See 48. Covariance of the Dirac Equation.) Not all interactions conserve parity, but for those that do, it is necessary in relativistic interactions to take into account the intrinsic parities of particles. Then the total parity of a state is the product of the intrinsic parities of the particles times the parities of the spatial parts of the wave functions. For example, the photon turns out to have negative intrinsic parity. This is related to Laporte’s rule, which we discuss presently, which states that the parity of an atom changes when a photon is emitted in an electric dipole transition. The parity of the atom changes, but the parity of the entire system, including the electromagnetic field, is conserved, as we see when we include the intrinsic parity of the photon. The intrinsic parity of particles cannot be understood very well outside of the relativistic theory, but we will have more to say about it 48. Covariance of the Dirac Equation.

Returning to nonrelativistic quantum mechanics, if we include the spatial degrees of freedom of a particle, so that the Hilbert space is spanned by the basis kets

∣x⟩⊗∣sm⟩=∣x,m⟩,\ket{\xvec}\otimes \ket{sm} = \ket{\xvec, m},

then the parity operator is defined by

π∣x,m⟩=∣−x,m⟩.\pi\ket{\xvec,m} = \ket{-\xvec,m}.

This is equivalent to the action on the wave function,

ψm(x)text⟶πψm(−x).\psi_m(\xvec) \mathrel{\mathop{\\text{\rightarrowfill}}^{\pi}} \psi_m(-\xvec).

As we say, parity is a purely spatial operator, and does not affect the spin (in the nonrelativistic theory).

Transformation Properties Under Parity

Returning to the conjugation relations (18), (21) and (22), we see that some vectors (for example, x\xvec and p\pvec) change sign when conjugated by parity, while some (for example, L\Lvec and S\Svec) do not. All these vectors are vector operators according to the definitions (15) or (21), that is, they all transform in the same way under conjugation by proper rotations, but their transformation laws under improper rotations differ. Thus we can take the set of vector operators and break it down according to how the vectors transform under parity. If V\Vvec is a vector operator and

πVπ†=±V,\pi \Vvec \pi^\hc = \pm \Vvec,

then we will say that V\Vvec is a true vector or polar vector if the sign is -1, or a pseudovector or axial vector if the sign is +1. Thus, x\xvec and p\pvec are true vectors, while L\Lvec and S\Svec (more generally, any angular momentum J\Jvec) are pseudovectors. Other vectors in physics can be classified similarly, for example, the electric field E\Evec is a true vector, while the magnetic field B\Bvec is a pseudovector.

Similarly, we can classify scalars. If KK is a scalar operator according to the definitions (8) or (10), and if

πKπ†=±K,\pi K \pi^\hc = \pm K,

then we will say that KK is a true scalar if the sign is +1 or a pseudoscalar if the sign is -1 (notice the signs are reversed compared to the terminology for a vector). An example of a true scalar is x⋅p\xvec\cdot\pvec, while an example of a pseudoscalar is p⋅S\pvec\cdot\Svec.

Hamiltonians Invariant Under Parity

We argued previously that the Hamiltonian for an isolated system must be a scalar, that is, invariant under proper rotations, because the energy could not depend on the orientation. See Energy Eigenstates in Isolated Systems. But are these Hamiltonians scalars or pseudoscalars or a combination of the two? That is, are they also invariant under parity?

As an example, consider a central force Hamiltonian,

H=p22m+V(r).H={p^2\over 2m} + V(r).

This certainly commutes with parity, for although p\pvec is mapped into −p-\pvec, the kinetic energy contains p⋅p\pvec\cdot\pvec, which is invariant. Also, the potential energy is a function of the radius rr, the magnitude of the vector x\xvec, which does not change when x\xvec changes its sign. Parity is a good quantum number for all central force Hamiltonians.

As we know, central force Hamiltonians usually represent a two-particle system interacting by means of forces derivable from a potential that is a function only of the distance between the two particles. See Two-Body Central Force Motion. All these Hamiltonians commute with parity. This generalizes to any number of particles interacting by means of forces derivable from potentials depending only on the distances between the particles, such as the electrostatic Coulomb approximation to the interactions of charged particles. This covers a wide variety of approximate Hamiltonians for atoms, molecules, and other systems. What happens when we include electromagnetic and relativistic effects? As it turns out, such Hamiltonians still commute with parity, as long as the electromagnetic fields are internally generated (that is, the system is still isolated). For example, one of the magnetic and relativistic corrections in a single electron atom is the spin-orbit term, which is proportional to L⋅S\Lvec\cdot\Svec. But by the rules above, this is a true scalar (the dot product of two pseudovectors). What about when we include the emission and absorption of radiation, that is, of photons? In this case, too, parity is conserved, as long as we include the parity of the electromagnetic field (that is, of the photons) in our account. We summarize these facts by saying that the electromagnetic forces conserve parity.

The strong forces, responsible for the interactions of nucleons inside a nucleus as well as interactions of other hadronic particles such as the π\pi-meson, also conserve parity. This fact was originally established by experiments involving nuclei and particles, which showed that in strong interactions the parity of the initial state is always equal to the parity of the final state. In reckoning the parity of states, however, it is necessary to take into account the intrinsic parities of various particles, for example, the π\pi-meson (all three charge states of it) has negative intrinsic parity (just like the photon).

For these reasons, it used to be believed that parity conservation was a fundamental law of nature, much like the rotational invariance of isolated systems. This belief ended abruptly in 1956, however, when Lee and Yang suggested that parity might be violated in weak interactions. They also suggested experiments to test the idea that were quickly carried out. These showed that, indeed, parity is violated in β\beta-decay, a prime example of a weak interaction.

The weak forces are present in ordinary atomic and nuclear interactions, but, as the name implies, they are very small compared to the electromagnetic and strong interactions, at least at low energy. The weak forces become relatively stronger as the energy is increased, and it is currently believed that at energies much higher than now available in accelerators the strong, electromagnetic and weak forces all become comparable in strength. Considering only low energy interactions, however, we can say that the weak interactions are extremely small compared to the electromagnetic and nuclear interactions, in fact, they are so small that they are quite difficult to detect even if one is looking for them. There has been some effort in the Berkeley Physics department to detect the effects of weak interactions in atomic physics, and these experiments are difficult. For most ordinary purposes at ordinary energies, the weak interactions can be ignored and all Hamiltonians for isolated systems conserve parity to an extremely high degree of accuracy.

It is, however, easy to create a system that does not conserve parity, simply by placing it in an external field. Such a system is no longer isolated, and parity need not be conserved. For example, suppose a particle of charge qq moving in a central force field is placed in an external, uniform electric field E0\Evec_0. Then the Hamiltonian is

H=p22m+V(r)−qx⋅E0.H={p^2 \over 2m} + V(r) - q\xvec\cdot\Evec_0.

This does not commute with parity because the last term changes sign upon conjugation by parity.

In such cases, however, parity conservation is restored if we enlarge our definition of “the system” to include the charges creating the electric field. For example, if E0\Evec_0 is created by positive and negative charges arrayed along the parallel plates of a capacitor, and if these charges are included in “the system,” then π\pi will swap the positive and negative charges and reverse the direction of E0\Evec_0. Since x\xvec also changes sign, the whole Hamiltonian will once again be invariant under parity.

Parity and Energy Eigenstates in Isolated Systems

Suppose we have an isolated system with Hamiltonian HH. As we have explained, in most applications HH commutes with parity to an extremely high degree of accuracy, so let us assume that πHπ†=H\pi H\pi^\dagger =H, that is, [π,H]=0[\pi,H]=0.

Since the system is isolated, HH commutes with all proper rotations, and is a scalar operator according to the definition (8). Then Theorem Theorem states that the energy eigenspaces consist of one or more irreducible subspaces under rotations, and the addendum says that, with few exceptions, they consist of precisely one irreducible subspace. So let us consider this generic case, for which each energy eigenvalue is associated with a definite jj value and has an energy eigenspace that is 2j+12j+1-fold degenerate.

As in 13. Representations of the Angular Momentum Operators and Rotations{} and \transfop{} let Sjm\SS_{jm} be the eigenspace of operators J2J^2 and J3J_3 with quantum numbers jj and mm, and let ∣ψ⟩\ket{\psi} be a bound energy eigenstate in the stretched subspace Sjj\SS_{jj} with eigenvalue EE, so that

J2∣ψ⟩=j(j+1)ℏ2 ∣ψ⟩,J3∣ψ⟩=jℏ ∣ψ⟩,H∣ψ⟩=E∣ψ⟩,J^2\ket{\psi}=j(j+1)\hbar^2\,\ket{\psi}, \quad J_3\ket{\psi}=j\hbar\,\ket{\psi}, \quad H\ket{\psi}=E\ket{\psi},

Under our assumptions of genericity, HH restricted to S\SS is nondegenerate, that is, ∣ψ⟩\ket{\psi} is the only state in Sjj\SS_{jj} with the energy EE.

Now π\pi commutes with J\Jvec and hence with J2J^2 and JzJ_z, and so it can be restricted to Sjj\SS_{jj}. And since ∣ψ⟩\ket{\psi} is a nondegenerate eigenstate of HH and since π\pi commutes with HH, according to Theorem Theorem, ∣ψ⟩\ket{\psi} must also be an eigenstate of parity,

π∣ψ⟩=e∣ψ⟩,\pi\ket{\psi}=e\ket{\psi},

where e=±1e=\pm1. Theorem Theorem does not tell us what the value of ee is.

Now it is easy to show that all the states obtained by applying lowering operators to ∣ψ⟩\ket{\psi} are also eigenstates of parity with the same eigenvalue ee as the stretched state. That is,

π(J−)r∣ψ⟩=(J−)rπ∣ψ⟩=e(J−)r∣ψ⟩,\pi (J_-)^r \ket{\psi} = (J_-)^r \pi\ket{\psi} = e(J_-)^r\ket{\psi},

for r=1,2,…,2jr=1,2,\ldots,2j. That is, all 2j+12j+1 degenerate energy eigenstates in the irreducible subspace have the same parity.

Thus the parity value can be associated with the energy level itself, just like the jj value. This is why the three nuclear states shown in Fig.~\figr\transfop.1 are assigned definite values of parity (all three, 57Fe\ironfiftyseven, 57Fe†\ironfiftyseven^\cc and 57Fe††\ironfiftyseven^{\cc\cc}, have negative parity).

In the case of nuclei, the energies, angular momenta (or spins) and parities of energy levels are determined experimentally, since theoretical calculations of these quantities are difficult. In atoms and molecules, theory is in better shape, and these quantities can be calculated with greater or lesser difficulty.

In the nongeneric case of the electrostatic model of hydrogen, it is not true that all energy eigenstates are eigenstates of parity. It is true that in central force motion each angular momentum value ℓ\ell is associated with the definite parity (−1)ℓ(-1)^\ell [see (52)], but since there are degeneracies of energy that cross angular momentum values, it is possible to form linear combinations of energy eigenstates of the same energy but opposite parity. For example, the 2s2s and 2p2p levels in hydrogen have the same energy but opposite parities. These facts play an important role in understanding the Stark effect, which is discussed in 24. The Stark Effect in Hydrogen and Alkali Atoms.

Parity in One-Dimensional Systems

Consider a one-dimensional Hamiltonian,

H=p22m+V(x).H={p^2\over 2m} + V(x).

This commutes with parity if and only if V(x)=V(−x)V(x)=V(-x), that is, if the potential is symmetric about x=0x=0. Now in one-dimensional systems with a localized potential, the bound energy eigenstates are nondegenerate because the wave functions must vanish at x=±∞x=\pm\infty. See Degeneracies in One Dimension. Therefore by Theorem Theorem, the energy eigenstates must also be eigenstates of parity. That is, the bound energy eigenfunctions satisfy ψ(x)=±ψ(−x)\psi(x) = \pm \psi(-x), for some choice of sign that depends on the eigenfunction. For example, in the harmonic oscillator, which has the symmetric potential V(x)=mω2x2/2V(x) = m\omega^2 x^2/2, the nn-th energy eigenfunction has parity (−1)n(-1)^n. See \harmosc.6.

Considerations of Symmetry Raised by Parity

Parity illustrates some issues that arise whenever a Hamiltonian possesses a symmetry. Some of these have already been discussed in connection with rotations, but parity is a simpler symmetry so some of the ideas are easier to appreciate.

Whenever a Hamiltonian possesses a symmetry and we wish to diagonalize it or otherwise study its properties, it is a good idea to diagonalize the symmetry operators first. This is because symmetry operators are usually easier to diagonalize than Hamiltonians, and their eigenbases, sometimes called symmetry-adapted bases, are especially convenient for further study of the Hamiltonian. For example, in the case of rotational invariance, the symmetry adapted basis is the standard angular momentum basis ∣γjm⟩\ket{\gamma jm}, an eigenbasis of J2J^2 and JzJ_z.

Similarly, if a Hamiltonian commutes with parity, it may be convenient to study it in a basis which is an eigenbasis of parity. Since π2=1\pi^2=1, the eigenvalues ee satisfy e2=1e^2=1, or e=±1e=\pm1. Thus the Hilbert space is decomposed into two subspaces, the eigenspaces of parity with e=±1e=\pm1. We can call these the even and odd subspaces,

E=Eeven⊕Eodd.\ES = \ES_even \oplus \ES_odd.

Then a symmetry adapted basis consists of an orthonormal basis inside Eeven\ES_even plus an orthonormal basis inside Eodd\ES_odd. We can write the basis vectors as ∣γe⟩\ket{\gamma e}, where

π∣γe⟩=e∣γe⟩,\pi \ket{\gamma e} = e \ket{\gamma e},

so that ee indicates which subspace the vectors lie in, and where γ\gamma is a label of the orthonormal basis vectors inside Eeven\ES_even (for e=+1e=+1) or inside Eodd\ES_odd (for e=−1e=-1). This notation is similar to the notation ∣γjm⟩\ket{\gamma jm} which we used for a standard angular momentum basis, where the index γ\gamma plays a similar role.

It is possible to write down projection operators that project onto the even and odd subspaces in terms of the parity operator itself. These are

P±=1±π2.P_\pm = {1 \pm \pi \over 2}.

It is easy to show that the operators P±P_\pm satisfy the requirements of a projection operator (see Projection Operators), and that when applied to an arbitrary state ∣ψ⟩\ket{\psi} (not necessarily an eigenstate of parity), they project it onto one or the other of the two subspaces, producing an eigenstate of parity of the eigenvalue indicated. For example, in a one-dimensional problem, if ψ(x)\psi(x) is an arbitrary wave function (not necessarily even or odd), then

(P+ψ)(x)=12[ψ(x)+ψ(−x)](P_+\psi)(x) = {1\over 2}[\psi(x)+\psi(-x)]

and

(P−ψ)(x)=12[ψ(x)−ψ(−x)](P_-\psi)(x) = {1\over 2}[\psi(x)-\psi(-x)]

are the even and odd parts of ψ(x)\psi(x), respectively.

The matrix elements of the Hamiltonian in the symmetry adapted basis are noteworthy. Assuming the Hamiltonian commutes with parity, we have

⟨γ′e′∣H∣γe⟩=⟨γ′e′∣π†Hπ∣γe⟩=e′e ⟨γ′e′∣H∣γe⟩,\matrixelement{\gamma'e'}{H}{\gamma e} =\matrixelement{\gamma'e'}{\pi^\hc H \pi}{\gamma e} = e'e \,\matrixelement{\gamma'e'}{H}{\gamma e},

which implies that the matrix element vanishes unless e′e=1e'e=1. But this is equivalent to e′=ee'=e, so the Hamiltonian is diagonal in the parity quantum number ee, and we have

⟨γ′e′∣H∣γe⟩=δe′e Cγ′γe,\matrixelement{\gamma'e'}{H}{\gamma e} = \delta_{e'e} \, C^e_{\gamma'\gamma},

where Cγ′γeC^e_{\gamma'\gamma} just indicates the indices the nonvanishing matrix elements can depend on. This may be compared to (55), which applies to Hamiltonians that are invariant under proper rotations. The fact that the matrix element vanishes if e′≠ee' \ne e is sometimes stated in words by saying that the Hamiltonian does not “connect” states of different parity.

In many cases it is necessary to diagonalize a Hamiltonian numerically, by setting up a matrix for it in some basis. There is the question, however, of what basis to use. If the Hamiltonian respects a symmetry then there is great advantage in using a symmetry adapted basis, since it reduces the size of the matrices that must be diagonalized. For example, in the parity-adapted basis ∣γe⟩\ket{\gamma e}, to diagonalize a Hamiltonian that commutes with parity we must diagonalize the two matrices Cγ′γeC^e_{\gamma'\gamma} with rows and columns indexed by γ′\gamma' and γ\gamma, for e=±1e=\pm1. Effectively we are diagonalizing the Hamiltonian separately in each of the subspaces Eeven\ES_even and Eodd\ES_odd. If we arrange the basis vectors so that all the even vectors ∣γ,+1⟩\ket{\gamma,+1} come first, and all the odd vectors ∣γ,−1⟩\ket{\gamma,-1} come second, then the Hamiltonian matrix has the form,

where the upper (or even-even) block is the matrix Cγ′γ+1C^{+1}_{\gamma'\gamma}, and the lower (or odd-odd) block is the matrix Cγ′γ−1C^{-1}_{\gamma'\gamma}, and where XX indicates a non-zero block.

There are many numerical algorithms for diagonalizing a matrix, but a common one (the Householder algorithm) involves a computational effort that scales as N3N^3, where NN is the size of the matrix. Thus the effort to diagonalize two matrices of size N/2N/2 is 2(N/2)3=N3/42(N/2)^3 = N^3/4, or about 1/41/4 the effort of diagonalizing a single N×NN\times N matrix. This is the improvement in computational efficiency when using a symmetry adapted basis for a Hamiltonian that commutes with parity, as compared to diagonalizing the Hamiltonian in an arbitrary (nonsymmetric) basis.

Parity and Selection Rules

Parity gives rise to important selection rules which we will illustrate for the case of dipole transitions in an electrostatic model of a single-electron atom. The energy eigenfunctions are the usual ones in central force motion,

ψnℓm(r,θ,ϕ)=Rnℓ(r)Yℓm(θ,ϕ),\psi_{n\ell m}(r,\theta,\phi) = R_{n\ell}(r) Y_{\ell m}(\theta,\phi),

which we shall write in ket language as ∣nℓm⟩\ket{n\ell m}. The matrix elements relevant to electric dipole transitions are

⟨n′ℓ′m′∣x∣nℓm⟩.\matrixelement{n'\ell'm'}{\xvec}{n\ell m}.

We need to find the action of parity on the eigenstates ∣nℓm⟩\ket{n\ell m}. Since the radius rr is invariant under the replacement x→−x\xvec \to -\xvec, only the YℓmY_{\ell m} part of the wave function is affected by parity. The YℓmY_{\ell m}'s have the property that if they are multiplied by rℓr^\ell, where rr is the radius, then the result is a homogeneous polynomial in the Cartesian coordinates xx, yy and zz of degree ℓ\ell. See YℓmY_{\ell m}'s in Cartesian Coordinates. But since xx, yy and zz all change sign under parity, and since the factor rℓr^\ell is invariant, the YℓmY_{\ell m} changes by (−1)ℓ(-1)^\ell under parity. We summarize this by writing

π∣nℓm⟩=(−1)ℓ∣nℓm⟩.\pi\ket{n\ell m} = (-1)^\ell \ket{n\ell m}.

The energy eigenfunctions (50) are also eigenfunctions of parity, with eigenvalue (−1)ℓ(-1)^\ell. This is often a useful fact. Notice that the eigenvalues do not depend on mm. This was to be expected, since π\pi is a scalar operator, and according to the Wigner-Eckart theorem, the eigenvalues of a scalar operator do not depend on mm. See Energy Eigenstates in Isolated Systems.

Now using the fact that x\xvec is odd under parity, we can make some statements about the matrix element (51). That is,

⟨n′ℓ′m′∣x∣nℓm⟩=−⟨n′ℓ′m′∣π†xπ∣nℓm⟩=−(−1)ℓ+ℓ′⟨n′ℓ′m′∣x∣nℓm⟩.\matrixelement{n'\ell'm'}{\xvec}{n\ell m} = -\matrixelement{n'\ell'm'}{\pi^\hc\xvec\pi}{n\ell m} = -(-1)^{\ell+\ell'} \matrixelement{n'\ell'm'}{\xvec}{n\ell m}.

Thus, the matrix element (51) vanishes unless −(−1)ℓ+ℓ′=1-(-1)^{\ell+\ell'} = 1, or

Δℓ=textodd,\Delta \ell = \\textodd,

where Δℓ=ℓ−ℓ′\Delta \ell = \ell-\ell'.

We see that the parity of the atomic state must change under an electric dipole transition. This selection rule is known as Laporte’s rule. Initially it was observed experimentally, and later it was explained by Wigner on the basis of parity. As noted above, Laporte’s rule is equivalent to the statement that the intrinsic parity of the photon is negative.

In the case of the matrix element (51), the Wigner-Eckart theorem gives the selection rule Δℓ=0,±1\Delta\ell = 0,\pm1, that is, these are the values of Δℓ\Delta \ell allowed by rotational invariance alone. But parity excludes the value Δℓ=0\Delta\ell =0, so the final selection rule is Δℓ=±1\Delta\ell=\pm1.

Problems

Problem 1.

This is essentially Sakurai’s problem 4.4, p.~282 (revised edition), or problem~4.4, p.~300 (second edition), but I’ve rewritten it to make it more clear.

We did not discuss the spin angular functions in class, but these are spinor functions on the unit sphere that arise when we combine orbital and spin angular momentum for a central force problem for a spinning particle. Here we will take the case s=12s=\frac{1}{2} (for example, in hydrogen). Let ∣nℓmℓ⟩\ket{n\ell m_\ell} be ket language for the wave function Rnℓ(r)Yℓmℓ(θ,ϕ)R_{n\ell}(r) Y_{\ell m_\ell}(\theta,\phi), the solution of the Schr"odinger equation for a spinless particle in a central force field, and let ∣sms⟩\ket{s m_s} be the usual spin states (here s=12s=\frac{1}{2} and ms=±12m_s = \pm \frac{1}{2}). We distinguish between mℓm_\ell and msm_s, the two types of magnetic quantum numbers. We multiply the wave functions times the spin functions and form linear combinations with the Clebsch-Gordan coefficients to get eigenstates of J2J^2 and JzJ_z, where J=L+S\Jvec=\Lvec+\Svec. These are

∣nℓjmj⟩=∑mℓ,ms∣nℓmℓ⟩∣sms⟩⟨ℓsmℓms∣jmj⟩.\ket{n\ell j m_j} = \sum_{m_\ell,m_s} \ket{n\ell m_\ell} \ket{s m_s} \braket{\ell s m_\ell m_s}{j m_j}.

The spatial wave functions ψnℓmℓ(x)=Rnℓ(r)Yℓmℓ(θ,ϕ)\psi_{n\ell m_\ell}(\xvec) =R_{n\ell}(r) Y_{\ell m_\ell}(\theta,\phi) factor into a radial part times an angular part. Let us write this in ket language as ∣nℓmℓ⟩=∣nℓ⟩∣ℓmℓ⟩\ket{n\ell m_\ell} = \ket{n \ell} \ket{\ell m_\ell}. The factor Rnℓ(r)R_{n\ell}(r) or ∣nℓ⟩\ket{n\ell} is the same for all terms in the sum above, so it can be taken out and what is left is a two-component spinor that depends only on the angles. This is what in wave function language Sakurai calls Yℓjmj{\cal Y}^{j m_j}_\ell; in ket language we will write

∣ℓjmj⟩=∑mℓms∣ℓmℓ⟩∣sms⟩⟨ℓsmℓms∣jmj⟩.\ket{\ell j m_j} = \sum_{m_\ell m_s} \ket{\ell m_\ell} \ket{s m_s} \braket{\ell s m_\ell m_s}{j m_j}.

(a) For the case ℓ=0\ell=0, j=12j=\frac{1}{2}, mj=12m_j=\frac{1}{2}, write out the two-component spinor Yℓjmj{\cal Y}^{j m_j}_\ell as functions of (θ,ϕ)(\theta,\phi).

(b) Multiply this by σ⋅x\sigmavec\cdot\xvec, and express the result as a linear combination of other spin angular functions Yℓjmj{\cal Y}^{j m_j}_\ell.

(c) Certain values of jj, mjm_j and ℓ\ell occur in the sum, and certain others do not. Use symmetry principles to explain why the values that occur are allowed and the others are not.

Hint: It may help to think of three kinds of rotations: spin rotations, orbital rotations, and total (spin plus orbital) rotations.

Problem 2.

This is Sakurai, problem 4.5, p.~282 (revised edition), or problem~4.5, p.~300 (second edition).

Because of weak (neutral-current) interactions, there is a parity-violating potential between the atomic electron and the nucleus as follows:

V=λ[δ(x) S⋅p+S⋅p δ(x)],V=\lambda[\delta(\xvec)\,\Svec\cdot\pvec + \Svec\cdot\pvec\, \delta(\xvec)],

where S\Svec and p\pvec are the spin and momentum operators of the electron, and the nucleus is assumed to be situated at the origin. As a result, the ground state of an alkali atom, usually characterized by ∣nℓjm⟩\ket{n\ell jm}, actually contains very tiny contributions from other eigenstates as follows:

∣nℓjm⟩→∣nℓjm⟩+∑n′ℓ′j′m′Cn′ℓ′j′m′ ∣n′ℓ′j′m′⟩,\ket{n\ell jm} \to \ket{n\ell jm} + \sum_{n'\ell'j'm'} C_{n'\ell'j'm'}\, \ket{n'\ell'j'm'},

On the basis of symmetry considerations alone, what can you say about the values of (n′ℓ′j′m′)(n'\ell'j'm') that give rise to nonvanishing contributions? Suppose the radial wave functions and energy levels are all known. Indicate how you can calculate Cn′ℓ′j′m′C_{n'\ell'j'm'}. Are there any further restrictions on (n′ℓ′j′m′)(n'\ell'j'm') that go beyond those due to symmetry considerations?

Problem 3.

This problem concerns a system whose symmetry is an improper rotation.

A collection of N≥3N\ge3 nuclei is distributed in the xx-yy plane. They do not lie on a line. Nucleus ii has position Ri=(Rix,Riy,0)\Rvec_i=(R_{ix},R_{iy},0). The positions of the nuclei are fixed. An electron moves in the field of these nuclei.

Let AA be the reflection in the xx-yy plane. It is represented by the matrix

A=(10001000−1),\begin{aligned} \Amat=\begin{pmatrix}1 & 0 & 0\\ 0 & 1 & 0\\ 0 & 0 & -1\\\end{pmatrix}, \end{aligned}

which maps (x,y,z)(x,y,z) into (x,y,−z)(x,y,-z). The electrostatic potential Φ(x)\Phi(\xvec) of the nuclei is invariant under reflection in the xx-yy plane, that is,

Φ(x,y,z)=Φ(x,y,−z),\Phi(x,y,z) = \Phi(x,y,-z),

or

Φ(Ax)=Φ(x)\Phi(\Amat\xvec) = \Phi(\xvec)

for short.

The x-y plane, viewed from the side. The spots in the x-y plane are the nuclei. \xvec is a typical point, and \Amat\xvec is the reflection of this point through the x-y plane. The figure illustrates Eqs.  and .

Figure 1:The xx-yy plane, viewed from the side. The spots in the xx-yy plane are the nuclei. x\xvec is a typical point, and Ax\Amat\xvec is the reflection of this point through the xx-yy plane. The figure illustrates Eqs. (62) and (63).

From E=−∇Φ\Evec=-\del\Phi and (60) or (61) it follows that

(Ex(x,y,−z)Ey(x,y,−z)Ez(x,y,−z))=(Ex(x,y,z)Ey(x,y,z)−Ez(x,y,z)),\begin{pmatrix}E_x(x,y,-z)\\ E_y(x,y,-z)\\ E_z(x,y,-z)\\\end{pmatrix} = \begin{pmatrix}E_x(x,y,z)\\ E_y(x,y,z)\\ -E_z(x,y,z)\\\end{pmatrix},

or

E(Ax)=AE(x)\Evec(\Amat\xvec)=\Amat\Evec(\xvec)

for short. This also follows by inspection of Figure 1.

(a) Every improper rotation can be written as PR(n^,θ)\Pmat \Rmat({\hat{\mathbf{n}}},\theta), where P\Pmat is given by (1) and R(n^,θ)\Rmat({\hat{\mathbf{n}}},\theta) is a proper rotation in axis-angle form. Write

A=PR(n^,θ)\Amat=\Pmat\Rmat({\hat{\mathbf{n}}},\theta)

and identify n^{\hat{\mathbf{n}}} and θ\theta.

(b) Is Φ(x)\Phi(\xvec) a scalar operator? Is E(x)\Evec(\xvec) a vector operator?

(c) In the electrostatic approximation the Hamiltonian for the electron is

H0=p22m+V(x),H_0={\pvec^2\over 2m} + V(\xvec),

where

V(x)=−eΦ(x).V(\xvec) = -e\Phi(\xvec).

Define a unitary operator WW that acts on the Hilbert space of the electron and that corresponds to A\Amat above in some reasonable way. I suggest you denote the usual parity operator by Π\Pi, to avoid confusion with the angle π\pi. You may find it convenient to refer to different types of rotation operators,

Uo(n^,θ)=e−iθn^⋅L/ℏ(orbital),Us(n^,θ)=e−iθn^⋅S/ℏ(spin),Ut(n^,θ)=e−iθn^⋅J/ℏ(total),\begin{aligned} U_o({\hat{\mathbf{n}}},\theta) &= e^{-i\theta{\hat{\mathbf{n}}}\cdot\Lvec/\hbar} \qquad {\rm(orbital)}, \\ U_s({\hat{\mathbf{n}}},\theta) &= e^{-i\theta{\hat{\mathbf{n}}}\cdot\Svec/\hbar} \qquad {\rm(spin)}, \\ U_t({\hat{\mathbf{n}}},\theta) &= e^{-i\theta{\hat{\mathbf{n}}}\cdot\Jvec/\hbar} \qquad {\rm(total)}, \\ \end{aligned}

as needed, where J=L+S\Jvec=\Lvec+\Svec.

Work out the conjugation relations,

W†(xyz)W,andW†(pxpypz)W,W^\dagger \begin{pmatrix}x\\ y\\ z\\\end{pmatrix} W, \qquad and \qquad W^\dagger \begin{pmatrix}p_x \\ p_y\\ p_z\\\end{pmatrix} W,

or simply W†xWW^\dagger \xvec W and W†pWW^\dagger \pvec W for short. Note that for any function f(x)f(\xvec),

W†f(x)W=f(W†xW).W^\dagger f(\xvec)W =f\bigl(W^\dagger \xvec W\bigr).

The WW you define should satisfy

W†H0W=H0.W^\dagger H_0 W=H_0.

If it does not, go back and change the definition so that it does. Notice that since you defined WW to be unitary, this is equivalent to [W,H0]=0[W,H_0]=0.

(d) So the system is invariant under reflection in the plane, and we have found an operator corresponding to that reflection that commutes with the electrostatic Hamiltonian H0H_0. Now we will see if the symmetry persists when we include relativistic corrections.

The spin-orbit Hamiltonian worked out in 25. Fine Structure in Hydrogen and Alkali Atoms{} [see (13)] applies only in the case of a radial electric field. The more general expression [see (49)] is

HSO=e4m2c2[E(x)×p−p×E(x)]⋅S,H_{SO} = {e\over 4m^2c^2} \bigl[\Evec(\xvec)\cross\pvec - \pvec\cross\Evec(\xvec)\bigr] \cdot\Svec,

where the quantity in the square brackets is the sum of an operator and its Hermitian conjugate, to make HSOH_{SO} Hermitian.

Does the operator WW you defined in part~(c) commute with HSOH_{SO}? If not, modify it so that it does, and call the new operator XX. Make sure it is still unitary. Depending on your choice of WW it is possible that no modification is necessary, then X=WX=W. If X≠WX\ne W, reconsider the commutation relations (68). Also work out the commutation relations

X†(Ex(x)Ey(x)Ez(x))XandX†(SxSySz)X,X^\dagger\begin{pmatrix}E_x(\xvec)\\ E_y(\xvec)\\ E_z(\xvec)\\\end{pmatrix}X \qquad and \qquad X^\dagger\begin{pmatrix}S_x \\ S_y\\ S_z\\\end{pmatrix}X,

or X†E(x)XX^\dagger \Evec(\xvec)X and X†SXX^\dagger\Svec X for short. Hint: notice that [(RA)×(RB)]⋅(RC)=(A×B)⋅C[(\Rmat\Avec)\cross(\Rmat\Bvec)]\cdot(\Rmat\Cvec) = (\Avec\cross\Bvec)\cdot\Cvec, for all proper rotations R\Rmat and all vectors A\Avec, B\Bvec, C\Cvec. See Prob. Paragraph.

(e) In quantum mechanics we like Hermitian operators that commute with the Hamiltonian, because they possess simultaneous eigenbases, and for other reasons. The operator XX you found in part~(d) should be unitary. Find a simple expression for X2X^2.

Is the XX you found in part~(d) Hermitian? If not, make a simple modification to it to create a new operator, call it KK, that is both unitary and Hermitian and that commutes with HH. If no modification is necessary, then KK will be the same as XX obtained in part~(d). If KK is both unitary and Hermitian, then its its eigenvalues must be ±1\pm1. Show that this is so.

Let the electron wave function be

ψ=(ψ+(x,y,z)ψ−(x,y,z)).\psi = \begin{pmatrix}\psi_+(x,y,z)\\ \psi_-(x,y,z)\\\end{pmatrix}.

Find the ± components of the wave function KψK\psi, and find the conditions that ψ+\psi_+ and ψ−\psi_- must satisfy in order that ψ\psi be an eigenfunction of KK with eigenvalues ±1\pm1.

Thus the Hilbert space breaks up into two orthogonal eigenspaces of KK, which we denote by E±{\cal E}_\pm, where k=±1k=\pm1 is the eigenvalue of KK. Write this as

E=E+⊕E−.{\cal E} = {\cal E}_+ \oplus {\cal E}_-.

This system and its symmetry become more interesting when time reversal is added to the mix of symmetries. See the continuation of this problem in Prob. Paragraph.