San José State University 

appletmagic.com Thayer Watkins Silicon Valley & Tornado Alley USA 


Consider a particle moving in a central field whose force function is F(r), where r is the distance of the particle from the center of the field. The force is directed along the radial vector from the particle to the center of the field. Let p be the linear momentum of the particle.
The potential energy function V(r) for the field is given by
The Hamiltonian function H for the particle is its total energy; i.e., the sum of its kinetic energy K and its potential energy V. For a single particle in a potential field V it is then given by
where m is the mass of the particle.
The trajectory for a point particle can be determined using classical mechanics methods. The wave function for a particle can be determined using quantum mechanical methods.
Although the problem may be analyzed using a general central force function F(r) it should be noted that there are only about two relevant cases. The field is carried by some particles such as photons in the case of the electrostatic force and pi mesons in the case of the nuclear strong force. This means there will be a (1/r²) dependence. If the forcecarrying particles do not decay, as is the case for photons, then F(r) is of the form K/r² where K is a constant. If they do decay, as is the case with the pi mesons, then the surviving proportions is a negative exponential function of time which translates into a negative exponential function of distance from the source; i.e.,
where H and r_{0} are constants. Generally there should be constant factors Qq in the equations, where Q is the charge of the force field and q is the charge of the particle with respect to the field.
The Hamiltonian operator is constructed from the Hamiltonian function for a system by substituting
−hi∇ for the momentum p, where ∇ is the gradient operator for the coordinate system, i is
the imaginary unit and h is Planck's constant divided by 2π.
Thus for a particle in a central field with a potential energy function V(r)
The expresssion ∇² is called the Laplacian operator of the system.
The wave function ψ(r,θ,φ,t) is such that the product of it times its complex conjugate ψ^{c} is the the probability density of finding the particle at (r,θ,φ) at time t. This means that the integral of ψ*ψ^{c} over all space is always equal to unity.
The wave function satisfies the equation
Suppose ψ(r,θ,φ,t) is of the form Ψ(r,θ,φ)T(t). Then
The lefthand side (LHS) of the above equation is a function of the space variables r, φ and θ only. The righthand side (RHS) is a function of t alone. Obviously the common value of the two sides is a constant, which will be denoted as E.
Thus, from the right side of the equation,
The solution to the above equation is
This is a periodic solution with a frequency equal to (E/h).
From the LHS of the previous equation
Thus Ψ is an eigenfunction of the Hamiltonian operator and E is an eigenvalue. The above equation is called the timeindependent Schrödinger equation for the system.
It is assumed that Ψ(r,θ,φ)=R(r)Θ(θ)Φ(φ). When this form is substituted into the the previous equations and the differentiation carried out, the result is
If this equation is first divided by RΘΦ and then multiplied by
−(2mr²sin²(θ)/h² the result can be expressed as
This is an equation in which the LHS depends only upon φ and the RHS does not depend upon φ. Therefore the common value is a constant. Let this constant be denoted as M. Thus
This equation has the solution
Because φ is an angle it must be that
This cannot occur for M^{½} being real. Therefore M^{½} must be imaginary, say im. This means that
which requires that m be an integer, positive, negative or zero. Therefor M=−m².
An equation that results from the RHS of the previous equation being equal to M=−m² is
The LHS of the above equation depends only upon r and the RHS only upon θ. Let the common value be denoted as Λ. There are then the two equations
It turns out that for the second equation to have a physically acceptable solution it must be that Λ is of the for form l(l+1) where l is of the form m+k, for k being a nonnegative integer. For the derivation of this result see Quantum Eigenvalues.
(To be continued.)
HOME PAGE OF Thayer Watkins 