Stellar astrophysics

Stellar astrophysics studies the structure, formation, energy production, observable properties, and evolution of stars.

Classification of stellar spectra by features controlled mainly by temperature and also by surface gravity and composition. The O-B-A-F-G-K-M temperature sequence does not mean a corresponding sequence in bulk hydrogen abundance.

A stellar spectral category based on gravity-sensitive line widths and ratios, distinguishing supergiants, giants and dwarfs at comparable spectral temperature. Examples are I for supergiants, III for giants and V for main-sequence dwarfs. It is inferred spectroscopically and is not a direct substitute for measuring absolute luminosity in every individual star.

The stellar temperature sequence O, B, A, F, G, K, M, with numerical subclasses approximately from hotter 0 to cooler 9 within each letter class. Spectral features follow temperature-dependent excitation, ionization and molecule survival, not a corresponding change in the bulk hydrogen abundance. A separate luminosity class describes gravity-sensitive spectral differences.

A plot of stellar magnitude against a colour index. It is an observational counterpart of the Hertzsprung-Russell diagram, but converting colour to effective temperature and magnitude to luminosity requires calibration, extinction and distance information. A coeval population's main-sequence turnoff helps estimate its stellar age.

An RR Lyrae star is a short-period pulsating evolved star used as a distance tracer in old populations. Its calibrated luminosity depends on composition and photometric band; infrared period-luminosity-metallicity relations are especially useful. extinction, crowding, mode identification and calibration remain important.

Time-dependent stellar phenomena such as starspots and chromospheric emission that can change brightness and spectral-line shapes, complicating radial-velocity planet searches.

Stellar variability that contaminates a planetary radial-velocity or photometric signal. Rotating active regions, changes in convective line shifts, oscillations and granulation can produce apparent velocities. A common spectral shift, line-shape diagnostics, activity indicators and simultaneous photometry provide different tests; none is an automatic numerical probability.

The midpoint of the two wings of a stellar absorption line at each depth. A changing line bisector can reveal profile distortions caused by active regions or blends rather than a rigid Doppler shift of the whole spectrum.

Neutrino astronomy infers properties of astronomical sources from detected neutrino reactions. Solar neutrino detection probes ongoing stellar nuclear fusion; the SN 1987A neutrino burst probes a brief episode of core-collapse supernova evolution. A detector observes reaction products, so source inference depends on cross sections, neutrino flavor, thresholds and detection efficiency.

Kamiokande, associated with Masatoshi Koshiba, was an underground water Cherenkov detector viewed by photomultipliers. Improved water purity and reduced background enabled solar neutrino observations. A concentration of reconstructed recoil directions toward the Sun distinguished a source signal above an approximately direction-independent background.^[ref]

For solar neutrinos, electron elastic scattering supplies a directional recoil, with different cross sections for active neutrino flavors. Neither the Homestake experiment deficit nor the SN 1987A neutrino burst alone established neutrino oscillation; later evidence compared complementary channels and propagation effects.^[ref]

Water Cherenkov neutrino reconstruction infers a vertex, time and charged particle direction from photomultiplier pulses. The 1987 analysis recorded charge and timing, calibrated with decaying muons and source-produced gamma rays, and rejected outer-detector activity and events following muons.^[ref]

In an ideal point-flash approximation, a light source at ${\mathbf{r}}_0$ and time $t_0$ gives the displayed travel-time relation at sensor ${\mathbf{r}}_i$ in a medium of refractive index $n$. Subtracting two times eliminates $t_0$. With a chosen constant $n=1.33$ and path lengths $3\mathrm{m}$ and $6\mathrm{m}$, the arrival-time difference is $n(3\mathrm{m})/c\simeq 13.31\mathrm{ns}$. These are synthetic inputs, not detector measurements. Several sensors constrain the vertex; an extended moving track requires the light-emission position and time along that track as well.

For a charged particle moving at $v=\beta c$, the Cherenkov cone angle satisfies $\cos {\theta }_C=1/(n\beta )$ when $n\beta >1$. With $\beta \simeq 1$ and the illustrative $n=1.33$, ${\theta }_C\simeq 41.25^{\circ }$. The geometrical cone carries track-direction information; electron deflections and the spread of emission positions broaden the observed pattern. The neutrino itself does not draw the light cone.

A Cherenkov radiation threshold and an analysis threshold are different constraints. For a fixed $n$, an electron must exceed $T_C=m_e{c}^2[(1-n^{-2}{)}^{-1/2}-1]$. With $m_e{c}^2=0.511\mathrm{MeV}$ and the chosen $n=1.33$, this gives about $0.264\mathrm{MeV}$. Recording and selecting a useful event requires much more light: the 1987 detector's trigger accepted $8.5\mathrm{MeV}$ electrons with 50% efficiency.^[ref]

Orbital dynamics gives the most direct stellar masses: a resolved binary orbit and distance determine the relative semimajor axis, and Kepler's third law gives the total mass. The barycentric size or radial-velocity amplitude ratio gives individual masses. An eclipsing binary inclination removes the sine-cubed projection factor from spectroscopic masses; a single-lined binary mass function alone does not uniquely determine both masses. Surface-gravity/radius and stellar-track estimates are model-dependent and are calibrated against dynamical binary measurements.

Loss of matter from a star or other astronomical system, for example through a stellar wind or explosive ejection. Mass loss rate conventions must distinguish the positive ejected rate $|\dot{M}|$ from the negative change of the retained mass, $\dot{M}$. In a binary star, the fate of an orbit also depends on the specific angular momentum carried away.

The spin of a star about its axis, which can be solid-body rotation or vary through the interior. Angular momentum transport and rotational mixing can change the stellar composition profile and available fuel. An observed surface rotation rate does not determine the entire interior profile.

In a steady axisymmetric star, a torque due to Reynolds stress, Maxwell stress tensor or viscosity must balance advection of specific angular momentum by stellar meridional circulation. Writing the nonadvective angular momentum transport flux as ${\mathbf{F}}_J$ gives the displayed balance. Thus stresses that maintain differential rotation also drive a meridional response. The stellar thermal-wind balance supplies a separate leading force constraint, not a replacement for this torque equation.

Taking the curl of axisymmetric hydrostatic equilibrium with centrifugal acceleration gives $\varpi {\partial }_z{\Omega }^2=(\nabla p\times \nabla \rho {)}_{\varphi }/{\rho }^2$. For an approximately spherical ideal gas with constant composition, use $\nabla \rho =\nabla p/c_s^2-\rho \nabla S/c_p$ to obtain the displayed relation. Latitudinal entropy differences can support noncylindrical differential rotation. A barotropic fluid gives cylindrical rotation in this leading balance; the relation does not determine the transport stresses or circulation rate.

Outside a spherical star, its gravitational acceleration is $-G{M}_{\star }\mathbf{r}/r^3$. It provides the dominant central force for Kepler orbits in a star-dominated planetary system.

A stellar population is an ensemble of stars described by its formation times, initial masses, composition and present stellar evolution states. Number fractions must be distinguished from mass-weighted fractions. Shared formation times can correlate the states of components of a binary star even when their initial masses have independence.

The event rate per unit initial stellar mass and per unit delay after star formation. An extended star formation history produces a rate $R(t)=\int \mathrm{SFR}(t-\tau )\Psi (\tau )d\tau$. Measuring rates versus age in old populations constrains progenitor channels only after accounting for stellar lifetimes, age spreads and selection.

A Population III star is a first-generation star formed from essentially metal-free primordial gas. Limited cooling generally raises its characteristic fragmentation mass and favors a more top-heavy initial mass function than in enriched star formation.

A Population II star is an old metal-poor star formed from gas enriched by earlier stellar generations. Metal and dust cooling let enriched gas fragment to lower masses than primordial gas under otherwise comparable conditions.

Integrated spectra can confuse older stellar populations with more metal-rich ones because both can redden the light and change absorption lines. Line strengths depend on stellar metallicity, temperatures, ages, and abundance ratios. Weak metal lines in a young population dominated by hot bright stars need not indicate low stellar metallicity.

Given a birth-mass distribution and age distribution, an evolutionary state occupies the mass-age set where its starting lifetime has elapsed and its ending lifetime has not. The state number fraction is the probability measure of this set. Initial mass function tail coordinates and the probability integral transform can turn these weighted regions into simple geometric areas.

A phase bounded by mass-dependent lifetimes ${\tau }_a(M)$ and ${\tau }_b(M)$ occupies ${\tau }_a(M)<t<{\tau }_b(M)$, clipped to the population age range. Its boundaries can be mapped into uniform probability coordinates to evaluate evolutionary-state fractions.

If a population forms at a constant number rate for duration $T$ and every birth object is retained in the census, its present ages are uniform on $[0,T]$. A fixed initial mass function with finite mean birth mass also makes constant mass-based star formation rate equivalent to constant number rate.

A stellar atmosphere is the outer region that determines the escaping spectrum and matches the interior stellar structure to its radiation boundary conditions. A photosphere is the layer near unit optical depth, not a zero-temperature boundary. The inward intensity at the outer surface vanishes when there is no external illumination.

The distribution of radiation from a star with wavelength or frequency, including its continuum and absorption lines. Continuum colors and temperature-dependent line strengths provide complementary estimates of effective temperature. Spectra of an unresolved binary star may need to be separated by eclipse observations or their different velocity shifts before the component effective temperatures can be inferred.

A stellar abundance can be expressed as an elemental or isotopic mass fraction, or as a number ratio to hydrogen. Inferring it from a spectral line requires a model for level populations, opacity, radiation and velocity fields. Agreement between suitable ionization stages and weak photospheric lines tests whether the adopted stellar atmosphere is adequate.

The abundance of chemical elements heavier than helium in astronomical material. A metal mass fraction measures their combined mass; spectroscopic abundance ratios are other measures of metallicity and need not be numerically identical to that fraction.

The fraction of total mass in chemical elements heavier than helium, so the hydrogen mass fraction, helium mass fraction and metal mass fraction satisfy $X+Y+Z=1$. Despite the astronomical name, this includes elements such as carbon and oxygen that are not chemical metals.

When radiative rates compete with or dominate collisions, level populations need statistical equilibrium coupled to radiative transfer, rather than thermal population laws at the local temperature. Radiation from other depths can make the local radiation field differ from the Planck function. This changes ionization and spectral-line strengths, so a geometrically thin atmosphere may still require non-LTE abundance modelling.

A heated layer above a photosphere can have an outward temperature increase. Strong lines formed there may show emission where an outward-cooling photosphere would produce absorption. The energy maintaining the temperature inversion is additional atmospheric heating; the exact height and transfer regime vary between stars.

A narrow spectral feature records enhanced emission or absorption near a transition frequency. In a stellar atmosphere, its strength depends on atomic populations, opacity, velocity fields and the source function. Lines therefore diagnose composition as well as the physical atmosphere in which they form.

Two nearby resolved spectral components. A two-line pattern can reflect fine-structure splitting, but counting lines is distinct from assigning a term's spin multiplicity.

The finite width of a spectral line can arise from finite optical coherence or a distribution of resonance frequencies. lifetime broadening and impact-regime pressure broadening often give Lorentzian line shapes; thermal Doppler broadening gives a Gaussian component. The instrument response must also be included when interpreting a measured spectrum.

Broadening for which all participating emitters share the same spectral response, as in a uniform lifetime-limited line. A resonant field can saturate the common transition response. This differs from a distribution of transition frequencies across independently resonant subensembles.

Finite interaction time with a spatially confined optical beam gives a finite spectral response width. Its detailed shape depends on the spatial beam profile, particle trajectory and interaction envelope; a duration alone does not imply a universal Lorentzian profile.

An ensemble-wide distribution of resonance frequencies, such as velocity-dependent Doppler shifts or spatially varying static-field shifts. The observed profile averages the individual responses. A velocity class need not have the full width of the ensemble's thermal Doppler broadening.

Broadening shared by otherwise identical emitters, often described by a common optical coherence time $T_2$. Exponential field-coherence decay produces a Lorentzian line shape with full angular-frequency width $2/T_2$. Population decay, dephasing and collisions can all contribute, with their time constants distinguished.

In the absence of additional dephasing, a transition whose excited-state population has lifetime $T_1$ has a field-coherence decay time $2T_1$. Its Lorentzian line shape has angular-frequency half-width $1/(2T_1)$ and full width $1/T_1$. Population decay time, amplitude decay time and linewidth conventions must therefore be distinguished. Extra dephasing shortens the optical coherence time without necessarily changing the population lifetime.

The locus of midpoints between the two sides of a spectral-line profile at equal intensities. Variations with depth or time diagnose changes in line asymmetry; a null result has finite sensitivity.

A pair of neighbouring spectral lines.

Collisions and fluctuating neighboring-particle electric fields alter atomic energy levels and coherence times, broadening a spectral line. In a simple impact treatment the coherence-decay rate gives a Lorentzian width of order $\Gamma /(2\pi )$; the precise definition of full or half width must be stated. Unlike Doppler broadening, this depends primarily on collision rates and perturber density. In stellar opacity calculations it alters line overlap and the low-opacity windows entering the Rosseland mean opacity.

If phase-randomizing collisions form a Poisson process of mean interval ${\tau }_c$, field correlation survives only when no such collision occurs during the delay. Thus its magnitude decays as $e^{-|\tau |/{\tau }_c}$. The Wiener-Khinchin theorem gives the normalized Lorentzian line shapeIts full width at half maximum is $2/{\tau }_c$ in angular frequency or $1/(\pi {\tau }_c)$ in ordinary frequency. A population lifetime is a different time constant: lifetime-limited amplitude decay at half the population rate introduces the familiar factor of two. The definition of the supplied lifetime must therefore accompany any numerical linewidth.

A spectral line profile with the shape of a Cauchy distribution. In the displayed peak-normalized convention, $\nu$ is frequency, ${\nu }_0$ is the line centre, and $\gamma >0$ is the half width at half maximum: $L({\nu }_0\pm \gamma )=1/2$, so the full width is $2\gamma$. The normalization used for a probability density function differs from this unit-peak convention.

A Voigt profile is the convolution of a Gaussian distribution and a Cauchy distribution. For an absorption line, thermal Doppler broadening contributes the Gaussian core while finite excited-state lifetime contributes the Lorentzian wings. Instrumental broadening and unresolved motions must also be modeled when inferring column density and temperature.

An emission line is a localized excess of specific intensity above the continuum. Heated chromospheric material can give a line source above the photospheric continuum source. Optically thin circumstellar gas can also add line photons through recombination or collisional excitation. Such line formation need not obey a thermal photospheric source prescription.

A forbidden emission line is suppressed by the electric-dipole selection rules but can occur through weaker radiative channels such as magnetic-dipole or electric-quadrupole emission. Its relatively long-lived upper atomic level is vulnerable to collisional de-excitation in dense gas. Low-density astrophysical plasmas can therefore display strong forbidden-line emission even when laboratory dense sources do not.

In a three-level ion, statistical equilibrium must include collisional coupling between both upper levels as well as excitations from the ground. The upper populations solve a two-by-two linear system when expressed relative to the ground population. Their line ratio is a rational function of electron density, interpolating between a collisionally excited low-density limit and a collisional Boltzmann limit. Closely spaced upper energies reduce temperature sensitivity, but neglecting upper-level collisional mixing can bias the inferred density.

For low-density gas with Maxwellian electrons, the excitation-energy difference between the ${}^1{S}_0$ and ${}^1{D}_2$ levels of doubly ionized oxygen makes $(I_{4959}+I_{5007})/I_{4363}$ temperature sensitive. A constant-collision-strength approximation is $R\simeq 7.90e^{32900\mathrm{K}/T_e}$. Increasing temperature strengthens the higher-excitation auroral line and lowers $R$. Density, extinction and nonthermal electrons require separate treatment; the $5007/4959$ pair shares one upper level and is primarily a branching ratio rather than this thermometer.

An absorption line is a deficit of observed spectral specific intensity relative to the surrounding continuum. In a photosphere with temperature falling outward, increased line opacity makes photons emerge from cooler layers than the continuum. For an absorption-dominated LTE source, the Planck function then supplies a smaller line specific intensity. A foreground layer with source below the incident continuum gives another direct transfer description.

A highly luminous evolved star, associated with luminosity class I. Its location in a Hertzsprung-Russell diagram is distinct from the low-luminosity, compact white dwarf region.

A cool, highly luminous evolved massive star with an extended envelope, generally later than its initial hot main sequence stage. Its eventual core-collapse outcome depends on the core, winds, mixing and binary history, not on its red colour alone.

An evolved star with substantially larger stellar radius and luminosity than a main sequence star of similar effective temperature. A red giant is a cool example; giant does not mean merely any star of large mass.

A star with broadly Sun-like properties; the range depends on the stated classification criterion. In elementary planetary estimates, a solar mass, solar radius, and solar effective temperature are common approximate inputs, rather than consequences of the name alone.

The physical radius assigned to a star, with a specified surface convention.

A low-mass star maintained by sustained hydrogen burning, rather than the continuing cooling of a brown dwarf. Near the hydrogen-burning minimum mass its radius is of order a Jovian radius. Stars below roughly $0.35M_{\odot }$ are approximately fully convective, with a radiative photosphere; mass, composition and age affect the precise transition.

The solar radius is the conventional unit of length based on the radius of the Sun. The nominal value is $R_{\odot }=6.957\times 10^8\mathrm{m}$.

The solar corona is the Sun's hot, dilute outer atmosphere. Its optically thin atomic lines diagnose temperature, density and composition through collisional-radiative models and ionization balance. Its plasma extends into the solar wind; the coronal approximation describes useful low-density equilibrium limits rather than every coronal condition.

A larger-scale horizontal cellular flow pattern on the Sun's surface. Transport by these flows helps arrange small magnetic elements into a solar magnetic network at cell boundaries. The observed pattern does not establish a unique thermal driving depth or an exact hierarchy of isolated steady eddies.

A solar granule is a photospheric convection cell, with relatively hot rising material in its interior and cooler descending material near its edges. The pattern formed by their bright interiors and darker edges is solar granulation. Their evolving horizontal outflows can sweep weak magnetic flux toward surrounding lanes.

The evolving cellular intensity pattern on the Sun's photosphere, produced by solar granules. Its radiative appearance and flow pattern are related but are not identical observables. Near-surface horizontal flow can redistribute weak magnetic field into the surrounding intergranular lanes.

An intergranular lane is a region between solar granules, commonly associated with cooler material and descending flow. Horizontal convergence can concentrate approximately vertical magnetic flux there. The field's geometry, diffusivity and feedback determine the actual concentration rather than the lane label alone.

A solar neutrino is a neutrino produced by a nuclear reaction in the Sun, principally through the proton–proton chain and CNO cycle. Its weak interactions allow it to escape directly from the stellar core. Accounting for neutrino oscillation allows measurements of the different detected flavors to test the nuclear reaction power and rates independently of escaping photons.

Solar neutrino detection asks how an invisible electron neutrino from the Sun becomes a measurable laboratory record. The Homestake experiment stores a reaction product for later counting; Kamiokande records light promptly. The former measures an exposure-integrated rate, whereas the latter can retain event times and statistical source-direction information. These are complementary ways to connect stellar nuclear fusion with neutrino astronomy.

A detector samples a weighted part of a neutrino flux, rather than counting every passing neutrino. For a target population $N_T$, an ideal reaction rate iswhere ${\varphi }_e(E)$ is the differential electron neutrino neutrino flux and ${\sigma }_e(E)$ is the relevant cross section. A detection efficiency belongs inside the integral if it depends on energy. Two detectors can therefore disagree in raw count rate while observing the same source. Interpreting a deficit requires both a source prediction and a detector response.

The Homestake experiment, led by Raymond Davis Jr., used electron neutrino capture on chlorine to produce argon subject to radioactive decay. The 1968 paper reported an upper limit, rather than a resolved positive signal.^[ref]

The underground target contained tetrachloroethylene. Circulated helium swept out argon, which was trapped on cooled charcoal, purified, and put into tiny proportional counters. Electron capture back to chlorine supplies the delayed readout. Added stable argon checked extraction recovery; background counting checked spurious pulses. The capture threshold is $0.814\mathrm{MeV}$ and the product's half-life is about 35 days. The chemical extraction retains neither the original neutrino's arrival time nor its direction.^[ref]

The reaction selects electron neutrinos, and its threshold excludes the lowest-energy part of the solar neutrino spectrum. A low count alone cannot distinguish a low source neutrino flux, a changed neutrino flavor mixture, or an underestimated loss in extraction or counting. This is why chemical recovery tests and an independent water Cherenkov neutrino reconstruction are scientifically valuable.

An idealized radiochemical neutrino counting model separates production, survival, recovery and readout. Suppose a freshly cleared target produces atoms at a constant rate $R$ during exposure $T$. Each atom has nuclear decay constant $\lambda =\log 2/h$, where $h$ is its half-life. The expected surviving population obeysThe proof is to multiply the differential equation by $e^{\lambda t}$ and apply integration. The expression also follows by summing the exponential decay survival probabilities of atoms produced at every earlier time. Production does not imply eventual recovery: atoms can decay during exposure.

Let $\eta$ be the fraction recovered instantaneously at extraction, $d$ the subsequent delay before counting, $ϵ$ the probability of recording a decay, and $W$ the counting duration. The expected signal pulses areWith constant background rate $b$, the expected total is $\mu ={\mu }_s+bW$. During counting, the expected instantaneous pulse rate is $s(t)+b$, with $s(t)=ϵ\eta \lambda N(T)e^{-\lambda (d+t)}$. Equal time windows therefore lose signal by the factor $e^{-\lambda \Delta t}$ while a constant background stays unchanged. Fitting the time dependence helps distinguish radioactive decay from that background; it does not recover individual capture times.

For an original worked example choose $Rh=20$, $T=W=h$, $d=0$, $\eta =0.8$, $ϵ=0.5$ and $bh=0.2$. These are illustrative assumptions. Then $N(T)=10/\log 2$, andDoubling the exposure does not double the surviving inventory: $N(2h)/N(h)=3/2$. A delay of one half-life halves the signal. At $\lambda T\ll 1$, $N(T)\simeq RT$; at long exposure it saturates at $R/\lambda$. These limits explain why a large target and repeated extraction can matter more than simply waiting longer. The model neglects time-dependent production, finite extraction duration, residual atoms between runs and time-dependent backgrounds.