\documentclass[11pt,reqno]{amsart} \usepackage{amsmath} \usepackage{amssymb} \usepackage{amsthm} \usepackage{amscd, amsfonts, mathrsfs} \usepackage{cases} \usepackage[dvips]{epsfig} \usepackage{verbatim} \usepackage{epsf} %\usepackage[usenames]{color} %\usepackage{hyperref} \usepackage[bookmarksnumbered,pdfpagelabels=true,plainpages=false,colorlinks=true, linkcolor=black,citecolor=black,urlcolor=blue]{hyperref} \theoremstyle{plain} \newtheorem{theorem}{Theorem}[section] \newtheorem{lemma}[theorem]{Lemma} \newtheorem{prop}[theorem]{Proposition} \newtheorem{coro}[theorem]{Corollary} \newtheorem{question}[theorem]{Question} \newtheorem{conjecture}[theorem]{Conjecture} \newtheorem{assump}[theorem]{Assumption} \theoremstyle{definition} \newtheorem{rema}[theorem]{Remark} \newtheorem{notation}[theorem]{Notation} \newtheorem{defi}[theorem]{Definition} \newtheorem{hyp}[theorem]{Hypothesis} \numberwithin{equation}{section} \newenvironment{dem}{{\noindent \bf Proof: }}{\hfill $\qed$\\ \medskip} \newenvironment{sketch}{{\noindent \bf Sketch proof: }}{\hfill $\qed$\\ \medskip} \newenvironment{ sys_eq }{\ left \ lbrace \ begin { array }{ @ {} l@ {}}}{\ end { array }\ right .} \setlength{\textwidth}{6.6in} \setlength{\textheight}{8.6in} \hoffset=-0.8truein \voffset=-0.1truein \baselineskip=7.0mm %\renewcommand{\baselinestretch}{2.2} \setlength{\baselineskip}{1.09\baselineskip} \gdef\dxd{\displaystyle} \gdef\x{\xi} \newcommand{\cA}{\mathcal A} \newcommand{\cB}{\mathcal B} \newcommand{\cC}{\mathcal C} \newcommand{\cD}{\mathcal D} \newcommand{\cE}{\mathcal E} \newcommand{\cF}{\mathcal F} \newcommand{\cG}{\mathcal G} \newcommand{\cH}{\mathcal H} \newcommand{\cI}{\mathcal I} \newcommand{\cJ}{\mathcal J} \newcommand{\cK}{\mathcal K} \newcommand{\cL}{\mathcal L} \newcommand{\cM}{\mathcal M} \newcommand{\cN}{\mathcal N} \newcommand{\cO}{\mathcal O} \newcommand{\cP}{\mathcal P} \newcommand{\cQ}{\mathcal Q} \newcommand{\cR}{\mathcal R} \newcommand{\cS}{\mathcal S} \newcommand{\cT}{\mathcal T} \newcommand{\cU}{\mathcal U} \newcommand{\cV}{\mathcal V} \newcommand{\cW}{\mathcal W} \newcommand{\cX}{\mathcal X} \newcommand{\cY}{\mathcal Y} \newcommand{\cZ}{\mathcal Z} \newcommand{\cBV}{\cB \cV} \newcommand{\ccB}{\mathscr{B}} \newcommand{\ccE}{\mathscr{E}} \newcommand{\ccD}{\mathscr{D}} \newcommand{\ccL}{\mathscr{L}} \newcommand{\ccW}{\mathscr{W}} \newcommand{\al}{\alpha} \newcommand{\be}{\beta} \newcommand{\ga}{\gamma} \newcommand{\Ga}{\Gamma} \newcommand{\de}{\delta} \newcommand{\ep}{\epsilon} \newcommand{\la}{\lambda} \newcommand{\La}{\Lambda} \newcommand{\si}{\sigma} \newcommand{\Si}{\Sigma} \newcommand{\vfi}{\varphi} \newcommand{\om}{\omega} \newcommand{\Om}{\Omega} \newcommand{\TT}{\mathbb T} \newcommand{\RR}{\mathbb R} \newcommand{\CC}{\mathbb C} \newcommand{\HH}{\mathbb H} \newcommand{\CP}{\mathbb C \mathbb P} \newcommand{\ZZ}{\mathbb Z} \newcommand{\NN}{\mathbb N} \newcommand{\Min}{\mathbb R^{1,3}} \newcommand{\cqd}{\hfill $\qed$\\ \medskip} \newcommand{\rar}{\rightarrow} \newcommand{\imp}{\Rightarrow} \newcommand{\tr}{\operatorname{tr}} \newcommand{\ve}{\varepsilon} \newcommand{\vol}{\operatorname{vol}} \newcommand{\id}{\operatorname{id}} \newcommand{\dive}{\operatorname{div}} \newcommand{\p}{\parallel} %\usepackage{xfrac} %\newcommand{\half}{\sfrac{1}{2}} %\newcommand{\half}{\frac{1}{2}} %The following allows to use footnotes in section headers \usepackage[stable]{footmisc} \title[MATH 234 Spring 14]{Vanderbilt University, MATH 234 Spring 14: The Schr\"odinger equation.} %\date{\today} %\author[Disconzi]{Marcelo M. Disconzi} %\address{Department of Mathematics\\ %Vanderbilt University, Nashville, TN 37240, USA} %\email{marcelo.disconzi@vanderbilt.edu} \urladdr{http://www.disconzi.net/Teaching/MAT234-Spring-14/MAT234-Spring-14.html} %\date{\today} %\pagestyle{empty} \begin{document} %\begin{abstract} %Abstract. %\end{abstract} \maketitle \tableofcontents \section{Introduction.} Our goal is to investigate solutions to the \emph{Schr\"odinger equation}, \begin{gather} i \hbar \frac{\partial \Psi}{\partial t} = -\frac{\hbar^2}{2\mu} \Delta \Psi + V \Psi, \label{Sch_eq} \end{gather} where $i$ is the imaginary number $i^2 = -1$; $\hbar = 1.51 \times 10^{-27}\, erg \, s$ is Planck's constant; $\mu$ is a positive constant called the mass; $V = V(t,x) : \RR\times \RR^3 \rar \RR$ is called the potential function; and the unknown is the complex-valued function $\Psi = \Psi(t,x): \RR \times \RR^3 \rar \CC$ called the wave-function. The variables $t$ and $x$ represent, respectively, the time and space variables. The Schr\"odinger equation describes the dynamics of a particle of mass $\mu$ interacting with a potential $V$, according to the laws of Quantum Mechanics. The physical interpretation of $\Psi$ is as follows. If $U \subseteq \RR^3$, then \begin{gather} \int_U \left| \Psi(t,x) \right|^2 \, dx \nonumber \end{gather} represents the probability of finding the particle in the region $U$ at a time $t$. In particular, one must have \begin{gather} \int_{\RR^3} \left| \Psi(t,x) \right|^2 \, dx = 1. \label{psi_normalized} \end{gather} Notice that, upon multiplying $\Psi$ by a suitable constant, condition (\ref{psi_normalized}) can always be fulfilled as long as \begin{gather} \int_{\RR^3} \left| \Psi(t,x) \right|^2 \, dx < \infty. \label{normalizable} \end{gather} Our treatment will be based on \cite{FY, T, W}, to which the student is referred for more details. \section{Separation of variables for a time-independent potential.} We shall assume that $V$ does not depend on time, i.e., $V(t,x) = V(x)$. We will have to divide several expressions by $\Psi$. In order to make this sensible, it will be assumed that $\Psi$ does not vanish (or, at least, does not vanish on an open set). Look for solutions of the form \begin{gather} \Psi(t,x) = T(t) \psi(x), \label{separation_t_x} \end{gather} Plugging (\ref{separation_t_x}) into (\ref{Sch_eq}) gives \begin{gather} i \hbar \frac{T^\prime}{T} = -\frac{\hbar^2}{2\mu} \frac{1}{\psi} \Delta \psi + V, \nonumber \end{gather} The left-hand side depends only on $t$, whereas the right-hand side depends only on $x$. Thus, both sides have to be equal to a constant, which we denote be $E$. Therefore \begin{gather} i\hbar T^\prime = E T, \label{T_eq} \end{gather} and \begin{gather} -\frac{\hbar^2}{2\mu} \Delta \psi + V \psi = E \psi. \label{time_independent} \end{gather} Equation (\ref{T_eq}) is easily solved. Its solution is \begin{gather} T(t) = e^{-\frac{i E}{\hbar} t}, \label{T_sol} \end{gather} where we ignored an arbitrary constant of integration (such constants will be neglected throughout, as an overall constant of integration can be fixed at the very end via condition (\ref{psi_normalized}) ). \subsection{The time-independent Schr\"odinger equation.} We now focus on (\ref{time_independent}), known as the \emph{time-independent Schr\"odinger equation.} To solve it, we assume further that $V$ is \emph{radially symmetric}, i.e., that $V(x) = V(\sqrt{x_1^2 + x_2^2 + x_3^2})$ or, in spherical coordinates, that $V = V(r)$. This assumption suffices to treat many physical systems of interest. Recall the expression for the Laplacian in spherical coordinates, \begin{gather} \Delta = \partial_r^2 + \frac{2}{r} \partial_r + \frac{1}{r^2} \Delta_{S^2}, \label{Laplacian_spherical} \end{gather} where \begin{gather} \Delta_{S^2} = \partial_\phi^2 + \frac{\cos \phi}{\sin \phi} \partial_\phi +\frac{1}{\sin^2 \phi} \partial^2_\theta \label{Laplacian_unit_sphere} \end{gather} is the \emph{Laplacian on the unit sphere}, $r \in [0,\infty)$, $\phi \in [0,\pi)$, and $\theta \in [0,2\pi)$. From now on we shall work in spherical coordinates. We suppose that \begin{gather} \psi(r,\phi, \theta) = R(r) Y(\phi,\theta). \label{separation_radial_angular} \end{gather} Plugging (\ref{separation_radial_angular}) into (\ref{time_independent}), and using (\ref{Laplacian_spherical}), \begin{gather} -\frac{\hbar^2}{2\mu} \frac{r^2}{R}\left( R^{\prime\prime} + \frac{2}{r} R^\prime \right) + (V - E) r^2 = \frac{\hbar^2}{2\mu} \frac{1}{Y} \Delta_{S^2} Y. \nonumber \end{gather} Since the left-hand side depends only on $r$ and the right-hand side only on $(\phi,\theta)$, both sides must be equal to a constant, which we denote by $-a$. Thus, \begin{gather} -\frac{\hbar^2}{2\mu}\left( R^{\prime\prime} + \frac{2}{r} R^\prime \right) + \left(V + \frac{a}{r^2} \right) R = ER, \label{radial_eq} \end{gather} and \begin{gather} \frac{\hbar^2}{2\mu} \Delta_{S^2} Y = - a Y. \label{angular_eq} \end{gather} \subsection{The angular equation.\label{angular_section}} We first investigate (\ref{angular_eq}), which, in light of (\ref{Laplacian_unit_sphere}), becomes \begin{gather} \partial^2_\phi Y + \frac{\cos \phi}{\sin \phi} \partial_\phi Y + \frac{1}{\sin^2\phi} \partial^2_\theta Y = -\frac{2a\mu}{\hbar^2} Y. \nonumber \end{gather} Supposing \begin{gather} Y(\phi, \theta) = \Phi(\phi) \Theta(\theta), \label{Y_Phi_Theta} \end{gather} one gets \begin{gather} - \frac{\Theta^{\prime\prime}}{\Theta} = \frac{\sin^2 \phi}{\Phi} \Phi^{\prime\prime} + \frac{\sin\phi \cos\phi}{\Phi} \Phi^\prime + \frac{ 2 a \mu \sin^2 \phi}{\hbar^2}. \label{eqs_angular} \end{gather} Once more, both sides ought to be equal to a constant, which we denote by $b$. One equation becomes \begin{gather} \Theta^{\prime\prime} = - b \Theta. \label{Theta_eq} \end{gather} To solve (\ref{Theta_eq}), we need to analyze the cases $b>0$, $b=0$, and $b< 0$. Notice the following boundary condition: the points with coordinates $\theta$ and $\theta + 2\pi$ must be identified as they correspond to the same point in $\RR^3$. Thus, \begin{gather} \Theta ( \theta + 2\pi) = \Theta(\theta). \label{Theta_bry} \end{gather} We immediately see that the case $b < 0$ does not yield a solution satisfying (\ref{Theta_bry}); $b=0$ and (\ref{Theta_bry}) give $\Theta=$constant; and $b > 0$ along with (\ref{Theta_bry}) give that $\Theta$ is a linear combination of $\cos( \sqrt{ b } \theta)$ and $\sin( \sqrt{ b } \theta)$. All this cases can be summarized by setting \begin{gather} b = m^2, \, m \in \ZZ, \label{b_m} \end{gather} and writing \begin{gather} \Theta( \theta ) = e^{im\theta}. \label{Theta_sol} \end{gather} Next, we move to the $\Phi$-equation. From (\ref{eqs_angular}) and (\ref{b_m}), one has \begin{gather} \frac{\sin \phi}{\Phi} \frac{d}{d \phi} \left( \sin\phi \frac{ d\Phi}{d\phi} \right) - m^2 = -\la \sin^2\phi, \label{Phi_eq} \end{gather} where \begin{gather} \la = \frac{ 2 \mu }{\hbar^2 } a, \label{la_a} \end{gather} and we used the product rule to rewrite the terms involving derivatives. In order to solve (\ref{Phi_eq}), let us make the following change of variables, \begin{gather} x = \cos \phi, \, 0 \leq \phi \leq \pi. \nonumber \end{gather} (strictly speaking, $\phi \in [0,\pi)$, but it is convenient to include $\phi = \pi$. Notice that this change of variables is well-defined since $\cos$ is one-to-one for $0 \leq \phi \leq \pi$). The chain rule gives \begin{gather} \sin \phi \frac{d}{d \phi} = \sin \phi \frac{dx}{d\phi} \frac{d}{dx} = - \sin^2\phi \frac{d}{dx} = (\cos^2\phi - 1 ) \frac{d}{dx} = (x^2 - 1 ) \frac{d}{dx}, \nonumber \end{gather} so that (\ref{Phi_eq}) becomes \begin{gather} \frac{d}{dx} \left( (1-x^2) \frac{d\Phi}{dx} \right) + \left( \la - \frac{m^2}{1-x^2}\right) \Phi = 0. \label{Phi_x_eq} \end{gather} To solve (\ref{Phi_x_eq}), we seek for a solution of the form \begin{gather} \Phi(x) = (1 - x^2)^\frac{|m|}{2} \frac{ d^{|m|} P(x) }{dx^{|m|}}, \label{Phi_P} \end{gather} where $P$ solves \begin{gather} (1-x^2) \frac{d^2 P}{dx^2} - 2 x \frac{d P }{dx} + \la P = 0. \label{Legengre} \end{gather} To see that this works, differentiate (\ref{Legengre}) $|m|$ times, obtaining \begin{gather} (1-x^2) \frac{d^{|m|+2} P}{dx^{|m|+2}} - 2 (|m|+1) x \frac{d^{|m|+1} P }{dx^{|m|+1}} + \left(\la - |m|(|m|+1) \right) \frac{d^{|m|} P}{dx^{|m|}} = 0. \label{Legengre_m} \end{gather} Students are encouraged to verify (\ref{Legengre_m}) (compute the first few derivatives to see that a pattern as (\ref{Legengre_m}) emerges). On the other hand, let $\widetilde{\Phi}$ be defined by \begin{gather} \Phi(x) = (1-x^2)^\frac{|m|}{2} \widetilde{\Phi}(x) \label{Phi_Phi_tilde} \end{gather} and plug this into (\ref{Phi_x_eq}). Computing the derivative terms, \begin{align} \begin{split} \frac{d}{dx} & \left( (1-x^2) \frac{d}{dx} \left( (1-x^2)^\frac{|m|}{2} \widetilde{\Phi} \right) \right) = \frac{d}{dx} \left( \frac{|m|}{2}(-2x)(1-x^2)^\frac{|m|}{2} \widetilde{\Phi} + (1-x^2)^{\frac{|m|}{2} + 1} \frac{d\widetilde{\Phi}}{dx} \right) \\ & = (1-x^2)^{\frac{|m|}{2}+1} \frac{d^2 \widetilde{\Phi}}{dx^2} + (1-x^2)^\frac{|m|}{2} \frac{d\Phi}{dx} \left( \left(\frac{|m|}{2} + 1 \right)(-2x) + \frac{|m|}{2}(-2x)\right)\\ &+ \frac{|m|}{2}\left( (-2x)\frac{|m|}{2} (1-x^2)^{\frac{|m|}{2}-1}(-2x) - 2(1-x^2)^\frac{|m|}{2}\right) \widetilde{\Phi} \\ & = (1-x^2)^{\frac{|m|}{2}+1} \frac{d^2 \widetilde{\Phi}}{dx^2} - 2x(1-x^2)^\frac{|m|}{2} \left(|m|+1\right) \frac{d\widetilde{\Phi}}{dx} + \frac{|m|}{2} (1-x^2)^\frac{|m|}{2}\left( \frac{2|m|x^2}{1-x^2} - 2 \right) \widetilde{\Phi} \\ & = (1-x^2)^\frac{|m|}{2} \left( (1-x^2) \frac{d^2 \widetilde{\Phi}}{dx^2} - 2x \left(|m|+1\right) \frac{d\widetilde{\Phi}}{dx} + |m|\left( \frac{|m|x^2}{1-x^2} - 1 \right) \widetilde{\Phi} \right). \end{split} \nonumber \end{align} By (\ref{Phi_x_eq}), this has to equal \begin{gather} - \left( \la - \frac{m^2}{1-x^2}\right) \Phi = \left( \la - \frac{m^2}{1-x^2}\right) (1-x^2)^\frac{|m|}{2} \widetilde{\Phi}(x), \nonumber \end{gather} what gives, after canceling $(1-x^2)^\frac{|m|}{2}$, \begin{gather} (1-x^2) \frac{d^2 \widetilde{\Phi}}{dx^2} - 2x \left(|m|+1\right) \frac{d\widetilde{\Phi}}{dx} + \la \widetilde{\Phi} + \left( \frac{|m|^2 x^2}{1-x^2} - |m| - \frac{|m|^2}{1-x^2} \right) \widetilde{\Phi} = 0. \nonumber \end{gather} But \begin{gather} \frac{|m|^2 x^2}{1-x^2} - |m| - \frac{|m|^2}{1-x^2} = \frac{|m|(|m|+1)x^2 - |m|(|m|+1)}{1-x^2} = -|m|(|m|+1), \nonumber \end{gather} and therefore \begin{gather} (1-x^2) \frac{d^2 \widetilde{\Phi}}{dx^2} - 2x \left(|m|+1\right) \frac{d\widetilde{\Phi}}{dx} + \left( \la - |m|(|m|+1) \right) \widetilde{\Phi} = 0. \label{Phi_tilde_eq} \end{gather} Comparing (\ref{Phi_tilde_eq}) with (\ref{Legengre_m}), we see that if $P$ solves (\ref{Legengre}), then (\ref{Phi_Phi_tilde}) solves (\ref{Phi_x_eq}), as claimed. Therefore, it suffices to solve (\ref{Legengre}). We seek a power series solution of the form \begin{gather} P(x) = \sum_{k=0}^\infty a_k x^k. \label{P_ps} \end{gather} Plugging (\ref{P_ps}) into (\ref{Legengre}) gives \begin{gather} (1-x^2)\sum_{k=0}^\infty k(k-1)a_k x^{k-2} - 2x \sum_{k=0}^\infty k a_k x^{k-1} + \la \sum_{k=0}^\infty a_k x^k = 0, \nonumber \end{gather} or yet, after rearranging some terms, \begin{gather} \sum_{k=0}^\infty \left( (k+2)(k+1) a_{k+2} - (k(k+1) - \la ) a_k \right) x^k = 0, \nonumber \end{gather} which implies the following recurrence relation, \begin{gather} a_{k+2} = \frac{ k(k+1) - \la }{ (k+1)(k+2) } a_k, \, k = 0, 1, 2, \dots. \label{recurrence} \end{gather} (\ref{recurrence}) determines all coefficients $a_k$ except $a_0$ and $a_1$, which remain arbitrary (this is consistent with the fact that we are solving a second order ODE). Furthermore, $a_0$ determines all even coefficients, giving rise to an even power series, while while $a_1$ determines all odd coefficients, giving rise to an odd power series. This two power series, even and odd, are two linearly independent solutions of (\ref{Legengre}). Next, we investigate the convergence of (\ref{P_ps}). Since it suffices to investigate the convergence of the even and odd expansions separately, as these are two linearly independent solutions, the ratio between two consecutive terms in the expansion is obtained from (\ref{recurrence}), yielding \begin{gather} \lim_{k\rar \infty} \left| \frac{ a_{k+2}x^{x+2} }{ a_k x^k } \right| = |x|^2, \nonumber \end{gather} and thus (\ref{P_ps}) converges for $|x| < 1$ by the ratio test. We need to investigate the case $|x| = 1$, in particular $x=1$ (which corresponds to $\phi = 0$). Plugging $x=1$ into (\ref{P_ps}) gives \begin{gather} P(1) = \sum_{k=0}^\infty a_k. \label{P_1} \end{gather} From (\ref{recurrence}) we have \begin{align} \begin{split} a_{k+2} = \frac{k^2 + O(k) }{k^2 + O(k) }a_k = \frac{k^2 + O(k) }{k^2 + O(k) }\frac{k^2 + O(k) }{k^2 + O(k) }a_{k-2} = \cdots = \begin{cases} \frac{k^{k+2} + O(k^{k+1} ) }{k^{k+2} + O(k^{k+1}) } a_0, & k \, \text{even}, \\ \frac{k^{k+1} + O(k^{k} ) }{k^{k+1} + O(k^{k}) } a_1, & k \, \text{odd}. \end{cases} \end{split} \nonumber \end{align} It follows that \begin{gather} \lim_{k \rar \infty} a_k \neq 0, \nonumber \end{gather} and therefore (\ref{P_1}) diverges by the divergence test, \emph{unless} (\ref{P_ps}) is in fact a \emph{finite sum}; i.e., \emph{unless} $a_k = 0$ \emph{for all $k$ greater than a certain $\ell$}. Hence, we must have, form some non-negative integer $\ell$, \begin{gather} a_{\ell+2} = 0 = \frac{ \ell(\ell+1) - \la }{ (\ell+1)(\ell+2) } a_\ell, \nonumber \end{gather} which implies \begin{gather} \la = \ell(\ell + 1), \label{la_ell} \end{gather} provided that $a_\ell \neq 0$. (\ref{la_ell}) determines $\la$, and hence the separation constant $a$ in view of (\ref{la_a}). The conclusion is that there is a family $\{ P_\ell\}$ of solutions to (\ref{Legengre}) parametrized by $\ell = 0, 1, 2, \dots$. After conveniently choosing $a_0$ and $a_1$ to obtain integer coefficients, the first few $P$'s are \begin{gather} P_0(x) = 1,\, P_1(x) = x, \, P_2(x) = 1 - 3x^2, \, P_3(x) = 3x - 5x^3. \nonumber \end{gather} Since $P_\ell$ is a polynomial of degree $\ell$, from (\ref{Phi_P}) it follows that $\Phi = 0$ for $|m| > \ell$. Thus, the values of $m$ are restricted to $|m| \leq \ell$, i.e., the allowed $m$-values depend on $\ell$ and satisfy \begin{gather} m \in \big \{ -\ell, -\ell + 1, \dots, -1, 0, 1, \dots, \ell -1, \ell \big \}. \label{m_restriction} \end{gather} We write $m = m_\ell$ when we want to stress this dependence of $m$ on $\ell$. One obtains a family of solutions $\{ \Phi_{\ell m_\ell} \}$ to (\ref{Phi_x_eq}) parametrized by $\ell$ and $m_\ell$, where $\ell=0,1,2,\dots$ and $m_\ell$ satisfies (\ref{m_restriction}). The first few $\Phi$'s are \begin{align} \begin{split} & \Phi_{00}(x) = 1, \\ & \Phi_{10}(x) = x, \, \Phi_{1,\pm 1}(x) = (1-x^2)^\frac{1}{2}, \\ & \Phi_{20}(x) = 1 - 3x^2, \, \Phi_{2 \pm 1}(x) = (1-x^2)^\frac{1}{2} x, \, \Phi_{2\pm 2}(x) = 1-x^2, \\ & \Phi_{30}(x) = 3x - 5x^3, \, \Phi_{3 \pm 1}(x) = (1-x^2)^\frac{3}{2}(1-5x^2), \, \Phi_{3\pm 2}(x) = (1-x^2)x, \, \Phi_{3 \pm 3} (1-x^2)^\frac{3}{2}. \end{split} \nonumber \end{align} Finally, it is necessary to rewrite our solutions in terms of the $\phi$ variable. Denoting $F_{\ell m_\ell} = \widetilde{\Phi}_{\ell m_\ell}$, and using $1-x^2 = \sin^2 \phi$, \begin{gather} \Phi_{\ell m_\ell}(\phi) = \sin^{|m_\ell|} F_{\ell m_\ell}(\cos\phi), \, \ell = 0, 1, 2, \dots, \, |m_\ell| \leq \ell. \label{Phi_sol} \end{gather} Combining (\ref{Y_Phi_Theta}), (\ref{Theta_sol}), and (\ref{Phi_sol}) gives \begin{gather} Y_{\ell m_\ell}(\phi, \theta) = e^{i m_\ell \theta} \sin^{|m_\ell|} \phi F_{\ell m_\ell}(\cos\phi), \, \ell = 0, 1, 2, \dots, \, |m_\ell| \leq \ell. \label{spherical_harmonics} \end{gather} Notice that, in view of (\ref{angular_eq}) and (\ref{la_a}), $Y_{\ell, m_\ell}$ solves \begin{gather} \Delta_{S^2} Y_{\ell m_\ell} = - \ell(\ell + 1) Y_{\ell m_\ell}. \nonumber \end{gather} We finish this section with some terminology. Equation (\ref{Legengre}) is known as Legendre equation, and its solutions $P_\ell$ are known as Legendre polynomials. The functions $F_{\ell m_\ell}$ are known as associated Legendre functions. The functions $Y_{\ell m_\ell}$ are called spherical harmonics. Legendre functions and spherical harmonics have many important applications in Physics. The interested reader is referred to \cite{B} for details. \subsection{The radial equation.} We now turn our attention to equation (\ref{radial_eq}). Using (\ref{la_a}) and (\ref{la_ell}), equation (\ref{radial_eq}) can be written as \begin{gather} \frac{1}{r^2} \frac{d}{dr}\left( r^2 \frac{d R}{dr}\right) + \frac{2\mu}{\hbar^2} \left( E - V(r) \right) R = \ell (\ell + 1) \frac{R}{r^2}. \label{radial_eq_ell} \end{gather} It is important to stress that the results of section \ref{angular_section} are general, i.e., they apply to separation of variables to any radially symmetric potential $V=V(r)$. To solve (\ref{radial_eq_ell}), on the other hand, we need to specify the function $V(r)$. We shall assume that $V$ is the potential describing the electromagnetic interaction of an electron with a proton. This covers the important case when one is solving the Schr\"odinger equation describing the evolution of an electron on a hydrogen atom. In this situation, $V$ takes the form \begin{gather} V(r) = - \frac{ Z e^2}{4\pi \varepsilon_0 r}, \label{Coulomb_potential} \end{gather} where $Z$ is the nuclear charge (for example, $Z=1$ for the hydrogen and $Z=2$ for an ionized helium atom), $-e$ is the electron charge, where $e=1.6 \times 10^{-19} \,C$, and $\varepsilon_0$ is the vacuum permitivity whose values is $\varepsilon_0 = 8.85 \times 10^{-12} \, F/m$ (farads per meters). In order to investigate solutions to (\ref{radial_eq_ell}) with $V$ given by (\ref{Coulomb_potential}), one needs more information about the separation constant $E$. We claim that $E$ must be real and negative. To see this, multiply equation (\ref{radial_eq_ell}) by $r^2 R^*$, where $R^*$ is the complex conjugate of $R$, and integrate from $0$ to $\infty$: \begin{gather} \int_0^\infty R^* \frac{d}{dr}\left( r^2 \frac{d R}{dr}\right) \, dr - \frac{2\mu}{\hbar^2} \int_0^\infty V |R|^2 r^2 \,dr -\ell (\ell + 1) \int_0^\infty |R|^2\, dr = -\frac{2\mu}{\hbar^2} E \int_0^\infty |R|^2 r^2 \, dr, \label{E_real_multiply} \end{gather} where we used that $|R|^2 = R^* R$. Integrating by parts the first term, \begin{align} \begin{split} \int_0^\infty R^* \frac{d}{dr}\left( r^2 \frac{d R}{dr}\right) \, dr = -\int_0^\infty \frac{d R^*}{dr} \frac{d R}{dr} r^2 \, dr + R^* r^2 \frac{d R}{dr}\Big|_0^\infty = -\int_0^\infty \frac{d R^*}{dr} \frac{d R}{dr} r^2 \, dr \end{split} \label{E_real_parts} \end{align} where it has been assumed that $R^*$ and $\frac{dR}{dr}$ vanish sufficiently fast at $\infty$. Writing \begin{gather} R = R_R + i R_C, \nonumber \end{gather} where $R_R$ and $R_C$ are real-valued, it comes \begin{align} \begin{split} \frac{dR^*}{dr} \frac{dR}{dr} = (\frac{d R_R}{dr} - i \frac{d R_C}{dr})(\frac{d R_R}{dr} + i \frac{d R_C}{dr}) = \left( \frac{d R_R}{dr} \right)^2 + \left( \frac{d R_C}{dr}\right)^2, \end{split} \nonumber \end{align} and we conclude that $\frac{dR^*}{dr} \frac{dR}{dr}$ is real-valued. But from (\ref{E_real_multiply}) and (\ref{E_real_parts}) we have \begin{gather} E = \frac{\int_0^\infty \frac{d R^*}{dr} \frac{d R}{dr} r^2 \, dr +\frac{2\mu}{\hbar^2} \int_0^\infty V |R|^2 r^2 \,dr +\ell (\ell + 1) \int_0^\infty |R|^2\, dr }{ \frac{2\mu}{\hbar^2} \int_0^\infty |R|^2 r^2 \, dr}. \label{E_int_R} \end{gather} Therefore, since all terms on the right-hand side are real, we conclude that $E$ is real as well. Students should notice that (\ref{E_int_R}) gives an explicit expression for $E$ in terms of (the integral of) $R$ and other data of the problem (although we shall derive a much more explicit expression for $E$, see below). Now that we know that $E$ is real, let us show that it is negative\footnote{It is possible to obtain $E<0$ by a more delicate analysis of (\ref{E_int_R}), but here we employ a simpler argument.}. Let us investigate the behavior of (\ref{radial_eq_ell}) for large values of $r$, i.e., $r \gg 1$. Then we can neglect the terms that contain $\frac{1}{r}$ and (\ref{radial_eq_ell}) gives, after expanding the terms in $\frac{d}{dr}$, \begin{gather} \frac{d^2R}{dr^2} \approx -\frac{2\mu E}{\hbar^2} R. \label{aprox_1} \end{gather} But for $r \gg 1$ we also have the approximation \begin{gather} r \frac{d^2R}{dr^2} + \frac{d R}{dr} \approx r\frac{d^2R}{dr^2}, \nonumber \end{gather} so that \begin{gather} \frac{d^2 (rR) }{dr^2} = r \frac{d^2R}{dr^2} + 2\frac{d R}{dr} \approx r\frac{d^2R}{dr^2}. \label{aprox_2} \end{gather} Hence, multiplying (\ref{aprox_1}) by $r$ and using (\ref{aprox_2}), \begin{gather} \frac{d^2( rR) }{dr^2} \approx -\frac{2\mu E}{\hbar^2}(rR). \nonumber \end{gather} This approximate equation can be easily solved, producing \begin{gather} rR \approx e^{\pm \frac{ \sqrt{-2\mu E }}{\hbar} r}. \nonumber \end{gather} If $E \geq 0$, then $R$ is a complex function which satisfies \begin{gather} | r R | \approx 1 \text{ for } r \gg 1. \nonumber \end{gather} Then the integral \begin{gather} \int_{\RR^3} \left| \Psi(t,x) \right|^2 \, dx = \left( \int_0^{2\pi} \int_0^\pi |Y(\phi,\theta)|^2 \sin \phi \, d\phi d\theta \right) \left( \int_0^\infty |R(r)|^2 r^2 \, dr \right) \nonumber \end{gather} diverges since $|R(r)|^2 r^2 \approx 1$ for large $r$. Consequently, condition (\ref{normalizable}) fails, and this does not produce a physically sensible solution. In light of the above arguments, we assume, once and for all, that $E< 0$. In this case, we can define the real constants \begin{gather} \be^2 = - \frac{2\mu E}{\hbar^2}, \label{be_def} \end{gather} and \begin{gather} \ga = \frac{\mu Z e^2}{4 \pi \varepsilon_0 \hbar^2 \beta}, \label{ga_def} \end{gather} and make the real change of variables \begin{gather} \varrho = 2\beta r. \nonumber \end{gather} With these definitions, equation (\ref{radial_eq_ell}), with $V$ given by (\ref{Coulomb_potential}), becomes \begin{gather} \frac{1}{\varrho^2} \frac{d}{d\varrho}\left(\varrho^2 \frac{d R}{d \varrho }\right) + \left( -\frac{1}{4} - \frac{\ell(\ell+1)}{\varrho^2} + \frac{\ga}{\varrho} \right) R = 0. \label{radial_rho} \end{gather} Equation (\ref{radial_rho}) will be solved using a power series expansion, but direct application of the method does not work. To see this, try plugging \begin{gather} R(\varrho) = \sum_{k=0}^\infty a_k \varrho^k \nonumber \end{gather} into (\ref{radial_rho}), obtaining \begin{gather} \sum_{k=0}^\infty k(k+1) a_k \varrho^{k-2} -\frac{1}{4} \sum_{k=0}^\infty a_k \varrho^k - \ell(\ell+1)\sum_{k=0}^\infty a_k \varrho^{k-2} + \ga \sum_{k=0}^\infty a_k \varrho^{k-1} = 0. \nonumber \end{gather} This can be rewritten as \begin{gather} -\ell(\ell+1) a_0\varrho^{-2} + \left( (2-\ell(\ell+1))a_1 + \ga a_0 \right) \varrho^{-1} \nonumber \\ + \sum_{k=0}^\infty \left( \left( (k+3)(k+2) - \ell(\ell + 1) \right) a_{k+2} + \ga a_{k+1} - \frac{1}{4} a_k \right) \varrho^k = 0. \label{attempt_R} \end{gather} Vanishing of each term order by order implies that $a_0 = 0$, then $a_1 = 0$, and subsequently $a_k = 0$ for any $k$, so $R = 0$. We need, therefore, to try a different approach. We shall focus on the behavior of (\ref{radial_rho}) when $\varrho \gg 1$, in which case the equation simplifies to \begin{gather} \frac{1}{\varrho^2} \frac{d}{d\varrho}\left(\varrho^2 \frac{d R}{d \varrho }\right) \approx \frac{R}{4}. \label{approx_rho} \end{gather} This (approximate) equation can be solved as follows. Look for a solution of the form $e^{A \varrho}$. Plugging into the equation we find $A = -\frac{1}{2}$, i.e., $e^{-\frac{\varrho}{2}}$ is a (approximate) solution of (\ref{approx_rho}). This suggests\footnote{The reader may remember that when one solves second order ODEs with constant coefficients, sometimes we have to multiply a solution by a suitable power of the variable in order to produce a particular solution or a second linearly independent solution. What it is being done here resembles that: we have some information about solutions, i.e., that $e^{-\frac{\varrho}{2}}$ solves the equation (in an approximate sense) for large values of $\varrho$. Thus, we try multiplying by $e^{-\frac{\varrho}{2}}$ to construct the full, exact solution. } looking for solutions of (\ref{radial_rho}) in the form \begin{gather} R(\varrho) = e^{-\frac{\varrho}{2}} G(\varrho). \label{an_G} \end{gather} Plugging (\ref{an_G}) into (\ref{radial_rho}) gives an equation for $G$, \begin{gather} \frac{ d^2 G}{d\varrho^2} + \left(\frac{2}{\varrho} - 1\right) \frac{\partial G}{\partial \varrho} + \left( \frac{\ga-1}{\varrho} - \frac{\ell(\ell+1)}{\varrho^2}\right) G = 0. \label{G_eq} \end{gather} We seek a solution of the form \begin{gather} G(\varrho) = \varrho^s \sum_{k=0}^\infty a_k \varrho^k = \sum_{k=0}^\infty a_k \varrho^{k+s}, \label{G_ps} \end{gather} where $s$ is to be determined. The term $\varrho^s$ has been included due to the $\frac{1}{\varrho}$ terms in the equation, as these may lead to singular terms that do not fit into a general recurrence relation, as it occurred in (\ref{attempt_R}). Notice that the traditional procedure is included in this approach by simply setting $s=0$. Plugging (\ref{G_ps}) into (\ref{G_eq}) gives, after some algebra, \begin{align} \begin{split} & \left(s(s+1) - \ell(\ell+1) \right) a_0 \varrho^{s-2} \\ + \sum_{k=0}^\infty \Big ( ( (s+k+1)&(s+k+2) - \ell(\ell+1) ) a_{k+1} - (s+k+1 -\ga )a_k \Big ) \varrho^{s + k-1} = 0. \end{split} \label{series_G_eq} \end{align} The vanishing of the term in $\varrho^{s-2}$ requires \begin{gather} s(s+1) - \ell(\ell+1) = 0, \nonumber \end{gather} which has roots $s = \ell$ and $s=-(\ell+1)$. This latter root is rejected on the basis that it does not yield a finite solution when $\varrho \rar 0^+$, i.e., $G(\varrho)$ blows up at the origin when $s=-(\ell+1)$ (recall that $\ell$ is non-negative). Using $s=\ell$, one finds from (\ref{series_G_eq}) the following recurrence relation, \begin{gather} a_{k+1} = \frac{ k + \ell + 1 - \ga }{(k+\ell + 1)(k + \ell + 2) - \ell(\ell + 1) } a_k. \label{recurrence_G} \end{gather} From (\ref{recurrence_G}) and the ratio test, we see at once that (\ref{G_ps}), with $s=\ell$, converges for all values of $\varrho$. In order for (\ref{G_ps}) to be an acceptable solution, we also must verify (\ref{normalizable}). From (\ref{recurrence_G}), it follows that \begin{gather} a_{k+1} = \frac{ k + \cdots }{k^2 + \cdots} a_k = \frac{ 1 + \cdots }{k + \cdots} a_k, \nonumber \end{gather} and \begin{gather} a_{k} = \frac{ k - 1 + \cdots }{(k-1)^2 + \cdots} a_{k-1} = \frac{ 1 + \cdots }{(k-1) + \cdots} a_{k-1}, \nonumber \end{gather} so that \begin{align} \begin{split} a_{k+1} & = \frac{ 1 + \cdots }{k + \cdots} a_k = \frac{ 1 + \cdots }{k + \cdots} \frac{ 1 + \cdots }{(k-1) + \cdots} a_{k-1}\\ & = \frac{ 1 + \cdots }{k(k-1) + \cdots} a_{k-1}. \end{split} \nonumber \end{align} Continuing this way, \begin{gather} a_{k+1} = \frac{1+\cdots}{k(k-1)(k-2) \cdots (k-j) + \cdots }a_{k-j}. \nonumber \end{gather} Remembering that \begin{gather} e^\varrho = \sum_{k=0}^\infty \frac{1}{k!} \varrho^k, \nonumber \end{gather} we see that $G(\varrho)$ is asymptotic to $\varrho^s e^\varrho$, i.e., its series expansion behaves very much like the series of $\varrho^\ell e^\varrho$ (recall that $s=\ell)$: \begin{gather} G(\varrho) \sim \varrho^\ell e^\varrho, \nonumber \end{gather} which implies, upon recalling (\ref{an_G}), \begin{gather} R(\varrho) = e^{-\frac{\varrho}{2}} G(\varrho) \sim e^{-\frac{\varrho}{2}} \varrho^\ell e^\varrho = \varrho^\ell e^{\frac{\varrho}{2}}, \nonumber \end{gather} which diverges when $\varrho \rar \infty$. As a consequence, (\ref{normalizable}) is not satisfied. This will be the case, \emph{unless} the series (\ref{G_ps}) terminates, i.e., unless $a_k = 0$ for all $k$ greater than a certain $n$. From (\ref{recurrence_G}), this means \begin{gather} k + \ell + 1 - \ga = 0, \nonumber \end{gather} i.e., \begin{gather} \ga = k + \ell + 1. \nonumber \end{gather} In particular, $\ga$ has to be an integer, \begin{gather} \ga = n, \, n = \ell+1, \ell +2, \dots. \nonumber \end{gather} With this, the series terminates at the $(n - (\ell+1))^{\text{th}}$ term, and $G$ is a polynomial of degree $n-1$. Recalling (\ref{be_def}) and (\ref{ga_def}), we have found the possible values for the separation constant $E = E_n$, namely, \begin{gather} E_n = -\frac{ \mu Z^2 e^4}{2(4\pi\varepsilon_0)^2 \hbar^2 n^2}, \, n=1, 2, 3, \dots \label{E_n} \end{gather} We write $R_{n\ell}$ to indicate that $R$ is parametrized by the integers $n$ and $\ell$, with $n=\ell,\ell+1, \dots$. We can now write, for each $n$, the corresponding $R_{n\ell}$ by using (\ref{recurrence_G}) to find the polynomial $G = G_{n\ell}$, and then $R_{n\ell}$ via (\ref{an_G}). Unwrapping all our definitions, \begin{gather} R_{n\ell}(r) = e^{-\frac{Zr}{na_0}} \left( \frac{Zr}{a_0}\right)^\ell G_{n\ell}\left(\frac{Zr}{a_0}\right), \nonumber \end{gather} where \begin{gather} a_0 = \frac{4\pi\varepsilon_0\hbar^2}{\mu e^2}. \nonumber \end{gather} In light of (\ref{separation_radial_angular}), we see that $\psi$ is also parametrized by $n, \ell,$ and $m_\ell$. Instead of thinking of $n$ varying according to $n = \ell, \ell +1, \dots$, we can equivalently think of $\ell$ as constrained by $\ell=0,1,\dots, n-1$, for each given $n=1,2,\dots$, what is more convenient in order to organize the parameters $n, \ell, m_\ell$. We obtain, therefore, a family of solutions to (\ref{time_independent}), \begin{gather} \psi_{n\ell m_\ell} = R_{n\ell} Y_{\ell m_\ell}, \label{psi_sol} \end{gather} where \begin{align} \begin{split} &n=1,2,3, \dots,\\ &\ell = 0, 1, 2, \dots, n-1,\\ &m_\ell = -\ell, -\ell+1, \dots, 0, \dots, \ell -1, \ell. \end{split} \label{relations_quantum_numbers} \end{align} Our final solution is then given, in view of (\ref{separation_t_x}) and (\ref{T_sol}), by \begin{gather} \Psi(t,x) = A_{n\ell m_\ell} e^{-\frac{i E_n}{\hbar} t} \psi_{n \ell m_\ell}, \nonumber \end{gather} where $n,\ell$, and $m_\ell$ satisfy (\ref{relations_quantum_numbers}), $E_n$ and $\psi_{n \ell m_\ell}$ are given by (\ref{E_n}) and (\ref{psi_sol}), respectively, and $A_{n\ell m_\ell}$ is a constant (depending on $n$, $\ell$, and $m_\ell$) that ensures (\ref{psi_normalized}), i.e., $A_{n\ell m_\ell}$ is given by \begin{gather} A_{n\ell m_\ell} = \left( \int_{\RR^3} |\psi_{n \ell m_\ell}|^2 \right)^{-\frac{1}{2}}. \nonumber \end{gather} The term $e^{-\frac{i E_n}{\hbar} t}$ does not contribute to $|\Psi|^2$ (since $(e^{-\frac{i E_n}{\hbar} t})^*(e^{-\frac{i E_n}{\hbar} t}) = (e^{+\frac{i E_n}{\hbar} t})(e^{-\frac{i E_n}{\hbar} t}) = 1$). It is customary to absorb the constant $A_{n\ell m_\ell}$ into $\psi_{n\ell m_\ell}$, in which case \begin{gather} \int_{\RR^3} |\psi_{n \ell m_\ell}|^2 =1. \nonumber \end{gather} Or course, (\ref{psi_normalized}) is automatically satisfied in this case. \section{Final comments.} We close with some remarks about the physical meaning of the problem we just described. Readers are referred to \cite{W} for a more thorough physical discussion. Below, we list some the first few $\psi_{n \ell m_\ell}$. \begin{center} \begin{tabular}{c|c|c|c} $n$ & $\ell$ & $m_\ell$ & $\psi_{n \ell m_\ell}$ \\ \hline 1 & 0 & 0 & $\psi_{100} = \frac{1}{\sqrt{\pi}}\left( \frac{Z}{a_0}\right)^\frac{3}{2} e^{-\frac{Zr}{a_0}}$\\ \hline 2 & 0 & 0 & $\psi_{200} = \frac{1}{4\sqrt{2\pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\left(2-\frac{Zr}{a_0}\right) e^{-\frac{Zr}{2a_0}}$\\ \hline 2 & 1 & 0 & $\psi_{210} = \frac{1}{4\sqrt{2\pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\frac{Zr}{a_0} e^{-\frac{Zr}{2a_0}} \cos \phi$\\ \hline 2 & 1 & $\pm$ 1 & $\psi_{21\pm1} = \frac{1}{8\sqrt{2\pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\frac{Zr}{a_0} e^{-\frac{Zr}{2a_0}} \sin \phi e^{\pm i \theta}$\\ \hline 3& 0 & 0 & $\psi_{300} = \frac{1}{81\sqrt{3\pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\left( 27-18\frac{Zr}{a_0}+2\frac{Z^2r^2}{z_0^2}\right) e^{-\frac{Zr}{3a_0}} $\\ \hline 3& 1 & 0 & $\psi_{310} = \frac{\sqrt{2}}{81\sqrt{\pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\left( 6 - \frac{Zr}{a_0}\right)\frac{Zr}{a_0} e^{-\frac{Zr}{3a_0}} \cos \phi$ \\ \hline 3 & 1 & $\pm$ 1 & $\psi_{31 \pm 1} = \frac{\sqrt{2}}{81\sqrt{\pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\left( 6 - \frac{Zr}{a_0}\right)\frac{Zr}{a_0} e^{-\frac{Zr}{3a_0}} \cos \phi e^{\pm i\theta}$\\ \hline 3& 2 & 0 & $\psi_{320} = \frac{1}{81\sqrt{6 \pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\frac{Z^2r^2}{a_0^2} e^{-\frac{Zr}{3a_0}} (3\cos^2 \phi - 1)$ \\ \hline 3& 2 & $\pm$ 1 & $\psi_{32\pm 1} = \frac{1}{81\sqrt{ \pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\frac{Z^2r^2}{a_0^2} e^{-\frac{Zr}{3a_0}} \sin\phi \cos \phi e^{\pm i \theta}$ \\ \hline 3& 2 & $\pm$ 2 & $\psi_{32\pm 2} = \frac{1}{262\sqrt{ \pi}} \left( \frac{Z}{a_0}\right)^\frac{3}{2}\frac{Z^2r^2}{a_0^2} e^{-\frac{Zr}{3a_0}} \sin^2\phi e^{\pm 2 i \theta}$ \\ \hline \end{tabular} \end{center} It is possible to show that the constants $E_n$, $\ell(\ell+1)$, and $m_\ell$ have important physical interpretation: $E_n$ corresponds to the electron energy, $\ell(\ell+1)$ to the magnitude of its orbital angular momentum, and $m_\ell$ to the projection of the orbital angular momentum onto the $z$-axis. The reader should notice that these quantities cannot be arbitrary, being allowed to take values only on a countable set of multiples of integers. This is a distinctive feature of Quantum Mechanics (we say that the energy and orbital angular momentum are ``quantized"). The indices $n$, $\ell$, and $m_\ell$ are called \emph{quantum numbers}. One-electron atoms with $\ell=0,1,2,3$ are labeled $s,p,d,f$. In hydrogen and hydrogen-like atoms, this letter is preceded by a number giving the energy level $n$. Thus, the lowest energy state of the hydrogen atom is $1s$; the next to the lowest are $2s$ and $2p$; the next $3s, 3p, 3d$ and so on. These are the so-called ``orbitals" that the student is likely to have learned in Chemistry. Remembering that $|\Psi|^2$ is a probability density, what these orbitals represent are ``clouds of probability," highlighting the regions of three-dimensional space where it is more likely to find the electron. We finish mentioning that in a more detailed treatment of the problem, $\mu$ is not exactly the mass of the particle being described, but rather the reduced mass of the system. This is because, strictly speaking, the electron does not orbit the proton, but both orbit the center of mass of the system electron--proton. This is very much like the situation of the Earth orbiting the Sun: both bodies move due to their reciprocal gravitational attraction, although the Sun, begin much more massive, barely feels the pull caused by Earth's gravitational field, and that is why one usually thinks of the Earth orbiting an standing-still Sun. A similar situation occurs for the proton and the electron. We remark, however, that the calculations we presented apply, with no change, to this more accurate situation: we only have to change the value of $\mu$ to be the reduced mass. \begin{thebibliography}{ZZZZ} \bibitem[B]{B} Butkov, E. \emph{Mathematical Physics.} Addison-Wesley Educational Publishers Inc; New edition (1973). \bibitem[FY]{FY} Faddeev, L. D., and Yakubovskii, O. A. \emph{Lectures on Quantum Mechanics for Mathematics Students.} American Mathematical Society (2009). \bibitem[T]{T} Takhtajan, L. A. \emph{Quantum Mechanics for Mathematicians.} American Mathematical Society (2008). \bibitem[W]{W} Weinberg, S. \emph{Lectures on Quantum Mechanics.} Cambridge University Press (2012). \end{thebibliography} \vskip 1cm \end{document}