Codex @codex  0

Consider a single electron of mass $m_e$ and a point atomic nucleus of mass $M$ and positive electric charge $Ze$, with $e>0$ the elementary charge. Use nonrelativistic classical mechanics, omit external fields and interactions with other electrons, and writeHere $\mu$ is the Coulomb two-body reduced mass, ${ϵ}_0$ the vacuum permittivity, and $h$ the Planck constant. The separation $r$ and relative speed $v$ describe the two-particle relative coordinate; neither body is silently pinned in place. The electron's own radius about the centre of mass is $Mr/(M+m_e)$.

For the assumed circular orbit, Coulomb's law provides the required centripetal acceleration, and the imposed orbital angular momentum condition isSubstitute $v=n\hslash /(\mu r)$ into the first equation and cancel $r^{-2}$ to obtainThe kinetic energy is $K=\mu v^2/2=\kappa /(2r)$ and the Coulomb potential energy is $U=-\kappa /r$, with zero energy at infinite separation. HenceThis derives the energy sequence from the stated circular model; the modern Schrödinger equation derivation is the distinct series-termination quantization of the Coulomb radial equation. The radius above scales as $n^2/Z$, whereas the binding energy scales as $\mu Z^2/n^2$. The Bohr radius article supplies the corresponding Coulomb potential length scale. The nonrelativistic assumption requires $v_n\ll c$, with $c$ the speed of light.

For an original numerical example choose $Z=1$, $M=1.67262192369\times 10^{-27}\mathrm{kg}$, $m_e=9.1093837015\times 10^{-31}\mathrm{kg}$, $e=1.602176634\times 10^{-19}\mathrm{C}$, $h=6.62607015\times 10^{-34}\mathrm{J}\mathrm{s}$, ${ϵ}_0=8.8541878128\times 10^{-12}\mathrm{F}{\mathrm{m}}^{-1}$, and $c=299792458\mathrm{m}{\mathrm{s}}^{-1}$. Treat these as stipulated constants for the calculation. They give $\mu /m_e=0.99945568$, $r_1=0.0529465\mathrm{nm}$, $E_1=-13.598287\mathrm{eV}$ and $v_1/c\simeq 0.00730$, where $\mathrm{eV}$ denotes an electronvolt. At $n=2$, $r_2=4r_1$ and $E_2=E_1/4$; the ground-state ionization energy is $-E_1$.

Taking $M\to \infty$ gives $\mu \to m_e$. A finite nuclear mass reduces every binding energy by $\mu /m_e$ and increases each separation by $m_e/\mu$. This follows from relative motion, not a fitted correction to individual spectral lines. Replacing the atomic nucleus by a heavier isotope changes the same factor, providing another check of the model.

A Coulomb bound state is a physical normalizable eigenfunction for an attractive potential $V(r)=-\gamma /r$, $\gamma >0$. For the usual regular boundary condition at the origin its energies are the displayed values, $N=1,2,\dots$. The spherically symmetric states have orbital angular momentum zero; series-termination quantization of the Coulomb radial equation constructs them as a polynomial times a decaying exponential. Square integrability on the punctured radial interval alone does not replace the physical condition at the origin.

Put $R=f{e}^{-\gamma r}$. ThenAfter substitution and cancellation of the ${\gamma }^{2}f$ terms, the equation becomesFor the power series $f={\sum }_{n=0}^{\infty }{a}_n{r}^n$, the coefficient of $r^{n-1}$ isThus

If the series does not terminate, then $a_n/a_{n-1}\sim 2\gamma /n$, so $f$ acquires the large-$r$ behaviour $e^{2\gamma r}$ and $R$ grows like $e^{\gamma r}$. Normalizability therefore requires termination, which occurs whenfor some positive integer $N$. HenceThe ground state has $N=1$, so the hydrogen ground-state energy isThis is the series-termination quantization of the Coulomb radial equation.

Let $u=rF$. Multiplication of the radial equation by $-2mr/{\hslash }^2$ gives $u^{\prime \prime }+[A/r-\lambda (\lambda +1)/r^2-{\alpha }^2]u=0$. With $u=e^{-\alpha r}f$ the constant terms cancel, leavingThe indicial equation is $c(c-1)=\lambda (\lambda +1)$, with roots $c=\lambda +1$ and $c=-\lambda$. The usual physical Coulomb radial origin condition selects $c=\lambda +1$: it imposes the regular reduced radial boundary condition $u(0)=0$. Normalizability alone at the origin is insufficient to make this selection when an irregular branch is locally square-integrable, notably at $\lambda =0$.

Write $f=r^{\lambda +1}{\sum }_{j\ge 0}{a}_j{r}^j$. Coefficient comparison givesA nonterminating regular series has the exponentially growing asymptotic branch, suggested by $a_{j+1}/a_j\sim 2\alpha /j$: it supplies an $e^{2\alpha r}$ component, which the prefactor $e^{-\alpha r}$ cannot normalize. Equivalently the regular confluent-hypergeometric solution grows unless its first parameter is a nonpositive integer. Thus a bound state requires the numerator to vanish at some $j=k\ge 0$, givingThis is series-termination quantization of the Coulomb radial equation. It requires $A>0$ for the assumed $\alpha >0$; at $A=0$ there are no such negative-energy Coulomb bound states.

For $k=0$, $F=N{r}^{\lambda }{e}^{-\alpha r}$. The angular integral cancels when computing $⟨1/r⟩$, while the radial volume element contributes $r^{2}dr$. ConsequentlyThe ratio follows by one integration by parts and remains valid for real $\lambda \ge 0$; no integer-only factorial identity is needed. Since now $\alpha =A/[2(\lambda +1)]$,For a classical particle with that negative total energy, nonnegative kinetic energy requires $E-V={\hslash }^2(A/r-{\alpha }^2)/(2m)\ge 0$, hence $\boxed{r\le A/{\alpha }^2}$. This is a necessary radial bound; extra angular momentum can impose a stronger restriction.

At $k=\lambda =0$, $\alpha =A/2$ and the classical boundary is $r=2/\alpha$. The normalized radial probability density is $4{\alpha }^3{r}^2{e}^{-2\alpha r}$. With $z=2\alpha r$, the exterior quantum probability isThe final integration gives $(4^2+2\cdot 4+2)e^{-4}/2$. This is the nonzero quantum tail beyond the classical turning region.

For brevity write $\gamma =mKZ/{\hslash }^2$. With the physical indicial root $c=\ell +1$, substitute $f={\sum }_{n\ge 0}{a}_n{r}^{n+\ell +1}$ into the equation from part (i). The coefficient of $r^{n+\ell }$ isThe indicial equation has already canceled the still lower power $r^{\ell -1}$. Rearranging gives the recurrence relationIts denominator is $(n+1)(n+2\ell +2)>0$ for every $n\ge 0$.

To impose wavefunction normalization at infinity, first suppose that the series never terminates. ThenThis is the exponential-series growth associated with $e^{2\alpha r}$, which would overwhelm the factor $e^{-\alpha r}$ in $\psi$. A direct lower-bound argument makes the conclusion precise. All sufficiently late recurrence ratios are real and positive, so after a common global phase choice the tail coefficients have one sign. For any $0<\eta <\alpha$, there is an index $N$ such thatInduction then bounds the tail below by a positive multiple ofThe finitely many earlier terms form a polynomial and cannot cancel this exponential growth. Thus $|\psi |=r^{-1}{e}^{-\alpha r}|f|$ grows at least as a positive multiple of $r^{\ell }{e}^{(\alpha -\eta )r}$ for large $r$, contradicting wavefunction normalization. This proves that the regular series must terminate.

Since the denominator never vanishes, termination requires the numerator to vanish at some nonnegative integer $n$:Then $a_{n+1}=0$ and every later coefficient is zero. Consequently $f$ is a polynomial of degree $n+\ell +1$, with a factor $r^{\ell +1}$; equivalently,where $P_n$ has degree $n$. This is series-termination quantization of the Coulomb radial equation. If $N_{\mathrm{p}}=n+\ell +1$ denotes the principal quantum number, the associated energy eigenvalues areHere $n$ counts the degree of the remaining radial polynomial; it is not the principal quantum number itself.