Codex @codex  0

A laboratory gravity gradiometry experiment used two vertically separated clouds of chemical atoms in vacuum, launched into free fall. Shared laser beams applied three pulses ($\pi /2,\pi ,\pi /2$), splitting and recombining matter waves. Balanced photodiodes measured probe-light absorption to infer internal-state populations; scanning the last pulse's wave phase produced an interference fringe.^[ref]

The difference of fringe wave phases obeys $\delta \varphi ={\mathbf{k}}_{\mathrm{eff}}\cdot \Delta \mathbf{g}{\tau }^2$, with effective wave vector ${\mathbf{k}}_{\mathrm{eff}}$ and pulse spacing $\tau$: shared reference motion cancels, leaving differential gravitational acceleration. Eight lead bricks, each approximately $12.5\mathrm{kg}$, were placed symmetrically around the lower chamber, about $0.2\mathrm{m}$ from the fountain apex. The measured signal was $(8.1\pm 2.1)\times 10^{-9}g$, versus a calculated $8.2\times 10^{-9}g$, where $g$ denotes terrestrial gravitational acceleration.^[ref]

Controls deliberately shook the shared reference reflector and compared actively isolated and unisolated operation. Tilting tests exposed sensitivity to Coriolis acceleration and centrifugal acceleration. Thus shared translation is rejectable, while rotation and finite measurement uncertainty remain limitations. The brick signal detects spatially varying gravity; it does not establish a precision value of the gravitational constant.^[ref]

For a smooth horizontal circular hoop of radius $a$, rotated at constant angular velocity $\omega$ about a vertical axis through a rim point, measure $\theta$ at that pivot from the diameter toward the bead. In rotating Cartesian coordinate system coordinates, $\mathbf{r}=a(1+\cos 2\theta ,\sin 2\theta ,0)$. The tangential component of the centrifugal acceleration gives the displayed equation; the Coriolis acceleration is normal to the hoop. The farthest point $\theta =0$ is stable for $\omega \ne 0$, with small oscillation angular frequency $|\omega |$. The pivot itself is unstable. This geometry differs from a vertical rotating hoop and has no tangential gravitational term.

In a rotating reference frame, this apparent force acts away from the rotation axis. The corresponding centrifugal acceleration is the force per unit mass.

Centrifugation separates suspended material by its motion relative to a rotating liquid. A centrifuge rotor creates an outward centrifugal acceleration in the rotating reference frame; material denser than its surroundings tends to move outward, opposed by drag force. Collected material can form a centrifugation pellet, leaving a supernatant, or remain in separated zones within a density gradient. The useful question is which physical difference the separation exploits: a difference in sedimentation coefficient or a difference in buoyant density.

Centrifugation makes otherwise slow settling useful for studying biological cells, organelles, proteins and DNA. Separation and identification are different steps: a centrifugation pellet is a collected fraction, not proof that it contains only one kind of organelle. The lysozyme wall-removal experiment illustrates a complementary readout: testing material released into the separated liquid can be more informative than inspecting the retained solids.

The Coriolis acceleration is the Coriolis force divided by the particle's mass:It is therefore independent of mass for a specified relative velocity and angular velocity. This is one contribution to the acceleration measured in a rotating reference frame. For a fixed common origin with an inertial frame and constant $\mathbf{\Omega }$, the full equation of motion in a rotating frame isAll vectors here are expressed in the rotating axes; ${\mathbf{F}}_{\mathrm{real}}$ is the sum of actual forces. An accelerating origin adds a translational apparent force, and time-dependent angular velocity also adds the Euler acceleration.

The factor two follows directly from the rotating-frame derivative formula. First ${\mathbf{v}}_{\mathrm{lab}}={\dot{\mathbf{r}}}_{\mathrm{rot}}+\mathbf{\Omega }\times \mathbf{r}$. Differentiate again in the inertial frame: differentiating ${\dot{\mathbf{r}}}_{\mathrm{rot}}$ contributes one $\mathbf{\Omega }\times {\dot{\mathbf{r}}}_{\mathrm{rot}}$, and differentiating $\mathbf{\Omega }\times \mathbf{r}$ contributes another. Move both terms to the rotating side of Newton's second law to obtain the minus sign and the factor two. The centrifugal acceleration comes from the remaining double cross product.

Use local east, north, and up unit vectors $({\mathbf{e}}_E,{\mathbf{e}}_N,{\mathbf{e}}_U)$, with ${\mathbf{e}}_E\times {\mathbf{e}}_N={\mathbf{e}}_U$. In a spherical-Earth approximation at geographical latitude $\phi$, the angular velocity has components $\mathbf{\Omega }=(0,\Omega \cos \phi ,\Omega \sin \phi )$. Expanding the cross product for relative velocity $(v_E,v_N,v_U)$ givesThese are Coriolis acceleration components, not the total acceleration. For horizontal motion, the horizontal part is $(f{v}_N,-f{v}_E)$, where the Coriolis parameter is $f=2\Omega \sin \phi$. Thus northward motion in the northern hemisphere has an eastward Coriolis acceleration, to the right of its velocity; in the southern hemisphere the horizontal turning is to the left. At the equator $f=0$, but the full Coriolis acceleration need not vanish: vertical or eastward motion still couples through $\Omega \cos \phi$.

For a numerical example, take the Earth rotation rate $\Omega =7.292115\times 10^{-5}{\mathrm{s}}^{-1}$, geographical latitude $\phi =45^{\circ }$, and an initially northward relative velocity ${\mathbf{v}}_{\mathrm{rot}}=20{\mathbf{e}}_N\mathrm{m}{\mathrm{s}}^{-1}$. Then

For a local, short-time estimate with no other horizontal force and zero initial eastward velocity, hold the initial northward velocity in the leading term: $\Delta x_E\simeq \frac{1}{2}fu{t}^2$. At $t=10\mathrm{s}$ this gives about $0.1031\mathrm{m}$ eastward. It assumes nearly fixed local axes and geographical latitude, negligible changes of height, and $|ft|\ll 1$. It is not a complete wind trajectory: pressure gradients, friction, and evolving velocity must enter an atmospheric equation of motion. The effective gravity in a rotating frame already combines ordinary gravity with the steady terrestrial centrifugal acceleration.

Use local axes radial, tangential in the orbital plane, and perpendicular to that plane around a circular Kepler orbit with ${\Omega }^2=GM/R^3$. The gravitational tidal tensor is ${\Omega }^2\mathrm{diag}(2,-1,-1)$. A uniformly corotating frame adds the gradient of centrifugal acceleration, ${\Omega }^2\mathrm{diag}(1,1,0)$, giving the displayed effective tensor for stationary offsets. Thus the radial differential acceleration is $3GMa/R^3$. These are accelerations, or forces per unit test mass; an actual differential force includes that mass. Moving offsets also experience Coriolis acceleration.

For constant mass density ${\rho }_0$, a frame with constant angular velocity $\mathbf{\Omega }$ absorbs centrifugal acceleration into reduced pressure $\Pi =p/{\rho }_0-|\mathbf{\Omega }\times \mathbf{r}{|}^2/2$. Indeed $\mathbf{\Omega }\times (\mathbf{\Omega }\times \mathbf{r})=-\nabla |\mathbf{\Omega }\times \mathbf{r}{|}^2/2$. Transforming the material acceleration in the Euler equations gives the displayed equation. The Coriolis acceleration remains velocity dependent and cannot be absorbed into a scalar potential. A conservative physical body force can be included by adding its potential to $\Pi$.

Consider two dilute populations of spheres with the same buoyant density and radii $b$ and $2b$. Adopt the fixed linear density gradient and constant local centrifugal acceleration of linear-gradient isopycnic focusing, with constant dynamic viscosity and negligible diffusion. Stokes's law makes $V/f$ proportional to $b^2$, so doubling radius multiplies the migration rate parameter by four, although the mass increases eightfold.

Use $u=z/L$, $\tau =t/T$ and $q=({\rho }_b-{\rho }_0)/(GL)=2$. The liquid occupies $0\le u\le 1$, so the matching density position lies beyond its outer boundary. Both populations begin at $u(0)=0.1$; choose ${\lambda }_A=0.15$ and ${\lambda }_B=0.6$. Their band-center trajectories areAt $\tau =0.5$, they lie at approximately $0.2373$ and $0.5924$; at $\tau =1$, at $0.3647$ and $0.9573$. Equal buoyant density has not prevented a useful finite-time separation. Solving $u_B(\tau )=1$ using the natural logarithm gives first boundary arrival at $\tau =\ln (1.9)/0.6\simeq 1.0698$. The corresponding arrival for A is $\ln (1.9)/0.15\simeq 4.2790$. Trajectories describe migration only until contact with the boundary; subsequent collection into a centrifugation pellet is not represented by extending the curve beyond $u=1$.

This is the useful contrast with isopycnic centrifugation: stopping time helps determine the recovered positions in rate-zonal centrifugation, whereas a stationary matching density position sets the limiting band center in isopycnic centrifugation. Real resolution also depends on initial zone width, diffusion and variations in particle properties; distinct calculated centers alone do not prove that fractions will be pure.

Consider a ball moving horizontally with no net horizontal force in the laboratory. Treat its horizontal motion as that of a free particle: ${\mathbf{r}}_{\mathrm{lab}}=(ut,0)$ with constant $u>0$. Axes rotating counterclockwise at constant angular velocity $\mathbf{\Omega }=\Omega {\mathbf{e}}_z$ about the launch point assign it the positionThis is an exact rotation matrix transformation, with the axes coinciding at $t=0$. Taking $u=0.25\mathrm{m}{\mathrm{s}}^{-1}$, $\Omega =0.8\mathrm{rad}{\mathrm{s}}^{-1}$, and $0\le t\le 3\mathrm{s}$ gives the two views below. At $t=3\mathrm{s}$, the laboratory position is $(0.75,0)\mathrm{m}$, but the rotating position is approximately $(-0.5530,-0.5066)\mathrm{m}$. Both assign the same distance $0.75\mathrm{m}$ from the origin.

The entire rotating acceleration must include the centrifugal acceleration as well as the Coriolis acceleration. To see the balance explicitly, write $c=\cos \Omega t$ and $s=\sin \Omega t$. Direct differentiation givesThe two apparent contributions areAt launch $\mathbf{r}=0$, so only the Coriolis acceleration contributes to the initial bending. Later the centrifugal acceleration also contributes; the full curved trajectory cannot be explained by keeping only the Coriolis force. In particular, $|{\mathbf{v}}_{\mathrm{rot}}{|}^2=u^2(1+{\Omega }^2{t}^2)$ grows even though the Coriolis force does no work: the rotating-frame kinetic energy increase is supplied by work from the centrifugal force. The laboratory kinetic energy remains $m{u}^2/2$.

The gravitational constant $G$ sets the strength of Newton's law of universal gravitation. Two point masses $M$ and $m$ a center-to-center distance $r$ apart attract with force magnitudeFor fixed masses and separation, doubling $G$ would double the force. For bodies with spherical symmetry, the same exterior formula follows from the spherical shell theorem; for arbitrary extended bodies, add the attractions of their mass density elements instead. The gravitational constant connects a laboratory attraction to the Newtonian gravitational field controlling a Kepler orbit.

The gravitational constant $G$ and local gravitational acceleration $g$ are different quantities. Applying Newton's second law to a test particle gives $g=F/m=GM/r^2$ for one spherical source. The source's mass and the observer's position change $g$; $G$ remains the same parameter of the law. A local falling-body experiment therefore measures a field, not $G$ alone. On rotating Earth, an instrument fixed to the surface also responds to centrifugal acceleration, while surrounding mass density and altitude change the gravitational contribution.

Dimensional analysis explains the units. Denote the base physical dimensions of length, mass and time by $\mathsf{L},\mathsf{M},\mathsf{T}$. Since force has dimension $\mathsf{ML}{\mathsf{T}}^{-2}$,By contrast, $[g]=\mathsf{L}{\mathsf{T}}^{-2}$ and the standard gravitational parameter $\mu =GM$ has dimension ${\mathsf{L}}^3{\mathsf{T}}^{-2}$. A small numerical value of a dimensional constant depends on the units chosen: it is the resulting force for specified physical masses and separation that conveys the laboratory scale.

NIST's current complete listing givesfor the 2022 CODATA adjustment, the latest adjustment listed when checked on 14 October 2026. Parentheses indicate a standard uncertainty $u(G)=0.00015\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$, approximately $22$ parts per million relative to $G$. The adjustment year is 2022; NIST's May 2024 data update is a different date.^{[ref] [ref]}

For an original scale calculation, stipulate two $1\mathrm{kg}$ point masses separated by $1\mathrm{m}$ and insert this recommended $G$: $F=6.67430\times 10^{-11}\mathrm{N}$. The corresponding acceleration of either mass is $6.67430\times 10^{-11}\mathrm{m}{\mathrm{s}}^{-2}$ at that instant. Reducing the separation to $0.10\mathrm{m}$ multiplies both by $100$. These are ideal model predictions, not observations; their tiny scale explains why stray forces matter.

Being a constant of a physical law does not make the numerical value exact. The SI fixes the numerical values of seven defining constants, including the speed of light and Planck constant; $G$ is determined experimentally. Exact unit definitions still require physical realizations with measurement uncertainty: calibrating a length or a mass does not become error-free.^[ref]

Use local east, north, and up unit vectors, with $x$ eastward, $y$ northward, and positive orientation $\psi$ anticlockwise as viewed from above. Let $\phi$ be signed geographical latitude, positive north, and $\Omega >0$ the magnitude of Earth's angular velocity. In the spherical-Earth approximation,Here $g$ is the local effective gravity in a rotating frame, including the steady terrestrial centrifugal acceleration. For small horizontal displacements compared with $\ell$, the horizontal restoring acceleration from tension is $-{\omega }_0^2(x,y)$. Vertical displacement is second order in the swing amplitude, so vertical velocity can be omitted at this order. The horizontal Coriolis acceleration in local Earth coordinates then gives the coupled linear differential equationsThese equations neglect damping, finite-amplitude corrections, support anisotropy, and spatial changes of the centrifugal force of order ${\Omega }^2$; they retain the leading terrestrial Coriolis force. A related derivation uses south/east axes instead.^[ref]

Introduce the complex number $z=x+iy$. Multiplying the second equation by $i$ and adding givesSet $z=e^{-i\alpha t}u$. Direct differentiation cancels the first derivative of $u$ and yieldsThe ${\alpha }^2$ in $\beta$ comes from solving these reduced equations exactly; it is not a complete prediction of terrestrial corrections of order ${\Omega }^2$. Thus the horizontal trajectory is an ellipse in coordinates rotating at angular velocity $-\alpha$ relative to the building. When $|\alpha |\ll {\omega }_0$, each rapid swing samples an almost fixed long axis; its slow precession isThe sidereal day $T_{\mathrm{sid}}=2\pi /\Omega$ is a rotation relative to distant directions. $T_{\mathrm{p}}$ denotes a full $2\pi$ turn of a continuously followed orientation; an unlabelled swing line repeats after $T_{\mathrm{p}}/2$. At the equator there is no finite precession period in this model. The Coriolis parameter is $f=2\alpha$, so the plane's rate is $-f/2$, not $-f$.

The limiting cases check the sign and scale. At the north pole $\alpha =\Omega$: the swing direction turns clockwise once per sidereal day. At the equator $\alpha =0$: both horizontal equations reduce to independent simple harmonic motions and a line launch stays on its original line. At the south pole $\alpha =-\Omega$: the direction turns anticlockwise once per sidereal day. Vanishing horizontal precession at the equator does not mean that every component of Coriolis acceleration vanishes there.

For a stipulated $\Omega =7.292115\times 10^{-5}{\mathrm{s}}^{-1}$, the sidereal day is approximately $23.9345\mathrm{h}$. At $45^{\circ }$ north the predicted rate is $-10.6356^{\circ }/\mathrm{h}$ and $T_{\mathrm{p}}\simeq 33.8485\mathrm{h}$; at $45^{\circ }$ south the rate changes sign. At the illustrative geographical latitude $34^{\circ }$ south, the prediction is $+8.4109^{\circ }/\mathrm{h}$. These calculations are model predictions, not new observations of the UNSW apparatus.

A swing plane fixed among the stars is generally only a polar picture. A north–south swing at the equator turns with Earth in an inertial frame, while having no precession against the building.^[ref]

Why does an isopycnic centrifugation band return after displacement? Choose a dilute mathematical model with outward coordinate $z$, constant local centrifugal acceleration $a>0$, fixed displaced volume $V>0$, linear friction coefficient $f>0$, and an imposed liquid density ${\rho }_f(z)={\rho }_0+Gz$ with $G>0$. Treat the species' buoyant density ${\rho }_b$ as constant in this medium. The sedimentation coefficient force balance gives the ordinary differential equationIts exponential function solution is $z(t)=z_{\ast }+(z(0)-z_{\ast })e^{-kt}$. The sign of the velocity reverses across $z_{\ast }$, proving stable equilibrium when $z_{\ast }$ lies inside the liquid. Changing $V/f$ changes the approach time $1/k$, while $z_{\ast }$ depends on buoyant density and the imposed density gradient. A real centrifuge rotor has radially varying centrifugal acceleration; the constant-$a$ approximation describes a short radial interval. Here the density gradient is held fixed, rather than modeling its formation.

For an original dimensionless example, let $u=z/L$ and $\tau =t/T$, where $L$ is column length and $T$ a chosen time scale. Write $q=z_{\ast }/L$ and $\lambda =kT$, so $u(\tau )=q+(u(0)-q)e^{-\lambda \tau }$. Species A and B have the same $q=0.45$ but different friction-dependent rates: $({\lambda }_A,u_A(0))=(0.4,0.1)$ and $({\lambda }_B,u_B(0))=(1.6,0.95)$. A moves outward and B inward, both toward $u=0.45$. Species C has $(q_C,{\lambda }_C,u_C(0))=(0.8,0.8,0.1)$ and approaches a separate band at $u=0.8$. At $\tau =4$, their positions are approximately $0.3793$, $0.4508$ and $0.7715$: a finite-time photograph would still show A short of its limiting position.

An equilibrium band has finite width because diffusion opposes focusing. In the same local model, let $c(z)$ be concentration and let the constant diffusion coefficient be $D>0$. Combining the restoring drift with Fick's laws, zero net flux requiresDividing by $c$ and integrating proves the Gaussian function profile, whose standard deviation is ${\sigma }_z=\sqrt{D/k}$. This local result assumes that the band is narrow relative to its distance from the boundaries and that $D$, $k$ and the density gradient vary negligibly across it. Equal buoyant densities therefore imply coincident centers in this model, not identical widths. A band outside the available density range has no such interior restoring point.

Use $D=(d/dt{)}_S$, $D^{\prime }=(d/dt{)}_{S^{\prime }}$, and ${\mathbf{v}}^{\prime }=D^{\prime }\mathbf{r}$. By the rotating-frame derivative formula, $D\mathbf{r}={\mathbf{v}}^{\prime }+\mathbf{\omega }\times \mathbf{r}$. Apply the product rule for the cross product and the same rotating-frame derivative formula again:We have $D{\mathbf{v}}^{\prime }=D^{\prime }{\mathbf{v}}^{\prime }+\mathbf{\omega }\times {\mathbf{v}}^{\prime }$ and $D\mathbf{\omega }=D^{\prime }\mathbf{\omega }$, since $\mathbf{\omega }\times \mathbf{\omega }=0$. ThereforeRearranging this identity into an equation of motion in a rotating frame produces the apparent Coriolis acceleration, Euler acceleration and centrifugal acceleration, with signs opposite to their corresponding terms on the right of this inertial acceleration identity. No translational term is needed because the two origins coincide.

Let $\Omega$ be the Earth's angular velocity magnitude. Use the local east, north and upward unit vectors ${\mathbf{e}}_E,{\mathbf{e}}_N,{\mathbf{e}}_U$, with ${\mathbf{e}}_E\times {\mathbf{e}}_N={\mathbf{e}}_U$. At latitude $\lambda$, the Earth's angular velocity and the train's ground-relative velocity areThus the east component of the Coriolis force isFor a north–south track, the train has no eastward acceleration relative to Earth. Its path curvature lies in the local north–vertical plane. The centrifugal acceleration and Newtonian gravity also lie in that plane, and the Earth's constant angular velocity gives no Euler acceleration. Hence the sideways force from the track must cancel the eastward Coriolis force:Its magnitude is $2m\Omega V\sin \lambda$, and its direction is west in the Northern hemisphere. The eastward tendency is the apparent Coriolis force; the requested track force has the opposite direction.

At $P$, use the local right-handed basis ${\mathbf{e}}_E,{\mathbf{e}}_N,{\mathbf{e}}_U$ pointing east, north and vertically outward. The Earth's rotation vector has componentsInitially $\dot{\mathbf{r}}=0$ in the rotating reference frame, so the Coriolis acceleration is zero. Evaluating the centrifugal acceleration givesThe gravitational acceleration is $-g{\mathbf{e}}_U$. HenceThe small angle $\delta$ from the downward vertical therefore satisfiesFor the northern hemisphere, $0\le \theta \le \pi /2$, this is the printed expression without absolute values. More generally, the printed signed meridional expression specifies the direction of the tilt; an unsigned angle uses its absolute value. The horizontal component is toward the equator and vanishes at both poles and at the equator.

Neglect the centrifugal acceleration, which is of order ${\omega }^2$, and use the local flat-surface approximation. To zeroth order in $\omega$, the height above $P$ and flight time areInserting this vertical velocity into the Coriolis acceleration gives the eastward componentWith $x(0)=\dot{x}(0)=0$, two integrations yieldto first order. Horizontal velocity is already of order $\omega$, so its effect on vertical acceleration is of order ${\omega }^2$. Thus the flight time remains $T_0$ at this order. The Coriolis deflection of a vertical projectile at landing isThe displacement is westward, and its magnitude is $4\omega v^3\sin \theta /(3g^2)$. The horizontal speed remains westward during the flight: $\dot{x}=-2\omega \sin \theta z_0(t)$. This explains why the eastward acceleration on descent does not cancel the earlier displacement.

Apply the rotating-frame derivative formula first to position and then to its inertial velocity:The angular velocity is constant, so there is no term involving $\dot{\mathbf{\omega }}$. By Newton's second law in the inertial frame, the inertial acceleration is $\mathbf{F}/m$. Rearrangement gives the equation of motion in a rotating frame:The two additional accelerations are the Coriolis acceleration and centrifugal acceleration.

Apply the rotating-frame derivative formula first to $r$, then to $v^{\prime }=[dr/dt{]}_{S^{\prime }}$. Since the angular velocity is constant in either frame,The two apparent forces in the rotating equation are the Coriolis force and the centrifugal acceleration multiplied by mass. ThusSet ${\varphi }_{\mathrm{eff}}=\varphi -\frac{m}{2}|\omega \times r{|}^2$. The supplied gradient identity gives $\nabla {\varphi }_{\mathrm{eff}}=\nabla \varphi +m\omega \times (\omega \times r)$. Taking the scalar product with $v^{\prime }$ removes $G$, since it is orthogonal to $v^{\prime }$, and removes the Coriolis force. Therefore the energy balance in a uniformly rotating frame isThe displayed quantity is conserved when the potential is stationary in the rotating coordinates. This includes a static inertial potential invariant under rotations about the chosen axis, such as vertical gravity in part (b). If the word conservative only means a time-independent inertial potential, that extra rotational invariance is needed; a conservative force alone does not imply the claimed conservation.

For an explicit counterexample to that unrestricted interpretation, take unit mass, an inertial force $F=e_x$, $\varphi =-x$, no $G$, rotation $\omega =\omega e_z$, and inertial trajectory $r(t)=(t^2/2,1,0)$ with initial velocity zero. It satisfies Newton's second law. Substituting $v^{\prime }=\dot{r}-\omega \times r$ into the proposed energy gives $E=\omega t$, which is not constant for nonzero rotation. The qualified conservation law above is the one used for the hoop.

In the inertial frame the steady convective acceleration is $-R{\Omega }^2{\mathbf{e}}_R$. For a barotropic fluid, $\nabla p/\rho =\nabla h$, where $h$ is the barotropic enthalpy function. The equilibrium momentum equation therefore giveson each connected fluid region. Uniform rotation is essential: it makes the centrifugal acceleration derivable from the centrifugal potential $-{\Omega }^2{R}^2/2$.

The Cowling approximation sets the gravitational-potential perturbation to zero while retaining background gravity. The barotropic pressure force linearizes as $\delta (\nabla p/\rho )=\nabla \delta h$, with $\delta h=p^{\prime }/\rho =c_s^2{\rho }^{\prime }/\rho$. Denote its mode amplitude by $W$. The advective time derivative on a scalar mode is ${\partial }_t+\Omega {\partial }_{\varphi }$, so the stated positive-frequency exponential produces $i\sigma$ with $\sigma =\omega +m\Omega$. Linearizing the radial centrifugal term produces $-2\Omega v_{\varphi }$, while linearizing the azimuthal convective term produces $2\Omega v_R$. ThusThese are the pressure-gradient and Coriolis acceleration terms in the corotating-frame form. Independently, linearizing continuity equation givesThe equilibrium mass density is axisymmetric, so no additional azimuthal background-density derivative appears.

For $\sigma \ne 0$ and $D=4{\Omega }^2-{\sigma }^2\ne 0$, invert the horizontal momentum system:Inserting these into the continuity equation, dividing by $i\sigma$ and multiplying by $D$ givesThe two terms proportional to $\rho W_R$ cancel; this cancellation leaves precisely the background-density derivative in the last term. Since $-D/{\sigma }^2=1-4{\Omega }^2/{\sigma }^2$, this is exactlyThe singular frequency cases require the original velocity equations rather than this inversion.

In the low-frequency anelastic approximation for a rotating barotropic star, neglect the left side. Spherical background mass density satisfies ${\rho }_R/R={\rho }_z/z$, understood by its smooth limiting form on the equatorial plane. For $W=z{R}^m$,The remaining equation reduces toWriting $q=2\Omega /\sigma$, the bracket factors as $(m+1-q)(q+1)$. The regular desired branch is thereforeThe other algebraic factor corresponds to $\sigma =-2\Omega$, where the eliminated horizontal system is singular; it is not classified by this pressure-equation inversion. A constant-density background makes the bulk factor identically zero but does not invalidate the displayed solution.

For the regular branch, an explicit velocity check is particularly informative. Up to a common mode normalization,It satisfies $R{v}_R+z{v}_z=0$ and $\nabla \cdot \mathbf{v}=0$. Thus $\nabla \cdot (\rho \mathbf{v})=0$ for every spherical mass density profile, directly verifying the anelastic continuity condition. The motion is tangential to spherical shells and is the sectoral inertial mode of a slowly rotating barotropic star. For single-valued azimuthal modes, $m$ is a positive integer. With the source's exponential convention,For $m\ge 2$ the pattern is retrograde relative to the star but prograde in the inertial frame. This is a leading slow-rotation anelastic mode, not an exact solution of the compressible equations whose left-hand term was discarded.