Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

American Journal of Physics

The “Radial Kinetic Energy” Term in the Schrödinger Equation for a Central Force

Frank S. Crawford

Published Aug 1, 1964• 2 citations
Publisher

Abstract

A straightforward derivation is given for the Hermitian operator pr ≡ 1r ħi δδrr corresponding to the radial component of linear momentum in the central force problem. The purpose is purely pedagogical; we believe that the radial term of the Schrödinger equation in spherical coordinates is more comprehensible when it appears in the form (pr2/2m)ψ = [(1r ħi δδrr)2/2m]ψ than when it appears, as it does in most textbooks, in the form −ħ22mr2 δδr[r2δδr]ψ. The latter expression is usually arrived at by a tedious transformation of variable done “in the appendix”; instead one may derive pr with a minimum of algebra and with the focus where it belongs, on the concepts of quantum mechanics.

DOI
10.1119/1.1970875

Related papers

From Primitives to the Schrödinger Equation

Allen Proxmire

The Schrödinger equation iℏ ∂ₜψ = Ĥψ, with non‑relativistic Hamiltonian Ĥ = ℏ²|p̂|²/(2m) + V(x̂), has been quantum mechanics’s central dynamical statement for nearly a century. It governs atomic spectra, molecular bonding, semiconductors, neutron interferometry, and countless other empirical settings. The question is why it has the form it does. In standard QM, the Schrödinger equation is a postulate. The Hilbert space is given. The momentum operator p̂ = −iℏ∇ is given. The kinetic form |p̂|²/(2m) is inherited directly from classical mechanics — quantum theory simply adopts the classical T = p²/(2m) structure without deeper justification. The first‑order time derivative is chosen because empirically it works. Each piece is accepted because the resulting equation matches experiment. This walkthrough shows how the Event Density framework derives the Schrödinger equation, in its full standard form, from a smaller set of structural commitments about what reality is made of at the substrate level. Every element — the Hilbert space, the momentum operator, linearity, first‑order time structure, the kinetic‑plus‑potential form, the factor of 1/(2m) — emerges as a forced consequence rather than a separate postulate. The derivation runs through Stone’s theorem applied twice — once on spatial translation (forcing p̂ = −iℏ∇) and once on time translation (forcing existence and uniqueness of a self‑adjoint generator with automatic linearity and first‑order structure) — closed by Galilean Lie algebra integration that identifies the time‑translation generator with the kinetic‑plus‑potential operator. The factor of 1/(2m) emerges as the chain‑rule Jacobian of the Galilean commutator [Ĥ, K̂ᵢ] = −iℏ p̂ᵢ, not as a borrowing from classical mechanics. The first‑order time derivative is a consequence of operator exponentials, not a chosen feature.

•May 12, 2026•0 citations

Unexpected Delta-Function Term in the Radial Schrodinger Equation

Anzor Khelashvili, Teimuraz Nadareishvili

Careful exploration of the idea that equation for radial wave function must be compatible with the full Schrodinger equation shows appearance of the delta-function while reduction of full Schrodinger equation in spherical coordinates. Elimination of this extra term produces a boundary condition for the radial wave function, which is the same both for regular and singular potentials.

arXiv:1002.1278•Feb 5, 2010•0 citations

Non-Hermitian Radial Momentum Operator and Path Integrals in Polar Coordinates

Kazuo Fujikawa

A salient feature of the Schrödinger equation is that the classical radial momentum term p2r in polar coordinates is replaced by the operator , where the operator is not hermitian in general. This fact has important implications for the path integral and semi-classical approximations. When one defines a formal hermitian radial momentum operator , the relation holds in d-dimensional space and this extra potential appears in the path integral formulated in polar coordinates. The extra potential, which influences the classical solutions in the semi-classical treatment such as in the analysis of solitons and collective modes, vanishes for d = 3 and attractive for d = 2 and repulsive for all other cases d ≥ 4. This extra term induced by the non-hermitian operator is a purely quantum effect, and it is somewhat analogous to the quantum anomaly in chiral gauge theory.

arXiv:0805.3879•Progress of Theoretical Physics•Aug 1, 2008•9 citations

Vector operators and spherical harmonics in quantum mechanics

M. Andrews

The usual explicit expressions for spherical harmonics are derived by a method which does not involve the solution of differential equations or the use of angular momentum operators in spherical polar form. Instead a raising operator constructed from an arbitrary vector is used together with the usual angular momentum raising and lowering operators. The method is first presented as simply as possible in Schrödinger representation using position as the vector. Then it is generalized to an arbitrary vector and arbitrary integral angular momentum without using a particular representation. Finally it is shown how a conserved vector (other than angular momentum) of a central system can be used to reveal properties of the energy eigenstates. As an example, eigenstates of the hydrogen atom are found using the Runge-Lenz vector.

American Journal of Physics•Mar 1, 1979•5 citations

ORBITAL ANGULAR MOMENTUM

Mircea S. Rogalski, Stuart B. Palmer

The operator is contained in the expression (31.1) of the central field Hamiltonian, and this can be shown using Eq.(31.6), where the component operators are replaced by their forms (29.16) written in Cartesian coordinates, and this gives ⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−∂ ∂⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−∂ ∂−⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂−= z x x z z x x z y z z y y z z yL 222ˆ == ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂− x y y x x y y x2= It is a straightforward matter to show that this relation may also be written as ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂+∂ ∂+⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂+∂ ∂+⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂+∂ ∂−= 2 2 yx z xz y zy xL = ⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+∂ ∂+∂ ∂+∂∂ ∂+∂∂ ∂+∂∂ ∂+ z z y y x x xz zx zy yz yx xy 222222 222 2= We now use the expansion ( )2 2ˆ ˆr p x y z x y zx y z x y z⎛ ⎞⎛ ⎞∂ ∂ ∂ ∂ ∂ ∂⋅ = − + + + +⎜ ⎟⎜ ⎟∂ ∂ ∂ ∂ ∂ ∂⎝ ⎠⎝ ⎠G G = ⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+∂ ∂+∂ ∂+∂∂ ∂+∂∂ ∂+∂∂ ∂+∂ ∂+∂ ∂+∂ ∂−= z z y y x x xz zx zy yz yx xy z z y y x x so that the addition of the last two relations yields ( ) ( ) 2 2 222 2 2 2 2 22 2 2ˆ ˆLˆ r p x y z x y zx y zx y z⎛ ⎞ ⎛ ⎞∂ ∂ ∂ ∂ ∂ ∂+ ⋅ = − + + + + + + +⎜ ⎟ ⎜ ⎟∂ ∂ ∂∂ ∂ ∂ ⎝ ⎠⎝ ⎠G G = = ( )pripr ˆˆˆ 22 GG=G ⋅+= (31.10) It follows that the momentum can be expressed as ( ) ( ) ( ) rrr r r L r priprL r p ∂ ∂−⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−=⎥⎦ ⎤⎢⎣ ⎡ ⋅−⋅+= 2 ˆ1ˆˆˆˆˆ1ˆ ==GG=GGG (31.11) where the radial part results from the substitution of Eq. (I.20), Appendix I, which gives 1 1ˆ sinr p i i e e e r r rθ ϕθ θ ϕ ⎛ ⎞∂ ∂ ∂= − ∇ = − + +⎜ ⎟∂ ∂ ∂⎝ ⎠ G G G G= = and hence ˆ ˆr p i r r ∂⋅ = − ∂ G G = Therefore the Hamiltonian (31.1) becomes ( )rUL rrrr r rm H + ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ −∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−= 222 2 ˆ111 2 ˆ = = and reduces to the equivalent form ( )rUL rr r rrm H +⎥⎦ ⎤⎢⎣ ⎡ −⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ 2 ˆ11 2 ˆ = = (31.12) used for a particle bound by a central force field. 31.2. EIGENVALUE EQUATIONS IN A CENTRAL FIELD Spherically symmetric systems, where the potential energy ( )rU is independent of the direction of ,rG are usually treated using the spherical coordinates ( )ϕθ ,,r related to the Cartesian coordinates as referred to in Appendix I, Eq.(I.19). As the angle θ is defined by the inclination to the z-axis, it is mathematically convenient to take the direction of the component of ( zyx ,, ) L G that is compatible with (in the sense introduced in Chapter 34) to be the z-axis. Therefore we consider the eigenvalue equations of the operators and given by ˆ ˆ, zH L 2Lˆ ( ) ( )ˆ , , , ,H r rθ ϕ ε θ ϕΨ = Ψ ( ) ( )ϕθϕθ ,,,,ˆ rLrL zz Ψ=Ψ (31.13) ( ) ( )ϕθϕθ ,,,,ˆ 22 rLrL Ψ=Ψ The appropriate form of the Hamiltonian Hˆ is conveniently written using the standard expression (I.23) of the Laplace operator in spherical coordinates as 2∇ ( )rU rrr r rrm H +⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂−= 2 sin 1sin sin 11 2 ˆ ϕθθθθθ = (31.14) Using the transformation (I.19) from Cartesian to spherical coordinates and the chain rule for differentiation, we find sin sin sin cosx y z r r x y x y z x y y θ ϕ θ ϕϕ ϕ ϕ ϕ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂= + + = − + = −∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂x so that the operator takes the form zLˆ ϕ∂ ∂=⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂= ix y y x i Lz ==ˆ (31.15) Furthermore, from Eqs.(31.12) and (31.14) it follows that is in fact expressed by that part of Hamiltonian which operates on the angular variables 2Lˆ ,θ ϕ only, that is, ⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂−= 2 sin 1sin sin 1ˆ ϕθθθθθ=L (31.16) The energy eigenfunctions in central field, which are simultaneously eigenfunctions of and are subjected, in addition to the usual bound state requirement (30.21), to the single-valueness condition zLˆ 2ˆ ,L ( ) ( )πϕθϕθ 2,,,, +Ψ=Ψ rr since the points ϕ and πϕ 2+ are in fact the same physical point. The solution to the time-independent Schrödinger equation (31.13) in central field can be found by a process of separation of variables, yielding three ordinary differential equations in the variables ,r θ and ,ϕ respectively, to which a physical significance will be assigned. It is common practice to separate the variables in two stages, first substituting the trial solution ( ) ( ) ( ), , ,r R r Yθ ϕΨ = θ ϕ (31.17) into the energy eigenvalue equation. If the Hamiltonian is expanded in spherical coordinates, as given by Eq.(31.14), Eq.(31.13) becomes ( )2 22 2 2 21 1 1 sin2 2 sin sin 21R Y Yr U r r m R r r m Y ε θθ θ θ θ ϕ ⎡ ⎤⎡ ⎤∂ ∂ ∂ ∂ ∂⎛ ⎞ ⎛ ⎞⎡ ⎤+ − = − +⎢ ⎥⎜ ⎟ ⎜ ⎟⎢ ⎥ ⎣ ⎦∂ ∂ ∂ ∂ ∂⎝ ⎠ ⎝ ⎠⎣ ⎦ ⎣ ⎦ = = (31.18) One side is independent of θ and ,ϕ and the other side is independent of r, so they must be separately equal to some constant, which is conveniently chosen to be equal to In this way we obtain a radial equation, given by .2/2 mL ( ) ( ) ( )rRrR mr LrU dr dr dr d rm ε=⎥⎦ ⎤⎢⎣ ⎡ ++⎟⎠ ⎞⎜⎝ ⎛− 2 2 1 2 = (31.19) where was replaced by for obvious reasons, and an angular equation, having the so-called spherical harmonics r∂∂ / drd / ( ),Y θ ϕ as solutions, which reads 2 2 ˆor sin 1sin sin 1 ==⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂− ϕθθθθθ= (31.20) The form of the spherical harmonics is obtained as a result of the second stage of the separation process, where we assume that ( ) ( ) ( )ϕθϕθ Φ= PY , so that Eq.(31.20) becomes 2 2 1sinsinsin ϕθθθθ θ ∂ Φ∂ Φ−=⎥⎦ ⎤⎢⎣ ⎡ +⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂ === LP P (31.21) Since one side is independent of ϕ while the other is independent of ,θ they must each be equal to a constant, which we call We thus obtain the differential equations .2zL d L dϕ Φ 2− = Φ= (31.22) where again ϕ∂∂ / has been replaced by ϕdd / , and 0 sin sin sin 1 22 ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ −+⎟⎠ ⎞⎜⎝ ⎛ PLL d dP d d z θθθθθ === (31.23) where θ∂∂ / has been replaced by θdd / . Therefore the energy equation (31.13) completely separates, in central field, into a radial equation (31.19) and two angular equations (31.22) and (31.23). The radial equation, giving the energy levels ε of the system, takes the form ( ) ( ) ( ) ( )rRrR mr LrUrR dr d rdr d m ε=⎥⎦ ⎤⎢⎣ ⎡ ++⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ +− 2 2 2 2 = (31.24) which can be simplified by making the substitution ( ) ( )rrRr =ξ (31.25) and using the identity 2 1r dd d r dr r rdr dr ξ ξ⎛ ⎞+ =⎜ ⎟⎝ ⎠ r Thus Eq.(31.24) reduces to ( ) ( ) ( ) ( )rr mr LrU dr rd m εξξξ =⎥⎦ ⎤⎢⎣ ⎡ ++− 2 22 = (31.26) which is formally identical to the time-independent Schrödinger equation in one dimension (30.17), provided the central force potential ( )rU is replaced by an effective potential energy of the form ( ) ( ) 2 2 2mr LrUrUeff += The repulsive potential, associated with the angular momentum, modifies the central potential as illustrated before, see Figure 3.3, for ( ) rCrU /−= . Further progress with the solution of the radial equation depends on the particular form of the potential function . We shall consider the significant example of one-electron atoms in the next chapter.

CRC Press eBooks•Jan 8, 2014•1 citation

THE $sub 4$RADIAL KINETIC ENERGY$sub 4$ TERM IN THE SCHRODINGER EQUATION FOR A CENTRAL FORCE

Frank S. Crawford

•Aug 1, 1964•0 citations
A straightforward derivation is given for the Hermitian operator pr ≡ 1r ħi δδrr corresponding to the radial component of linear momentum in the central force problem. The purpose is purely pedagogical; we believe that the radial term of the Schrödinger equation in spherical coordinates is more comprehensible when it appears in the form (pr2/2m)ψ = [(1r ħi δδrr)2/2m]ψ than when it appears, as it does in most textbooks, in the form −ħ22mr2 δδr[r2δδr]ψ. The latter expression is usually arrived at by a tedious transformation of variable done “in the appendix”; instead one may derive pr with a minimum of algebra and with the focus where it belongs, on the concepts of quantum mechanics.
The Schrödinger equation iℏ ∂ₜψ = Ĥψ, with non‑relativistic Hamiltonian Ĥ = ℏ²|p̂|²/(2m) + V(x̂), has been quantum mechanics’s central dynamical statement for nearly a century. It governs atomic spectra, molecular bonding, semiconductors, neutron interferometry, and countless other empirical settings. The question is why it has the form it does. In standard QM, the Schrödinger equation is a postulate. The Hilbert space is given. The momentum operator p̂ = −iℏ∇ is given. The kinetic form |p̂|²/(2m) is inherited directly from classical mechanics — quantum theory simply adopts the classical T = p²/(2m) structure without deeper justification. The first‑order time derivative is chosen because empirically it works. Each piece is accepted because the resulting equation matches experiment. This walkthrough shows how the Event Density framework derives the Schrödinger equation, in its full standard form, from a smaller set of structural commitments about what reality is made of at the substrate level. Every element — the Hilbert space, the momentum operator, linearity, first‑order time structure, the kinetic‑plus‑potential form, the factor of 1/(2m) — emerges as a forced consequence rather than a separate postulate. The derivation runs through Stone’s theorem applied twice — once on spatial translation (forcing p̂ = −iℏ∇) and once on time translation (forcing existence and uniqueness of a self‑adjoint generator with automatic linearity and first‑order structure) — closed by Galilean Lie algebra integration that identifies the time‑translation generator with the kinetic‑plus‑potential operator. The factor of 1/(2m) emerges as the chain‑rule Jacobian of the Galilean commutator [Ĥ, K̂ᵢ] = −iℏ p̂ᵢ, not as a borrowing from classical mechanics. The first‑order time derivative is a consequence of operator exponentials, not a chosen feature.
A salient feature of the Schrödinger equation is that the classical radial momentum term p2r in polar coordinates is replaced by the operator , where the operator is not hermitian in general. This fact has important implications for the path integral and semi-classical approximations. When one defines a formal hermitian radial momentum operator , the relation holds in d-dimensional space and this extra potential appears in the path integral formulated in polar coordinates. The extra potential, which influences the classical solutions in the semi-classical treatment such as in the analysis of solitons and collective modes, vanishes for d = 3 and attractive for d = 2 and repulsive for all other cases d ≥ 4. This extra term induced by the non-hermitian operator is a purely quantum effect, and it is somewhat analogous to the quantum anomaly in chiral gauge theory.
The usual explicit expressions for spherical harmonics are derived by a method which does not involve the solution of differential equations or the use of angular momentum operators in spherical polar form. Instead a raising operator constructed from an arbitrary vector is used together with the usual angular momentum raising and lowering operators. The method is first presented as simply as possible in Schrödinger representation using position as the vector. Then it is generalized to an arbitrary vector and arbitrary integral angular momentum without using a particular representation. Finally it is shown how a conserved vector (other than angular momentum) of a central system can be used to reveal properties of the energy eigenstates. As an example, eigenstates of the hydrogen atom are found using the Runge-Lenz vector.
The operator is contained in the expression (31.1) of the central field Hamiltonian, and this can be shown using Eq.(31.6), where the component operators are replaced by their forms (29.16) written in Cartesian coordinates, and this gives ⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−∂ ∂⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−∂ ∂−⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂−= z x x z z x x z y z z y y z z yL 222ˆ == ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂− x y y x x y y x2= It is a straightforward matter to show that this relation may also be written as ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂+∂ ∂+⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂+∂ ∂+⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂+∂ ∂−= 2 2 yx z xz y zy xL = ⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+∂ ∂+∂ ∂+∂∂ ∂+∂∂ ∂+∂∂ ∂+ z z y y x x xz zx zy yz yx xy 222222 222 2= We now use the expansion ( )2 2ˆ ˆr p x y z x y zx y z x y z⎛ ⎞⎛ ⎞∂ ∂ ∂ ∂ ∂ ∂⋅ = − + + + +⎜ ⎟⎜ ⎟∂ ∂ ∂ ∂ ∂ ∂⎝ ⎠⎝ ⎠G G = ⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+∂ ∂+∂ ∂+∂∂ ∂+∂∂ ∂+∂∂ ∂+∂ ∂+∂ ∂+∂ ∂−= z z y y x x xz zx zy yz yx xy z z y y x x so that the addition of the last two relations yields ( ) ( ) 2 2 222 2 2 2 2 22 2 2ˆ ˆLˆ r p x y z x y zx y zx y z⎛ ⎞ ⎛ ⎞∂ ∂ ∂ ∂ ∂ ∂+ ⋅ = − + + + + + + +⎜ ⎟ ⎜ ⎟∂ ∂ ∂∂ ∂ ∂ ⎝ ⎠⎝ ⎠G G = = ( )pripr ˆˆˆ 22 GG=G ⋅+= (31.10) It follows that the momentum can be expressed as ( ) ( ) ( ) rrr r r L r priprL r p ∂ ∂−⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−=⎥⎦ ⎤⎢⎣ ⎡ ⋅−⋅+= 2 ˆ1ˆˆˆˆˆ1ˆ ==GG=GGG (31.11) where the radial part results from the substitution of Eq. (I.20), Appendix I, which gives 1 1ˆ sinr p i i e e e r r rθ ϕθ θ ϕ ⎛ ⎞∂ ∂ ∂= − ∇ = − + +⎜ ⎟∂ ∂ ∂⎝ ⎠ G G G G= = and hence ˆ ˆr p i r r ∂⋅ = − ∂ G G = Therefore the Hamiltonian (31.1) becomes ( )rUL rrrr r rm H + ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ −∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂−= 222 2 ˆ111 2 ˆ = = and reduces to the equivalent form ( )rUL rr r rrm H +⎥⎦ ⎤⎢⎣ ⎡ −⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ 2 ˆ11 2 ˆ = = (31.12) used for a particle bound by a central force field. 31.2. EIGENVALUE EQUATIONS IN A CENTRAL FIELD Spherically symmetric systems, where the potential energy ( )rU is independent of the direction of ,rG are usually treated using the spherical coordinates ( )ϕθ ,,r related to the Cartesian coordinates as referred to in Appendix I, Eq.(I.19). As the angle θ is defined by the inclination to the z-axis, it is mathematically convenient to take the direction of the component of ( zyx ,, ) L G that is compatible with (in the sense introduced in Chapter 34) to be the z-axis. Therefore we consider the eigenvalue equations of the operators and given by ˆ ˆ, zH L 2Lˆ ( ) ( )ˆ , , , ,H r rθ ϕ ε θ ϕΨ = Ψ ( ) ( )ϕθϕθ ,,,,ˆ rLrL zz Ψ=Ψ (31.13) ( ) ( )ϕθϕθ ,,,,ˆ 22 rLrL Ψ=Ψ The appropriate form of the Hamiltonian Hˆ is conveniently written using the standard expression (I.23) of the Laplace operator in spherical coordinates as 2∇ ( )rU rrr r rrm H +⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂−= 2 sin 1sin sin 11 2 ˆ ϕθθθθθ = (31.14) Using the transformation (I.19) from Cartesian to spherical coordinates and the chain rule for differentiation, we find sin sin sin cosx y z r r x y x y z x y y θ ϕ θ ϕϕ ϕ ϕ ϕ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂= + + = − + = −∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂x so that the operator takes the form zLˆ ϕ∂ ∂=⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ ∂ ∂−∂ ∂= ix y y x i Lz ==ˆ (31.15) Furthermore, from Eqs.(31.12) and (31.14) it follows that is in fact expressed by that part of Hamiltonian which operates on the angular variables 2Lˆ ,θ ϕ only, that is, ⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂−= 2 sin 1sin sin 1ˆ ϕθθθθθ=L (31.16) The energy eigenfunctions in central field, which are simultaneously eigenfunctions of and are subjected, in addition to the usual bound state requirement (30.21), to the single-valueness condition zLˆ 2ˆ ,L ( ) ( )πϕθϕθ 2,,,, +Ψ=Ψ rr since the points ϕ and πϕ 2+ are in fact the same physical point. The solution to the time-independent Schrödinger equation (31.13) in central field can be found by a process of separation of variables, yielding three ordinary differential equations in the variables ,r θ and ,ϕ respectively, to which a physical significance will be assigned. It is common practice to separate the variables in two stages, first substituting the trial solution ( ) ( ) ( ), , ,r R r Yθ ϕΨ = θ ϕ (31.17) into the energy eigenvalue equation. If the Hamiltonian is expanded in spherical coordinates, as given by Eq.(31.14), Eq.(31.13) becomes ( )2 22 2 2 21 1 1 sin2 2 sin sin 21R Y Yr U r r m R r r m Y ε θθ θ θ θ ϕ ⎡ ⎤⎡ ⎤∂ ∂ ∂ ∂ ∂⎛ ⎞ ⎛ ⎞⎡ ⎤+ − = − +⎢ ⎥⎜ ⎟ ⎜ ⎟⎢ ⎥ ⎣ ⎦∂ ∂ ∂ ∂ ∂⎝ ⎠ ⎝ ⎠⎣ ⎦ ⎣ ⎦ = = (31.18) One side is independent of θ and ,ϕ and the other side is independent of r, so they must be separately equal to some constant, which is conveniently chosen to be equal to In this way we obtain a radial equation, given by .2/2 mL ( ) ( ) ( )rRrR mr LrU dr dr dr d rm ε=⎥⎦ ⎤⎢⎣ ⎡ ++⎟⎠ ⎞⎜⎝ ⎛− 2 2 1 2 = (31.19) where was replaced by for obvious reasons, and an angular equation, having the so-called spherical harmonics r∂∂ / drd / ( ),Y θ ϕ as solutions, which reads 2 2 ˆor sin 1sin sin 1 ==⎥⎦ ⎤⎢⎣ ⎡ ∂ ∂+⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂− ϕθθθθθ= (31.20) The form of the spherical harmonics is obtained as a result of the second stage of the separation process, where we assume that ( ) ( ) ( )ϕθϕθ Φ= PY , so that Eq.(31.20) becomes 2 2 1sinsinsin ϕθθθθ θ ∂ Φ∂ Φ−=⎥⎦ ⎤⎢⎣ ⎡ +⎟⎠ ⎞⎜⎝ ⎛ ∂ ∂ ∂ ∂ === LP P (31.21) Since one side is independent of ϕ while the other is independent of ,θ they must each be equal to a constant, which we call We thus obtain the differential equations .2zL d L dϕ Φ 2− = Φ= (31.22) where again ϕ∂∂ / has been replaced by ϕdd / , and 0 sin sin sin 1 22 ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡ ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ −+⎟⎠ ⎞⎜⎝ ⎛ PLL d dP d d z θθθθθ === (31.23) where θ∂∂ / has been replaced by θdd / . Therefore the energy equation (31.13) completely separates, in central field, into a radial equation (31.19) and two angular equations (31.22) and (31.23). The radial equation, giving the energy levels ε of the system, takes the form ( ) ( ) ( ) ( )rRrR mr LrUrR dr d rdr d m ε=⎥⎦ ⎤⎢⎣ ⎡ ++⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛ +− 2 2 2 2 = (31.24) which can be simplified by making the substitution ( ) ( )rrRr =ξ (31.25) and using the identity 2 1r dd d r dr r rdr dr ξ ξ⎛ ⎞+ =⎜ ⎟⎝ ⎠ r Thus Eq.(31.24) reduces to ( ) ( ) ( ) ( )rr mr LrU dr rd m εξξξ =⎥⎦ ⎤⎢⎣ ⎡ ++− 2 22 = (31.26) which is formally identical to the time-independent Schrödinger equation in one dimension (30.17), provided the central force potential ( )rU is replaced by an effective potential energy of the form ( ) ( ) 2 2 2mr LrUrUeff += The repulsive potential, associated with the angular momentum, modifies the central potential as illustrated before, see Figure 3.3, for ( ) rCrU /−= . Further progress with the solution of the radial equation depends on the particular form of the potential function . We shall consider the significant example of one-electron atoms in the next chapter.
THE sub4sub 4sub4RADIAL KINETIC ENERGYsub4sub 4sub4 TERM IN THE SCHRODINGER EQUATION FOR A CENTRAL FORCE