In short

Nuclear fusion is the process in which two light nuclei combine to form a heavier nucleus plus a particle or two, and in so doing release energy. The key reactions:

  • Proton–proton chain (the Sun): net 4\,^1\text{H} \to {}^4\text{He} + 2 e^+ + 2\nu_e + \gamma, releasing about 26.7 MeV per helium nucleus made.
  • Deuterium–tritium (D–T) fusion (the target of every terrestrial reactor): ^2\text{D} + {}^3\text{T} \to {}^4\text{He} + n + 17.6\,\text{MeV}. The fusion reaction with the largest cross-section at achievable temperatures.

Fusion energy comes from exactly the same bookkeeping as fission, just running in the opposite direction on the binding-energy curve. Fission starts from heavy nuclei (uranium, B/A \approx 7.6 MeV) and moves toward the peak at iron (8.79 MeV). Fusion starts from light nuclei (hydrogen, B/A = 0 for a single proton; deuterium, 1.11 MeV) and moves toward the peak from below. The climb on the fusion side is steeper — which is why fusion releases roughly seven times more energy per nucleon than fission.

Fusion requires extreme conditions because every fusion pair is two positive nuclei that have to touch, and they repel each other hard. The Coulomb barrier for proton–proton is about 0.7\,\text{MeV}, corresponding to a classical temperature of 8 \times 10^9\,\text{K} — far hotter than the Sun. The Sun fuses at 1.5 \times 10^7\,\text{K} only because a small quantum-mechanical tail of its Maxwell–Boltzmann distribution tunnels through the barrier. This tunnelling rate is ferociously sensitive to temperature.

Magnetic confinement (the tokamak, doughnut-shaped magnetic field) and inertial confinement (laser compression of a frozen D–T pellet) are the two routes being pursued to turn fusion into electricity. India operates the Aditya-U and SST-1 tokamaks at the Institute for Plasma Research in Gandhinagar, and contributes 9.1% of the cost of ITER — the international tokamak under construction at Cadarache, France.

Look up on a clear day. The Sun is about 150 million kilometres away, which is roughly 8 light-minutes. Every photon hitting the cricket pitch at the Wankhede in Mumbai right now left the Sun eight and a half minutes ago. If you could follow any one of those photons backwards through its last few million years — photons in the Sun's dense interior bounce off charged particles so relentlessly that it takes something like 10^5 years for a photon of gamma-ray energy to diffuse from the core to the surface — you would find it in the centre of the Sun, in a region about 20% of the solar radius across, where four hydrogen nuclei were squeezing themselves into one helium nucleus. The 26.7 MeV released by that fusion event is what the photon carried to your face today.

Every second, the Sun's core converts roughly 6.2 \times 10^{11} kilograms of hydrogen into 6.16 \times 10^{11} kilograms of helium. The missing 4 \times 10^9 kilograms per second — four million tonnes — has become energy, via E = mc^2. That is the Sun's luminosity: 3.85 \times 10^{26} watts. It has been doing this for 4.6 billion years and has another 5 billion or so left before the hydrogen in the core runs out.

Everything about that description is physics. This article derives the pieces — why fusion releases so much energy, why the Coulomb barrier is the obstacle, why the Sun fuses at a temperature that should be classically impossible, and what the machines at Gandhinagar and Cadarache are trying to do to bring the same process to a power grid near you.

Where the energy comes from — the binding-energy curve, again

Plot the binding energy per nucleon B/A against mass number A. The curve rises very steeply on the light side — from B/A = 0 for a free proton (there is nothing to bind), past B/A = 1.11 MeV for deuterium, 2.83 MeV for tritium, a sharp spike to 7.07 MeV for ^4He, through a smooth rise to the peak at ^{56}Fe (8.79 MeV per nucleon) — then falls slowly on the heavy side all the way to uranium (7.6 MeV).

Binding energy per nucleon is "how much glue is holding each nucleon in place." A higher B/A means a more tightly bound, more stable nucleus. Any nuclear process that moves the constituents up the curve — to a higher B/A — releases energy; any process that moves them down consumes energy. Fission is "heavy \to middle," moving up the curve from the right. Fusion is "light \to middle," moving up the curve from the left.

Binding energy per nucleon showing where fusion and fission sitA horizontal axis labelled A from zero to two hundred fifty, a vertical axis labelled B over A in MeV from zero to ten. A red curve starts at zero for A equals one, rises sharply through deuterium, helium-four at seven MeV, carbon, oxygen, peaks at iron near A equals fifty-six at eight point seven nine MeV, then gently falls to uranium at about seven point six MeV. An arrow labelled fusion points from hydrogen rightward up the curve; an arrow labelled fission points from uranium leftward up the curve. Both arrows meet near the iron peak. A (mass number) B/A (MeV) 0 3 6 9 50 100 150 200 250 ⁴He (7.07) ⁵⁶Fe (8.79, peak) ²³⁸U (7.57) fusion → ← fission
Binding energy per nucleon $B/A$ versus mass number $A$. Fusion moves from the far left (free hydrogen, $B/A = 0$) toward the peak at iron. The steep rise on the light side is why fusion releases so much energy per nucleon: from $^2$D ($B/A = 1.11$) to $^4$He ($B/A = 7.07$) is an increase of nearly 6 MeV per nucleon, versus roughly 0.9 MeV per nucleon for the uranium-to-middle-mass fission jump.

The fusion side of the curve is the steeper side. Climbing from hydrogen to helium is worth about 6 MeV per nucleon; climbing from uranium down toward middle masses is worth only 0.9 MeV per nucleon. This is the fundamental statement that fusion should, gram for gram, release roughly seven times more energy than fission. Everything else is engineering.

Energy of D–T fusion from the mass defect

Compute the energy release of the central terrestrial fusion reaction, deuterium + tritium:

^2_1\text{D} + {}^3_1\text{T} \;\longrightarrow\; {}^4_2\text{He} + {}^1_0 n + Q.

The atomic masses (in unified atomic mass units, u):

Compute the mass defect:

\Delta m \;=\; [2.014102 + 3.016049] - [4.002602 + 1.008665] \;=\; 5.030151 - 5.011267 \;=\; 0.018884\ \text{u}.

Why: atomic masses cancel the electron count (D and T each have one electron, ^4He has two, and that pairing works out on both sides). So you can use atomic masses directly here without tracking electrons, which is the conventional shortcut for fusion-reaction energy calculations.

Convert to MeV using 1\ \text{u}\,c^2 = 931.494\ \text{MeV}:

Q \;=\; \Delta m \cdot c^2 \;=\; 0.018884 \times 931.494 \;=\; 17.59\ \text{MeV}. \tag{1}

Why: this is the energy released per D–T fusion event. Most of it (14.1 MeV) goes to the neutron; the remaining 3.5 MeV goes to the alpha particle. The mass of the products is smaller than the mass of the reactants, and the shortfall has left the reaction as kinetic energy of the \alpha and the n. This is the Einstein bookkeeping.

Per nucleon: 17.59\ \text{MeV}/5 = 3.52\ \text{MeV} per nucleon. Compare with uranium fission, which releases 200/236 \approx 0.85\ \text{MeV} per nucleon. D–T fusion gives roughly four times more energy per nucleon than uranium fission — and for pure hydrogen fusion (as in the Sun) it is even higher.

Why fusion needs extreme temperatures — the Coulomb barrier

Fusion looks easy on paper. You just put two hydrogen nuclei close enough that the short-range nuclear force takes over and glues them together. The problem is "close enough." Two nuclei are both positively charged; they repel each other electrostatically. To get them within the nuclear-force range (about r_0 \approx 2\ \text{fm} = 2 \times 10^{-15}\ \text{m}), you have to push them against a mountain of Coulomb energy.

Deriving the barrier height

The Coulomb potential energy of two point charges Z_1 e and Z_2 e at separation r is

U_\text{C}(r) \;=\; \frac{1}{4\pi\varepsilon_0}\,\frac{Z_1 Z_2 e^2}{r}. \tag{2}

Why: the standard Coulomb potential-energy formula, applied to the two nuclei. Positive energy (repulsive) since both charges have the same sign.

The nuclei stop being point-like and the nuclear (strong) force takes over at r \approx r_0 \approx R_1 + R_2, the sum of the two nuclear radii. For proton–proton, both radii are about 1.0\ \text{fm}, so r_0 \approx 2\ \text{fm}.

Plug in numbers for proton–proton (Z_1 = Z_2 = 1):

U_\text{C}(r_0) \;=\; \frac{(9 \times 10^9)(1)(1.6 \times 10^{-19})^2}{2 \times 10^{-15}} \;=\; \frac{9 \times 10^9 \times 2.56 \times 10^{-38}}{2 \times 10^{-15}}
U_\text{C}(r_0) \;=\; \frac{2.304 \times 10^{-28}}{2 \times 10^{-15}} \;=\; 1.152 \times 10^{-13}\ \text{J}. \tag{3}

Why: this is the height of the Coulomb barrier in SI units. Plug in the two fundamental charges and the relevant distance, do the arithmetic.

Convert to electron-volts using 1\,\text{eV} = 1.6 \times 10^{-19}\ \text{J}:

U_\text{C}(r_0) \;=\; \frac{1.152 \times 10^{-13}}{1.6 \times 10^{-19}}\ \text{eV} \;\approx\; 720\ \text{keV} \;=\; 0.72\ \text{MeV}. \tag{4}

Why: convert to the energy unit of nuclear physics. About 0.7 MeV is the height of the Coulomb barrier for proton–proton — a number worth memorising.

For D–T fusion, both charges are 1 (both are hydrogen isotopes) and r_0 \approx 3.5\ \text{fm}, giving a slightly lower barrier of about 0.4\ \text{MeV}. For heavier nuclei the barrier grows as Z_1 Z_2/r_0 \propto Z_1 Z_2/(A_1^{1/3} + A_2^{1/3}), so fusion between two carbon nuclei (as in advanced stellar burning) sees a barrier of about 8\ \text{MeV} — ten times higher.

Classical temperature estimate

If every particle needs kinetic energy E \sim U_\text{C}(r_0) to climb the barrier, and the average kinetic energy in a thermal plasma is \tfrac{3}{2} k_B T, then the required temperature is

T \;\sim\; \frac{2 U_\text{C}}{3 k_B} \;=\; \frac{2 \times 1.15 \times 10^{-13}}{3 \times 1.38 \times 10^{-23}} \;\approx\; 5.6 \times 10^{9}\ \text{K}. \tag{5}

Why: set the average thermal kinetic energy equal to the barrier height, solve for temperature. This is a direct consequence of the Maxwell–Boltzmann distribution's mean.

That is five billion kelvin. The Sun's core is only 1.5 \times 10^7 kelvin — four hundred times cooler than this estimate. If fusion required each individual nucleus pair to climb the barrier classically, the Sun would not shine, and neither would any other star. Something else is going on.

The Sun cheats — quantum tunnelling through the barrier

A proton approaching another proton does not behave like a ball rolling toward a hill. It behaves like a quantum wavefunction, spread out over a region about one de Broglie wavelength across. When that wavefunction overlaps the far side of the barrier, there is a finite probability the proton will "tunnel through" and end up on the other side without ever having had the classical energy to climb.

This is the same tunnelling that alpha decay exploits to let alpha particles escape heavy nuclei. The probability of tunnelling through a Coulomb barrier at centre-of-mass energy E is the Gamow factor:

P(E) \;=\; \exp\!\left(-\frac{b}{\sqrt{E}}\right), \qquad b \;=\; \pi Z_1 Z_2 e^2 \sqrt{2\mu}/(4\pi\varepsilon_0 \hbar), \tag{6}

where \mu is the reduced mass of the two nuclei. The constant b encapsulates the barrier. Its important feature is that it is large — so the tunnelling probability is very small. For proton–proton at the Sun's core temperature, P \sim 10^{-9} per collision.

Fusion rate has two competing factors. Higher E gives a bigger tunnelling probability (the Gamow factor \exp(-b/\sqrt{E}) grows sharply with E). But fewer particles sit at higher E (the Maxwell–Boltzmann distribution is \propto \exp(-E/k_B T), falling with E). Their product,

\text{reaction rate} \;\propto\; \exp\!\left(-\frac{E}{k_B T} - \frac{b}{\sqrt{E}}\right),

peaks at a specific energy E_\text{G} called the Gamow peak — the energy at which most stellar fusion actually happens. For the Sun's core (k_B T \approx 1.3\ \text{keV}), E_\text{G} \approx 5-6 keV, still two orders of magnitude below the 720 keV classical barrier. The Sun fuses almost entirely on the high-energy tail of its thermal distribution, tunnelling through a barrier whose top it cannot classically reach.

The full calculation is in the going-deeper section. Here is the qualitative picture: the reaction rate depends on temperature through a steep exponential in T^{-1/3}. At T = 1.5 \times 10^7 K the Sun's core makes just enough fusions per second to power the Sun. Raise the temperature by a factor of two and the fusion rate would go up by a factor of about 2 \times 10^4. Halve the temperature and fusion would essentially stop. This extreme sensitivity is what makes stellar cores self-regulating — any local increase in temperature accelerates fusion and pushes the gas outward (cooling it), any local decrease chokes fusion and lets gravity recompress the core.

Explore — reaction rate versus temperature

Interactive: relative D–T fusion reactivity versus plasma temperature A curve showing relative reactivity on the vertical axis versus temperature in keV on the horizontal axis from zero to thirty. The curve starts near zero at low T, rises steeply past ten keV, and reaches its maximum around sixty-five keV before slowly declining. Readouts display current T, relative reactivity, and the thermal energy k T in kelvin equivalent. T (keV) ⟨σv⟩ (rel. to peak) 0 0.5 1 10 20 30 peak near 60 keV drag red dot to change plasma temperature
The D–T reactivity $\langle\sigma v\rangle$ — the velocity-averaged fusion cross-section — plotted against plasma temperature in keV. For reference, $1\ \text{keV}$ corresponds to $1.16 \times 10^7\ \text{K}$. The Sun's core is at about $1.3\ \text{keV}$ (off the bottom of this plot); a fusion reactor needs $T \approx 10$–$20\ \text{keV}$ ($100$–$200$ million kelvin). The steep rise past 5 keV is the Gamow factor; the peak near 60 keV and the slow decline at very high $T$ reflect the decrease of cross-section when nuclei pass each other too quickly to tunnel.

Drag the slider. The take-away: terrestrial fusion reactors aim for the steep rise, not the peak. The plasma density and confinement time at 10 keV are already as hard as any engineering challenge on the planet; pushing to 60 keV for a few extra units of reactivity is not worth it.

The proton–proton chain — how the Sun makes helium

The Sun fuses four protons into one helium nucleus. It does not do this in one step — four-body collisions are too rare. Instead, it takes a sequence of two-body steps called the proton–proton chain (pp-I):

Step 1. Two protons fuse into a deuteron (with the emission of a positron and a neutrino):

^1\text{H} + {}^1\text{H} \;\longrightarrow\; {}^2\text{D} + e^+ + \nu_e + 0.42\ \text{MeV}.

One of the protons has to beta-plus-decay into a neutron during the collision (a proton becomes a neutron + positron + neutrino). This is a weak interaction process, which is why it is rare. The mean lifetime of a proton in the Sun's core before it undergoes this fusion is about 10^{10} years. This is the bottleneck that controls the Sun's entire luminosity.

Why: a single fast process in a chain does not set the rate — the slowest step does. The pp fusion of two protons is the slowest step by an enormous margin, which is why the Sun burns steadily over billions of years rather than in a flash.

Step 2. The deuteron captures another proton to form ^3He:

^2\text{D} + {}^1\text{H} \;\longrightarrow\; {}^3\text{He} + \gamma + 5.49\ \text{MeV}.

This is a strong-interaction capture, so it is fast — a deuteron in the Sun's core lives only about one second before being captured.

Step 3. Two ^3He nuclei fuse into a ^4He plus two protons:

^3\text{He} + {}^3\text{He} \;\longrightarrow\; {}^4\text{He} + 2\,^1\text{H} + 12.86\ \text{MeV}.
Proton-proton chain showing the three stepsA diagram showing four protons entering, two fusing to deuterons through beta-plus emission, each deuteron combining with another proton to make helium-3, and two helium-3 nuclei combining to make helium-4 plus two protons returned. Step 1 (slow, weak) p p → ²D + e⁺ + νₑ 0.42 MeV Step 2 (fast, strong) ²D p → ³He + γ 5.49 MeV Step 3 (strong) ³He ³He → ⁴He + 2 p 12.86 MeV
The three steps of the pp-I branch. Step 1 is rate-limiting — a proton in the Sun's core waits about $10^{10}$ years before it fuses. Steps 2 and 3 are fast. Net result of the sequence (Step 1 done twice + Step 2 done twice + Step 3): four protons become one $^4$He plus two positrons, two neutrinos, and two gammas.

Net reaction (Step 1 happens twice, Step 2 happens twice, Step 3 happens once):

4\,^1\text{H} \;\longrightarrow\; {}^4\text{He} + 2 e^+ + 2\nu_e + 2\gamma + 26.73\ \text{MeV}.

The positrons annihilate with ambient electrons, adding 4 m_e c^2 \approx 2.04\ \text{MeV} of extra gamma rays. The neutrinos escape the Sun almost immediately carrying about 0.59\ \text{MeV} between them. The remaining 26.14\ \text{MeV} is deposited as thermal energy in the core — this is what heats the Sun and eventually leaks out as the light hitting your face.

Worked examples

Example 1: Solar luminosity from the pp chain

The Sun radiates L_\odot = 3.85 \times 10^{26}\ \text{W}. If every bit of that luminosity comes from the pp chain (which deposits \approx 26.7\ \text{MeV} per ^4He made, minus the neutrinos, so about 26.1\ \text{MeV} to the core), how many helium nuclei must the Sun's core make per second? And how many kilograms of hydrogen does that consume per second?

Budget diagram for solar luminosityA flow diagram. Sixty-two trillion kilograms of hydrogen per second flow in from the left. An arrow points to a central box labelled pp chain. From the box, an arrow labelled helium output points to the right showing four trillion kilograms of mass becomes energy at three point eight five times ten to the twenty-six watts. H fuel 6.2×10¹¹ kg/s pp chain 4 p → ⁴He 26.7 MeV L = 3.85×10²⁶ W Δm = 4×10⁹ kg/s
Every second, $6.2 \times 10^{11}$ kg of hydrogen enters the pp chain in the Sun's core. A fraction ($4 \times 10^9$ kg/s) is converted into pure energy at the luminosity $L_\odot = 3.85 \times 10^{26}$ W. The rest comes out as helium.

Step 1. Convert the energy per fusion event to joules.

26.1\ \text{MeV} \;=\; 26.1 \times 10^6 \times 1.6 \times 10^{-19}\ \text{J} \;=\; 4.18 \times 10^{-12}\ \text{J}.

Why: the luminosity is in SI (watts = J/s), so the per-event energy must be in joules to compare.

Step 2. Divide the luminosity by the energy per event.

N_{\text{He}} \;=\; \frac{L_\odot}{E_{\text{per He}}} \;=\; \frac{3.85 \times 10^{26}}{4.18 \times 10^{-12}} \;\approx\; 9.2 \times 10^{37}\ \text{helium nuclei per second}.

Why: this is just "power equals number of events times energy per event." Rearrange to get the rate.

Step 3. Each ^4He uses 4 protons. So the rate of hydrogen consumption is

N_H \;=\; 4 N_{\text{He}} \;=\; 3.68 \times 10^{38}\ \text{protons per second}.

Convert to kilograms using the proton mass m_p = 1.673 \times 10^{-27} kg:

\text{mass rate} \;=\; 3.68 \times 10^{38} \times 1.673 \times 10^{-27} \;=\; 6.16 \times 10^{11}\ \text{kg/s}.

Why: multiply the number of protons per second by the mass per proton to get kilograms of hydrogen consumed per second. The answer is 620 billion kilograms — six hundred megatonnes of hydrogen fused every second.

Step 4. Mass-to-energy conversion rate. The luminosity divided by c^2 gives the rate at which mass becomes energy:

\dot m \;=\; \frac{L_\odot}{c^2} \;=\; \frac{3.85 \times 10^{26}}{(3 \times 10^8)^2} \;=\; \frac{3.85 \times 10^{26}}{9 \times 10^{16}} \;=\; 4.3 \times 10^{9}\ \text{kg/s}.

Why: this is the amount of mass that actually disappears per second. It is far less than the 6.2 \times 10^{11} kg of hydrogen consumed — because most of the hydrogen mass is retained as helium; only a 0.7% fraction of each fused hydrogen batch is converted to energy.

Result: The Sun makes about 9.2 \times 10^{37} helium nuclei per second, consuming 6.2 \times 10^{11} kg of hydrogen per second. About 4.3 \times 10^9 kg of that becomes energy; the remaining 6.16 \times 10^{11} kg stays behind as helium. Doing this for 4.6 billion years, the Sun has converted roughly 9 \times 10^{28} kg of hydrogen into helium — about 5% of its initial hydrogen inventory.

What this shows: The Sun is a very inefficient mass-to-energy engine (0.7%) but its enormous hydrogen reservoir keeps it burning for billions of years. A fusion power reactor running at 1000\ \text{MW} electrical and 40% conversion efficiency would burn about 110\ \mu\text{g} of D–T mixture per second — the Sun's mass-burn rate scaled down by a factor of 4 \times 10^{13}.

Example 2: Energy from one gram of D–T

One gram of a 50–50 (by atom count) mixture of deuterium and tritium — 3 \times 10^{23}/5 = 6 \times 10^{22} pairs, approximately — undergoes complete fusion. What is the total energy released, expressed in joules and in kilowatt-hours?

Energy yield from one gram of D-T fusion fuelA bar chart comparing three fuels by energy released per gram. One gram of coal gives 24 kilojoules. One gram of uranium-235 fission gives 82 billion joules. One gram of D-T fusion gives about 340 billion joules. The fusion bar is the tallest. Energy per gram (log scale) coal 24 kJ ²³⁵U fission 82 GJ D–T fusion 340 GJ 14× fission 10⁷× coal
Energy per gram of fuel. Coal: $\sim 24\ \text{kJ/g}$ (chemical combustion). $^{235}\text{U}$ fission: $\sim 82\ \text{GJ/g}$. D–T fusion: $\sim 340\ \text{GJ/g}$. The fusion bar is four times taller than the fission bar, ten million times taller than the coal bar.

Step 1. Count D–T pairs in one gram.

The mixture is a 50–50 atom mixture of ^2D (mass 2.014 u) and ^3T (mass 3.016 u). The average atomic mass is (2.014 + 3.016)/2 = 2.515\ \text{u}. So one gram contains

N_\text{atoms} \;=\; \frac{1\ \text{g}}{2.515\ \text{u} \times 1.66 \times 10^{-24}\ \text{g/u}} \;=\; \frac{1}{4.175 \times 10^{-24}} \;\approx\; 2.40 \times 10^{23}\ \text{atoms}.

Half of these are D and half are T, so the number of D–T pairs is N_\text{atoms}/2 = 1.20 \times 10^{23}.

Why: each fusion event takes one deuteron plus one triton. The number of possible events is the number of pairs, not the number of atoms.

Step 2. Multiply by energy per fusion from equation (1).

E_\text{total} \;=\; N_\text{pairs} \times 17.59\ \text{MeV} \;=\; 1.20 \times 10^{23} \times 17.59 \times 10^6 \times 1.6 \times 10^{-19}\ \text{J}.

Compute:

E_\text{total} \;=\; 1.20 \times 10^{23} \times 2.814 \times 10^{-12} \;=\; 3.38 \times 10^{11}\ \text{J} \;\approx\; 340\ \text{GJ}.

Why: 340 gigajoules from one gram of fuel. For comparison, burning one gram of petrol gives about 46 kilojoules — a factor of 7 \times 10^6 less.

Step 3. Convert to kilowatt-hours (1\ \text{kWh} = 3.6 \times 10^6\ \text{J}):

E_\text{total} \;=\; \frac{3.38 \times 10^{11}}{3.6 \times 10^6} \;\approx\; 9.4 \times 10^{4}\ \text{kWh}.

Why: convert to a unit your electricity bill speaks. Ninety-four thousand kilowatt-hours from one gram of D–T, assuming you could capture 100% of the energy (in practice a fusion reactor would capture about 35–40% as electricity).

Step 4. Household interpretation.

A typical Indian household consumes about 300\ \text{kWh} per month. One gram of D–T fusion fuel would therefore supply that household for

\frac{9.4 \times 10^4}{300} \;\approx\; 310\ \text{months} \;\approx\; 26\ \text{years}.

Result: One gram of D–T fusion releases 3.4 \times 10^{11} J, equivalent to 9.4 \times 10^4 kWh — enough to run a typical Indian household for 26 years, or (at 40% conversion efficiency) about 10 years of real grid power.

What this shows: The energy density of fusion fuel is astonishing. The hard problems are not the fuel — seawater contains 33\ \text{mg} of deuterium per litre, a millennia-scale resource — but the plasma confinement. The Aditya-U tokamak at IPR Gandhinagar typically holds its plasma for a few tenths of a second, at densities and temperatures around 100 million kelvin; scaling that to a commercial reactor needs confinement of seconds to minutes at the same temperatures.

Controlled fusion — magnetic and inertial confinement

A fusion reactor has three numbers to maximise simultaneously: plasma density n, confinement time \tau_E, and temperature T. Their product, the Lawson triple product nT\tau_E, must exceed about 3 \times 10^{21}\ \text{keV}\cdot\text{s}/\text{m}^3 for a D–T reactor to produce more power than it consumes. This is called the ignition criterion, and no reactor on Earth has yet achieved it in a sustained way (though the ignition milestone at the National Ignition Facility in December 2022 crossed it for a single shot, and ITER is projected to do so in the 2030s).

Two strategies compete:

Magnetic confinement. A plasma at 10^8\ \text{K} will instantly vaporise any solid container. Instead, use a magnetic bottle — charged particles (nuclei and electrons) spiral along magnetic field lines, so a magnetic field that closes on itself can trap them. The tokamak is a doughnut (torus) shaped vessel with a helical magnetic field, and it is the leading candidate reactor design. The plasma is kept in the middle of the doughnut, never touching the walls.

India operates two tokamaks at the Institute for Plasma Research, Gandhinagar:

And India is a 9.1%-stakes member of ITER at Cadarache, France, supplying (among other things) the cryostat — the world's largest stainless-steel vacuum vessel, built by Larsen & Toubro at their Hazira facility in Gujarat.

Inertial confinement. A different strategy: compress a millimetre-scale pellet of frozen D–T with 192 ultra-high-power lasers (as at NIF in the United States), so briefly and intensely that the fuel fuses before its own inertia lets it fly apart. You don't need to contain the plasma; you just need it to fuse in the nanosecond before it disassembles itself. India has a small inertial-fusion programme at BARC using lasers at the Raja Ramanna Centre for Advanced Technology in Indore, focused on inertial-confinement physics studies.

Tokamak schematic showing toroidal vessel and magnetic field linesA cross section of a torus-shaped tokamak. A rectangular outer boundary marks the vacuum vessel. Inside, a plasma ring is shown. Arrows going around the torus represent the toroidal magnetic field. A helical arrow represents the net field lines that trap the plasma. vacuum vessel plasma (10⁸ K) helical B field toroidal B
Tokamak schematic. A torus-shaped vacuum vessel (dashed rectangle in cross-section) contains a hot plasma (red). The magnetic field inside is helical — the *toroidal* component (red arrows) runs around the big loop of the doughnut; an additional *poloidal* component (not shown) wraps around the small loop. Together they twist the field lines into helices that trap charged particles. The field is produced by external coils (also not shown) that carry hundreds of kiloamperes of current.

Common confusions

If you came here to understand the Sun and the tokamak at a conceptual level, you have what you need. What follows is the quantitative version of the Gamow tunnelling calculation, the Lawson criterion, and the CNO cycle.

The Gamow factor from tunnelling

The probability that a particle of energy E tunnels through a Coulomb barrier U_\text{C}(r) = (Z_1 Z_2 e^2)/(4\pi\varepsilon_0 r) is, in the WKB approximation,

P \;=\; \exp\!\left(-\frac{2}{\hbar}\int_{r_0}^{b}\sqrt{2\mu[U_\text{C}(r) - E]}\, dr\right),

where r_0 is the nuclear-contact radius and b = Z_1 Z_2 e^2 /(4\pi\varepsilon_0 E) is the classical turning point (where U_\text{C} = E). For E \ll U_\text{C}(r_0) (always the case at stellar temperatures), r_0 contributes negligibly to the integral, and the result reduces to

P \;=\; \exp\!\left(-\frac{\pi Z_1 Z_2 e^2}{4\pi\varepsilon_0\, \hbar}\sqrt{\frac{2\mu}{E}}\right) \;=\; \exp\!\left(-\frac{b_G}{\sqrt{E}}\right),

with

b_G \;=\; \pi Z_1 Z_2 e^2 \sqrt{2\mu}/(4\pi\varepsilon_0\, \hbar).

Why: the WKB integral evaluates to a beta-function-like expression that reduces, in the low-energy limit, to \pi times the square-root ratio. This is the Gamow factor; its key feature is the 1/\sqrt{E} in the exponent, which makes tunnelling rise very steeply with energy.

Gamow peak. The fusion rate per unit volume at temperature T is the Maxwell–Boltzmann distribution times the tunnelling probability:

\text{rate}(E) \;\propto\; E\, \exp\!\left(-\frac{E}{k_B T}\right) \exp\!\left(-\frac{b_G}{\sqrt{E}}\right).

Maximise the argument of the combined exponential — differentiate -E/k_B T - b_G/\sqrt{E} with respect to E, set to zero:

-\frac{1}{k_B T} + \frac{b_G}{2 E^{3/2}} \;=\; 0 \quad\Longrightarrow\quad E_\text{G} \;=\; \left(\frac{b_G k_B T}{2}\right)^{2/3}.

Why: this is the Gamow peak — the specific energy at which most fusion events occur. For the Sun's core (k_B T = 1.3\ \text{keV} and b_G for p-p), E_\text{G} \approx 5.9\ \text{keV}, roughly 5 k_B T — well into the high-energy tail of the thermal distribution, but still far below the 720 keV classical barrier.

Temperature dependence. Substitute E_\text{G} back to find the fusion rate's scaling with T:

\text{rate}(T) \;\propto\; T^{-2/3}\exp\!\left(-\frac{3 E_\text{G}}{k_B T}\right) \;=\; T^{-2/3}\exp\!\left(-\frac{3}{(2)^{2/3}}\frac{b_G^{2/3}}{(k_B T)^{1/3}}\right).

This is the famous \exp(-\text{const}/T^{1/3}) temperature dependence of thermonuclear fusion. Differentiating \ln(\text{rate}) with respect to \ln T gives the local exponent n in \text{rate} \propto T^n:

n \;=\; \frac{d\ln(\text{rate})}{d\ln T} \;=\; \frac{E_\text{G}}{k_B T} - \frac{2}{3}.

For p-p at T = 1.5 \times 10^7\ \text{K} (k_B T \approx 1.3\ \text{keV}, E_\text{G} \approx 5.9\ \text{keV}): n \approx 5.9/1.3 - 2/3 \approx 3.9. The Sun's pp-chain luminosity goes as T^4. For the CNO cycle (which dominates in heavier stars), E_\text{G} is larger because the barrier is higher, and n can reach 15–20 — so CNO burning is much more temperature-sensitive than pp burning. This is why the Sun burns pp but a star 50% more massive burns CNO.

The Lawson criterion

For a D–T reactor to break even, the fusion power must exceed the power lost via bremsstrahlung, conduction, and convection. A simple bookkeeping gives the Lawson criterion:

n\tau_E \;\geq\; \frac{12 k_B T}{\langle\sigma v\rangle\, E_f}, \tag{7}

where n is the ion density, \tau_E the energy confinement time, \langle\sigma v\rangle the reactivity at temperature T, and E_f the fusion energy per event. At T = 1020\ \text{keV} (the optimum for D–T), \langle\sigma v\rangle \approx 10^{-22}\ \text{m}^3/\text{s} and E_f = 17.6\ \text{MeV}, giving

n\tau_E \;\geq\; 10^{20}\ \text{s}/\text{m}^3 \quad\text{(ignition)}.

The triple product n T \tau_E \geq 3 \times 10^{21}\ \text{keV}\cdot\text{s}/\text{m}^3 is the more commonly cited form. Tokamak progress over the last fifty years is tracked on a single plot of this triple product versus year — it has risen by a factor of 10^5 since the 1970s and is projected to cross ignition at ITER in the 2030s.

Two extreme limits of the Lawson criterion:

  • Magnetic confinement uses low density (n \sim 10^{20}\ \text{m}^{-3}, about 10^{-6} of atmospheric) and long confinement time (\tau_E \sim 1 s).
  • Inertial confinement uses enormous density (n \sim 10^{31}\ \text{m}^{-3}, 10^5 times solid D–T density) and brief confinement (\tau_E \sim 10^{-11} s).

Both can in principle reach the same n\tau_E product. Neither has yet sustained it.

The CNO cycle

In stars more massive than about 1.3\,M_\odot, the carbon–nitrogen–oxygen (CNO) cycle takes over from the pp chain as the dominant hydrogen-burning process. It uses pre-existing ^{12}C as a catalyst: the catalyst is reused in a six-step loop, not consumed. The net reaction is still 4\,^1\text{H} \to {}^4\text{He} with the same 26.7 MeV energy release, but the rate is much more temperature-sensitive (\propto T^{17} at T = 1.5 \times 10^7\ \text{K}). This higher temperature sensitivity is why massive stars burn out quickly — their higher core temperatures make fusion fast, exhausting the fuel in tens of millions of years rather than the Sun's ten billion.

Why D–T and not D–D

The D–D reaction has two branches of roughly equal probability:

  • ^2\text{D} + {}^2\text{D} \to {}^3\text{T} + p + 4.03\ \text{MeV}
  • ^2\text{D} + {}^2\text{D} \to {}^3\text{He} + n + 3.27\ \text{MeV}

D–D fuel is cheaper (tritium has to be bred from lithium, an extra engineering burden), but the cross-section at 10–20 keV is about a hundred times smaller than D–T, putting its Lawson criterion an order of magnitude further out of reach. First-generation reactors will run D–T; if the fusion industry matures, later reactors may switch to D–D, or to aneutronic fuels like D–^3He or p–^{11}B. Both of the advanced fuels have vanishing neutron emission — making the reactor much easier to maintain — but vastly higher temperature requirements.

Where this leads next