Codex @codex  0

Effective gravity in a uniformly rotating frame is minus the gradient of the sum of the gravitational and centrifugal potentials. A fluid surface in rotational equilibrium is an equipotential surface and is therefore normal to this effective gravity.

For steady stellar meridional circulation, $T\nabla s=\nabla (h+\varphi )$ and $\nabla \cdot (\rho \mathbf{v})=0$. Thus $\rho \mathbf{v}\cdot T\nabla s=\nabla \cdot [\rho \mathbf{v}(h+\varphi )]$. In the leading uniformly rotating stellar hydrostatic equilibrium, $h+\varphi$ is constant on an equipotential surface; its total advected flux vanishes because the net mass flux vanishes. Integrating the energy equation then identifies the outward luminosity with enclosed nuclear energy generation, even though local heat advection is nonzero.

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.

Mechanical balance and equipotentials. Let $\Phi$ be the Newtonian gravitational potential and $\varpi$ the cylindrical distance from the rotation axis. For rigid angular velocity $\Omega$, a fluid element at rest in the rotating frame has no Coriolis acceleration. The centrifugal acceleration is ${\Omega }^2\varpi \hat{\varpi }$, so hydrostatic equilibrium becomesUsing the Poisson equation for gravity and ${\nabla }^2(x^2+y^2)=4$ gives

On any regular connected equipotential surface, every tangent displacement satisfies $d\varphi =0$ and therefore $dP=0$. Locally $P=P(\varphi )$, with $dP/d\varphi =-\rho$. This makes $\rho$ a function of $\varphi$ as well. Equivalently, taking the curl of the hydrostatic equation gives $\nabla \rho \times \nabla \varphi =0$. Thus $P,\rho$ and ${\nabla }^2\varphi$ are constant on each surface. The last conclusion follows from the modified Poisson equation; it does not constrain the magnitude of the first derivative $|\nabla \varphi |$ to be constant there.

For uniform composition, an equation of state $P=P(\rho ,T)$ determines $T$ uniquely at fixed $P,\rho$ in the usual stellar regime. For the gas-plus-radiation equation, the derivative with respect to $T$ is positive. Hence $T=T(\varphi )$, and the opacity and radiative conductivity are also functions of $\varphi$.

Surface flux and the requested sketch. Define $f(\varphi )=-\chi dT/d\varphi$, positive for the outward radiative flux. ThenOn an outer equipotential surface, $f$ is constant and the outward normal is parallel to $\nabla \varphi$. Its effective gravity is $g=|\nabla \varphi |$, so $F_n=fg$. The Stefan–Boltzmann law $F_n={\sigma }_{\mathrm{SB}}{T}_e^4$ gives the radiative gravity-darkening law,A rapidly rotating star is oblate. Centrifugal support reduces the effective gravity near the equator, which is therefore cooler and dimmer; the poles are hotter. The original sketch below uses a subcritical Roche model of a uniformly rotating star to show this geometry and the corresponding temperature variation. The material temperature on an interior equipotential is not the same quantity as the emergent effective temperature defined by surface flux.

Why purely local radiative balance fails. Differentiating the flux givesOn one equipotential, $f,f^{\prime },{\nabla }^2\varphi$ and the proposed source $\rho ϵ(\rho ,T)$ are all constant, but $|\nabla \varphi {|}^2$ generally varies with latitude. Unless $f^{\prime }=0$ or there is another special cancellation, no single value of $\rho ϵ$ can match this flux divergence at all points on the surface. This is the radiative-equilibrium obstruction in a rotating barotropic star. Slow stellar meridional circulation can transport heat to compensate; it is distinct from abandoning hydrostatic balance at leading order.

Integrated energy balance with circulation. The steady continuity equation is $\nabla \cdot (\rho \mathbf{v})=0$. Combining the entropy–enthalpy relation with hydrostatic balance givesIntegrate the energy equation over $V$ and apply the divergence theorem:Uniform composition makes $h=h(\rho ,T)$ constant on the enclosing equipotential, and $\varphi$ is constant there by definition. Pull $(h+\varphi )$ out of the last integral; the remaining net mass flux is zero by steady continuity. Thus the equipotential luminosity conservation with stellar circulation isThis does not require $\mathbf{v}\cdot \hat{n}=0$ pointwise on $S$: inward and outward circulation may cancel in the net mass flux.

Finally, $f$ is constant on $S$ and the divergence theorem applied to the effective potential givesTherefore $-\chi dT/d\varphi =L/(4\pi Gm-2{\Omega }^{2}V)$. Since $dP/d\varphi =-\rho$,The mass $m$ is the mass enclosed by this surface, not necessarily the total stellar mass. The denominator is positive for the regular hydrostatic surfaces used here; near loss of centrifugal confinement the assumed stellar model must be reconsidered. The derivation uses slow circulation, retaining hydrostatic/equipotential thermodynamics to leading order while allowing a nonzero advective heat flux.

Let $\Psi$ be the gravitational potential, with ${\nabla }^2\Psi =4\pi G\rho$, and let ${\varpi }^2=x^2+y^2$ measure distance from the rotation axis. Uniform solid-body rotation contributes the centrifugal potential $-{\Omega }^2{\varpi }^2/2$. In the corotating frame, static force balance therefore gives the uniformly rotating stellar hydrostatic equilibrium equationsThe final coefficient follows because ${\nabla }^2(x^2+y^2)=4$.

A tangent displacement along a regular connected equipotential surface has zero dot product with $\nabla \varphi$, and hence with $\nabla P$. Thus $P=P(\varphi )$ locally. Taking the curl of the hydrostatic equilibrium equation gives $\nabla \rho \times \nabla \varphi =0$, so $\rho =\rho (\varphi )$ on such surfaces as well. The Poisson equation then shows that ${\nabla }^2\varphi$ is constant on each equipotential surface.

The magnitude $g=|\nabla \varphi |$ need not be constant. For instance $\varphi =A(x^2+y^2)+B{z}^2$, with positive $A\ne B$, has constant Laplacian but different pole and equator gradient magnitudes on an ellipsoidal equipotential. In a rotating star, flattening and centrifugal acceleration likewise make the effective gravity latitude dependent.

With uniform composition, an equation of state determines $T$ from $P$ and $\rho$; for the supplied ideal gas plus radiation pressure law its right side increases monotonically with $T$. Hence $T=T(\varphi )$, and the radiative conductivity $\chi (\rho ,T)$ is also constant on an equipotential surface. Put $f(\varphi )=-\chi dT/d\varphi$. Then $\mathbf{F}=f(\varphi )\nabla \varphi$. On the stellar surface, below mass shedding, the outward normal is $\nabla \varphi /g$, soThis is the radiative gravity-darkening law. The effective temperature measures emergent flux; it is not the same quantity as the interior material temperature that was constant on an equipotential.

A rapidly rotating star is oblate, with reduced effective gravity and cooler effective temperature around the equator. The poles are hottest. The original sketch uses a Roche model of a uniformly rotating star only to illustrate the shape and surface-flux variation:

The obstruction to purely local radiative equilibrium is visible by taking the divergence:On an equipotential surface, $\rho ϵ$, $f$, $f^{\prime }$ and the second term are constant, but $g^2$ generally varies. Unless $f^{\prime }=0$ or another special circumstance removes the mismatch, $\nabla \cdot \mathbf{F}=\rho ϵ$ cannot hold at every latitude. This is the radiative-equilibrium obstruction in a rotating barotropic star; heat advection by stellar meridional circulation resolves the missing energy transport in the leading hydrostatic treatment.

For the steady circulation, mass conservation gives $\nabla \cdot (\rho \mathbf{v})=0$. The specific enthalpy identity and hydrostatic equilibrium giveConsequentlyBoth $h$ and $\varphi$ are constant on the bounding equipotential surface. By the divergence theorem and steady continuity equation,Integrating the given thermal energy equation therefore gives equipotential luminosity conservation with stellar circulation:Individual circulation streams can cross the surface; only their net steady mass flux has to vanish.

Since $f$ is constant on that surface, integrating the Poisson equation givesThus $dT/d\varphi =-L/[\chi (4\pi Gm-2{\Omega }^{2}V)]$, while $dP/d\varphi =-\rho$. Dividing these derivatives and inserting $\chi =4ac{T}^3/(3\kappa \rho )$ yields the radiative gradient in a uniformly rotating star:The bracket must be positive for the regular outward-gravity configuration being used.

For an interior equipotential, both mass factors must be the enclosed mass $m$. The printed final bracket uses $M$, but its derivation gives $m$; at the outer surface $m=M$ and the printed bracket is recovered. For a direct check, a uniform-density model has $m=\rho V$, so the true correction ${\Omega }^{2}V/(2\pi Gm)={\Omega }^2/(2\pi G\rho )$ is independent of the chosen enclosed volume. Replacing $m$ by a fixed total mass would incorrectly make it volume dependent.

Assume nonrelativistic pressureless matter, an initial irrotational vector field of growing-mode linear cosmological density perturbations, a homogeneous expanding background, and no shell crossing. All gradients below refer to initial comoving coordinates. Let ${\rho }_0={\bar{\rho }}_m(t)a^3(t)={\bar{\rho }}_m(t_i)$, with $a(t_i)=1$, and normalize the linear growth factor by $D(t_i)=1$. The paper calls $D$ a growth rate, but it is the dimensionless growth factor; $\dot{D}$ is its time derivative.

Writing ${\mathbf{x}}_i$ for the actual position at $t_i$, the Zeldovich approximation is, to first order,This explicitly gives $\mathbf{x}(t_i)={\mathbf{x}}_i$. The usual notation $\mathbf{x}=\mathbf{q}+D\mathbf{s}(\mathbf{q})$ instead uses the uniform reference Lagrangian coordinate $\mathbf{q}$, with ${\mathbf{x}}_i=\mathbf{q}+\mathbf{s}(\mathbf{q})$. The two labels can be interchanged inside a first-order displacement, but not in the unperturbed position without the initial displacement correction.

To determine $\mathbf{s}$, mass conservation between the initial and current positions givesDemanding $\delta (t)=D(t){\delta }_i$ gives $\nabla \cdot \mathbf{s}=-{\delta }_i$. The irrotational vector field assumption makes $\mathbf{s}$ a gradient; the initial cosmological Poisson equation is ${\nabla }^2{\Phi }_i=4\pi G{\rho }_0{\delta }_i$. With the same boundary conditions for these potential equations, and after removing an irrelevant uniform translation,The evolution of $D$ follows from the linearized peculiar-motion equation $\ddot{\mathbf{x}}+2H\dot{\mathbf{x}}=-\nabla \Phi /a^2$. Since cosmological Poisson equation gives $\Phi (t)=(D/a){\Phi }_i$ to this order, substitution yields the linear growth equationThus the approximation extrapolates the linear growing displacement along a fixed initial direction, even as its density mapping becomes nonlinear.

The early spin of a dark matter halo comes from an external gravitational torque on its nonspherical initial mass region. A uniform external acceleration moves its centre of mass without spinning it; the spatially varying tidal tensor exerts different forces on different parts. The tidal torque theory requires misalignment of the region's shape and the surrounding tidal field. A local irrotational vector field of velocity does not imply zero integrated angular momentum for a nonspherical region.

Put $\mathbf{y}={\mathbf{x}}_i-{\hat{\mathbf{x}}}_i$ and use the physical peculiar velocity $\mathbf{v}=a\dot{\mathbf{x}}$. The physical lever arm is $a(\mathbf{x}-\hat{\mathbf{x}})$, while a leading-order mass element is ${\rho }_0{d}^3{x}_i$. Substituting the velocity formula into the supplied definition of angular momentum about the centre of mass givesTerms from the perturbed mass measure are higher order. Subtracting the barycentre's velocity also changes nothing because ${\int }_{V_L}\mathbf{y}{d}^3{x}_i=0$ at this order.

Apply the divergence theorem componentwise:The antisymmetry of the Levi-Civita symbol kills the last term. Consequently, for the ordinary oriented comoving surface element $d\mathbf{S}$,The original PDF prints a different coefficient, $-\dot{D}/(4\pi G{\bar{\rho }}_m{a}^2)$. That coefficient does not follow from its own mass measure and velocity equation and has the wrong dimensions for total angular momentum. The corrected coefficient above also gives the time dependence consistent with the paper's later tensor expression; this is a source typo, not a TeX transcription issue.

For a spherical $V_L$ centered on its barycentre, the outward normal is parallel to $\mathbf{y}$, so $\mathbf{y}\times d\mathbf{S}=0$ pointwise. HenceFor an equipotential surface, ${\Phi }_i={\Phi }_{\ast }$ on ${\Sigma }_L$, andThereforeThe second conclusion holds for any boundary shape at this order, not just a sphere.

To extract the leading tidal torque, expand the initial peculiar gravitational potential about the barycentre using its Taylor series:The constant has no force; $\mathbf{g}$ has zero integrated torque because the first mass moment vanishes. With the mass second-moment tensorwe obtainThe second equality swaps the indices $k,l$ and uses the symmetry of both tensors. Define the initial acceleration tidal tensor byThen the formula in the question has exactly the stated sign:If one defines $T$ as the positive Hessian matrix of ${\Phi }_i$ instead, the displayed contraction has a plus sign. Defining the convention is essential. Also, the paper's $I$ is a mass second-moment tensor; the mechanical inertia tensor is $\mathrm{tr}(I)1-I$.

Only the anisotropic parts contribute: a multiple of the identity produces zero contraction with the Levi-Civita symbol. In axes where $I=\mathrm{diag}(I_1,I_2,I_3)$, for example,Thus a spherical region or tensors sharing principal axes have zero leading torque; unequal principal moments together with off-diagonal tidal components create spin. The Taylor series truncation assumes that higher spatial derivatives of the tidal field are sufficiently small across the region.

In an Einstein-de Sitter universe, $a\propto t^{2/3}$, $D\propto a$, and ${\bar{\rho }}_m\propto a^{-3}$. The initial $I$ and $T$ are time independent, soThis is early-time growth for fixed initial matter, before collapse invalidates the extrapolation. It does not predict indefinite linear spin growth for a virialized halo.

Two reasons for only approximate agreement with simulations are that nonlinear collapse and shell crossing change the trajectories and the tidal field, invalidating the first-order displacement and fixed leading tidal tensor; and that real halos undergo dark-matter halo mergers, anisotropic accretion, and exchange of material and angular momentum, so the halo identified at a later time need not be the same isolated collection of initial particles. These processes can change both the magnitude and the direction of the spin.

Use cylindrical coordinates $(R,\varphi ,z)$, put $s=(R^2+z^2{)}^{1/2}$, and write the purely azimuthal velocity as $\mathbf{u}=R\Omega {\mathbf{e}}_{\varphi }$. The steady radial, azimuthal, and vertical momentum equations areThus $\nabla p=\rho {\mathbf{g}}_{\mathrm{eff}}$, where the effective gravity in a rotating frame isWhen radiation pressure supplies the pressure support, force balance gives

On the upper conical surface $z=R\tan \alpha$, the tangent and outward normal areBecause a pressure surface is an equipotential surface, ${\mathbf{g}}_{\mathrm{eff}}\cdot \mathbf{t}=0$. Since $s=R/\cos \alpha$, this fixesThe lower face gives the reflected result. One face has area element $dA=2\pi R,dR/\cos \alpha$. Adding both faces and using $F=c|{\mathbf{g}}_{\mathrm{eff}}|/\kappa$ yieldsConsequentlyFor a broad radial range, each logarithmic interval can radiate an Eddington-scale contribution because rotation reduces the surface-normal gravity. The factor $\sin \alpha$ makes a thin wedge faint and lets a geometrically thick funnel expose a larger effective gravity and emitting area.

In the synchronously rotating reference frame, stationary fluid obeyswhere $\Psi$ is the sum of the two gravitational potentials and the centrifugal potential. The circular binary's Kepler third law gives ${\Omega }^2=G(M_1+M_2)/a^3=GM/a^3$, hence

Uniform composition and the assumed central condensation make star 1 approximately barotropic, so one may define the specific enthalpyHydrostatic balance becomes $\nabla (h+\Psi )=0$, or $h+\Psi =\text{constant}$ throughout the star. Since $h$ is a monotone function of $P$ and $\rho$, constant-$P$ and constant-$\rho$ surfaces are equipotential surfaces of $\Psi$.

On a regular equipotential surface, $f,f^{\prime },\rho$ and ${\nabla }^2\varphi$ are constant, but $|\nabla \varphi {|}^2$ generally varies. Consequently a flux $\mathbf{F}=f(\varphi )\nabla \varphi$ cannot ordinarily satisfy local radiative equilibrium $\nabla \cdot \mathbf{F}=\rho ϵ$ at all latitudes. Special cases such as $f^{\prime }=0$ or spherical symmetry evade this obstruction. Stellar meridional circulation supplies heat advection in the leading hydrostatic treatment.

The divergence theorem and Poisson equation give ${\int }_S\nabla \varphi \cdot d\mathbf{S}=4\pi Gm-2{\Omega }^{2}V$ on a bounding equipotential surface. Equipotential luminosity conservation with stellar circulation then fixes $f=L/(4\pi Gm-2{\Omega }^{2}V)$. Combining $dT/d\varphi =-f/\chi$, $dP/d\varphi =-\rho$ and radiative conductivity gives the displayed logarithmic temperature-pressure gradient. Both appearances of mass mean the enclosed mass $m$; only on the outer boundary can both be replaced by the total mass.

Uniform solid-body rotation adds the centrifugal potential $-{\Omega }^2{\varpi }^2/2$ to the gravitational potential. In the corotating frame, hydrostatic equilibrium becomes $\nabla P=-\rho \nabla \varphi$. The Poisson equation gives ${\nabla }^2\varphi =4\pi G\rho -2{\Omega }^2$. On regular connected equipotential surfaces, $P$ is constant; taking the curl shows that $\rho$ is constant too. Uniform composition and an invertible equation of state then make $T$ constant on each such surface. These statements do not require $|\nabla \varphi |$ to be constant there.