The Quantum Harmonic Oscillator Stephen Webb

Transcription

1 The Quantum Harmonic Oscillator Stephen Webb The Importance of the Harmonic Oscillator The quantum harmonic oscillator holds a unique importance in quantum mechanics, as it is both one of the few problems that can really be solved in closed form, and is a very generally useful solution, both in approximations and in exact solutions of various problems. The harmonic oscillator is characterized by the Hamiltonian: H = P2 2 m 1 2 m 2 X 2 This Hamiltonian appears in various applications, and in fact the approximation of the harmonic oscillator is valid near the minimum of any potential function. Expanded around a minimum point x*, any potential can then be Taylor expanded as: V x = V 0 x x * V x x = x * 1 2 x x 2 V * 2 x 2 x = x * In this case, V 0 is an arbitrary constant that can be taken to zero, and the second term is zero since the derivative is zero at critical points. If we assume that the variation around the critical point is small, then the higher order terms may be neglected, and we may readjust the variables to place the critical point at 0 and approximate it as a parabola. The harmonic oscillator also gives the exact solution for a particle in a uniform magnetic field of a given vector potential, as that vector potential merely takes the form of a two-dimensional harmonic oscillator. Finally, it serves as an excellent pedagogical tool. It introduces people to the methods of analytically solving the differential equations frequently encountered in quantum mechanics, and also provides a good introduction to the use of raising and lowering operators, and using the abstract vectors that are frequently used in quantum mechanics to solve problems by knowing the action of operators upon state vectors, rather than using integration to evaluate expectation values. The Classical Oscillator: A Review The classical harmonic oscillator is most frequently introduced as a mass on an undamped spring. The Hamiltonian for such a system is: and from the canonical relations we find that: H p = ẋ = p m H = p2 2 m 1 2 k x2 H x = ṗ = kx Solving this set of differential equations then gives the solution: x t = x 0 cos k m t It is this solution that we should approximately get in the high energy limit solution to the quantum harmonic oscillator, and this will be our test that we have found the solution to the problem.

2 The Quantum Harmonic Oscillator: Analytical Solution With the quantum harmonic oscillator we are presented with the problem of finding the eigenfunctions of the given Hamiltonian, which, in the position representation, is: H = ħ2 2 2 m x m 2 x 2 The Schrodinger equation then reads: ħ m x m 2 x 2 = i ħ t This is a second order, nonhomogeneous differential equation, for which we may apply separation of variables to turn it into the standard time independent Schrodinger equation: ħ m x m 2 x 2 = E Our first instinct in solving this problem should be to use a power series solution, as we know that the wave functions must be well behaved enough to be expressed as a power series, and the problem itself demands the use of it. First, we get the problem in unitless dimensions, doing so by making the following substitutions and dividing through by constants: = E ħ y = m x ħ Making these substitutions, we then get a unitless differential equation: ' ' 2 y 2 = 0 To find the solution to this differential equation, we first analyze its asymptotics, that is, the extreme positions of the differential equation. In this case, we look at small energy, when the y^2 term dominates. In this form, the differential equation is approximately: ' ' = y 2 which suggests that the wave function may be a Gaussian, or something of similar form. Thus, we guess that the wave function is of the form: x = u y exp[ y 2 /2] Plugging this into the differential equation we get: u ' ' y exp[ y 2 /2] 2 u ' y y exp[ y 2 /2] 2 1 u y exp[ y 2 /2] = 0 We can get rid of the Gaussian terms, as they represent a general multiplying factor, thus: u ' ' 2 y u ' 2 1 u = 0 If we assume that the equation for u is in the form of a power series, we make the substitution: u y = C n y n C n [n n 1 y n 2 2 ny n 2 1 y n ] = 0 n=0 n=0 Analyzing the resulting series solution, we arrive at the recursion relation: C n 2 = C n 2 n 1 2 n 2 n 1 In theory we've solved the problem, except for the normalization condition. In the limit of large n, we find that the ratio between the succeeding terms is n/2, which grows faster than the terms in the exponential. This series diverges faster than our Gaussian converges, and the only way for the series to converge, and thus for the wave function to have any physical meaning, is for the series to truncate. The best way to do this is to make the numerator zero, as any successive terms to that would be zero as well, and from this we get: 2 = 2 n 1 = n 1 2 E = ħ n 1 2 From the requirement that the power series converge, and therefore truncate, we find that the energy levels for the harmonic oscillator are quantized.

3 This power series solution is not complete, as the individual eigenstates of the Hamiltonian are not yet made orthogonal. When this is carried out, what results is a solution that involves the Hermite polynomials, and is of the form: n x = m 1/4 ħ ħ x] x2 exp m 2 ħ 1 2 n n! n[ H m The computation of expectation values with this involves knowledge of the Hermite polynomials, and involves integrals of Gaussians that become tedious and leave open the possibility for algebraic error. It also does not take much advantage of the tools of linear algebra that can be applied to quantum mechanics. An alternative method, developed by Paul A.M. Dirac, utilizes a different method, in which the Hamiltonian is factored into raising and lowering operators to acquire the different eigenstates of the Hamiltonian, and will be discussed now. The Quantum Harmonic Oscillator: Energy Basis Method The first thing to do is to generate a unitless form of the Hamiltonian, which we do by introducing a new operator: H ' ħ =H This makes the new operator measure the energy states in terms of the characteristic energies, and removes several constants that would otherwise appear in the process of solving the problem, and such a substitution is made with a priori knowledge of the results that are to follow. The substitution the makes the new Hamiltonian operator: H '= P2 2 m ħ 1 m 2 ħ X 2 We now introduce two new operators, what we will eventually recognize as the raising and lowering operators: a = m 1 X P 2 ħ i 2 m ħ a = m 1 X P 2 ħ i 2 m ħ The new Hamiltonian operator we created can be expressed in terms of these operators as: a a = m 2 ħ X m ħ P 2 i 2 ħ [ X, P ] a a = H ' 1 2 H ' = a a 1 2 These operators then satisfy the commutation relations [a,a ] = 1 [a, H ' ] = a [a, H ' ] = a Such commutation relations we will find useful, and these relations indicate that these operators produce new eigenstates of the Hamiltonian: H ' a n = a H ' [a, H ' ] n H ' a n = a H ' a n = n 1 a n The indication of this is that the operator a, which we will now refer to as the lowering operator, generates some new eigenstate for the Hamiltonian, which is one energy level below the previous, to within a normalization constant. A similar procedure can be done to show that a has a similar effect, only raising the eigenenergy of the state. So here we have a mechanism for determining all future eigenstates of the Hamiltonian given only one of them. Theoretically we have the capacity to generate an infinite number of energy states just by

4 repeatedly applying these raising and lowering operators, but we require that the energies be positive just because of our choice of potential and the fact that the state vector must be zero below the minimum potential energy. Therefore, we conclude that for some eigenstate of the Hamiltonian: a 0 = 0 The next question to ask is what the minimum energy of this eigenstate is. We can do this simply by exploiting the linearity of the raising and lowering operators: a 0 = 0 a a 0 =0 a a = H ' 1 2 H ' =0 H ' 0 = H 0 = 1 2 ħ 0 The minimum energy of the harmonic oscillator is ħ, which is exactly what we predicted using the power series method to solving the oscillator. Now we need to find the energy level of the n th eigenstate for the Hamiltonian. To do this, we use the fact that any eigenstate can be represented as the ground state acted upon n times by the raising operator. Doing so, we find that the energy difference between any two energy eigenstates is ħ which then gives us that the energy levels are E n = n ħ, which is exactly the result we gained using the power series method, but without the messy algebra. Now all that is left for using the raising and lowering operators is to find the normalization constant sitting in front of the newly created eigenstates generated by the raising and lowering operators. First, we consider the equation a n = C n n 1 and its complex conjugate: n a = n 1 C n * We combine these two equations to get: n a a n = n 1 n 1 C n * C n n H ' 1 2 n = C * nc n n n n = C n 2 C n 2 = n C n = n e i where the φ in the above equation is an undetermined phase constant usually taken as zero by convention, as it is an unobservable rotation in the energy eigenspace. A similar procedure gives us that a n = n 1 n 1 From all this we can then find that any eigenstate of our Hamiltonian can be expressed as: n = a n n! 0 To recap what we have done so far, we can express everything about the energy of the eigenstates of the harmonic oscillator as follows: H n = n 1 ħ n 2 a n = n n 1 a n = n 1 n 1 We then have a nice mechanism for computing the eigenstates of the Hamiltonian, but we can also express expectation values using the raising and lowering operators. By simply making the substitutions:

5 ħ a a 2 m P = i m ħ 2 a a we can compute any arbitrary expectation values that depend upon these quantities, merely by knowing the effects of the raising and lowering operators upon the eigenstates of the Hamiltonian. Remarks on the Oscillator X = There are several items of note that should be verified by the reader. First of all, the ground state wave function of the harmonic oscillator represents the minimum uncertainty state, for which the Heisenberg uncertainty inequality for momentum and position x p ħ 2 is minimized, that is, the above expression is a strict equality rather than an inequality. This can be easily confirmed using the raising and lowering operators to find that x = 0 p = 0 x 2 = ħ 2 m p2 = m ħ 2 The harmonic oscillator also has its uses in a uniform magnetic field, in which the vector potential and revised Hamiltonian can be written: A = B 0 y x x y 2 H = 1 2 m p q A 2 When the new Hamiltonian is expanded, you get a two dimensional harmonic oscillator in the x-y coordinates and a free particle in the z coordinates. Einstein utilized the 3N uncoupled harmonic oscillators to approximate the heat capacity of a crystal lattice, and approximation far more accurate than the classical one but corrected upon by Debye to deal with the fact that all the oscillators are actually coupled, and provides a good solution for low temperatures. These last two examples can be worked out on your own without much difficulty; the magnetic field example can be done by simply plugging into the new Hamiltonian, whereas the heat capacity problem is an example that can be found in several textbooks. The Three Dimensional Oscillator For a three dimensional oscillator, the solutions are the same as three one-dimensional oscillators added together. Each raising and lowering operator has to be specified to coordinates, and the eigenstates are given by H n x n y n z = n x n y n z 3/2 ħ n x n y n z for an isotropic oscillator. This oscillator experiences various degeneracies on the level of the number of solutions to n x n y n z = k which is given by the formula D k = k 3 1 k Such degeneracies require a study in identical particles to resolve further. Like the one-dimensional harmonic oscillator, we have raising and lowering operators, but these

6 operators act specifically on the different spatial coordinates, thus, we have the raising and lowering operators a x, a x, a y, a y, a z, a z so that their individual actions on an eigenstate are given by a i n i n j n k = n i n i 1 n j n k a i n i n j n k = n i 1 n i 1 n j n k a i a i n i n j n k = n i n i n j n k while the individual operators satisfy their own commutation relations for the same spatial coordinate, but commute with the operators of different coordinates: [a i, a j ] = i, j This three dimensional quantum oscillator was used by Albert Einstein in 1907 when he provided a very primitive and crude explaination of the heat capacity of crystalline solids from the viewpoint of each solid being held in place as a quantum harmonic oscillator. Although this description lacks the sophistication required to provide a legitimate theory explaining the phenomena (it neglects the fact that all the oscillators are coupled), it still provided a qualitative approach which Debye followed with greater success. The Harmonic Oscillator Approximation for the Minima of Potentials As mentioned, in the introduction, the quantum harmonic oscillator is useful as it provides an approximation of the minima of any potential well. To illustrate this, we will analyze the quantum states given by the Hamiltonian H = p2 2 m V 0cosh x/a Our first job is to find the critical point of the potential x V 0cosh x/a = V 0 sinh x/a = 0 x = 0 a Now we have to Taylor expand the potential function about zero, and use zero to find the approximate value for. The Taylor expansion is given by: V 0 cosh x/a = V 0 1 sinh0 x a cosh 0 2 x2 O x 3 The first constant can be disregarded as a scaling factor that can be taken into account with the energy states, finding the total energy by adding back one V 0. Therefore we find that the Hamiltonian can be approximately expressed as: H p2 2 m 1 2 V 0 a 2 x2 V 0 From this and keeping in mind that the scalar in front of the potential Taylor expansion gives: E n n 1 2 ħ V 0 m a 2 V 0 Keeping this in mind, we therefore see that we can use the harmonic oscillator to give a partial approximation to the first few energy states for more esoteric potentials for which we have no hope of finding a closed form solution for.

Chapter 3 (Lecture 4-5) Postulates of Quantum Mechanics Now we turn to an application of the preceding material, and move into the foundations of quantum mechanics. Quantum mechanics is based on a series

The Essentials of Quantum Mechanics Prof. Mark Alford v7, 2008-Oct-22 In classical mechanics, a particle has an exact, sharply defined position and an exact, sharply defined momentum at all times. Quantum

Quantum Mechanics Bohr's model: BUT: In 1925-26: by 1930s: - one of the first ones to use idea of matter waves to solve a problem - gives good explanation of spectrum of single electron atoms, like hydrogen

CHAPTER 5 THE HARMONIC OSCILLATOR The harmonic oscillator is a model which has several important applications in both classical and quantum mechanics. It serves as a prototype in the mathematical treatment

CHAPTER 4 PRINCIPLES OF QUANTUM MECHANICS In this Chapter we will continue to develop the mathematical formalism of quantum mechanics, using heuristic arguments as necessary. This will lead to a system

The Schrödinger Equation When we talked about the axioms of quantum mechanics, we gave a reduced list. We did not talk about how to determine the eigenfunctions for a given situation, or the time development

PHYS 1624 Physics I An introduction to mechanics, heat, and wave motion. This is a calculus- based course for Scientists and Engineers. 4 hours (3 lecture/3 lab) Prerequisites: Credit for MATH 2413 (Calculus

Lecture 2 Observables 13 14 LECTURE 2. OBSERVABLES 2.1 Observing observables We have seen at the end of the previous lecture that each dynamical variable is associated to a linear operator Ô, and its expectation

ecture VI Abstract Before learning to solve partial differential equations, it is necessary to know how to approximate arbitrary functions by infinite series, using special families of functions This process

Math 2280 Section 002 [SPRING 2013] 1 Today well learn about a method for solving systems of differential equations, the method of elimination, that is very similar to the elimination methods we learned

Chapter 41. One Dimensional Quantum Mechanics Quantum effects are important in nanostructures such as this tiny sign built by scientists at IBM s research laboratory by moving xenon atoms around on a metal

An Introduction to Separation of Variables with Fourier Series Math 391w, Spring 2010 Tim McCrossen Professor Haessig Abstract: This paper aims to give students who have not yet taken a course in partial

Fourier Series A Fourier series is an infinite series of the form a b n cos(nωx) c n sin(nωx). Virtually any periodic function that arises in applications can be represented as the sum of a Fourier series.

arxiv:physics/0004029v1 [physics.ed-ph] 14 Apr 2000 Lagrangians and Hamiltonians for High School Students John W. Norbury Physics Department and Center for Science Education, University of Wisconsin-Milwaukee,

4. The Infinite Square Well Copyright c 215 216, Daniel V. Schroeder In the previous lesson I emphasized the free particle, for which V (x) =, because its energy eigenfunctions are so simple: they re the

5.6 Physical Chemistry 5 Helium Atom page HELIUM ATOM Now that we have treated the Hydrogen like atoms in some detail, we now proceed to discuss the next simplest system: the Helium atom. In this situation,

Chapter 5 Harmonic Oscillator and Coherent States 5. Harmonic Oscillator In this chapter we will study the features of one of the most important potentials in physics, it s the harmonic oscillator potential

APPLIED MATHEMATICS ADVANCED LEVEL INTRODUCTION This syllabus serves to examine candidates knowledge and skills in introductory mathematical and statistical methods, and their applications. For applications

Average Particle Energy in an Ideal Gas - the total energy of the system is found by summing up (integrating) over all particles n(ε) at different energies ε - with the integral - we find - note: - the

Second Order Linear Nonhomogeneous Differential Equations; Method of Undetermined Coefficients We will now turn our attention to nonhomogeneous second order linear equations, equations with the standard

MASTER OF SCIENCE IN PHYSICS Admission Requirements 1. Possession of a BS degree from a reputable institution or, for non-physics majors, a GPA of 2.5 or better in at least 15 units in the following advanced

Dynamics Any physical system, such as neurons or muscles, will not respond instantaneously in time but will have a time-varying response termed the dynamics. The dynamics of neurons are an inevitable constraint

Page 1 Bra-ket notation FromWikipedia,thefreeencyclopedia Bra-ket notation is the standard notation for describing quantum states in the theory of quantum mechanics. It can also be used to denote abstract

Taylor and Maclaurin Series In the preceding section we were able to find power series representations for a certain restricted class of functions. Here we investigate more general problems: Which functions

CHAPTER I VECTOR QUANTITIES Quantities are anything which can be measured, and stated with number. Quantities in physics are divided into two types; scalar and vector quantities. Scalar quantities have

2. Spin Chemistry and the Vector Model The story of magnetic resonance spectroscopy and intersystem crossing is essentially a choreography of the twisting motion which causes reorientation or rephasing

Chapter 7 ket and bra notation Introduces the bra and ket notation and gives some examples of its use. When you change the description of the world from the inutitive and everyday classical mechanics to

Chapter 3 Operator methods in quantum mechanics While the wave mechanical formulation has proved successful in describing the quantum mechanics of bound and unbound particles, some properties can not be

Math 5311 Gateaux differentials and Frechet derivatives Kevin Long January 26, 2009 1 Differentiation in vector spaces Thus far, we ve developed the theory of minimization without reference to derivatives.

Chapter 3 Fourier Series 3.1 Some Properties of Functions 3.1.1 Goal We review some results about functions which play an important role in the development of the theory of Fourier series. These results

Oscillations Vern Lindberg June 10, 2010 You have discussed oscillations in Vibs and Waves: we will therefore touch lightly on Chapter 3, mainly trying to refresh your memory and extend the concepts. 1

Lecture-XXIV Quantum Mechanics Expectation values and uncertainty Expectation values We are looking for expectation values of position and momentum knowing the state of the particle, i,e., the wave function

Quantum Mechanical Operators and Commutation C I. Bra-Ket Notation It is conventional to represent integrals that occur in quantum mechanics in a notation that is independent of the number of coordinates

Math Department Student Learning Objectives Updated April, 2014 Institutional Level Outcomes: Victor Valley College has adopted the following institutional outcomes to define the learning that all students

Basic Quantum Mechanics Postulates of QM - The state of a system with n position variables q, q, qn is specified by a state (or wave) function Ψ(q, q, qn) - To every observable (physical magnitude) there

Introduction to quantum mechanics Lecture 3 MTX9100 Nanomaterjalid OUTLINE -What is electron particle or wave? - How large is a potential well? -What happens at nanoscale? What is inside? Matter Molecule

State of Stress at Point Einstein Notation The basic idea of Einstein notation is that a covector and a vector can form a scalar: This is typically written as an explicit sum: According to this convention,

Lecture Relevant sections in text: 3.1, 3. Rotations in quantum mechanics Now we will discuss what the preceding considerations have to do with quantum mechanics. In quantum mechanics transformations in

Chapter 9 Unitary Groups and SU(N) The irreducible representations of SO(3) are appropriate for describing the degeneracies of states of quantum mechanical systems which have rotational symmetry in three

Chapter 18 Statistical mechanics of classical systems 18.1 The classical canonical partition function When quantum effects are not significant, we can approximate the behavior of a system using a classical

From Einstein to Klein-Gordon Quantum Mechanics and Relativity Aline Ribeiro Department of Mathematics University of Toronto March 24, 2002 Abstract We study the development from Einstein s relativistic

STRAND 1: QUANTITATIVE LITERACY AND LOGIC STANDARD L1: REASONING ABOUT NUMBERS, SYSTEMS, AND QUANTITATIVE SITUATIONS Based on their knowledge of the properties of arithmetic, students understand and reason

Accuracy of the coherent potential approximation for a onedimensional array with a Gaussian distribution of fluctuations in the on-site potential I. Avgin Department of Electrical and Electronics Engineering,

UNIVERSITY OF LONDON BSc/MSci EXAMINATION June 2007 for Internal Students of Imperial College of Science, Technology and Medicine This paper is also taken for the relevant Examination for the Associateship

Advanced Higher Mathematics Course Assessment Specification (C747 77) Valid from August 2015 This edition: April 2016, version 2.4 This specification may be reproduced in whole or in part for educational

Angular Velocity and Momentum Kurt Nalty September 2, 2008 Introduction The presentation of angular momentum in the classical physics world is dominated by rigid body mechanical concepts. This, in turn,

Algebra 1 Chapter 7 Notes Name Section 7.1 Solving Linear Systems by Graphing System of Linear Equations: Two or more equations in the same variables, also called a. Solution of a System of Linear Equations:

Big Ideas in Mathematics which are important to all mathematics learning. (Adapted from the NCTM Curriculum Focal Points, 2006) The Mathematics Big Ideas are organized using the PA Mathematics Standards

A Theoretical Model for Mutual Interaction between Coaxial Cylindrical Coils Lukas Heinzle Page 1 of 15 Abstract: The wireless power transfer link between two coils is determined by the properties of the

MATH022 Linear Algebra Brief lecture notes 48 Similarity and Diagonalization Similar Matrices Let A and B be n n matrices. We say that A is similar to B if there is an invertible n n matrix P such that