courseshub.world · Lives in Mathematics and Physics

Hans Bethe

Hans Albrecht Bethe

A course in two strands: the physics of the calculator who explained why the stars shine, tamed the infinities of quantum electrodynamics and wrote nuclear physics its first textbook, and the life of a refugee who built the atomic bomb and then spent fifty years trying to limit it.

Born 2 July 1906, Strasbourg Died 6 March 2005, Ithaca, New York Nobel Prize Physics, 1967 Level advanced undergraduate → graduate
The CN cycle of 1938–39. Carbon-12 is a catalyst: four protons go in, one helium nucleus comes out, and the carbon is returned.
Lesson 00

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.

Prerequisites for Part II
  • 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).

Reference

Timeline

1906Born 2 July in Strasbourg, then in the German Empire, only child of the physiologist Albrecht Bethe and Anna Kuhn.
1912–15The family moves to Kiel, then to Frankfurt.
1924Enters the University of Frankfurt.
1926Moves to Munich to study with Arnold Sommerfeld.
1928July: doctorate on the diffraction of electrons by crystals. Short posts at Frankfurt and Stuttgart (with Paul Ewald).
1929Crystal field theory: the splitting of atomic levels in crystals.
1930Theory of the passage of fast particles through matter. Privatdozent at Munich. Rockefeller fellowship: Cambridge.
1931Rome with Enrico Fermi; exact solution of the one-dimensional Heisenberg chain.
1932–33Acting assistant professor at Tübingen. Dismissed in 1933 under the Nazi civil-service law.
1933–34Lecturer at Manchester; with Rudolf Peierls explains the photodisintegration of the deuteron. Fellow at Bristol in autumn 1934.
1935February: assistant professor at Cornell University. Professor from 1937.
1936–37The three Reviews of Modern Physics articles known as “Bethe’s Bible”.
1938Washington conference on stellar energy; proton–proton reaction with Charles Critchfield; the CN cycle.
1939Energy Production in Stars, Physical Review. 13 September: marries Rose Ewald.
1941–42American citizen. Radar work at the MIT Radiation Laboratory.
1943–46Head of the Theoretical Division at Los Alamos.
1947Shelter Island conference; nonrelativistic calculation of the Lamb shift.
1948Name added by George Gamow to the Alpher–Bethe–Gamow paper.
1951Bethe–Salpeter equation. Returns to thermonuclear weapons work after the Teller–Ulam design.
1954Testifies for J. Robert Oppenheimer at his security hearing.
1957–58President’s Science Advisory Committee; chairs a panel on verifying a test ban; Geneva conference of experts.
1961Enrico Fermi Award.
1967Nobel Prize in Physics.
1968With Richard Garwin, Scientific American article against anti-ballistic-missile systems.
1979Begins collaboration with Gerald Brown on stellar collapse.
1985–86Delayed neutrino-heating supernova mechanism with James Wilson; MSW explanation of the solar-neutrino deficit.
1995Calls on scientists to stop work on new nuclear weapons.
2005Dies 6 March at home in Ithaca, aged 98.
Part I · Formation
Lesson 01

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.

Lesson 02

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.

Lesson 03

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).

Part II · The physics
Lesson 04

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.

Goals
  • 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.

Principle (Bethe, 1929)

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.

Bethe formula (1930; relativistic form 1932)

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
\(\chi^{(2)}(\theta)=\sin(5\theta/2)/\sin(\theta/2)\) gives \(5,\ 1,\ -1,\ 1,\ -1\) on the five classes. The sum of squares weighted by class size is \(25+3+6+6+8=48=2|O|\), so \(D^{(2)}\) contains two distinct irreducibles. Taking inner products with the character table of \(O\cong S_4\) gives \(E\) (dimension 2, character \(2,2,0,0,-1\)) and \(T_2\) (dimension 3, character \(3,-1,-1,1,0\)): \(\tfrac1{24}(10+6+0+0+8)=1\) and \(\tfrac1{24}(15-3+6+6+0)=1\). So the level splits into a doublet and a triplet.

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
The formula depends only on \(z\) and \(v\). Kinetic energy \(\tfrac12Mv^2\) at equal speed scales with mass; the alpha particle has about four times the proton mass, so 16 MeV alpha and 4 MeV proton have nearly the same speed. With \(z=2\) versus \(z=1\), the alpha loses energy about four times faster per unit length.

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
\(dE=Mv\,dv\) and \(-dE/dx\approx C z^2/v^2\). So \(dx=-\frac{Mv^3}{Cz^2}dv\) and \(R=\frac{M v_0^4}{4Cz^2}\). In terms of \(E_0=\tfrac12Mv_0^2\), \(R=\frac{E_0^2}{CMz^2}\): range grows as the square of the initial energy. Including the logarithm gives the empirical power closer to 1.7–1.8.
Lesson 05

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.

Goals
  • 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}\).

Theorem (Bethe, 1931)

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
With \(\lambda=\tfrac12\cot(k/2)\), \(\frac{\lambda+i/2}{\lambda-i/2}=\frac{\cos(k/2)+i\sin(k/2)}{\cos(k/2)-i\sin(k/2)}=e^{ik}\). Then \(\lambda^2+\tfrac14=\frac{1}{4\sin^2(k/2)}\), so \(\frac{1/2}{\lambda^2+1/4}=2\sin^2(k/2)=1-\cos k\).

Exercise 5.2For \(M=1\), what do the Bethe equations say? Interpret the solution with \(k=0\).

Show solution
The product is empty, so \(e^{ikN}=1\) and \(k=2\pi m/N\). The solution \(k=0\) (\(\lambda=\infty\)) has energy 0: it is the uniform superposition, obtained from the all-up state by the total lowering operator. It belongs to the same spin multiplet, which is why highest-weight states use finite rapidities.

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
For a Néel state \(\langle\mathbf S_n\cdot\mathbf S_{n+1}\rangle=\langle S^z_nS^z_{n+1}\rangle=-\tfrac14\), since the transverse terms have zero expectation. The exact ground state has \(\tfrac14-\ln2\approx-0.443\) per bond, much lower. Quantum fluctuations destroy long-range order in one dimension and lower the energy.
Lesson 06

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.

Goals
  • 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.

Effective-range expansion (Bethe, 1949)

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
\(\hbar c\,\kappa=\sqrt{Mc^2B}=\sqrt{939\times2.2246}\approx45.7\) MeV, so \(1/\kappa\approx197.3/45.7\approx4.3\) fm. The wavefunction extends well beyond the 1–2 fm range of the force, which is why the zero-range approximation works.

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
Maximize \(-a_CZ^2A^{-1/3}-a_A(A-2Z)^2/A\): \(-2a_CZA^{-1/3}+4a_A(A-2Z)/A=0\), so \(Z(8a_A/A+2a_CA^{-1/3})=4a_A\), giving \(Z=\frac{A}{2+(a_C/2a_A)A^{2/3}}\). With \(a_C/2a_A\approx0.0154\) and \(56^{2/3}\approx14.6\), \(Z\approx56/2.23\approx25\). The stable isobar is iron, \(Z=26\); including the neutron–proton mass difference moves the estimate up slightly.

Exercise 6.3Estimate the binding energy of \({}^{56}\)Fe with the coefficients above. The measured value is 492.3 MeV.

Show solution
Volume \(15.8\times56=884.8\); surface \(18.3\times14.64=267.9\); Coulomb \(0.714\times26\times25/3.826=121.3\); asymmetry \(23.2\times16/56=6.6\); pairing \(+12/\sqrt{56}=+1.6\). Total \(\approx490.6\) MeV, about \(8.76\) MeV per nucleon, within 0.4% of experiment.
Lesson 07

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.

Goals
  • 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 pp chain (Bethe and Critchfield, 1938) \[ \begin{aligned} p+p&\to d+e^++\nu_e\\ d+p&\to{}^3\mathrm{He}+\gamma\\ {}^3\mathrm{He}+{}^3\mathrm{He}&\to{}^4\mathrm{He}+2p \end{aligned} \]

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.

The CN cycle (Bethe, 1938–39; independently von Weizsäcker) \[ \begin{aligned} {}^{12}\mathrm C+p&\to{}^{13}\mathrm N+\gamma, & {}^{13}\mathrm N&\to{}^{13}\mathrm C+e^++\nu_e,\\ {}^{13}\mathrm C+p&\to{}^{14}\mathrm N+\gamma, & {}^{14}\mathrm N+p&\to{}^{15}\mathrm O+\gamma,\\ {}^{15}\mathrm O&\to{}^{15}\mathrm N+e^++\nu_e, & {}^{15}\mathrm N+p&\to{}^{12}\mathrm C+{}^4\mathrm{He}. \end{aligned} \]

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. \]
Result (Gamow peak)

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
p + p: \(E_0=(493\times1.35^2/4)^{1/3}\approx6.1\) keV, \(\tau=3\times6.1/1.35\approx13.5\), \(\nu\approx3.8\). p + \({}^{14}\)N: \(E_0=(45000\times1.35^2/4)^{1/3}\approx27\) keV, \(\tau\approx61\), \(\nu\approx20\). Reactions occur at energies several times \(kT\), and the CN rate rises like \(T^{20}\).

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
26.7 MeV \(=4.28\times10^{-12}\) J, so the Sun performs \(\approx8.9\times10^{37}\) conversions per second and emits \(\approx1.8\times10^{38}\) neutrinos per second. Dividing by \(4\pi(1.5\times10^{11})^2\approx2.8\times10^{23}\) m\(^2\) gives \(\approx6.4\times10^{14}\) m\(^{-2}\)s\(^{-1}\), or about \(6\times10^{10}\) per cm\(^2\) per second.

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
In equilibrium every step proceeds at the same rate, so each abundance is inversely proportional to its destruction rate. \({}^{14}\)N + p is the slowest step, so \({}^{14}\)N accumulates: most carbon is converted into nitrogen. This is the origin of much of the nitrogen in the universe.
Lesson 08

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.

Goals
  • 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.

Bethe’s formula (1947)

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
\(\alpha^5mc^2\approx1.06\times10^{-5}\) eV and \(\frac{4}{3\pi\cdot8}\approx0.0531\), giving \(5.6\times10^{-7}\) eV. The logarithm is \(\ln(511000/242)\approx7.66\). Product \(\approx4.29\times10^{-6}\) eV, or about 1040 MHz.

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
Expand \(\langle n|[p_x,[H,p_x]]|n\rangle=2\langle n|p_xHp_x|n\rangle-\langle n|Hp_x^2+p_x^2H|n\rangle=2\sum_m|\langle m|p_x|n\rangle|^2(E_m-E_n)\). With \([H,p_x]=i\hbar\,\partial_xV\), \([p_x,i\hbar\partial_xV]=\hbar^2\partial_x^2V\). Summing over \(x,y,z\) gives \(\tfrac{\hbar^2}{2}\langle\nabla^2V\rangle\), and for the Coulomb potential \(\nabla^2V\) is a delta function at the origin.

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
The shift is proportional to \(|\psi_n(0)|^2\propto1/n^3\), so the 1s shift is about 8 times larger (about 8.2 GHz measured). For \(l>0\), \(\psi(0)=0\) and the leading term vanishes; only a much smaller logarithm-free contribution remains.
Lesson 09

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.

Goals
  • 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.

Bethe–Salpeter equation (Salpeter and Bethe, 1951)

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.

Bethe–Goldstone equation (1957) \[ G(\omega)=V+V\,\frac{Q}{\omega-H_0}\,G(\omega), \]

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
\(\rho=4\cdot\frac{k_F^3}{6\pi^2}=\frac{2k_F^3}{3\pi^2}\), so \(k_F^3=0.16\times3\pi^2/2\approx2.37\) fm\(^{-3}\) and \(k_F\approx1.33\) fm\(^{-1}\). Then \(E_F=(\hbar ck_F)^2/(2Mc^2)=(263)^2/1878\approx37\) MeV. The average kinetic energy is \(\tfrac35E_F\approx22\) MeV per nucleon, against a binding of about 16 MeV per nucleon.

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
Substituting the equation into itself repeatedly gives the series; each term is a ladder with \(n\) rungs \(V\) joined by intermediate two-particle propagators \(G_0\). For a hard core, \(V=\infty\) inside the core, so individual terms diverge, but the sum \(G\) is the scattering matrix of a well-defined problem (the wavefunction vanishes inside the core), and \(G\) itself is finite.
Lesson 10

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.

Goals
  • 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.

Delayed neutrino heating (Bethe and Wilson, 1985)

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}}. \]
Resonance condition

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
\(M=2.79\times10^{30}\) kg, so \(GM^2/R=6.67\times10^{-11}\times7.8\times10^{60}/1.2\times10^4\approx4.3\times10^{46}\) J and \(\tfrac35\) of it is \(2.6\times10^{46}\) J \(=2.6\times10^{53}\) erg. The explosion takes only about 0.4% of this; neutrinos carry the rest.

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
\(\cos2\theta=1-2(0.31)=0.38\). Then \(n_e=\frac{\Delta m^2\cos2\theta}{2\sqrt2G_FE}\) \(\approx2.85\times10^{-5}/3.30\times10^{-16}\approx8.6\times10^{10}\) eV\(^3\) in natural units; dividing by \((\hbar c)^3\) gives \(n_e\approx1.1\times10^{25}\) cm\(^{-3}\), or about 18.7 moles of electrons per cm\(^3\). Mass density \(\rho\approx18.7/Y_e\approx27\) g cm\(^{-3}\). This lies inside the Sun, so high-energy \({}^8\)B neutrinos pass through resonance.
Part III · Life and legacy
Lesson 11

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?

Lesson 12

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.

Reference

Sources and further reading

Biography

  1. Hans Bethe, Wikipedia (with extensive references).
  2. Hans Bethe, Biographical, NobelPrize.org.
  3. J. J. O’Connor and E. F. Robertson, Hans Albrecht Bethe, MacTutor History of Mathematics.
  4. In memoriam: Hans Bethe, 1906–2005, Arms Control Today, April 2005.
  5. S. S. Schweber, Nuclear Forces: The Making of the Physicist Hans Bethe, Harvard University Press, 2012.
  6. J. Bernstein, Hans Bethe, Prophet of Energy, Basic Books, 1980.
  7. G. E. Brown and S. Lee, Hans Albrecht Bethe 1906–2005, Biographical Memoirs of the National Academy of Sciences, 2009.

Physics

  1. H. Bethe, Theorie der Beugung von Elektronen an Kristallen, Ann. Phys. 87 (1928).
  2. H. Bethe, Termaufspaltung in Kristallen, Ann. Phys. 3 (1929).
  3. H. Bethe, Zur Theorie des Durchgangs schneller Korpuskularstrahlen durch Materie, Ann. Phys. 5 (1930).
  4. H. Bethe, Zur Theorie der Metalle. I. Eigenwerte und Eigenfunktionen der linearen Atomkette, Z. Phys. 71 (1931).
  5. H. Bethe and R. Peierls, Quantum theory of the diplon, Proc. Roy. Soc. A 148 (1935).
  6. H. Bethe, R. F. Bacher and M. S. Livingston, Nuclear physics, Rev. Mod. Phys. 8 (1936) and 9 (1937).
  7. H. Bethe and C. L. Critchfield, The formation of deuterons by proton combination, Phys. Rev. 54 (1938).
  8. H. Bethe, Energy production in stars, Phys. Rev. 55 (1939).
  9. H. Bethe, The electromagnetic shift of energy levels, Phys. Rev. 72 (1947).
  10. H. Bethe, Theory of the effective range in nuclear scattering, Phys. Rev. 76 (1949).
  11. E. E. Salpeter and H. A. Bethe, A relativistic equation for bound-state problems, Phys. Rev. 84 (1951).
  12. H. A. Bethe and J. Goldstone, Effect of a repulsive core in the theory of complex nuclei, Proc. Roy. Soc. A 238 (1957).
  13. H. A. Bethe and J. R. Wilson, Revival of a stalled supernova shock by neutrino heating, Astrophys. J. 295 (1985).
  14. H. A. Bethe, Possible explanation of the solar-neutrino puzzle, Phys. Rev. Lett. 56 (1986).
  15. H. A. Bethe, Supernova mechanisms, Rev. Mod. Phys. 62 (1990).
  16. H. A. Bethe, Nobel Lecture: Energy production in stars, 1967, NobelPrize.org.
  17. M. Karbach and G. Müller, Introduction to the Bethe ansatz I, arXiv:cond-mat/9809162.