Nuclear physics 2026-10-06
The study of atomic nuclei, their structure, reactions and interactions. Nuclear forces between nucleons and nuclear binding energies connect it to the low-energy consequences of Quantum chromodynamics.
QCD supplies the underlying strong interaction; Skyrmions provide a mesonic effective description of baryons, and quantized multi-Skyrmions can model nuclei. These are related descriptions at different scales, not three identical theories.
In QCD, quarks carry color charge and interact through gluons, the gauge fields of the color special unitary group . Its Lagrangian density has the form
Asymptotic freedom makes short-distance processes accessible through small-coupling expansions, but nuclear scales involve strongly coupled, confined dynamics. Observable hadrons are color singlets. A nucleon, either a proton or a neutron, is a baryon with baryon number one; a nucleus contains such units of baryon number. Directly extracting all nuclear binding energies, spectra and interactions from QCD is difficult, motivating low-energy effective field theories that preserve its symmetries and relevant degrees of freedom.
For the two light quark flavours, the small-mass limit has approximate chiral symmetry . Chiral symmetry breaking leaves its vector subgroup , the approximate isospin symmetry. The three pions are the associated Goldstone bosons in the massless limit and pseudo-Goldstone bosons when the light quark masses are retained. Package these pions into a special unitary group field
where is the pion decay constant in this normalization and are the Pauli matrices. The nonlinear sigma model is the leading two-derivative mesonic theory. The Skyrme model adds a specific four-derivative stabilizing interaction. One conventional normalization is
Here is a dimensionless model coupling, not electric charge. The last term accounts for a common pion mass and preserves vector isospin. It vanishes in the chiral massless limit. This effective field theory uses color-singlet mesonic fields and does not resolve constituent quarks or gluons inside a baryon.
The condition at spatial infinity compactifies physical space to . Since is itself a three-sphere, the field defines a map with integer topological charge in . This is identified with the topological baryon number in the Skyrme model:
The associated topological current is identically conserved. A single Skyrmion has , and a multi-Skyrmion with is a candidate intrinsic configuration for an ordinary nucleus; negative charge describes antibaryonic sectors. Integer topology prevents a smooth finite-energy unwinding into the classical vacuum, but it does not by itself guarantee a nonzero-size energy minimum.
The energetic reason for the Skyrme term is Derrick scaling. For the rescaled field in three dimensions, let be the quadratic-gradient, quartic-gradient, and potential energies. Their scale dependence is
The two-derivative nonlinear sigma model alone can lower its static energy by shrinking. The positive quartic Skyrme term instead grows under shrinking, permitting a balance and a stable soliton size. Without the mass term, this balance gives . The displayed scaling convention uses , so it is the inverse of the equally common convention.
The connection with QCD is strengthened by large-Nc baryon scaling. Generalize the number of colors to while keeping fixed. Meson masses remain of order one, their interactions weaken, and an effective mesonic action has an overall scale of order . In the Skyrme model this corresponds to and of order , so the soliton mass and rotational moment of inertia are also of order , whereas rotational level splittings are of order . These are the expected baryon scaling properties of large- QCD. A massive, semiclassical soliton built from meson fields is therefore consistent with the underlying theory, even though physical is only a finite value and the simplest Skyrme model is not uniquely determined by this argument.
To represent a nucleon, a classical Skyrmion must be quantized. The unit Skyrmion hedgehog ansatz ties spatial rotations to isospin rotations. Its collective coordinates include its position and orientation; rotational quantization of a unit Skyrmion gives the rotor spectrum
in units with . The Finkelstein-Rubinstein constraints impose the correct fermionic sign under a nontrivial configuration-space loop. In particular a spatial rotation acts on a charge- state by in the physical odd-color theory: odd admits half-integer spin, while even has integer spin. For , the lowest allowed doublet represents the proton and neutron; the rotor state represents the Delta baryon resonance. A bosonic pion field can therefore describe fermionic baryons because the quantum wavefunction carries this nontrivial topological sign.
For nuclei, minimize the classical energy in a fixed baryon number sector, then quantize the permitted rotations, isospin rotations, and relevant vibrations or relative motions. The toroidal two-Skyrmion has a lowest nuclear state with , identifying it with the deuteron. The cubic four-Skyrmion has an allowed state appropriate to the alpha particle. The rational map approximation for Skyrmions makes these intrinsic symmetries easier to construct, while collective-rotation constraints for a Skyrmion select allowed nuclear quantum numbers. A spin-zero state has rotationally invariant laboratory expectation values; a classical cubic intrinsic field should not be interpreted as a fixed cube visible in every orientation. Collective-coordinate quantization restores this distinction between intrinsic shape and a physical quantum state.
The nuclear force also has a mesonic interpretation. At large separation the tails of Skyrmions are weak pion fields; with nonzero mass their multipole falloff derives from derivatives of the Yukawa potential. Their interaction depends on relative orientation, and after quantization generates the familiar spin- and isospin-dependent pion-exchange structure of the nuclear force. Attractive channels allow several unit Skyrmions to form a lower-energy multi-Skyrmion. In nuclear language the positive nuclear binding energy is the difference between the separated nucleon masses and the mass of the quantized bound state, not just a count of topological units.
The limitations remain physical. The simplest Skyrme model retains only selected terms in a derivative expansion, and finite solitons probe gradients where omitted terms can matter. Its parameters require matching or calibration; predicted binding can be too strong, and masses, radii and spectra are not all fixed correctly by topology. Rotational quantization alone neglects quantum and vibrational corrections, especially when clustering or breakup channels are important. More general mesonic interactions, additional meson fields, and less restrictive classical ansätze can improve the description, but they introduce further low-energy information. The organizing relation is therefore
It links underlying quark and gluon dynamics to a geometric, symmetry-based account of baryons and nuclei, while keeping the distinction between an effective approximation and a full derivation from QCD.
For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate . Let be specific thermal internal energy and follow a mass element. The stellar energy balance equation is
The outward luminosity is energy transported through the mass shell, is nuclear rest-mass energy released per unit mass and time, and is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by . At fixed composition it is by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in .
In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives , so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is . After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.
Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is , including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:
For nonresonant charged-particle reactions the energy weighting contains ; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as or are local approximations.
For the proton–proton chain, first produce deuterium and helium-3:
The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:
In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:
In pp III, the beryllium-7 instead captures a proton:
Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, is a useful local approximation.
The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists of
The carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through , , , and , after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is with a much steeper exponent than the PP chain; the preceding homology problem specifies . The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.
Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:
The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is , releasing about . Because three alpha particles are required, the specific rate scales as ; a schematic resonant rate is
Its local temperature exponent is , of order forty near . A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing reaction determines much of the eventual carbon-oxygen mixture.
If the core becomes hot enough, carbon burning follows, for example through and . Neon burning begins with photodisintegration, , followed by alpha capture such as . The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.
The increasing nuclear binding energy per nucleon supplies progressively less energy as products approach the iron group. Further fusion past that region cannot provide sustained net heat to support a core. Massive stars can therefore end with an unstable iron-group core; lower-mass stars do not attain all these burning stages and instead leave remnants such as carbon-oxygen white dwarfs. Heavy-element neutron captures can synthesize nuclei beyond the iron group without being the principal hydrostatic power source.
Finally, stellar neutrino energy loss competes with every heat source. Reaction neutrinos accompany the proton–proton and CNO chains. Hot or dense material also emits neutrino pairs through , plasmon decay, , and electron-ion bremsstrahlung. These neutrinos usually escape, unlike the photons whose transport is diffusive; neutrino trapping requires much more extreme collapse conditions. Late burning has smaller fuel-energy reservoirs and strong neutrino cooling, so its duration is much shorter than core hydrogen burning. The surface luminosity is the escaping photon power, , not the sum of photons and unobserved neutrino power. Radiative diffusion in a star, convection and electron conduction redistribute energy, while the energy-balance equation distinguishes true sources, losses and storage. Accretion or mass loss adds boundary energy and mechanical work when present, rather than changing the nuclear reaction budget.
Silicon burning 2026-10-06
At very high temperatures, photodisintegration and proton, neutron and alpha captures rearrange silicon-group material through networks approaching quasi-statistical equilibrium. The result is an iron-group mixture determined by electron fraction and weak-reaction timescales. Direct silicon-plus-silicon fusion is not the usual mechanism. The approach to maximal nuclear binding energy per nucleon ends sustained hydrostatic fusion power.