Part IV · Quantum Mechanics — Chapter 4.8
The Oscillator, and the Ladder
One potential, solved twice over. The differential equation was done in Chapter 4.5 and this chapter does not repeat it. What it does instead is take the same operator apart into two factors, discover that the factorisation fails by exactly one commutator, and get the whole spectrum out of that failure. The answer is famous. The method is what Parts V and VII are made of.
Chapter 4.7 solved four potentials and then said what it could not do. Its method was to write the general solution on each interval where the potential is constant and join the pieces where it changes, and that method needs pieces. The harmonic oscillator has none. Its potential is nowhere constant, there is no interval to solve on and no step to match across, so nothing from that chapter's toolkit touches it. Chapter 4.7's closing brick named the two ways forward from that position: solve one problem exactly by an algebra, or approximate all of them. This chapter is the first of the two, and Chapter 4.10 is the second.
The oscillator is worth the exact treatment for a reason Chapter 0.8 §4 already gave. Expand any smooth potential about a stable minimum and the constant is irrelevant, the linear term vanishes because the point was a minimum, and what survives is quadratic. So this is not one problem among many. It is what every stable system looks like close enough to equilibrium, which is why Chapter 0.8 called it the equation worth solving thoroughly, and why Chapters 5.3 and 7.4 are this single calculation with an index attached to it.
Here is the route. Section 1 says plainly that there are two ways to solve this problem, that one of them is already done and sitting in Chapter 4.5, and why the other one is what this chapter spends its pages on. Section 2 tries to factorise the energy into two pieces the way you would factorise a sum of two squares, watches the factorisation fail, and finds that the failure is exactly a commutator and exactly the ground-state energy. Section 3 turns that leftover into a machine that moves between eigenstates, and names the machine, because you will meet it three more times before the book is over. Section 4 is the step most books skip: showing that the machine has to stop, and that what stops it is the non-negativity of a length. Out of that comes , from the algebra and nothing else, once §5 has supplied the one state the algebra needs to exist. Section 5 finds the wavefunctions by solving a single first-order equation and climbing. Section 6 collects the phase-space area that Chapters 0.8 and 1.3 have been promising since Part 0, and checks the whole spectrum against a computation that shares no step with it. Section 7 builds a state that is not an energy eigenstate at all, moves exactly like a classical oscillator, and never changes its width. Sections 8 and 9 are worked examples and problems.
Conventions. The particle has mass and the oscillator angular frequency , and every length in the chapter is measured against , which Chapter 4.5 §5.1 introduced and which is the only length , and can be made into. Its partner is the only momentum, and . The letter always counts rungs and is always a non-negative integer. One result in this chapter is quoted rather than derived and it is experimental. It is the vibrational band positions of two isotopologues of hydrogen chloride, together with a bond length and three atomic masses, all used in Worked example 1. Its mark is raised at the head of §8, and nothing in §§1 to 7 is quoted at all. Worked example 1(c) additionally imports two premises that no chapter of this book derives: that the two isotopologues share a force constant, and that a reaction rate carries an Arrhenius factor. It names both where it uses them. Several marks standing in earlier chapters are leaned on and cited rather than raised again, and the closing brick lists all of it.
Tools you'll need — Chapter 4.6 §4.2 for the identification of with , §5 for in the position representation, §9 for stationary states, and §10 for a packet that spreads. Chapter 4.5 §5 for the Hermite functions, their generating function and their completeness, §5.5 for which operator actually is, and §8.2 for when an expansion in eigenstates is legitimate. Chapter 4.4 §3 for an operator as a formula together with a domain. Chapter 4.2 §8 for the canonical commutator, and §5.3 for the spread of an observable in a state. Chapter 4.7 §2 for parity and §3.5 for a zero-point energy read as an uncertainty. Chapter 1.3 §4.4 for and the Bohr–Sommerfeld guess, and §6.3 for the algebra a commutator obeys. Chapter 1.1's Problem 2 for the reduced mass. Chapter 0.9 §6 for the bandwidth theorem and §6.5 for the functions that saturate it. Chapter 0.8 §2.1 for a first-order linear equation solved by separation, §3.1 for the dimension of a solution space, §4.2 for the classical complex amplitude, and §4.4 for the phase-space ellipse and its area. Chapter 0.5 §1.1 for the inner-product axioms and §4 for the adjoint. Chapter 0.2 §4 for the Gaussian integral.
1 · Two routes, and why only one is taken
Here is where this section is going. There are exactly two standard ways to solve the quantum oscillator, one of them is already in this book, and this chapter takes the other. By the end of this section you will know what the route not taken looks like, why skipping it costs nothing, and what the route taken is for, which is not the answer but the technique.
1.1 · The operator, and what is already settled about it
Chapter 4.6 §4.2 identified the Hamiltonian of a particle in a potential as , and putting the harmonic potential into it gives the operator this chapter is about:
Chapter 4.4 spent a chapter on the fact that a formula like this is not yet an operator, because an operator is a formula together with a domain, and that for the particle in a box the missing information was four real numbers. Here that question is already answered, and the answer is not something this chapter has to work for. Chapter 4.5 §5.5 exhibited a complete orthonormal family of eigenfunctions for (4.8.1) and used it to write the operator as multiplication by a real function in a different description, which makes it self-adjoint by construction. That section said in so many words that it was settling the question in order to save an argument here. So we may treat (4.8.1) as a genuine observable and get straight to its spectrum.
What that means has a name, and two earlier chapters gave it. Chapter 0.5's closing list said the spectral theorem would return in this part of the book as the statement that solving a system means diagonalising its Hamiltonian, and Chapter 4.2's closing list picked out this chapter as one of the places it happens. Everything below is that diagonalisation: a complete orthonormal family of eigenvectors of (4.8.1), together with the number attached to each. Section 6.4 then does the same job a second time on a computer, and the two answers are compared.
1.2 · The route already taken
The first route is the differential equation. Write (4.8.1) in the position representation, strip out the constants with the length , and the eigenvalue problem becomes in the dimensionless variable , with the energy measured in half-quanta. Peel off the Gaussian that the large- behaviour demands, write what is left as a power series, and the equation turns into a recurrence between successive coefficients. The series then has to terminate, because if it does not the coefficients eventually behave like those of and the solution is not normalisable. Termination happens only at particular , and those are the eigenvalues.
That route exists, it is standard, and it is not being taken here. Two reasons, and the first is short. The series route does two separate jobs, and this book has both of them covered elsewhere. Exhibiting the solutions is Chapter 4.5 §5's job, and it is done. That section defined the Hermite polynomials by a generating function and extracted their recurrences. It then showed that the resulting functions satisfy the equation above and read the eigenvalues off, without a single power series. Excluding every other energy is the second job, the one termination performs in the series argument, and §4 below does it with a norm instead. The guarantee that no state has been missed is Chapter 4.5 §5.4's completeness proof, which §5.5 collects. Doing a Frobenius expansion here would produce the same number from the same side of the problem, and this book does not have a theory of special functions to spend it on.
The second reason is the one that matters. The algebraic route is not a cheaper way to the same answer. It is a different technique, and the technique is what the rest of the book runs on. Chapter 5.3 quantises a field by attaching a copy of this chapter's algebra to every mode, and the integer that labels a rung becomes the number of particles present. Chapter 7.4 does the same for a vibrating string, where the rungs are the particle spectrum. Neither of those uses a Hermite polynomial anywhere. What they use is and what follows from it, and that is what §§2 to 4 build. This book builds every special function it needs algebraically and never by a series, here and in Chapters 4.12 and 4.13 alike, and the decision was taken once for all three.
1.3 · What Chapter 4.7 hands over
Three things were promised to this chapter by name. Here they are, collected before starting rather than after finishing.
- The fifth exactly solvable potential. Chapter 4.7 solved four and said this is the one its method cannot reach, since the potential is nowhere constant and there are no pieces to join. What replaces matching is an algebra, and §2 is where the replacement happens.
- One of exactly two routes onward. Once the potential stops being piecewise constant you either solve one case exactly by algebra or approximate all cases. Chapter 4.10 does the second, and this chapter's exact answer is what its approximation gets scored against.
- The alternating parities, by a third route. Chapter 4.7 found that the states of a symmetric potential alternate even and odd up the ladder, and got it two ways: from a non-degeneracy argument in its §2.4, and from the roots of a matching condition in its §4.3. Section 5.6 below gets the same alternation for the oscillator out of the algebra alone, without looking at a single wavefunction.
What this chapter needs from Chapter 4.7 is only the shape of the answer: confinement gives a discrete ladder, the parities alternate, and the ground state does not sit at the bottom of the well. All three turn up again here, and none of them is assumed.
2 · Factorising , and what the leftover is
Here is where this section is going. We are going to try to write the energy as a product rather than a sum, because a product of two first-order objects is a far easier thing to handle than a second-order differential operator. The attempt does not quite work, and the interesting part is the amount by which it fails. That amount is the commutator of position and momentum, it is a pure number, and it is the ground-state energy of the oscillator arriving before any eigenvalue problem has been solved.
2.1 · One length, and the operator with the constants removed
Constants are in the way, and there is exactly one length that removes them. Out of , and you can build precisely one quantity with the dimensions of length, and one with the dimensions of momentum:
Chapter 4.5 §5.1 introduced the first of these for the same reason. Now define dimensionless operators by measuring against them, and . Two things happen at once. The Hamiltonian loses its constants, since and , and the canonical commutator of Chapter 4.2 §8 loses its , because dividing by leaves
Everything in this chapter comes out of those two statements. The second one is Chapter 4.2's postulate P6, the one relation about position and momentum that was assumed rather than derived, wearing different units. Nothing else about the operators will be used.
2.2 · The factorisation, and the piece that will not go
We want the energy as a product, and the shape of says what product to try. For two ordinary numbers, , and the reason that identity works is that the two cross terms and cancel. So write the same product with operators and watch the cancellation fail:
The cross terms do not cancel. What is left over in (4.8.4) is , and (4.8.3) makes that . Rearranging (4.8.4) to put the energy on the left carries that term across the equals sign and flips its sign, so the term standing beside the product below is , which is :
The sum of two squares has been written as a product, and the price of doing so is one extra term. That term is not an approximation, not a correction and not a small quantity. It is the commutator itself, and it is there for exactly the reason that position and momentum do not commute. Multiply (4.8.5) by and read what the price is in units of energy:
The is already here. It has arrived before any eigenvalue has been computed, before any wavefunction has been written down, and before any boundary condition has been imposed. It is the arithmetic cost of factorising a sum of two squares whose factors do not commute. Anything that later calls itself the zero-point energy of this oscillator is this term, and it can no more be removed than the commutator can. Two things are still owed before the term can be called the bottom of the spectrum rather than an offset above which the levels might start anywhere. One is that the product standing beside it in (4.8.6) is never negative, which is §4.2. The other is that a state the right-hand factor sends to zero exists at all, which is §5.2. Both are paid inside this chapter, and neither is assumed here.
2.3 · Naming the two factors
The two factors are going to be used constantly, so give them names, and choose the numerical factor in front so that the algebra afterwards is as clean as possible. Define
The second form is the first one with (4.8.2) substituted back, and it is the form to use when a real mass and a real frequency are in play. The factor is a convention and §2.4 is where it earns its keep.
Now the other factor. Chapter 0.5 §4 defined the adjoint by how it moves across the inner product, and one of the three rules it derived on the spot is that . Both and are self-adjoint, so taking the adjoint of leaves the two operators alone and conjugates the :
These two are the only new objects in the chapter's first half, and inverting the pair gives position and momentum back whenever they are wanted:
2.4 · The one relation everything else comes from
With the factor of in place, the commutator of the two new objects is as simple as it can be. Expanding using (4.8.7) and (4.8.8), the with and with terms commute and drop out, and what survives is
That is what the bought: without it the answer would be , and every square root later in the chapter would carry a stray factor. And rewriting (4.8.6) with the new names, since ,
What has happened is a genuine change of problem. Equation (4.8.1) was a second-order differential operator and its eigenvalue problem was a second-order differential equation. Equation (4.8.11) is a product of two first-order objects and a constant, and the only property of those objects that the next two sections use is (4.8.10). Not their form, not the space they act on, not what and mean. One line of algebra.
Two cautions attach to (4.8.7) and they will both be needed later, so they are worth stating at the point of definition rather than at the point of trouble.
is not self-adjoint, since its adjoint is and the two are different. So it is not an observable in the sense of Chapter 4.2's postulate P2, no measurement returns a value of it, and none of the theorems that chapter proved about Hermitian operators applies to it. In particular there is no reason for its eigenvalues to be real, and §7 will build states whose eigenvalue is a complex number. That is not a paradox and it is not an exception. It is what happens when an operator that was never Hermitian turns out to have eigenvectors.
What is an observable is the combination , whose adjoint is itself, which is why (4.8.11) is a legitimate Hamiltonian and on its own is not a legitimate anything.
Both are unbounded, which is Chapter 4.4's subject and not a technicality. Neither operator is defined on all of the space, and every manipulation below that moves across an inner product is using the adjoint relation of Chapter 0.5 §4 on vectors where both sides make sense. The states this chapter actually works with are the Hermite functions of §5 and finite combinations of them, and on those everything is finite, because each is a polynomial times a Gaussian and the ladder maps that set to itself. Section 7 then builds states that are infinite combinations, and §7.2 checks the condition for those separately rather than assuming it carries over. The restriction is stated out loud in §4.2 below, where it does real work, and §7.2 is the only other place it has to be met.
The energy of an oscillator is a sum of two squares, one for the motion and one for the stretch. For ordinary numbers a sum of two squares can be written as a product of two factors, and the reason is that the two cross terms in the product cancel each other. That trick is worth wanting here, because a product of simple things is far easier to work with than a sum of squares of complicated things.
Try it with position and momentum and the cancellation almost happens. The two cross terms do not quite kill each other, because the order in which you write position and momentum matters, and the difference between the two orders is the one thing quantum mechanics insisted on at the start. So the factorisation is correct apart from a single leftover term, and the leftover is precisely that difference.
Now look at the size of the leftover in units of energy. It is half of Planck's constant times the frequency, and it is exactly the amount by which the lowest state of an oscillator sits above the bottom of its well. That number has not been calculated yet. It has fallen out of an attempt to write a sum of two squares as a product, which is the closest this book comes to saying that the zero-point energy is a piece of arithmetic rather than a piece of physics.
3 · The ladder
Here is where this section is going. We have a Hamiltonian built out of two objects satisfying one commutation relation. What we want from them is motion between eigenstates: an operation that takes a state of definite energy and hands back another one, higher up. That operation exists, one commutator proves it in two lines, and the manoeuvre that does it runs three more times before this book ends.
3.1 · The number operator
The combination that appears in (4.8.11) deserves its own name, because every question about the energy is now a question about it. Define
Its adjoint is , using Chapter 0.5 §4's rule that taking the adjoint of a product reverses it, so is symmetric. And because is multiplied by a positive constant and shifted by another, the two have exactly the same eigenvectors, with eigenvalues related by . So the entire problem is now to find the eigenvalues of . We do not yet know that they are integers, and it would be circular to call them before showing it, so write for a generic eigenvalue and let the algebra decide what can be.
3.2 · Two commutators
To move between eigenstates of we need to know how reacts when is pushed past it, which is exactly what a commutator records. Compute it from (4.8.10) alone. Since commutes with itself, the only thing to move is the :
Neither line needs an identity quoted from anywhere. Writing the first one out, has a common factor of on the left, and what remains inside is , which is (4.8.10). The second goes the same way with the common factor on the right. Chapter 1.3 §6.3 noted that a commutator obeys the Leibniz rule exactly as a Poisson bracket does, and this is that rule used once. The two relations differ only in sign, so it is enough to follow the first and read the second off at the end.
3.3 · What a commutator like that does
A commutator is a difference of two orderings, and the use of one is that it lets you move an operator through another at the cost of a known remainder. Here the remainder is the operator itself, and that is what makes the relation powerful. Suppose is an eigenvector of with eigenvalue . Ask what does to the vector , and use (4.8.13) to move past so that it lands on the eigenvector, where we know what it does:
That line is the whole of the method. The vector is again an eigenvector of , and its eigenvalue is one greater. Running the same two lines with gives the mirror statement, that has eigenvalue . So the two operators move a state up and down a ladder of eigenvalues spaced exactly one apart, which in energy is a ladder spaced apart.
One qualification has to be attached now rather than later, because §4 turns on it. Equation (4.8.14) is an identity between vectors, and the zero vector satisfies it as readily as any other. So the honest reading is that is either an eigenvector with eigenvalue or the zero vector, which is not a state at all. The same disjunction holds downward. Everything interesting about the oscillator's spectrum comes from asking when the second alternative happens.
3.4 · The move, named once, because it happens three more times
Nothing in §3.3 used the oscillator. Strip the physics out and what remains is a pattern with three ingredients: an operator whose eigenvalues you want, a second operator , and a commutation relation for some number . Whenever those three hold, carries eigenvectors of to eigenvectors of with the eigenvalue shifted by , and the proof is (4.8.14) with the letters changed. Call it a commutator that shifts an eigenvalue. It runs three more times in this book, and each time the algebra around it is heavier and the shape harder to see, so here are the three in advance.
- Chapter 4.11 §4 runs it on angular momentum. There is , the shifting operators are , and the relation is . The extra work in that chapter is that the ladder has to close at both ends rather than one, and the argument that closes it is §4 below run twice.
- Chapter 4.13 §4 runs it on hydrogen. There the shifting operator carries you between different values of the orbital angular momentum at fixed energy, and it is built by factorising the radial Hamiltonian the way §2 factorised this one. It is the same manoeuvre, and that chapter says so.
- Chapter 7.3 builds the Virasoro algebra out of oscillator commutators, and Chapter 7.4 quantises the string with it, where the operator counting rungs is the one whose eigenvalue fixes a particle's mass. The algebra there is this chapter's with an index running over the string's modes.
So the move appears twice more inside the next five chapters and once more in Part VII. The reason to name it here, on the smallest example there is, is that here the whole content is (4.8.10) and a two-line calculation. By Chapter 4.11 there are three operators, a Casimir and a sign convention in the way, and a reader who has not been told what to look for will read the algebra instead of the shape.
4 · Why it stops
Here is where this section is going. Section 3 built a machine that lowers the eigenvalue by one every time it is applied, and nothing so far stops it running forever into negative numbers. Something must, and the something is not a physical argument, an appeal to a boundary condition, or a statement that the energy ought to be positive. It is that a length cannot be negative. By the end of this section the spectrum will be in hand, derived from (4.8.10) and one axiom of the inner product.
4.1 · The problem, stated sharply
Suppose with a genuine non-zero state. Apply once and, by §3.3, you have a state with eigenvalue or the zero vector. Apply it times and you have a state with eigenvalue or the zero vector. If the zero vector never turns up, then choosing larger than produces an eigenvector of with a negative eigenvalue, and choosing larger still produces one as negative as you please. Nothing in §3 forbids that. The spectrum would run down without end, the energies with it, and the oscillator would have no ground state.
So the question is not whether the ladder terminates. The question is what could possibly terminate it, and most books answer with a sentence to the effect that it must stop somewhere. It does not have to stop somewhere. It stops because of the following.
4.2 · A norm is not negative, and that is the whole argument
We want a lower bound on , so look for a quantity that is manifestly non-negative and that has in it. There is one, and it is the length of the vector . Chapter 0.5 §1.1 laid down the axioms of an inner product, and the last of them is that for every vector, with equality only for the zero vector. That is not a theorem about quantum mechanics. It is part of what the word length means. Apply it to and move one factor of across the inner product using Chapter 0.5 §4's defining property of the adjoint:
The middle step is the only one carrying a condition, and the condition is the domain question flagged in §2.4: it needs to be in the domain of and in the domain of . For the vectors this section is about, the eigenvectors of , both hold automatically. The operator is a product, and the domain of a product is by definition the set of vectors the two factors can be applied to in turn. That is what Chapter 4.4 §3 means by an operator being a formula together with a domain. Calling an eigenvector of already says that lies in the domain of , so it is one of those vectors, and nothing about which it belongs to has entered. That matters for §4.3, which applies (4.8.15) to a hypothetical eigenvector that is deliberately not one of §5's states. Section 5's states satisfy both conditions as well, being polynomials times a Gaussian, so (4.8.15) also says that the expectation value of is non-negative in every state this chapter writes down. Now take to be the eigenvector itself. The right-hand side is , the norm on the right is a positive number because is not the zero vector, and therefore
No eigenvalue of is negative. That is the bound the ladder needs, and notice what produced it. Not the form of the potential, not the behaviour of a wavefunction at infinity, not any claim about energies being positive. The positivity of the spectrum of is the positivity of , which is an axiom of the inner product and was written down in Chapter 0.5 before quantum mechanics appeared in this book at all.
4.3 · Closing the argument
Now put §4.1 and §4.2 together. Suppose is an eigenvalue and is not a non-negative integer. A negative is excluded by (4.8.16) at once, so the case to rule out is a positive that is not a whole number. Descending gives eigenvalues , and the descent can only be stopped by hitting the zero vector, which by (4.8.16) read backwards happens only when the eigenvalue being lowered from is zero. A number that is not a whole number never becomes zero by having whole numbers subtracted from it, so the descent is never stopped and eventually produces a negative eigenvalue. That contradicts (4.8.16). Hence every eigenvalue of is a non-negative integer, and we may finally call it .
The converse direction needs one more line, since so far we know only which values are permitted, not which occur. Suppose a state exists with . Then , so it has eigenvalue zero, and climbing from it produces eigenvalue at the -th step provided the climb never dies. It never does: by the same move as (4.8.15), with from (4.8.10),
So each rung is genuinely occupied, and the ladder is unbounded above and stops dead below. All of it rests on a state annihilated by existing at all, which is not yet proved. Section 5 proves it by solving one first-order differential equation, and finds exactly one solution.
4.4 · The spectrum
Substituting the permitted values of into (4.8.12) gives the result the chapter was built to reach:
Everything that went into it was: the canonical commutator, the definition of the adjoint, and the non-negativity of a norm. No series was summed and no boundary condition was imposed, and the only differential equation anywhere in the chapter is §5.2's, which is first order and settles nothing but whether the bottom rung exists. That is the sense in which this is a different technique and not a shorter road to the same one.
Three debts fall due here. Chapter 0.8's closing brick promised ladder operators, and the that will not go away, and it does not go away because §2.2 showed it is the commutator. Chapter 1.3 §4.4 imposed the Bohr–Sommerfeld condition on the classical phase-space area, got , and marked the condition itself as quoted forward to this chapter, saying the result would be derived here exactly, with ladder operators and no semiclassical approximation. It has been. Section 6 closes that loop from the other end by computing the area. And Chapter 4.5 §5.2, having produced the same eigenvalues from the Hermite side, said that this chapter would produce them a second way and get the same number.
The last of the three is the reason for doing a calculation twice, and it takes a paragraph. The two derivations share no step. Chapter 4.5 defined a family of polynomials by a generating function, differentiated it twice to get recurrences, assembled a differential equation and read the eigenvalue off the coefficient. This chapter wrote one commutator and used the fact that a length is not negative. They have the operator in common and nothing else, and they agree. An answer reached twice by disjoint routes carries more weight than an answer reached once by either, and more still when one of the routes generalises to field theory and the other does not.
4.5 · The rungs, normalised
Before going further we need to know not only that moves up a rung but by how much it stretches the vector, because §§5 to 7 all compute with these. Equation (4.8.17) already has the answer for the upward move, and running the same line downward with (4.8.15) gives when is normalised. So has length and has length . Each is a multiple of the neighbouring normalised state, and only one of the two multiples is free. Define to be , which fixes that phase to at every rung. The downward relation then follows rather than being chosen a second time, since , with the same sign:
The second relation at reads , which is the terminating condition of §4.3 written in the new notation. Iterating the first from the bottom expresses every state in terms of the lowest one, with the factorials coming from the product of the square roots:
One thing this does not yet establish is that each rung holds exactly one state. The algebra shows that a state at rung descends to a state at rung , so the number of independent states at rung is at most the number at rung , which is the dimension of the space of solutions of . That dimension is shown to be one in §5.2 below, and non-degeneracy follows.
Two operations have been built. One raises the energy of a state by a fixed step and the other lowers it by the same step, and the fixed step is Planck's constant times the frequency. That much is a matter of pushing one operator past another and collecting the remainder, and it involves no picture of the system at all.
The obvious worry is that lowering can be repeated forever, which would give states of arbitrarily negative energy and no lowest one. What rules that out is a fact about lengths. Apply the lowering operation to a state that sits on a definite rung, and the squared length of the result turns out to be exactly the rung number you started from. A squared length cannot be negative. So no rung is below zero, and the only way the descent can be stopped is by reaching a state the lowering operation sends to nothing at all.
That is the whole reason the energies of an oscillator are a discrete evenly spaced ladder starting half a step above the bottom of the well. The spacing comes from a commutator, the floor comes from a length being non-negative, and the half step comes from the failed factorisation of §2. It is worth noticing how little was used. Nothing here refers to a spring, a wavefunction, or the shape of the potential beyond the two lines it took to set the problem up.
a natural place to stop · the spectrum is derived and the algebra is finished; what follows turns the same two operators on the wavefunctions, the classical orbit, and a state that moves
5 · The wavefunctions, from
Here is where this section is going. We have a spectrum but no functions, and one thing still owed from §4.3: that a state killed by exists and is unique. Both are settled by the same short calculation, and the calculation is the payoff of the whole method. A second-order eigenvalue problem with infinitely many solutions to find has become one first-order equation, solved once, plus an algebra that generates the rest. By the end of the section the states will be identified with Chapter 4.5's Hermite functions, and their parities will have been read off without looking at any of them.
5.1 · One first-order equation
To turn into something solvable we need as a differential operator, which means going to the position representation. Chapter 4.6 §5 established that there acts as and as multiplication by . Substituting both into (4.8.7) and writing everything in the dimensionless turns the operator into
The relative sign is the one place to be careful, and it comes from the in (4.8.7) meeting the in . Checking (4.8.10) against these two expressions directly costs one line: applying and its reverse to a test function and subtracting leaves , and the halves in front make it .
5.2 · Solving it, and why the solution is unique
Now write out . The multiplies zero and can be dropped, leaving
This is a first-order linear equation with separable variables, which Chapter 0.8 §2.1 solved in general: divide by , integrate both sides, exponentiate. Doing that gives , and fixing the constant with Chapter 0.2 §4's Gaussian integral,
The is the whole difference between the two forms, and every wavefunction after this one is written in the second of them. The left-hand form is normalised against and therefore carries the dimensions of an inverse square root of length. The right-hand form pulls that dimension out into and leaves a dimensionless function of , normalised against . From here on the states are the dimensionless factors, and getting back to a function of means restoring the and reading as .
Two things follow at once. First, the ground state is a Gaussian, which is going to matter twice more in this chapter, at §6.1 and again at §7. Second, and this is what §4.5 was waiting for, the solution is unique up to a constant multiple. Chapter 0.8 §3.1 showed that the solution space of a homogeneous linear differential equation has dimension equal to its order, and (4.8.22) is of order one. So kills a one-dimensional space of functions, there is exactly one normalised ground state, and therefore exactly one state at every rung. The oscillator's spectrum is non-degenerate, and the reason is that a first-order equation has a one-parameter family of solutions.
Chapter 0.8 opened by promising that its central sentence, that a linear differential equation is an eigenvalue problem, would appear three times in different costumes, and named this chapter as one of them. Here is the costume. The eigenvalue problem that does the work is not the second-order one that the Hamiltonian looks like. It is (4.8.22), first order, with the eigenvalue zero, and an algebra supplies everything else.
The effort saved is easy to measure. The eigenvalue problem for (4.8.1) is a second-order equation, one for each of infinitely many eigenvalues, with normalisability deciding which of them are allowed. What replaced it was (4.8.22), solved once in three lines, plus (4.8.20).
5.3 · Climbing
With the ground state in hand, every other state follows from (4.8.20). In the position representation that reads
Doing the first two by hand shows the pattern. Applying the operator once, the derivative of the Gaussian brings down a factor , which adds to the already there: . Applying it again and dividing by gives . In every case a polynomial multiplies the same Gaussian, and the polynomial gains one degree at each step, because the operator either multiplies by or differentiates.
5.4 · These are the Hermite functions, and that is a derived fact
The polynomials just generated look like the Hermite polynomials of Chapter 4.5, and they are, but saying so needs an argument rather than an appeal to appearance. The argument is one line and it uses a result that chapter proved about polynomials, not about eigenvalues.
Chapter 4.5 §5.2 defined by a generating function and derived, in its grind box, the identity . Write and apply the raising operator of (4.8.21) to it. The derivative hits both factors, and the Gaussian contributes , which cancels one of the two terms:
The raising operator is the Hermite recurrence. That is the identification, and the normalisation follows from it. Chapter 4.5 §5.1 normalised its functions as , and dividing (4.8.25) by that constant and by the in (4.8.21) turns the coefficient into , which is precisely (4.8.19). Since agrees with (4.8.23), induction gives for every , including the sign and the constant.
What has and has not been imported needs stating precisely. The recurrence used above is a statement about a family of polynomials. It says nothing about energies. The eigenvalues of this chapter came out of (4.8.10) and (4.8.16) and would be what they are if Chapter 4.5 had never been written. What §5.4 adds is that the functions produced by the two routes are the same functions, which is a check the calculation could have failed and did not.
5.5 · Completeness, cited and not reproved
One question remains and this chapter does not answer it, because it has already been answered. Are these all the states? The ladder produces one state per rung and the algebra says nothing about whether some other vector of the space lies outside all of them. Chapter 4.5 §5.4 proved that no non-zero vector is orthogonal to every Hermite function, its §5.5 used that to complete the spectral verification for this exact operator, and its §5.6 stated how far the verification reaches. Chapter 4.3's closing brick promised in writing that this would happen, and said what would follow from it, which is that the present chapter expands in the Hermite functions freely.
So we do, from here on, and the licence is Chapter 4.5's rather than this chapter's. Any state of the oscillator may be written with finite, the sum converges in the norm, and Chapter 4.5 §8.2 lists what may then safely be done with it. Section 7 uses that licence hard, and it is the only place in the chapter that needs it. This also collects a promise made a long way back: Chapter 0.5 said that applying its orthogonalisation procedure with a weight on the whole line produces the Hermite polynomials, and that those are the quantum harmonic oscillator states of this chapter. They are, and (4.8.24) is how they got here without an orthogonalisation being run.
5.6 · The parities, by a third route
Chapter 4.7 §2 built parity, the operator that reflects a function through the origin, and showed it commutes with the Hamiltonian whenever the potential is even. The oscillator's is, so its states must each be even or odd, and Chapter 4.7 §3.3 said this chapter would find the alternation by a third route. Here it is, and the route uses no wavefunction.
Reflection reverses position and momentum, which in operator terms is and , using . Both appear linearly in (4.8.8), so the raising operator reverses too:
the second because (4.8.23) is a Gaussian in and is therefore unchanged by . Now apply to (4.8.20) and slide it through the copies of , inserting between each pair. Each passage costs a minus sign and the last lands on and does nothing, so
The parities alternate up the ladder starting even, exactly as they did in the infinite well and in the finite well. Three routes to the same fact now exist: Chapter 4.7 §2.4 got it from the non-degeneracy of bound states in one dimension, its §4.3 got it from which branch of a transcendental condition each root lies on, and this got it from moving one operator past another. Chapter 4.7 §2.5 listed the places parity gets spent, and Chapters 4.16 and 4.17 will spend it on states of this chapter's ladder in particular.
One corollary is free. The degree of is exactly , which Chapter 4.5 §5.4 established from the recurrence, so has at most zeros. Combined with (4.8.27), those zeros come in pairs at , apart from a possible one at the origin, which every odd state has. The numerically computed states in §6.5's figure carry exactly of them. The infinite well gave the same count, but the two chapters start their counters in different places, so the agreement needs stating in a form neither convention owns. In both problems, the state rungs above the lowest has exactly interior zeros. Chapter 4.7 numbers that state and this chapter numbers it .
Finding the states of an oscillator would ordinarily mean solving a second-order differential equation once for every energy level, with the requirement that the answer stay finite doing the work of picking out which energies are allowed. The method of this chapter replaces all of that with one first-order equation. Ask which function the lowering operation destroys, and you are asking for a function whose slope is minus its position times itself, which separates and integrates in three lines. The answer is a bell curve.
Every other state is then made by applying the raising operation to that bell curve as many times as you like. Each application multiplies by the coordinate and subtracts a derivative, so each one adds one power to a polynomial sitting in front of the same bell curve. Those polynomials are exactly the ones Chapter 4.5 built from a generating function, and the raising operation turns out to be their recurrence relation, which is a fact worth checking rather than assuming, since the two constructions had no reason to agree.
Two further things come out for nothing. Because a first-order equation has a one-parameter family of solutions, there is exactly one lowest state, and therefore exactly one state on every rung. And because reflecting through the origin reverses both position and momentum, it reverses the raising operation as well, so each rung is the opposite of the one below in its behaviour under reflection. Even, odd, even, odd, all the way up, established without drawing a single one of them.
6 · The phase-space area, collected
Here is where this section is going. The classical oscillator of Chapter 0.8 traced an ellipse in phase space, and the area of that ellipse was computed there and left with an explicit promise attached to it. This section pays the promise. First we measure how large a quantum state is in position and in momentum, which turns out to place it exactly on the classical ellipse of the same energy. Then we compute the enclosed area and find that it is , which is what Chapter 1.3 guessed semiclassically and what this chapter has now derived without guessing. The section ends with a numerical check of the entire spectrum that shares no step with any of the above.
6.1 · How big is a state
We want the spread of position and of momentum in the state , and (4.8.9) makes both easy, because it turns a question about into a question about which rungs the operators connect. Start with the means. Since lowers a rung and raises one, both send to a state orthogonal to , so and are both zero, at every rung. No energy eigenstate of an oscillator is anywhere in particular, and none of them is moving. Section 7 exists to fix that.
For the squares, expand from (4.8.9). Four products appear. The terms and move two rungs and contribute nothing to a diagonal expectation value, and the two surviving terms are and , whose sum in the state is :
The second comes out of the same expansion. Squaring the momentum in (4.8.9) puts an overall minus sign in front, from the , and also flips the sign of the two terms that survive on the diagonal, so the two changes cancel and the answer differs from the first only in carrying rather than . Two consequences are immediate. Since the means vanish, the spreads of Chapter 4.2 §5.3 are the square roots of these, and their product is
At the bottom rung that is , and that number has already been met. Chapter 0.9 §6's bandwidth theorem bounds the product of a function's width and its transform's width from below, and Chapter 4.6 §10.2 turned that into using the single substitution . So the ground state sits exactly on a bound proved before quantum mechanics entered the argument. Chapter 0.9 §6.5 went further and showed that the only functions attaining it are Gaussians, and (4.8.23) is a Gaussian, so the two statements agree and neither was arranged to. The general inequality for an arbitrary pair of observables belongs to Chapter 4.9 §2 and is not needed here. What this section has is one instance, computed exactly, before the theorem arrives, which is the order this book prefers wherever it can get it.
The other consequence is a check on the arithmetic that costs nothing. The mean kinetic energy is and the mean potential energy is , which is the same number. So the energy divides equally between the two, at every rung, and the sum is as it must be.
6.2 · The ellipse, rationed
Now bring in the classical picture, because the numbers in (4.8.28) are about to land on it exactly. Chapter 0.8 §4.4 plotted the classical oscillator's state as a point and found that energy conservation confines it to an ellipse, . Substitute the two quantum averages into the left-hand side, at energy . Each term contributes exactly one half, and
The root-mean-square point of the quantum state sits on the classical orbit of the same energy. That is not an approximation and it does not require to be large. It holds at , where the state has no classical counterpart at all.
Chapter 0.8 §4.4 then computed the area the ellipse encloses, getting , and Chapter 1.3 §4.4 wrote the same quantity as the loop integral and made it the central object of its phase-space section. Verified symbolically: integrating around the closed orbit returns exactly. Now put the derived spectrum (4.8.18) into it, and use :
Three promises are paid by that line. Chapter 0.8 §4.4 said that the classical ellipse is not replaced by quantum mechanics but rationed, and its closing brick called the area the area that Chapter 4.8 will quantise. Chapter 1.3 §4.4 marked the Bohr–Sommerfeld condition as quoted forward to here. And Chapter 1.3 said of its three observations that this chapter would confirm them by an exact operator calculation that uses none of this reasoning, which is a strong claim. Sections 2 to 4 used a commutator, an adjoint and the non-negativity of a norm. No orbit appeared, no adiabatic invariant, no correspondence argument and no appeal to large quantum numbers. The claim holds.
Equation (4.8.31) reads two ways. Successive allowed orbits enclose areas differing by exactly , so the ring between one orbit and the next has area equal to Planck's constant, always, everywhere on the ladder. And the innermost orbit encloses rather than nothing, which is the zero-point energy of §2.2 drawn as a picture: the state cannot be a point in phase space, and the smallest region it can occupy is half a unit of .
6.3 · What is this chapter's, and what is Chapter 4.10's
A boundary needs drawing here, so that you know exactly what has been established and what has not. Chapter 1.3's Bohr–Sommerfeld condition is a general statement about any bound one-dimensional motion, and it was a guess when it was made, supported by the argument that an adiabatic invariant is the right thing to quantise. This chapter has derived one case of it, exactly, and has said nothing about any other case.
The general statement is Chapter 4.10 §6, where it comes out of a genuine approximation scheme rather than a guess, and where it is scored: for potentials that are not quadratic the condition is approximate, the error is largest at the bottom of the ladder, and that chapter prints the numbers. So the oscillator is the case where the semiclassical answer is exact, and the reason it is exact is that the answer came from an algebra with no small parameter in it. Chapter 4.10 also supplies a count that Chapter 4.1 needed and could not get. That chapter's classical partition function had to divide phase space into cells before it could count states, and classical mechanics offered no size for a cell, so the sum was defined only up to a constant. The count that closes it is that an orbit enclosing area holds about states, which is (4.8.31) read as a density rather than as a level, and 4.10 is where it is proved.
6.4 · The numerical confirmation
Everything above came out of one commutation relation, so the whole of it can be held against a computation that has never heard of that relation. The check used here is a direct diagonalisation of (4.8.1) in position space, with no ladder operator, no Hermite polynomial and no generating function anywhere in it.
The construction is short to describe, and the one choice in it needs a reason. Work in units and put down a uniform grid of points at spacing , spanning . What has to be represented on that grid is the kinetic energy . An ordinary finite difference would represent it only to some power of the spacing, and the comparison below would then be measuring the difference scheme rather than the spectrum. So use instead the matrix that is exact on every function whose Fourier content fits between the grid's own limits. Its entries are on the diagonal and off it, all divided by twice the squared spacing. The is already inside them, so what those entries define is the kinetic-energy operator and not the second derivative, which is why the next step adds the potential rather than subtracting it. Add down the diagonal and diagonalise the resulting symmetric matrix. The whole computation is a discretisation of the differential operator and a rotation to its eigenbasis.
The lowest sixteen eigenvalues come out as , with the largest departure from anywhere in that range equal to and the largest departure of any gap from equal to . Both are at the level of double-precision rounding for a matrix of this size, so the levels are equally spaced as far as the arithmetic can see. Taking the same eigenvectors and computing , which is read straight off the grid, reproduces (4.8.28) to a worst relative departure of over the same sixteen states, and the corresponding reproduces (4.8.29) to twelve figures at every one of them.
Two remarks about what that does and does not test. The accuracy is not an accident of a fine grid: it comes from the eigenfunctions decaying faster than any exponential, so a grid that comfortably contains them and comfortably resolves them commits an error that falls off faster than any power of the spacing. The grid has to do both, and it fails at each end for its own reason. Shrinking the span to at the same points degrades the sixteenth level to , because the highest states no longer fit inside it. Stretching it to degrades the same level to , because the spacing is then too coarse to follow their oscillations. Between those two the window is wide, and inside it the answer is right to the level of double-precision rounding quoted above. The test is also genuinely independent, since a grid Laplacian shares no step with §§2 to 4. What it does not test is completeness, which is Chapter 4.5's and cannot be settled by any finite matrix.
6.5 · The ladder, operable
The figure below is the ladder with a button on it. Press and the state climbs a rung, press and it falls, and at the bottom pressing returns nothing at all, which is (4.8.19) at and is the whole reason the spectrum has a floor. The levels, the wavefunctions and the expectation values are all taken from §6.4's diagonalisation rather than from any formula in this chapter, so the readouts are testing the algebra rather than displaying it.
Chapter 0.8 drew the classical oscillator as a point going round and round an ellipse, with position on one axis and momentum on the other, and worked out the area the ellipse encloses. It came to the energy divided by the frequency, times two pi. That chapter then said, in writing, that a later one would ration those ellipses, and this is the later one.
Feed the energies just derived into that area and the answer is Planck's constant multiplied by a half-integer. So the allowed orbits are not a continuum. They are a set of nested rings, each enclosing exactly one more unit of Planck's constant than the one inside it, and the innermost enclosing half a unit rather than none. That last half unit is the same zero-point energy that appeared when the factorisation failed, now wearing its geometric costume: the state cannot shrink to a point, and the smallest patch of phase space it can occupy has a definite size.
Chapter 1.3 guessed this rule and said so honestly at the time. What has changed is that the rule is now a consequence rather than a guess, in this one case, and the derivation used nothing about orbits. The general version, and an honest account of how badly it does on potentials that are not parabolas, is Chapter 4.10's business.
a natural place to stop · the energy eigenstates are finished; what follows builds a state that is not one of them, and it is the state that behaves classically
7 · Coherent states
Here is where this section is going. Six sections have produced energy eigenstates, and §6.1 noted a strange fact about all of them: in every one, the mean position and the mean momentum are zero, and they stay zero forever, because a stationary state's expectation values do not move. Nothing built so far swings. This section builds the state that does, finds that its width never changes while it swings, and observes that it is not an energy eigenstate at all.
7.1 · What to ask for
The wanted state should behave as much like a classical oscillator as quantum mechanics permits, and the natural way to look for it is to ask what the classical variables become. Put Chapter 0.8's solution into the combination , using , and the two terms assemble into . Its modulus is the amplitude and its argument turns at the oscillator's own rate, which is that chapter's §4.2 complex amplitude set rotating and its §4.4 circle written as one number. That combination is precisely what (4.8.7) was built out of, so ask for a state on which acts as a number:
The warning in §2.4 is what makes that request legitimate. Because is not self-adjoint, nothing requires to be real, and none of Chapter 4.2's theorems about the eigenvalues of an observable applies. There is no measurement whose outcome is . What there is, if such states exist, is a family labelled by a complex number, and the label should turn out to be the classical amplitude.
7.2 · Building it
To find the states, expand in the basis we have. Chapter 4.5 §5.4's completeness is what permits the expansion, and §5.5 above collected it, so write and apply (4.8.19) term by term. Lowering shifts every index down one, so
Now compare the coefficient of on the two sides, which is legitimate because the are orthonormal. On the left the term with comes from , so . That recurrence unwinds immediately, giving , and the constant is fixed by demanding , where the sum is the exponential series. So the state exists, for every complex , and it is
Those coefficients also settle the domain question of §2.4 for this state, which is not a finite combination and so is not covered by what that box said. Since falls off faster than any power of , the sum converges for every , so lies in the domain of every polynomial in and . That is what licenses moving across the inner product in §7.3 and §7.5, and again in Problem 3(a).
7.3 · Its position, its momentum, and its width
What we want to know first is where this state is and how wide it is, and both come from (4.8.9) with (4.8.32) used on the right and its adjoint on the left. Since and , every expectation value reduces to arithmetic in :
So the label is the position and momentum of the state, in the natural units, and the state can be put anywhere. For the widths, use the same trick on , being careful with the one term that is not already in normal order: write as before evaluating, which is (4.8.10) earning its keep for the last time. The cross terms combine into exactly, and the survives on its own:
What is missing from those expressions is the point: the widths do not depend on . Every state of this family, however far from the origin it is placed and however fast it is moving, has exactly the width of the ground state, and saturates Chapter 0.9 §6's bound. That is consistent with §6.1, where the ground state saturated it too, and indeed in (4.8.34) gives back. The coherent states are the ground state, moved.
7.4 · The state that does not spread
Now let it evolve, which is where the point of the whole section is. Chapter 4.6 §9 established that a stationary state acquires the phase and nothing else, so each term of (4.8.34) picks up . Pull the half out front, where it is a global phase and affects nothing measurable, and what is left multiplies by , which is the same as replacing by throughout:
A coherent state evolves into a coherent state, with its label rotating at the classical frequency. Everything follows from that one line. Since never changes, the widths in (4.8.36) never change either, and the packet does not spread. And feeding the rotating label into (4.8.35), with ,
which is Chapter 0.8's classical solution exactly, amplitude and phase and all. The state is a rigid Gaussian of fixed width sliding back and forth along the classical trajectory.
This settles a debt from the chapter before last. Chapter 4.6 §10 built a free Gaussian packet and watched it spread, and its closing brick said that this chapter would do the same for the oscillator and find the state whose width does not change. That state is (4.8.34), its width is (4.8.36), and in the units that chapter's figure uses, where , the number is . That figure's oscillator panel was run at exactly that width for exactly this reason.
Equation (4.8.38) is easy to over-read, and two boundaries keep it in place. Chapter 4.9 §5 will show that the mean position of any state follows the classical trajectory exactly whenever the potential is at most quadratic, so the motion of the centre is not special to coherent states. What is special is the width. And no state of any system is genuinely classical: (4.8.36) is a fixed non-zero spread, and §7.5 says what else is spread.
7.5 · Not an energy eigenstate, and that is the point
After six sections of eigenstates, (4.8.34) is not one of them, and that needs saying plainly. It is a superposition of all of them, with weights that never vanish, and its energy is not definite. The mean rung follows from , and the spread from the same normal-ordering step used in §7.3:
The probability of finding rung is , which is a Poisson distribution with mean , and (4.8.39) is its standard deviation written out. That is the reason a coherent state looks classical when it is large. The fractional uncertainty in the energy falls like , so a state carrying a million quanta has its energy defined to a part in a thousand, and a state carrying one does not have a defined energy in any useful sense.
7.6 · Where this goes
Chapter 0.8 §7 closed with a sentence this chapter can now honour. A field in a box is a collection of independent modes, each one a harmonic oscillator with its own frequency, so quantising a field means doing the calculation of §§2 to 4 once and attaching a copy of it to every mode. The integer labelling a rung then counts how many quanta of that mode holds, and those quanta are what we call particles. A photon is one rung of one mode of the electromagnetic field, and Chapter 5.3 is where that sentence becomes a construction. Chapter 7.4 runs the same machine on the vibrational modes of a string, where the rungs are the particle spectrum. Chapter 0.3's closing list named all three destinations in order, and this chapter is the first of them.
Coherent states are the piece that travels furthest in one respect. A field mode in a coherent state has a definite amplitude and phase and an indefinite number of quanta, which is what a classical electromagnetic wave is, and it is why the classical limit of a quantum field is a coherent state rather than a state of many photons with definite number. The half rung per mode travels too, and less comfortably: every mode of every field carries whether anything is in it or not, and adding that up over infinitely many modes gives infinity. Chapter 5.3 has to say what to do about it before it can say anything else.
This box picks up the second paragraph of "Where we are" rather than anything in §§1 to 7. That paragraph is where the chapter said why one potential deserves an exact treatment, and the reason it gave was Chapter 0.8 §4's. Expand a smooth potential about a stable minimum: the constant is irrelevant, the linear term vanishes because the point was a minimum, and what survives is quadratic. That expansion, and not the algebra of this chapter, is the move you make every time you fit a model. A log-likelihood is maximised at , the first derivative vanishes there because it is a maximum, and expanding to second order leaves with the observed information. Everything you then do with a Wald interval treats that parabola as the whole of the likelihood near the estimate, and the width that comes out, , is built from the curvature and from nothing else. Chapter 0.8's expansion is the identical algebra with the sign flipped, and plays the part of .
The box sits here rather than beside that paragraph because both of the places the analogy breaks need equipment assembled later: the ground state's actual width, which is §6.1's, and the discrete ladder, which is §4's. It waits until the end of the main line so that both are in hand. Where the analogy stops is more interesting than where it holds, and there are two places.
The width needs a second constant. In the likelihood case the curvature alone fixes the standard error, because the units of are already the units of the answer. Here it does not: has the dimensions of a stiffness and no length can be made from it. The length of (4.8.2) needs as well, and the ground state's width is that length. Quantum mechanics is the extra constant, and every quantitative statement in this chapter is a statement about it. Set to zero and the state shrinks to a point sitting at the bottom of the well, which is the classical answer.
There is no ladder in the statistical case, and there cannot be. A likelihood surface has a continuum of parameter values near its maximum and no notion of an allowed set. What produces the discrete rungs here is not the parabola but the commutator, which has no counterpart in a log-likelihood at all, and which is what §3 turned into an integer. So the shared structure is the expansion and only the expansion. Read further than that and you would be predicting a discrete set of permitted hazard ratios, which is exactly the kind of nearly-right analogy that costs more than it buys.
Every state built so far has a definite energy, and every one of them has a curious property: on average the particle is at the centre and not moving, and it stays that way forever. Nothing so far in the chapter oscillates, which is odd for a chapter about an oscillator.
The state that does oscillate is found by asking for something different. Instead of a state on which the energy takes a definite value, ask for one on which the lowering operation acts as multiplication by a number. That is a legitimate thing to ask for, because the lowering operation is not an observable, so nothing forces the number to be real, and such a state does exist for every complex number you name. Expanding it over the rungs gives coefficients that make its probabilities a Poisson distribution, of the kind that describes counts of independent rare events.
What this state does is the point. Its centre moves exactly along the classical path, with the right amplitude and the right phase, and its width never changes at all, which is the opposite of what a free packet does. It is a bell curve of fixed shape sliding back and forth in the well. The price is that it has no definite energy: its rung is uncertain by the square root of its mean rung, so a large one is sharply defined in relative terms and a small one is not. That is the sense in which a classical oscillation is a quantum state with a great many quanta in it and no idea how many.
8 · Worked examples
Everything in this chapter up to here is a consequence of (4.8.10) and the axioms of an inner product. Worked example 1 tests it against a molecule, and to do that it needs numbers that no calculation in this book produces. They are quoted, which is what this mark means.
The vibrational constants of hydrogen chloride. Infrared spectroscopy of the two isotopologues gives harmonic vibrational wavenumbers of for and for . The observed fundamental band of the first, meaning the transition from the lowest rung to the next, is at , and its first overtone, from the lowest rung to the one above next, is at . These are measurements.
Two more quoted numbers, used once each. The equilibrium bond length of is , which part (b) uses only to say what fraction of the bond the ground state's motion amounts to. And the three atomic masses, for , for and for , are measured quantities like any other.
What is done with them. The example computes a force constant from each isotopologue separately and asks whether the two agree. It also predicts the second molecule's frequency from the first's, using the masses together with one premise about the electronic problem, which part (c) states where it uses it. Neither of those uses the measurements to fix a parameter of the theory, because the theory of this chapter has no free parameter to fix. No result in §§1 to 7 depends on any of this, and removing this box would change nothing above it.
Two premises, imported and not derived. Part (c) also leans on two statements that no chapter of this book proves, and both are premises the example brings in from outside rather than results it inherits. The first is that two isotopologues of the same molecule share a force constant, which is a claim about the electronic energy curve. Part (c) treats that one as a hypothesis under test rather than as a fact, and reads the residual as its failure rather than the oscillator's. The second is that a reaction rate carries an Arrhenius factor , which is a claim about chemistry, and that one is assumed rather than tested: the kinetic-isotope factor of rests on it entirely. Both are named again where part (c) uses them. Nothing in §§1 to 7, and nothing in the other worked examples or problems, depends on either.
Hydrogen chloride is two atoms joined by a bond, and near its equilibrium separation the bond behaves like a spring. (a) Treat the vibration as a one-dimensional oscillator and say what plays the part of the mass. (b) From 's wavenumber find the force constant, the zero-point energy and the amplitude of the ground state. (c) Predict 's wavenumber from 's alone, and test it. (d) Say where the harmonic ladder is visibly wrong.
(a) Two masses joined by a spring is not literally the problem of this chapter, since there are two coordinates rather than one. Chapter 1.1's Problem 2 changed variables to the centre of mass and the separation, and found that the kinetic energy splits into one term for each with no cross term at all, the second carrying the reduced mass . That change of variables is linear, so it does the same thing to the momentum operators, and the internal motion is (4.8.1) with replaced by . For the atomic masses and give , which is . The reduced mass is almost the hydrogen mass, which is why hydrogen stretches sit at the high-frequency end of an infrared spectrum.
(b) Spectroscopy quotes a wavenumber , so first convert: , which for gives . Then the force constant is , which comes to . That is a stiff spring by everyday standards and a completely ordinary one for a chemical bond.
The zero-point energy is , by (4.8.18) at , and it is , or , or half the quoted wavenumber, . Compare that with at room temperature, which is . The molecule's vibration carries seven times the thermal energy scale in a state from which no energy can be removed, and the fraction of molecules in the first excited rung at is . Essentially every hydrogen chloride molecule in a room is on the bottom rung, and is nonetheless vibrating.
The amplitude follows from (4.8.28). The oscillator length is , so the root-mean-square displacement of the ground state is , against a bond length of . The bond breathes by about six percent of its length while in its lowest state, which is small enough for the quadratic expansion to be reasonable and large enough that part (d) can detect its failure.
(c) This test uses no fitted quantity at all. The force constant is a property of the electronic energy curve, which depends on the charges and not on the nuclear masses, so the two isotopologues should share it. That is a premise imported from the electronic problem, and neither this chapter nor any other in this book derives it. The flag box above says so, and the test below is as much a test of that premise as of the ladder. What this chapter supplies is the consequence: if is shared then , and the ratio of the two frequencies is fixed by the masses alone,
With from the deuterium mass , the predicted ratio is . The measured ratio is . Predicting 's wavenumber from 's alone gives against measured, so the wavenumbers depart by two parts in a hundred thousand. Read the other way, the force constants extracted separately from the two molecules are and , and those agree to four parts in a hundred thousand. The two residuals are the same discrepancy counted twice, since goes as and the force constants therefore disagree by twice as much as the wavenumbers.
What that residual is, is the failure of the shared-force-constant hypothesis rather than of the oscillator, since nothing in (4.8.18) knows about isotopes. The zero-point energies differ by between the two molecules, and that difference is a real and measurable thing. Suppose a reaction rate goes as in the energy needed to reach the top of a barrier. That is the second imported premise, a statement about chemistry that this book derives nowhere, and unlike the shared force constant it is assumed rather than tested. Granted it, the heavier isotope starts lower and has that much further to climb, which at room temperature is a factor of in the rate. That factor is the primary kinetic isotope effect. Its size comes out of (4.8.18) and the two masses once the Arrhenius form is granted, and it does not survive without it.
(d) The ladder of (4.8.18) is exactly evenly spaced, so the transition from rung to rung must cost exactly twice the transition from to . Measured, the ratio is , low by . The fundamental itself sits below the harmonic wavenumber quoted in the flag box, for the same reason. A real bond is not a parabola: it softens as it stretches and it dissociates, so the true levels crowd together going up rather than staying evenly spaced. Two percent is the size of that effect at the second rung of a stiff bond, and it is what Chapter 4.15 computes by treating the departure from a parabola as a perturbation. Problem 4 below sets up the first step of that calculation with this chapter's algebra.
(a) Compute for all and . (b) Do the same for . (c) Say what the answers mean for which transitions a light wave can drive. (d) Check both against direct integration of the Hermite functions.
(a) The point of the exercise is that no integral is needed. Chapter 4.5 would have you compute ; the ladder computes it by asking which rung the operator moves you to. Put (4.8.9) between the two states and act with (4.8.19) on the right:
So connects a rung only to its two neighbours, and to nothing else, ever. The whole infinite matrix has two diagonals and is zero everywhere else, and the entries grow like .
(b) Square the same expression, keeping the operators in order. Four products appear. takes to with coefficient , its adjoint goes the other way with , and the two mixed terms both stay put and add to :
Three diagonals, and the middle one is (4.8.28) recovered. The pattern behind it is that changes by one, by nought or two, and by an amount of the same parity as . That is (4.8.27) restated, since an integral of an odd function over a symmetric interval vanishes, and it is why parity kills matrix elements without their being computed.
(c) A light wave couples to a charge through a term proportional to , which Chapter 4.17 derives. Part (a) then says the coupling vanishes unless the rung changes by exactly one, which is the selection rule . So an oscillator absorbs and emits at one frequency, , no matter which rung it starts on, and that is why a harmonic vibrational spectrum is a single line rather than a series. Worked example 1(d) is the measurement of how badly that fails, and it fails by two percent, which is the size of the departure from a parabola.
(d) Numerical integration of on the line, using the Hermite functions of Chapter 4.5 §5.1 directly, returns , and for , and , against part (a)'s at , and , with the entries at and vanishing to eleven decimal places. For the diagonal entries return , and at , and the two-step entries at and at , all matching the formulae to ten figures. The integrals were never needed and they agree.
A charge in an oscillator sits in a uniform field , adding to the potential. (a) Solve it exactly with the algebra. (b) Identify the new ground state. (c) Put numbers on it for a trapped ion. (d) Say what Chapter 4.15 will do with the answer.
(a) The new Hamiltonian is , and the fastest route is to notice that the added term is linear in , hence linear in and by (4.8.9). A term linear in and sitting beside can be absorbed by shifting both, which is completing a square with operators in place of numbers. Define with a real constant. Subtracting a number changes no commutator, so , and every consequence of §§3 to 5 holds for with no work at all. Expanding and comparing with fixes , and what is left over is a constant:
The spectrum is the old one pushed down bodily. The spacing does not change, so the spectroscopy does not change either, and the shift is quadratic in the field rather than linear even though the added term was linear. In position space the same statement is the completion of the square with , verified symbolically, which is the classical displacement of the equilibrium point.
(b) The new ground state is the state killed by , which means . By (4.8.32) that is a coherent state of the original ladder, with label . So §7's states are not an ornament: they are the ground states of displaced oscillators, and a field applied to an oscillator in its ground state leaves it in a coherent state of the original one. Writing the shift in terms of the label, , so the depression of the ladder is one quantum for every unit of , which is also the mean rung of the displaced state measured on the old ladder by (4.8.39).
(c) Take a ion in a trap of frequency , a standard laboratory setting. Its oscillator length is and its level spacing is . A stray field of one volt per metre, which is very small, displaces the equilibrium by , which is oscillator lengths, so and the ladder drops by quanta. The same field applied to an electron bound at moves nothing measurable, and it takes to shift that ladder by one percent of its spacing. The difference is entirely in , which is for the electron and for the ion, twelve orders of magnitude apart. This is why trapped-ion experiments compensate stray fields as their first act of calibration.
(d) Chapter 4.15 develops perturbation theory, in which a small addition to a Hamiltonian is handled by an expansion. Applied to , its first-order term is , which Worked example 2(a) says is zero, and its second-order term is a sum over the two neighbouring rungs. This example says in advance what that sum has to produce, namely at every rung and nothing at any higher order. An exact answer waiting for an approximate method is the most useful thing a chapter can hand forward, and it is why this example is here rather than there.
9 · Your turn
Problem 1 — the oscillator in three dimensions, and half a quantum per direction
(a) An isotropic three-dimensional oscillator has . Define one ladder per direction and give the spectrum. (b) Show that the level is degenerate and count the states at each . (c) Say what the zero-point energy is and what it is per degree of freedom. (d) Chapter 4.11 will find that a rotationally symmetric problem has levels labelled by an angular momentum with states each. Check that your count in (b) is consistent with that, and say which occur at each .
Solution
(a) The three directions do not talk to each other: by Chapter 4.2 §8, so operators built from different directions commute. Define by (4.8.7) in each direction. Then and . Every argument of §§3 and 4 applies to each separately, so each has spectrum and . Writing for the total, .
(b) The count is the number of ways to write as an ordered sum of three non-negative integers. Fix , and then ranges over to , which is choices with forced. Summing over from to gives . Checked by direct enumeration: for . Unlike the one-dimensional case, the levels are degenerate, and the degeneracy grows quadratically.
(c) The ground-state energy is , which is for each of the three directions. That is the general rule, and §2.2 is why: each independent quadratic pair contributes one failed factorisation, hence one commutator, hence one half quantum. A system with quadratic degrees of freedom has zero-point energy if they share a frequency, and if they do not. Chapter 5.3 takes to infinity and has to deal with the consequences.
(d) The degeneracies must decompose into odd numbers . They do: , , , , . So level contains down to or , and only of the same parity as appears, which is (4.8.27) again in three dimensions. Levels with different share an energy, which is a degeneracy rotation alone does not explain. Chapter 4.14 meets the same phenomenon in hydrogen and identifies the larger symmetry responsible.
Problem 2 — the argument of §4, run in both directions
(a) Suppose had an eigenvalue . Follow the descent explicitly and say at which step the contradiction with (4.8.16) appears. (b) Where exactly does the same argument fail to bound the ladder from above? Identify the step that has no upward counterpart. (c) Chapter 4.11 §4 runs this machinery on angular momentum and finds a ladder bounded at both ends. What must be true of that problem which is not true here? (d) Suppose someone proposes an operator with , together with a state that annihilates. Compute and say what it rules out. Then drop the second hypothesis and say why the sign on its own settles nothing.
Solution
(a) Descending gives eigenvalues , , . None of the intermediate states is the zero vector, because (4.8.15) gives , which is strictly positive at , and . So the state at genuinely exists, and (4.8.16) applied to it says . The contradiction arrives at the third application, and it would arrive eventually for any that is not a non-negative integer.
(b) The bound came from . The upward counterpart is (4.8.17), , and non-negativity of that gives , which is weaker than what is already known and constrains nothing above. The asymmetry is in the commutator: but , so the downward norm can vanish and the upward one cannot. Physically, the Hamiltonian is bounded below and not above, and the algebra knows it.
(c) There the operator whose eigenvalues are being bounded is , and the quantity playing the part of is , which is bounded above and below once is fixed, because it is a fixed number minus a square. Here has no such ceiling. The general lesson is that the ladder terminates wherever a non-negative quantity would be forced negative, and how many ends it has is decided by how many such quantities there are.
(d) The relation gives , so with ,
A vector of squared length cannot exist in a space where the inner product is positive definite, which is Chapter 0.5 §1.1's last axiom. So the two hypotheses together are impossible, and one of them has to go.
Now drop the second. Without a state that annihilates, the sign alone rules out nothing, because satisfies and the algebra is this chapter's with the two names exchanged. What the sign does is decide which of the pair is the lowering operator, and the trouble in the display above comes from the physics insisting on the other one. That is exactly the situation Chapter 7.4 meets when it quantises the timelike component of a string coordinate: the metric puts a minus sign in the commutator while the definition of the vacuum is fixed independently, so states of negative norm appear, and removing them honestly is a substantial part of that chapter.
Problem 3 — counting statistics in a coherent state
(a) Verify from (4.8.34) that the probability of rung is Poisson with mean , and derive without summing a series. (b) A laser beam of power at wavelength is a coherent state of one field mode. Estimate for one nanosecond of a one-milliwatt beam at , and state the fractional uncertainty in the number of photons. (c) Two coherent states with different labels are not orthogonal. Compute and say when it is negligible. (d) Why does (c) not contradict Chapter 4.2's theorem that eigenvectors of an observable belonging to distinct eigenvalues are orthogonal?
Solution
(a) Squaring the coefficient in (4.8.34) gives , which is the Poisson probability with mean . For the spread without a series, use the operators. . For the second moment, normal-order first: , whose expectation is . Subtracting the square of the mean leaves , so . Variance equal to mean is the Poisson signature, and it arrives here from one application of (4.8.10).
(b) A photon at carries , and a milliwatt for a nanosecond delivers , so . The fractional uncertainty is . The beam has a well-defined amplitude and phase and an indefinite photon number, and the indefiniteness is a twentieth of a percent, which is why nothing about a laser looks quantised. It is also measurable, and it is the shot noise floor of any photodetector.
(c) Using (4.8.34) twice and the orthonormality of the , , so . Two coherent states overlap unless their labels are separated by more than about one, which in position is more than about one oscillator length. They are a family that is complete and not orthogonal, which is unusual and useful.
(d) That theorem is about eigenvectors of a self-adjoint operator, and its proof uses self-adjointness at the one step where the eigenvalue is moved across the inner product. The warning in §2.4 was that is not self-adjoint. The mirror statement is that is false; the correct relation is , and the proof of orthogonality does not run.
Problem 4 — the first anharmonic correction
A real bond is not a parabola, and Worked example 1(d) measured the departure. Model it by adding to (4.8.1). The letter is being reused: this is the size of the quartic term and has nothing to do with the displacement label of Worked example 3. (a) Compute with the ladder. (b) Chapter 4.15 will show that the leading shift of level is . Using that, say how the spacing between neighbouring levels changes with , and compare the sign with what Worked example 1(d) measured. (c) Why does an term contribute nothing at this order, and what does that have to do with §5.6? (d) Estimate for hydrogen chloride from Worked example 1(d)'s measurement, and say whether the perturbation is small.
Solution
(a) Write and use Worked example 2(b): the state has components on rungs , and with coefficients , and . Since is self-adjoint, is the squared length of that vector, which is the sum of the three squares:
Checked against the values , , , and in units of at , computed from the matrix directly. No integral of a fourth power against a Gaussian was needed anywhere.
(b) The shift is , so the gap between and changes by , growing linearly in . For the levels spread out going up, and for they crowd together. The measurement in Worked example 1(d) puts the gap below the gap, so the levels crowd, so the effective is negative. That is what a bond that softens as it stretches should give, and it is the correct sign.
(c) An term connects rungs differing by or , by the same counting as Worked example 2(b), and never connects a rung to itself. So its diagonal element vanishes and it contributes nothing at leading order. The reason is (4.8.27): an odd power of is odd under reflection, the states have definite parity, and the integrand is odd. The cubic term does contribute at the next order, where it enters squared, and Chapter 4.15 handles it there.
(d) Part (b) says the overtone falls short of twice the fundamental by . Measured, it falls short by times the fundamental, which is about of , so . The quartic term is therefore about one percent of the level spacing where the ladder is being used, which is what makes perturbation theory the right tool and Chapter 4.15 the right chapter. At high rungs the correction grows like and the expansion fails, which is the same statement as a real bond having only finitely many bound vibrational levels before it dissociates.
The factorisation, and where the half came from. The energy of an oscillator is a sum of two squares, and a sum of two squares wants to be a product. Writing as fails by exactly one term, that term is , and in units of energy it is . So the zero-point energy is visible in (4.8.5), before any eigenvalue problem has been posed, and it is the arithmetic price of factorising something whose factors do not commute. Naming the two factors and turns the Hamiltonian into and the canonical commutator into , and from that point onward nothing about position, momentum, or the space they act on is used again.
The ladder, and the norm that stops it. From follow and , and a commutator of that shape shifts an eigenvalue by a fixed step. That is the whole mechanism, and it is why the levels are evenly spaced. What bounds the ladder below is not a physical argument but , together with the fact that a squared length cannot be negative, which is the positive-definiteness axiom of Chapter 0.5 §1.1 and nothing else. A non-integer eigenvalue would descend past zero and produce a state of negative norm, so the eigenvalues are the non-negative integers and . Chapter 4.5 §5.2 reached the same numbers from a generating function and a differential equation, sharing not one step with this, and said in place that it expected them to be reproduced. They were.
The states, from one first-order equation. The ground state is what annihilates, and in position space that is , solved by separation in three lines and giving a Gaussian. Because a first-order equation has a one-dimensional solution space, the ground state is unique and every level is non-degenerate. Climbing with generates the rest, and the raising operator turns out to be the Hermite recurrence , so the states produced here are Chapter 4.5's functions including their normalisation, which is a check that could have failed. Completeness is Chapter 4.5 §5.4's and is cited rather than reproved, which is what that section was built to license. Parity comes out of the algebra alone, since reflection reverses and fixes the Gaussian, giving and the alternation Chapter 4.7 predicted by two other routes.
The classical connection, and the numbers. The ladder gives and , so the ground state saturates Chapter 0.9 §6's bound and does so as a Gaussian, which is the only thing that can. Those two averages place the state exactly on Chapter 0.8's classical ellipse of the same energy, at every rung including the lowest. The enclosed area is , verified symbolically, so successive orbits differ by exactly one unit of Planck's constant and the innermost encloses half a unit. That collects Chapter 0.8's the area that Chapter 4.8 will quantise in full, and Chapter 1.3 §4.4's Bohr–Sommerfeld mark in one case only. Section 6.3 draws that second line: 1.3's condition is a general statement about bound one-dimensional motion, the general version is Chapter 4.10 §6's, and this chapter has derived the quadratic case exactly and said nothing about any other. It also collects Chapter 1.3's claim that an exact operator calculation using none of that reasoning would confirm it. Verified: a -point grid Laplacian with on the diagonal, containing no ladder operator and no Hermite function, returns the sixteen lowest levels as with the worst departure and the worst gap error , and its eigenvectors reproduce to .
One mark, four leaned on, and two premises carrying none. The single flag in this chapter is experimental and is raised at the head of §8: the vibrational band positions of two isotopologues of hydrogen chloride, together with a bond length and three atomic masses, used in Worked example 1 to test a mass-ratio prediction that has no free parameter. Nothing in §§1 to 7 depends on any of it. Every result in §§1 to 7 is derived here or in a chapter that derived it. Worked example 1(c) is the one place where that stops being true, and it says so twice in its own text. It imports two premises: that two isotopologues share a force constant, which is a statement about the electronic problem, and that a reaction rate carries the Arrhenius factor , which is a statement about chemistry. Neither is derived in this book and neither will be. The first is under test rather than assumed, and the second carries the factor of entirely. Four marks standing elsewhere are leaned on and cited rather than raised again. Chapter 4.2 §8's postulate P6 is the commutator everything is built from. Chapter 4.6 §4.2's identification of with , which that chapter marks as chosen rather than derived, is the operator being solved. Chapter 4.5 §3's spectral theorem stands behind the claim that the ladder exhausts the spectrum, reaching this chapter through its §5.4's completeness proof. And Chapter 0.8 §1.1's Picard–Lindelöf mark stands behind §5.2's non-degeneracy, since the dimension count cited there is Chapter 0.8 §3.1's, and that is the section where 0.8 says its quoted theorem gets spent. One further debt is not a mark: Chapter 4.5 §5.5 settles which self-adjoint operator the formula (4.8.1) denotes, a question Chapter 4.4 showed is never automatic, and it settles it by construction rather than by quotation. Two statements are borrowed from chapters that come after this one, and both sit outside the derivations: Worked example 2(c) uses the fact that light couples to a charge through a term proportional to , which Chapter 4.17 derives, and Problem 4 is set up on Chapter 4.15's first-order perturbation formula, which it states as the problem's premise. Nothing in §§1 to 7 rests on either.
Where this gets spent. Chapter 4.9 §1 turns §6.1's one instance into the bound for every state, in a single multiplication, and its §2 then proves the general relation for an arbitrary pair of observables. Its §5 explains why the mean of any state in a quadratic potential follows the classical path, of which §7.4's coherent state is the case where the width holds still as well. Chapter 4.10 §6 recovers the Bohr–Sommerfeld condition from a real approximation scheme and scores it, and (4.8.31) is the one case where its answer is exact. Chapter 4.11 §4 runs §3's commutator on angular momentum and closes the ladder at both ends with §4's argument used twice, and Chapter 4.13 §4 runs it a third time on hydrogen's radial equation. Chapter 4.15 perturbs this spectrum, using Worked example 2's matrix elements and Worked example 3's exact answer as its score. Chapters 4.16 and 4.17 spend §5.6's parity. Chapter 5.3 attaches one copy of §§2 to 4 to every mode of a field, and the rung number becomes the number of particles, which is the sentence Chapter 0.8 wrote and which needed this chapter's algebra before it could be made good. Chapter 7.3 builds the Virasoro algebra out of the same commutators and Chapter 7.4 quantises a string with it. What this chapter does not do is solve any potential that is not quadratic, and the honest scale of that limitation is in Worked example 1(d), where a real molecular bond departs from the ladder by two percent at the second rung.