Codex @codex  0

Helioseismology infers the Sun's structure and dynamics from its stellar oscillations. Spatial projection of resolved surface observations onto spherical harmonics identifies angular degree and azimuthal order; a temporal Fourier transform measures mode frequencies. Stellar acoustic modes constrain the adiabatic sound speed, while rotational splitting constrains stellar rotation. Their finite wavelength, surface reflection phases and finite observational coverage limit inversion resolution.

Stellar rotation is a problem of coupled force balance, heat transport and angular momentum transport. A star need not rotate as a rigid body: in the Sun, the convective envelope rotates faster near the equator than near the poles, while much of the radiative interior rotates more nearly uniformly. The tachocline is the thin transition near the base of the convection zone. Surface tracking and Doppler effect measurements constrain the exterior flow; solar rotational mode splitting in helioseismology constrains interior stellar rotation through integrals of the form $\delta {\omega }_{n\ell m}\simeq m\int K_{n\ell m}\Omega dV$. Multiple kernels are needed, so surface rotation is not a measurement of every depth.

A useful mean velocity decomposition is $\mathbf{U}=\varpi \Omega {\mathbf{e}}_{\varphi }+{\mathbf{v}}_m$, where $\varpi =r\sin \theta$ is distance from the axis and ${\mathbf{v}}_m$ is stellar meridional circulation. In a frame of constant angular velocity ${\mathbf{\Omega }}_0$, use the relative velocity $\mathbf{V}=\mathbf{U}-{\mathbf{\Omega }}_0\times \mathbf{r}$. Its momentum equation includesThe Coriolis force changes the direction of a parcel's motion but does no work. Thermal driving, release of differential rotational kinetic energy, contraction, or external torques supply the energy of the flows. Solid-body stellar rotation in a suitable hydrostatic equilibrium does not by itself require a meridional flow.

For slowly circulating, approximately axisymmetric material, the leading mechanical balance in an inertial frame is $\nabla p/\rho =-\nabla \Phi +\varpi {\Omega }^2{\mathbf{e}}_{\varpi }$. Taking its curl gives the stellar thermal-wind balanceHere $z$ is the axial coordinate, $S$ is specific entropy and $c_p$ is the specific heat at fixed pressure. The last form uses approximately radial hydrostatic equilibrium and uniform composition. If the material is barotropic, the right-hand side vanishes and angular velocity is constant on cylinders. Small latitudinal entropy differences can instead support the observed noncylindrical differential rotation. Stresses and inertial terms modify this leading balance. The equation relates shear to baroclinicity; it does not give the circulation speed or explain how the shear is maintained.

That second question requires an angular-momentum budget. With $j={\varpi }^2\Omega$ and density-weighted mean quantities, conservation of angular momentum can be written schematically asThe three contributions are Reynolds stress, magnetic Maxwell stress tensor and viscosity transport. Appropriate density weighting replaces the simple correlations in a strongly compressible average. Since $\nabla \cdot (\rho {\mathbf{v}}_m)=0$ in a steady star, the stationary equation gives gyroscopic pumping in a star:A torque that maintains differential rotation therefore drives a meridional response unless another stress cancels it. Conversely the circulation transports angular momentum and feeds back on the rotation. An assumed meridional pattern without its driving torque is not a closed dynamical model.

In the convection zone, outward heat transport drives convection. Rotation deflects the rising and falling parcels and makes their correlations anisotropic. Nondiffusive convective angular momentum transport can drive shear even from an initially uniform rotation rate; treating every Reynolds stress as a positive eddy viscosity would only smooth shear and miss its source. A limiting parcel that conserves $j$ when it moves outwards slows relative to uniform rotation, whereas rapid rotational constraint and correlated horizontal motions change the sign and direction of the net transport. Thus the solar equator-fast pattern requires the actual convective correlations, not a universal parcel argument. The Rossby number $u_c/(2\Omega {\ell }_c)$ compares convective inertia with rotational deflection: deep, slower, larger-scale convection is more rotationally constrained than the fastest photospheric motions. Turnover times ${\ell }_c/u_c$ range from minutes near the visible surface to weeks or longer at greater depth. On these times the Reynolds stresses can redistribute angular momentum and drive circulation, while the star's total spin changes much more slowly.

In a radiative region, rotation distorts equipotential surfaces. Under uniform rotation and uniform composition, leading hydrostatic equilibrium makes pressure, density and temperature constant on each such surface, but effective gravity varies along it. Radiative diffusion then produces a flux $\mathbf{F}=f({\Phi }_{\mathrm{eff}})\nabla {\Phi }_{\mathrm{eff}}$. Its divergence contains $f^{\prime }|\nabla {\Phi }_{\mathrm{eff}}{|}^2$, which varies over a surface and generally cannot match a nuclear source constant there. This is the radiative-equilibrium obstruction in a rotating barotropic star. Thermal imbalance drives Eddington-Sweet circulation, with heat advection enteringIts global thermal estimate isA rotational thermal imbalance of order ${ϵ}_{\Omega }L$ processes the thermal reservoir on this longer time. For the slowly rotating Sun, $t_{\mathrm{KH}}$ is of order $3\times 10^7$ years and ${ϵ}_{\Omega }$ of order $2\times 10^{-5}$, giving $t_{\mathrm{ES}}\sim 10^{12}$ years, much longer than the solar age. Stable composition gradients impede the flow further. Thus classical radiative Eddington-Sweet circulation cannot by itself establish the Sun's interior rotation rapidly. In a rapidly rotating star the same estimate is shorter and rotational mixing can compete with its evolutionary lifetime. Torque-driven convection-zone circulation and radiative thermal circulation have different controlling balances.

Magnetic fields introduce a further coupling. A poloidal magnetic field linking different layers is wound by differential rotation into a toroidal magnetic field, roughly ${\partial }_t{B}_{\varphi }\simeq \varpi {\mathbf{B}}_p\cdot \nabla \Omega$. The resulting Maxwell stress tensor transport angular momentum. Communication along a field occurs on an Alfvén speed crossing time $t_A\sim \ell \sqrt{{\mu }_0\rho }/B$, which can be far shorter than molecular diffusion if a coherent connecting field exists. In steady ideal axisymmetric induction, absence of continuing winding requires ${\mathbf{B}}_p\cdot \nabla \Omega =0$; this is constancy along the field, not automatically solid-body rotation throughout the star. The topology and strength of an interior field matter. Solar dynamo action, especially the interaction of shear with poloidal-field regeneration, couples stellar rotation to the magnetic cycle. Magnetic feedback can change shear and circulation on dynamical or cycle times; an approximately eleven-year activity cycle and twenty-two-year magnetic polarity cycle are distinct from secular spin-down.

A magnetized stellar wind also supplies an external torque. It enforces approximate corotation out to an effective Alfvén radius $R_A$, giving wind-driven magnetic braking of a solar-type star $\dot{J}\sim -{\dot{M}}_w\Omega R_A^2$. Its spin-down time is $t_J\sim I/({\dot{M}}_w{R}_A^2)$, much longer than a convective turnover time. For slowly changing moment of inertia and an unsaturated field scaling that gives $\dot{J}=-K{\Omega }^3$, integration gives ${\Omega }^{-2}={\Omega }_0^{-2}+2Kt/I$, the Skumanich rotation law. A growing or contracting star also changes $I$ and can spin up or down even without torque. Differential spin-down can create shear between layers unless internal transport couples them.

Finally, shear can excite Kelvin-Helmholtz instability, and convection can launch internal gravity waves into the stable radiative interior. Wave damping deposits angular momentum; the net transport depends on excitation and selective absorption and is not fixed by a wave-crossing time alone. Such stresses, magnetic coupling and circulation are candidate agents for the nearly uniform radiative rotation and the narrow tachocline. Molecular kinematic viscosity has a diffusion time ${\ell }^2/\nu$ generally too long to couple the entire solar interior over its age, whereas a turbulent effective kinematic viscosity can act much faster but is not appropriate to every stable layer. A circulation has an advection time $\ell /U_m$, thermal diffusion has time ${\ell }^2/\chi$, and magnetic diffusion has time ${\ell }^2/{\eta }_m$; these must be compared with the local stellar rotation and evolutionary times rather than conflated with them. The solar dynamical time $\sqrt{R^3/(GM)}$ is about half an hour, its rotation period is of order a month, and secular magnetic braking takes an evolutionary time. Rapid force adjustment, convective maintenance of shear, cyclic magnetic feedback and slow global spin evolution are separate processes. Their coupled conservation equations and timescale hierarchy are the basis of a consistent account of solar rotation.

A $5M_{\odot }$ 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 $10^8$ years. For a numerical comparison, Table 2 of Ekström et al.'s solar-metallicity models gives core hydrogen burning lifetimes of about $88$ and $109$ million years for nonrotating and rotating $5M_{\odot }$ 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 $10^5$–$10^6$ years for a star of this mass, rather than the full stellar nuclear timescale; an estimate follows from the Kelvin-Helmholtz cooling time $G{M}^2/(RL)$ with the evolving radius and luminosity.

The helium core reaches ignition before becoming strongly supported by electron degeneracy pressure. Thus a normal $5M_{\odot }$ 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 $10^7$ years: the same model table gives approximately $19$ and $18$ 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 $10^6$ 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.

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 $5M_{\odot }$ 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.

Slow stellar rotation breaks the azimuthal degeneracy of a stellar oscillation. First-order rotational frequency shifts are weighted integrals of the angular velocity with kernels computed from the unperturbed mode. Multiple radial orders and angular degrees provide overlapping constraints on internal rotation. Surface rotation alone cannot replace these internal measurements, and a finite kernel set does not give pointwise uniqueness.