Just over thirty years ago, a new generation of quarks was discovered when Fermilab announced they had found the bottom quark [1], adding to the known up, down, strange, and charm quarks. The discovery was indirect—the actual detection involved finding bottom-antibottom quark pairs ($b\phantom{\rule{0}{0ex}}\overline{b}$) that form bound states via strong interactions and have a rich spectroscopy analogous to that of the hydrogen atom [2]. These composite particles are called bottomonium, an analogy to the well-known electron-positron pairs called positronium. The first two $b\phantom{\rule{0}{0ex}}\overline{b}$ states that were discovered are named upsilon particles ($\mathrm{{\rm Y}}$ and $\mathrm{{\rm Y}}\prime $) and were found in 1977 during experiments with collisions of 400-GeV protons on nuclear targets at Fermilab [1]. Subsequently, a variety of other excited states (all spin triplets) have been observed.

However, no spin-singlet state had been seen until the observation of the ground state called ${\eta}_{b}$, now reported in *Physical Review Letters* [3] by the BABAR collaboration. The difference in mass between the $\mathrm{{\rm Y}}$ and the ${\eta}_{b}$ is important in understanding quark-antiquark states (generally called quarkonia) by testing existing models, the applicability of perturbative quantum chromodynamics to the $b\phantom{\rule{0}{0ex}}\overline{b}$ system, and the results of the numerical approach, known as lattice quantum chromodynamics (lattice QCD), to calculate hadron properties [4]. More importantly, having a measured value will challenge theorists to perform more precise calculations that can be compared to experiment.

Heavy quarkonia, which are bound states of a heavy quark and antiquark, are well described by nonrelativistic potential models originally derived to describe charm-anticharm ($c\phantom{\rule{0}{0ex}}\overline{c}$) states [5, 6]. The potentials incorporate general features of quantum chromodynamics (QCD)—the theory of quarks and gluons describing the strong interactions. At short distances, these QCD-motivated potentials take the form of a one-gluon exchange potential, analogous to the photon exchange that is responsible for the Coulomb interaction in quantum electrodynamics (QED). Added to this are relativistic corrections, such as spin-spin and spin-orbit terms, all with “color” factors reflecting the more complicated group structure of QCD compared to QED. The spin-spin term, for example, is analogous to the hyperfine interaction that gives rise to the 21-cm line in hydrogen. At large separation the potential is described by a linearly rising interaction that confines the quarks. The QCD-motivated phenomenological potential is in good agreement with results obtained using numerical lattice-QCD methods [4]. Lattice QCD is a nonperturbative approach that deals with the nonlinear nature of the strong interaction by dividing space and time into discrete grid points and then integrating over quark and gluon configurations.

In these potential models, quarkonium energy levels are found by solving a nonrelativistic Schrödinger equation, although more sophisticated calculations take into account relativistic corrections [7]. The calculations yield energy levels that are characterized by the radial quantum number $n$, which is equal to one plus the number of nodes of the radial wave function, and $L$, the relative orbital angular momentum between the quark and antiquark. In fact, much of the nomenclature is familiar from atomic physics. The orbital levels are labeled by $S$, $P$, $D$ (corresponding to $L=0,1,2$). The spins of the quark and antiquark couple to give total spin $S=0$ (spin-singlet) or $S=1$ (spin-triplet) states. $S$ and $L$ couple to give the total angular momentum of the state $J$, which can take on values $J=L-1$, $L$, or $L+1$. Thus the $L=0$ states are ${}^{1}{S}_{0}$ and ${}^{3}{S}_{1}$; the $L=1$ states are ${}^{1}{P}_{1}$ and ${}^{3}{P}_{0}$, ${}^{3}{P}_{1}$, ${}^{3}{P}_{2}$, etc.

In addition to the spin-independent potential, there are spin-dependent interactions that give rise to splittings within multiplets [7]. With these, we can predict $\mathrm{{\rm Y}}-{\eta}_{b}$ splittings in bottomonium and similar splittings in charmonium that are analogous to the hyperfine splittings in hydrogen. Splittings within $P$-wave and higher $L$-state multiplets are due to spin-orbit and spin-spin interactions arising from one-gluon exchange and a relativistic spin-orbit precession term. The contact spin-spin splitting between the singlet and triplet $P$-wave states is predicted to be small due to its short range and because the wavefunction at the origin for $P$-wave states is zero.

The observations of the ${\eta}_{b}$ by BABAR [3] and the charmonium state ${h}_{c}$ by CLEO [8] are important validations of this picture. In the experiment reported by Aubert *et al.*, electrons and positrons from the PEP-II storage ring at SLAC collide with a center-of-mass energy of 10.355 GeV. This energy is selected so that the collisions create $\mathrm{{\rm Y}}(3S)$ particles, some of which then decay radiatively to the ${\eta}_{b}(1S)$ state. Figure 1 shows the $b\phantom{\rule{0}{0ex}}\overline{b}$ spectrum of observed states along with predictions for missing states [9] by Isgur and myself. The commonly used names of observed levels are shown. Note that $b\phantom{\rule{0}{0ex}}\overline{b}$ states with mass greater than two times the mass of a $B$ meson (the ground state of a meson made up of a bottom-quark and a light up or down quark), will have a large decay rate into $B-\overline{B}$ pairs so the branching ratio for radiative decays will be small.

Electromagnetic transitions between the levels can be calculated in the quark model and provide an important tool in understanding the quarkonium internal structure [10]. The theory and terminology of electromagnetic transitions between quarkonium states closely follows the treatment given for transitions in the hydrogen atom in undergraduate quantum mechanics textbooks with the replacement of the electric charge of the electron with that of the quark charge and one has to include both the quark and antiquark transition amplitudes. The leading-order transition amplitudes are due to electric dipole transitions ($E1$) between states with the same total spin and magnetic dipole transitions ($M1$) which flip the quark or antiquark spin and are inversely proportional to the constituent quark mass. The predictions for $E1$ transitions, ${}^{3}{P}_{J}{\leftrightarrow}^{3}{S}_{1}$, in the bottomonium system are in good agreement with experimental data [10]. Recently, the CLEO experiment observed a $1D$ bottomonium state in a cascade of $E1$ transitions with a mass of $10161.1\pm 0.6\pm 1.6\phantom{\rule{0.333em}{0ex}}\text{MeV}/{\text{c}}^{2}$ [11], which is in good agreement with theoretical predictions [7].

In the nonrelativistic limit the spatial overlap integrals for $M1$ transitions equal one between $S$-wave states within the same multiplet (that is, they are favored transitions) and zero for transitions between states with different radial quantum numbers (that is, these transitions are hindered). Relativistic corrections leads to small overlaps for these hindered transitions, which can be compensated by large phase-space factors [12]. Until the observation of the hindered $\mathrm{{\rm Y}}(3S)$ to ${\eta}_{b}$ transition, no $M1$ transitions had been observed in the bottomonium system.

Until now, all of the observed states in the bottomonium system were spin-triplet states but quark models predict the existence of spin-singlet partners including the ground state. As mentioned above, while the decay amplitudes for hindered transitions are much smaller than those for favored transitions, this can be compensated with the larger available phase space in transitions such as $\mathrm{{\rm Y}}(3S)\to {\eta}_{b}(1S)$. BABAR collected a large data set by tuning the ${e}^{+}\phantom{\rule{0}{0ex}}{e}^{-}$ energy to the mass of the $\mathrm{{\rm Y}}(3S)$ and observed a signal in the photon energy with ${E}_{\gamma}=921.2+2.1/-2.8(\text{stat})\pm 2.4(\text{syst})\text{MeV}$ where the first error is statistical and the second systematic, which they interpreted as an $M1$ transition to the ${\eta}_{b}(1S)$ [3]. This corresponds to an ${\eta}_{b}(1S)$ mass of $9388.9+3.1/-2.3(\text{stat})\pm 2.7(\text{syst})\text{MeV}/{\text{c}}^{2}$ with corresponding $\mathrm{{\rm Y}}(1S)-{\eta}_{b}$ hyperfine mass splitting of $71.4+2.3/-3.1(\text{stat})\pm 2.7(\text{syst})\text{MeV}/{\text{c}}^{2}$.

The measured $\mathrm{{\rm Y}}(3S)-{\eta}_{b}$ splitting is consistent with potential model predictions although a significant subset of predictions lie outside experimental one-sigma error bounds [12]. A recent lattice-QCD calculation predicts a value of $61\pm 14\phantom{\rule{0.333em}{0ex}}\text{MeV}/{\text{c}}^{2}$, which is consistent within the large errors [4]. Two recent calculations using a perturbative-QCD approach predict splittings of $39\pm 14\phantom{\rule{0.333em}{0ex}}\text{MeV}/{\text{c}}^{2}$ [13] and $44\pm 11\phantom{\rule{0.333em}{0ex}}\text{MeV}/{\text{c}}^{2}$ [14], both being over two standard deviations away from the BABAR measurement. One can see that the recent BABAR result poses a serious challenge to theorists, which should spur renewed effort to improve calculations. More precise measurement of the ${\eta}_{b}$ mass would allow precision tests of lattice-QCD and perturbative-QCD calculations of the $\mathrm{{\rm Y}}-{\eta}_{b}$ splitting.

The large amount of data that BABAR has accumulated on the $\mathrm{{\rm Y}}(3S)$ state should allow searches for other missing $b\phantom{\rule{0}{0ex}}\overline{b}$ states. In particular, it may be possible to observe the ${\eta}_{b}(2S)$ state via $M1$ transitions. Many models predict the branching ratio to the ${\eta}_{b}(2S)$ to be only a factor of 2 or 3 smaller than that to the ${\eta}_{b}(1S)$ and therefore possibly observable. Other interesting possibilities consist of searching for the ${h}_{b}({1}^{1}{P}_{1}$) in the processes $\mathrm{{\rm Y}}(3S)\to {\pi}^{0}\phantom{\rule{0}{0ex}}{h}_{b}({1}^{1}{P}_{1})\to {\pi}^{0}\phantom{\rule{0}{0ex}}\gamma \phantom{\rule{0}{0ex}}{\eta}_{b}$ and the sequential process $\mathrm{{\rm Y}}(3S)\to {\pi}^{+}\phantom{\rule{0}{0ex}}{\pi}^{-}\phantom{\rule{0}{0ex}}{h}_{b}({1}^{1}{P}_{1})\to {\pi}^{+}\phantom{\rule{0}{0ex}}{\pi}^{-}\phantom{\rule{0}{0ex}}\gamma \phantom{\rule{0}{0ex}}{\eta}_{b}$ [15]. The discovery of these states would represent an important step in completing the bottomonium spectrum and provide an important test of QCD-based models and calculations. Measurement of the hyperfine mass splittings between the triplet and singlet quarkonium states is crucial to understanding the role of spin-spin interactions in quarkonium models and in testing QCD calculations [7].

## Acknowledgments

The author gratefully acknowledges Nathan Isgur and Jon Rosner for teaching him much of what he knows about this subject. This research was supported in part by the Natural Sciences and Engineering Research Council of Canada.

### References

- S. W. Herb et al., Phys. Rev. Lett.
**39**, 252 (1977). - W. Kwong, J. L. Rosner, and C. Quigg, Annu. Rev. Nucl. Part. Sci.
**37**, 325 (1987). - B. Aubert et al. (), Phys. Rev. Lett.
**101**, 071801 (2008). - A. Gray, I. Allison, C. T. H. Davies, E. Gulez, G. P. Lepage, J. Shigemitsu, and M. Wingate (), Phys. Rev. D
**72**, 094507 (2005). - E. Eichten and K. Gottfried, Physics Letters B
**66**, 286 (1977). - W. Celmaster, H. Georgi, and M. Machacek, Phys. Rev. D
**17**, 886 (1978). - N. Brambilla, et al. (), arXiv:hep-ph/0412158.
- J. L. Rosner et al. (), Phys. Rev. Lett.
**95**, 102003 (2005). - S. Godfrey and N. Isgur, Phys. Rev. D
**32**, 189 (1985). - E. Eichten, S. Godfrey, H. Mahlke, and J. L. Rosner, arXiv:hep-ph/0701208.
- G. Bonvicini et al. (), Phys. Rev. D
**70**, 032001 (2004). - S. Godfrey and J. L. Rosner, Phys. Rev. D
**64**, 074011 (2001);**65**, 039901(E) (2002). - B. A. Kniehl, A. A. Penin, A. Pineda, V. A. Smirnov, and M. Steinhauser, Phys. Rev. Lett.
**92**, 242001 (2004). - S Recksiegel and Y. Sumino, Phys. Lett. B
**578**, 369 (2004). - S. Godfrey and J. L. Rosner, Phys. Rev. D
**66**, 014012 (2002). - W.-M. Yao et al., J. Phys. G: Nucl. Part. Phys.
**33**, 1 (2006).