Particle in a spherically symmetric potential

(Redirected from Spherical potential well)


In quantum mechanics, a spherically symmetric potential is a system of which the potential only depends on the radial distance from the spherical center and a location in space. A particle in a spherically symmetric potential will behave accordingly to said potential and can therefore be used as an approximation, for example, of the electron in a hydrogen atom or of the formation of chemical bonds.[1]

Hydrogen atomic orbitals of different energy levels. The more opaque areas are where one is most likely to find an electron at any given time.

In the general time-independent case, the dynamics of a particle in a spherically symmetric potential are governed by a Hamiltonian of the following form: Here, is the mass of the particle, is the momentum operator, and the potential depends only on the vector magnitude of the position vector, that is, the radial distance from the origin (hence the spherical symmetry of the problem).

To describe a particle in a spherically symmetric system, it is convenient to use spherical coordinates; denoted by , and . The time-independent Schrödinger equation for the system is then a separable, partial differential equation. This means solutions to the angular dimensions of the equation can be found independently of the radial dimension. This leaves an ordinary differential equation in terms only of the radius, , which determines the eigenstates for the particular potential, .

Structure of the eigenfunctions

edit

If solved by separation of variables, the eigenstates of the system will have the form:  in which the spherical angles   and   represent the polar and azimuthal angle, respectively. Those two factors of   are often grouped together as spherical harmonics, so that the eigenfunctions take the form:  

The differential equation which characterises the function   is called the radial equation.

Derivation of the radial equation

edit

The kinetic energy operator in spherical polar coordinates is: The spherical harmonics satisfy  Substituting this into the Schrödinger equation we get a one-dimensional eigenvalue equation, This equation can be reduced to an equivalent 1-D Schrödinger equation by substituting  , where   satisfies which is precisely the one-dimensional Schrödinger equation with an effective potential given by where  . The correction to the potential V(r) is called the centrifugal barrier term.

If  , then near the origin,  .

Spherically symmetric Hamiltonians

edit

Since the Hamiltonian is spherically symmetric, it is said to be invariant under rotation, ie:

 

Since angular momentum operators are generators of rotation, applying the Baker-Campbell-Hausdorff Lemma we get:

 

Since this equation holds for all values of  , we get that  , or that every angular momentum component commutes with the Hamiltonian.

Since   and   are such mutually commuting operators that also commute with the Hamiltonian, the wavefunctions can be expressed as   or   where   is used to label different wavefunctions.

Since   also commutes with the Hamiltonian, the energy eigenvalues in such cases are always independent of  .

 

Combined with the fact that   differential operators only act on the functions of   and  , it shows that if the solutions are assumed to be separable as  , the radial wavefunction   can always be chosen independent of   values. Thus the wavefunction is expressed as:[2]

 

Solutions for potentials of interest

edit

There are five cases of special importance:

  1.  , or solving the vacuum in the basis of spherical harmonics, which serves as the basis for other cases.
  2.   (finite) for   and zero elsewhere.
  3.   for   and infinite elsewhere, the spherical equivalent of the square well, useful to describe bound states in a nucleus or quantum dot.
  4.  for the three-dimensional isotropic harmonic oscillator.
  5.   to describe bound states of hydrogen-like atoms.

The solutions are outlined in these cases, which should be compared to their counterparts in cartesian coordinates, cf. particle in a box.

Vacuum case states

edit

Let us now consider  . Introducing the dimensionless variables the equation becomes a Bessel equation for  : where regular solutions for positive energies are given by so-called Bessel functions of the first kind   so that the solutions written for   are the so-called spherical Bessel function

 .

The solutions of the Schrödinger equation in polar coordinates in vacuum are thus labelled by three quantum numbers: discrete indices and m, and k varying continuously in  : These solutions represent states of definite angular momentum, rather than of definite (linear) momentum, which are provided by plane waves  .

Sphere with finite "square" potential

edit

Consider the potential   for   and   elsewhere - that is, inside a sphere of radius   the potential is equal to   and it is zero outside the sphere. A potential with such a finite discontinuity is called a square potential.[3]

We first consider bound states, i.e. states which display the particle mostly inside the box (confined states). Those have an energy   less than the potential outside the sphere, i.e., they have negative energy. Also worth noticing is that unlike Coulomb potential, featuring an infinite number of discrete bound states, the spherical square well has only a finite (if any) number because of its finite range.

The resolution essentially follows that of the vacuum case above with normalization of the total wavefunction added, solving two Schrödinger equations — inside and outside the sphere — of the previous kind, i.e., with constant potential. The following constraints must hold for a normalizable, physical wavefunction:

  1. The wavefunction must be regular at the origin.
  2. The wavefunction and its derivative must be continuous at the potential discontinuity.
  3. The wavefunction must converge at infinity.

The first constraint comes from the fact that Neumann and Hankel functions are singular at the origin. The physical requirement that   must be defined everywhere selected Bessel function of the first kind over the other possibilities in the vacuum case. For the same reason, the solution will be of this kind inside the sphere: Note that for bound states,  . Bound states bring the novelty as compared to the vacuum case now that  . This, along with the third constraint, selects the Hankel function of the first kind as the only converging solution at infinity (the singularity at the origin of these functions does not matter since we are now outside the sphere): The second constraint on continuity of   at   along with normalization allows the determination of constants   and  . Continuity of the derivative (or logarithmic derivative for convenience) requires quantization of energy.

Sphere with infinite "square" potential

edit

In case where the potential well is infinitely deep, so that we can take   inside the sphere and   outside, the problem becomes that of matching the wavefunction inside the sphere (the spherical Bessel functions) with identically zero wavefunction outside the sphere. Allowed energies are those for which the radial wavefunction vanishes at the boundary. Thus, we use the zeros of the spherical Bessel functions to find the energy spectrum and wavefunctions. Calling   the kth zero of  , we have: so that the problem is reduced to the computations of these zeros  , typically by using a table or calculator, as these zeros are not solvable for the general case.

In the special case   (spherical symmetric orbitals), the spherical Bessel function is  , which zeros can be easily given as  . Their energy eigenvalues are thus:  

3D isotropic harmonic oscillator

edit

The potential of a 3D isotropic harmonic oscillator is   An N-dimensional isotropic harmonic oscillator has the energies   i.e.,   is a non-negative integral number;   is the (same) fundamental frequency of the   modes of the oscillator. In this case  , so that the radial Schrödinger equation becomes,  

Introducing   and recalling that  , we will show that the radial Schrödinger equation has the normalized solution,   where the function   is a generalized Laguerre polynomial in   of order  .

The normalization constant   is,  

The eigenfunction   is associated with energy  , where   This is the same result as the quantum harmonic oscillator, with  .

Derivation

edit

First we transform the radial equation by a few successive substitutions to the generalized Laguerre differential equation, which has known solutions: the generalized Laguerre functions. Then we normalize the generalized Laguerre functions to unity. This normalization is with the usual volume element r2 dr.

First we scale the radial coordinate   and then the equation becomes   with  .

Consideration of the limiting behavior of v(y) at the origin and at infinity suggests the following substitution for v(y),   This substitution transforms the differential equation to   where we divided through with  , which can be done so long as y is not zero.

Transformation to Laguerre polynomials
edit

If the substitution   is used,  , and the differential operators become   and  

The expression between the square brackets multiplying   becomes the differential equation characterizing the generalized Laguerre equation (see also Kummer's equation):   with  .

Provided   is a non-negative integral number, the solutions of this equations are generalized (associated) Laguerre polynomials  

From the conditions on   follows: (i)   and (ii)   and   are either both odd or both even. This leads to the condition on   given above.

Recovery of the normalized radial wavefunction
edit

Remembering that  , we get the normalized radial solution:  

The normalization condition for the radial wave function is:  

Substituting  , gives   and the equation becomes:  

By making use of the orthogonality properties of the generalized Laguerre polynomials, this equation simplifies to:  

Hence, the normalization constant can be expressed as:  

Other forms of the normalization constant can be derived by using properties of the gamma function, while noting that   and   are both of the same parity. This means that   is always even, so that the gamma function becomes:   where we used the definition of the double factorial. Hence, the normalization constant is also given by:  

Hydrogen-like atoms

edit

A hydrogenic (hydrogen-like) atom is a two-particle system consisting of a nucleus and an electron. The two particles interact through the potential given by Coulomb's law:   where

In order to simplify the Schrödinger equation, we introduce the following constants that define the atomic unit of energy and length:

 

where   is the reduced mass in the   limit. Substitute   and   into the radial Schrödinger equation given above. This gives an equation in which all natural constants are hidden,   Two classes of solutions of this equation exist:

(i)   is negative, the corresponding eigenfunctions are square-integrable and the values of   are quantized (discrete spectrum).

(ii)   is non-negative, every real non-negative value of   is physically allowed (continuous spectrum), the corresponding eigenfunctions are non-square integrable. Considering only class (i) solutions restricts the solutions to wavefunctions which are bound states, in contrast to the class (ii) solutions that are known as scattering states.

For class (i) solutions with negative W the quantity   is real and positive. The scaling of  , i.e., substitution of   gives the Schrödinger equation:  

For   the inverse powers of x are negligible and the normalizable (and therefore, physical) solution for large   is  . Similarly, for   the inverse square power dominates and the physical solution for small   is x+1. Hence, to obtain a full range solution we substitute  

The equation for   becomes,  

Provided   is a non-negative integer, this equation has polynomial solutions written as   which are generalized Laguerre polynomials of order  . The energy becomes 

The principal quantum number   satisfies  . Since  , the total radial wavefunction is   with normalization which absorbs extra terms from     via[4]

 

The corresponding energy is  

References

edit
  1. ^ Ruedenberg, Klaus; Schmidt, Michael W. (2009-03-12). "Physical Understanding through Variational Reasoning: Electron Sharing and Covalent Bonding". The Journal of Physical Chemistry A. 113 (10): 1954–1968. doi:10.1021/jp807973x. ISSN 1089-5639. PMID 19228050.
  2. ^ Littlejohn, Robert G. "Physics 221A: Central Force Motion" (PDF). Archived (PDF) from the original on 8 December 2023. Retrieved 18 February 2024.
  3. ^ A. Messiah, Quantum Mechanics, vol. I, p. 78, North Holland Publishing Company, Amsterdam (1967). Translation from the French by G.M. Temmer
  4. ^ H. Margenau and G. M. Murphy, The Mathematics of Physics and Chemistry, Van Nostrand, 2nd edition (1956), p. 130. Note that convention of the Laguerre polynomial in this book differs from the present one. If we indicate the Laguerre in the definition of Margenau and Murphy with a bar on top, we have  .