Why can a star shine for so long, and why do hotter stellar cores favor different nuclear reactions? A useful theory must supply both an adequate energy reservoir and the right stellar energy-generation rate. Nuclear binding energy supplies the reservoir; the thermonuclear reaction rate, stellar composition and stellar structure determine how quickly it is spent. The proton–proton chain and CNO cycle are two routes for hydrogen burning with very different temperature sensitivities.
In his 1967 lecture, Hans Albrecht Bethe contrasted the short Kelvin-Helmholtz cooling time with long ages inferred from terrestrial and stellar evidence.[ref]
For an original solar-scale estimate, choose mass , radius and luminosity . The characteristic Kelvin-Helmholtz cooling time is about ; the precise gravitational reservoir has a structure-dependent coefficient. If a fraction of the mass can burn its initial hydrogen mass fraction , with a chosen rest-energy conversion fraction , the fuel estimate is
This assumes fixed luminosity and neglects neutrino losses. It demonstrates the reservoir advantage, rather than predicting the Sun's age or its full stellar evolution. Having enough fuel still does not identify a reaction: the fuel must release energy at a rate compatible with hydrostatic equilibrium and radiative transfer.
Hans Albrecht Bethe and Charles Louis Critchfield calculated deuterium formation from two protons in 1938, combining quantum tunnelling with a weak conversion of a proton into a neutron. Their calculated nuclear reaction cross-section was a theoretical input to hydrogen burning.[ref] In 1939, Hans Albrecht Bethe analysed the catalytic carbon–nitrogen route and compared required stellar core temperatures with stellar structure estimates in Table VIII. His reported solar agreement used historical abundance and structure assumptions; it should not be read as a modern demonstration that the CNO cycle dominates the Sun.[ref]
Hans Albrecht Bethe's 1967 lecture credits Carl Friedrich von Weizsäcker's independent 1938 carbon–nitrogen cycle discovery, distinguishing it from calculation of rates and temperature dependence.[ref] The 1967 Nobel Prize in Physics recognized Hans Albrecht Bethe's contributions to nuclear reaction theory, especially stellar energy production.[ref]
Laboratory nuclear reaction cross-sections were inputs, not measurements of an artificial star. The 1939 paper converted thick-target yields to nuclear reaction cross-sections using a range–energy relation and tabulated capture data.[ref] A concrete apparatus example is the Cockcroft-Walton lithium disintegration experiment, with an accelerated proton beam, lithium target, zinc sulfide screen and beam-blocking control. Its nuclear reaction differs from both stellar burning routes. The 1967 lecture stresses extrapolation below laboratory access and tests for nearby nuclear levels.[ref]
If all mass of a spherical body lies within radius , then in the shell-assembly integral for its gravitational binding energy. Therefore . The inequality is strict for an extended interior mass distribution; a thin shell gives equality. Dividing by a stipulated luminosity gives a lower bound on that reservoir-to-luminosity ratio. A Kelvin-Helmholtz cooling time also needs a structure model and the accessible fraction of gravitational release.
Core accretion builds a solid core and then a massive gaseous envelope before protoplanetary disk dispersal. The solid and gas stages overlap; this is not direct collapse of a whole disk patch.
In a protoplanetary disk, refractory grains and, beyond a snow line, condensed ices provide solids. Settling and protoplanetary dust coagulation build aggregates, but fragmentation, bouncing and radial drift of solids in a gas disk can interrupt simple sticking. Concentration processes such as the streaming instability can produce dense particle structures that collapse into planetesimals. Subsequent collisions build planetary embryos; alternatively, embryos can efficiently capture drifting particles through gas-assisted pebble accretion. Neither the initial planetesimal population nor unrestricted growth through every size should be assumed without a mechanism.
For gravitationally focused growth of a planetary core, the two-body capture cross-section is
An illustrative thin-disk estimate is , up to geometry and velocity-distribution factors. At fixed solid density and random speed, strong gravitational focusing gives : runaway planetesimal growth. As large embryos stir the smaller bodies, encounter speeds increase and focusing weakens, giving oligarchic growth. Consuming the local feeding zone leads toward an isolation mass; radial supply, migration and mergers can extend growth. Snow lines can increase the solid reservoir, but a distant core also grows more slowly because orbital times are longer.[ref]
A growing core binds gas over a scale limited by the smaller of its Bondi radius and Hill radius. The envelope approximately satisfies
An equation of state, disk boundary conditions, and radiative or convective energy transport close the structure. In an optically thick radiative region,
Here is opacity; an unstable radiative gradient instead drives convection. Solid accretion releases roughly , heating and supporting the envelope. Gas capture consequently depends on cooling, not just whether the potential can bind gas instantaneously.
The envelope loses energy on a Kelvin-Helmholtz cooling time, . Lower grain opacity, higher core mass, or reduced solid-accretion heating can accelerate contraction. The critical core mass for gas accretion is therefore model dependent; “of order ten Earth masses” is a common scale, not a universal switch. As envelope self-gravity becomes important and the gas and solid masses become comparable, contraction permits still more gas capture, producing runaway gas accretion. The classical sequence has a rapid solid-growth phase, a slower envelope-growth phase, and accelerated gas accumulation.[ref]
Figure 1.
Core accretion requires both solid growth and envelope cooling while the disk still supplies gas
. Source. Details.
The decisive race is . Failure to cool or obtain a sufficiently large core before gas removal leaves a smaller solid or ice-rich planet. Successful growth is ultimately limited by disk supply, gap opening by a planet, and dispersal; a gap does not completely seal off accretion. Angular-momentum exchange drives planetary migration, so the observed orbit need not be the formation radius. Accretion-shock radiation and subsequent contraction set the young giant's entropy and luminosity; final core mass, envelope enrichment and orbital location depend on the whole growth history, not on one fixed threshold. The necessary cooling, supply and orbital constraints make the model testable rather than an assertion that every embryo must become a gas giant.
A dense cloud core first collapses when self-gravity overwhelms its effective support. While radiation escapes readily, compression can remain approximately isothermal and the Jeans mass decreases as density rises, encouraging gravitational fragmentation. Angular momentum usually leads to a disk rather than unrestricted spherical infall.
As the central optical depth increases, cooling becomes slower than compression and pressure produces a first hydrostatic core. molecular hydrogen dissociation near absorbs compressional energy and lowers the effective adiabatic response, permitting a second collapse of a protostar. After dissociation the gas stiffens and a much more compact stellar core forms. An accretion shock and surrounding disk then convert infall energy into radiation. A small central core forms before most of the envelope arrives, as demonstrated in radiation-hydrodynamic collapse calculations.[ref]
The isothermal inside-out collapse of a singular isothermal sphere model provides one idealized supply law, , but rotation, magnetic fields and an unsteady reservoir modify real accretion. Its luminosity scale is . Dust can absorb the emerging short-wavelength radiation and re-emit it in the infrared, so an embedded source need not have a directly visible stellar photosphere.
As envelope accretion declines, gravitational contraction becomes important on the Kelvin-Helmholtz cooling time . deuterium-burning energy-generation rate can temporarily influence the radius; sustained core hydrogen burning eventually establishes the main sequence if the final mass is sufficient. A low-mass convective object approaches along a Hayashi-like track, while more massive objects need not share that track. accretion, contraction and nuclear burning are distinct energy sources, and a protostar is not generally a fixed-mass object.
A cold self-gravitating core initially contracts on a dynamical scale if cooling prevents pressure from halting collapse. Increasing optical depth eventually traps radiation and forms a pressure-supported first hydrostatic core. Dissociation of molecular hydrogen absorbs compressional energy, softens the equation of state and triggers a second protostellar collapse to a compact protostar. Rotation produces a disk; angular-momentum transport permits continuing accretion, while outflows remove some mass and angular momentum.
An accreting object's luminosity includes a contribution of order , with a coefficient set by where accretion energy is radiated or retained. After accretion subsides, gravitational contraction supplies internal luminosity over the Kelvin-Helmholtz cooling time . Hydrostatic contraction heats the interior even while energy is radiated, as follows from the virial theorem and the negative heat capacity of a virialized gravitational system. accretion luminosity must not be substituted for in the contraction time without an energy-budget justification.
Low-mass, nearly fully convective pre-main-sequence stars contract down a Hayashi track at roughly constant effective temperature and decreasing radius and luminosity. A developing radiative core can lead onto a Henyey track toward higher effective temperature; the route depends on mass and opacity. The deuterium-burning energy-generation rate can temporarily moderate contraction, but sustained core hydrogen burning marks arrival on the main sequence. Massive objects may begin hydrogen burning while still accreting, and substellar objects never establish long-lived hydrogen burning. An infrared spectral class reflects envelope and disk emission as well as inclination and interstellar extinction, so it is not a unique stellar age or interior evolutionary track.
Nuclear burning explains both the large available energy and its release under stellar-core conditions. Converting hydrogen into helium releases approximately of the fuel rest-mass energy. The net proton–proton chain reaction can be written , with about MeV available per helium nucleus, including the neutrino energy loss. The CNO cycle achieves the same net hydrogen-to-helium conversion with catalysts. Measured nuclear reaction rates and stellar-core temperatures give rates capable of supplying the observed luminosity.
The energy budget distinguishes stellar nuclear fusion from alternatives. Kelvin-Helmholtz contraction has a characteristic Kelvin-Helmholtz cooling time , only tens of millions of years for a solar-type star. Geological ages and stellar-population ages are much longer. Chemical binding energy per unit mass is smaller still. Nuclear fuel in a substantial fraction of the core gives a stellar nuclear timescale of billions of years, explaining long-lived main-sequence stars.
Main sequence luminosities, the correlation of lifetime with mass, and cluster turnoffs are reproduced by models in which central hydrogen is consumed and helium accumulates. Spectroscopic abundances and the appearance of helium-, carbon- and heavier-element-rich evolved material agree with sequential nuclear processing and stellar dredge-up. These observations support the reaction network and evolutionary consequences, rather than independently measuring every reaction in every star.
Solar neutrinos provide direct information from the burning core because they escape without the long radiative diffusion delay of photons. Their measured energies and fluxes test weak reactions in the hydrogen-burning network. The historical deficit of electron neutrinos did not require abandoning nuclear power; flavor conversion changes the measured electron component. Later flavor-sensitive measurements confirmed the total flux expected from solar nuclear burning, as described in the SNO collaboration's first-results report. Thus the evidence combines an adequate energy source, quantitatively successful stellar models and a direct nuclear by-product. Cooling remnants and brief gravitationally powered evolutionary phases are exceptions to the statement that a star's instantaneous luminosity is nuclear.
The first argument is the energy budget. The gravitational energy reservoir is of order , so the Kelvin-Helmholtz cooling time is . For the Sun, , and give years. This is much shorter than the approximately -year solar-system age. Chemical binding energies are only of order electronvolts per atom and provide an even smaller reservoir. Thus neither ordinary chemical energy nor continuing gravitational contraction can supply the long-term solar luminosity.
Stellar nuclear fusion provides an adequate reservoir. Conversion of four hydrogen atoms into one helium-4 atom releases about , corresponding to a mass defect of about of the initial mass; some energy leaves as neutrinos. If a fraction of stellar mass with hydrogen fraction is available for core burning, the energy scale is . Taking and supplies roughly , enough for several billion years at the solar luminosity, consistent with a stellar nuclear timescale of order years. The core temperature and density inferred from stellar structure permit the proton–proton chain; hotter stellar cores favor the CNO cycle.
The most direct evidence is the detection of solar neutrinos. They are created in the weak reactions of the proton–proton chain and can escape from the core, unlike the photons whose transport through the interior is slow. Measurements accounting for neutrino oscillation recover a total neutrino flux consistent with nuclear-powered solar models. Their spectra and reaction-dependent fluxes test the core reactions rather than merely the surface effective temperature.
There is also independent evidence from stellar evolution: the main sequence is a long-lived core-hydrogen-burning phase, more massive luminous stars consume fuel more quickly, and the main-sequence turnoff of a star cluster records its age. Subsequent evolution and the production of heavier elements fit further nuclear reactions. Agreement among these observations, nuclear energy release and stellar-model lifetimes is much stronger than the energy-budget argument alone. The combination of a sufficient nuclear reservoir, measured core neutrinos and the observed evolutionary sequence identifies nuclear reactions as the main long-term power source of ordinary stars.
A star of roughly solar stellar composition begins on the zero-age main sequence with a hot convective core, a mostly radiative envelope, and hydrogen burning dominated by the CNO cycle. The strong temperature dependence of the CNO cycle concentrates the stellar energy-generation rate near the centre, helping produce the convective core. Core mixing replenishes hydrogen within that region; it does not homogenize the entire star. The nuclear fuel supply gives a stellar nuclear timescale of order years. For a numerical comparison, Table 2 of Ekström et al.'s solar-metallicity models gives core hydrogen burning lifetimes of about and million years for nonrotating and rotating models, respectively.
After central hydrogen exhaustion, an inert helium core contracts and heats through Kelvin-Helmholtz contraction, while a hydrogen-burning shell supplies the growing luminosity. The envelope expands and cools: the track crosses the Hertzsprung gap toward the red giant region. Its deepening convective envelope causes first dredge-up, bringing material altered by hydrogen burning to the surface. The rapid structural crossing has a thermal timescale, roughly – years for a star of this mass, rather than the full stellar nuclear timescale; an estimate follows from the Kelvin-Helmholtz cooling time with the evolving radius and luminosity.
The helium core reaches ignition before becoming strongly supported by electron degeneracy pressure. Thus a normal star starts core helium burning without the helium flash characteristic of lower-mass stars. The Triple-alpha process produces carbon-12, and carbon-12 alpha capture produces oxygen-16; these reactions build a carbon-oxygen core. There is a convective core and continued hydrogen-burning shell activity. Core helium burning lasts order years: the same model table gives approximately and million years. The changing core and envelope structure can produce a blue loop in the Hertzsprung-Russell diagram; its extent depends on stellar composition, mixing and stellar rotation.
After central helium exhaustion, the carbon-oxygen core contracts and becomes increasingly supported by electron degeneracy pressure. A helium-burning shell forms above it, beneath the hydrogen-burning shell. The star returns to a cool, luminous Asymptotic giant branch configuration. On this early branch, helium-burning shell activity is an important energy source; the overlying hydrogen-burning shell can weaken during readjustment and resume later. At this mass, a deepening convective envelope can cause second dredge-up. This shell-supported evolution to the first AGB thermal pulse takes order years, with model-dependent duration. As the helium-burning shell becomes sufficiently thin, the thin-shell instability ends steady shell burning and marks the onset of AGB thermal pulses.
Figure 1.
Schematic evolution of a five-solar-mass star to the first AGB thermal pulse
. Details.
The original sketch shows the direction of stellar evolution, with higher effective temperature to the left. It illustrates a possible blue loop; its coordinates and loop shape are schematic, rather than a computed stellar model. The stages before the first AGB thermal pulse are the requested track; the later remnant evolution is described below.
In the usual final evolution, stellar winds and an AGB superwind remove the envelope before the carbon-oxygen core can ignite substantial carbon burning. Once shell fuel is exhausted, the exposed core contracts and heats, potentially ionizing the expelled gas as a planetary nebula, and finally becomes a carbon-oxygen white dwarf that loses its stored thermal energy by cooling. The relevant competition is between the AGB envelope-loss and core-growth timescales, not just the total initial mass.
Much higher late-stage mass loss strips the envelope sooner, shortens the Asymptotic giant branch phase, and allows less core growth, usually producing a smaller carbon-oxygen white dwarf. Much lower mass loss prolongs shell burning and permits greater core growth. In an extreme hypothetical case, retention of the envelope could let a degenerate carbon-oxygen core approach the Chandrasekhar mass and undergo explosive carbon burning: the proposed type 1.5 supernova channel. This is not an inevitable fate for a star. Doherty et al.'s evolutionary calculations emphasize the uncertainty of that channel and find that their adopted mass loss removes carbon-oxygen envelopes before such ignition. If carbon burning instead proceeds without disruption, the possible further evolution depends on the resulting core and continued envelope retention; lowering mass loss alone does not establish an electron-capture supernova outcome.
The conventional Kelvin-Helmholtz cooling time is of order . For the monatomic uniform-density stellar model, its precise accessible total-energy reservoir is , so
Using , and , the usual order-of-magnitude normalization without the structural factor is . The full gravitational binding reservoir is about , but that overcounts the radiatable energy in virial equilibrium.
The Sun and the Solar system are approximately old, much older than any of these contraction estimates. Gravitational contraction cannot sustain the Sun's long-lived present power. The dominant long-term source is stellar nuclear fusion. Contraction can still supply transient luminosity, especially in a pre-main-sequence star. The timescale comparison rules out a contraction-only explanation over the observed age; it does not assert that all gravitational energy release is absent.
For a quasi-static self-gravitating object of fixed mass, negligible surface pressure, and no accretion, nuclear source or additional deep heating, conservation of energy gives the intrinsic planetary luminosity
Here is internal energy, and is negative gravitational potential energy. Gravitational contraction makes more negative and releases power . Some of that release raises the internal energy; it cannot all be radiated while the object maintains hydrostatic equilibrium.
To quantify the split assume thermal gas pressure with constant specific-heat ratio and negligible rotation or magnetic support. This excludes the dynamically unstable gas-supported regime. The stellar virial theorem gives
It follows that
For a monatomic ideal gas, , so : half the released gravitational energy heats the gas and half escapes. If the density profile remains homologous, write , with a fixed structure constant. Then
up to the structure and virial coefficients, where is the Kelvin-Helmholtz cooling time.
For a diatomic gas with the idealized fraction is , already demonstrating that the half-factor is not universal. Partial electron degeneracy pressure, dissociation, changing structure, accretion and deep heating require the full energy equation. The general link between luminosity and Kelvin-Helmholtz contraction is , with any genuine source terms added explicitly.
Let be Jupiter's current total emitted luminosity. Under the stipulated split, its intrinsic luminosity after removal of sunlight is approximately . If the accessible energy reservoir down to is , the constant-luminosity estimate gives .
Assume the more massive giant has approximately the same radius as Jupiter, a comparable structure factor, and an accessible cooling reservoir scaling as , as in the Kelvin-Helmholtz cooling time. Assume also that the final cold state contributes negligibly to that reservoir and that the giant radiates at the constant stipulated luminosity . Then
Therefore
This is a characteristic scaling under the stated assumptions, not a detailed evolutionary age. A different mass-radius relation, a different accessible internal-energy fraction, or time-dependent luminosity changes the result; the mass and luminosity alone do not uniquely determine a cooling time.
A cataclysmic variable is a close interacting binary star with a white dwarf accretor and a low-mass donor star, usually a red dwarf, filling its Roche lobe. Matter crosses by Roche-lobe overflow and commonly forms an accretion disk; sufficiently strong white-dwarf magnetism instead permits magnetically channelled accretion. Accretion produces strong variability, while classical novae result from unstable burning of accumulated hydrogen.
The cataclysmic-variable orbital-period distribution contains a short-period population between about 82 minutes and two hours, a pronounced cataclysmic-variable period gap around two to three hours, and a longer-period population above the gap. There is an accumulation near the cataclysmic-variable period minimum. These describe ordinary hydrogen-rich systems; evolved donors and helium-rich systems can have shorter periods. Selection by luminosity and outburst activity affects the observed counts, so is not simply the intrinsic distribution.
For the contact red dwarf, . Using its supplied thermal-equilibrium mass-radius relation, , gives . Kepler's third law then gives
independent of . This is a special case of the Roche-lobe-filling period-density relation.
To examine rapid mass loss, use the stellar radius response exponent . During conservative binary mass transfer, and are fixed, and . Therefore
The change in overfill is . Since , overfill grows, and mass loss runs away, when . The dynamical stability of binary mass transfer criterion here is consequently
The printed instability inequality is reversed. With the defined , the specified radius laws imply the inequality above; is the stable side. Equality is marginal, and the supplied Roche-radius approximation restricts this calculation to its stated mass-ratio range.
Long-term contact is maintained by angular momentum transport out of the orbit. Gravitational-wave emission from a binary system provides one loss mechanism. Magnetic braking of a binary star uses the donor's stellar wind: its magnetic field enforces approximate corotation out to the Alfvén radius, so the wind carries much more specific angular momentum than a nonmagnetic surface outflow. At that radius, roughly
The geometry contributes order-one factors to the torque. Even when , the large lever arm can remove substantial spin angular momentum. Tidal locking couples donor spin to orbital motion, so the wind torque ultimately brakes the orbit. In practice classical novae can also eject matter; the conservative calculations here isolate the stated idealization.
For the prescribed , now use the supplied equilibrium law for the slowly evolving donor. Contact implies , assuming negligible net wind mass loss. Taking the logarithmic derivative of gives
Since along this equilibrium sequence,
For the dynamically stable mass ratios, this is negative. The result uses thermal equilibrium as well as slow hydrostatic evolution; evolution slow compared only with the dynamical time does not by itself guarantee the radius law .
In the disrupted magnetic braking model, the donor above the gap loses mass fast enough to be inflated relative to thermal equilibrium. When the donor becomes a fully convective star, the model assumes a substantial reduction in magnetic braking. The red dwarf then contracts toward equilibrium on its Kelvin-Helmholtz cooling time, becoming smaller than its Roche lobe; accretion largely stops. Gravitational-wave emission from a binary system continues to reduce the separation, carrying the detached binary through the gap until contact is restored near its lower edge.
This explanation requires strong braking immediately above the gap: the donor must already be out of equilibrium and inflated, so there is room to contract and detach when braking drops. With nearly fixed component masses during detachment, , giving the illustrative inflation factor
A weak braking torque that kept the donor in thermal equilibrium just above the gap would not produce this detachment. The model requires a decrease in the torque, not disappearance of the donor's magnetic field.
A free-fall timescale estimate follows from :
Falling to a central point mass gives the more precise factor , or ; a distributed star changes an order-one coefficient. For any spherical mass distribution with radius ,
The inequality is strict for an extended star; equality is approached by a thin shell at . For an monatomic gas the virial theorem is , so half the gravitational potential energy released by contraction heats the star and half can be radiated. The Kelvin-Helmholtz cooling time is therefore of order
A uniform-density model gives . The given bound fixes a scale, rather than an exact contraction history; both estimates are far shorter than the Sun's main-sequence lifetime.
A free-fall timescale estimate follows from :
Falling to a central point mass gives the more precise factor , or ; a distributed star changes an order-one coefficient. For any spherical mass distribution with radius ,
The inequality is strict for an extended star; equality is approached by a thin shell at . For an monatomic gas the virial theorem is , so half the gravitational potential energy released by contraction heats the star and half can be radiated. The Kelvin-Helmholtz cooling time is therefore of order
A uniform-density model gives . The given bound fixes a scale, rather than an exact contraction history; both estimates are far shorter than the Sun's main-sequence lifetime.
Build the star by adding spherical shells. With the enclosed mass, its gravitational binding energy is
For a nondegenerate extended star, on the interior, so
This is a gravitational binding-energy lower bound for a spherical star. The Kelvin-Helmholtz cooling time is of order
The displayed energy bound gives years; it supplies a lower bound on that reservoir-to-luminosity ratio. A stellar structure model sets its coefficient. For an ideal monatomic, gas-supported star the virial theorem puts about half the released gravitational energy into heat, leaving half available for radiation; a uniform-density model would give about years. These contraction estimates are all of order years, much shorter than the nuclear lifetime.
Write with a structure factor of order unity. Then the thermal radiation timescale is
For a uniform sphere , this is about years. The conventional Kelvin-Helmholtz cooling time is about years; a centrally concentrated solar model changes the coefficient. Both estimates are far shorter than the nuclear lifetime.
Ionization and dissociation can store compression energy without the full temperature rise of a monatomic gas, lowering the effective adiabatic response in partially ionized outer layers. Radiation pressure tends towards an adiabatic exponent . These effects, changing composition and degeneracy make the effective exponent position dependent. The local dynamical exponent is ; it equals in the simple fixed-composition ideal-gas model, but need not do so for an ionizing mixture.
If an envelope has an effective exponent at or below , its pressure response is soft and sufficiently extensive such regions can promote radial pulsation, envelope instability or mass loss. In the ideal homologous global model, is marginal and is unstable to contraction or expansion. A low value confined to a thin outer layer does not by itself prove collapse of the entire star: global stability depends on the pressure-weighted structure. Its location and mass extent matter.
Stellar nuclear timescale 2026-10-14
The time for stellar nuclear fusion to consume the accessible fuel at characteristic luminosity . If a fraction of a star's mass supplies fuel abundance and released energy per unit fuel mass , hydrogen burning gives . Accessible core fuel, rather than the entire star, sets the lifetime unless mixing is complete. This is longer than a structural Kelvin-Helmholtz cooling time when the nuclear reservoir greatly exceeds the gravitational contraction reservoir.