How to use this course
Hans Bethe worked as a theoretical physicist for more than seventy-five years. His first paper appeared in 1928 and his last research papers in the 2000s. In between he contributed to solid-state physics, atomic physics, nuclear physics, quantum electrodynamics, astrophysics and neutrino physics. His style was distinctive. He picked a concrete problem, chose the simplest model that captured it, and pushed the calculation through to a number that could be compared with experiment. Colleagues called the method “Bethe’s way”, and several of his papers became the standard reference for decades.
- Quantum mechanics at the level of a first graduate course: perturbation theory, the hydrogen atom, scattering, spin.
- Special relativity and basic nuclear physics (binding energy, beta decay).
- For Lessons 08 and 09, some exposure to quantum field theory or to second-order perturbation theory with a continuum of states.
- For Lessons 07 and 10, the elements of stellar structure are helpful but not required.
Each technical lesson opens with a paragraph for the general reader, then gives precise statements, derivations and exercises with folded solutions. The lessons follow the order of his work: crystals and stopping power (1928–1932), exactly solvable chains (1931), the nucleus (1934–1937), stellar energy (1938–1939), quantum electrodynamics (1947), bound states and nuclear matter (1951–1970), and supernovae and neutrinos (1979–2000).
Timeline
| 1906 | Born 2 July in Strasbourg, then in the German Empire, only child of the physiologist Albrecht Bethe and Anna Kuhn. |
| 1912–15 | The family moves to Kiel, then to Frankfurt. |
| 1924 | Enters the University of Frankfurt. |
| 1926 | Moves to Munich to study with Arnold Sommerfeld. |
| 1928 | July: doctorate on the diffraction of electrons by crystals. Short posts at Frankfurt and Stuttgart (with Paul Ewald). |
| 1929 | Crystal field theory: the splitting of atomic levels in crystals. |
| 1930 | Theory of the passage of fast particles through matter. Privatdozent at Munich. Rockefeller fellowship: Cambridge. |
| 1931 | Rome with Enrico Fermi; exact solution of the one-dimensional Heisenberg chain. |
| 1932–33 | Acting assistant professor at Tübingen. Dismissed in 1933 under the Nazi civil-service law. |
| 1933–34 | Lecturer at Manchester; with Rudolf Peierls explains the photodisintegration of the deuteron. Fellow at Bristol in autumn 1934. |
| 1935 | February: assistant professor at Cornell University. Professor from 1937. |
| 1936–37 | The three Reviews of Modern Physics articles known as “Bethe’s Bible”. |
| 1938 | Washington conference on stellar energy; proton–proton reaction with Charles Critchfield; the CN cycle. |
| 1939 | Energy Production in Stars, Physical Review. 13 September: marries Rose Ewald. |
| 1941–42 | American citizen. Radar work at the MIT Radiation Laboratory. |
| 1943–46 | Head of the Theoretical Division at Los Alamos. |
| 1947 | Shelter Island conference; nonrelativistic calculation of the Lamb shift. |
| 1948 | Name added by George Gamow to the Alpher–Bethe–Gamow paper. |
| 1951 | Bethe–Salpeter equation. Returns to thermonuclear weapons work after the Teller–Ulam design. |
| 1954 | Testifies for J. Robert Oppenheimer at his security hearing. |
| 1957–58 | President’s Science Advisory Committee; chairs a panel on verifying a test ban; Geneva conference of experts. |
| 1961 | Enrico Fermi Award. |
| 1967 | Nobel Prize in Physics. |
| 1968 | With Richard Garwin, Scientific American article against anti-ballistic-missile systems. |
| 1979 | Begins collaboration with Gerald Brown on stellar collapse. |
| 1985–86 | Delayed neutrino-heating supernova mechanism with James Wilson; MSW explanation of the solar-neutrino deficit. |
| 1995 | Calls on scientists to stop work on new nuclear weapons. |
| 2005 | Dies 6 March at home in Ithaca, aged 98. |
Strasbourg, Kiel, Frankfurt, 1906–1926
Hans Albrecht Bethe was born on 2 July 1906 in Strasbourg, then the capital of the German imperial territory of Alsace-Lorraine. His father, Albrecht Bethe, was a physiologist at the university. His mother, Anna Kuhn, was the daughter of a Strasbourg professor of medicine. Her family was Jewish; Hans was raised as a Protestant, like his father. He was an only child.
The boy loved numbers. Family stories describe him computing with fractions and powers at an age when other children were learning to read, and filling notebooks with tables. In 1912 the family moved to Kiel, where his father directed the Institute of Physiology, and in 1915 to Frankfurt, where Albrecht Bethe took a chair at the new university. Hans entered the Goethe-Gymnasium, but tuberculosis in 1916 interrupted his schooling. He spent periods at a school in Bad Kreuznach and at the progressive Odenwaldschule before returning to the Goethe-Gymnasium, from which he graduated in 1924. By his own account he taught himself calculus as a teenager.
At the University of Frankfurt he began in chemistry and mathematics. He found the experimental work clumsy and the physics teaching weak, apart from the lectures of Walter Gerlach, who left in 1925. A lecturer, Karl Meissner, advised him that the place to learn theoretical physics was Munich. In April 1926 he moved there.
Sommerfeld’s Munich and the wander years, 1926–1933
Arnold Sommerfeld’s institute in Munich had already trained Werner Heisenberg, Wolfgang Pauli and Peter Debye. Bethe arrived in the year Erwin Schrödinger published wave mechanics, and he learned the new quantum mechanics as it was being made. Sommerfeld set him a problem prompted by the experiments of Clinton Davisson and Lester Germer on electron diffraction. Bethe’s thesis, Theorie der Beugung von Elektronen an Kristallen, earned a summa cum laude in July 1928. It introduced the dynamical theory of electron diffraction and showed that the mean electrostatic potential inside a crystal shifts the Bragg peaks.
He then held one-semester posts at Frankfurt and at the Technische Hochschule Stuttgart, where he worked with Paul Ewald, the X-ray crystallographer whose daughter Rose he would later marry. Back in Munich he wrote two papers that would each found a subject: crystal field theory (1929) and the theory of stopping power (1930). He became a Privatdozent at Munich in 1930.
A Rockefeller Foundation fellowship took him to Cambridge in the autumn of 1930, where he worked in the circle of Ralph Fowler, and to Rome in the spring terms of 1931 and 1932. In Rome he learned from Enrico Fermi a lighter, more intuitive style, one that looked for the essential physics before the mathematics. He later said that Sommerfeld taught him rigour and Fermi taught him simplicity. In Rome he also solved the one-dimensional Heisenberg chain exactly (Lesson 05). While in Cambridge he wrote, with Guido Beck and Wolfgang Riezler, a mock paper deriving absolute zero from the fine-structure constant, a parody of Arthur Eddington’s numerology. Die Naturwissenschaften printed it, and the authors later had to publish a retraction.
With Sommerfeld he wrote the long Handbuch der Physik article on the electron theory of metals (1933), and alone the article on one- and two-electron atoms, which remained a standard reference for decades. In the winter of 1932–33 he became acting assistant professor at Tübingen, under Hans Geiger.
Exile and Cornell, 1933–1938
In April 1933 the Nazi Law for the Restoration of the Professional Civil Service removed from state employment anyone with a Jewish grandparent. Bethe’s mother came from a Jewish family, and he lost his post at Tübingen. Sommerfeld found him a fellowship in Munich for the summer. William Lawrence Bragg offered a temporary lectureship at Manchester, and in October 1933 Bethe left Germany.
In Manchester he shared lodgings with Rudolf Peierls, another refugee theorist. When James Chadwick and Maurice Goldhaber disintegrated the deuteron with gamma rays in 1934, Bethe and Peierls worked out the theory within weeks (Lesson 06). In the same year Bethe and Walter Heitler published the quantum theory of bremsstrahlung and pair production, which became the basis of the theory of cosmic-ray showers. After a semester at Bristol with Nevill Mott, he accepted an offer from Cornell University and arrived in Ithaca in February 1935. He remained at Cornell for the rest of his life.
Cornell had little theoretical physics. Bethe built a programme by teaching, and by writing. With Robert Bacher and Milton Stanley Livingston he produced three long articles in Reviews of Modern Physics (1936–37) that surveyed all of nuclear physics, much of it rederived by Bethe himself. A generation of physicists learned the subject from “Bethe’s Bible”. In 1938 George Gamow and Edward Teller persuaded a reluctant Bethe to attend a small Washington conference on the energy of stars. Within months he had answered the question (Lesson 07).
Crystals and fast particles
For the general reader. Put an atom inside a crystal and its energy levels split, because the surrounding ions break the symmetry of empty space. Shoot a fast charged particle into a block of matter and it slows down by kicking electrons out of their atoms. Bethe gave the theory of both effects between 1929 and 1932. The first underlies the colours of gemstones and the magnetism of transition-metal compounds. The second is still used every day in radiation therapy, detector design and nuclear physics.
- See how group theory predicts the splitting of atomic levels in a crystal.
- Derive the Bethe formula for stopping power in the nonrelativistic limit.
- Understand the role of the mean excitation energy \(I\).
Crystal field theory (1929)
An ion in a crystal feels an electrostatic potential with the point symmetry \(G\) of its site, a finite subgroup of \(\mathrm{O}(3)\). A free-ion level of angular momentum \(l\) carries the \((2l+1)\)-dimensional representation \(D^{(l)}\) of \(\mathrm{SO}(3)\). Restricted to \(G\), it is in general reducible.
If the crystal potential is a perturbation of the free-ion Hamiltonian with site symmetry \(G\), a level of angular momentum \(l\) splits into sublevels labelled by the irreducible representations \(\Gamma_i\) of \(G\) that occur in \(D^{(l)}\!\restriction_G=\bigoplus_i m_i\Gamma_i\). Each sublevel has degeneracy \(\dim\Gamma_i\) (barring accidental degeneracy), and the multiplicities follow from characters: \(m_i=\frac{1}{|G|}\sum_{g\in G}\overline{\chi_{\Gamma_i}(g)}\,\chi^{(l)}(g)\), with \(\chi^{(l)}(\text{rotation by }\theta)=\dfrac{\sin\big((l+\tfrac12)\theta\big)}{\sin(\theta/2)}\).
For a \(d\) level (\(l=2\)) in octahedral symmetry the five orbitals split into a triplet \(t_{2g}\) and a doublet \(e_g\). This is the starting point of ligand-field theory in chemistry. Bethe’s paper also treated half-integer \(j\) through double groups.
Stopping power (1930–1932)
A heavy particle of charge \(ze\) and speed \(v\) crosses a material with \(n\) electrons per unit volume. Classically, Niels Bohr had computed the momentum given to an electron at impact parameter \(b\): \(\Delta p = 2ze^2/(bv)\) in Gaussian units. The energy lost per unit path is then
\[ -\frac{dE}{dx} = \int_{b_{\min}}^{b_{\max}} \frac{(\Delta p)^2}{2m}\, n\, 2\pi b\, db = \frac{4\pi n z^2 e^4}{m v^2}\,\ln\frac{b_{\max}}{b_{\min}}. \]The logarithm depends on cutoffs. Bethe replaced the classical collision by quantum transitions in the first Born approximation, summed over all final states of the atom using a sum rule, and found that the cutoffs combine into a single property of the material.
For a particle much heavier than the electron, fast compared with atomic electrons but without radiative losses, the mean energy loss per unit length is
\[ -\frac{dE}{dx} = \frac{4\pi n z^2}{m c^2\beta^2}\Big(\frac{e^2}{4\pi\varepsilon_0}\Big)^2\left[\ln\frac{2mc^2\beta^2}{I\,(1-\beta^2)}-\beta^2\right], \]where \(\beta=v/c\), \(m\) is the electron mass and \(I\) is the mean excitation energy, defined by \(\ln I=\sum_k f_k\ln E_k\big/\sum_k f_k\) over the oscillator strengths \(f_k\) of the atom. The nonrelativistic limit is \(\frac{4\pi n z^2e^4}{mv^2}\ln\frac{2mv^2}{I}\) (Gaussian units).
The formula depends on the projectile only through \(z\) and \(v\), and on the material only through \(n\) and \(I\). Empirically \(I\approx 10\,\text{eV}\times Z\) for heavier elements. Later refinements add shell corrections at low speed and the density-effect correction at very high energy. Because \(-dE/dx\propto1/v^2\), the loss rises sharply as the particle slows, producing the Bragg peak near the end of the track. Proton therapy uses that peak to deposit dose in a tumour while sparing tissue beyond it.
Exercise 4.1Using the character formula, show that an \(l=2\) level splits into a 3 + 2 pattern under the rotation group \(O\) of the cube, whose classes are the identity (1 element), rotations by \(\pi\) about face axes (3), by \(\pm\pi/2\) about face axes (6), by \(\pi\) about edge axes (6), and by \(\pm2\pi/3\) about body diagonals (8).
Show solution
Exercise 4.2Show from the nonrelativistic Bethe formula that particles of the same speed but different charge and mass have stopping powers in the ratio \(z_1^2:z_2^2\). Use this to compare a 4 MeV proton with a 16 MeV alpha particle.
Show solution
Exercise 4.3For a fixed material, show that the range \(R\) of a nonrelativistic particle of mass \(M\), charge \(z\), and initial speed \(v_0\) scales as \(R\propto \dfrac{M}{z^2}\,F(v_0)\) for some function \(F\), ignoring the slow variation of the logarithm. How does \(R\) depend on initial energy?
Show solution
The Bethe ansatz
For the general reader. A chain of quantum spins, each coupled to its neighbours, is one of the simplest models of magnetism. In 1931 Bethe guessed the exact form of every energy eigenstate: a superposition of plane waves, with the waves scattering off each other only in pairs. The guess, now called the Bethe ansatz, was the first exact solution of an interacting quantum many-body system. Its descendants run through statistical mechanics, cold-atom physics and string theory.
- Find the one-magnon spectrum of the Heisenberg chain.
- State the Bethe equations and the energy formula.
- Understand why two-body scattering suffices.
The model
Take \(N\) spin-\(\tfrac12\) sites on a ring and the Heisenberg Hamiltonian
\[ H = J\sum_{n=1}^{N}\Big(\mathbf S_n\cdot\mathbf S_{n+1}-\tfrac14\Big) = \frac J2\sum_{n=1}^N\big(P_{n,n+1}-1\big),\qquad \mathbf S_{N+1}\equiv\mathbf S_1, \]where \(P_{n,n+1}\) exchanges the two spins. The state with all spins up has energy 0. Since \(H\) commutes with \(S^z_{\text{tot}}\), we work in sectors with \(M\) flipped spins at positions \(n_1<\dots<n_M\).
The ansatz
For \(M=1\), \(H|n\rangle=\frac J2\big(|n+1\rangle+|n-1\rangle-2|n\rangle\big)\), so the plane wave \(\sum_n e^{ikn}|n\rangle\) has energy \(J(\cos k-1)\), with \(e^{ikN}=1\). Bethe proposed that for general \(M\) the amplitudes are
\[ \psi(n_1,\dots,n_M)=\sum_{P\in S_M} A_P\,\exp\Big(i\sum_{j=1}^M k_{P(j)}\,n_j\Big), \]with amplitudes for permutations differing by a transposition related by a two-body phase. Parametrize momenta by rapidities, \(e^{ik}=\dfrac{\lambda+i/2}{\lambda-i/2}\).
The eigenstates of \(H\) with \(M\le N/2\) flipped spins that are highest-weight for total spin are given by the ansatz, with rapidities \(\lambda_1,\dots,\lambda_M\) solving the Bethe equations
\[ \left(\frac{\lambda_j+i/2}{\lambda_j-i/2}\right)^{\!N}=\prod_{l\ne j}\frac{\lambda_j-\lambda_l+i}{\lambda_j-\lambda_l-i},\qquad j=1,\dots,M, \]and energy \(E=-J\sum_{j=1}^M\dfrac{1/2}{\lambda_j^2+1/4}=J\sum_j(\cos k_j-1)\). Bethe also showed how complex solutions (“strings”) describe bound states and argued that the count of solutions gives all \(2^N\) states.
Why it works. When two flipped spins are not adjacent, \(H\) acts on each separately and any superposition of plane waves solves the equation. Adjacency imposes one extra condition, which fixes the ratio \(A_{\dots k_l k_j\dots}/A_{\dots k_jk_l\dots}\) to the two-body S-matrix \(S(\lambda_j-\lambda_l)=\frac{\lambda_j-\lambda_l-i}{\lambda_j-\lambda_l+i}\) (up to convention). With three or more flipped spins, consistency requires the S-matrices to compose in either order. This is the Yang–Baxter equation, made explicit by Chen Ning Yang and Rodney Baxter decades later. Periodicity then gives the Bethe equations.
For the antiferromagnet (\(J>0\)) the ground state has \(M=N/2\) real rapidities. In 1938 Lamek Hulthén solved the equations in the limit \(N\to\infty\) and found the ground-state energy per site \(E_0/N=-J\ln2\), equivalently \(J(\tfrac14-\ln2)\) for \(J\sum\mathbf S_n\cdot\mathbf S_{n+1}\).
Exercise 5.1Verify that \(e^{ik}=\frac{\lambda+i/2}{\lambda-i/2}\) is solved by \(\lambda=\tfrac12\cot(k/2)\), and that \(J(\cos k-1)=-J\,\dfrac{1/2}{\lambda^2+1/4}\).
Show solution
Exercise 5.2For \(M=1\), what do the Bethe equations say? Interpret the solution with \(k=0\).
Show solution
Exercise 5.3A classical Néel state alternates up and down spins. Compute its energy per bond for \(J\sum\mathbf S_n\cdot\mathbf S_{n+1}\) and compare with Hulthén’s exact value.
Show solution
The nucleus: deuteron, mass formula, the Bible
For the general reader. After the neutron was discovered in 1932, physicists knew what nuclei were made of but not what held them together. Bethe’s approach was to extract as much as possible from the simplest nucleus, the deuteron, and from smooth trends across all nuclei. His three review articles of 1936–37 organized the whole field so thoroughly that it was taught from them for a decade.
- Solve the zero-range deuteron and estimate its size.
- State the effective-range expansion.
- Use the semi-empirical mass formula to find the valley of stability.
The deuteron (Bethe and Peierls, 1935)
The nuclear force has a range of about \(1\) fm, while the deuteron binding energy \(B=2.2246\) MeV is small. Outside the force range the radial wavefunction obeys the free equation with reduced mass \(\mu\approx M/2\):
\[ u(r)\propto e^{-\kappa r},\qquad \kappa=\frac{\sqrt{2\mu B}}{\hbar}=\frac{\sqrt{MB}}{\hbar}. \]Bethe and Peierls realized that most low-energy properties depend only on \(\kappa\), not on the shape of the potential. They computed the photodisintegration cross-section and neutron–proton scattering in this “zero-range” approximation.
For s-wave scattering by a short-range potential, the phase shift \(\delta_0(k)\) satisfies
\[ k\cot\delta_0(k)=-\frac1a+\frac12 r_0k^2+O(k^4), \]where the scattering length \(a\) and effective range \(r_0\) are independent of the details of the potential to this order. A shallow bound state corresponds to a pole at \(k=i\kappa\), so \(\kappa=\frac1a+\frac12r_0\kappa^2+\dots\).
Bethe’s 1949 paper gave a short derivation from the radial equation and became the standard treatment of low-energy two-nucleon scattering. The same expansion now organizes cold-atom physics near Feshbach resonances.
The semi-empirical mass formula
Carl Friedrich von Weizsäcker proposed in 1935 that nuclear binding energies follow from a liquid-drop picture plus quantum corrections. Bethe and Bacher simplified and refitted the formula in the 1936 Bible, and it is often called the Bethe–Weizsäcker formula:
\[ B(A,Z)=a_VA-a_SA^{2/3}-a_C\frac{Z(Z-1)}{A^{1/3}}-a_A\frac{(A-2Z)^2}{A}+\delta(A,Z). \]Typical modern fits give \(a_V\approx15.8\), \(a_S\approx18.3\), \(a_C\approx0.714\), \(a_A\approx23.2\) MeV, and a pairing term \(\delta\approx+12/\sqrt A\) MeV for even–even nuclei, \(-12/\sqrt A\) for odd–odd, and 0 for odd \(A\).
Exercise 6.1Compute \(1/\kappa\) for the deuteron, with \(Mc^2\approx939\) MeV and \(\hbar c=197.3\) MeV fm. Compare with the range of the nuclear force.
Show solution
Exercise 6.2Treating \(Z(Z-1)\approx Z^2\), minimize the mass at fixed \(A\) using only the Coulomb and asymmetry terms (ignore the neutron–proton mass difference). Show that the most stable \(Z\) satisfies \(Z\approx\dfrac{A}{2+(a_C/2a_A)A^{2/3}}\) and evaluate for \(A=56\).
Show solution
Exercise 6.3Estimate the binding energy of \({}^{56}\)Fe with the coefficients above. The measured value is 492.3 MeV.
Show solution
Energy production in stars
For the general reader. Arthur Eddington argued in 1920 that the Sun shines by turning hydrogen into helium, but no one knew how. In 1938 Bethe found the two answers: a chain of reactions starting with two protons fusing, and a cycle in which carbon acts as a catalyst. The Nobel Prize of 1967 was given for this work and for his contributions to the theory of nuclear reactions.
- Write out the proton–proton chain and the CN cycle.
- Derive the Gamow peak and the temperature sensitivity of a thermonuclear rate.
- See why the CN cycle dominates in stars more massive than the Sun.
The two routes
Both routes convert four protons into one helium nucleus, releasing \(26.7\) MeV, part of it carried away by neutrinos.
The first step is a weak interaction and is extremely slow: a given proton in the solar core waits billions of years. Other branches through \({}^7\)Be and \({}^8\)B were added later.
Bethe reached the cycle by elimination. He checked systematically which light nuclei could react with protons fast enough at stellar temperatures without being destroyed, and found that only a closed loop through carbon and nitrogen survived. He also argued, correctly, that no stable nucleus of mass 5 or 8 exists, so helium could not be built up further in main-sequence stars by adding protons or alpha particles.
The Gamow peak
Nuclei in a star at temperature \(T\) have Maxwell–Boltzmann energies, and must tunnel through their Coulomb barrier with probability \(e^{-\sqrt{E_G/E}}\), the Gamow factor, where \(E_G=2\mu c^2(\pi\alpha Z_1Z_2)^2\). Writing the cross-section as \(\sigma(E)=\frac{S(E)}{E}e^{-\sqrt{E_G/E}}\), with \(S(E)\) slowly varying, the rate per pair is
\[ \langle\sigma v\rangle=\sqrt{\frac{8}{\pi\mu}}\,\frac{1}{(kT)^{3/2}}\int_0^\infty S(E)\,\exp\!\Big(-\frac{E}{kT}-\sqrt{\frac{E_G}{E}}\Big)\,dE. \]The integrand peaks at \(E_0=\big(\tfrac14E_G(kT)^2\big)^{1/3}\). In a Gaussian (saddle-point) approximation,
\[ \langle\sigma v\rangle\propto T^{-2/3}\exp\!\Big(-\frac{3E_0}{kT}\Big),\qquad \frac{3E_0}{kT}=3\Big(\frac{E_G}{4kT}\Big)^{1/3}\equiv\tau, \]and near a given \(T\) the rate behaves like \(T^{\nu}\) with \(\nu=\dfrac{\tau-2}{3}\).
Because \(E_G\) grows as \(Z_1^2Z_2^2\), reactions with nitrogen (\(Z=7\)) are far more temperature-sensitive than proton–proton fusion. At low temperatures the pp chain wins; above about 17 million kelvin the CN cycle does. Bethe concluded that the Sun runs mainly on the pp chain and more massive stars on the CN cycle, which is the modern picture.
Exercise 7.1At the solar centre, \(T\approx15.7\times10^6\) K, so \(kT\approx1.35\) keV. For p + p, \(E_G\approx493\) keV. Find \(E_0\) and \(\nu\). Repeat for p + \({}^{14}\)N, with \(E_G\approx45\) MeV.
Show solution
Exercise 7.2The solar luminosity is \(3.83\times10^{26}\) W. Each conversion of four protons to helium releases about 26.7 MeV and two neutrinos. Estimate the solar neutrino flux at Earth (distance \(1.50\times10^{11}\) m).
Show solution
Exercise 7.3Show that the CN cycle changes the relative abundances of C, N and O isotopes even though it conserves their total. Which isotope accumulates, and why?
Show solution
The Lamb shift and renormalization
For the general reader. Paul Dirac’s theory predicted that two excited levels of hydrogen have exactly the same energy. In 1947 Willis Lamb and Robert Retherford measured a small difference. Days later Bethe explained it, with a calculation he did largely on a train. The interaction of the electron with the quantum vacuum gives infinite answers, but the difference between a bound electron and a free one is finite. This idea, mass renormalization, made quantum electrodynamics a working theory.
- Set up the self-energy of a bound electron in second-order perturbation theory.
- See how subtracting the free-electron self-energy removes the linear divergence.
- Evaluate the 2s–2p splitting.
The calculation
At the Shelter Island conference (2–4 June 1947), Lamb reported that the \(2s_{1/2}\) level of hydrogen lies about 1000 MHz above \(2p_{1/2}\). Hendrik Kramers had suggested that the electromagnetic mass of the electron should be absorbed into the observed mass. Bethe took the idea literally. In nonrelativistic second-order perturbation theory, the emission and reabsorption of a transverse photon shifts level \(n\) by
\[ W_n=-\frac{2\alpha}{3\pi}\frac{1}{m^2c^2}\int_0^{K}k\,dk\sum_m\frac{|\langle m|\mathbf p|n\rangle|^2}{E_m-E_n+k}, \]with photon energy \(k\) and a cutoff \(K\). This diverges linearly. A free electron has a similar term, which only renormalizes its mass. Subtracting it, Bethe obtained an expression that diverges only logarithmically:
\[ W_n-W_n^{\text{free}}=\frac{2\alpha}{3\pi}\frac{1}{m^2c^2}\sum_m|\langle m|\mathbf p|n\rangle|^2(E_m-E_n)\ln\frac{K}{|E_m-E_n|}. \]He argued that relativity would cut off the integral at \(K\approx mc^2\). For the sum he used a closure relation, \(\sum_m|\langle m|\mathbf p|n\rangle|^2(E_m-E_n)=\frac{\hbar^2}{2}\langle n|\nabla^2V|n\rangle=2\pi\hbar^2\alpha\hbar c\,Z|\psi_n(0)|^2\), which vanishes unless the state has \(l=0\), and he replaced the energy in the logarithm by a numerically computed average.
For an s state of hydrogen-like ions, to leading order,
\[ \Delta E_{n,l=0}=\frac{4\alpha^5Z^4}{3\pi n^3}\,mc^2\,\ln\frac{mc^2}{\langle E_m-E_n\rangle_{\text{av}}}, \]with \(\langle E_m-E_n\rangle_{\text{av}}\approx17.8\) Ry for the 2s state. States with \(l>0\) are shifted much less. The result for the 2s level is about 1040 MHz.
Bethe’s three-page paper, The electromagnetic shift of energy levels, appeared in Physical Review in August 1947. It used no relativistic field theory at all. Within two years Julian Schwinger, Richard Feynman, Sin-Itiro Tomonaga and Freeman Dyson, who had come to Cornell as Bethe’s student, gave the fully relativistic, covariant renormalization of QED. The modern value, including vacuum polarization and higher orders, is 1057.8 MHz.
Exercise 8.1Evaluate Bethe’s formula for the 2s state of hydrogen, with \(\alpha=1/137.036\), \(mc^2=511.0\) keV and 1 Ry \(=13.606\) eV. Convert to MHz using \(h=4.1357\times10^{-15}\) eV s.
Show solution
Exercise 8.2Show that \(\sum_m|\langle m|p_x|n\rangle|^2(E_m-E_n)=\tfrac12\langle n|[p_x,[H,p_x]]|n\rangle\), and evaluate the double commutator for \(H=p^2/2m+V\).
Show solution
Exercise 8.3Ignoring the change in the logarithm, how does the Lamb shift of the 1s level compare with that of the 2s level? Why is the shift of the 2p level small?
Show solution
Bound states and nuclear matter
For the general reader. After the Lamb shift, Bethe turned to the hardest problems of bound systems: two relativistic particles held together by exchanging quanta, and the dense sea of nucleons inside a heavy nucleus. In both cases he and his collaborators found equations that sum infinite classes of interactions, and that are still the starting points of their fields.
- State the Bethe–Salpeter equation and its ladder approximation.
- Understand the Bethe–Goldstone equation and Pauli blocking.
- Compute basic properties of nuclear matter.
The Bethe–Salpeter equation (1951)
Edwin Salpeter and Bethe wrote a covariant equation for the two-body bound state in field theory. Let \(G\) be the full two-particle Green’s function and \(K\) the two-particle-irreducible interaction kernel. Then \(G=G_0+G_0KG\), where \(G_0\) is the product of the two dressed one-particle propagators.
A bound state of total four-momentum \(P\) appears as a pole of \(G\). Its amplitude \(\chi_P(p)\) satisfies the homogeneous equation
\[ \chi_P(p)=S_1\big(\tfrac P2+p\big)\,S_2\big(\tfrac P2-p\big)\int\frac{d^4p'}{(2\pi)^4}\,K(p,p';P)\,\chi_P(p'). \]Keeping only single-quantum exchange in \(K\) gives the ladder approximation, which sums all ladder diagrams.
The equation includes relative time (or energy) as a variable, which makes it hard to solve and to interpret. It became the basis for the precise theory of positronium and for relativistic models of mesons. Murray Gell-Mann and Francis Low derived the same equation independently.
Nuclear matter and the Bethe–Goldstone equation
Nuclear forces have a strongly repulsive core at short distance, so ordinary perturbation theory fails. Keith Brueckner proposed in the mid-1950s to replace the bare potential \(V\) by an effective interaction summing repeated two-body scatterings inside the medium. Bethe and his student Jeffrey Goldstone (1957) solved the resulting two-body problem.
where \(Q\) projects onto two-particle states with both momenta above the Fermi momentum \(k_F\). Because the Pauli principle forbids scattering into occupied states, the correlated pair wavefunction “heals” to the uncorrelated one within about 1 fm, even for an infinite hard core.
The healing property explains why the independent-particle shell model works despite strong forces. Over the following fifteen years, with Baird Brandow, Albert Petschek and others, Bethe pushed this theory to a quantitative account of nuclear saturation, summarized in his 1971 review Theory of nuclear matter. The three-body extension led to the Bethe–Faddeev equations. The same machinery later fed his work on the equation of state of neutron stars.
Exercise 9.1Nuclear matter has density \(\rho_0\approx0.16\) fm\(^{-3}\) with four nucleon species (spin and isospin). Find \(k_F\) and the Fermi kinetic energy, with \(\hbar c=197.3\) MeV fm and \(Mc^2=939\) MeV.
Show solution
Exercise 9.2Iterate \(G=V+VG_0G\) to show that \(G=V+VG_0V+VG_0VG_0V+\cdots\). Draw the corresponding diagrams and explain why the series makes sense for a hard core even though each term is infinite.
Show solution
Supernovae and solar neutrinos
For the general reader. In his seventies Bethe started a third career in astrophysics. With Gerald Brown he worked on why the collapsing core of a massive star bounces and explodes as a supernova. With James Wilson he identified the mechanism still favoured today: neutrinos from the hot new neutron star reheat the stalled shock. He also gave a clear explanation of why detectors saw only a third of the expected solar neutrinos: the neutrinos change type on their way out of the Sun.
- Estimate the energy released in core collapse.
- Outline the delayed neutrino-heating mechanism.
- Derive the MSW resonance condition.
Core collapse
The iron core of a massive star, about 1.4 solar masses, collapses when electron captures and photodisintegration remove its pressure support. In 1979 Bethe, Brown, Jim Applegate and James Lattimer showed that the collapsing matter has low entropy, so nuclei survive until nuclear density, and the core bounces stiffly. The bounce launches a shock. Detailed simulations showed that the shock loses energy dissociating iron and stalls a few hundred kilometres out.
Behind the stalled shock, matter absorbs a small fraction of the neutrinos streaming from the proto-neutron star, mainly through \(\nu_e+n\to p+e^-\) and \(\bar\nu_e+p\to n+e^+\). If heating exceeds cooling in a “gain region” over a time of order several hundred milliseconds, the shock is revived and the star explodes.
The neutron star’s binding energy, about \(3\times10^{53}\) erg, is radiated almost entirely as neutrinos of all flavours over roughly ten seconds. About 1% reaches the ejecta as kinetic energy. The detection of about two dozen neutrinos from SN 1987A confirmed the energy scale and the duration.
The MSW effect (Bethe, 1986)
Lincoln Wolfenstein (1978) showed that electron neutrinos moving through matter acquire an extra potential \(V=\sqrt2G_Fn_e\) from forward scattering on electrons. Stanislav Mikheyev and Alexei Smirnov (1985) found that this can cause resonant conversion. Bethe’s short 1986 paper gave the adiabatic picture in a form that made it widely understood. For two flavours with mixing angle \(\theta\) and mass-squared difference \(\Delta m^2\), the effective mixing angle in matter satisfies
\[ \tan2\theta_m=\frac{\sin2\theta}{\cos2\theta-\dfrac{2\sqrt2\,G_Fn_eE}{\Delta m^2}}. \]Maximal mixing occurs where \(\sqrt2\,G_Fn_e=\dfrac{\Delta m^2\cos2\theta}{2E}\). If the density decreases slowly enough along the neutrino path (adiabatic condition), an electron neutrino created above resonance density emerges as the heavier mass eigenstate, mostly of another flavour.
Bethe used the deficit to estimate \(\Delta m^2\) and argued that solar neutrino data required physics beyond the Standard Model. The Sudbury Neutrino Observatory confirmed flavour conversion in 2001–02; the large-mixing-angle solution is the one realized in nature.
Exercise 10.1Estimate the gravitational binding energy \(\tfrac35GM^2/R\) of a neutron star with \(M=1.4M_\odot\) and \(R=12\) km. Compare with the kinetic energy of a typical supernova, \(10^{51}\) erg.
Show solution
Exercise 10.2For \(\Delta m^2=7.5\times10^{-5}\) eV\(^2\), \(\sin^2\theta=0.31\) and \(E=10\) MeV, find the resonance electron density and the corresponding mass density for \(Y_e\approx0.7\) (electrons per nucleon). Use \(G_F=1.166\times10^{-23}\) eV\(^{-2}\) and \(\hbar c=1.973\times10^{-5}\) eV cm. Is the resonance inside the Sun, whose central density is about 150 g cm\(^{-3}\)?
Show solution
Los Alamos, the H-bomb and arms control, 1941–1995
Bethe married Rose Ewald, Paul Ewald’s daughter, on 13 September 1939; they had two children, Henry and Monica. When war came he looked for a way to help his new country. He became an American citizen in 1941 and joined the MIT Radiation Laboratory, where he worked on microwave radar and invented a waveguide directional coupler still called the Bethe-hole coupler. In the summer of 1942 he joined Oppenheimer’s study group in Berkeley on the design of a fission bomb. When Edward Teller raised the possibility that a nuclear explosion could ignite the atmosphere, Bethe’s calculations helped show that it could not.
In 1943 Oppenheimer made Bethe head of the Theoretical Division at Los Alamos. The appointment hurt Teller, who had hoped for the post, and the rift between them never fully healed. Bethe’s division calculated critical masses, the efficiency of the explosion and the hydrodynamics of implosion. With Richard Feynman he derived a formula for the yield of a fission weapon. He witnessed the Trinity test in July 1945.
After the war Bethe helped found the Federation of American Scientists. In 1950 he argued in Scientific American against a crash programme to build a hydrogen bomb. After the Korean War began and Truman ordered the programme, and after Teller and Stanislaw Ulam found a workable design in 1951, he went back to Los Alamos to help, saying later that he had hoped to show it could not be done. In 1954 he testified for Oppenheimer at the security hearing that removed Oppenheimer’s clearance.
From 1957 he served on President Eisenhower’s Science Advisory Committee. He chaired an interagency panel on whether a nuclear test ban could be verified and took part in the 1958 Geneva conference of experts on detection. These efforts fed into the Partial Test Ban Treaty of 1963. In 1968 he and Richard Garwin argued in Scientific American that anti-ballistic-missile defences could be defeated cheaply, a case that helped shape the ABM Treaty of 1972. In the 1980s he opposed the Strategic Defense Initiative. In 1995, on the fiftieth anniversary of Hiroshima, he called on scientists everywhere to stop working on new nuclear weapons. He remained a strong supporter of civil nuclear power.
DiscussionBethe worked on the hydrogen bomb while opposing its construction, and argued that his presence made the programme and later arms control more sensible. Assess this position. What responsibilities does a scientist have when the state asks for a weapon?
Cornell, honours and legacy
Bethe was Cornell’s theoretical physics programme for decades. He brought Feynman there after the war, and he attracted students who became leaders in many fields. His doctoral students included Robert Marshak, Francis Low, Jun John Sakurai, Jeffrey Goldstone, David Thouless, Roman Jackiw, Michel Baranger and Peter Carruthers. Freeman Dyson came to work with him in 1947 without ever taking a doctorate. Bethe was famous for his ability to estimate anything quickly and for writing carefully worked lecture notes. He continued to publish into his nineties, much of it with Gerald Brown on supernovae, neutron stars and black-hole binaries.
The Nobel Prize in Physics for 1967 cited his contributions to the theory of nuclear reactions, especially concerning energy production in stars. Other honours included the Henry Draper Medal (1947), the Max Planck Medal (1955), Foreign Membership of the Royal Society (1957), the Eddington Medal (1961), the Enrico Fermi Award (1961), the Rumford Prize (1963), the National Medal of Science (1975), the Lomonosov Gold Medal (1989), the Oersted Medal for teaching (1993) and the Bruce Medal (2001).
Bethe died of congestive heart failure at his home in Ithaca on 6 March 2005, aged 98. Rose Bethe survived him and lived until 2019.
Legacy
- Methods: the Bethe ansatz and its descendants, the Bethe lattice and Bethe–Peierls approximation in statistical mechanics, effective-range theory, and the Bethe–Salpeter equation.
- Applications: the Bethe formula underlies particle detectors, shielding and proton therapy; crystal field theory underlies ligand-field chemistry.
- Astrophysics: the theory of stellar energy and the delayed neutrino mechanism for supernovae.
- Public role: five decades of technical advice on nuclear weapons, test bans and missile defence.
Final projectChoose one of Bethe’s papers from Part II. Reproduce its central estimate with modern data, identify which approximation in the original carries the largest error, and explain why the answer was nevertheless close to the truth.
Sources and further reading
Biography
- Hans Bethe, Wikipedia (with extensive references).
- Hans Bethe, Biographical, NobelPrize.org.
- J. J. O’Connor and E. F. Robertson, Hans Albrecht Bethe, MacTutor History of Mathematics.
- In memoriam: Hans Bethe, 1906–2005, Arms Control Today, April 2005.
- S. S. Schweber, Nuclear Forces: The Making of the Physicist Hans Bethe, Harvard University Press, 2012.
- J. Bernstein, Hans Bethe, Prophet of Energy, Basic Books, 1980.
- G. E. Brown and S. Lee, Hans Albrecht Bethe 1906–2005, Biographical Memoirs of the National Academy of Sciences, 2009.
Physics
- H. Bethe, Theorie der Beugung von Elektronen an Kristallen, Ann. Phys. 87 (1928).
- H. Bethe, Termaufspaltung in Kristallen, Ann. Phys. 3 (1929).
- H. Bethe, Zur Theorie des Durchgangs schneller Korpuskularstrahlen durch Materie, Ann. Phys. 5 (1930).
- H. Bethe, Zur Theorie der Metalle. I. Eigenwerte und Eigenfunktionen der linearen Atomkette, Z. Phys. 71 (1931).
- H. Bethe and R. Peierls, Quantum theory of the diplon, Proc. Roy. Soc. A 148 (1935).
- H. Bethe, R. F. Bacher and M. S. Livingston, Nuclear physics, Rev. Mod. Phys. 8 (1936) and 9 (1937).
- H. Bethe and C. L. Critchfield, The formation of deuterons by proton combination, Phys. Rev. 54 (1938).
- H. Bethe, Energy production in stars, Phys. Rev. 55 (1939).
- H. Bethe, The electromagnetic shift of energy levels, Phys. Rev. 72 (1947).
- H. Bethe, Theory of the effective range in nuclear scattering, Phys. Rev. 76 (1949).
- E. E. Salpeter and H. A. Bethe, A relativistic equation for bound-state problems, Phys. Rev. 84 (1951).
- H. A. Bethe and J. Goldstone, Effect of a repulsive core in the theory of complex nuclei, Proc. Roy. Soc. A 238 (1957).
- H. A. Bethe and J. R. Wilson, Revival of a stalled supernova shock by neutrino heating, Astrophys. J. 295 (1985).
- H. A. Bethe, Possible explanation of the solar-neutrino puzzle, Phys. Rev. Lett. 56 (1986).
- H. A. Bethe, Supernova mechanisms, Rev. Mod. Phys. 62 (1990).
- H. A. Bethe, Nobel Lecture: Energy production in stars, 1967, NobelPrize.org.
- M. Karbach and G. Müller, Introduction to the Bethe ansatz I, arXiv:cond-mat/9809162.