Astrophysical S factor 2026-10-18
The astrophysical S factor separates the dominant Coulomb barrier suppression from a charged-particle nuclear reaction cross-section :
Here is relative kinetic energy, is reduced mass, is the fine-structure constant, and are atomic numbers. The astrophysical S factor has units of energy times area. Factoring out a barrier does not remove all nuclear physics: nuclear levels and reaction mechanisms can still make vary strongly.
For a dilute nonrelativistic Maxwell-Boltzmann velocity distribution of relative motion, multiplying the relative-energy density by speed and nuclear reaction cross-section cancels the explicit factors of , giving
Thus the Gamow peak describes a rate-integral kernel, not a histogram of reacting stars or even the unweighted particle population. Plasma screening and departures from the assumed Maxwell-Boltzmann velocity distribution require additional physics.
A hotter stellar core can favor a higher-barrier route despite fewer catalyst nuclei. To isolate this competition, prescribe a local hydrogen burning toy model at fixed mass density , with reference and
Choose and , rounded barrier parameters for a proton pair and proton capture on nitrogen-14. These prescribe the temperature dependence, while the following freely chosen heating normalization incorporates effective nuclear factors and deposited energy:
Here is hydrogen mass fraction, is CNO mass fraction, and is an arbitrary positive unit of stellar energy-generation rate. For and , the reference ratio is by construction. It is not a measured solar ratio. The CNO cycle bottleneck is represented by one effective rate, without evolving the existing proton-proton chain branches or catalyst network.
At , the logarithmic temperature slopes are and . The larger Coulomb barrier gives a steeper fractional response. Because these slopes change with temperature, this example retains the exponential rather than extending two fixed power laws. Writing and gives
provided the bracket is positive. With the chosen composition, . Reducing to at unchanged shifts it to . For , this catalytic route supplies no heat in this model. These are local toy crossovers: a whole star's luminosity requires integrating its stellar energy-generation rate through its stellar composition profile and stellar structure. Nuclear features, plasma screening and changes of the dominant bottleneck limit the model's extrapolation.
Figure 1.
Chosen pp and CNO heating rates and their composition-dependent crossover
. Original local toy rates at fixed density and hydrogen mass fraction 0.70. The CNO-to-pp ratio at 15 MK is chosen as 0.01 for CNO mass fraction 0.01. Exponential barrier responses give crossovers at 20.30 and 23.91 MK for the two catalyst abundances. These are model curves, not stellar measurements or a fitted solar model. Details.
Debye length 2026-10-18
For mobile thermal Electrons with a fixed neutralizing background, the Debye length is the characteristic distance over which plasma screening suppresses an electrostatic perturbation. Linearizing the Boltzmann density response and using Poisson's equation gives the displayed scale. Other mobile species contribute additively to through .
A plasma contains mobile charged species whose collective electromagnetic response matters. On scales large compared with the Debye length, positive and negative charge densities nearly cancel: this is quasineutrality, not the absence of electric currents or electric fields. For electrons of temperature and number density , the characteristic plasma screening and oscillation scales are
A collective plasma description requires many particles within a screening volume; long-wavelength, sufficiently slow behavior can often be described by magnetohydrodynamics with an appropriate fluid closure.
An individual particle feels the Lorentz force . In a uniform magnetic field, its transverse motion gyrates while its parallel motion remains free, giving a helical path when both components are present. The nonrelativistic cyclotron frequency and gyroradius are
When the system scale is much larger than and the relevant time scale longer than , the field strongly organizes transport across it, provided collisions do not interrupt the gyromotion. A perpendicular electric field also gives the common guiding-centre drift .
For a conducting fluid, the macroscopic Lorentz force density has the useful decomposition
These terms are magnetic pressure and magnetic tension. Magnetic tension resists bending, while magnetic-pressure gradients resist nonuniform compression. The plasma beta compares thermal and magnetic stresses.
With uniform magnetic diffusivity , the resistive induction equation is
If the magnetic Reynolds number , ideal magnetohydrodynamics applies approximately and magnetic flux freezing links the field geometry to fluid motion. Finite diffusivity permits field slippage and magnetic reconnection. Field tension supports Alfvén waves with ; coupled magnetic and gas compression also supports magnetosonic waves. These effects allow a field to transmit forces and angular momentum between distant parts of an astrophysical plasma.
For a cold, unmagnetized, collisionless plasma, take an equilibrium Electron density , Electrons of mass and charge , and approximately stationary ions. Use fields proportional to . The Electron displacement satisfies
Its electric polarization is , hence the cold collisionless plasma dielectric response
The negative susceptibility arises from the inertial Electron response, rather than a restoring force binding an Electron to a particular atom.
A transverse electromagnetic wave then obeys
For , the plasma supports propagation, with refractive index . The phase velocity is greater than , while the group velocity is ; their product is , and no superluminal signal is implied. Below the plasma frequency, the bulk wavenumber is imaginary and the field is evanescent. A sufficiently thick plasma largely reflects such radiation in the lossless model.
Including an Electron momentum-relaxation rate replaces the denominator by :
Its imaginary part describes absorption with the stated time convention. Finite Electron temperature also produces plasma screening of slowly varying electrostatic fields, with . That spatial, longitudinal screening problem is distinct from simply taking in the cold transverse-wave formula. Magnetic bias, ion motion and kinetic resonances require more general dielectric responses.
Electromagnetic waves drive mobile charges in a conducting medium. Their current changes both the phase and attenuation of the wave. A description in terms of a real, frequency-independent refractive index is therefore generally insufficient: the material response and the propagation constant are complex and frequency dependent.
Use phasors proportional to , with and for a passive medium. In a homogeneous, isotropic, linear material, let represent the permittivity associated with bound-charge electric polarization, the magnetic permeability and the mobile-charge conductivity. The macroscopic Maxwell equations give
For a transverse plane wave, taking another curl yields
This division of the response avoids counting mobile Electrons both in and in . Define the complex refractive index . For a nonmagnetic medium, . Its real part determines the phase variation, while the positive optical extinction coefficient gives amplitude attenuation and intensity attenuation . The amplitude penetration length is
A simple microscopic explanation comes from the Drude model. Let Electrons of charge , mass and number density move against an effectively stationary positive background. Collisions relax their mean velocity at rate :
Consequently
With a vacuum-like bound-charge background, this gives
Its imaginary part is positive for positive , consistent with dissipation under the chosen phasor convention. Real materials can require an effective Electron mass, a nontrivial bound-charge response, or additional absorption associated with the electronic band structure.
A plasma is a conducting ionized medium with mobile Electrons and ions, approximately neutral on large scales because of plasma screening. For the electron-frequency transverse waves considered here, the heavy ions can often be treated as stationary. In a dilute weakly collisional plasma, , so
For , is real in this collisionless model and the wave propagates. Its phase velocity is , while its group velocity is ; their product is . A phase crest exceeding does not imply superluminal information transport. For , the refractive-index square is negative and
The transmitted field is evanescent rather than a traveling bulk wave. A sufficiently thick collisionless plasma reflects such radiation without Ohmic absorption. Collisions make the response complex and permit absorption. This accounts for a plasma cutoff and, qualitatively, the reflection of some radio frequencies by ionized atmospheric layers. The simple model describes transverse waves; plasmas also support longitudinal collective modes.
In a metal, positive ions usually form a crystal lattice. Mobile conduction Electrons occupy states described by the electronic band structure and normally form a degenerate Fermi gas. The conduction-electron density is usually much higher than in a dilute gas plasma. Thus its plasma frequency is commonly in the optical or ultraviolet range, rather than the lower radio-frequency values possible in dilute plasmas. Electron scattering from the lattice, defects and other excitations gives a finite relaxation rate. The distinction is not that every plasma is collisionless and every metal collisional: the controlling ratio for either material is , together with its carrier density and bound-charge response.
At sufficiently low frequency, , the metal's conductivity is approximately real, . In a good conductor the conduction current dominates the displacement current, . The wave equation then gives
For example, a tangential electric field varies as
The current is therefore concentrated within a few skin depths of the surface: this is the skin effect. As frequency increases within this regime, . In a wire whose radius greatly exceeds , only a thin outer region carries appreciable alternating current, increasing the effective resistance. In a conducting shield, the same penetration suppresses the field transmitted through a sufficiently thick wall.
The Poynting theorem identifies the loss as Joule heating. The cycle-averaged heating density for a complex conductivity is
The wave impedance is . For incidence from vacuum at normal incidence, the electric reflection amplitude is , where . A good conductor has small and therefore strong reflection, while the part of the incident power entering it is dissipated over the skin depth.
At high enough frequency that , a metal's free-electron response also approaches the collisionless plasma form. The low-frequency formula must not simply be extrapolated into that regime. Real optical behavior additionally depends on the electronic band structure and transitions between bands. Both plasmas and metals are therefore described by the same electromagnetic response framework, but their relevant densities, scattering rates and frequency regimes differ.
In a dilute thermal helium gas with rapid formation–decay balance of beryllium-8, write its small equilibrium abundance schematically as . A narrow isolated Hoyle state gives a second-step coefficient proportional to . Multiplication yields
with the stabilizing decay strength absorbed into the coefficient. Here is the Boltzmann constant. The two positive offsets count energy relative to their respective entrance thresholds, not relative to the carbon-12 ground state. A resonance enhances capture at the relevant energy; the Boltzmann factor still selects how many encounters reach it. This is a sequential narrow-level approximation with fixed nuclear widths, excluding plasma screening and changes of composition or mass density.
William Alfred Fowler lists the Hoyle state above the carbon-12 ground state, above beryllium-8 plus an alpha particle, and above three alpha particles.[ref] Subtraction gives for the diagram.
Normalize the calculated rate at and . Holding alpha particle number density and the decay strength fixed gives the original sensitivity calculation
At , a deliberately chosen shift multiplies the rate by ; gives . These offsets are hypothetical perturbations, not alternative measured nuclear levels. The reference curve has local temperature exponent of a thermonuclear reaction : raising temperature by gives a factor in this fixed-environment calculation. A real star changes its stellar structure and composition as it burns, so these factors do not directly predict final carbon abundance or a constraint on life.
Figure 1.
Triple-alpha energy thresholds and calculated resonance sensitivity
. Original energy bookkeeping using Fowler's rounded offsets (1983 lecture, p. 179). The normalized narrow-resonance rate calculation holds alpha density and decay strength fixed. The plus/minus 100 keV shifts are deliberately chosen sensitivity tests. No plotted curve is laboratory data. Bound carbon levels lie below the diagram's energy window; the downward arrow indicates stabilization rather than a measured lifetime or branching fraction. Details.