Next Article in Journal
The Entropic Dynamics of Quantum Scalar Fields Coupled to Gravity
Next Article in Special Issue
Symmetry, Transactions, and the Mechanism of Wave Function Collapse
Previous Article in Journal
Asymmetrical Pseudo-Extinction Phenomenon in the Illusory Line Motion
Previous Article in Special Issue
Inflation with Scalar-Tensor Theory of Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Dynamical Symmetries of the H Atom, One of the Most Important Tools of Modern Physics: SO(4) to SO(4,2), Background, Theory, and Use in Calculating Radiative Shifts

Quantum Fields LLC, St. Charles, IL 60174, USA
Symmetry 2020, 12(8), 1323; https://doi.org/10.3390/sym12081323
Submission received: 18 June 2020 / Revised: 13 July 2020 / Accepted: 31 July 2020 / Published: 7 August 2020
(This article belongs to the Special Issue Symmetries in Quantum Mechanics)

Abstract

:
Understanding the hydrogen atom has been at the heart of modern physics. Exploring the symmetry of the most fundamental two body system has led to advances in atomic physics, quantum mechanics, quantum electrodynamics, and elementary particle physics. In this pedagogic review, we present an integrated treatment of the symmetries of the Schrodinger hydrogen atom, including the classical atom, the SO(4) degeneracy group, the non-invariance group or spectrum generating group SO(4,1), and the expanded group SO(4,2). After giving a brief history of these discoveries, most of which took place from 1935–1975, we focus on the physics of the hydrogen atom, providing a background discussion of the symmetries, providing explicit expressions for all of the manifestly Hermitian generators in terms of position and momenta operators in a Cartesian space, explaining the action of the generators on the basis states, and giving a unified treatment of the bound and continuum states in terms of eigenfunctions that have the same quantum numbers as the ordinary bound states. We present some new results from SO(4,2) group theory that are useful in a practical application, the computation of the first order Lamb shift in the hydrogen atom. By using SO(4,2) methods, we are able to obtain a generating function for the radiative shift for all levels. Students, non-experts, and the new generation of scientists may find the clearer, integrated presentation of the symmetries of the hydrogen atom helpful and illuminating. Experts will find new perspectives, even some surprises.

1. Introduction

1.1. Objective of This Paper

This pedagogic review is focused on the symmetries of the Schrodinger nonrelativistic hydrogen atom exclusively to give it the attention that we believe it deserves. The fundamental results of the early work are known and do not need to be derived again. However, having this knowledge permits us to use the modern language of group theory to do a clearer, more focused presentation, and to use arguments from physics to develop the proper forms for the generators, rather than dealing with detailed, mathematical derivations to prove results we know are correct.
There are numerous articles about the symmetry of the Schrodinger hydrogen atom, particularly the SO(4) group of the degenerate energy eigenstates, including discussions from classical perspectives. The spectrum generating group SO(4,1) and the non-invariance group SO(4,2) have been discussed, but, in many fewer articles, often in appendices, with different bases for the representations. For example, numerous papers employ Schrodinger wave functions in parabolic coordinates, not with the familiar n l m quantum numbers, often with the emphasis on the details of the mathematical structures, not the physics, and with the emphasis on the potential role of the symmetry in elementary particle physics. Generators may be expressed in complex and unfamiliar terms, for example, in terms of the raising and lowering operators for quantum numbers characteristic of parabolic coordinates. Other approaches involve regularizing Schrodinger’s equation by, for example, multiplying by r. The approach results in generators that are not always Hermitian or manifestly Hermitian, and the need for nonstandard inner products. Indeed, most of the seminal articles do not have the words “hydrogen atom” in the title but are focused on the prize down the road, understanding what was called at the time “the elementary particle zoo”.
Yet, these foundational articles and books taken together present the information we have about the symmetry of the H atom, of which experts in the field are aware. On the other hand, to the non-expert, the student, or a researcher new to the field, it does not appear that the relevant information is in a form that is conveniently accessible. Since the hydrogen atom is the most fundamental physical system with an interaction, whose exploration and understanding has led to much of the progress in atomic physics, quantum physics, and quantum electrodynamics, we believe that a comprehensive treatment is warranted and, since most of the relevant papers were published four decades ago, is timely. Many younger physicists may not be acquainted with these results.
Unlike in a number of the foundational papers, here the operators are all Hermitian, and given in terms of the canonical position and momentum variables in the simplest forms. The transformations they generate are clearly explained, and we provide brief explanations for the group theory used in the derivations.
In most papers, a separate treatment for bound and scattering states in needed. In contrast we are able to clarify and simplify the exposition since we use a set of basis states which are eigenfunctions of the inverse of the coupling constant [1], that include both the bound states and the scattering states in a uniform way, and that employ the usual Cartesian position and momenta, with the usual inner product, with the exact same quantum numbers n l m as the ordinary bound states; a separate treatment for bound and scattering states is not required. In addition, we have two equivalent varieties of this uniform basis, one that is more suitable for momentum space calculations and one more suitable for configuration space calculations. This advantage again allows us to simplify the exposition.
We focus on the utility of group theoretic methods using our representation and derive expressions for the unitary transformation of group elements and some new results that allow for us to readily compute the first order radiative shift (Lamb shift) of a spinless electron, which accounts for about 95% of the total shift. This approach allows for us to obtain a generating function for the shifts for all energy levels. For comparison, we derive an expression for the Bethe log.
In summary, we present a unified treatment of the symmetries of the Schrodinger hydrogen atom, from the classical atom to SO(4,2) that focuses on the physics of the hydrogen atom, that gives explicit expressions for all the manifestly Hermitian generators in terms of position and momenta operators in a Cartesian space, that explains the action of the generators on the basis states, which evaluates the Casimir operators characterizing the group representations, and that gives a unified treatment of the bound and continuum states in terms of wave functions that have the same quantum numbers as the ordinary bound states. We give an example of the use of SO(4,2) in a practical application, the computation of the first order radiative shift in the hydrogen atom.
Hopefully, students and non-experts and the new generation of scientists will find this review helpful and illuminating, perhaps motivating some to use these methods in various new contexts. Senior researchers will find new perspectives, even some surprises and encouragements.

1.2. Outline of This Paper

In the remainder of Section 1 we give a brief historical account of the role of symmetry in quantum mechanics and of the work done in order to explore the symmetries of the Schrodinger hydrogen atom.
In Section 2, we provide some general background observations regarding symmetry groups and non-invariance groups and discuss the degeneracy groups for the Schrodinger, Dirac, and Klein-Gordon equations. We also introduce the uncommon ( Z α ) 1 eigenstates that allow us to treat the bound and scattering states in a uniform way, using the usual quantum numbers. In Section 3, the classical equations of motion of the nonrelativistic hydrogen atom in configuration and momentum space are derived from symmetry considerations. The physical meaning of the symmetry transformations and the structure of the degeneracy group SO(4) is discussed. In Section 4, we discuss the symmetries using the language of quantum mechanics. In order to display the symmetries in quantum mechanics in the most elegant and uniform way we use a basis of eigenstates of the inverse of the coupling constant ( Z α ) 1 . In Section 5 we discuss these wave functions in momentum and configuration space, how they transform and their classical limit for Rydberg states.
In Section 6, we discuss the noninvariance or spectrum generating group of the hydrogen atom SO(4,l) and relate it to the conformal group in momentum space. In Section 7, the enlarged spectrum generating group SO(4,2) is introduced, with a discussion of the physical meaning of the generators. All of the physical states together form a basis for a unitary irreducible representation of these noninvariance groups. We derive manifestly Hermitian expressions in terms of the momentum and position canonical variables for the generators of the group transformations and obtain the values possible for the Casimir operators. We discuss the important subgroups of SO(4,2).
In Section 8, we use the group theory of SO(4,2) to determine the radiative shifts in energy levels due to the interaction of a spinless electron with its own radiation field, or equivalently with the quantum vacuum. In the nonrelativistic or dipole approximation the level shift contains a matrix element of a rotation operator of an O(1,2) subgroup of the group SO(4,2). We can sum this over all states, obtaining the character of the representation, yielding a single integral that is a generating function for the radiative shift for any level in the nonrelativistic or dipole approximation. A brief conclusion follows.

1.3. Brief History of Symmetry in Quantum Mechanics and Its Role in Understanding the Schrodinger Hydrogen Atom

The hydrogen atom is the fundamental two-body system and perhaps the most important tool of atomic physics and the continual challenge is to continually improve our understanding of the hydrogen atom and to calculate its properties to the highest accuracy possible. The current QED theory is the most precise of any physical theory [2]:
The study of the hydrogen atom has been at the heart of the development of modern physics...theoretical calculations reach precision up to the 12th decimal place...high resolution laser spectroscopy experiments...reach to the 15th decimal place for the 1S–2S transition...The Rydberg constant is known to six parts in 10 12 [2,3]. Today, the precision is so great that measurement of the energy levels in the H atom has been used to determine the radius of the proton.
Continual progress in understanding the properties of the hydrogen atom has been central to progress in quantum physics [4]. Understanding the atomic spectra of the hydrogen atom drove the discovery of quantum mechanics in the 1920’s. The measurement of the Lamb shift in 1947 and its explanation by Bethe in terms of atom’s interaction with the quantum vacuum fluctuations ushered in a revolution in quantum electrodynamics [5,6,7]. Exploring the symmetries of the hydrogen atom has been an essential part of this progress. Symmetry is a concept that has played a broader role in physics in general; for example, in understanding the dynamics of the planets, atomic, and molecular spectra, and the masses of elementary particles.
When applied to an isolated system, Newton’s equations of motion imply the conservation of momentum, angular momentum and energy. But the significance of these conservation laws was not really understood until 1911 when Emily Nother established the connection between symmetry and conservation laws [8]. Rotational invariance in a system results in the conservation of angular momentum; translational invariance in space results in conservation of momentum; and translational invariance in time results in the conservation of energy. We will discuss Nother’s Theorem in more detail in Section 2.
Another critical ingredient of knowledge, on which Nother based her proof, was the idea of an infinitesimal transformation, such as a infinitesimal rotation generated by the angular momentum operators in quantum mechanics. These ideas of infinitesimal transformations originated with the Norwegian mathematician Sophus Lie who was studying differential equations in the latter half of the nineteenth century. He studied the collection of infinitesimal transformations that would leave a differential equation invariant [9]. In 1918, German physicist and mathematician Hermann Weyl, in his classic book with the translated title “The Theory of Groups and Quantum Mechanics”, would refer to this collection of differential generators leaving an operator invariant as a linear algebra, ushering in a little of the terminology of modern group theory [10]. Still this was a very early stage in understanding the role of symmetry in the language of quantum theory. When he introduced the new idea of a commutator on page 264, he put the word “commutator” in quotes. In the preface Weyl made a prescient observation: “.. the essence of the new Heisenberg-Schrodinger-Dirac quantum mechanics is to be found in the fact that there is associated with each physical system a set of quantities, constituting a non-commutative algebra in the technical mathematical sense, the elements of which are the physical quantities themselves”.
A few years later Eugene Wigner published in German, “Group Theory and Its Application to the Quantum Mechanics of Atomic Spectra [11]”. One might ask why was this classic not translated into English until 1959. In the preface to the English edition, Prof. Wigner recalled: “When the first edition was published in 1931, there was a great reluctance among physicists toward accepting group theoretical arguments and the group theoretical point of view. It pleases the author that this reluctance has virtually vanished..” It was the application of group theory in particle physics in the early sixties, such as SU(3) and chiral symmetry, which reinvigorated interest in Wigner’s book and the field in general. In the 1940’s, Wigner and Bargmann developed the representation theory of the Poincare group that later provided an infrastructure for the development of relativistic quantum mechanics [11,12].
The progress in understanding the symmetries of the hydrogen atom, in particular, has some parallels to the history of symmetry in general: there were some decades of interest but after the 1930’s interest waned for about three decades in both fields, until stimulated by the work on symmetry in particle physics.
Probably, the first major advance in understanding the role of symmetry in the classical treatment of the Kepler problem after Newton’s discovery of universal gravitation, elliptical orbits, and Kepler’s Laws, was made two centuries ago by Laplace when he discovered the existence of three new constants of the motion in addition to the components of the angular momentum [13]. These additional conserved quantities are the components of a vector which determines the direction of the perihelion of the motion (point closest to the focus) and whose magnitude is the eccentricity of the orbit. The Laplace vector was later rediscovered by Jacobi and has since been rediscovered numerous times under different names. Today it is generally referred to as the Runge-Lenz vector. However, the significance of this conserved quantity was not well understood until the nineteen thirties.
In 1924, Pauli made the next major step forward in understanding the role of symmetry in the hydrogen atom [14]. He used the conserved Runge-Lenz vector A and the conserved angular momentum vector L to solve for the energy spectrum of the hydrogen atom by purely algebraic means, a beautiful result, yet he did not explicitly identify that L and A formed the symmetry group SO(4) corresponding to the degeneracy. At this time, the degree of degeneracy in the hydrogen energy levels was believed to be n 2 for a state with principal quantum number n, clearly greater than the degeneracy due to rotational symmetry which is ( 2 l + 1 ) . The n 2 degeneracy arises from the possible values of the angular momentum l = 0 , 1 , 2 , n 1 , and the 2 l + 1 values of the angular momentum along the azimuthal axis m = l , l + 1 , 0 , 1 , 2 , l + 1 . The additional degeneracy was referred to as “accidental degeneracy [15].”
Six years after Pauli’s paper, Hulthen used the new Heisenberg matrix mechanics to simplify the derivation of the energy eigenvalues of Pauli by showing that the sum of the squares L 2 + A 2 could be used to express the Hamiltonian and so could be used to find the energy eigenvalues [16]. In a one sentence footnote in this three page paper, Hulthen gives probably the most important information in his paper: Prof. Otto Klein, who had collaborated for years with Sophus Lie, had noticed that the two conserved vectors formed the generators of the Lorentz group, which we can describe as rotations in four dimensions, the fourth dimension being time. This is the non-compact group SO(3,1), the special orthogonal group in four dimension whose transformations leave the magnitude g μ ν z μ z ν = t 2 + x 2 + y 2 + z 2 unchanged [17]. Klein’s perceptive observation triggered the introduction of group theory to understanding the hydrogen atom.
About a decade later, in 1935, the Russian physicist Vladimir Fock published a major article in Zeitshrift fur Physik, the journal in which all the key articles about the hydrogen atom cited were published [18]. He transformed Schrodinger’s equation for a given energy eigenvalue from configuration space to momentum space, and did a stereographic projection onto a unit sphere, and showed that the bound state momentum space wave functions were spherical harmonics in four dimensions. He stated that this showed that rotations in four dimensions corresponded to the symmetry of the degenerate bound state energy levels in momentum space, realizing the group SO(4), the group of special orthogonal transformations which leaves the norm of a four-vector U 0 2 + U 1 2 + U 2 2 + U 3 2 constant. By counting the number of four-dimensional spherical harmonics Y n l m in momentum space ( m = l , l + 1 0 , 1 , l , where the angular momentum l can equal l = n 1 , n 2 , 0 ) , he determined that the degree of degeneracy for the energy level characterized by the principal quantum number n was n 2 . It is interesting that Fock did not cite the work by Pauli, implying the four dimensional rotational symmetry in configuration space. Fock also presented some ideas about using this symmetry in calculating form factors for atoms.
A year later, the German-American mathematician and physicist Valentine Bargmann showed that for bound states (E < 0) Pauli’s conserved operators, the angular momentum L and the Runge–Lenz vector A, obeyed the commutation rules of the SO(4) [12]. His use of commutators was so early in the field of quantum mechanics, that Bargmann explained the square bracket notation he used for a commutator in a footnote [19]. He gave differential expression for the operators, adapting the approach of Lie generators in the calculation of the commutators. He linked solutions to Schrodinger’s equation in parabolic coordinates to the existence of the conserved Runge–Lenz vector and was thereby able to establish the relationship of Fock’s results to the algebraic representation of SO(4) for bound states implied by Fock and Pauli [12]. He also pointed out that the scattering states (E > 0) could provide a representation of the group SO(3,1). In a note at the end of the paper, Bargmann, who was at the University in Zurich, thanked Pauli for pointing out the paper of Hulthen and the observation by Klein that the Lie algebra of L and A was the same as the infinitesimal Lorentz group, which is how he referred to a Lie algebra. Bargmann’s work was a milestone demonstrating the relationship of symmetry to conserved quantities and it clearly showed that to fully understand a physical system one needed to go beyond the usual ideas of geometrical symmetry. This work was the birth, in 1936, without much fanfare, of the idea of dynamical symmetry.
Little attention was paid to these developments until the 1960’s, when interest arose primarily because of the applications of group theory in particle physics, particularly modeling the mass spectra of hadrons. Particle physicists were faced with the challenge of achieving a quantitative description of hadron properties, particularly the mass spectra and form factors, in terms of quark models. Since little was know about quark dynamics they turned to group-theoretical arguments, exploring groups like SU(3), chiral U(3)xU(3), U(6)xU(6), etc. The success of the eight-fold way of SU(3) (special unitary group in three dimensions) of American physicist Murray Gell-Mann in 1962 brought attention to the use of symmetry considerations and group theory as tools for exploring systems in which one was unsure of the exact dynamics [20].
In 1964, three decades after Fock’s work, American physicist Julian Schwinger published a paper using SO(4) symmetry to construct a Green function for the Coulomb potential, which he noted was based on a class he taught at Harvard in 1949 [21]. The publication was a response to the then current emphasis on group theory and symmetry, which led to, as Israeli physicist Yuval Ne’eman described it, “‘the great-leap-forward’ in particle physics during the years 1961–1966 [22]”. Some of the principal researchers leading this effort were Ne’eman [23], Gell-Mann [20,24,25], Israeli physicist Y. Dothan [26], Japanese physicist Yochiri Nambu [27], and English-American Freeman Dyson [28]. Advantage was taken of the mathematical infrastructures of group theory developed years earlier [10,11,12,29,30].
Interest was particularly strong in systems with wave equations with an infinite number of components, which characterize non-compact groups. In about 1965, this interest in particle physics gave birth to the identification of SO(4,1) and SO(4,2) as Spectrum Generating Algebras that might serve as models for hadronic masses. The hydrogen atom was seen as a model to explore the infinite dimensional representations of non-compact groups. The first mention of SO(4,1) was by Barut, Budini, and Fronsdal [31], where the H atom was presented as an illustration of a system characterized by non-compact representation, and so comprising an infinite number of states. The first mention of a six dimensional symmetry, referred to as the “non-compact group O(6)”, appears to be by the Russian physicists I. Malin and V. Man’ko of the Moscow Physico-technical Institute [32]. In a careful three page paper, they showed that all of the bound states of the H atom energy spectrum in Fock coordinates provided a representation of this group, and they calculated the Casimir operators for their symmetric tensor representation in parabolic coordinates.
Very shortly thereafter, Turkish-American theoretical physicist Asim Barut and his student at University of Colorado, German theoretical physicist Hagen Kleinert, showed that including the dipole operator e r as a generator led to the expansion of SO(4,1) to SO(4,2), and that all the bound states of the H atom formed a representation of SO(4,2) [33]. This allowed them to calculate dipole transition matrix elements algebraically. They give a position representation of the generators based on the use of parabolic coordinates. The generators of the transformations are given in terms of the raising and lowering operators for the quantum numbers for solutions to the H atom in parabolic coordinates. The dilation operator is used to go from one SO(4) subspace with one energy to a SO(4) subspace with different energy and it has a rather complicated form. They also used SO(4,2) symmetry to compute form factors [34].
The papers of the Polish-American physicist Myron Bander and French physicist Claude Itzakson published in 1966, when both were working at SLAC (Stanford Linear Accelerator in California) provide the first mathematically rigorous and “succinct” review of the O(4) symmetry of the H atom and provide an introduction to SO(4,1) [35,36], which is referred to as a spectrum generating algebra SGA, meaning that it includes generators that take the basis states from one energy level to another. They use two approaches in their mathematical analysis, the first is referred to as “the infinitesimal method,” based on the two symmetry operators, L and A and the O(4) group they form, and the other, referred to as the “global method”, first done by Fock, converts the Schrodinger equation to an integral equation with a manifest four dimensional symmetry in momentum space. They establish the equivalence of the two approaches by appealing to the solutions of the H atom in parabolic coordinates, and demonstrate that the symmetry operators in the momentum space correspond to the symmetry operators in the configuration space. As they note, the stereographic projection depends on the energy, so the statements for a SO(4) subgroup are valid only in a subspace of constant energy. They then explore the expansion of the SO(4) group to include scale changes so the energy can be changed, transforming between states of different principal quantum number, which correspond to different subspaces of SO(4). To insure that this expansion results in a group, they include other transformations, which results in the the generators forming the conformal group O(4,1). Their mathematical analysis introducing SO(4,1) is based on the projection of a p dimensional space (4 in the case of interest) on a parabaloid in p + 1 dimensions (5 dimensions). In their derivation they treat bound states in their first paper [35] and scattering states in the second paper [36].
As we have indicated the interest in the SO(4,2) symmetry of the Schrodinger equation was driven by a program focused on developing equations for composite systems that had infinite multiplets of energy solutions and ultimately could lead to equations that could be used to predict masses of elementary particles, perhaps using other than four dimensions [27,35,36,37,38,39,40]. In 1969 Jordan and Pratt showed that one could add spin to the generators A and L , and still form a SO(4) degeneracy group. By defining J = 1 2 ( L + A ) + S , they showed one could obtain a representation of O(4,1) for any spin s [41].
In their review of the symmetry properties of the hydrogen atom, Bander and Itzakson emphasize this purpose for exploring the group theory of the hydrogen atom [35]:
The construction of unitary representations of non-compact groups that have the property that the irreducible representations of their maximal subgroup appear at most with multiplicity one is of certain interest for physical applications. The method of construction used here in the Coulomb potential case can be extended to various other cases. The geometrical emphasis may help to visualize things and provide a global form of the transformations.
Special attention was also given to solutions for the hydrogen atom from the two body Bethe-Salpeter equation for a proton and electron interacting by a Coulomb potential, since the symmetry was that of a relativistic non-compact group [27,40,42,43].
Finally in 1969, five years after it was published, Schwinger’s form of the Coulomb Green’s function based on the SO(4) symmetry was used to calculate the Lamb shift by Michael Lieber, one of Schwinger’s students at Harvard [21,44]. A year later, Robert Huff, a student of Christian Fronsdal at UCLA, focused on the use of the results from SO(4,2) group theory to compute the Lamb shift [45]. He converted the conventional expression for the Lamb shift into a matrix element containing generators of SO(4,2), and was able to perform rotations and scale changes to simplify and evaluate the matrix elements. After clever mathematical manipulation, he obtained an expression for the Bethe log in terms of a rapidly terminating series for the level shifts. He provided an appendix with a brief discussion of the fundamental of SO(4,2) representations for the H atom, showing the expressions for the three generators needed to express the Schrodinger equation.
In the next few years, the researchers published a few mathematically oriented papers [41,46,47,48,49,50], a short book [51] dealing with the symmetries of the Coulomb problem, and a paper by Barut presenting a SO(4,2) formulation of symmetry breaking in relativistic Kepler problems, with a 1 page summary of the application of SO(4,2) for the non-relativistic hydrogen atom [34,52]. Bednar published a paper applying group theory to a variety of modified Coulomb potentials, which included some matrix elements of SO(4,2) using hydrogen atom basis states with quantum numbers n l m [53]. There also was interest in application of the symmetry methods and dynamical groups in molecular chemistry [54] and atomic spectroscopy [55].
In the 1970’s, researchers focused on developing methods of group theory and on understanding dynamical symmetries in diverse systems [56,57,58]. A book on group theory and its applications appeared in 1971 [59]. Barut and his collaborators published a series of papers dealing with the hydrogen atom as a relativistic elementary particle, leading to an infinite component wave equation and mass formula [60,61,62,63].
Papers on the classical Kepler problem, the Runge–Lenz vector, and SO(4) for the hydrogen atom have continued to appear over the years, from 1959 to today. Many were published in the 1970’s [64,65,66,67,68,69,70] and some since 1980, including [71,72,73,74,75]. Papers dealing with SO(4,2) are much less frequent. In 1986 Barut, A. Bohm, and Ne’eman published a book on dynamical symmetries that included some material on the hydrogen atom [76]. In 1986, Greiner and Muller published the second edition of Quantum Mechanics Symmetries, which had six pages on the Hydrogen atom, covering only the SO(4) symmetry [77]. The 2005 book by Gilmore on Lie algebras has 4 pages of homework problems on the H atom to duplicate results in early papers [78]. The last papers I am aware of that used SO(4,2) were applications in molecular physics [79,80] and more general in scope [81]. Carl Wulfman published a book on dynamical symmetries in 2011, which provides a helpful discussion of dynamical symmetries for the hydrogen atom [82]. He regularizes the Schrodinger equation, essentially multiplying by r, obtaining Sturmian wave functions in parabolic coordinates. This approach allows for him to treat bound and scattering states for SO(4,2) at one time, but requires redefining the inner product, and leads to a non-Hermitian position operator.

1.4. The Dirac Hydrogen Atom

We have focused our discussion on the symmetries of the non-relativistic hydrogen atom described by the Schrodinger equation. Quantum mechanics also describes the hydrogen atom in terms of the relativistic Dirac equation, which we will only discuss briefly in this paper.
The gradual understanding of the dynamical symmetry of the Dirac atom parallels that of the Schrodinger atom, but it has received much less attention, probably because the system has less relevance for particle physics and for other applications. It was known that the rotational symmetry was present and that the equation predicted that the energy depended on the principal quantum number and the quantum number for the total angular momentum j, but not the spin s or orbital angular momentum l separately. This remarkable fact meant that, in some sense, angular momentum contributed the same to the total energy no matter whether it was intrinsic or orbital in origin. This degeneracy is lifted if we include the radiative interactions which leads to the Lamb shift.
In order to understand the symmetry group for the Dirac equation consider that for a given total angular momentum quantum number J > 0 there are two degenerate levels for each energy level of the Dirac hydrogen atom: one level has l = J + 1/2 and the other has l = J − 1/2. Since the l values differ by unity, the two levels have opposite parity. Dirac described a generalized parity operator K, which was conserved. For an operator Λ to transform one degenerate state into the other, it follows that the operator has to commute with J and have parity −1. This means it has to anticommute with K, and so it is a conserved pseudoscalar operator.
The parity ( 1 ) l + J 1 / 2 is conserved in time, so the states are parity eigenstates. Using the two symmetry operators Λ and K, one can build a SU(2) algebra. If we include the O(3) symmetry due to the conservation of angular momentum, we obtain the full symmetry group S U ( 2 ) x O ( 3 ) which is isomorphic to S O ( 4 ) for the degeneracy of the Dirac hydrogen atom.
In 1950, M.Johnson and B. Lippman discovered the operator Λ [83]. Further work was done on understanding Λ by Biedenharn [84]. The Johnson–Lippman operator has been rediscovered and reviewed several times over the decades [85,86,87]. It has been interpreted in the non-relativistic limit as the projection of the Runge–Lenz vector onto the spin angular momentum [87,88,89]. The SO(4) group can be expanded to include all states, and then the spectrum generating group is SO(4,1) or SO(4,2) depending on the assumptions regarding relativistic properties and the charges present [33,38]. We will not discuss the symmetries of the Dirac H atom further.

2. Background

2.1. The Relationship between Symmetry and Conserved Quantities

The nature of the relationship between symmetry, degeneracy, and conserved operators is implicit in the equation
[ H , S ] = 0
where H is the Hamiltonian of our system, S is a Hermitian operator, and the brackets signify a commutator if we are discussing a quantum mechanical system, or i times a Poisson bracket if we are discussing a classical system. If S is viewed as the generator of a transformation on H, then Equation (1) says the transformation leaves H unchanged. Therefore, we say S is a symmetry operator of H and leaves the energy invariant. The fact that a non- trivial S exists means that there is a degeneracy. To show this, consider the action of the commutator on an energy eigenstate | E :
[ H , S ] | E = 0
or
H ( S | E ) = E ( S | E ) .
If S is nontrivial then S | E is a different state from | E but has the same energy eigenvalue. If we label all such degenerate states by
| E , m , m = 1 , , N
then clearly S | E , m is a linear combination of degenerate states:
S | E , m = S m s | E , s .
S m s is a matrix representation of S in the subspace of degenerate states. In a classical Kepler system S generates an orbit deformation that leaves H invariant. The existence of a nontrivial S therefore implies degeneracy, in which multiple states have the same energy eigenvalue. We can show that the complete set of symmetry operators for H forms a Lie algebra by applying Jacobi’s identity to our set of Hermitian operators S i :
H , S i = 0 , i = 1 , , L .
[ S j , [ H , S i ] ] + [ S i , [ S j , H ] ] + [ H , [ S i , S j ] ] = 0
so
[ H , [ S i , S j ] ] = 0 .
Therefore, the commutator of S i and S j is a symmetry operator of H. Either the commutator is a linear combination of all the symmetry operators S i , i = 1 , , L :
[ S i , S j ] = C i j k S k
or the commutator defines a new symmetry operator which we label S L + 1 . We repeat this procedure until the Lie algebra closes as in Equation (9).
By exponentiation, we assume that we can locally associate a group of unitary transformations
e i S i a i
for real a i with our Lie algebra and so conclude that a group of transformations exists under which the Hamiltonian is invariant [90]. We call this the symmetry or degeneracy group of H. Our energy eigenstates states form a realization of this group.
It is possible to form scalar operators, called Casimir operators, from the generators of the group that commute with all the generators of the group, and, therefore, have numerical values. The values of the Casimir operators characterize the particular representation of the group. For example, for the rotation group in three dimensions, the generators are L = ( L 1 , L 2 , L 3 ) and the quantity L 2 = L ( L + 1 ) commutes with all of the generators. L can have any positive integer value for a particular representation. The Casimir operator for O(3) is L 2 . The number of Casimir operators that characterize a group is called the rank of the group. O(3) is rank 1 and SO(4,2) is rank 3.
Now let us consider Equation (1) in a different way. If we view H as the generator of translations in time, then we recall that the total time derivative of an operator S i is
d S i d t = i H , S i ˙ + S i t
where the commutator and the partial derivative give the implicit and explicit time dependence respectively. Provided that the symmetry operators have no explicit time dependence ( S i d t = 0 ) , then Equation (11) implies that Equation (1) means that the symmetry operators S are conserved in time and d S i d t = 0 . Conversely, we can say that conserved Hermitian operators with no explicit time dependence are symmetry operators of H. This very important relationship between conserved Hermitian operators and symmetry was first discovered by German mathematician Emmy Nother in 1917, and is called Nother’s Theorem [8,91,92,93,94,95].

2.2. Non-Invariance Groups and Spectrum Generating Group

As we have discussed, the symmetry algebra contains conserved generators S i that transform one energy eigenstate into a linear combination of eigenstates all with the same energy. In order to illustrate with hydrogen atom eigenstates:
S i | n l m = l , m S n l m n l m | n l m
where | n l m refers to a state with energy E n , angular momentum l ( l + 1 ) and l z = m .
A non-invariance algebra contains generators D i that can be used to transform one energy eigenstate | n l m into a linear combination of other eigenstates, with the same or a different energy, different angular momentum l, and different azimuthal angular momentum m:
D | n l m = n , l , m D n l m n l m | n l m .
Because the set of energy eigenstates is complete, the action of the most general operator would be identical to that shown to Equation (13). Therefore this requirement alone is not sufficient to select the generators needed.
The goal is to expand the degeneracy group with its generators S i into a larger group, so that some or all of the eigenstates form a representation of the larger group with the degeneracy group as a subgroup. Thus, we need to add generators G i , such that the combined set of generators
{ S i , G j ; for all i , j } { D k ; for all k }
forms an algebra that closes
[ D i , D j ] = i ϵ i j k D k .
This is the Lie algebra for the expanded group. To illustrate with a specific example, consider the O(4) degeneracy group with six generators [96]. One can expand the group to O(5) or O(4,1) which has ten generators by adding a four-vector of generators. One component might be a scalar and the other three a three-vector. The question then is can some or all of the energy eigenstates provide a representation of O(5)? If so, then this would be considered a non-invariance group. The group might be expanded further in order to obtain generators of a certain type or to include all states in the representation. For the H atom the generators D i can transform between different energy eigenvalues meaning between eigenfunctions with different principal quantum numbers.
Another way to view the expansion of the Lie algebra of the symmetry group is to consider additional generators D i that are constants in time [97] but do not commute with the Hamiltonian so
d D i d t = 0 = i H , D i ˙ + D i t .
If we make the additional assumption that the time dependence of the generators is harmonic
2 D i ( t ) t 2 = ω i n D n ( t ) .
then generators D i , and the first and second partial derivatives with respect to time could close under commutation, forming an algebra. This approach does not tell us what generators to add, but, as we demonstrate in Section 7.5, it does reflect the behavior of the generators that have been added to form the spectrum generating group in the case of the hydrogen atom.
We may look for the largest set of generators D i , which can transform the set of solutions into itself in an irreducible fashion (meaning no more generators than necessary). These generators form the Lie algebra for the non-invariance or spectrum generating algebra [32,98]. If the generators for the spectrum generating algebra can be exponentiated, then we have a group of transformations for the spectrum generating group. The corresponding wave functions form the basis for a single irreducible representation of this group. This group generates transformations among all the solutions for all energy eigenvalues and it is called the Spectrum Generating Group [99]. For the H atom, SO(4,1) is a spectrum generating group or non-invariance group, which can be reduced to contain one separate SO(4) subgroup for each value of n.
To get a representation of SO(4,1), we need an infinite number of states, which we have for the H atom. This group has been expanded by adding a five vector to form SO(4,2) because the additional generators can be used to express the Hamiltonian and the dipole transition operator. The group SO(p,q) is the group of orthogonal transformations that preserve the quantity X = x 1 2 + x 2 2 + + x p 2 x p + q 2 , which may be viewed as the norm or a p+q-dimensional vector in a space that has a metric with p plus signs and q minus signs. The letters SO stand for special orthogonal, meaning the orthogonal transformations have determinant equal to +1.
In terms of group theory, there is a significant difference between a group like SO(4) and SO(4,1). SO(4) and SO(3) are both compact groups, while SO(4,1) and SO(4,2) are non-compact groups. A continuous group G is compact if each function f(g), continuous for all elements g of the group G, is bounded. The rotation group in three dimensions O(3), which conserves the quantity r 2 = x 1 2 + x 2 2 + x 3 2 , is an example of a compact group.
For a non-compact group, consider the Lorentz group O(3,1) of transformations to a coordinate system moving with a velocity v. The transformations preserve the quantity r 2 c 2 t 2 . The matrix elements of the Lorentz transformations are proportional to 1 / 1 β 2 , where β = v / c , and are not bounded as β 1 . Therefore, r and c t may increase without bound, while the difference of the squares remains constant. Unitary representations of non-compact groups are infinite dimensional, for example, the representation of the non-invariance group SO(4,1) has an infinite number of states. Unitary representations of compact groups can be finite dimensional, for example, our representation of SO(4) for an energy level E n has dimension n 2 .
In the nineteen sixties and later, the spectrum generating group was of special interest in particle physics, because it was believed it could provide guidance where the precise particle dynamics were not known. The hydrogen atom provided a physical system as a model. Because the application was in particle physics, there was less interest in exploring representations in terms of the dynamical variables for position and momentum.
The expansion of the group from SO(4,1) to SO(4,2) was motivated by the fact that the additional generators could be used to write Schrodinger’s equation entirely in terms of the generators, and to express the dipole transition operator. This allowed for algebraic techniques and group theoretical methods to be used to obtain solutions, calculate matrix elements, and other quantities [33,38].

2.3. Basic Idea of Eigenstates of ( Z α ) 1

We briefly introduce the idea behind these states, since they are unfamiliar [1]. The full derivation is given in Section 4. Schrodinger’s equation in momentum space for bound states can be written as
p 2 + a 2 2 m Z α r | a = 0 .
where a 2 = 2 m E > 0 and Z α is the coupling constant, which we will now view as a parameter. This equation has well behaved solutions for certain discrete eigenvalues of the energy or a 2 , namely
a n 2 = 2 m E n
where E n = 1 2 m ( Z α ) 2 n 2 . We can write the eigenvalue condition equivalently as
( a n m Z α ) = 1 n .
This last equation shows that solutions exist for certain values a n of the RMS momentum a. To introduce eigenstates of ( Z α ) 1 we simply take a different view of this last equation and say that instead of quantizing a and obtaining a n , we imagine that we quantize ( Z α ) 1 , let a remain unchanged, obtaining
a m ( Z α ) n = 1 n .
Accordingly, now we can interpret Schrodinger’s equation as an eigenvalue equation that has solutions for certain values of ( Z α ) 1 namely
( Z α ) n 1 = m a n .
We have the same equation but can view the eigenvalues differently but equivalently. Instead of quantizing a we quantize ( Z α ) 1 .
This roughly conveys the basic idea of eigenstates of the inverse of ( Z α ) , but this simplified version does at all reveal the advantages of our reformulation because we have left the Hamiltonian unchanged. In Section 4, we transform Schrodinger’s equation to an eigenvalue equation in a, so that the kernel is bounded, which means that there are no states with E > 0, no scattering states, and all states have the usual quantum numbers. Other important advantages to this approach will also be discussed.

2.4. Degeneracy Groups for Schrodinger, Dirac and Klein-Gordon Equations

The degeneracy groups for the bound states described by the different equations of the hydrogen atom are summarized in Table 1. The degeneracy (column 2) is due to the presence of conserved operators which are also symmetry operators (column 3), forming a degeneracy symmetry group (column 4). For example, The symmetry operators for the degeneracy group in the Schrodinger hydrogen atom are the angular momentum L and the Runge–Lenz vector A . In Section 3, it will be shown that together these are the generators for the direct product SO(3)xSO(3) which is isomorphic to SO(4). Column 5 gives the particular representations present. These numbers are the allowed values of the Casimir operators for the group and they determine the degree of degeneracy (last column) and the corresponding allowed values of the quantum numbers for the degenerate states.
The Casimir operators, which are made from generators of the group, have to commute with all of the members of the group, and the only way this can happen is if they are actually constants for the representation. The generators are formed from the dynamical variables of the H atom, so the Casimir operators are invariants under the group composed of the generators, and their allowed numerical values reflect the underlying physics of the system and determine the appropriate representations of the group [9,82,100]. For example, L 2 is the Casimir operator for the group O(3) and can have the values l ( l + 1 ) . The relationship between Casimir operators and group representations is true for all irreducible group representations, including the SO(4) degeneracy group, as well as the spectrum generating group SO(4,2) [43,53].
For the Schrodinger equation, there are n 2 states | n l m that form a representation of the degeneracy group SO(4) formed by L and A . These states correspond to the principal quantum number n, the n different values of the angular momentum quantum number l , and 2 l + 1 different values of the z component of the angular momentum l Z = m .
For the Dirac equation, the 2 ( 2 J + 1 ) dimensional degeneracy group for bound states is realized by the total angular momentum operator J , the generalized parity operator K, and the Johnson–Lippman operator Λ , which together form the Lie algebra for SO(4).
For the fully relativistic Klein–Gordon equation, only the symmetry from rotational symmetry survives, leading to the degeneracy group O(3). If the V 2 term, the four-potential term squared is dropped in a semi-relativistic approximation as we describe in Section 4.3, then the equation can be rewritten in the same form as the non-relativistic Schrodinger equation, so a Runge–Lenz vector can be defined and the degeneracy group is again SO(4).

3. Classical Theory of the H Atom

In order to discuss orbital motion and the continuous deformation or orbits we give this discussion in terms of classical mechanics, but much of it is valid in terms of the Heisenberg representation of quantum mechanics if the Poisson brackets are converted to commutators, as will be discussed in Section 4.
For a charged particle in a Coulomb potential, there are two classical conserved vectors: the angular momentum L, which is perpendicular to the plane of the orbit, and the Runge–Lenz vector A, which goes from the focus corresponding to the center of mass and force along the semi-major axis to the perihelion (closest point) of the elliptical orbit. The conservation of A is related to the fact that non-relativistically the orbits do not precess. The Hamiltonian of our bound state classical system with an energy E < 0 is [17]
H = p 2 2 m Z α r = E
where m=mass of the electron, r is the location of the electron, p is its momentum, α is the fine structure constant, and E is the total non-relativistic energy.
The Runge–Lenz vector is
A = 1 2 m H ( p × L m Z α r r )
where L is the angular momentum. From Hamilton’s equation, H = E , so
A = p × L a m Z α a r r
where a is defined by
a = 2 m E .
From the virial theorem, the average momentum p 2 = 2 m E so a is the root mean square momentum. We are discussing bound states so E < 0. It is straightforward to verify that A is conserved in time:
[ A , H ] = d A d t = 0 .
From the definition of A and the definition of angular momentum
L = r × p
if follows that A is orthogonal to the angular momentum vector
A · L = 0 .
Using the fact that A and L are conserved, we can easily obtain equations for the orbits in configuration and momentum space and the eccentricity, and other quantities, all usually derived by directly solving the equations of motion.
We will show that A and L are the generators of the group O(4). If we introduce the linear combinations N = 1 2 ( L + A ) and M = 1 2 ( L A ) , we find that N and M commute reducing the nonsimple group O(4) to SU(2)x SU(2), which we will discuss in Section 4.2 in the language of quantum mechanics.

3.1. Orbit in Configuration Space

In order to obtain the equation of the orbit one computes
r · A = r A c o s ϕ r = r m Z α a + r · p × L .
Noting that r · p × l = L 2 , we can solve for r
r = L 2 / m Z α ( a / m Z α ) A c o s ϕ r + 1 .
This is the equation of an ellipse with eccentricity e = ( a / m Z α ) A and a focus at the origin (Figure 1). To find e we calculate A · A using the identity p × L · p × L = p 2 L 2 and obtain
A 2 = p 2 L 2 a 2 2 m Z α a 2 L 2 r + m Z α a 2 .
Substituting E for the Hamiltonian Equation (15) gives the usual result
a m Z α A = e = 2 E L 2 m ( Z α ) 2 + 1 .
The length r c of the semi-major axis is the average of the radii at the turning points
r c = r 1 + r 2 2 .
Using the orbit equation we find
r c = L 2 m Z α 1 1 e 2
or
r c = Z α 2 E = m Z α a 2 .
The energy depends only on the length of the semi-major axis r c , not on the eccentricity. This important result is a consequence of the symmetry of the problem. It is convenient to parameterize the eccentricity in terms of the angle ν (see Figure 1) where
e = sin ν
From this definition and from Equations (25), (27), and (28) follow the useful results
L = r c a cos ν , A = r c a sin ν
which immediately imply
L 2 + A 2 = r c a 2 = m Z α a 2 .
This equation is the classical analogue of an important quantum mechanical result first obtained by Pauli and Hulthen allowing us to determine the energy levels from symmetry properties alone [14,16]. From Figure 1, it is apparent that this equation is a statement of Pythagoras’s theorem for right triangles.
The energy equation (Equation (15)) and the orbit equation (Equation (23)) respectively may be rewritten in terms of a , r, and ν :
r c r = p 2 + a 2 2 a 2
r = r c cos 2 ν 1 + sin ν cos ϕ r .

3.2. The Period

To obtain the period, we use the geometrical definition of the eccentricity
e = 1 ( b / r c ) 2
where b is the semi-minor axis. Using e = sin ν we find
b = r c cos ν
so from Equation (30), we obtain
L = a b .
From classical mechanics, we know the magnitude of the angular momentum is equal to twice the mass times the area swept out by the radius vector per unit time. The area of the ellipse is π b r c . It the period of the classical motion is T, then L = 2 m π b r c / T . Therefore, the classical period is
T = 2 π m r c a = 2 π m ( Z α ) 2 8 E 3 .
and the classical frequency ω c l = 2 π / T is
ω c l = a m r c .

3.3. Group Structure SO(4)

The generators of our symmetry operations form the closed Poisson bracket algebra of O ( 4 ) :
[ L i , L j ] = i ϵ i j k L k , [ L i , A j ] = i ϵ i j k A k , [ A i , A j ] = i ϵ i j k L k .
The brackets mean i times the Poisson bracket, which is the classical limit of a commutator. The first bracket says that the angular momentum generates rotations and forms a closed Lie algebra corresponding to O(3). The second bracket says that the Runge-Lenz vector transforms as a vector under rotations generated by the angular momentum. The last commutator says that the multiple transformations generated by the Runge-Lenz vector are equivalent to a rotation. Taken together the commutators form the Lie algebra of O(4). The connected symmetry group for the classical bound state Kepler problem is obtained by exponentiating our algebra giving the symmetry group SO(4). The scattering states with E > 0 form an infinite dimensional representation of the non-compact group SO(3,1).
We now want to determine the nature of the transformations generated by A i and L i . Clearly, L · δ ω generates a rotation of the elliptical orbit about the axis δ ω by an amount δ ω . To investigate the transformations generated by A · δ ν we assume a particular orientation of the orbit, namely that it is in the 1–2 (or x-y) plane and that A is along the 1-axis (see Figure 1). The more general problem is obtained by a rotation generated by L i . For an example, we choose a transformation with δ ν pointing along the 2-axis, so that A · δ ν = A 2 δ ν . The change in A is defined by δ A , where
δ A = i [ A · δ ν , A ] .
From the Poisson bracket relations, we find for this particular case:
δ A 1 = L 3 δ ν , δ A 2 = 0 , δ A 3 = L 1 δ ν
For our orbit, L 1 = 0 so δ A 3 = 0 . We perform a similar computation to find δ L . We find we can characterize the transformation by
δ A 1 = L 3 δ v or δ A = L δ ν δ L 3 = A 1 δ ν or δ L = A δ ν δ e = 1 e 2 δ ν or δ ( sin ν ) = cos ν δ ν .
Recalling e = sin ν and Equation (30), we see that these transformations are equivalent to the substitution
ν ν + δ ν .
In other words, the eccentricity of the orbit, and therefore A and L are all changed in such a way that the energy, a and r c (length of the semi-major axis) remain constant. In our example, both L and A are changed in length but not direction, so the plane and orientation of the orbit are unchanged. The general transformation A · δ μ will also rotate the plane of the orbit or the semi-major axis.
Figure 2 shows a set of orbits in configurations space with different values of the eccentricity e = sin ν , but the same total energy and the same semimajor axis r c , which is the bold hypotenuse. The bold vertical and horizontal legs are A / a and L / a and they are related to the hypotenuse r c by Pythagoras’s theorem. The generator A 2 ν produces a deformation of the circular orbit into the various elliptical orbits shown. This classical degeneracy corresponds to the quantum mechanical degeneracy in energy levels that occurs for different eigenvalues of the angular momentum with a fixed principal quantum number.
We can visualize all possible elliptical orbits for a fixed total energy or semi-major axis through a simple device. It is possible to produce an elliptical orbit with eccentricity sin ν as the shadow of a circle of radius r c which is rotated an amount ν about an axis perpendicular to the illuminating light. With a complete rotation of the circle, we will see all possible classical elliptical orbits corresponding to a given total energy. In quantum mechanics only certain angles of rotation would be possible corresponding to the quantized values of L. As the circle is rotated, we must imagine that the force center shifts as the sine of the angle of rotation, so that it always remains at the focus [101].

3.4. The Classical Hydrogen Atom in Momentum Space

We can derive the equation for the classical orbit in momentum space of a particle bound in a Coulomb potential using the conserved operators L and A . For convenience, we assume that we have rotated our axes, so that L lies along the 3-axis and A the 1-axis, as shown in Figure 1. We compute
p · A = p 1 A = m Z α a p · r r m Z α a p r
and we employ Equations (28) and (30) to show A = m Z α a sin ν , to obtain [102]
p r = sin ν p 1
which we substitute into the identity
p r 2 + L 2 r 2 = p 2 = p 1 2 + p 2 2
Using Equations (30) and (32) we find
1 = 2 a p 1 p 2 + a 2 2 + 2 a p 2 p 2 + a 2 2 1 cos 2 v
and
p 2 a 2 = 2 a p 2 tan ν
which may also be written as
p 1 2 + p 2 a tan ν 2 = a 2 cos 2 ν .
From Equation (49), we see the orbit in momentum space is a circle of radius a / cos ν with its center displaced from the origin a distance a tan ν along the 2-axis. Figure 3 shows the momentum space orbit that corresponds to the configuration space orbit in Figure 1.
As an alternative method of showing the momentum space orbits are circular, we can compute [103]
p a L × A L 2 2 = C 2 .
Using the lemma
p × A = L 2 a ( p 2 a 2 ) ,
the fact that L · A = 0 , and Equation (30), we find C = a cos ν . The orbit is a circle of radius a c o s ν whose center lies at a L × A L 2 , in agreement with the previous result.
We now consider what the generators A i and L i do to the orbit in momentum space. Clearly, L generates a rotation of the axes. For an A i transformation, consider the same situation that we considered in our discussion of the configuration space orbit (see Figure 1 and Figure 3). Because the generator A 2 δ ν changes ν to ν + δ ν , we conclude that in momentum space this shifts the center of the orbit along the 2-axis and changes the radius of the circle. However, the distance a from the 2-axis to the intersection of the orbit with the 1-axis remains unchanged. Figure 4 shows a set of momentum space orbits for a fixed energy that correspond to the set of orbits in configuration space shown in Figure 2.

3.5. Four-Dimensional Stereographic Projection in Momentum Space

It is interesting that in classical mechanics the bound state orbits in a Coulomb potential are simpler in momentum space than in configuration space. In quantum mechanics the momentum space wave functions become simply four-dimensional spherical harmonics if one normalizes the momentum p by dividing by the RMS momentum a = 2 m E and performs a stereographic projection onto a unit hyper-sphere in a four-dimensional space [18,35]. We will do the analogous projection procedure for the classical orbits. The three-dimensional momentum space hyperplane passes through the center of the four-dimensional hypersphere, as shown in Figure 5. The unit vector in the fourth direction is n ^ = ( 1 , 0 , 0 , 0 ) . The unit vector U ^ goes from the center of the sphere to the surface of the hypersphere, where it is intersected by the line connecting the vector p / a to the north pole of the sphere.
We find
U i = 2 a p i p 2 + a 2 i = 1 , 2 , 3 . U 4 = p 2 a 2 p 2 + a 2
or inverting,
p i = a U i 1 U 4 p 2 = a 2 1 + U 4 1 U 4 .
Momentum space vectors for which p / a < 1 are mapped onto the lower hyperhemisphere. The advantage of this projection over one in which the hypersphere is tangent to the hyperplane is that we may have | n ^ | = | U ^ | = 1 . At times it is convenient to describe U ^ in terms of spherical polar coordinates in four dimensions. Because U ^ is a unit vector we define
U 4 = cos θ 4 U 3 = sin θ 4 cos θ U 2 = sin θ 4 sin θ sin ϕ U 1 = sin θ 4 sin θ cos ϕ
where θ and ϕ are the usual coordinates in three dimensions. By comparison to Equation (52), we have
θ 4 = 2 cot 1 p a θ = cos 1 p 3 p ϕ = tan 1 p 2 p 1

3.6. Orbit in U space

We want to find the trajectory of the particle on the surface of the hypersphere corresponding to the Kepler orbits in configuration space or the displaced circles in momentum space. We assume we have rotated the axes in configuration space so that L is along the 3-axis and A is along the 1-axis, as shown in Figure 1. The equation for the orbit in three-dimensional momentum space is given by Equations (48) or (50). Dividing Equation (48) by p 2 + a 2 immediately gives a parametric equation for the projected orbit in U space:
U 4 = U 2 tan ν .
Because the orbit is in the 1–2 plane in configurations space, U 3 = 0 . The orbit lies in a hyperplane perpendicular to the 2–4 plane that goes through the origin and makes an angle π / 2 ν with the 4-axis, as shown in Figure 6 [104]. The orbit is the intersection of this plane with the hypersphere and it is therefore a great circle. To derive the exact equation for the projected orbit we express p in Cartesian components p 1 and p 2 in Equation (48) and substitute Equation (52) obtaining
U l 2 + U 2 2 2 tan ν U 2 ( l U 4 ) = ( 1 U 4 ) 2 .
To interpret this equation, we consider it in a rotated coordinate system. If we perform a rotation by an amount δ ν about the 1–3 plane ( A 2 δ ν is the generator of this rotation), the equations of transformation may be written [105]
U 2 = U 2 cos δ ν + U 4 sin δ ν , U 3 = U 3 U 4 = U 4 cos δ ν U 2 sin δ ν , U 1 = U 1 .
This transformation is equivalent to making the substitution ν ν + δ ν in the equations relating to the orbit. For example, Equation (56) becomes
U 4 = U 2 tan ( ν + δ ν ) .
We choose δ ν = ν , which means the orbital plane becomes U 4 = 0 . Writing Equation (57) in terms of the primed coordinates, we find
U 1 2 + U 2 2 = 1
which in the original system is
U 1 2 + ( U 2 cos ν + U 4 sin ν ) 2 = 1 .
This is the equation of a great hypercircle ( ν , 0 ) centered at the origin and lying in a hyperplane making an angle π / 2 ν with the 4-axis and an angle π / 2 0 with the 3-axis. If L did not lie along the 3-axis, but, for example, was in the 1–3 plane, at an angle Θ from the 3-axis, then Equation (61) would be modified by the substitution
U 1 U 1 c o s Θ + U 3 sin Θ
which follows, since U i transforms as a three-vector. The corresponding great circle ( ν , Θ ) lies in a hyperplane making an angle π / 2 ν with the 4-axis and π / 2 Θ with the 3-axis.
The motion of the orbiting particle corresponds to a dot moving along the great circle ( ν , 0 or Θ ) with a period T given by the classical period Equation (37). The velocity in configuration space can be expressed in terms of U 4 by using its definition in terms of p 2 Equation (53) or in terms of θ 4 Equation (55). The particle is moving at maximum velocity when θ 4 is a minimum, which occurs at the perihelion when θ 4 = π / 2 ν :
m a x ( p a ) = 1 + e 1 e = r 2 r 1
and at a minimum velocity when θ 4 = π / 2 + ν :
m i n ( p a ) = 1 e 1 + e .
These values of θ 4 correspond to turning points, at which r and p have extreme values. This is apparent when we use Equation (32) for the total energy to show
U 4 = r c r r c .
When r > r c then p 2 < a 2 , so the particle is moving more slowly than the RMS velocity. Applying the virial theorem to any orbit we find p 2 = a 2 so as expected a is the RMS momentum and 1 1 U 4 = 1 = r c / r .
Figure 7 is a picture of a simple device illustrating the stereographic projection of the orbit in p/a-space onto the four-dimensional hypersphere in U-space. We assume that the orbit is in the 1–2 plane and that A lies along the 1-axis, so p 3 = 0 , U 3 = 0 . Because of this trivial dependence on p 3 , we have omitted the 3-axis. The vertical pin or rod represents the unit vector n ^ lying along the 4-axis. The circumference of the larger circle perpendicular to the 4-axis represents the orbit in p / a -space. One can see that it is displaced from the origin along the 2-axis. Centered at the origin, we must imagine a hypersphere of unit radius U ^ 2 = 1 . The stereographic projection U ^ of the vector p / a is obtained by placing the string coming from the top of n directly at the head of the vector p / a . The intersection of the string with the unit hypersphere defines U ^ . As the string is moved along the orbit in p / a -space, it intersects the hypersphere along a great circle shown by the circumference of the unit circle making an angle ν with the 1–2 plane. We can see for example, that at the closest approach, θ 4 is a minimum and U4 is a maximum, U ^ · A = U 1 is a minimum, and p / a is a maximum.

3.7. Classical Time Dependence of Orbital Motion

We can determine the time dependence of the orbital motion by integrating the expression for the angular momentum L = m r 2 d ϕ r / d t
d t = m r 2 L d ϕ r
where r is given by the orbit equation Equation (33) and we are assuming the orbit is in the 1–2 plane. After integrating, we can use the equations relating the momentum space and configuration space variables to obtain the time dependence in p-space and U-space. We obtain
1 r c 2 1 cos ν L m 0 t d t = cos 3 ν 0 ϕ r ( t ) d ϕ r 1 + sin ν cos ϕ r 2 .
The left-hand side of this equation is equal ω c l t , where ω c l is the classical frequency. This follows by substituting Equations (30) and (38)
L = a r c cos ν ω c l = a m r c .
The integral on the right side gives [106]
ω c l t = sin ν cos ν sin ϕ r 1 + sin v cos ϕ r tan 1 cos ν sin ϕ r sin ν + cos ϕ r
which may be simplified as
ω c l t = A L y r c tan 1 A L y r c r
where y = r sin ϕ r .
The relationship between the angle ϕ ϕ p in momentum space and ϕ r in configuration space follows by either differentiating the orbit equation Equation (33) with respect to time and using L = m r 2 ϕ r ˙ or by simultaneously solving the configuration space orbit, the momentum space orbit equation Equation (48) and the energy equation Equation (32). We find
sin ϕ r = p cos ϕ cos ν a cos ϕ r = p sin ϕ cos v a sin ν .
From these equations, the definitions of the U i , Equation (52), and the orbit and energy equations, it follows that for the classical orbit in the 1–2 plane
U 1 = r sin ϕ r r c cos ν U 2 = r r c sin ν + cos ϕ r cos v = U 4 cot ν U 3 = 0 U 4 = r c r r c .
Using these results in Equation (69) gives
ω c l t = U 1 sin ν + tan 1 U 2 U 1 cos ν
which gives the time dependence in U space, and it agrees with the results of [97,107]. We could do rotations to generalize this result [102]. We can also rewrite the inverse tangent as cos 1 U 1 using
U 1 2 + U 2 2 c o s 2 ν = 1 .
If we consider the last two equation for circular orbits with e = sin ν = 0 we obtain tan 1 ( U 2 / U 1 ) = ϕ ( t ) = ω c l t = ϕ r ( t ) + π / 2 , and our familiar circle U 1 2 + U 2 2 = 1 .

Remark on Harmonic Oscillator

We can find a conserved Runge–Lenz vector for the non-relativistic hydrogen atom because the elliptical orbit does not precess, as it does for the relativistic atom. The only central force laws that yield classical elliptical orbits that do not precess are the inverse Kepler force and the linear harmonic oscillator force [108]. Thus, it seems reasonable that one could construct a constant vector similar to A for the oscillator, although the force center for the atom is at a focus and for the oscillator it is at the center of the ellipse. However, it is not readily possible [109]. Instead, one can construct a constant Hermitian second rank tensor T i j :
T i j = 1 m ω 0 p i p j + m ω 0 x i x j .
This constant tensor is analogous to the moment of inertia tensor for rigid body motion. The eigenvectors of the tensor will be constant vectors along the principal axes for the particular orbit being considered. The existence of the conserved tensor leads to the U(3) symmetry algebra of the oscillator. The generators are λ i j a T i j , where the λ ’s are are the usual U ( 3 ) matrices [22]. The spectrum generating algebra is SU(3,1).
In another approach, the Schrodinger equation for the hydrogen atom has been transformed into an equation for a four dimensional harmonic oscillator or two two dimensional harmonic oscillators. This approach which fits well with parabolic coordinates was used especially in the 1980’s to analyze the group structure of the atom and relate it to SU(3) [110,111,112,113,114,115,116,117,118,119]. We will not discuss this approach further.

4. The Hydrogenlike Atom in Quantum Mechanics; Eigenstates of the Inverse of the Coupling Constant

In this section we switch from classical dynamics to quantum mechanics and discuss the group structure and exploit it to determine the bound state energy spectrum directly, as Pauli and followers did almost a century ago [14,16]. In Section 4.3 we introduce a new set of basis states for the hydrogenlike atom, eigenstates of the coupling constant. Using these states allows us to display the symmetries in the most convenient manner and to treat bound and scattering states uniformly.

4.1. The Degeneracy Group SO(4)

The quantum mechanical Hamiltonian is
H = p 2 2 m Z α r = E .
The classical expression for the Runge–Lenz vector needs to be symmetrized to insure the corresponding quantum mechanical operator A is Hermitian:
A = 1 2 m H p × L L × p 2 m Z α r r .
We may verify that A and L = r × p both commute with the Hamiltonian H. The commutation relations of L i and A i are the same as the corresponding classical Poisson bracket relations for bound states:
[ L i , L j ] = i ϵ i j k L k , [ L i , A j ] = i ϵ i j k A k , [ A i , A j ] = i ϵ i j k L k
and form the algebra of O(4) [35]. We can write the commutation relations in a single equation which makes the 0(4) symmetry explicit. If we define
S i j = ϵ i j k L k S i 4 = A i
then
[ S a b , S c d ] = i ( δ a c S b d + δ b d S a c δ a d S b c δ b c S a d ) a , b = 1 , 2 , 3 , 4 .
The Kronecker delta function δ a b acts like a metric tensor.

4.2. Derivation of the Energy Levels

We can obtain the energy levels by determining which representations of the group SO(4) are realized by the degenerate eigenstates of the hydrogenlike atom [12,14,35]. The representations of SO(4) may be characterized by the numerical values of the two Casimir operators for SO(4):
C 1 = L · A C 2 = L 2 + A 2
Once we know the value of C 2 , then the eigenvalues of H follow from the quantum mechanical form of Equation (31), namely
L 2 + A 2 + 1 = ( m Z α ) 2 2 m H .
In order to determine the possible values of C 2 , we factor the 0(4) algebra into two disjoint SU(2) algebras [120], each of which has the same commutation relations as the ordinary angular momentum operators,
N = 1 2 ( L + A ) M = 1 2 ( L A ) .
The commutation relations are
[ M i , N j ] = 0 [ M i , M j ] = i ϵ i j k M k [ N i , N j ] = i ϵ i j k N k .
In analogy with the results for the ordinary angular momentum operators, the Casimir operators are
M 2 = j 1 j 2 + 1 , j 1 = 0 , 1 2 , 1 , N 2 = j 2 j 2 + 1 , j 2 = 0 , 1 2 , 1 ,
The numbers j 1 and j 2 , which may have half-integral values for SU(2) but not O(3), define the ( j 1 , j 2 ) representation of SO(4). From the definitions of A and L in terms of the canonical variables, it follows that C 1 = L · A = 0 which means j 1 = j 2 = j , as in the classical case. For our representations, we find
M 2 = N 2 = 1 4 ( L 2 + A 2 ) = j ( j + 1 ) , j = 0 , 1 2 , 1 , . .
and therefore
L 2 + A 2 + 1 = ( 2 j + 1 ) 2 .
Substituting this result into Equation (82) gives the usual formula for the bound state energy levels of the hydrogen atom:
H = m ( Z α ) 2 2 n 2 = E n
where the principal quantum number n = 2j+1 = 1,2,.. and the prime on H signifies an eigenvalue of the operator H.
Within a subspace of energy E n , the Runge–Lenz vector is
A = 1 a n p × L L × p 2 m Z α r r .
where a n = 2 m E n = m Z α n .
Our considerations of the Casimir operators have shown that the hydrogen atom provides completely symmetrical tensor representations of SO(4), namely, ( j , j ) = ( n 1 2 , n 1 2 ) , n = 1 , 2 , The dimensionality is n 2 , corresponding to the n 2 degenerate states. The appearance of only symmetrical tensor representations ( j 1 = j 2 ) may be traced to L · A vanishing, which is a consequence of the structure of L and A in terms of the dynamical variables for the hydrogenlike atom. For systems other than the hydrogenlike atom, it is not generally possible to find the expression for the energy levels in terms of all the different quantum numbers alone. It worked here, since we could express the Hamiltonian as a function of the Casimir operators that contained all quantum numbers explicitly.
There are a variety of possible basis states. We could choose basis states for the SO(4) representation that reflect the SU(2) decomposition, namely eigenstates of M 2 , N 2 , M 3 and N 3 [120]. Another possibility is to have a basis with eigenstates of the Casimir operator C 2 , and A 3 , and M 3 . This choice fits well with the use of parabolic coordinates [60]. A more physically understandable choice is to choose the common basis states | n l m that are eigenstates of C 2 , L 2 , and L 3 . For this set of basis states, we have
( L 2 + A 2 + 1 ) | n l m = n | n l m L 2 | n l m = l ( l + 1 ) | n l m L 3 | n l m = m | n l m
We can define raising and lowering operators for m:
L ± = L 1 ± i L 2
which obeys the commutation relations
[ L 2 , L ± ] = 0 ] [ L 3 , L ± ] = ± L ± .
Therefore, we can use L ± to change the value of m for the basis states
L ± | n l m = ( l ( l + 1 ) m ( m ± 1 ) | n l m ± 1
for l 1 . We can also use the generators A to change the angular momentum. A general SO(4) transformation can be expressed as a rotation induced by L , followed by a rotation induced by A 3 , followed by another rotation generated by L [121]. Our interest is primarily in changing the angular momentum l, which is most directly done while using A 3 , which commutes with L 3 and C 2 , and so only changes l:
A 3 | n l m = ( n 2 ( l + 1 ) 2 ) ( ( l + 1 ) 2 m 2 ) 4 ( l + 1 ) 2 1 1 2 | n l + 1 m + ( n 2 l 2 ) ( l 2 m 2 ) 4 l 2 1 1 2 | n l 1 m .
for l 1 .

4.3. Relativistic and Semi-Relativistic Spinless Particles in the Coulomb Potential and Klein–Gordon Equation

The Klein–Gordon equation
p 2 ( E ˜ V ) 2 + m 2 ψ ˜ = 0
where E ˜ is the relativistic total energy, may be solved exactly for a Coulomb potential, V = ( Z α ) / r [122]. The energy levels depend on a principal quantum number and on the magnitude of the angular momentum but not its direction. The only degeneracy present is associated with the O(3) symmetry of the Hamilton. For a relativistic scalar particle, there is no degeneracy to be lifted by a Lamb shift.
If we neglect the V 2 term the resulting equation can be written in the form
p 2 2 E ˜ + V E ˜ 2 m 2 2 E ˜ ψ ˜ = 0
This is exactly the same as the nonrelativistic Schrodinger equation with the substitutions
m E ˜ E E ˜ 2 m 2 2 E ˜ .
Thus, we regain the O(4) symmetry of the nonrelativistic hydrogen atom, and can define two conserved vectors, as indicated in Table 1. It is possible to take the “square root” of this approximate Klein-Gordon equation (in the same sense that the Dirac equation is the square root of the Klein–Gordon equation) and get an approximate Dirac equation whose energy eigenvalues are independent of the orbital angular momentum [123].

4.4. Eigenstates of the Inverse Coupling Constant ( Z α ) 1

Solutions to Schrodinger’s equation for a particle of energy E = a 2 2 m in a Coulomb potential
p 2 + a 2 2 m Z α r | a = 0
may be found for certain critical values of the energy E n = a n 2 2 m where a n = m Z α n . The corresponding eigenstates of the Hamiltonian are | n l m which satisfy Equation (95) with a replaced by a n . In addition to the bound states, because there is no upper bound on p 2 in the Hamiltonian, we also have the continuum of scattering states that have E > 0.
Because the quantity that must have discrete values for a solution to exist is actually a m Z α , as noted in Section 2.3, we might ask if eigensolutions to Equation (95) exist for certain critical values of Z α while keeping a and the energy fixed [1]. To investigate such solutions it is convenient to algebraically transform Equation (95):
1 ρ ( a ) ( p 2 + a 2 a ) 1 ρ ( a ) 1 ρ ( a ) ( 2 m Z α a r ) 1 ρ ( a ) ρ ( a ) | a = 0
where
ρ ( a ) = p 2 + a 2 2 a 2 .
Because ρ ( a ) commutes with p 2 , we obtain the eigenvalue equation
a m Z α K ( a ) | a ) = 0
where the totally symmetric and real kernel is
K ( a ) = 2 a 2 p 2 + a 2 1 a r 2 a 2 p 2 + a 2
and
| a ) = ( ρ ( a ) ) 1 2 | a .
As before, solutions to this transformed equation may found for the eigenvalues
K ( a ) = ( a m Z α ) = 1 n .
If we hold Z α constant and let a vary, we obtain the usual spectrum a n = 2 m E n = m Z α n . For a = a n , Equation (99) reduces to the equation for the eigenstates
ρ ( a n ) | n l m = | n l m ; a n ) .
Alternatively, if we hold a constant, then Z α has the spectrum
( Z α ) n = n a m
with the corresponding eigenstates of ( Z α ) 1 being
| n l m ; a ) .
The relationship between the usual energy eigenstates of the energy | n l m and the eigenstates | n l m ; a ) of ( Z α ) 1 is
| n l m = 1 ρ ( a n ) | n l m ; a n ) ,
which requires that both sets of states have the same quantum numbers. Note that the magnitude ( a | K ( a ) | a ) is proportional to 1 / a r n = ( 1 / a n ) ( 1 / n 2 a 0 ) where a 0 is the Bohr radius for the ground state, and so is positive and bounded. The kernel K(a) is real and symmetric in p and r and manifestly Hermitean. Because the kernel K in Equation (98) is bounded, definite, and Hermitian with respect to the eigenstates | n l m ; a n ) [124] the set of normalized eigenstates
| n l m ; a ) n = 0 , 1 , 2 . ; l = 0 , 1 , . . n 2 , n 1 ; m = l , l + 1 , l 1 , l .
where
1 n K ( a ) | n l m ; a ) = 0
is a complete orthonormal basis for the hydrogenlike atom:
n l m ; a | n l m ; a = δ n n δ l l δ m m
n l m | n l m ; a ) ( n l m ; a | = 1 .
There are several important points to notice with regard to these eigenstates of the inverse of the coupling constant:
(1) Because of the boundedness of K, there is no continuum portion in the eigenvalue spectrum of ( Z α ) 1 , the eigenvalues are discrete. Because K is a positive definite Hermitian operator, all eigenvalues are positive, real numbers. This feature leads to a unified treatment of all states of the hydrogenlike atom as opposed to the treatment in terms of energy eigenstates in which we must consider separately the bound states and the continuum of scattering states.
(2) It follows from Equation (104) that the quantum numbers, multiplicities, and degeneracies of these states | n l m ; a ) are precisely the same as those of the usual bound energy eigenstates. For example, there are n 2 eigenstates of ( Z α ) 1 that have the principal quantum number equal to n or ( Z α ) equal to n a m .
(3) A single value of the RMS momentum a or the energy E = a 2 2 m applies to all the states in our complete basis, as opposed to the usual energy eigenstates where each nondegenerate state has a different value of a. We have made this explicit by including a in the notation for the states: | n l m ; a ) . Sometimes we will write the states as | n l m ) , provided the value of a has been specified. This behavior in which a single value of a applies to all states will prove to be very useful. In essence, it permits us to generalize from statements applicable in a subspace of Hilbert space with energy E n or energy parameter a n to the entire Hilbert space.
(4) By a suitable scale change or dilation, we may give the quantity a any positive value we desire. This is affected by the unitary operator
D ( λ ) = e i 1 2 ( p · r + r · p ) λ
which transforms the canonical variables
D ( λ ) p D 1 ( λ ) = e λ p D ( λ ) r D 1 ( λ ) = e λ r
and the kernel K(a)
D ( λ ) K ( a ) D 1 ( λ ) = K ( a e λ ) .
and the eigenvalue equation
1 n K ( a e λ ) D ( λ ) | n l m ; a ) = 0 .
Therefore the states transform as
D ( λ ) | n l m ; a ) = | n l m ; a e λ ) .
These transformed states form a new basis corresponding to the new value e λ a of the RMS momentum.
The relationship between the energy eigenstates and the ( Z α ) 1 eigenstates can be written while using the dilation operator:
| n l m = 1 ρ a n D λ n | n m ; a ) where e λ n = a n a .
The usual energy eigenstates | n l m are obtained from the eigenstates of ( Z α ) 1 by first performing a scale change to ensure that the energy parameter a has the value a n and then multiplying by a factor ρ 1 / 2 . The need for the scale change is apparent from dimensional considerations: from the ( Z α ) 1 eigenvalue equation, we see that the eigenfunctions are functions of p / a or a r while from the energy eigenvalue equation the eigenfunctions are functions of p / a n and a n r . The factor ρ 1 / 2 was required in order to convert Schrodinger’s equation to one involving a bounded Hermitean operator.
Using the eigenstates of ( Z α ) 1 as our basis allows for us to analyze the mathematical and physical structure of the hydrogenlike atom in the easiest and clearest way.

4.5. Another Set of Eigenstates of ( Z α ) 1

We may transform Schrodinger’s equation Equation (95) to an eigenvalue equation for ( Z α ) 1 that differs from Equation (96) by similar methods:
1 n K 1 ( a ) | n l m ; a ) = 0
where
K 1 ( a ) = 1 a r 2 a 2 p 2 + a 2 1 a r ρ ( a ) = n / a r ρ ( a n ) | n l m = | n l m ; a n )
This kernel, like K(a), is a bounded, positive definite, Hermitian operator so the eigenstates form a complete basis. The relationship of these basis states to the energy eigenstates is the same as that of the previously discussed eigenstates of ( Z α ) 1 Equation (114) but with ρ ( a ) = n / a r . The n insures that the two sets of eigenstates have consistent normalization, which may be checked by means of the virial theorem. The n cancels out when similarity transforming from the basis of energy eigenstates to the basis of ( Z α ) 1 . Note that, classically, both kernels equal 1 / a r c .
The first set of basis states of ( Z α ) 1 with ρ ( a ) = p 2 + a 2 2 a 2 is more convenient to use when working in momentum space and the second set with ρ ( a ) = n / a r is more convenient in configuration space.
Other researchers have used other approaches to secure a bounded kernel for the Schrodinger hydrogen atom, for example, by multiplying the equation from the left by r to regularize it [82]. However, the methods used have not symmetrized the kernels to make them Hermitian, nor are all the generators of the corresponding groups Hermitian, and they have to redefine the inner product [38,82].

4.6. Transformation of A and L to the New Basis States

We must transform A as given in Equation (89) and L = r × p when we change our basis states from eigenstates of the energy to eigenstates of the inverse coupling constant. The correct transformation may be derived by requiring that the transformed generators produce the same linear combination of new states as the original generators produced of the old states. Thus, because
A | n l m = l , m | n l m A l m l m
where the coefficients A l m l m are the matrix elements of A , we require that the transformed generator a satisfies the equation
a | n l m ) = l m | n l m ) A l m l m .
In other words, since the Runge–Lenz vector A is a symmetry operator of the original energy eigenstates, a will be a symmetry operator of the new states with precisely the same properties and matrix elements. Because A is Hermitian, a is Hermitian.
To obtain a differential expression for a acting on the new states, we need to transform the generator using Equation (114):
a = D 1 ( λ n ) ρ ( a n ) A 1 ρ ( a n ) D ( λ n ) .
The effect of the scale change on the quantity in large parenthesis is to replace a n everywhere by a. By explicit calculation, we find
a = 1 2 a r p 2 + p 2 r 2 r · p p p p · r a r 2
for ρ ( a ) = ( p 2 + a 2 ) / 2 a 2 . And we obtain
a = 1 2 a r p 2 + p 2 r 2 r · p p p p · r r 4 r 2 a r 2
for ρ ( a ) = n / a r .
Both of these expressions for a are manifestly Hermitian. In addition, since there is no dependence on the principal quantum number these expressions are valid in the entire Hilbert space, and not just in a subspace spanned by the degenerate states, as was the case when we used the energy eigenstates as a basis (Equation (89)).
The angular momentum operator is invariant under scale changes and it commutes with scalar operators. Therefore L is invariant under the similarity transformation 1 ρ ( a n ) D ( λ n ) and the expression for the angular momentum operator with respect to the eigenstates of ( Z α ) 1 is the same as the expression with respect to the eigenstates of the energy.

4.7. The U | Representation

The U coordinates provide the natural representation for the investigation of the symmetries of the hydrogenlike atom in quantum mechanics, as in classical mechanics [125]. Therefore, we briefly consider the relevant features of this representation and, in particular, its relationship to the momentum representation. The eigenstate < U | of U b , b = 1,2,3,4, is defined by
U | U b = U b U |
These states are complete on the unit hypersphere in four dimensions:
| U U | d 3 Ω = 1
where Ω refers to the angles ( θ 4 , θ , ϕ ) defined in Equation (54). The U variables are defined in terms of the momentum variables and the quantity a in Equation (52). Therefore, the momentum and the U operators commute
[ p i , U b ] = 0
and the state U | is proportional to a momentum eigenstate p | :
U | = p | J ( p )
where the momentum eigenstate is defined by p | p = p p | and
d 3 p | p p | = 1 .
The function J ( p ) may be determined by equating the completeness conditions and substituting Equation (125):
1 = d 3 p | p p | = d 3 Ω | U U | = d 3 Ω J ( p ) | p p |
which leads to the identification of the differential quantities
d 3 p = d 3 Ω J ( p )
demonstrating that J ( p ) is the Jacobian of the transformation from the p- to the U- space. Noting that on the unit sphere
U 4 2 = 1 U i U i
we can compute the Jacobian
J ( p ) = p 2 + a 2 2 a 3 .
Therefore from Equation (125) we have the important result
U | = p | p 2 + a 2 2 a 3 / 2 .
We can use this result to compute the action of r on U | in terms of the differential operators. Using the equation
p | r = i p p |
we obtain
U | r = i p 3 i p p 2 + a 2 U | .

Action of a and L on U |

Using Equation (133) for the action of r on U | and using the expression Equation (120) for a since U | is proportional to p | , we immediately find that when acting on U | , a has the differential representation
a = i 2 a ( p 2 a 2 ) p 2 p p · p
where
U | a = a U | .
We can also write a in terms of the U variables by using the relationship Equation (52) between the p and U variables:
a = U 4 i U U i / 4
where the spatial part of the four vector U is U = ( U 1 , U 2 , U 3 ) and U 4 is the fourth component. This is the differential representation of a rotation operator mixing the spatial and the fourth components of U a . When acting on the state U | , clearly e i a · ν generates a four-dimensional rotation that produces a new eigenstate U | . To explicitly derive the form of the finite transformation, we compute
[ a j , U j ] = i U 4 δ i j [ a j , U 4 ] = i U i .
For a finite transformation a i a · n ν with n 2 = 1 , we have
U = e i a · n ν U e i a · n ν = U n n · U + n n · U cos ν n U 4 sin ν
and
U 4 = e i a · n ν U 4 e i a · n ν = U 4 cos ν + n · U sin ν
These equations of transformation are like those for a Lorentz transformation of a four-vector (r, it). We can illustrate the equations for e i a 2 ν (cf Equation (58)), which mixes the two and four components of U :
U 1 = U 1 U 3 = U 3 U 2 = U 2 cos ν U 4 sin ν U 4 = U 2 sin ν + U 4 cos ν .
When L acts on U | , it has the differential representation
L = U × i U
This result follows directly, since U i equals p i times a factor that is a scalar under rotations in three dimensions. When e i L · ω acts on U | it produces a new state U | , where the spatial components of U have been rotated to produce U .
In summary, we see that U is a four-vector under rotations generated by a and L . Therefore the states U | provide a vector representation of the group of rotations in four dimensions SO(4), with the generators a and L .

5. Wave Functions for the Hydrogenlike Atom

In this section, we analyze the wave functions of the hydrogenlike atom, working primarily in the U | representation and using eigenstates of the inverse of the coupling constant ( Z α ) 1 for the basis states. In this representation the wave functions are spherical harmonics in four dimensions. We derive the relationship of the usual energy eigenfunctions in momentum space to the spherical harmonics and discuss the classical limits in momentum and configuration space.

5.1. Transformation Properties of the Wave Functions under the Symmetry Operations

We can show that the wave functions Y n l m ( U ) in the U | representation with respect to the eigenstates of ( Z α ) 1
Y n l m ( U ) U | n l m )
transform as four-dimensional spherical harmonics under the four-dimensional rotations generated by the Runge-Lenz vector a and the angular momentum L . We note that the quantity a is implicit in both the bra and the ket in Equation (142). For our basis states we employ the set of ( Z α ) 1 eigenstates | n l m ) of the inverse coupling constant that are convenient for momentum space calculations ( ρ = p 2 + a 2 2 a 2 ). We choose these states rather than those convenient for configuration space calculations, because the U | eigenstates are proportional to the p | eigenstates.
If we transform our system by the unitary operator e i θ where θ = L · ω + a · ν , then the wave function in the new system is
Y n l m ( U ) = U | e i θ | n l m ) .
There are two ways in which we may interpret this transformation, corresponding to what have been called the active and the passive interpretations. In the passive interpretation we let e i θ act on the coordinate eigenstate U | . As we have seen in Section 3.6, this produces a new eigenstate U | , where the four-vector U is obtained by a four-dimensional rotation of U (Equations (138) and (139)). Thus, we have
Y n l m ( U ) = U | n l m ) = Y n l m ( U ) .
In the active interpretation, we let e i θ act on the basis state | n l m ) . Because L and a are symmetry operators of the system, transforming degenerate states into each other, it follows that e i θ | n l m ) must be a linear combination of states with principal quantum number equal to n. Therefore, we have
Y n l m ( U ) = l m U | R n l m n l m | n l m ) = l m R n l m n l m Y n l m ( U ) .
The wave functions for degenerate states with a given n transform irreducibly among themselves under the four-dimensional rotations, forming a basis for an irreducible representation of SO(4) of dimensions n 2 . Equating the results of the two different interpretations gives
Y n l m ( U ) = l m R n l m n l m Y n l m ( U ) .
The transformation properties Equation (146) of Y n l m are precisely analogous to those of the three-dimensional spherical harmonic functions. It follows that the Y n l m are four-dimensional spherical harmonics [35,124].

5.2. Differential Equation for the Four Dimensional Spherical Harmonics Y n l m ( U )

The differential equation for the Y n l m ( U ) may be obtained from the equation
( L 2 + a 2 ) Y n l m ( U ) = ( n 2 1 ) Y n l m ( U )
which follows from C 2 = n 2 1 and the definition of C 2 , Equation (81). Substituting in the differential expressions Equations (134) and (141) for a and L we find that L 2 + a 2 equals U 2 ( U · U ) 2 , which is the angular part of the Laplacian operator in four dimensions (cf in three dimensions, L 2 / r 2 = p 2 p r 2 ) . Thus, Equation (147) is the differential equation for four-dimensional spherical harmonics with the degree of homogeneity equal to n 1 , which means n 2 such functions exist, in agreement with the know degree of degeneracy.

5.3. Energy Eigenfunctions in Momentum Space

We want to determine the relationship between the usual energy eigenfunctions in momentum space ψ n l m ( p ) p | n l m (with a = a n ) and the four-dimensional spherical harmonic eigenfunctions Y n l m ( U ; a ) = U | n l m ; a ) .
We choose the RMS momentum a to have the value a n . If we use the expression Equation (131) for U | in terms of p |
U | = p | p 2 + a n 2 2 a n 3 / 2
and the expression Equation (104) for the eigenstates of ( Z α ) 1 in terms of the energy eigenstates
| n l m ; a n ) = p 2 + a n 2 2 a n 2 | n l m
we find the desired result
Y n l m ( U ; a n ) = p 2 + a n 2 2 a n 2 1 a n ψ n l m ( p ) .
The usual method of deriving this relationship between the wave function in momentum space and the corresponding spherical harmonics in four dimensions involves transforming the Schrodinger wave equation to an integral equation in momentum space [18,35]. As in the classical case, to do this we first replace p by p / a and perform a stereographic projection from the hyperplane corresponding to the three- dimensional momentum space to a unit hypersphere in a four-dimensional space. The resulting integral equation manifests a four-dimensional invariance. When the wave functions are normalized as in Equation (150), the solutions are spherical harmonics in four dimensions. As another alternative to this procedure, we can Fourier transform the configuration space wave functions directly [126].

5.4. Explicit Form for the Spherical Harmonics

The spherical harmonics in four dimensions can be expressed as [127]:
Y n l m ( Ω ) = N 1 ( n , l ) ( sin θ 4 ) l C n 1 l l + 1 ( cos θ 4 ) · N 2 ( l , m ) ( sin θ ) m C l m m + 1 / 2 ( cos θ ) e i m ϕ 2 π .
The factor in brackets is equal to Y l m ( θ , ϕ ) , the usual spherical harmonic in three-dimensions [127]. The Gegenbauer polynomials C n λ of degree n and order λ are defined in terms of a generating function:
1 ( 1 2 t x + t 2 ) λ = n = 0 t n C n λ .
N 1 ( n , l ) and N 2 ( l , m ) are chosen to normalize the Y n l m on the surface of the unit sphere:
| Y n l m ( Ω ) | 2 d 3 Ω U = 1
where d 3 Ω U = sin 2 θ 4 sin θ d θ d ϕ . We find
N 1 ( n , l ) = 2 2 l + 1 π n ( n l 1 ) ! ( l ! ) 2 ( n + l ) !
N 2 ( l , M ) = 2 2 m π ( l + 1 2 ) ( l m ) ! ( l + m ) ! [ Γ ( m + 1 2 ) ] 2 .
In the next section, we discuss the asymptotic behavior of Y n l m for large quantum numbers and compare it to the classical results of Section 3.

5.5. Wave Functions in the Classical Limit

5.5.1. Rydberg Atoms

Advances in quantum optics, such as the development of ultra short laser pulses, microwave spectroscopy, and atom inteferometery, have opened new possibilities for experiments with atoms and Rydberg states, meaning hydrogenlike atoms in states with very large principal quantum numbers and correspondingly large diameter electron orbits. The pulsed electromagnetic fields can be used to modify the behavior of the orbital electrons. Semi-classical electron wave packets in hydrogenlike atoms were first generated in 1988 by ultrashort laser pulses, and today are often generated by unipolar teraherz pulses [128,129,130]. Over the last few decades, there has been interest in the classical limit of the hydrogenlike atom for n very large, Rydberg states, for a number of reasons [131]: 1. Rydberg states are at the border between bound states and the continuum, and any process which leads to excited bound states, ions or free electrons usually leads to the production of Rydberg states. This includes, for example, photo-ionization or the application of microwave fields. The very large cross section for scattering is unique. 2. Rydberg states can be used to model atoms with a higher atomic number that have an excited valence electron that orbits beyond the core. 3. In Rydberg states, the application of electric and magnetic fields breaks the symmetry of the atom and allows for the study of different phenomena, including the transition from classical to quantum chaos [132]. 4. Rydberg atoms can be used to study coherent transient excitation and relaxation, for example, the response to short laser pulses creating coherent quantum wave packets that behave like a classical particle.
The square of the wave function for a given quantum state gives a probability distribution for the electron that is independent of time. If we want to describe the movement of an electron in a semiclassical state, with a large radius, going around the nucleus with a classical time dependence, then we need to form a wave packet. The wave packet is built as a superposition of many wave functions with a band of principal quantum numbers.
A variety of theoretical methods have been used to derive expressions for the hydrogen atom wave functions and wave packets for highly excited states. There is general agreement on the wave functions for large n, and that the wave functions display the expected classical behavior, elliptical orbits in configurations space, and great circles in the four dimensional momentum space [133,134,135,136].
Researchers have proposed a variety of wave packets to describe Rydberg states. There are general similarities in the wave packets that describe electrons going in circular or elliptical orbits with a classical time dependence for some characteristic number of orbits, and it is maintained that the quantum mechanical wave packets provide results that agree with the classical results [129,130,133,134,135,136,137,138,139]. Most of the approaches exploit the SO(4) or SO(4,2) symmetry of the hydrogen atom which is used to rotate a circular orbital to an elliptical orbit. The starting orbital is often taken as a coherent state, which is usually considered a classical like state with minimum uncertainty. The most familiar example of a coherent state is for a one dimensional harmonic oscillator characterized by creation and annihilation operators a and a. The coherent state | α is a superposition of energy eigenstates that is an eigenstate of a where a | α = α | α for a complex α . This coherent state will execute harmonic motion like a classical particle [140]. To obtain a coherent state for the hydrogenlike atom, eigenstates of the operator than lowers the principal quantum number n (which will be discussed in Section 7.4) have been used [141], as well as lowering operators based on the equivalence of the four dimensional harmonic oscillator representation of the hydrogen atom [134,135,142].
In either case, this coherent eigenstate is characterized by a complex eigenvalue, which needs to be specified. Several constraints have been used to obtain the classical wave packet that presumably obeys Kepler’s Laws, such as requiring that the orbit lie in a plane so z = 0 for the orbital, or that r r c l a s s i c a l be a minimum, or that some minimum uncertainty relationship is obtained. In addition, there are issues regarding the approximations used, in particular, those that relate to time. For times characteristic of the classical hydrogen atom, the wavepackets act like a classical system. For longer times, the wavepacket spreads in the azimuthal direction and after some number of classical revolutions of order 10 to 100 the spread is 2 π , so the electron is uniformly spread over the entire orbit. The spread arises because the component wave functions forming the wavepacket have different momenta. In two derivations, still longer times were considered, and recoherence was predicted to occur after about n / 3 (where n is the approximate principal quantum number) revolutions, although there is some difference in the predicted amount of recoherence [131,136]. Because of the conservation of L and A , the spread of the wave packets is inhibited, except in the azimuthal direction.
This is a system with very interesting physics. For example, one can view the transition from a well defined wave packet representing the electron to a 2 π spread and back again as a dynamical illustration of the wave-particle duality.
Brown took a different approach to develop a wavepacket for a circular orbit [139]. He first developed the asymptotic wave function for large n and then optimized the coefficients in a Gaussian superposition to minimize the spread in ϕ , obtaining a predicted characteristic decoherence time of about 10 minutes, considerable longer than any other predicted decoherence time.
Other authors have explored the problem from the perspective of classical physics and the correspondence principle [133,137,143,144]. The results from the different methods are similar with the basic conclusion that the wave functions are peaked on the corresponding classical Kepler trajectories: “atomic elliptic states sew the wave flesh on the classical bones” [129].
With the variety of experimental methods used to generate Rydberg states, a variety of Rydberg wave packets are created, and it is not clear which theoretical model, if any, is preferred [131]. We take a very simple approach to forming a wavepacket and simply use a Gaussian weight for the different frequency components. This does not give an intentionally optimized wave packet, but it is a much simpler approach and the result has all of the expected classical behavior that is very similar to that obtained from much more complicated derivations. We start with a circular orbit and then do a SO(4) rotation to secure an elliptical orbit. We show that it has the classical period of rotation.

5.5.2. Wave Functions in the Semi-Classical Limit

We need to derive the semiclassical limit of the wave functions that correspond to circular orbits in configuration space. For this case, sin ν , which we interpret as the expectation value of the eccentricity, vanishes. We derive expressions for the wave functions in momentum space and then form a wavepacket. To obtain corresponding expressions for elliptical orbits, we perform a rotation by e i a · ν , which does not alter the energy but changes the eccentricity and the angular momentum.
Case 1: Circular orbits, sin ν = e = 0
We derive the asymptotic form of Y n l m for large quantum numbers, where for simplicity we choose the quantum numbers n 1 = l = m corresponding to a circular orbit in the 1–2 plane. From Equation (151) we see we encounter Gegenbauer polynomials of the form C 0 λ , which, by Equation (152), are unity. For a very large l, sin l θ will have a very strong peak at θ = π / 2 , so we make the expansion [103]
sin θ = sin π 2 + θ π 2 = 1 1 2 θ π 2 2 + e ( 1 / 2 ) ( θ π / 2 ) 2
to obtain
sin l θ e ( 1 / 2 ) l ( θ π / 2 ) 2 .
The asymptotic forms for N 1 and N 2 can be computed using the properties of Γ functions:
Γ ( 2 z ) = 1 2 π 2 2 z 1 2 Γ ( z ) Γ z 1 2 lim z Γ ( a z + b ) 2 π e a z ( a z ) a z + b 1 2 .
We finally obtain [145]
Y n , n 1 , n 1 ( Ω ) = n 2 π 2 e 1 2 n θ 4 π 2 2 · e 1 2 n θ π 2 2 e i ( n 1 ) ϕ ,
which gives the probability density
Y n , n 1 , n 1 ( Ω ) 2 = n 2 π 2 e n θ 4 π 2 2 · e n θ π 2 2 .
We have Gaussian probability distributions in θ 4 and θ about the value π / 2 . The distributions are quite narrow with widths Δ θ 4 Δ θ 1 / n and the spherical harmonic essentially describes a circle ( θ 4 = θ = π / 2 ) on the unit sphere in the 1–2 plane. As n becomes very large, U 4 = cos θ 4 ( r r c ) / r (Equations (54) and (72)) and U 3 = sin θ 4 cos θ , which is proportional to p 3 , both go to zero as 1 / n . The distribution approaches the great circle U 1 2 + U 2 2 = 1 that we found in Section 3.6 for a classical particle moving in a circular orbit in the 1–2 plane in configuration space. Note that this state is a stationary state with a constant probability density. To get the classical time dependence we need to form a wavepacket.
Forming a Wavepacket
We form a time dependent wavepacket for circular orbits by superposing circular energy eigenstates:
χ ( Ω , t ) = n e i t E n Y n , n 1 , n 1 A n N
where A n N is an amplitude peaked about n = N > > 1 . For n > > 1 we expand E n about E N :
E n = E N + E n N s + 2 E n 2 N s 2 +
where s = n N . From the equation for the energy levels, E = m ( Z α ) 2 / ( 2 n 2 ) we compute
E n n | N = m ( Z α ) 2 N 3 = 8 E N 3 m ( Z α ) 2 .
In agreement with the Bohr Correspondence Principle, the right-hand side of this equation is just the classical frequency ω c l as given in Equation (38). For the second order derivative, we have
2 E n 2 | N = 3 N ω c l β
which gives
χ ( Ω , t ) = e i t E N e i ϕ ( N 1 ) s = N + 1 e i ω c l s t β 2 s 2 t s ϕ · A s | Y N + s , N + s 1 , N + s 1 | .
We choose a Gaussian form for A s
A s = 1 2 π N e s 2 / ( 2 N ) .
Brown used A s = C e s 2 3 ω c l t / N which minimizes the diffusion in ϕ at time t [139]. Since | Y N + s , N + s 1 , N + s 1 | varies slowly with s for N > > 1 we can take it outside the summation in Equation (165). We now replace the sum by an integral over s. Because A s is peaked about N we can integrate from s = to s = + . We perform the integral by completing the square in the usual way. The final result for the probability amplitude for a circular orbital wave packet is
| x ( Ω , t ) | 2 = Y N , N 1 , N 1 2 1 + β 2 t 2 N 2 1 · exp ϕ ω c l t 2 N 1 + ( β t N ) 2 .
This represents a Gaussian distribution in ϕ that is centered about the classical value ϕ = ω c l t , meaning that the wavepacket is traveling in the classical trajectory with the classical time dependence. The width of the ϕ distribution is
Δ ϕ = ( N ) ( 1 / 2 ) ( 1 + β 2 t 2 N 2 ) 1 / 2 = ( N ) ( 1 / 2 ) ( 1 + 9 ω c l 2 t 2 ) 1 / 2 .
The distribution in ϕ at t = 0 is very narrow, proportional to 1 / N , but after several orbits Δ ϕ is increasing linearly with time.
The distributions in θ 4 and θ are Gaussian and centered about π / 2 in each case as for the circular wave function (cf. Equation (160)) with widths equal to ( N ) 1 / 2 . The spreading of these distributions in time is inhibited because of the conservation of angular momentum and energy. The detailed behavior of the widths depend on our use of the Gaussian distribution. Other distributions will give different widths, although the general behavior is expected to be similar.
As a numerical example, consider a hydrogen atom that is in the semiclassical region when the orbital diameter is about 1 cm. The corresponding principal quantum number is about 10 4 , the mean velocity is about 2.2 × 10 4 cm/s and the period about 1.5 × 10 4 s. After about 34 revolutions or 5 × 10 3 s, the spread in ϕ is about 2 π , meaning the electron is spread uniformly about the entire circular orbit. This characteristic spreading time can be compared to 1.6 × 10 3 s for a fully optimized wave packets formed from coherent SO(4,2) states [136,146]. In order to make predictions about significantly longer times, we would need to retain more terms in the expansion Equation (162) of E n .
Case 2: Elliptical orbits sin ν = e 0
We can obtain the classical limit of the wave function for elliptical orbits by first writing our asymptotic form Equation (159) for Y n , n 1 , n 1 in terms of the U variables instead of the angular variables by using definitions Equation (52), and setting a = a n . Retaining only the lowest order terms in ( θ 4 π / 2 ) and ( θ π / 2 ) , we find
Y n , n 1 , n 1 ( U ) = n 2 π 2 1 2 e i ( n 1 ) tan 1 U 2 U 1 · e 1 2 n U 4 2 e 1 2 n U 3 2
For large n, this represents a circular orbit in the 1–2 plane. We now perform a rotation by A 2 ν . which will change the eccentricity to sin ν , and change the angular momentum, but will not change the energy or the orbital plane. Using Equation (140) to express the old coordinates in terms of the new coordinates, we find to lowest order
Y n , n 1 , n 1 ( U ) = n 2 π 2 1 / 2 e i ( n 1 ) tan 1 U 2 U 1 cos v · e 1 2 n U 2 sin ν U 4 cos ν 2 · e 1 2 n U 3 2 .
In Section 3.6, we found that the vanishing of the term in braces 0 = U 2 sin ν U 4 cos ν specifies the classical great hypercircle orbit (Equation (56)) corresponding to an ellipse in configuration space with eccentricity e = sin ν and lying in the 1–2 plane. The probability density | Y n , n 1 , n 1 ( U ) | 2 vanishes except within a hypertorus with a narrow cross section of radius approximately 1 n , which is centered about the classical distribution. Because the width 1 n of the distribution is constant in U space, it will not be constant when projected onto p space.
In terms of the original momentum space variables, the asymptotic spherical harmonic is
Y n , n 1 , n 1 ( p ) = n 2 π 2 1 2 e i ( n 1 ) tan 1 p 2 p 1 cos ν · exp n 2 p 1 2 + p 2 a tan ν 2 a 2 cos 2 ν 2 cos ν p 2 + a 2 2 · exp n 2 2 p 3 a p 2 + a 2 2 | a = a n .
The expression in brackets corresponds to the momentum space classical orbit equation we found previously (Equation (49)). As we expect, p 3 is Gaussian about zero since the classical orbit is in the 1–2 plane. We can simplify the expressions for the widths by observing that to lowest order we can use Equation (48), which implies p 2 + a 2 = 2 a 2 + 2 a p 2 tan ν in the exponentials. The widths of both distributions therefore increase linearly with p 2 . We also note that, since classically there exists a one-to-one correspondence between each point of the trajectory in momentum space and each point in configuration space, we may interpret the widths of the distributions using Equation (32) p 2 + a 2 a 2 = 2 r c r . Accordingly, the widths increase as the momentum increases or as the distance to the force center decreases (Figure 8).
Forming a Wave packet for Elliptical Motion
We may form a time dependent wavepacket superposing the wave functions of Equation (170). Care must be taken to include the first order dependence (through a n ) of tan 1 ( U 2 / U 1 cos ν ) on the principal quantum number when integrating over the Gaussian weight function. The result for the probability density is the same as before (Equation (167)), except | Y N , N 1 , N 1 | 2 (given in Equation (170)) replaces | Y N , N 1 , N 1 | 2 and
ω c l t = tan 1 U 2 U 1 cos ν + sin ν U 1
where replaces ω c l t = ϕ . The result is exactly the same as the classical time dependence Equation (73). The spreading of the wave packet will be controlled by the same factor as for the circular wave packet.
Remark on the Semiclassical Limit in Configuration Space
The time dependent quantum mechanical probability density follows the classical trajectory in momentum space meaning that the probability is greatest at the classical location of the particle in momentum space. Because the configuration space wave function is the Fourier transform of the momentum space wave function, the classical limit must also be obtained in configuration space. That this limit is obtained is made explicit by observing that the momentum space probability density is large when
( U 2 sin ν U 4 cos ν ) 2 0 .
However from Section 3 we can show that
U 2 sin ν U 4 cos ν = cos ν r r c l a s s i c a l r c l a s s i c a l
where r c l a s s i c a l is given by the classical orbit Equation (33). Accordingly, we see that the configurations space probability will be large when
r r c l a s s i c a l r c l a s s i c a l 2 0 .

5.6. Quantized Semiclassical Orbits

It is convenient at times to have a semiclassical model for the orbitals of the hydrogenlike atom. Historically this was first done by Pauling and Wilson [147]. We can obtain a model by interpreting the classical formulae for the geometrical properties of the orbits as corresponding to the expectation values of the appropriate quantum mechanical expressions. Thus, when the energy E = a 2 / 2 m appears in a classical formula, we employ the expression for a for the quantized energy levels a = 1 n r 0 where r 0 = ( m Z α ) 1 , which is the radius 0.53 Angstrom of the ground state. Similarly if L 2 appears in a classical formula, we substitute l ( l + 1 ) , where l is quantized l = 0 , 1 , 2 . . n 2 , n 1 ; and m, the component of L along the 3-axis is quantized: m= l , l + 1 , l + 2 , . . l .
Orbits in Configuration Space
Recalling Equation (28) r c = m Z α a 2 , we see a r c = n , which gives a semimajor axis of length r c = n 2 r 0 , where r 0 is the radius for the circular orbit of the ground state and r c is for a circular orbit for a state with principal quantum number n. The equations for the magnitude of L and A are
L = r c a cos v = n cos v = l ( l + 1 )
A = r c a sin v = n sin v = n 2 l ( l + 1 )
This gives an eccentricity sin ν equal to
e = sin v = 1 l ( l + 1 ) n 2
and a semi-major axis equal to
b = r c sin ν = n l ( l + 1 ) .
Note that the expression for e is limited in its meaning. For an s state, it always gives e = 1, and for states with l = n 1 it give e = 1 / n , not the classically expected 0 for a circular orbit.
Orbits in U-Space
The corresponding great hypercircle orbits ( ν , Θ ) in U-space are described by giving the quantized angle ν , between the three-dimensional hyperplane of the orbit and the 4-axis, and the quantized angle Θ , between the hyperplane of the orbit and the 3-axis:
cos ν = l ( l + 1 ) n 2
cos Θ = m 2 l ( l + 1 ) .
Note the similarity in these two equations, suggesting that m relates to l the same way that l relates to n, which suggests a generalization of the usual vector model of the atom which only describes the precession of L about the z axis.
The results for orbits in configuration and momentum space illustrate some interesting features:
1. The equation a r = n illustrates that the characteristic dimensions of an orbit in configuration space and the corresponding orbit in momentum space are inversely proportional, as expected, since they are related by a Fourier transform, consistent with the Heisenberg Uncertainty Principle.
2. If l = 0 , then no classical state exists. The orbit in configuration space degenerates into a line passing through the origin while the corresponding circular orbit in momentum space attains an infinite radius and an infinite displacement from the origin. Although this seems peculiar from the pure classical viewpoint, quantum mechanically it follows, since for S states there is a nonvanishing probability of finding the electron within the nucleus.
In order to interpret these statements about quantized semiclassical elliptical orbits we observe that for quantum mechanical state of the hydrogenlike atom with definite n , l , m , the probability density is (1) independent of ϕ r or ϕ p and (2) it does not confine the electron to some orbital plane. Because the quantum mechanical distribution for such a state specifies no preferred direction in the 1–2 plane, we must imagine this distribution as corresponding in some way to an average over all possible orientations of the semiclassical elliptical orbit. This interpretation is supported by the fact that the region within which the quantum mechanical radial distribution function differs largely from zero is included between the values of r corresponding to the semiclassical turning points r c ( 1 ± sin ν ) .

5.7. Four-Dimensional Vector Model of the Atom

In configuration space or momentum space, the angle between the classical plane of the orbit and the 3-axis is Θ , which is usually interpreted in terms of the vector model of the atom in which we imagine L to be a vector of magnitude l ( l + 1 ) precessing about the 3-axis, with m as the component along the 3-axis. This precession may be linked to the ϕ r independence of the probability and the absence of an orbital plane as mentioned at the end of the preceding section. The precession constitutes a classical mechanism which yields the desired average over all possible orientations of the semiclassical elliptical orbit. Because the angle Θ is restricted to have only certain discrete values one can say that there is a quantization of space.
The expression for cos ν = l ( l + 1 ) / n 2 is quite analogous to that for Θ and so suggests a generalization of the vector model of the atom to four dimensions. The projection of the four-dimensional vector model onto the physical three-dimensional subspace must give the usual vector model. We can achieve this by imagining that a four-dimensional vector of length n, where n is the principal quantum number, is precessing in such a way that its three and four components are constants, while the one and two components vary periodically. The projection onto the 1–2–3 hyperplane is a vector of constant magnitude l ( l + 1 ) precessing about the 3-axis. The component along the 3-axis is m. The component along the 4-axis is A = n 2 l ( l + 1 ) the magnitude of the vector A . The vectors L and A are perpendicular to each other. Thus, the precessing n vector makes a constant angle Θ with the 3-axis and a constant angle π / 2 ν with the 4-axis. Because both angles are restricted to certain values, we may say that we have a quantization of four-dimensional space.

6. The Spectrum Generating Group SO(4,1) for the Hydrogenlike Atom

We consider the Schrodinger hydrogen atom and its unitary “noninvariance” or spectrum generating operators e i D i β i where D i is a generator and β i is a real parameter, using eigenstates of ( Z α ) 1 for our basis of our representation. These operators transform an eigenstate of the kernel K with a definite value of the coupling constant (or principal quantum number) into a linear combination of eigenstates with different values of the coupling constant (or different principal quantum numbers), and different l and m. Unlike the invariance generators L and A , the noninvariance generators clearly do not generally commute with the kernel K, [ D i , K ] 0 , so they change the principal quantum number.
The set of all invariance and noninvariance operators forms a group with which we may generate all eigenstates in our complete set from a given eigenstate. We show that this group, called the Spectrum Generating Group of the hydrogenlike atom, is SO(4,1), the group of orthogonal transformations in a 5-dimensional space with a metric g A B = ( 1 , 1 , 1 , 1 , 1 ) , where A , B = 0 , 1 , 2 , 3 , 4 . The complete set of eigenstates of ( Z α ) 1 for the hydrogenlike atom forms a unitary, irreducible, infinite-dimensional representation of SO(4,1) which, we shall find, can be decomposed into an infinite sum of irreducible representations of SO(4), each corresponding to the degeneracy group for a particular principal quantum number. A unitary representation means all generators are unitary operators. An irreducible representation does not contain lower dimensional representations of the same group. In Section 6.3, we discuss the isomorphism between the spectrum generating group SO(4,l) and the group of conformal transformations in momentum space. An isomorphism means the groups have the same structure and can be mapped into each other.

6.1. Motivation for Introducing the Spectrum Generating Group Group SO(4,1)

We have examined the group structure for the degenerate eigenstates of ( Z α ) 1 for the Schrodinger hydrogenlike atom: the n 2 degenerate states form an irreducible representation of SO(4). The next question we might ask is: Do all or some of the states with different principal quantum numbers form an irreducible representation of some larger group which is reducible into SO(4) subgroups? If such a group exists then it clearly is not an invariance group of the kernel K (Equation (106)). If we want our noninvariance group to include just some of the states then it will be a compact group, since unitary representations of compact groups can be finite dimensional. If we want to include all states then it will be a noncompact group since there are an infinite number of eigenstates of ( Z α ) 1 and all unitary representations of noncompact groups are infinite dimensional [9].
We can find a compact noninvariance group for the first N levels of the coupling constant, n = 1 , 2 , N . The dimensionality of our representation is
n = 1 N n 2 = N ( N + 1 ) ( 2 N + 1 ) 6 .
Mathematical analysis of the group SO(5) shows that this is the dimensionality of the irreducible symmetrical tensor representation of SO(5), given by the tensor with five upper indices T a b c , where a , b , . . = 1 , 2 , 3 , 4 or 5 [11]. Reducing this representation of SO(5) into its SO(4) components gives
( symm · tensor N ) S O ( 5 ) = ( 0 , 0 ) 1 2 , 1 2 N 1 2 , N 1 2 = ( symm · tensor n = 1 ) S O ( 4 ) ( symm · tensor n = 2 ) S O ( 4 ) ( symm · tensor n = N ) S O ( 4 )
which is precisely the structure of the first N levels of a hydrogenlike atom. If we want to include all levels then we guess that the appropriate noncompact group is SO(4,1), whose maximal compact subgroup is SO(4). Thus, we conjecture that all states form a representation of SO(4,1).
Consider the Lie algebra of O(4,1) and the general structure of its generators in terms of the canonical variables. The algebra of O(n) has n ( n 1 ) 2 generators so to extend the algebra of O(n) to O(n+1) takes n generators, which can be taken as the components of a n-vector. To extend the Lie algebra from O(4) to O(5) or O(4,1) we can choose the additional generators G a to be components of a four-vector G under O(4):
[ S a b , G c ] = i ( G b δ a c G a δ b c ) a , b , c = 1 , 2 , 3 , 4 .
If we apply Jacobi’s identity to S a b , G a , and G b and use Equation (184) we find
[ S a b , [ G a , G b ] ] = 0 .
We require that the Lie algebra closes, so [ G a , G b ] must be a linear combination of the generators, clearly proportional to S a b and we choose the normalization, such that
[ G a , G b ] = i S a b .
If we define
G 4 = S 40 = S ; G i = S i 0 = B i
and recall Equation (79)
L i = e i j k S j k A i = S i 4
then the additional commutation relations that realize SO(4,1) may be written in terms of L , A , B , and S:
L i , B j = i ϵ i j k B k L i , S = 0 S , A j = i B j S , B j = i A j A j , B k = i δ j k S B i , B j = i ϵ i j k L k .
The top two commutators show that B transforms as a three-vector under O(3) rotations and that S is a scalar under rotations. Alternatively, we can write the commutation relations in terms of the generators S A B , A , B = 0 , 1 , 2 , 3 , 4 :
S A B , S C D = i g A C S B D + g B D S A C g A D S B C g B C S A D
where g 00 = 1 , g a a = 1 .
The commutators above follow directly from the mathematical theory of SO(4,1), but the theory does not tell us what these generators represent, just their commutation properties. We now investigate the general features of the representations of SO(4,1) provided by the hydrogenlike atom and how to represent the generators in terms of the canonical variables.

6.2. Casimir Operators

The two Casimir operators of SO(4,1) are [120]
Q 2 = 1 2 S A B S A B = S 2 + B 2 A 2 L 2
and
Q 4 = w A w A = ( S L A x B ) 2 1 4 [ L · ( A + B ) ( A + B ) · L ] 2
where w A = 1 8 ϵ A B C D E S B C S D E .
For SO(4), we recall that for the SO(4) representations the structure of the generators in terms of the canonical variables led to the vanishing of one Casimir operator C 1 = L · A and, consequently, the appearance of only symmetrical tensor representations. We will find Q 4 vanishes for analogous reasons.
If B is a pseudovector, it is proportional to L , which is the only independent pseudovector that can be constructed from the dynamical variables. The coefficient of proportionality, a scalar, X need not commute with H:
B = X L [ X , L ] = 0 [ X , H ] 0
Because [ B i , B j ] = i e i j k L k it follows that X 2 = 1 and B would therefore be a constant multiple of L and not an independent generator. Thus B must be a vector and expressible as
B = f r + h p
where f and h are scalar functions of r , p 2 , and r · p . Accordingly we find
B · L = L · B = 0
Further, since B is a vector and A is a vector, A × B is a pseudovector and therefore is proportional to L :
A × B = Y L , [ Y , L ] = 0
For this equation to be consistent with the SO(4,1) commutation relations we find Y = S and, therefore
A × B = S L .
It follows from substituting L · A = 0 and Equations (194) and (196) into Equation (191) that for the SO(4,1) representations realized by the hydrogenlike atom
Q 4 = 0 .
As with the SO(4) symmetry, the dynamics of the hydrogen atom require that only certain representations of SO(4,1) appear. From the mathematical theory of irreducible infinite dimensional unitary representations of SO(4,1), we have the following results:
Class I : Q 4 = 0 ; Q 2 real , > 0 SU ( 2 ) × SU ( 2 ) content : ( Q ) I = ( 0 , 0 ) 1 2 , 1 2 ( 1 , 1 )
Class II : Q 4 = 0 , Q 2 = ( s 1 ) ( s + 2 ) , s = integer > 0 S U ( 2 ) × S U ( 2 ) content : ( Q ) I I = s 2 , s 2 s + 1 2 , s + 1 2
The class I representations are realized by the complete set of eigenstates of ( Z α ) 1 for the hydrogenlike atom. Note, however, that we have an infinite number of such class I representations since Q 2 may have any positive real value. We shall find that for Q 2 = 2 we may extend our group from SO(4,l) to SO(4,2). The class II representations are realized by the eigenstates of ( Z α ) 1 with principal quantum numbers from n = s + 1 to n becomes infinite. The first s levels could, if we desire, be described by SO(5).
In this section, we have analyzed the group structure and the representations using the complete set of eigenstates of ( Z α ) 1 for our basis. We might ask: what if we used energy eigenstates instead as a basis for the representations? From Section 4.3 we know that the quantum numbers and multiplicities of the ( Z α ) 1 eigenstates are precisely the same as those of the bound energy eigenstates. Thus, with the energy eigenstates as our basis, we would reach the same conclusions about the group structure as before but we would be including only the bound states in our representations and we would be ignoring all scattering states.

6.3. Relationship of the Dynamical Group SO(4,1) to the Conformal Group in Momentum Space

We can give a more complete analysis of the hydrogenlike atom in terms of SO(4,1) by considering the relationship between the four-dimensional rotations of the four-vector U a , with a = l , 2 , 3 , 4 , which we discussed in Section 4.6, and the group of conformal transformations in momentum space. Conformal transformations preserve the angles between directed curves, but not necessarily lengths. The rotations generated by the Runge-Lenz vector a and the angular momentum L leave the scalar product U a V a of four-vectors invariant and, therefore, are conformal transformations. The stereographic projection we employed is also a conformal transformation. Since the product of two conformal transformations is itself a conformal transformation, we must conclude that a generates a conformal transformation of the momentum three-vector p .
We must introduce two additional operators that correspond to the operators B and S introduced in Section 6.1 in order to express the most general conformal transformation. By employing the isomorphism between the generators L , a , B , and S of SO(4,1) and the generators of conformal transformations in momentum space we can immediately obtain expressions for the additional generators B and S in terms of the canonical variables, which is our objective. We need these additional generators to complete our SO(4,1) group for the hydrogen atom.
To derive the isomorphism we use the most convenient representation, namely that based on eigenstates of ( Z α ) 1 convenient for momentum space calculations ( ρ = ( p 2 + a 2 ) 2 a 2 ). Once established, the isomorphism becomes a group theoretical statement and it is independent of the particular representation.
The Conformal Group in Momentum Space
An arbitrary infinitesimal conformal transformation in momentum three-space may be written as
δ p j = δ a j + δ ω j k p k + δ ρ p j + p 2 δ c j 2 p j p · δ c
where δ ω j k = δ ω k j .
The terms in δ p j arise as follows:
δ a j translation generated by R · δ a δ ω j k rotation generated by J · δ ω , J i j = ϵ i j k J k δ ρ dilation generated by D δ ρ δ c j special conformal transformation generated by K · δ c
This is a ten parameter group with the generators ( R , J , D , K ) which obey the following commutation relations:
D , R j · = i R j D , J i = 0 D , K j = i K j R i , J k = i ϵ i j k R k K n , R m = 2 i ϵ n m r J r 2 i δ m n D J i , J k = i ϵ i k m J m R i , R j = 0 [ K i , J k ] = i ϵ i k m K m K j , K j = 0
There is an isomorphism between the algebra of the generators of conformal transformations and the dynamical noninvariance algebra of SO(4,1) of the hydrogen atom. Because J i is the generator of spatial rotations we make the association L i = J i . Comparing the differential change in p i from a transformation generated by a · δ ν (in the representation with ρ = ( p 2 + A 2 ) 2 a 2 )
δ p i = i a · δ ν , p i = 1 2 a p 2 a 2 δ ν i 2 p · δ ν p i
to the differential change in p i from a conformal transformation leads to the association
a i = 1 2 K i a a R i .
To confirm the identification we can use the commutation relations of the conformal group to show that the O(4) algebra of L and a corresponds precisely to that of J and 1 2 ( K a a R ) . This result alone suggests that our SO(4) degeneracy group should be considered as a subgroup of the larger group SO(4,1). It suggests introducing the operators
B = 1 2 K a + a R S = D .
The commutation relations of S and B which follow from Equations (205) and the commutation relations Equation (202) are identical to the commutation relations given for S and B in Section 6.1. Thus, by considering the a and L transformations in momentum space as conformal transformations, we were led to introduce the generators B and S and obtain the dynamical algebra SO(4,l). Further, we are led to the expressions for these generators in terms of the canonical variables.
By comparing the expression for a in terms of the conformal generators with our known expressions for a , Equation (120) or Equation (121), we obtain expressions for K i and R i in terms of the canonical variables. If we use the eigenstates convenient for configuration space calculations ( ρ = n / a r ) we make the identifications
K = 1 2 ( r p 2 + p 2 r ) r · p p p p · r r 1 4 r 2 R = r .
Substituting these results in the equation for B we find
B = 1 2 a p 2 r + r p 2 2 r · p p p p · r r 4 r 2 + a r 2
which is a manifestly Hermitean operator valid throughout Hilbert space. From Equations (204) and (205) we note that
B a = a r .
To compute D, we substitute the expressions for K and R into the commutation relation
D = i 2 [ K i , R i ]
obtaining the result
S = 1 2 ( p · r + r · p ) = D
which is identical to the generator of the scale change transformation D ( λ ) defined in Equation (109) in Section 4.3.
The significance of the generator D = S of the scale change in terms of SO(4,1) is apparent if we compute
e i λ D K a ± a R e i λ D = K a ± a R
where a = e λ a .
The unitary transformation e i λ D may be viewed as generating an inner automorphism of SO(4,1) which is an equivalent representation of SO(4) that is characterized by a different value of the quantity a or the energy. In other words, under the scale change e i λ D , the basis states for our representation of SO(4,1), | n l m ; a ) , transform to a new set, | n l m ; e λ a ) in agreement with our discussion in Section 4.3.
Because the algebra of our generators closes, we may also view e i λ D as transforming a given generator into a linear combination of the generators. With the definitions of a and B (Equation (204) and (205)), we can easily show that Equation (210), with the upper sign, may also be written
e i λ D B e i λ D = B cosh λ + a sinh λ .

7. The Group SO(4,2)

7.1. Motivation for Introducing SO(4,2)

We would like to express Schrodinger’s equation as an algebraic equation in the generators of some group [37,38]. As we are unable to do this with our SO(4,l) generators S A B we again expand the group. To guide us, we recall that to expand SO(3) to SO(4) we added a three-vector of generators A , and to expand SO(4) to SO(4,1), we added a four-vector of generators ( S , B ) . In both cases, this type of expansion produced a set of generators convenient for the study of the hydrogenlike atom. We guess that the appropriate expansion of SO(4,1) is obtained by adding a five-vector (under SO(4,1)) of generators Γ A to obtain SO(4,2) [37,38]. We can provide additional motivation for this choice by considering Schrodinger’s equation. The generators in terms of which we want to express this equation must be scalars under L i rotations. Additionally we know S = S 40 (Equation (187)) generates scale changes of Schrodinger’s equation. The fact that S 40 mixes the zero and four components of a five-vector suggests that Schrodinger’s equation may be expressed in terms of the components Γ 0 and Γ 4 which are scalars under L i , of the five vector Γ A . Since Γ A is a five-vector under SO(4,1), it must satisfy the equation
S A B , Γ C = i Γ B g A C Γ A g B C .
The spatial components of Γ A which are ( Γ 1 , Γ 2 , Γ 3 ) = Γ transform as a vector under rotations generated by L .
To construct the Lie algebra of SO(4,2) we require that the set of operators { Γ A , S A B ; A , B = 0 , 1 , 2 , 3 , 4 } must close under the operations of commutation. By applying Jacobi’s identity to Γ A , Γ B , and S A B , and requiring that Γ A and Γ B do not commute, we find
[ S A B , [ Γ A , Γ B ] ] = 0 A , B = 0 , 1 , 2 , 3 , 4 .
Because we require that our Lie algebra closes, the commutator of Γ A and Γ B must be proportional to S A B . We normalize Γ , so
[ Γ A , Γ B ] = i S A B A , B = 0 , 1 , 2 , 3 , 4 .
If we define
S A 5 = Γ A = S 5 A A = 0 , 1 , 2 , 3 , 4 .
and recall
A i = S i 4 B i = S i o L i = e i j k S j k S = S 40
then we may unite all the commutations relations of Γ A and S A B in the single equation:
[ S AB , S CD ] = i ( g AC S BD + g BD S AC g AD S BC g BC S AD )
where A , B , . . = 0 , 1 , 2 , 3 , 4 , 5 and g 00 = g 55 = 1 ; g a a = 1 , a = 1 , 2 , 3 , 4 .
These are the commutation relation for the Lie algebra of SO(4,2). In terms of A , B , L , S and Γ A the additional commutation relations for the noncommuting generators are [53]:
B i , Γ j = i Γ 0 δ i j Γ i , Γ j = i ϵ i j k L k A i , Γ j = i Γ 4 δ i j Γ i , Γ 0 = i B i L i , Γ j = i ϵ i j k Γ k Γ i , Γ 4 = i A i Γ 4 , Γ 0 = i S B i , Γ 0 = i Γ i A i , Γ 4 = i Γ i S , Γ 0 = i Γ 4 [ S , Γ 4 ] = i Γ 0

7.2. Casimir Operators

The Lie algebra of SO(4,2) is rank three, so it has three Casimir operators W 2 , W 3 , and W 4 [43]:
W 2 = 1 2 S AB S AB = Q 2 + Γ A Γ A
where Q 2 is the nonvanishing SO(4,1) Casimir operator Equation (190) and
W 3 = ϵ A BCDEF S AB S CD S EF
W 4 = S AB S BC S CD S DA .
Computation of W 3
We can show that W 3 = 0 from dynamical considerations similar to those used in the discussion of SO(4,1) Casimir operators. The only terms that can be included in W 3 are scalars that are formed from products of three generators with different indices
B · A × Γ , A · Γ × B , Γ · B × A Γ 4 L · B , Γ 0 L · A , S Γ · L
It is interesting that these terms are actually all pseudoscalars. Terms like B · A × L are simply not possible because of the structure of W 3 . We know that Γ = ( Γ 1 , Γ 2 , Γ 3 ) must not be pseudovector, otherwise it would be proportional to L . Since it is a vector, it must equal a linear combination of r and p . Therefore, we conclude
Γ · L = L · Γ = 0 .
Because Γ and B are both vectors and L is the only pseudovector we have Γ × B = Λ L . In order to determine the scalar Λ we evaluate the commutators
[ B k , ( Γ × B ) k ] , and [ Γ k , ( Γ × B ) k ]
and find
Γ × B = Γ 0 L = B × L .
The analogous equations for A and Γ , and for A and B , are
Γ × A = Γ 4 L = A × Γ A × B = S L = B × A .
From Equations (223) and (224), we see that because of the dynamical structure of the generators each of the quantities in the first line of Equation (220) is proportional to the quantity directly below in the second line. We also have shown that (Equations (194) and (221))
L · B = L · A = Γ · L = 0 .
Accordingly, each scalar in our list vanishes and
W 3 = 0 .
Computation of W 2
In order to compute W 2 we need to evaluate
Γ 2 Γ A Γ A = Γ 4 2 + Γ i Γ i Γ 0 2 .
From the structure of W 2 as shown in Equation (217), we see that Γ 2 must be a number since W 2 and Q 2 are both Casimir operators and therefore equal numbers for a particular representation. Accordingly, we have
[ Γ 2 , Γ A ] = 0 .
From this equation, we can deduce a lemma allowing us to easily evaluate W 2 and W 4 in terms of the number Γ 2 . Using Equation (227) and the definition of S A B Equation (213) we find
Γ A S A B + S A B Γ A = 0 .
Contracting Equation (212) with g A C gives
S A B Γ A Γ A S A B = 4 i Γ B .
Consequently, it must follow that
S A B Γ A = 2 i Γ B = Γ A S A B .
We are now able to evaluate the quantity
S A B S C B = i S A B [ Γ B , Γ C ] = i ( S A B Γ B Γ C S A B Γ C Γ B ) .
Using Equation (212) for the commutator of S A B with Γ C and Equation (229) for the contraction S A B Γ B we prove the lemma
S A B S C B = 2 i S C A Γ A Γ C + Γ 2 g A C .
The value of the SO(4,1) Casimir operator Q 2 = 1 2 g A C S A B S C B follows directly from the lemma:
Q 2 = 2 Γ 2 .
Accordingly, we have from Equation (217)
W 2 = 3 Γ 2 .
Computation of W 4
The Casimir operator W 4 can be written as
W 4 = S A B S B C S C D S D A + S A B S B 5 S 5 D S D A + S 5 B S B C S C D S D 5 + S 5 B S B 5 S 5 D S D 5 .
where B , D = 0 , 1 , 2 , 3 , 4 , 5 and A , C = 0 , 1 , 2 , 3 , 4 .
In order to evaluate W 4 in terms of Γ 2 we compute S A B S B C . Recalling Γ A = S A 5 we see
S A B S B C = Γ A Γ C + S A B S B C .
Substituting the lemma Equation (230), we find
S A B S B C = 2 i S A C Γ 2 g A C .
From Equation (229), it follows that
S 5 B S B C = 2 i Γ C .
Substituting Equations (231), (232), (235), and (236) into Equation (233) for W 4 we find
W 4 = 6 ( Γ 2 ) 2 24 Γ 2 .
The fact that the nonvanishing Casimir operators ( Q 2 , W 2 , and W 4 ) for SO(4,1) and SO(4,2) are given in terms of Γ 2 implies that the representation of SO(4,2) determines the particular representation of SO(4,l) appropriate to the hydrogenlike atom. In turn the value of Γ 2 is determined by the structure of the Γ s in terms of the canonical variables. In Section 7.4, we derive these structures and find that
Γ 2 = 1 .
Therefore, the quadratic SO(4,1) Casimir operator Q 2 has the value
Q 2 = 2
and the SO(4,2) Casimir operators have the values:
W 2 = 3 W 3 = 0 W 4 = 18 .
The researchers that have published different representations of SO(4,2) based on the hydrogen atom that give their Casimir operators all have W 2 = 3 (or its equivalent) and W 3 =0 [34,79,82], however, two authors have representations with W 4 = 0 [34,82] and one [79] has W 4 = 12 , as compared to our value of −18.
From the mathematical theory of representations, it follows that our representations of SO(4,1) and SO(4,2) are both unitary and irreducible. This means there is no subset of basis vectors that transform among themselves as either SO(4,1) or as SO(4,2).

7.3. Some Group Theoretical Results

In this section, we derive the transformation properties of the generators of SO(4,2) and then a novel contraction formula that will prove useful for situations in which we want to employ perturbation theory, for example, in our calculation of the radiative shift for the hydrogen atom in Section 8. We will work primarily with the SO(4,2) generators expressed as the combination of the SO(4,1) generators S A B and the five-vector Γ , with g A B = ( 1 , 1 , 1 , 1 , 1 ) where A , B = 0 , 1 , 2 , 3 , 4 .
Transformation Properties of the Generators
We can evaluate quantities like
A B Γ B ( θ ) e i S A B θ Γ B e i S A B θ no sum over A or B
by expanding the exponentials in an infinite series and then using the commutation relations Equations (212) and (213) of the generators S A B and Γ A , Γ B repeatedly. However, it is easier to solve the differential equations satisfied by A B Γ B and to use the appropriate boundary conditions. Differentiating Equation (238) and using the commutation relations, we obtain the equations
d d θ A B Γ B = g B B A B Γ A d 2 d θ 2 A B Γ B = g A A g B B A B Γ B
which have the solution
A B Γ B = Γ B cos g A A g B B θ + g B B g A A g B B Γ A sin g A A g B B θ .
Using a similar procedure we find
e i Γ A θ S A B e i Γ A θ = S A B cosh g A A θ + g A A Γ B sinh g A A θ
e i Γ A θ Γ B e i Γ A θ = Γ B cosh g A A θ + 1 g A A S A B sinh g A A θ
where no summation over A or B is implied.
These formulae, Equation (240)–(242), give the SO(4,2) transformation properties of the SO(4,2) generators.
The Contraction Formula
If we multiply Equation (241) from the right by e i Γ A θ and then contract from the left with Γ B , we obtain
B Γ B e i Γ A θ Γ B = ( 1 g A A Γ A 2 ) cosh g A A θ + 2 i Γ A g A A sinh g A A θ e i Γ A θ + g A A Γ A 2 e i Γ A θ
where we have used Γ 2 = 1 and Equation (229). Expanding the hyperbolic functions in terms of exponentials and collecting terms gives
B Γ B e i Γ A θ Γ B = 1 2 ( 1 + i Γ A g A A ) 2 e i ( Γ A i g A A ) θ + 1 2 ( 1 i Γ A g A A ) 2 e i ( Γ A + i g A A ) θ + g A A Γ A 2 e i Γ A θ .
A Fourier decomposition of a function Γ A may be written
f ( Γ A ) = 1 2 π d θ h ( θ ) e i Γ A θ .
Consequently we have
B Γ B f ( Γ A ) Γ B = 1 2 ( 1 + i Γ A g A A ) 2 f ( Γ A i g A A ) + 1 2 ( 1 i Γ A g A A ) 2 f ( Γ A + i g A A ) + g A A Γ A 2 f ( Γ A ) .
By performing a suitable rotation we can generalize this formula from functions of Γ A to functions of Γ A n A where n A n A = ± 1 . For n 2 = 1 we start with Γ A = Γ 0 and rotate to obtain a very general result
B Γ B f ( n Γ ) Γ B = 1 2 ( n Γ + 1 ) 2 f ( n Γ + 1 ) + 1 2 ( n Γ 1 ) 2 f ( n Γ 1 ) ( n Γ ) 2 f ( n Γ ) .
We will have occasion to apply this formula for the special case
f ( n Γ ) = 1 Γ n ν .
Using the representation
1 Γ n ν = 0 d s e ν s e Γ n s
we obtain the result
Γ A 1 Γ n ν Γ A = 2 ν 0 d s e ν s d d s ( sinh 2 s 2 e Γ n s )
which is in a form convenient for perturbation calculations.
Derivation of the Γ A in terms of the Canonical Variables
For our basis states, we shall use eigenstates of ( Z α ) 1 convenient for configuration space calculations ( ρ = n a / r ). We choose these states rather than those convenient for momentum space calculations because they lead to simpler expressions for the Γ A in terms of the canonical variables, although the expression for a is slightly more complicated. Thus, our states obey the equation
1 K 1 ( a ) n | n l m ) = 0
where K 1 ( a ) is given by Equation (116). We know that K 1 1 must commute with the generators of the SO(4) symmetry group ( a ) i = S i 4 and S i j = ϵ i j k L k . This suggests that we choose
Γ 0 = [ K 1 ( a ) ] 1 = a r p 2 + a 2 2 a 2 a r = 1 2 r p 2 r a + a r
so that
Γ 0 n | n l m ) = 0 .
This last equation is the Schrodinger equation expressed in our language of SO(4,2): our states | n l m ) are eigenstates of Γ 0 with eigenvalue n.
To find Γ 4 , we calculate Γ 4 = i [ S , Γ 0 ] , using Equation (209) for S,
Γ 4 = a r p 2 a 2 2 a 2 a r = 1 2 r p 2 r a a r .
Sometimes it is convenient to use the linear combinations
Γ 0 Γ 4 = a r Γ 0 + Γ 4 = r p 2 r a
which can be used to express the dipole transition operator [33]. We can find Γ i from Equation (216), Γ i = i [ B i , Γ 0 ]
Γ i = r p i r
which we might have guessed initially since [ r p i , r p j ] L k . Every component of Γ A is Hermitean, consequently the generators S A B given by the commutators Equation (213) are also Hermitian. We may explicitly verify that these expressions for Γ A lead to a consistent representation of all generators in the SO (4,2) Lie algebra.
Under a scale change generated by S, Γ i is invariant and Γ 4 and Γ 0 transform in the same manner as a and B (Equation (210)): they retain their form but a is transformed into e λ a :
e i λ S Γ 0 Γ 4 e i λ S = 1 2 r p 2 r e λ a ± e λ a r
The scale change generates an inner automorphism of SO(4,2) characterized by a different value of the parameter a.

7.4. Subgroups of SO(4,2)

The two most significant subgroups are [53]:
1. L i , a i or S j k , S i 4 , forming an SO(4) subgroup. These generators commute with Γ 0 and therefore constitute the degeneracy group for states of energy a 2 / ( 2 m ) and fixed principal quantum number n (or fixed coupling constant n a / m ). The Casimir operator for this subgroup is
a 2 + L 2 = n 2 1 = Γ 0 2 1 .
We discussed this subgroup in Section 4.2 in terms of L and A and the states | n l m . The same results are obtained with the generators L and a with the states | n l m ) . For example, we have the raising and lowering operators for m and l (Equations (93) and (94)). With the definition
L ± = L 1 ± i L 2
it follows that
[ L 3 , L ± ] = ± L ±
which gives
L ± | n l m ) = ( l ( l + 1 ) m ( m ± 1 ) | n l m ± 1 ) L 3 | n l m ) = m | n l m ) .
for l 1 . In analogy to L ± one can define
a ± = a 1 ± i a 2
which obey the relations
[ a 3 , a ± ] = ± L 3 [ L 3 , a ± ] = ± a ±
and
a ± | n l m ) = ( n 2 ( l + 1 ) 2 ) ( l + 2 ± m ) ( l + 1 ± m ) 4 ( l + 1 ) 2 1 1 2 | n l + 1 m ± 1 ) ± ( n 2 l 2 ) ( l m ) ( l 1 m ) 4 l 2 1 1 2 | n l 1 m ± 1 )
for l 1 . The action of a ± is not directly analogous to that of L ± , because we are using | n l m ) as basis states. If we used | n a 3 l 3 = m ) as basis states, the action would be similar. An operator that only changes the angular momentum is a 3
a 3 | n l m = ( n 2 ( l + 1 ) 2 ) ( ( l + 1 ) 2 m 2 ) 4 ( l + 1 ) 2 1 1 2 | n l + 1 m + ( n 2 l 2 ) ( l 2 m 2 ) 4 l 2 1 1 2 | n l 1 m .
for l 1 . Since a 3 commutes with L 3 and Γ 0 , it does not change n or m.
2. Γ 4 , S = S 40 , Γ 0 , forming a SO(2,1) subgroup. These operators commute with L but not with Γ 0 , hence then can change n but not L or m. The Casimir operator for this subgroup is
Γ 0 2 Γ 4 2 S 2 = L 2 = l ( l + 1 ) .
We can define the operators [53]
j 1 = Γ 4 j 2 = S j 3 = Γ 0
with commutators
[ j 1 , j 2 ] = i j 3 [ j 2 , j 3 ] = i j 1 [ j 3 , j 1 ] = i j 2
We can define the raising and lowering operators
j ± = j 1 ± i j 2 = Γ 4 ± i S
which obey the commutation relations
[ j ± , j 3 ] = j ± .
We find (in analogy to Equation (261))
Γ 0 | n l m ) = n | n l m ) ( Γ 4 ± i S ) | n l m ) = n ( n ± 1 ) l ( l + 1 ) | n ± 1 l m )
We can express the action of Γ 0 Γ 4 = a r on our states
a r | n l m ) = 1 2 ( ( n ) ( n l ) l ( l + 1 ) ) 1 2 | n 1 l m ) + n | n l m ) + 1 2 ( ( n ) ( n + l ) l ( l + 1 ) ) 1 2 | n + 1 l m )
As previously mentioned, the operator S generates scale changes as shown in Equation (257), where the value of a is changed. We can also express the action of S equivalently as transforming Γ 0 into Γ 4
e i S λ Γ 0 e i S λ = Γ 0 cosh λ Γ 4 sinh λ e i S λ Γ 4 e i S λ = Γ 4 cosh λ Γ 0 sinh λ .

7.5. Time Dependence of SO(4,2) Generators

For a generator to be a constant it must commute with the Hamiltonian as discussed in Section 2.1. Because the SO(4,2) group is the non-invariance or spectrum generating group, the additional generators do not all commute with the Hamiltonian and may have a harmonic time dependence as discussed in Section 2.2. It is notable that as far as we know only one paper considers the time dependence of the generators of non-invariance groups in general and one considers SO(4,2) specifically [97,136]. Our results certainly clarify and make explicit the time dependence, and show that it is just a particular aspect of the SO(4,2) transformations. In our representation with basis states | n l m ; a ) , the Hamiltonian, which is the generator of translations in time, has been transformed into Γ 0 and the Schrodinger energy eigenvalue equation has become Γ 0 | n l m ) = n | n l m ) . Accordingly, all of the generators that commute with Γ 0 are constants of the motion, which includes a , L . The other operators, B , Γ , S , Γ 4 have a time dependence given by Equations (241) and (242), for example
S ( t ) = e i H t S ( 0 ) E i H t = e i Γ 0 t S e i Γ 0 t = S cos t + Γ 4 sin t .
Γ 4 ( t ) = e i H t Γ 4 ( 0 ) E i H t = e i Γ 0 t Γ 4 e i Γ 0 t = Γ 4 cos t S sin t .
Consequently, terms like j ± have a simple exponential time dependence
j ± ( t ) = j ± ( 0 ) e ± i t .
Similarly Γ ± i B has an exponential time dependence.

7.6. Expressing the Schrodinger Equation in Terms of the Generators of SO(4,2)

We can write the Schrodinger equation for the energy eigenstate E n = a n 2 / 2 m of a particle in a Coulomb potential in terms of the SO(4,2) generators, which are expressed in terms of the energy a 2 / 2 m , by making a scale change. From Section 4.3, Equation (114), the relationship between the Schrodinger energy eigenstate | n l m and the eigenstate of ( Z α ) 1 is:
| n l m ; a ) = e i S λ n ρ ( a n ) | n l m
where
e λ n = a n a ρ ( a n ) = n a n r .
Substituting Equation (277) into the eigenvalue equation Equation (253) for | n l m ; a ) and employing the transformation Equation (273), we find the usual Schrodinger equation can be expressed in SO(4,2) terms as
( Γ n n ) ρ ( a n ) | n l m = 0
where
Γ n Γ A n A = Γ 0 n 0 + Γ i n i + Γ 4 n 4
n o = cosh λ n = a 2 + a n 2 2 a a n , n i = 0 , n 4 = sinh λ n = a 2 a n 2 2 a a n
and n A n A = n 4 2 n 0 2 = 1 .
Equation (279) expresses Schrodinger’s equation for an ordinary energy eigenstate | n l m with energy E N = a n 2 / 2 m in the language of SO(4,2). It shows the relationship between these energy eigenstates and the basis states of ( Z α ) 1 used for the SO(4,2) representation,

8. SO(4,2) Calculation of the Radiative Shift for the Schrodinger Hydrogen Atom

In the 1930’s, it was generally believed that the Dirac equation predicted the energy levels of the hydrogen atom with excellent accuracy, but there were some questions about the prediction that the energy levels for a given principal quantum number and given total angular momentum were independent of the orbital angular momentum. To finally resolve this issue, in 1947, Willis Lamb and his student Robert Retherford at Columbia University in New York City employed rf spectroscopy and exploited the metastability of the hydrogen 2 s 1 / 2 level in a beautiful experiment and determined that the 2 s 1 / 2 and 2 p 1 / 2 levels were not degenerate and that the energy difference between the levels was about 1050 MHz, or 1 part in 10 6 of the 2 s 1 / 2 level [5,148]. Shortly thereafter Hans Bethe [6] published a ground breaking nonrelativistic quantum theoretical calculation of the shift assuming it was due to the interaction of the electron with the ground state electromagnetic field of the quantum vacuum field. This radiative shift accounted for about 96% of the measured shift. The insight that one needed to include the interaction of the atom with the vacuum fluctuations and how one could actually do it ushered in the modern world of quantum electrodynamics [7]. Here, we compute in the non-relativistic dipole approximation and to first order in the radiation field, as did Bethe, the radiative shift, but we use group theoretical methods based on the SO(4,2) symmetry of the non-relativistic hydrogen atom as developed in this paper. Bethe’s calculation required the numerical sum over intermediate states to obtain the average value of the energy of the states contributing to the shift. In our calculation, we do not use intermediate states, and we derive an integral equivalent to Bethe’s log, and more generally derive the shift for all levels in terms of a double integral.
An expression for the radiative shift Δ N L for energy level E N of a hydrogen atom in a state | N L can be readily obtained using second order perturbation theory (to first order in α the radiation field) [6,149,150,151]
Δ N L = 2 α c 3 π m 2 n s 0 ω c d ω ( E n E N ) N L | p i | n n | p i | N L E n E N + ω i ϵ ,
where ω C is a cutoff frequency for the integration that we will take as ω c = m .
This expression, which is the same as Bethe’s, has been derived by inserting a complete set of states | n n | , a step that we eliminate with our group theoretical approach:
Δ N L = 2 α 3 π m 2 0 ω c d ω N L | p i H E H ( E ω ) i ϵ p i | N L
If we add and subtract ω from the numerator, we find the real part of the shift is
R e Δ N L = 2 α 3 π m 2 R e 0 ω c d ω [ N L | p 2 | N L ω Ω N L ]
where
Ω N L = N L | p i 1 H E N + ω i ϵ p i | N L
and
H = p 2 2 m Z α r .
The imaginary part of the shift gives the width of the level [7].
The matrix element Ω N L can be converted to a matrix element of a function of the generators Γ A taken between eigenstates | n l m ) of ( Z α ) 1 . To do this we insert factors of 1 = r 1 r and use the definitions of the Γ A in terms of the canonical variables, Equations (254)–(256). Letting the parameter a take the value a N , we obtain the result
Ω N L = m ν N 2 ( N L | Γ i Γ n ( ξ ) ν Γ i | N L )
where
n 0 ( ξ ) = 2 + ξ 2 1 + ξ = cosh ϕ n i = 0 n 4 ( ξ ) = ξ 2 1 + ξ = sinh ϕ
and
ξ = ω | E N | ν = N 1 + ξ = N e ϕ .
From the definitions we see ϕ = 1 2 l n ( 1 + ξ ) > 0 and n A ( ξ ) n A ( ξ ) = 1 . The quantity
ν = m Z α 2 m ( E N ω )
may be considered the effective principal quantum number for a state of energy E N ω . The contraction over i in Ω N L may be evaluated using the group theoretical formula Equation (250):
Ω N L = 2 m ν 2 N 2 0 d s e ν s d d s sinh 2 s 2 M N L ( s ) m ν N 2 ( N L | Γ 4 1 Γ n ( ξ ) ν Γ 4 | N L ) + m ν N 2 ( N L | Γ 0 1 Γ n ( ξ ) ν Γ 0 | N L )
where
M N L ( s ) = ( N L | e Γ n ( ξ ) s | N L ) .
In order to evaluate the last-two terms in Ω N L , we can express the action of Γ 4 on our states in terms of Γ n ( ξ ) ν . Substituting the equation
Γ 0 | N L ) = N | N L )
into the expression for Γ n ( ξ ) ν , with n ( ξ ) given by Equation (288), gives
Γ 4 = N 1 sinh ϕ Γ n ( ξ ) ν
when acting on the state |NL). If we substitute Equation (293) into the expression for the R e Δ E N L Equation (284) and simplify using the virial theorem
( N L | p 2 | N L ) = a N 2
we find that the term in p 2 exactly cancels the last two terms in Ω N L , yielding the result
R e Δ E N L = 4 m α ( Z α ) 4 3 π N 4 0 ϕ c d ϕ sinh ϕ e ϕ 0 d s e ν s d d s sinh 2 s 2 M N L ( s )
where
ϕ c = 1 2 l n 1 + ω c | E N | = 1 2 l n 1 + 2 N 2 ( Z α ) 2
and ω c = m .
Comparison to the Bethe Logarithm
The first order non-relativistic radiative shift is commonly given in terms of the Bethe logarithm γ ( N , L ) , which is interpreted as the average over all states, including scattering states, of l n | E n E N | 1 2 m ( Z α ) 2 [149]:
γ ( N , L ) n S E n E N N 0 p i n n p i N 0 = n S E n E N N L p i n n p i N L ln E n E N 1 2 m ( Z α ) 2 .
We use the dipole sum rule [150]
2 n s E n E N N p i n n p i N = N 2 V N
and apply it for the Coulomb potential 2 V ( r ) = 4 π Z α δ ( r ) . The use of the Bethe log allowed Bethe to take the logarithmic expression obtained from the frequency integration outside the summation over the states, and replace it with the average value. Only the S states contribute to the expectation value in Equation (297), giving, from Equation (282), an expression for the shift
Re Δ E N L = 4 m 3 π α ( Z α ) 4 1 N 3 δ L 0 ln 2 ( Z α ) 2 γ ( N , L ) .
Comparing the shift in terms of M N L Equation (294) to the shift in terms of γ ( N , L ) we find that the Bethe log is
γ ( N , L ) = 0 ϕ c d ϕ sinh ϕ e ϕ 0 d s e ν s d d s sinh 2 s 2 M N L ( s ) δ L 0 l n 1 ( Z α ) 2

8.1. Generating Function for the Shifts

We can derive a generating function for the shifts for any eigenstate characterized by N and L if we multiply Equation (291) by N 4 e β N and sum over all N , N L + 1 . To simplify the right side of the resulting equation, we use the fact that the O(2,1) algebra of Γ 0 , Γ 4 , and S closes. We can compute the sum on the right hand side:
N = L + 1 e β N M N L = N = L + 1 ( N L | e j · ψ | N L ) .
where
e j · ψ e β Γ 0 e s Γ n ( ξ ) .
We perform a j transformation, such that
e j · ψ e j 3 ψ = e Γ 0 ψ
so
N = L + 1 e β N M N L = N = L + 1 ( N L | e j 3 ψ | N L ) = N = L + 1 e N ψ
= e ψ ( L + 1 ) 1 e ψ .
In order to find a particular M N L , we must expand the right hand side of the equation in powers of e β and equate the coefficients to those on the left hand side. First, we need an equation for e ψ . This can be obtained using the isomorphism between j and the Pauli σ matrices:
( Γ 4 , S , Γ 0 ) ( j 1 , j 2 , j 3 ) ( i 2 σ 1 , i 2 σ 2 , 1 2 σ 3 )
Using the formula
e i 2 s n · σ = cos s 2 + i n · σ sin s 2
where | n | = 1 , we find
cosh ψ 2 = cosh β 2 cosh s 2 + sinh β 2 sinh s 2 cosh ϕ .
We can rewrite this equation in a form easier for expansion
e + 1 2 ψ = d e 1 2 β + b e 1 2 β e 1 2 ψ
where
d = cosh s 2 + sinh s 2 cosh ϕ b = cosh s 2 sinh s 2 cosh ϕ .
Let β become very large and iterate the equation for e 1 2 ψ to obtain the result
e ψ = A e β 1 + A 1 e β + A 2 e 2 β +
where
A = A 0 = 1 d 2 A 1 = 2 d b d 1 A 2 = 3 d 2 b d 1 2 2 2 b d 1
Note b d 1 = d 1 sinh 2 s 2 sinh 2 ϕ .

8.2. The Shift between Degenerate Levels

Expressions for the energy shift between degenerate levels with the same value of N may be obtained directly from the generating function using Equations (294) and (304). We find
N = L + 1 e β N N 4 Re Δ E N L N = L + 1 e β N N 4 Re Δ E N L =
4 m α ( Z α ) 4 3 π 0 ϕ c d ϕ e ϕ sinh ϕ 0 d s e ν s d d s sinh 2 s 2 e ψ ( L + 1 ) e ψ L + 1 1 e ψ .
For an example, consider L = 1 , L = 0 . For the shifts between levels we obtain
N = 2 e β N N 4 Re Δ E N O Δ E N 1 + Re Δ E 10 e β =
4 m α ( z α ) 4 3 π 0 ϕ c d ϕ e ϕ sinh ϕ 0 d s e ν s d d s sinh 2 s 2 e ψ
Substituting Equation (310) for e ψ , using the coefficient A A N 1 of e N β , gives
R e ( Δ E N 0 Δ E N 1 ) = 4 m α ( Z α ) 4 3 π N 4 0 ϕ c d ϕ e ϕ sinh ϕ 0 d s e ν s d d s sinh 2 s 2 A A N 1 .
where A and A N 1 are given in Equation (311) in terms of the variables of integration s and ϕ .
General Expression for M N L
Once we have a general expression for M N L , we can use Equation (294) to calculate the shift for any level E N L . We can obtain expressions for the values of M N L by letting β become large, expanding the denominator in Equation (304) and equating coefficients of powers of e β . For large β , we have large ψ . We have
e ψ ( L + 1 ) 1 e ψ = m = 1 e ψ ( m + L )
and for large β it follows from Equation (310) that
N = L + 1 e β N M N L = m = 1 e β A 1 + A 1 e β + m + L .
Using the multinomial theorem [124], the right side of the equation becomes
m = 1 A m + L r , s , t , ( m + L ) ! r ! s ! t ! A 1 s A 2 t e β ( m + L + s + 2 t + ) .
where r + s + t + = m + L .
To obtain the expression for M N L , we note N is the coefficient of β so N = m + L + s + 2 t + = r + 2 s + 3 t + Accordingly we find
M N L = r , s , t , A ( r + s + t + ) ( r + s + t + ) ! r ! s ! t ! A 1 s A 2 t
where r + 2 s + 3 t + = N and r + s + t + >L.
By applying this formula, we obtain the results:
N = 1:
M 10 = A
N = 2:
M 20 = A 2 + A A 1 M 21 = A 2
Shifts for N = 1 and N = 2
To illustrate these results, we can calculate the shift for a given energy level using Equation (294). For N = 1 , we note from Equation (318) that M 10 = A , and from Equation (311) that A = 1 / d 2 . We find that the real part of the radiative shift for the 1 S ground state is
R e Δ E 10 = 4 m α ( Z α ) 4 3 π 0 ϕ c d ϕ e ϕ sinh ϕ 0 d s e s e ϕ d d s 1 coth s 2 + cosh ϕ 2
For the shift between two states Equation (314) can be used. For the N = 2 Lamb shift between 2S-2P states, the radiative shift to first order in α is
R e ( Δ E 20 Δ E 21 ) = m α ( Z α ) 4 6 π 0 ϕ c d ϕ e ϕ sinh 3 ϕ 0 d s e 2 s e ϕ d d s 1 coth s 2 + cosh ϕ 4
The s integral can be computed in terms of a Jacobi function of the second kind [127].
As one check on our group theoretical methods, we can compare our matrix elements ( 10 | e i S ϕ | n 0 ) with those of Huff [45]. To go from Equation (301) to Equation (302), we did a rotation R ( ϕ ) = e i ϕ S generated by S that transformed Γ n into Γ 0 . For N , L = 1 , 0 we have
M 10 = ( 10 | e Γ n s | 10 ) = ( 10 | R ( ϕ ) e Γ 0 s R 1 ( ϕ ) | 10 ) = 1 ( cosh s 2 + sinh s 2 cosh ϕ ) 2
Expanding the hyperbolic functions, we get
M 10 = 4 e s ( 1 + cosh ϕ ) 2 1 e s tanh 2 s 2 2
= 4 ( 1 + cosh ϕ ) 2 n = 1 n e n s tanh 2 ϕ 2 n 1 .
We can also compute M 10 by inserting a complete set of states and using Γ 0 | n 0 ) = n | n 0 ) in Equation (322). Because the generator S is a scalar, only states with L = 0 , m = 0 can contribute:
M 10 = n l m e n s | n ( 10 | R ( ϕ ) | n 0 ) | 2 .
Comparing this to Equation (323), we make the identification
| ( 10 | R ( ϕ ) | n 0 ) | 2 = 4 n ( 1 + cosh ϕ ) 2 tanh 2 ϕ 2 n 1 .
Huff computes this matrix element by analytically continuing the known O ( 3 ) matrix element of e i J y ϕ obtaining
| 10 | R ( ϕ ) | n 0 | 2 = 4 n cosh 2 ϕ 1 tanh 2 ϕ 2 n · [ 2 F 1 ( 0 , 1 ; n ; 1 2 ( 1 cosh ϕ ) ) ] 2 .
By algebraic manipulation and using 2 F 1 = 1 for the arguments here, we see that this result agrees with our much more simply expressed result from group theory.

9. Conclusions and Future Research

Measuring and explaining the properties of the hydrogen atom has been central to the development of modern physics over the last century. One of the most useful and profound ways to understand its properties is through its symmetries, which we have explored in this paper, beginning with the symmetry of the Hamiltonian, which reflects the symmetry of the degenerate levels, then the larger non-invariance and spectrum-generating groups, which include all of the states. The successes in using symmetry to explore the hydrogen atom led to use of symmetry to understand and model other physical systems, particularly elementary particles.
The hydrogen atom will doubtless continue to be one of testing grounds for fundamental physics. Researchers are exploring the relationship between the hydrogen atom and quantum information [152], the effect of non-commuting canonical variables [ x i , x j ] 0 on energy levels [153,154,155], muonic hydrogen spectra [156], and new physics using Rydberg states [157,158,159,160,161,162]. The ultra high precision of the measurement of the energy levels has led to new understanding of low Z two body systems, including muonium, positronium, and tritium [151]. As mentioned in the introduction, measurements of levels shifts are currently being used to measure the radius of the proton [2]. We can expect that further investigations of the hydrogen atom and hydrogenlike atoms will continue to reveal new vistas of physics and that symmetry considerations will play an important part.

Funding

This research received no external funding.

Acknowledgments

I would especially like to thank Lowell S. Brown very much for the time and energy he has spent on my education as a physicist. I also thank Peter Milonni for his encouragement and his helpful comments and many insightful discussions and MDPI for the invitation to be the Guest Editor for this special issue of Symmetry on Symmetries in Quantum Mechanics.

Conflicts of Interest

The author declares no conflict of interest.

References and Notes

  1. Brown, L. Bound on Screening Corrections in Beta Decay. Phys. Rev. 1964, 135, B314. [Google Scholar] [CrossRef]
  2. Beyer, A. The Rydberg constant and proton size from atomic hydrogen. Science 2017, 358, 79–85. [Google Scholar] [CrossRef] [Green Version]
  3. Mohr, P.D.B.N.; Taylor, B.N. CODATA recommended values of the fundamental physical constants: 2014. Rev. Mod. Phys. 2016, 88, 035009. [Google Scholar] [CrossRef] [Green Version]
  4. Rigden, J. Hydrogen, The Essential Element; Harvard University Press: Cambridge, MA, USA, 2002. [Google Scholar]
  5. Lamb, W.; Retherford, R. Fine Structure of the Hydrogen Atom by a Microwave Method. Phys. Rev. 1947, 72, 241. [Google Scholar] [CrossRef]
  6. Bethe, H. The Electromagnetic Shift of Energy Levels. Phys. Rev. 1947, 72, 339. [Google Scholar] [CrossRef]
  7. Maclay, J. History and Some Aspects of the Lamb Shift. Physics 2020, 2, 8. [Google Scholar] [CrossRef] [Green Version]
  8. Noether, E. Invariante Variationsprobleme. In Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen; Akademie der Wissenschaften zu Göttingen: Göttingen, Germany, 1918; pp. 235–257. [Google Scholar]
  9. Hamermesh, M. Group Theory; Adddison-Wesley Publishing Co.: Reading, MA, USA, 1962. [Google Scholar]
  10. Weyl, H. The Theory of Groups and Quantum Mechanics, 2nd ed.; Dover Reprint; Dover Publications: New York, NY, USA, 1928. [Google Scholar]
  11. Wigner, E. Group Theory and Its Application to the Quantum Mechanics of Atomic Spectra; Academic Press: New York, NY, USA, 1959. [Google Scholar]
  12. Bargmann, V. Zur Theorie des Wasserstffatoms. Z. Phys. 1936, 99, 576. [Google Scholar] [CrossRef]
  13. Laplace, P. A Treatise of Celestial Mechanics; Forgotten Books: Dublin, Ireland, 1827. [Google Scholar]
  14. Pauli, W. Uber das Wasserstoffspektrum vom Standpunkt der neuen Quantummechanik. Z. Phys. 1926, 36, 336–363. [Google Scholar] [CrossRef]
  15. McIntosh, H. On Accidental Degeneracy in Classical and Quantum Mechanics. Am. J. Phys. 1959, 27, 620–625. [Google Scholar] [CrossRef]
  16. Hulthen, E. Über die quantenmechanische Herleitung der Balmerterme. Z. Phys. 1933, 86, 21–23. [Google Scholar] [CrossRef]
  17. We employ natural Gaussian units so = 1, c = 1, and α = (e2/ℏc) ≈ 1/137. The notation for indices and vectors is μ,ν,.. = 0,1,2,3; i,j,. = 1,2,3; pμpμ = − p 0 2 + p 2 , p = (p1, p2, p3), gμν = (−1,1,1,1)
  18. Fock, V. Zur Theorie des Wasserstoffatoms. Z. Phys. 1935, 98, 145–154. [Google Scholar] [CrossRef]
  19. Dirac, P. Quantum Mechanics, 1st ed.; Oxford University Press: Oxford, UK, 1930. [Google Scholar]
  20. Gell-Mann, M. Symmetries of Baryons and Mesons. Phys. Rev. 1962, 125, 1067. [Google Scholar] [CrossRef] [Green Version]
  21. Schwinger, J. Coulomb’s Green’s Function. J. Math. Phys. 1964, 5, 1606–1608. [Google Scholar] [CrossRef]
  22. Ne’eman, Y. Algebraic Theory of Particle Physics; Benjamin: New York, NY, USA, 1967. [Google Scholar]
  23. Ne’eman, Y. Derivation of strong interactions from a gauge invariance. Nucl. Phys. 1961, 26, 222–229. [Google Scholar] [CrossRef]
  24. Gell-Mann, M. A schematic model of baryons and mesons. Phys. Lett. 1964, 8, 214–215. [Google Scholar] [CrossRef]
  25. Gell-Mann, M.; Ne’eman, Y. The Eightfold Way; Benjamin: New York, NY, USA, 1964. [Google Scholar]
  26. Dothan, Y.; Gell-Man, M.; Ne’eman, Y. Series of Hadron Energy Levels as Representations of Non-Compact Groups. Phys. Letters 1965, 17, 148. [Google Scholar] [CrossRef]
  27. Nambu, Y. Infinite-Component Wave Equations with Hydrogenlike Mass Spectra. Phys. Rev. 1967, 160, 1171. [Google Scholar] [CrossRef]
  28. Dyson, F. Symmetry Groups in Nuclear and Particle Physics; Benjamin: New York, NY, USA, 1966. [Google Scholar]
  29. Thomas, L. On the Unitary Representations of the Group of de Sitter Space. Ann. Math. 1941, 42, 113–126. [Google Scholar] [CrossRef]
  30. Harish-Chandra, Representations of Semisimple Lie Groups II. Trans. Am. Math. Soc. 1954, 76, 26. [CrossRef]
  31. Barut, A.; Budini, P.; Fronsdal, C. Two examples of covariant theories with internal symmetries involving spin. Proc. Roy. Soc. 1966, A291, 106–112. [Google Scholar]
  32. Malkin, I.; Man’ko, V. Symmetry of the Hydrogen Atom. Sov. Phys. Jetp Lett. 1966, 2, 146. [Google Scholar]
  33. Barut, A.; Kleinert, K. Transition Probabilities of the Hydrogen Atom from Noncompact Dynamical Groups. Phy. Rev. 1967, 156, 1541. [Google Scholar] [CrossRef]
  34. Barut, A.; Kleinert, H. Transition Form Factors in the H Atom. Phys. Rev. 1967, 160, 1149. [Google Scholar] [CrossRef]
  35. Bander, M.; Itzykson, C. Group Theory and the Hydrogen Atom (I). Rev. Mod. Phys. 1966, 38, 330. [Google Scholar] [CrossRef] [Green Version]
  36. Bander, M.; Itzykson, C. Group Theory and the Hydrogen Atom (II). Rev. Mod. Phys. 1966, 38, 346. [Google Scholar] [CrossRef] [Green Version]
  37. Fronsdal, C. Infinite Multiplets and Local Fields. Phys. Rev. 1967, 156, 1653. [Google Scholar] [CrossRef]
  38. Fronsdal, C. Infinite Multiplets and the Hydrogen Atom. Phys. Rev. 1967, 156, 1665. [Google Scholar] [CrossRef]
  39. Barut, A.; Fronsdal, C. On Non-Compact Groups. II Representations of the 2+1 Lorentz Group. Proc. R. Soc. 1965, A287, 532–548. [Google Scholar]
  40. Fronsdal, C. Relativistic Lagrangian Field Theory for Composite Systems. Phys. Rev. 1968, 171, 1811. [Google Scholar] [CrossRef]
  41. Pratt, R.; Jordan, T. Coulomb Group Theory for and Spin. Phys. Rev. 1969, 188, 2534. [Google Scholar]
  42. Fronsdal, C. Relativistic and Realistic Classical Mechanics of Two Interacting Point Particles. Phys. Rev. D 1971, 4, 1689. [Google Scholar] [CrossRef]
  43. Kyriakopoulos, R. Dynamical Groups and the Bethe-Salpeter Equation. Phys. Rev. 1968, 174, 1846. [Google Scholar] [CrossRef]
  44. Lieber, M. O(4) Symmetry of the Hydrogen Atom and the Lamb Shift. Phys. Rev. 1968, 174, 2037. [Google Scholar] [CrossRef]
  45. Huff, R. Simplified Calculation of Lamb Shift Using Algebraic Techniques. Phys. Rev. 1969, 186, 1367. [Google Scholar] [CrossRef]
  46. Musto, R. Generators of SO(4,1) for the Quantum Mechanical Hydrogen Atom. Phys. Rev. 1966, 148, 1274. [Google Scholar] [CrossRef]
  47. Barut, A.; Bornzin, G. SO(4,2)-Formulation of the Symmetry Breaking in Relativistic Kepler Problems with of without Magnetic Charge. J. Math. Phys. 1971, 12, 841–846. [Google Scholar] [CrossRef] [Green Version]
  48. Barut, A.; Kleinert, H. Current Operators and Majorana Equation for the Hydrogen Atom from Dynamical Groups. Phys. Rev. 1967, 157, 1180. [Google Scholar] [CrossRef]
  49. Mack, G.; Todorov, I. Irreducibility of the Ladder representations when restricted to the Poincare Subgroup. J. Math Phys. 1969, 10, 2078–2085. [Google Scholar] [CrossRef]
  50. Decoster, A. Realization of the Symmetry Groups of the Nonrelativistic Hydrogen Atom. Nuovo Cimento 1970, 68A, 105–117. [Google Scholar] [CrossRef]
  51. Englefield, M. Group Theory and the Coulomb Problem; Wiley-Interscience: New York, NY, USA, 1972. [Google Scholar]
  52. Barut, A. Dynamical Groups; University of Canterbury Press: Christchurch, New Zealand, 1972. [Google Scholar]
  53. Bednar, M. Algebraic Treatment of Quantum-Mechanical Models with Modified Coulomb Potentials. Ann. Phys. 1973, 75, 305–331. [Google Scholar] [CrossRef]
  54. Wulfman, C.; Takahata, Y. Noninvariance Groups in Molecular Quantum Mechanics. J. Chem. Phys. 1967, 47, 488–498. [Google Scholar] [CrossRef]
  55. Wybourne, B. Symmetry Principles in Atomic Spectroscopy. J. Phys. 1970, 31, C4-33. [Google Scholar] [CrossRef]
  56. Mariwalla, K. Dynamical Symmetries in Mechanics. Phys. Rep. 1975, 20, 287–362. [Google Scholar] [CrossRef]
  57. Akyildiz, Y. On the dynamical symmetries of the Kepler problem. J. Math. Phys. 1980, 21, 665–670. [Google Scholar] [CrossRef]
  58. Fronsdal, C.; Huff, R. Two-Body Problem in Quantum Field Theory. Phys. Rev. D 1971, 3, 933. [Google Scholar] [CrossRef]
  59. Loebl, E. (Ed.) Group Theory and Its Applications; Academic Press: New York, NY, USA, 1971. [Google Scholar]
  60. Barut, A.; Rasmussen, W. The hydrogen atom as a relativistic elementary particle I. The wave equation and mass formulae. J. Phys. 1973, B6, 1695. [Google Scholar] [CrossRef]
  61. Barut, A.; Rasmussen, W. The hydrogen atom as a relativistic elementary particle II. Relativistic scattering problems and photo-effect. J. Phys. 1973, B6, 1713. [Google Scholar] [CrossRef]
  62. Barut, A.; Bornzin, G. Unification of the external conformal symmetry group and the internal conformal dynamical group. J. Math. Phys. 1974, 15, 1000–1006. [Google Scholar] [CrossRef]
  63. Barut, A.; Schneider, C.; Wilson, R. Quantum theory of infinite component fields. J. Math. Phys. 1979, 20, 2244–2256. [Google Scholar] [CrossRef]
  64. Shibuya, T.; Wulfman, C. The Kepler Problem in Two-Dimensional Momentum Space. Am. J. Phys. 1965, 33, 570–574. [Google Scholar] [CrossRef]
  65. Dahl, J. Physical Interpretation of the Runge-Lenz Vector. Phys. Let. 1968, 27A, 62–63. [Google Scholar] [CrossRef]
  66. Collas, P. Algebraic Solution of the Kepler Problem Using the Runge-Lenz Vector. Am. J. Phys. 1970, 38, 253–255. [Google Scholar] [CrossRef]
  67. Rodgers, H. Symmetry transformations of the classical Kepler problem. J. Math. Phys. 1973, 14, 1125–1129. [Google Scholar] [CrossRef]
  68. Majumdar, S.; Basu, D. O(3,1) symmetry of the hydrogen atom. J. Phys. Math. Nuc. Gen. 1974, 7, 787. [Google Scholar] [CrossRef]
  69. Stickforth, J. The classical Kepler problem in momentum space. Am. J. Phys. 1978, 46, 74–75. [Google Scholar] [CrossRef]
  70. Ligon, T.; Schaaf, M. On the Global Symmetry of the Classical Kepler Problem. Rep. Math. Phys. 1976, 9, 281–300. [Google Scholar] [CrossRef]
  71. Lakshmanan, M.; Hasegawa, H. On the canonical equivalence of the Kepler problem in coordinate and momentum space. J. Phys. Math. Gen. 1984, 17, L889. [Google Scholar] [CrossRef]
  72. O’Connell, R.; Jagannathan, K. Illustrating dynamical symmetries in classical mechanics: The Laplace-Runge-Lenz vector revisited. Am. J. Phys. 2003, 71, 243–246. [Google Scholar] [CrossRef]
  73. Valent, G. The hydrogen atom in electric and magnetic fields: Pauli’s 1926 article. Am. J. Phys. 2003, 71, 171–175. [Google Scholar] [CrossRef] [Green Version]
  74. Morehead, J. Visualizing the extra symmetry of the Kepler problem. Am. J. Phys. 2005, 73, 234–239. [Google Scholar] [CrossRef]
  75. Huntington, L.; Nooijen, M. An SO(4) invariant Hamiltonian and the two-body bound state. I: Coulomb interaction between two spinless particles. Int. J. Quant. Chem. 2009, 109, 2885–2896. [Google Scholar] [CrossRef]
  76. Barut, A.; Bohm, A.; Neeman, Y. Dynamical Groups and Spectrum Generating Algebras; World Scientific: Singapore, 1986. [Google Scholar]
  77. Greiner, W.; Muller, B. Quantum Mechanics, Symmetries; Springer: Berlin, Germany, 1989. [Google Scholar]
  78. Gilmore, R. Lie Groups, Lie Algegras and Some of Their Applications; Dover Books on Mathmatics; Dover: Mineola, NY, USA, 2005. [Google Scholar]
  79. Kibler, M. On the use of the group SO(4,2) in atomic and molecular physics. Mol. Phys. 2004, 102, 1221–1229. [Google Scholar] [CrossRef] [Green Version]
  80. Hammond, I.; Chu, S. Irregular wavefunction behavior in dimagnetic Rydberg atoms:a dynamical SO(4,2) group study. Chem. Phys. Let. 1991, 182, 63. [Google Scholar] [CrossRef]
  81. Lev, F. Symmetries in Foundation of Quantum Theory and Mathematics. Symmetry 2020, 12, 409. [Google Scholar] [CrossRef] [Green Version]
  82. Wulfman, C. Dynamical Symmetry; World Scientific Publishing: Singapore, 2011. [Google Scholar]
  83. Johnson, M.; Lippmann, B.A. Relativistic Kepler problem. Phys. Rev. 1950, 78, 329. [Google Scholar]
  84. Biedenharn, L. Remarks on the relativistic Kepler problem. Phys. Rev. 1962, 126, 845–851. [Google Scholar] [CrossRef]
  85. Lanik, J. The Reformulations of the Klein-Gordon and Dirac Equations for the Hydrogen Atom to Algebraic Forms. Czech. J. Phys. 1969, B19, 1540–1548. [Google Scholar] [CrossRef]
  86. Stahlhofen, A. Algebraic solutions of relativistic Coulomb problems. Helv. Phys. Acta 1997, 70, 1141. [Google Scholar]
  87. Chen, J.; Deng, D.; Hu, M. SO(4) symmetry in the relativistic hydrogen atom. Phys. Rev. A 2008, 77, 034102. [Google Scholar] [CrossRef] [Green Version]
  88. Khachidze, T.; Khelashvili, A. The hidden symmetry of the Coulomb problem in relativistic quantum mechanics: From Pauli to Dirac. Am. J. Phys. 2006, 74, 628–632. [Google Scholar] [CrossRef] [Green Version]
  89. Zhang, F.; Fu, B.; Chen, J. Dynamical symmetry of Dirac hydrogen atom with spin symmetry and its connection to Ginocchio’s oscillator. Phys. Rev. A 2008, 78, 040101(R). [Google Scholar] [CrossRef] [Green Version]
  90. Heine, V. Group theory in Quantum Mechanics; Dover Publications: New York NY, USA, 1993. [Google Scholar]
  91. Noether, E.; Mort, T. Invariant Variation Problems. Transp. Theory Stat. Phys. 1971, 1, 186–207. [Google Scholar] [CrossRef] [Green Version]
  92. Neuenschwander, D.E. Emmy Noether’s Wonderful Theorem; Johns Hopkins University Press: Baltimore, MD, USA, 2010. [Google Scholar]
  93. Hanca, J.; Tulejab, S.; Hancova, M. Symmetries and conservation laws: Consequences of Noether’s theorem. Am. J. Phys. 2004, 72, 428–435. [Google Scholar] [CrossRef] [Green Version]
  94. Byers, N.E. Noether’s Discovery of the Deep Connection Between Symmetries and Conservation Laws. arXiv 1998, arXiv:physics/9807044. [Google Scholar]
  95. The daughter of a mathematician, she wanted to be a mathematician, but since women were not allowed to take classes at the University of Erlingen, she audited courses. She did so well in the exams, that she received a degree and was allowed to enroll in the university and received a PhD in 1907. She remained at the university, unpaid, in an unofficial status, for 8 years. Then she went to the University at Gottengen, where she worked for 8 years with no pay or status before being appointed as Lecturer with a modest salary. She was invited in 1915 by Felix Klein and David Hilbert, two of the most famous mathematicians in the world at the time, to work with them and address issues in Einstein’s theory of General Relativity about energy conservation. She discovered Nother’s First Theorem (and a second theorem also). She remained there until 1933 when she, as a Jew, lost her job. At Einstein’s suggestion, she went to Bryn Mawr College in Pennsylvania. She died from ovarian cysts two years later.
  96. A rotation in 4 dimension can be represented by an antisymmetric 4 × 4 matrix which has 3 + 2 + 1 = 6 independent non-diagonal elements corresponding to 6 generators. Similarly a rotation in 5 dimensions has 10 independent elements or 10 generators.
  97. Dothan, Y. Finite-Dimensional Spectrum-Generating Algebras. Phys. Rev. 1970, D2, 2944. [Google Scholar] [CrossRef]
  98. Mukanda, N.; Sudarshan, E. Characteristic Noninvariance Groups of Dynamical Systems. Phys. Rev. Lett. 1965, 15, 1041. [Google Scholar] [CrossRef]
  99. Kyriakopoulos, E. Algebraic Equations for Bethe-Salpeter and Coulomb Green’s Functions. J. Math. Phys. 1972, 13, 1729–1735. [Google Scholar] [CrossRef]
  100. Lipkin, H. Lie Groups for Pedestrians; Dover Publications: New York, NY, USA, 2001. [Google Scholar]
  101. Were it not for this displacement of the force center, the observation that a rotated circle projects onto a plane as an ellipse would manifest the four-dimensional symmetry of the hydrogenlike atom directly in configuration space. The elliptical orbits could be viewed as projections of a rotated hypercircle onto a three-dimensional hyperplane. These considerations can be applied with some modification to the three-dimensional harmonic oscillator for which the force center and the center of the ellipse coincide.
  102. This equation and any other equation written in this specific coordinate system can be generalized to an arbitrary coordinate system by noting that the Cartesian unit vectors may be written in a manner that is independent of the coordinate system: i = A A ,   j = L × A L A ,   k = L L
  103. Brown, L.; University of Washington, Seattle, WA, USA. Unpublished lecture notes. 1972.
  104. We define the angle between a three-dimensional hyper-plane and a line as π/2 minus the angle between the line and the normal to the hyperplane.
  105. It is desirable to first show that A (and of course L) generate rotations of the hypersphere or U ^ . However, as we prefer to do the necessary calculations in terms of commutators rather than Poisson brackets, we defer these considerations to Section 4. There we show that the generator Li rotates U ^ about the i–4 plane; the generator A1 rotates U ^ about the 2–3 plane, etc., thereby changing the orbit with respect to the 4-axis and changing the eccentricity.
  106. Bois, G. Tables of Indefinite Integrals; Dover Pub1ications: New York, NY, USA, 1961; p. 123. [Google Scholar]
  107. Using Equation (44) and [102], Equation (73) may be written as cos−1(U · A/A) = p · r/arc + ωclt. This agrees with the time dependent function ϕ = p · r/arcωclt Equation (70) defined in [97].
  108. Brown, L. Forces giving no orbit precession. Am. J. Phys. 1978, 46, 930–931. [Google Scholar] [CrossRef]
  109. Bacry, H. Lectures in Theoretical Physics; Brittin, W.E., Barut, A.O., Guenin, M., Eds.; Gordon and Breach: New York, NY, USA, 1967. [Google Scholar]
  110. Barut, A. Dynamics of a Broken SUN Symmetry for the Oscillator. Phys. Rev. 1965, 139, B1433. [Google Scholar] [CrossRef]
  111. Boiteux, M. The Three-Dimensional Hydrogen Atom as a Restricted Four-Dimensional Harmonic Oscillator. Physica 1972, 65, 381–395. [Google Scholar] [CrossRef]
  112. Hughes, J. The harmonic oscillator:values of the SU(3) invariants. J. Phys. A Math. Gen. 1973, 6, 453. [Google Scholar] [CrossRef]
  113. Chen, A. Hydrogen atom as a four-dimensional oscillator. Phys. Rev. A 1980, 22, 333. [Google Scholar] [CrossRef]
  114. Chen, A. Homomorphism between SO(4,2) and SU(2,2). Phys. Rev. A 1981, 23, 1653. [Google Scholar] [CrossRef]
  115. Kibler, M.; Negadi, T. Connection between the hydrogen atom and the harmonic oscillator: The zero-energy case. Phys. Rev. A 1984, 29, 2891. [Google Scholar] [CrossRef]
  116. Chen, A.; Kibler, M. Connection between the hydrogen atom and the four-dimensional oscillator. Phys. Rev. A 1985, 31, 3960. [Google Scholar] [CrossRef]
  117. Gerry, C. Coherent states and the Kepler-Coulomb problem. Phys. Rev. A 1986, 33, 6. [Google Scholar] [CrossRef]
  118. Chen, A. Coulomb–Kepler problem and the harmonic oscillator. Am. J. Phys. 1987, 55, 250–252. [Google Scholar] [CrossRef]
  119. Van der Meer, J. The Kepler system as a reduced 4D oscillator. J. Geom. Phys. 2015, 92, 181–193. [Google Scholar] [CrossRef]
  120. Bacry, H. The de Sitter Group L4,1 and the Bound States of the Hydrogen Atom. Nuovo Cimento 1966, 41, 222–234. [Google Scholar] [CrossRef] [Green Version]
  121. Biedenharn, L. Wigner Coefficients for the R4 Group and Some Applications. J. Math. Phys. 1961, 2, 433–441. [Google Scholar] [CrossRef]
  122. Shiff, L. Quantum Mechanics; McGraw Hill: New York, NY, USA, 1955. [Google Scholar]
  123. Biedenharn, L.; Swamy, N. Remarks on the Relativistic Kepler Problem. II. Approximate Dirac-Coulomb Hamiltonian Possessing Two Vector Invariants. Phys. Rev. 1964, 133, B1353. [Google Scholar] [CrossRef]
  124. Morse, P.; Feshbach, H. Methods of Theoretical Physics, Vol. 1; McGraw-Hill: New York, NY, USA, 1953. [Google Scholar]
  125. The primes indicates eigenvalues of operators, and unprimed quantities indicate abstract operators. The quantity x′ means the four-vector (t′, r )
  126. Morse, P.; Feshbach, H. Methods of Theoretical Physics, Vol. 2; McGraw-Hill: New York, NY, USA, 1953. [Google Scholar]
  127. Erdeli, A. (Ed.) Higher Transcendental Functions, Bateman Manuscript Project; McGraw-Hill Book Co.: New York, NY, USA, 1953. [Google Scholar]
  128. Makowski, A.; Pepłowski, P. Zero-energy wave packets that follow classical orbits. Phys. Rev. A 2012, 86, 042117. [Google Scholar] [CrossRef]
  129. Bellomo, P.; Stroud, C., Jr. Classical evolution of quantum elliptical orbits. Phys. Rev. A 1999, 59, 2139. [Google Scholar] [CrossRef] [Green Version]
  130. Berry, M.; Mount, K.E. Semiclassical approximations in wave mechanics. Rep. Prog. Phys. 1972, 35, 315. [Google Scholar] [CrossRef] [Green Version]
  131. Eberly, J.; Stroud, C. Chapters 14 (Rydberg Atoms) and Chapter 73 (Coherent Transients). In Springer Handbook of Atomic, Molecular, and Optical Physics; Drake, G., Ed.; Springer Science and Business Media: New York, NY, USA, 2006. [Google Scholar]
  132. Lakshmanan, M.; Ganesan, K. Rydberg atoms and molecules-Testing grounds for quantum manifestations of chaos. Curr. Sci. 1995, 68, 38–44. [Google Scholar]
  133. Kay, K. Exact Wave Functions for the Coulomb Problem from Classical Orbits. Phys. Rev. 1999, 25, 5190. [Google Scholar] [CrossRef]
  134. Lena, C.; Deland, D.; Gay, J. Wave functions of Atomic Elliptic States. Europhys. Lett. 1991, 15, 697. [Google Scholar] [CrossRef]
  135. Bhaumik, D.; Dutta-Roy, B.; Ghosh, G. Classical limit of the hydrogen atom. J. Phys. A Math. Gen. 1986, 19, 1355. [Google Scholar] [CrossRef]
  136. McAnally, D.; Bracken, A. Quasiclassical states of the Coulomb system and SO(4, 2). J. Phys. A Math. Gen. 1990, 23, 2027. [Google Scholar] [CrossRef]
  137. Pitak, A.; Mostowski, J. Classical limit of position and matrix elements for Rydberg atoms. Eur. J. Phys. 2018, 39, 025402. [Google Scholar] [CrossRef]
  138. Nauenberg, M. Quantum wavepackets on Kepler elliptical orbits. Phys. Rev. A 1989, 40l, 1133. [Google Scholar] [CrossRef] [PubMed]
  139. Brown, L.S. Classical limit of the hydrogen atom. Am. J. Phys. 1973, 41, 525–530. [Google Scholar] [CrossRef]
  140. Leonhardt, U. Measuring the Quantum State of Light; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  141. Barut, A.; Girardello, L. New “coherent” states associated with non-compact groups. Commun. Math. Phys. 1971, 21, 41–55. [Google Scholar] [CrossRef]
  142. Satyanarayana, M. Squeezed coherent states of the hydrogen atom. J. Phys. A Math. Gen. 1986, 19, 1973. [Google Scholar] [CrossRef]
  143. Liu, Q.; Hu, B. The hydrogen atom’s quantum-to-classical correspondence in Heisenberg’s correspondence principle. J. Phys. A Math. Gen. 2001, 34, 5713. [Google Scholar] [CrossRef]
  144. Zverev, V.; Rubinstein, B. Dynamical symmetries and well-localized hydrogenic wave packets. Proc. Inst. Math. Nas Ukr. 2004, 50, 1018. [Google Scholar]
  145. The wave function in momentum space ψ(p) is obtained by multiplying Ynlm by the normalizing factor ( a n ) 3 / 2 ( 1 U 4 ) 2 , cf Equation (150).
  146. Nandi, S.; Shastry, C. Classical limit of the two-dimensional and three-dimensional hydrogen atom. J. Phys. A Math. Gen. 1989, 22, 1005. [Google Scholar] [CrossRef] [Green Version]
  147. Pauling, L.; Wilson, E.B. Introduction to Quantum Mechanics; McGraw-Hill: New York, NY, USA, 1935; p. 36. [Google Scholar]
  148. Lamb, W.; Retherford, R. Fine Structure of the H Atom, Part I. Phys. Rev. 1950, 79, 549. [Google Scholar] [CrossRef]
  149. Bethe, H.; Salpeter, E. The Quantum Mechanics of One and Two Electron Atoms; Springer: Berlin, Germany, 1957. [Google Scholar]
  150. Milonni, P. The Quantum Vacuum; Academic Press: San Diego, CA, USA, 1994. [Google Scholar]
  151. Eides, M.; Grotch, H.; Shelyuto, V. Theory of Light Hydrogenic Bound States, Springer Tracts in Modern Physics 222; Springer: Berlin, Germany, 2007. [Google Scholar]
  152. Rau, A.; Alber, G. Shared symmetries of the hydrogen atom and the two-bit system. J. Phys. B At. Mol. Opt. 2017, 50, 242001. [Google Scholar] [CrossRef] [Green Version]
  153. Castro, P.; Kullock, R. Physics of the SOp(4) Hydrogen Atom. Theo. Math. Phys. 2015, 185, 1678. [Google Scholar] [CrossRef] [Green Version]
  154. Alavi, A.; Rezaei, N. Dirac equation, hydrogen atom spectrum and the Lamb shift in dynamical non-commutative spaces. Pramana-J. Phys. 2017, 88, 5. [Google Scholar] [CrossRef]
  155. Gnatenko, K.P.; Krynytskyi, Y.S.; Tkachuk, V.M. Perturbation of the ns levels of the hydrogen atom in rotationally invariant noncommutative space. Mod. Phys. Lett. 2015, 30, 8. [Google Scholar] [CrossRef] [Green Version]
  156. Haghighat, M.; Khorsandi, M. Hydrogen and muonic hydrogen atomic spectra in non-commutative space-time. Eur. Phys. J. 2015, 75, 1. [Google Scholar] [CrossRef] [Green Version]
  157. Praxmeyer, L. Hydrogen atom in phase space: The Wigner representation. J. Phys. A Math. Gen. 2006, 39, 14143. [Google Scholar] [CrossRef]
  158. Jones, M.; Potvliege, M.R.; Spannowsky, M. Probing new physics using Rydberg states of atomic hydrogen. Phys. Rev. Res. 2020, 2, 013244. [Google Scholar] [CrossRef] [Green Version]
  159. Jentschura, U.; Mohr, P. Calculation of hydrogenic Bethe logarithms for Rydberg States. Phys. Rev. A 2002, 72, 012110. [Google Scholar] [CrossRef] [Green Version]
  160. Jentschura, U.; LeBigot, E.; Evers, J.; Mohr, P.; Keitel, C. Relativistic and radiative shifts for Rydberg states. J. Phys. B At. Mol. Opt. Phys. 2005, 38, S97. [Google Scholar] [CrossRef] [Green Version]
  161. Jentschura, U.; Mohr, P.; Tan, J. Fundamental constants and tests of theory in Rydberg states of one-electron ions. J. Phys. B At. Mol. Opt. Phys. 2010, 43, 074002. [Google Scholar] [CrossRef]
  162. Cantu, S.H.; Venkatramani, A.V.; Xu, W. Repulsive photons in a quantum nonlinear medium. Nat. Phys. 2020. [Google Scholar] [CrossRef]
Figure 1. Classical Kepler orbit in configuration space. The orbit is in the 1–2 plane (plane of the paper). One focus, where the proton charge is located, is the origin. The semi-minor axis is b = r c sin ν . The semi-major axis is r c .
Figure 1. Classical Kepler orbit in configuration space. The orbit is in the 1–2 plane (plane of the paper). One focus, where the proton charge is located, is the origin. The semi-minor axis is b = r c sin ν . The semi-major axis is r c .
Symmetry 12 01323 g001
Figure 2. Kepler bound state orbits in configuration space for a fixed energy and different values of the eccentricity e = sin ν . The bold hypotenuse is the semi-major axis r c , which makes an angle ν with the vertical 2-axis.
Figure 2. Kepler bound state orbits in configuration space for a fixed energy and different values of the eccentricity e = sin ν . The bold hypotenuse is the semi-major axis r c , which makes an angle ν with the vertical 2-axis.
Symmetry 12 01323 g002
Figure 3. Kepler orbit in momentum space of radius a / cos ν , with its center at p 2 = a tan ν , corresponding to the orbit in configuration space shown in Figure 1. A circular orbit in configuration space corresponds to a circular orbit in momentum space centered on the origin with radius a.
Figure 3. Kepler orbit in momentum space of radius a / cos ν , with its center at p 2 = a tan ν , corresponding to the orbit in configuration space shown in Figure 1. A circular orbit in configuration space corresponds to a circular orbit in momentum space centered on the origin with radius a.
Symmetry 12 01323 g003
Figure 4. Kepler orbits in momentum space for a fixed energy and RMS momentum a with different values of the eccentricity e = sin ν corresponding to the orbits in configuration space shown in Figure 2.
Figure 4. Kepler orbits in momentum space for a fixed energy and RMS momentum a with different values of the eccentricity e = sin ν corresponding to the orbits in configuration space shown in Figure 2.
Symmetry 12 01323 g004
Figure 5. Stereographic projection in momentum space for a fixed energy, mapping p / a into U ^ . The unit vector in the four direction is n ^ and n ^ · U ^ = cos Θ 4 .
Figure 5. Stereographic projection in momentum space for a fixed energy, mapping p / a into U ^ . The unit vector in the four direction is n ^ and n ^ · U ^ = cos Θ 4 .
Symmetry 12 01323 g005
Figure 6. Showing the hyperplane containing the orbit making an angle ν with the U 2 plane. Notice tan ν = A / L as required by Equation (30).
Figure 6. Showing the hyperplane containing the orbit making an angle ν with the U 2 plane. Notice tan ν = A / L as required by Equation (30).
Symmetry 12 01323 g006
Figure 7. Model illustrating the stereographic projection from the 1–2 plane to a 4-d hypersphere. The pin represents the unit vector n ^ along the 4-axis, normal to the 1–2 plane.
Figure 7. Model illustrating the stereographic projection from the 1–2 plane to a 4-d hypersphere. The pin represents the unit vector n ^ along the 4-axis, normal to the 1–2 plane.
Symmetry 12 01323 g007
Figure 8. Wave function probability distribution | Y n , n 1 , n 1 ( p ) | 2 in momentum space for large n, showing the variation in the width of the momentum distribution about the classical circular orbit. The center of the distribution is at p 2 = a tan ν . The classical orbit is in the 1–2 plane.
Figure 8. Wave function probability distribution | Y n , n 1 , n 1 ( p ) | 2 in momentum space for large n, showing the variation in the width of the momentum distribution about the classical circular orbit. The center of the distribution is at p 2 = a tan ν . The classical orbit is in the 1–2 plane.
Symmetry 12 01323 g008
Table 1. In the table L = rxp is the orbital angular momentum; A is the Runge-Lenz vector; J = L +  σ /2 is the total angular momentum; K is the generalized parity operator; Λ is the conserved pseudoscalar operator.
Table 1. In the table L = rxp is the orbital angular momentum; A is the Runge-Lenz vector; J = L +  σ /2 is the total angular momentum; K is the generalized parity operator; Λ is the conserved pseudoscalar operator.
Degeneracy Groups for Bound States in a Coulomb Potential
EquationDegeneracyConserved QuantitiesDegeneracy GroupRepresentationDimension
SchrodingerE indep. of l , l z A , L S O ( 4 ) ( n 1 2 , n 1 2 ) n 2
Klein-GordonE indep. of l z L O ( 3 ) Casimir op. is l ( l + 1 ) 2 l + 1
Klein-Gordon without V 2 termE indep. of l , l z A , L S O ( 4 ) ( n 1 2 , n 1 2 ) n 2
DiracE depends on J, n only Λ , K , J S O ( 4 ) ( 1 / 2 , J ) 2(2J+1)

Share and Cite

MDPI and ACS Style

Maclay, G.J. Dynamical Symmetries of the H Atom, One of the Most Important Tools of Modern Physics: SO(4) to SO(4,2), Background, Theory, and Use in Calculating Radiative Shifts. Symmetry 2020, 12, 1323. https://doi.org/10.3390/sym12081323

AMA Style

Maclay GJ. Dynamical Symmetries of the H Atom, One of the Most Important Tools of Modern Physics: SO(4) to SO(4,2), Background, Theory, and Use in Calculating Radiative Shifts. Symmetry. 2020; 12(8):1323. https://doi.org/10.3390/sym12081323

Chicago/Turabian Style

Maclay, G. Jordan. 2020. "Dynamical Symmetries of the H Atom, One of the Most Important Tools of Modern Physics: SO(4) to SO(4,2), Background, Theory, and Use in Calculating Radiative Shifts" Symmetry 12, no. 8: 1323. https://doi.org/10.3390/sym12081323

APA Style

Maclay, G. J. (2020). Dynamical Symmetries of the H Atom, One of the Most Important Tools of Modern Physics: SO(4) to SO(4,2), Background, Theory, and Use in Calculating Radiative Shifts. Symmetry, 12(8), 1323. https://doi.org/10.3390/sym12081323

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop