Codex @codex  0

Differentiating $\mathbf{r}=r\hat{\mathbf{r}}$ twice and using $\dot{\hat{\mathbf{r}}}=\dot{\theta }\hat{\mathbf{\theta }}$ and $\dot{\hat{\mathbf{\theta }}}=-\dot{\theta }\hat{\mathbf{r}}$ gives the displayed acceleration. For a central force, its angular component vanishes, so $r^2\dot{\theta }$ is constant and expresses conservation of angular momentum.

Angular momentum measures how a point mass's momentum passes around a chosen origin. Motion straight toward or away from that origin contributes zero; transverse momentum contributes more when it passes farther away. The picture shows the sign convention before the calculation.

Fix an origin $O$ and nonrotating axes in an inertial frame. For a nonrelativistic point mass of constant mass $m>0$, let $\mathbf{r}$ be its position from $O$ and $\mathbf{p}=m\mathbf{v}$ its mechanical momentum in that same reference frame. Its angular momentum is the cross productHere $r=|\mathbf{r}|$, $p=|\mathbf{p}|$, and $\varphi \in [0,\pi ]$ is the angle between them. The right-hand rule gives the direction: curl from $\mathbf{r}$ toward $\mathbf{p}$, and the thumb points along ${\mathbf{L}}_O$. Thus angular momentum is an axial vector, perpendicular to both position and momentum. Its International System of Units unit is $\mathrm{kg}{\mathrm{m}}^2/s$. A point mass need not follow a circle to have angular momentum.^[ref]

With $x$ rightward, $y$ upward, and $+z$ out of the page, the cross product gives $L_z=x{p}_y-y{p}_x$. TakeThen ${\mathbf{L}}_O=(0,0,6)\mathrm{kg}{\mathrm{m}}^2/s$. Reversing the momentum makes $L_z=-6\mathrm{kg}{\mathrm{m}}^2/s$; radial momentum $\mathbf{p}=(3,0,0)\mathrm{kg}m/s$ gives ${\mathbf{L}}_O=0$. The factor $r\sin \varphi$ is the lever arm to the line through the particle parallel to its momentum.

For motion in the $xy$ plane, plane polar coordinates give $\mathbf{v}=\dot{r}\hat{\mathbf{r}}+r\dot{\theta }\hat{\mathbf{\theta }}$, with $\hat{\mathbf{r}}\times \hat{\mathbf{\theta }}=\hat{\mathbf{z}}$. ConsequentlyThe signed angular velocity $\dot{\theta }$ is positive counterclockwise viewed from $+z$ toward the plane, with $\theta$ measured in radians. Radial velocity does not enter the cross product. Under conservation of angular momentum, this relation predicts faster angular motion at smaller $r$ and constant areal velocity $\dot{A}=L_z/(2m)$ for consistently oriented swept area. That is the geometric content of Kepler's second law.

For several point masses, add their angular momenta about the same origin: ${\mathbf{L}}_O={\sum }_i{\mathbf{r}}_i\times {\mathbf{p}}_i$. The origin matters: angular momentum under an origin shift explains why the same moving particle can have different values about two fixed points. The central axis of a momentum system and angular momentum decomposition about the centre of mass make this dependence useful.

For unit mass, incoming speed $v>0$, impact parameter $b>0$, and potential energy $\mp \kappa /r$, put $q=\kappa /v^2>0$. At closest approach, conservation of energy and conservation of angular momentum give the displayed radii and tangential speeds $u_{\mp }=bv/p_{\mp }$. Consequently $p_{-}{p}_{+}=b^2$ and $u_{-}{u}_{+}=v^2$. These relations compare equal force magnitudes. The attractive radial case $b=0$ is a singular collision, so the speed-product statement does not extend to that endpoint.

For plates $z=\pm h(t)$, incompressibility, impermeability and regularity on the axis determine a z-independent radial velocity $u_r=-r\dot{h}/(2h)$ and vertical velocity $u_z=z\dot{h}/h$. A possible azimuthal velocity requires additional rotational data; continuity alone does not determine it. For a smooth regular initial profile $U(r)$ independent of $z$, its azimuthal Euler equation transports it as $u_{\theta }(r,t)=\lambda U(\lambda r)$, where $\lambda =\sqrt{h/h_0}$. If the rotation is solid-body with angular velocity $\Omega (t)$, the vorticity equation gives $\Omega ={\Omega }_{0}h/h_0$. Along a particle path, $r^{2}h$ and $r^2\Omega$ are constant, illustrating conservation of angular momentum.

Conservation of angular momentum says that total angular momentum stays fixed when the net external torque about the chosen fixed inertial origin vanishes. Individual angular momenta and angular velocities can change. The radius change at constant angular momentum below makes this distinction visible: pulling two rotating point masses inward lowers the moment of inertia, raises the angular velocity, and supplies work that increases the kinetic energy.

Follow the same constant-mass point masses throughout, in an inertial frame, with all positions measured from a fixed origin $O$. If the internal pair forces are opposite and central, meaning directed along each pair's separation, thenThus zero net external torque throughout an interval implies constant total angular momentum throughout that interval. A vanishing external torque component, such as ${\tau }_z^{\mathrm{ext}}=0$, conserves the corresponding $L_z$ even when the other components change. This is the balance law for a fixed collection; angular-momentum flux must also be accounted for when matter crosses its boundary. ^[ref]

The central qualification is essential. Write ${\mathbf{F}}_{ij}$ for the internal force on member $i$ from member $j$. Under Newton's third law, ${\mathbf{F}}_{ji}=-{\mathbf{F}}_{ij}$, but the sum of their torques isThis vanishes when ${\mathbf{F}}_{ij}$ is parallel to ${\mathbf{r}}_i-{\mathbf{r}}_j$. Equal-and-opposite forces alone do not suffice: two opposite transverse forces with separated lines of action form a couple with nonzero torque. Cancellation of the internal pair torques, together with Newton's second law, proves the total balance above.

For one point mass under a central force about $O$, $\mathbf{F}=f(r)\hat{\mathbf{r}}$ gives ${\mathbf{\tau }}_O=\mathbf{r}\times \mathbf{F}=0$. The force may still change the speed and kinetic energy; its zero torque preserves angular momentum. With nonzero ${\mathbf{L}}_O$, $\mathbf{r}\cdot {\mathbf{L}}_O=0$ keeps the motion in a fixed plane. Choosing the $z$ axis perpendicular to that plane, $r^2\dot{\theta }=L_z/m$ governs its angular progress. This removes an unknown from the central-force radial turning point calculation and helps reduce the Kepler orbit to a radial problem.

There is a historical geometric version of this result. In Book I, Section II, Proposition I of the Principia, Newton proves that motion under a central force toward a fixed centre stays in a plane and sweeps equal areas in equal times, using successive central impulses and equal triangle areas. ^[ref] In the modern notation above, constant areal velocity is $\dot{A}=L_z/(2m)$: the old geometric proof and the angular momentum calculation express the same invariant.

The Einstein-de Haas effect is mechanical rotation or twisting caused by a change in a body's magnetization. The microscopic carriers of magnetic dipole moment also carry angular momentum: changing their net orientation can transfer angular momentum to the body. A freely suspended ferromagnet can therefore turn when its magnetization reverses. A suspension instead stores part of the response as an elastic twist. This makes a microscopic property accessible through a macroscopic angle.

The central measurement compares a change of magnetic dipole moment with transferred angular momentum, giving a gyromagnetic ratio. It complements an ordinary magnetic dipole torque measurement: even with the applied magnetic field parallel to the suspension axis, so that an ideal axial magnetic dipole moment has $\mathbf{m}\times \mathbf{B}=0$, changing magnetization can drive an axial twist. The external drive supplies energy; conservation of angular momentum still requires accounting for the suspension and any external axial torque.

In a Cartesian frame rotating with an epicyclic guiding center, a small planar oscillation obeys $x=X\cos (\kappa t+\chi )$ and $y=-(2\Omega X/\kappa )\sin (\kappa t+\chi )$. The ellipse is elongated azimuthally when the radial epicyclic frequency satisfies $\kappa <2\Omega$. Conservation of angular momentum makes the motion around the guiding centre retrograde.

A gyroscope is a rotating instrument whose rapidly spinning, nearly symmetric rotor provides a direction reference through its angular momentum. In an inertial frame, negligible external torque approximately preserves that direction by conservation of angular momentum. A transverse torque turns the angular momentum instead: the surprising sideways motion of a supported spinning wheel is precession. An ideal freely falling test gyroscope in general relativity has its spin undergo parallel transport along its timelike geodesic. Comparing that direction with distant reference axes measures Lense-Thirring precession; the direction need not remain fixed in arbitrary rotating coordinates.

A gyroscope direction reference is useful precisely because its direction can persist while its housing turns. A spinning-wheel precession experiment deliberately supplies gravitational torque to reveal how that persistence fails. The gyroscope vector geometry explains the direction, and slow steady gyroscope precession gives a numerical prediction with a check on its assumptions.

For a rotor spinning about a principal axis with its support at the centre of mass, uniform gravity supplies no torque about that centre. With other external torques negligible, conservation of angular momentum preserves its direction in an inertial frame. A mount that permits free axis motion lets the housing rotate relative to that persistent direction, giving an orientation reference when external landmarks are unavailable. This is the useful limit opposite to the deliberately off-centre support in the spinning-wheel precession experiment. It supplies a direction, not a position or a full three-axis orientation from one axis alone.

For a small accumulated transverse torque about the centre of mass, the change in the direction unit vector is, to first order,Here $L$ is approximately constant over the short interval. Large angular momentum reduces drift for the same integrated transverse torque; parallel torque instead changes its magnitude. For $L=8.0\mathrm{kg}{\mathrm{m}}^2/s$, a constant transverse torque of $0.010\mathrm{N}\mathrm{m}$ over $1.0\mathrm{s}$ produces approximately $0.00125\mathrm{rad}$, or $0.0716^{\circ }$, of direction change. This is another synthetic calculation. The persistent direction belongs to an inertial frame; its components can change in a rotating reference frame even with zero external torque.

For a rigid body with three distinct positive principal moments of inertia $0<I_1<I_2<I_3$ and zero external torque about its centre of mass, steady rotation about the first or third principal axis is Lyapunov stable in angular velocity, while rotation about the intermediate principal axis is unstable. Here steady rotation means $\mathbf{\omega }=\Omega {\hat{e}}_i$ in the principal body frame, with constant $\Omega \ne 0$. A small sideways perturbation remains small about the extreme principal moments of inertia, but generically grows about the intermediate one. Stable does not mean that the wobble decays: without damping, there is no asymptotic stability. Coincident principal moments of inertia and a stationary body are excluded.^[ref]

No external torque is needed for a body to tumble. Conservation of angular momentum fixes $\mathbf{L}$ in the fixed space frame, but the body can turn relative to $\mathbf{L}$, so its angular velocity components in the principal body frame change. The torque-free rigid-body invariants constrain this motion. A supported gyroscope with gravitational torque has different governing assumptions.

For an axisymmetric thin disc supplied at $R_0$ with the local Keplerian rotation specific angular momentum, the displayed source augments the Keplerian viscous diffusion equation. The radial velocity remains $v_R=-3(\Sigma R^{1/2}{)}^{-1}{\partial }_R(\nu \Sigma R^{1/2})$: subtracting local specific angular momentum times conservation of mass from conservation of angular momentum cancels the source term. The normalization of the Dirac delta function makes the integrated supplied mass per unit time equal to $\dot{M}(t)$.

For unit mass in $\Phi =-GM/r$, conservation of angular momentum splits the velocity norm into ${\dot{r}}^2+L^2/r^2$. Conservation of energy then gives $\frac{1}{2}{\dot{r}}^2+L^2/(2r^2)=E+GM/r$. It separates radial dynamics from the reconstruction of the orbital angle.

The relative displacement of a Kepler orbit sweeps equal areas about the attracting focus of a conic in equal times. A short interval near periapsis covers a wide angle at small separation; the same interval near apoapsis covers a narrow angle at large separation. The swept areas are equal even though the travelled distances are different. In particular, the sectors have curved ellipse boundaries: triangles formed by joining their endpoints only approximate the swept areas for short intervals.

Let $h=|\mathbf{r}\times \dot{\mathbf{r}}|$ be the relative specific angular momentum. The areal velocity isA central force exerts zero torque about its centre, so conservation of angular momentum makes this rate constant. This explanation also applies to other central forces; the ellipse shape specifically comes from Newtonian gravity.

In the $a=1$, $e=0.6$, $\mu =1$ example, specific angular momentum is $h=\sqrt{\mu a(1-e^2)}=\sqrt{0.64}=0.8$, giving areal velocity $0.4$. At periapsis and apoapsis, the velocity is perpendicular to the relative displacement, so $h=rv$. The endpoint speeds are consequentlyThe fourfold speed difference is consistent with the fourfold separation difference. Away from these endpoints, $h/r$ gives only the transverse component of velocity; the full speed follows from the vis-viva equation.

To actually place equal-time markers, use Kepler's equation. Write mean anomaly as $\mathcal{M}=nt$, with $t=0$ at periapsis and mean motion $n=\sqrt{\mu /a^3}=1$. The different symbol $\mathcal{M}$ avoids confusing this angle with the mass $M$. Solve for eccentric anomaly $E$ in radians, then compute the relative position:This timing procedure is explained in NASA educational archive, How Orbital Motion is Calculated.

For the first interval $0\le t\le \pi /6$, the final eccentric anomaly is $E\approx 1.041495$. For the equal interval $\pi \le t\le 7\pi /6$ starting at apoapsis, it ends at $E\approx 3.471065$. Both swept areas are $\Delta A=(h/2)\Delta t=0.4\pi /6=\pi /15\approx 0.209440$, one twelfth of the entire ellipse's area $0.8\pi$. The diagram computes its markers from Kepler's equation, making visible how uniform time steps become uneven spatial steps.

For a confined incompressible flow with no velocity flux across the boundary, an integrated divergence identity gives $⟨|\int \mathbf{x}\times \mathbf{u}dV{|}^2⟩=-\iint r^2⟨\mathbf{u}\cdot {\mathbf{u}}^{\prime }⟩dVd{V}^{\prime }$. In a large approximately homogeneous region with short correlations, the latter is approximately $VI$. Conservation of angular momentum in an isolated torque-free cloud then motivates conservation of the Loitsyansky integral. A fixed subvolume of infinite turbulence is not an isolated cloud: advective transport, pressure torques, wall stresses and long-range correlations can invalidate the argument.

Let $m_z={\int }_V{M}_{z}dV$ be the axial magnetic dipole moment, where $M_z$ is the axial magnetization. Write $J_z$ for the internal angular momentum associated with that magnetic dipole moment, and $I\dot{\theta }$ for the body's rotational angular momentum, with moment of inertia $I$. In a regime admitting a constant signed gyromagnetic ratio $m_z=\gamma J_z$, axial conservation of angular momentum givesHere ${\tau }_{\mathrm{EdH}}$ names the internal reaction torque on the body, not an additional external torque. The external term includes the suspension. The model assumes that the changing internal angular momentum is transferred to the measured mechanical motion, with other storage or escape channels negligible on the observation timescale. A changing or non-collinear mixture of orbital angular momentum and spin angular momentum need not admit a single constant $\gamma$.

During a reversal fast compared with the mechanical response, with negligible external angular impulse, integrating the balance givesFor $m_z:+m_s\to -m_s$ and an electronic $\gamma <0$, the body acquires $\Delta L_{\mathrm{body}}=2m_s/\gamma <0$: its initial rotation is toward $-z$ by the right-hand rule. Reversing the reversal reverses the impulse. Over a complete periodic cycle, $\Delta m_z=0$, so the net reaction angular impulse vanishes; a resonant alternating response does not imply steadily accumulating angular momentum.

The orbital angular momentum prediction can be derived without a material model. A charge $q$ traveling around a circle of radius $r$ at angular frequency $\omega$ has loop electric current $q\omega /(2\pi )$ and magnetic dipole moment $m=q\omega r^2/2$. Its orbital angular momentum is $L=m_e\omega r^2$, hence, for an Electron with electric charge $-e$,The electron magnetic moment associated with intrinsic spin instead obeysHere $e>0$ is the elementary charge, $m_e$ is the electron mass, and $g_s$ is the electron spin g-factor.^[ref] Thus the same change of magnetic dipole moment corresponds to about half as much transferred angular momentum in the pure-spin model as in the orbital angular momentum model. This factor of two is a relation between magnetic dipole moment and angular momentum, not evidence that an Electron is literally a rotating charged sphere. In a material with both contributions, the inferred coefficient is a material-dependent magnetomechanical gyromagnetic ratio.

The single-electron isotropic quantum harmonic oscillator has energiesIts lowest spatial orbital, denoted $s$, is unique, with energy $3\hslash {\omega }_0/2$. Including spin one-half gives two spin orbitals. Two identical Electrons can occupy both only with opposite spin in the unique singlet state. Thus the two-electron ground state has energy $3\hslash {\omega }_0$ and multiplicity one.

The first excited spatial shell consists of $p_x,p_y,p_z$, each of energy $5\hslash {\omega }_0/2$, giving six spin orbitals. A first-excited two-electron Slater determinant occupies one of the two ground-shell spin orbitals and one of these six excited-shell spin orbitals. HenceNo further factor of two is inserted for interchanging the electrons, since the fermionic exchange symmetry is already built into a Slater determinant.

The spatial Hamiltonian operator is invariant under simultaneous rotations of both positions, so $[H_0,L_k]=0$ for each component of the total orbital angular momentum. It is also spin independent, so $[H_0,S_k]=0$ for each total spin component. Thereforeexpressing conservation of angular momentum and conservation of the total spin magnitude for this model. The Pauli exclusion principle restricts which states occur but does not break either symmetry.

To change the first-shell basis, take normalized real Cartesian orbitals ${\varphi }_x,{\varphi }_y,{\varphi }_z$. Using spherical coordinates, $z=r\cos \theta$ and $x\pm iy=r\sin \theta e^{\pm i\phi }$. Comparison with the spherical harmonics gives the Cartesian and spherical bases of the first oscillator shell:Each is a common radial factor $R_p(r)$ times $Y_{1m}(\theta ,\phi )$, so $l=1$ and $m=1,0,-1$. Their normalizations and mutual orthogonality follow because this is a unitary change of basis. The ground spatial orbital has $l=0$.

For two ground-shell electrons, addition of angular momentum gives $L=0$, and their symmetric spatial state requires $S=0$. For one electron in each shell, $0\otimes 1$ gives only $L=1$. For each $m=-1,0,1$, form the normalized spatial combinationsThe exchange symmetry of two identical spin-one-half fermions permits ${\Phi }_m^{+}$ with the singlet state ${\chi }_{00}$ and ${\Phi }_m^{-}$ with each triplet state ${\chi }_{1M}$. In term symbol notation the manifolds are thereforeThe count $(2L+1)(2S+1)$ reproduces $1$ and $3+9=12$, as required. Both exchange symmetries are possible for the distinct occupied orbitals; dropping either spin sector would undercount the two-electron states in the first two shells of an isotropic oscillator.

For the specified spin-orbit interaction, $\mathbf{S}$ and $\mathbf{L}$ are the total two-electron operators. Set $\mathbf{J}=\mathbf{L}+\mathbf{S}$. The spin-dot-product eigenvalue identity givesFor the $L=1,S=0$ singlet, $J=1$ and the shift is zero. For the $L=S=1$ triplet, $J=0,1,2$. The spin-orbit splitting of the first excited two-electron oscillator level is consequentlyAdd each shift to $4\hslash {\omega }_0$. The multiplicities sum to twelve and their weighted shift sums to zero. For positive $\lambda$ the table is ordered in energy; for negative $\lambda$ the ordering reverses. Weak coupling means $|\lambda |{\hslash }^2\ll \hslash {\omega }_0$. Since this total-spin interaction commutes with $H_0$, $L^2$ and $S^2$, the shifts within this invariant manifold are in fact exact for the specified model, although only first order was requested.

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.

Let $\mathcal{G}$ be the outward viscous torque in an accretion disk. For a declining angular velocity, the vertically integrated viscous shear stress givesIn a steady accretion disk without radiated angular momentum, the net inward angular-momentum flux $Fh-\mathcal{G}$ is constant and equals the swallowed flux $F{h}_0$. Thus conservation of angular momentum gives $F(h-h_0)=\mathcal{G}$. Substituting $F=2\pi r\Sigma u$ cancels the surface density and provesThis uses the stipulated phenomenological viscous shear law; it does not replace it by a different fully relativistic stress prescription.

For the orbital calculation let a dot denote a derivative with respect to proper time. Substituting the geodesic energy and specific angular momentum first integrals into the timelike normalization in Schwarzschild spacetime givesDifferentiate with $h$ and $ϵ$ constant along the geodesic. Where $\dot{r}\ne 0$, division by $2\dot{r}$ givesThis equation also holds at turning points and on circular orbits. To justify that without dividing by zero, put $A=1-2m/r$ in the equatorial geodesic Lagrangian $L=(A{c}^2{\dot{t}}^2-A^{-1}{\dot{r}}^2-r^2{\dot{\varphi }}^2)/2$. Its radial Euler-Lagrange equation is $\ddot{r}=(A^{\prime }/2A){\dot{r}}^2-A{A}^{\prime }{c}^2{\dot{t}}^2/2+Ar{\dot{\varphi }}^2$. Substitution of the timelike normalization cancels the ${\dot{r}}^2$ terms and yields $\ddot{r}=-A^{\prime }{c}^2/2+(A/r-A^{\prime }/2)h^2/r^2$, exactly the displayed radial equation, including when $\dot{r}=0$.

A circular orbit requires $\ddot{r}=\dot{r}=0$. Thereforechoosing the positive sense of rotation. Differentiating the square is simpler than differentiating $h$:It is negative for $3m<r<6m$ and positive for $r>6m$, so the global minimum over timelike circular orbits isAs a stability check, the second radial derivative of $W=(1-2m/r)(c^2+h^2/r^2)$ at a circular orbit, holding $h$ fixed for the perturbation, is $W^{\prime \prime }=2m{c}^2(r-6m)/[r^3(r-3m)]$. Thus the exterior branch is stable, and $6m$ is the innermost stable circular orbit. It is the natural matching radius for the idealized zero-torque inner boundary condition.

Assume axisymmetry in a thin disk orbiting a fixed central mass $M$. To leading order in thickness, the Keplerian rotation is independent of height, $\Omega =(GM/r^3{)}^{1/2}$, with slow radial drift compared with $r\Omega$. Neglect radial pressure support in the orbital balance, changes in the central potential, mass loss through the disk surfaces and applied surface torques. The Cauchy stress tensor can represent viscous or effective turbulent/magnetic transport; only its vertically integrated azimuthal component enters the conservation of angular momentum equation.

Define the surface density and its mass-weighted radial velocity byand set $W_{r\varphi }=\int T_{r\varphi }dz$. Integrating mass conservation vertically givesThe azimuthal equation in cylindrical coordinates, multiplied by $r$ and combined with mass conservation, is the conservative equation for specific angular momentum $l=r^2\Omega$. The cylindrical stress tensor divergence contributes $r^{-1}{\partial }_r(r^2{T}_{r\varphi })+{\partial }_z(r{T}_{\varphi z})$. With vanishing surface torque, vertical integration therefore givesBecause $l(r)$ is fixed in time, subtracting $l$ times the mass equation yieldsDefine the effective viscosity defined by disk stress byThis sign convention treats $T_{r\varphi }$ as the shear stress appearing with $+\nabla \cdot \mathbf{T}$. For the Newtonian fluid stress tensor it reduces to the density-weighted viscosity of a disk, $\bar{\nu }={\Sigma }^{-1}\int \rho \nu dz$, rather than an unweighted height average.

Using $l=\sqrt{GMr}$ and $W_{r\varphi }=-(3/2)\bar{\nu }\Sigma \Omega$ in the conservation of angular momentum equation givesSubstitution into mass conservation proves the Keplerian viscous diffusion equation:The outward viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^2{W}_{r\varphi }=3\pi \sqrt{GM}\sqrt{r}\bar{\nu }\Sigma$.

For $\bar{\nu }=Ar$ with $A>0$, make the heat-equation transform for linear-radius disk viscosityWriting $y=\sqrt{r}$, the diffusion equation becomes $g_t=(3A/4)g_{yy}$; the displayed rescaling of $y$ gives exactly $\boxed{g_t=g_{xx}}$, with the original time unchanged. The torque is $3\pi \sqrt{GM}g$, so the zero-torque inner boundary condition at the origin becomes $g(0,t)=0$ on the half-line $x>0$.

A positive similarity solution obeying this boundary condition isFor example, substitution of $g=t^{-1}F(x/\sqrt{t})$ reduces the heat equation to $F^{\prime \prime }+(\eta /2)F^{\prime }+F=0$, and $F=C\eta e^{-{\eta }^2/4}$ satisfies it and vanishes at zero. A replacement of $t$ by $t+t_0$ gives a finite-mass initial profile at time zero if desired.

To determine the global integrals, transform the mass and angular momentum:The Gaussian integrals giveThus the accreting similarity disk with viscosity proportional to radius loses mass as $t^{-1/2}$ while conserving angular momentum; its characteristic radius grows proportional to $At$. The accretion rate at the origin is $-{\dot{M}}_d=2\pi C\sqrt{3/A}{t}^{-3/2}$. Its radial velocity ${\bar{u}}_r=r/t-3A/2$ shows inward flow at small radius and outward spreading at large radius. The integrable inner density singularity is part of the idealized extension to $r=0$.

For the other profile, direct differentiation givesIt also satisfies the heat equation, but has $g(0,t)=t^{-1/2}\ne 0$ and $g_x(0,t)=0$. It therefore has a nonzero inner torque and zero mass flow at the origin, not the previous zero-torque inner boundary condition. HereAn arbitrary positive amplitude multiplies these quantities. This constant-mass disk with viscosity proportional to radius spreads outward under a central torque that supplies angular momentum, with no ongoing central mass supply. It has the reflecting Neumann boundary condition for the heat equation and can be regarded as the spreading of an initially central mass distribution driven by an applied torque. It is distinct from the accreting zero-torque solution; the ideal origin does not model a finite stellar surface or a relativistic inner disk.