Harmonic Oscillator Physics
Total Page:16
File Type:pdf, Size:1020Kb
Physics 342 Lecture 9 Harmonic Oscillator Physics Lecture 9 Physics 342 Quantum Mechanics I Friday, February 12th, 2010 For the harmonic oscillator potential in the time-independent Schr¨odinger equation: 2 1 2 d (x) 2 2 2 −~ + m ! x (x) = E (x); (9.1) 2 m dx2 we found a ground state 2 − m ! x 0(x) = A e 2 ~ (9.2) 1 with energy E0 = 2 ~ !. Using the raising and lowering operators 1 a+ = p (−i p + m ! x) 2 m ! ~ (9.3) 1 a− = p (i p + m ! x); 2 ~ m ! we found we could construct additional solutions with increasing energy using a+, and we could take a state at a particular energy E and construct solutions with lower energy using a−. The existence of a minimum energy state ensured that no solutions could have negative energy and was used to define 1: 0 n 1 n a− 0 = 0 H (a+ 0) = + n ~ ! a+ 0: (9.4) 2 The operators a+ and a− are Hermitian conjugates of one another { for any 1I am leaving the \hats" off, from here on { we understand that H represents a differ- ` ~ @ ´ ential operator given by H x; i @x for a classical Hamiltonian. 1 of 10 9.1. NORMALIZATION Lecture 9 f(x) and g(x) (vanishing at spatial infinity), the inner product: Z 1 Z 1 ∗ ∗ @g(x) f(x) a± g(x) dx = f(x) ∓~ + m ! x g(x) dx −∞ −∞ @x Z 1 ∗ @f(x) ∗ = ±~ g(x) + f(x) m ! x g(x) dx −∞ @x Z 1 ∗ = (a∓ f(x)) g(x) dx; −∞ (9.5) (integration by parts) or, in words, we can \act on g(x) with a± or act on f(x) with a∓" in the inner product. 9.1 Normalization n The state n that comes from n applications of a+ is not normalized, nor does the eigenvalue form of the time-independent Schr¨odingerequation de- mand that it be. Normalization is the manifestation of our probabilistic interpretation of jΨ(x; t)j2. Consider the ground state, that has an undeter- mined constant A. If we want jΨ(x; t)j2 to represent a probability density, then Z 1 2 2 − m ! x A e ~ dx = 1; (9.6) −∞ and the left-hand side is a Gaussian integral: Z 1 2 r 2 − m ! x 2 π ~ A e ~ dx = A ; (9.7) −∞ m ! 1=4 so the normalization is A = m ! . Just because we have normalized the π ~ ground state does not mean that 1 ∼ a+ 0(x) is normalized. Indeed, we have to normalize each of the n(x) separately. Fortunately, this can be done once (and for all). R 1 2 Suppose we have a normalized set of n, i.e. −∞ n dx = 1 for all n = 0 −! 1. We know that n±1 = α± a± n where the goal is to find the constants α± associated with raising and lowering while keeping the wavefunctions normalized. Take the norm of the resulting raised or lowered state: Z 1 Z 1 2 2 ∗ j n±1j dx = α± (a± n(x)) (a± n(x)) dx −∞ −∞ Z 1 (9.8) 2 ∗ = α± (a∓ a± n(x)) n(x) dx; −∞ 2 of 10 9.2. ORTHONORMALITY Lecture 9 and the operator a∓ a± is related to the Hamiltonian, as we saw last time: 1 H = ~ ! a± a∓ ± 2 . Then H 1 1 1 a∓ a± n = ± n = + n ± n; (9.9) ~ ! 2 2 2 so (9.8) becomes Z 1 Z 1 2 ∗ 2 1 1 2 α± (a∓ a± n(x)) n(x) dx = α± + n ± j n(x)j dx: −∞ 2 2 −∞ (9.10) R 1 2 The integral −∞ j n(x)j dx = 1 by assumption, so we have 1 1 α+ = p α− = p ; (9.11) n + 1 n and our final relation is 1 1 n+1 = p a+ n n−1 = p a− n: (9.12) n + 1 n Starting from the ground state, for which we know the normalization, 1 = a+ 0 1 1 2 2 = p a+ 1 = p a+ 0 2 1 × 2 (9.13) 1 1 3 3 = p a+ 2 = p a+ 0: 3 1 × 2 × 3 The general case is 1 n n = p a+ 0 (9.14) n! for the appropriately normalized 0, 1=4 2 m ! − m ! x 0(x) = e 2 ~ : (9.15) π ~ 9.2 Orthonormality The states described by n are complete (we assume) and orthonormal { take our usual inner product for m and n: Z 1 Z ∗ 1 1 n ∗ m n · m = n(x) m(x) dx = p p (a+ 0) (a+ 0) dx: −∞ m! n! −∞ (9.16) 3 of 10 9.2. ORTHONORMALITY Lecture 9 Now, suppose m > n, then using the fact that a+ and a− are Hermitian m conjugates, we can flip the a+ onto the other term: Z Z 1 1 n ∗ m 1 1 m n ∗ p p (a+ 0) (a+ 0) dx = p p (a− a+ 0) ( 0) dx; m! n! −∞ m! n! −∞ (9.17) n m and we know that a+ 0 ∼ n and a− n ∼ n−m { but for m > n, we m−n have a− 0 = 0, the defining property of the ground state. In the case n > m, we just use the conjugate in the other direction, and make the same argument: Z Z 1 1 n ∗ m 1 1 ∗ n m p p (a+ 0) (a+ 0) dx = p p ( 0) (a−a+ 0) dx = 0: m! n! −∞ m! n! −∞ (9.18) R 1 2 Only when m = n will we get a non-zero result, and of course, −∞ m dx = 1 by construction. So Z 1 ∗ n(x) m(x) dx = δmn: (9.19) −∞ 9.2.1 Hermite Polynomials The prescription for generating n does not provide a particularly easy way to obtain the functional form for an arbitrary n { we have to repeatedly apply the raising operator to the ground state. There is a connection between the Hermite polynomials and our procedure of \lifting up" the ground state. Using the Frobenius method, it is possible to solve Schr¨odinger'sequation as a power series expansion (described in Griffiths), and we won't re-live that argument (yet). But it is important to understand the connection between the algebraic, operational approach and the brute force series expansion. Starting from the ground state, let's act with a+ (call the normalization constant A again, just to make the expressions more compact): 2 1 @ − m ! x 1 = a+ 0 = p −~ + m ! x A e 2 ~ 2 m ! @x ~ (9.20) r 2 A m ! − m ! x = p 2 x e 2 ~ : 2 ~ 4 of 10 9.2. ORTHONORMALITY Lecture 9 Consider the second state, r 2 1 1 1 @ A m ! − m ! x 2 = p a+ 1 = p p −~ + m ! x p 2 x e 2 ~ 2 2 2 m ! ~ @x 2 ~ 2 A − m ! x @ 2 x = p e 2 ~ −~ 2 2 @x ~ r 2 1 m ! 1 @ A − m ! x + p 2 x p −~ + m ! x p e 2 ~ 2 ~ 2 m ! ~ @x 2 2 2 r 2 A − m ! x A m ! − m ! x = p e 2 ~ (−2) + p 2 x e 2 ~ 2 2 2 2 ~ 2 A m ! 2 − m ! x = p 4 x − 2 e 2 ~ : 2 2 ~ (9.21) The pattern continues { we always have some polynomial in x multiplying the exponential factor. That polynomial, for the nth wave function is called th p m ! Hn, the n Hermite polynomial. In the dimensionless variable ξ = x, 2 ~ we can read off the first two { H1(ξ) = 2 ξ, and H2(ξ) = 4 ξ −2. Normalized, we have the expression 1=4 2 m ! 1 − ξ n(x) = p Hn(ξ) e 2 (9.22) π ~ 2n n! with H0(ξ) = 1 H1(ξ) = 2 ξ 2 H2(ξ) = 4 ξ − 2 (9.23) 3 H3(ξ) = 8 ξ − 12 ξ . This set of polynomials is well-known, and they have a number of interesting recursion and orthogonality properties (many of which can be developed from the a+ and a− operators). One still needs a table of these in order to write down a particular n, but that's better than taking n successive derivatives of 0 { in essence, the Hermite polynomials have accomplished that procedure for you. Once again, we can plot the first few wavefunctions (see Figure 9.1), and as we increase in energy, we see a pattern similar to the infinite square well case (note that for the harmonic oscillator, we start with n = 0 as the ground state rather than 1). 5 of 10 9.3. EXPECTATION VALUES Lecture 9 n =3 n =2 Energy n =1 n =0 Figure 9.1: The first four stationary states: n(x) of the harmonic oscillator. 9.3 Expectation Values 9.3.1 Classical Case The classical motion for an oscillator that starts from rest at location x0 is x(t) = x0 cos(! t) : (9.24) The probability that the particle is at a particular x at a particular time t is given by ρ(x; t) = δ(x − x(t)), and we can perform the temporal average to get the spatial density. Our natural time scale for the averaging is a half π cycle, take t = 0 ! ! , π 1 Z ! ρ(x) = π δ(x − x0 cos(! t)) dt: (9.25) ! 0 We perform the change of variables to allow access to the δ, let y = x0 cos(! t) so that ! Z −x0 δ(x − y) ρ(x) = − dy π x0 x0 ! sin(! t) 1 Z x0 δ(x − y) = dy p 2 π −x0 x0 1 − cos (! t) (9.26) 1 Z x0 δ(x − y) = dy p 2 2 π −x0 x0 − y 1 = : p 2 2 π x0 − x 6 of 10 9.3.