Solutions of Schrödinger's equation for Rydberg states of hydrogen are shown to display patterns that can be understood semi-classically within the framework of Bohr-Sommerfeld planetary orbits and de Broglie's concept of matter waves. Using *n* = 10 and *n* = 100 quantum levels we show how de Broglie-wave-analyzed Kepler orbits can be matched to excited-state wavefunctions and charge distributions. Making the connection between the old and new quantum theories can help beginning students develop an intuition about quantum mechanics and allow them to develop an appreciation of how Schrödinger's wave mechanics was stimulated by the theories of his predecessors.

## I. CIRCULAR ORBITS AND “CIRCULAR STATES”

Niels Bohr's 1913 paper described an analytical method that built upon Ernest Rutherford's concept of an atom as a positively charged nucleus surrounded by a system of negatively charged electrons.^{1} An early success of Bohr's atomic theory was its ability to derive and explain the empirical formula of Johannes Rydberg for the frequencies of spectral lines of atomic hydrogen.^{2} Bohr analyzed the hydrogen atom as a single electron orbiting an infinitely massive positive nucleus (that is, a single proton). He simplified the problem by assuming further that the electron orbit was circular. Bohr used non-relativistic classical mechanics as the foundation for his analysis, but he also drew upon the work of Max Planck to understand radiation associated with an orbiting electron. Planck earlier had proposed the concept of the energy-quanta to explain blackbody radiation.^{3} In Bohr's model, rather than an orbiting electron radiating energy continuously, he suggested that radiation from excited hydrogen occurred in distinctly separated emissions (described today as jumps) between stationary states. Stationary states of hydrogen were those that satisfied an integer condition connected with Planck's idea of quantized action.

Bohr designated the integer factor in his analysis with *τ*—we will instead use *n*, which is known commonly as the principal quantum number. Bohr showed that the energy of the stationary states is proportional to $ \u2212 1 / n 2$, and the angular momentum of the stationary states is $ n h / 2 \pi $, where *h* is Planck's constant. Bohr identified a most tightly bound state (the ground state) corresponding to the smallest possible positive integer, *n* = 1. The radius of the circular orbits according to Bohr's model is $ n 2 a 0$, where *a*_{0} is the Bohr radius. The model also shows that the electron velocity in an orbit is $ v / c = ( 1 / n ) ( 2 \pi e 2 / ( 4 \pi \epsilon 0 h c ) )$. (Here, $ \epsilon 0$ is the permittivity of free space.^{4}) The numerical factor in the velocity relationship is Sommerfeld's fine structure constant, $ \alpha = 2 \pi e 2 / ( 4 \pi \epsilon 0 h c ) \u2248 1 / 137$. The velocity scales as $ 1 / n$. The energy of stationary states of the electron according to the Bohr model is $ \u2212 ( h c \u2009 R \u221e ) / n 2$, where $ R \u221e$ is the Rydberg constant and *c* is the speed of light. Bohr's derived formulas for both *a*_{0} and $ R \u221e$ were given in terms of fundamental constants: the electron mass, *m*_{0}; the electron charge, *e*; and Planck's constant, *h*. Bohr's value for $ R \u221e$ based on fundamental constants is very close to Rydberg's experimentally determined value.

The notion of a moving electron as a wave is often invoked in modern explanations of Bohr's quantum theory, but the concept did not emerge until about 10 years after Bohr's landmark paper. Bohr only introduced a stability condition. Louis de Broglie, drawing on the works of Maupertuis (principle of least action) and Fermat (principle of least time), reasoned that particles following a trajectory have a wave nature. In a paper published in 1925, de Broglie noted the connection between the concept of a particle wave and Bohr's stability condition:^{5} “Propagation is, therefore, analogue to a liquid wave in a channel closed on itself but of variable depth. It is physically obvious, that to have a stable regime, the length of the channel must be resonant with the wave; …The resonance condition is $ L = n \lambda $ if the wavelength is constant …” [*L* is de Broglie's value for the length of the channel and *n* is an integer.] Then he continued in a later paragraph, “In the particular case of closed circular Bohr orbits in an atom, one gets: $ m 0 \u222e v d L = 2 \pi R \u2009 m 0 v = n h$, where $ v = R \omega $, when *ω* is angular velocity, $ m 0 \omega \u2009 R 2 = n h / 2 \pi $. This is exactly Bohr's fundamental formula.” Since the circumference of the Bohr circle is $ 2 \pi R$, the wavelength extracted from Bohr's formula is $ \lambda = h / m 0 v$, that is, de Broglie's well known relationship. We see that the circular Bohr orbit contains *n* electron wavelengths, and the electron phase is a multiple of $ 2 \pi $ at *n* points equally spaced along the orbit.

Now we consider the Schrödinger wavefunction, which in a spherical basis $ ( r , \theta , \phi )$ is characterized by the quantum numbers, $ n , \u2113$, and *m*. We will focus on the “circular state”^{6} of quantum level *n*, which has maximal values of both $\u2113$ and *m*, namely, $ \u2113 = n \u2212 1$ and $ m = \xb1 \u2113$. Since we compare orbits to a slice of the wavefunction, it is helpful first to consider the bigger picture of what the electron charge distribution looks like for the *n* = 10 “circular state.” This is given in Fig. 1, which shows in 3D that the distribution appears as a torus centered on the origin. The charge density peaks at a radius somewhere between the inner and outer radii of the torus.

The sum of the squares of the real and the imaginary parts of the wavefunction is used to compute the probability density in Fig. 1. Figure 2 is a 2D contour plot of the real part of the wavefunction in the *x* – *y* plane. Red contours are for positive real values. Blue dashed contours are for negative real values. Black lines indicate where the real part of the wavefunction is zero.

The imaginary part of the wavefunction in the *x* – *y* plane has the same pattern as the real part shown in Fig. 2 but is rotated about the z-axis by $ + 10 \xb0$ (CCW) for the $ m = + 9$ state and by $ \u2212 10 \xb0$ (CW) for the *m* = – 9 state. Because the *n* = 10 “circular state” has 9-fold symmetry, a spatial rotation of $ 40 \xb0$ (that is, $ 360 \xb0 / 9$) in either the CCW or CW direction about *z* will leave the complex wavefunction unchanged. The $ \xb1 10 \xb0$ relative rotation between real and imaginary parts means that the azimuthal pattern of the real and imaginary parts oscillates and is in quadrature (that is, one quarter of the way to the symmetry rotation). This is fundamental.^{7} It allows the sum of the squares of real and imaginary parts to be constant as a function of the azimuthal angle $\varphi $. We are interested in using orbits to understand spatial patterns in wavefunctions. Since the patterns of real and imaginary parts differ only by a rotation, we restrict our comparison, without loss of generality, to the real part of the wavefunction. All subsequent mention of the wavefunction, therefore, will be the real part only.

We now compare Schrödinger's wavefunction for *n* = 10 “circular states” with Bohr's circular orbit for *n* = 9 and *n* = 10. Why consider *n* = 9? We include the *n* = 9 orbit because Schrödinger and Bohr analyses differ in their values for angular momentum. Schrödinger's value of angular momentum for “circular states” of quantum level *n* is $ ( n \u2212 1 ) h / 2 \pi $, whereas the circular orbits of Bohr for quantum level *n* have angular momentum $ n h / 2 \pi $, that is, $ h / 2 \pi $ larger. The Bohr method overestimates angular momentum of each quantum level—this well-known flaw of Bohr's theory^{8} was corrected by Schrödinger in 1926.^{9} This error leads Bohr to assume that the angular momentum of the *n* = 1 ground state is $ h / 2 \pi $, whereas Schrödinger's ground state has angular momentum zero. So, we will start by considering both *n* = 9 and *n* = 10 Bohr orbits in comparison with the *n* = 10 “circular state.” Plotted in Fig. 2 are contours of the Schrödinger wavefunction in the *x* – *y* plane as compared with these two Bohr orbits.

We see that while the *n* = 10 Bohr orbit does a good job at matching the anti-node peaks of the real part of the wavefunction, the red dots in that orbit miss the wavefunction phase (that is, the oscillations)—whereas the red dots of the *n* = 9 Bohr orbit do a good job at matching the phase – they do this because the angular momentum for the *n* = 9 Bohr case matches the angular momentum (and the symmetry) of the *n* = 10 circular Schrödinger case. Our approach then for comparing Schrödinger *n* = 10 results to the Bohr model is to use the *n* = 9 orbit but scaling it up in physical dimension by the factor, 100/81. Now both the size of the orbit and the phase will line up as shown in Fig. 4(a). For comparisons with states of higher values of *n*, for example, *n* = 100, as considered in Sec. III of this paper, the scaling is not necessary because the scale factor, $ ( n 2 ) / ( n \u2212 1 ) 2$, approaches 1 when n is very large, but we still will use a quantum number for Bohr orbits that is one less than that of the Schrödinger state.

## II. ELLIPTICAL ORBITS AND SPHERICAL STATES

In 1916, Sommerfeld expanded Bohr's approach by including the possibility of elliptical orbits.^{10–12} This possibility corrected a shortcoming of Bohr's model by providing an analysis that improved comparison with experimentally observed level splitting due to external electric and magnetic fields. To allow elliptical orbits, Sommerfeld postulated a new quantum stability condition involving multiple quantum numbers, that is, a quantum condition for each degree of freedom for an orbit of an electron. He assumed that the orbit was confined to a plane, and therefore, there were two quantum conditions and two quantum numbers.

Here, we compare the Bohr–Sommerfeld planar elliptical orbits with the solutions of atomic wavefunctions. To begin, we restrict our attention to the *x* – *y* plane and consider only Bohr–Sommerfeld orbits in that plane. The classical angular momentum is perpendicular to the orbital plane. Thus, for these orbits, the angular momentum points along the *z* axis. The Sommerfeld method, like that of Bohr, overestimates the total angular momentum by $ h / 2 \pi $. In a state with principal quantum number *n*, the maximal values for the quantum numbers $\u2113$ and *m* are *n* – 1. So, as in Sec. I, when we compare the Bohr–Sommerfeld model to the Schrödinger equation solutions, we compare cases of the same angular momentum with the same projection of angular momentum along the *z* axis. Therefore, we match Bohr–Sommerfeld orbits having quantum number, *n* – 1 (which have $\u2113$ values that range from 1 to *n* – 1) to Schrödinger wavefunctions for states having quantum number, *n* (with $\u2113$ values that also range from 1 to *n* – 1—the $ \u2113 = 0$ case is not considered). We also restrict the comparison to those *m* states that have the maximum angular momentum projection along *z—*that is, $ m = \xb1 \u2113$, for each case.

Sommerfeld noted that Bohr's circular orbits have only one degree of freedom and that made it relatively easy to apply Planck's quantized action condition. The elliptical case for motion in a plane involves two degrees of freedom, and, therefore, he postulated that the quantized action conditions for stability should also have two conditions: $ \u222e p \theta \u2009 d \theta = n \theta h$ and $ \u222e p r \u2009 d r = n r h$, where $ n \theta + n r = n$. Here, $ \u2009 p \theta = m 0 r 2 d \theta / d t$ is the moment of momentum, and $ p r = m 0 d r / d t$ is the radial component of momentum, which is *m*_{0} times the radial velocity. Following this approach, he concluded that there are *n* orbits that satisfy the stability condition for the quantum level *n*: the circular orbit and *n* – 1 additional elliptical orbits. The possible solutions have angular momenta ranging from $ n h / 2 \pi $ down to $ h / 2 \pi $. These correspond to $ n \theta $ ranging from *n* down to 1. The case $ n \theta = n$ is Bohr's circular orbit. Sommerfeld eliminated the case of $ n \theta = 0$ arbitrarily, arguing that this case would have the electron collide with the nucleus. Sommerfeld computed the *n* = 4 case as an example. He noted that the circle and all the ellipses share a common semi-major axis, $ n 2 a 0$. Figure 3(a) is an image taken directly from Sommerfeld's original paper.^{10}

^{13}Knowing how the velocity varies along the orbit will allow us then to determine the locations along the elliptical path where the phase is a multiple of $ 2 \pi $. mathematica is used to solve for the Kepler orbits according to the parametric equations

Here, *ε* is the eccentricity, *a* is the semi-major axis of the ellipse, *b* is the semi-minor axis of the ellipse, *ω* is the orbital frequency in radians/seconds, and *M* is the mean anomaly. The equation for *M* is Kepler's equation. Sommerfeld's model provides expressions for $ a , b , \epsilon $, and *ω* in terms of *n* and $\u2113$: $ a = n 2 a 0$; $ b = n \u2009 \u2113 \u2009 a 0$; $ \epsilon = ( 1 \u2212 ( b 2 / a 2 )$; and $ \omega = ( m 0 e 4 ) / ( ( 4 \pi \epsilon 0 ) 2 ( h / 2 \pi ) 3 n 3 )$. To solve for *x*(*t*) and *y*(*t*), we invert Kepler's equation to solve for $ \theta ( t )$. This cannot be done analytically—so we invert it numerically using the “InverseFunction” function in mathematica. Once we have $ \theta ( t )$, we make a table of $ \theta ( t )$, *x*(*t*), and *y*(*t*). This table is based on a time step of typically $ 10 \u2212 4$– $ 10 \u2212 6$ of the orbital period. The table entries let us solve for the electron velocity along the elliptical path.

To figure out the points along the path where the phase advances by a multiple of $ 2 \pi $, we rely on the action phase-integral, $ ( \u222e p d q ) / h$. We assume a single degree of freedom for the electron as it moves along the ellipse. To get *p* at time step $ t + \Delta t$, we use the chain rule. We differentiate Eqs. (1) and (2) with respect to *θ* and determine $ d \theta / d t$ by differentiating Eq. (3). Then $ p = m 0 ( d \theta / d t ) ( [ d x / d \theta ] 2 + [ d y / d \theta ] 2 )$, and $ d q = d \theta ( [ d x / d \theta ] 2 + [ d y / d \theta ] 2 )$. We determine $ d \theta $ by taking the differences in the table for $ \theta ( t )$ at times $ t + \Delta t$ and *t*. We start at the perihelion (the point of closest separation between the electron and the nucleus) and define the phase at that point to be equal to zero. The phase tabulation begins at *t* = 0. The electron moves along the elliptical path, and for each new time step, we record the accumulated phase in an additional column in the table. We plot the orbit using the *x*(*t*) and *y*(*t*) values in the table and place a red dot at those points where the accumulated phase reaches a multiple of $ 2 \pi $.

This was done for the *n* = 4 case as shown in Fig. 3(b). As expected, when the electron is moving quickly near the perihelion, the electron wavelength is smaller than it is when near the aphelion (the most distant separation between electron and the nucleus). For *n* = 4, there are four electron wavelengths in the circular orbit and in each of the ellipses.

In Fig. 3(c), we have done this also for the *n* = 9 case. Here, there are nine electron wavelengths in the circular orbit and in each of the eight additional ellipses. This family of orbits will be used to compare with the *n* = 10 Schrödinger wavefunctions.

Now we overlay the scaled *n* = 9 Bohr–Sommerfeld orbits on the *n* = 10 Schrödinger wavefunctions. We do this for each of the possible $\u2113$ values (leaving out $ \u2113 = 0$—so $ \u2113 = 9 , 8 , \u2026 , 2 , 1$). The comparisons are for equal values of $\u2113$ and only for the sub-levels that have a maximum projection of angular momentum on the *z* axis. Thus, if we are making comparisons for the $ \u2113 = 9$ state, we choose $ m = \xb1 9$. If we are comparing the $ \u2113 = 8$ state, we choose $ m = \xb1 8$ for comparison, and so on. In Fig. 4, *n*= 9 Bohr–Sommerfeld orbits are overlayed on *n* = 10 Schrödinger wavefunctions for four cases.

Figure 4(a) is the circular state, $ n = 10 , \u2113 = 9 , m = \xb1 9$, compared to the *n* = 9 Bohr and Sommerfeld model solution. The state $ n = 10 , \u2113 = 8 , m = \xb1 8$, shown in Fig. 4(b), displays several interesting features. It has 8-fold symmetry, one less than the $ n = 10 , \u2113 = 9 , m = \xb1 9$ state. Along the radial direction, a second ring of anti-nodes appears. Recall Sommerfeld's suggestion of two quantum numbers, $ n \theta $ and *n _{r}*—so the reduction of the circumferential symmetry by 1 and the increase in the radial rings by 1 is consistent with Sommerfeld's idea that $ n \theta + n r = n$. As we will see for lower $\u2113$ states, as the circumferential symmetry drops by 1, the ring of radial anti-nodes will increase by 1. Next, we see that the perihelion of the elliptical orbit aligns with the location of the inner ring. This trend will continue with lower $\u2113$ states. Each time a new inner ring is added by decreasing $\u2113$, the location of the inner-most ring of anti-nodes will match the location of the perihelion of the elliptical orbit. We see also that there is rough correspondence with the location of the positive circumferential anti-nodes with the red dots along the ellipse.

Note that the overall size of the Schrödinger $ n = 10 , \u2113 = 8 , m = \xb1 8$ wavefunction appears larger than the $ n = 10 , \u2113 = 9 , m = \xb1 9$ wavefunction. This can be seen as the result of sweeping the ellipse about the nucleus—recall that the semi-major axis of all of the ellipses and the circle is $ n 2 a 0$. According to Bethe and Salpeter,^{14} the average value of the radius for a Schrödinger state of angular momentum $\u2113$ is $ ( 1 / 2 ) ( 3 n 2 \u2212 \u2113 ( \u2113 + 1 ) ) a 0$. That is, the smaller $\u2113$, the larger the average size of the state—only when $\u2113$ is *n* – 1 does the size of the state approach the classical value of $ n 2 a 0$. So according to this relationship, the average radius of Schrödinger's $ n = 10 , \u2113 = 9$ state is $ 105 a 0$ (as compared to $ 100 a 0$), and the radius of the $ n = 10 , \u2113 = 8$ state is $ 114 a 0$ (a larger value than $ 105 a 0$). The expected average radius of the $ n = 10 , \u2113 = 1$ state should be $ 149 a 0$ (that is, 50% larger than $ 100 a 0$).

The next two wavefunctions in the *n* = 10 sequence are $ \u2113 = 7 , m = \xb1 7$, and $ \u2113 = 6 , m = \xb1 6$, shown in Figs. 4(c) and 4(d). These states also match up with the de Broglie-wave-analyzed orbits. The lower $\u2113$ states $ ( \u2113 = 4 , 3 , 2 ,$ and 1) continue to follow these trends—as the $\u2113$ value decreases, the circumferential symmetry decreases, and the number of radial rings increases. All this is consistent with Sommerfeld's hypothesis that $ n \theta + n r = n$.

One of the reviewers of this paper asked if the correlation of the de Broglie phase between the Bohr and Bohr–Sommerfeld models and the actual Schrödinger wavefunctions is significant or merely coincidence? This is a good question worthy of discussion.

We know that classical elliptical orbits have been used with much success to understand quantum mechanical behavior of hydrogen and hydrogen-like atoms. For example, Horbatsch, Hessels, and Horbatsch in their paper entitled, “Classical calculation of the lifetimes and branching ratios for radiative decays of hydrogenic atoms,”^{15} found that they could accurately estimate quantum mechanical radiative decay in hydrogen using semi-classical methods for all states having principal quantum numbers *n* and $ \u2113 \u2265 1$. They assumed that each $ n , \u2113$ state had a single well-defined eccentricity that depended on the energy of the state and its angular momentum. Their agreement for lifetimes was excellent (<100 parts per million error for $ \u2113 \u2265 30 , < 0.1 %$ for $ \u2113 \u2265 9 , < 1 %$ for $ \u2113 \u2265 3$), and their agreement for branching ratios was quite reasonable. Their calculation was based on Larmor radiation of a circulating electron in orbit around the nucleus. This same type of calculation is posed in a question at the end of Chapter 14 of Jackson's *Classical Electrodynamics* textbook.^{16} In addition to the study of lifetimes and branching ratios, there are other successful studies using classical orbits to analyze quantum phenomenon in hydrogen, and recent studies give us confidence that looking for elliptical orbits where the eccentricity is well defined is sensible and not likely due to coincidence.

However, perhaps the best test is whether the method works with other sublevels. So, in response to a reviewer's query, we examined the wavefunctions for another group of states, namely, *m* = 0 states having $ \u2113 = 9 , 8 , \u2026 , 2 , 1$. We found that all of the *m* = 0 wavefunctions follow the same patterns as the case discussed earlier, but now the patterns lie in any plane that contains the *z* axis. In the *m* = 0 case, the classical Bohr–Sommerfeld orbit used for comparison with the quantum result will lie in the same *z*-axis-containing plane. Since the angular momentum vector is perpendicular to the orbital plane, the projection of the angular momentum along *z* is also zero. Shown in Fig. 5 are two of these patterns for the case of $ \u2113 = 5$. The one on the left is the wavefunction pattern in the *x* – *y* plane for the $ m = \xb1 5$ sublevels, and the one on the right is the wavefunction pattern in the *x* – *z* plane for the *m* = 0 sublevel. The two are strikingly similar. With these additional states, we find now that 27 of the 100 sublevels of *n* = 10 quantum level have patterns of nodes and anti-nodes that can be explained using these semi-classical ideas. Observe that the charge distributions for the *m* = 0 case and the $ m = \xb1 5$ cases are very different from one another—the charge distribution of the *m* = 0 case is a collection of spherical rings, whereas the charge distributions of the $ m = \xb1 5$ cases are planar rings. These are shown in Fig. 6. It is evident that the semi-classical analysis can help one understand nodal patterns in charge distributions as well as in wavefunctions.

## III. ELLIPTICAL ORBITS AND SPHERICAL STATES—CHARGE DISTRIBUTIONS COMPARED RADIALLY

Using the method of separation of variables, the Schrödinger wavefunction for hydrogen can be factored into two parts, $ R n \u2113 ( r )$ and $ Y \u2113 m ( \theta , \phi )$. Here, $ R n \u2113 ( r )$ is the solution to the hydrogen radial equation. $ Y \u2113 m ( \theta , \phi )$ are the spherical harmonics. Because of this factorization, we can determine the charge probability density distribution as a function of *r*. The Schrödinger radial distribution is given by $ r 2 [ R n \u2113 ( r ) ] 2$. Our objective in this section is to compare this radial distribution to the radial distribution expected for de Broglie-wave-analyzed Bohr–Sommerfeld orbits. We will be able to do this because Sommerfeld's stability conditions are in two parts: one that is radial and another that is azimuthal. This quantum-classical comparison has two aspects. One involves matching the amplitude envelope of the quantum probability density to a classical analysis of the radial probability density for Kepler orbits. The other involves matching the phase.

To match the amplitude envelope, we note that the radial probability density for a classical electron in orbit is proportional to $ 1 / [ d r / d t ]$, which for Keplerian orbits is known to be $ r / ( \pi a ( \epsilon 2 a 2 \u2212 ( r \u2212 a ) 2 ) )$, where *ε* is the eccentricity and *a* is the size of the semi-major axis of the ellipse.^{17,18} Comparisons between the amplitude envelope and the classical radial probability distribution have been noted by authors in previous works.^{17,19,20}

To get the de Broglie phase, we apply Sommerfeld's radial stability condition: $ \u222e ( p r / h ) d r = n r$. Here, *p _{r}* is $ m 0 d r / d t$. We tabulate

*r*(

*t*) in a Bohr–Sommerfeld orbit in a manner similar to that in Sec. II, using $ r ( t ) = [ x ( t ) ] 2 + [ y ( t ) ] 2$, where $ x ( \theta ( t ) )$ and $ y ( \theta ( t ) )$ are given in Eqs. (1) and (2). We get

*dr*/

*dt*by applying the chain rule, $ d r / d t = ( d r / d \theta ) ( d \theta / d t )$, and $ d r = [ ( d r / d t ) ( d t ) ]$. $ d \theta / d t$ is determined by differentiating Eq. (3). Previously with the phase integral along the ellipse, we noted every time the phase advanced to a multiple of $ 2 \pi $, but in this case, we will note the values of

*r*where the phase advances to a multiple of

*π*because we are plotting the probability density rather than the wavefunction. Note that the method here considers only the radial path,

*r*, using the differential radial phase of $ ( p r / h ) d r$. In contrast, in Sec. II, the method for determining the phase was based a one-dimensional elliptical path,

*q*, using the differential electron phase, $ ( p / h ) d q$.

In Fig. 7, the radial charge distribution for the $ n = 100 , \u2113 = 80$ Schrödinger state is compared to the classical probability distribution for the n $ = 99 , \u2113 = 80$ Bohr–Sommerfeld orbit. The classical distribution is the solid black line. This figure also has red dots placed at the radii, where the phase advances by multiples of *π*. The semi-classical result helps one to understand aspects of the quantum result. One sees that the average amplitude of the oscillatory solution matches the classical value. One also sees that the anti-nodes of the oscillatory solution are in reasonable correspondence with the red dots. We are using the *n* = 99 Bohr–Sommerfeld orbit for comparison here for the same reasons that we did in Secs. I and II. In this case, the *n* = 99 Bohr–Sommerfeld orbital size is very close to the *n* = 100 orbital size, so there was no need to scale the results as we did previously.

## IV. ELLIPTICAL ORBITS AND PARABOLIC STATES

Sommerfeld used parabolic coordinates for his comparisons to the Stark effect. He used the Hamilton–Jacobi method of separation of variables to develop three quantum stability conditions according to the parabolic coordinate variables $ \xi , \eta $, and $\varphi $. These stability conditions gave rise to quantum numbers, *n*_{1}, *n*_{2}, and *m*. Schrödinger's equation can be analyzed using parabolic coordinates and also using the method of separation of variables. As pointed out by Bethe and Salpeter,^{14} the form of the perturbation potential for the electric field along the *z* axis preserves the method of separation of variables since the variable *z* is a linear combination of *ξ* and *η*, namely, $ z = 0.5 ( \xi \u2212 \eta )$. The other variables are $ x = \xi \eta \u2009 cos \u2009 ( \varphi ) ; y = \xi \eta \u2009 sin \u2009 ( \varphi )$; and $ \varphi = arctan ( y / x )$. Parabolic states are used often to analyze the Stark effect in atomic hydrogen. For a state to have a linear shift in the electric field, it means that in the absence of a field, there would be an electric dipole moment. That is true for parabolic states but not for spherical states, which are symmetric or anti-symmetric with respect to inversion. The energy shift of a state with an electric dipole moment is the product of the electric dipole moment and the external field. If the state has an electric dipole moment in the absence of a field, then application of an external field leads to a shift that is linear in the field—that is what happens in excited atomic hydrogen.

Figure 8 is a contour plot of the charge distribution of a parabolic state for the quantum level, *n* = 10 in the *x* – *z* plane. This state has the largest electric dipole moment of the 100 sub-levels of *n* = 10 and the largest linear shift in an external electric field–it also has *m* = 0. A Bohr–Sommerfeld ellipse in the *x* – *z* plane has its angular momentum vector along the *y* axis and, therefore, would also have *m* = 0. We know from Schrödinger's analysis that parabolic states do not have a well-defined angular momentum. Instead, they are made up of a superposition of multiple angular momentum states.^{19} The fractional component is determined by projecting the parabolic state onto spherical states of various values of $\u2113$. One can develop insight about the Schrödinger result if we compare a collection of Bohr–Sommerfeld elliptical orbits to an *x* – *z* slice of the Schrödinger charge distribution—we could have also used a *y* – *z* slice, or any slice that contains the *z* axis—this state has rotational symmetry about *z*. In Fig. 8, four Bohr–Sommerfeld orbits having $ \u2113 = 1 , 2 , 3$, and 4 are overlayed on the quantum charge distribution. For these four Bohr–Sommerfeld orbits, the angular momentum vector points along y. The contours shown are for a slice in the *x* – *z* plane. The red dots appear at multiples of *π* in phase along the orbits. You will note that the anti-nodes of the parabolic state roughly line up with the red dots. We see that the size and de Broglie phase of the classical elliptical orbit help one to understand both the oscillations and the amplitude of the quantum charge density.

## V. CONCLUSIONS

Wavefunctions and charge densities for atomic hydrogen display patterns that are understandable in terms of the planetary models of the old quantum theories. The seemingly complex myriad of nodes and anti-nodes of hydrogen are recognizable in comparison with de Broglie-wave-analyzed Bohr–Sommerfeld orbits. One need not dismiss the old quantum theories so quickly when teaching or studying quantum mechanics. Appreciation of physical ideas like the de Broglie phase wave can help students look for features in quantum solutions that trigger their imagination. How might a classical orbiting electron behave? If the wavefunction or the charge distribution is rapidly changing in space, then the classical interpretation is that the electron is moving quickly. Schrödinger's 1926 paper in the *Physical Review*^{9} reminds us of how important classical mechanics and de Broglie's ideas were in his developing wave mechanics. He started off his paper as follows: “The theory which is reported in the following pages is based on the very interesting and fundamental researches of L. de Broglie on what he called ‘phase wave’ (‘ondes de phase’) and thought to be associated with the motion of material points, especially with the motion of an electron or proton. The point of view taken here, which was first published in a series of German papers, is rather that material points consist of, or are nothing but, wave-systems.” Schrödinger was inspired by classical motion and the phase waves of de Broglie. So, too, can these ideas help the beginning student see the connections between quantum solutions and classical ideas and help the student to develop an understanding of why quantum wavefunctions and charge distributions look the way they do.

## ACKNOWLEDGMENTS

The authors thank Marko Horbatsch of York University for helpful discussions. The authors are also grateful to the reviewers for useful suggestions and corrections.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

## REFERENCES

*Atombau und Spektrallinien*

*Atomic Structure and Spectral Lines*

*Quantum Mechanics of One and Two Electron Atoms*

*Classical Electrodynamics*