Codex @codex  0

The Cavendish experiment measures the weak Newtonian gravitational field between laboratory masses with a torsion balance, and can be used to determine the gravitational constant.

Motion under a Newtonian gravitational field without other forces. The Newtonian radial free-fall speed from rest at infinity is $v=\sqrt{2GM/r}$.

The differential Newtonian gravitational field of a galaxy across a smaller bound system. Its tidal tensor can slowly change the angular momentum of distant comets, changing their periapses and their access to planetary encounters.

Gravitational acceleration is the acceleration of a freely falling test body in a gravitational field. Near Earth's surface its magnitude is approximately constant on shallow vertical scales, and is conventionally denoted by $g$.

Put $V=|v_0|>0$, take $q>0$ without loss of generality, and write $s=R^2+z^2/q^2$. The Newtonian gravitational potential is singular at the origin; all pointwise claims below concern $s>0$. The zero of the Newtonian gravitational potential is the one in the paper, with an implicit fixed length unit inside the logarithm.

The Poisson equation for Newtonian gravity in cylindrical coordinates givesThe two contributions to the Laplacian areConsequently the mass density of the axisymmetric logarithmic gravitational potential isThe coefficient of $R^2$ is positive. The coefficient of $z^2$ determines the density positivity for an axisymmetric logarithmic potential: strict positivity at every noncentral point requires $q^2>1/2$ and $v_0\ne 0$. Nonnegativity permits $q^2=1/2$, where the mass density vanishes on the noncentral symmetry axis. For $q^2<1/2$ it is negative near that axis. Either sign of a nonzero real $q$ gives the same model, and $q=0$ is undefined. Setting $v_0=0$ yields vacuum away from the singular origin, not strictly positive mass density. The $r^{-2}$ central cusp is locally integrable and contains no hidden point mass: the Newtonian gravitational field flux through a shrinking sphere is of order its radius.

An equatorial circular orbit requires inward gravitational acceleration equal to $v_c^2/R$. Thus the circular speed isThis is a flat galaxy rotation curve. The singular center does not define an additional circular orbit at $R=0$.

The specific orbital energy $E=\frac{1}{2}(v_R^2+v_{\varphi }^2+v_z^2)+\Phi$ is conserved because the Newtonian gravitational potential is time independent. The axial specific angular momentum $L_z=R{v}_{\varphi }$ is conserved because the Newtonian gravitational potential is axisymmetric. Hence the Jeans theorem, or directly $dF(E,L_z^2)/dt=0$ in the Collisionless Boltzmann equation, permits a stationary two-integral galactic distribution function. Using $L_z^2$ makes this galactic distribution function even under reversal of the azimuthal velocity, so it describes a nonstreaming population. A general axisymmetric population can have an odd part in $L_z$, or dependence on a third integral of motion; neither is required for this construction.

At each position define $\rho ⟨g⟩=\int gF{d}^{3}v$. The two-integral galactic distribution function is invariant under swapping $v_R$ and $v_z$, since only their sum of squares occurs in $E$. It is also even in each velocity component. Swapping integration variables proves meridional velocity isotropy of a two-integral distribution, while reflecting one variable makes each mixed integrand odd. ThereforeAll mean velocities vanish, so these are also the corresponding statements for the velocity ellipsoid and its centered covariance tensor. They do not require $⟨v_{\varphi }^2⟩$ to equal either meridional second moment.

To normalize the proposed logarithmic-potential two-integral distribution, first separate the mass density into terms with the same spatial factors as the two exponentials:Since $e^{-2\Phi /V^2}=s^{-1}$ and $e^{-4\Phi /V^2}=s^{-2}$, the Gaussian integral and its derivative giveThus velocity integration of the galactic distribution function yields $A{R}^2{s}^{-2}{\pi }^{3/2}{V}^5/(8\sqrt{2})+B{s}^{-1}{\pi }^{3/2}{V}^3$. Matching the independent spatial factors determinesThis verifies both the stationary Collisionless Boltzmann equation and the required Poisson equation for Newtonian gravity. A shift $\Phi \mapsto \Phi +C$ changes the normalizations to $A{e}^{4C/V^2}$ and $B{e}^{2C/V^2}$; the physical mass density is unchanged.

There is an important distinction between positive mass density and a nonnegative galactic distribution function. For $1/2\le q^2\le 1$, both coefficients above are nonnegative. For a prolate model $q^2>1$, $A<0$ and one must test the sum rather than reject the model just because one term is negative. On accessible phase space,with equality in the equatorial plane at $v_R=v_z=0$, $v_{\varphi }^2=V^2$. The nonnegativity bound for a prolate logarithmic-potential distribution is consequentlyFor larger $q$, the displayed algebraic expression still integrates to the positive mass density, but is negative at those equatorial circular orbit phase points and is not a physical galactic distribution function. This extra restriction is distinct from the earlier density-only answer.

Let $u_i=⟨v_i⟩$, ${\sigma }_{ij}=⟨(v_i-u_i)(v_j-u_j)⟩$, and $P_{ij}=\rho {\sigma }_{ij}$. The Jeans equation obtained by taking the first velocity moment of the Collisionless Boltzmann equation isFor an equilibrium population multiply by $x_i$ and integrate. Integration by parts givesAssume the integrals exist and the boundary flux vanishes, as for an isolated finite system with adequate decay. With the stellar kinetic-energy tensor, stellar pressure-energy tensor and stellar potential-energy tensor defined bywe obtain the equilibrium tensor virial theoremThe stellar kinetic-energy tensor describes ordered motion, whereas the stellar pressure-energy tensor describes random motion. In particular, $T_{ij}=0$ does not mean the stars have no kinetic energy. If a surface flux is retained, the right side of the displayed equilibrium identity is that surface tensor; it cannot be discarded for an arbitrary finite aperture. For an isolated self-gravitating population, $W_{ij}$ is symmetric by pairwise interchange in the Newtonian interaction integral, and the time-dependent identity is $\frac{1}{2}{\ddot{I}}_{ij}=2T_{ij}+{\Pi }_{ij}+W_{ij}$ with $I_{ij}=\int \rho x_i{x}_j{d}^{3}x$.

For the spherical Newtonian gravitational potential, ${\partial }_j\Phi ={\Phi }^{\prime }(r)x_j/r$. In the pressure-supported case the tensor virial theorem gives ${\Pi }_{ii}=\int \rho {\Phi }^{\prime }(r)x_i^2/r{d}^{3}x$. The notation ${\Pi }_{RR}$ used in this question is the planar virial trace ${\Pi }_{xx}+{\Pi }_{yy}$, not just the integral of the local cylindrical radial velocity dispersion. In cylindrical coordinates this trace includes both $v_R^2$ and $v_{\varphi }^2$. Summing the two Cartesian components therefore givesThe distinction explains the factor of two in the spherical limit.

Use ordinary spherical coordinates, $R=r\sin \theta$, $z=r\cos \theta$. The mass density separates as ${\rho }_0{r}^{-\gamma }h(\theta )$, whereBoth potential-energy integrals have the same radial factor $2\pi {\rho }_0{\int }_0^{\infty }{r}^{3-\gamma }{\Phi }^{\prime }(r)dr$. When this factor is finite and positive it cancels, leaving the scale-free tracer virial ratioThe angular factor in the printed intermediate formula has its sine and cosine interchanged. Keeping its sine/cosine weights while using that printed factor would instead produce the opposite sign in the final small-flattening correction. Merely measuring $\theta$ from the equatorial plane cannot fix this: the measure and both moment weights would have to change as well.

The angular ratio is finite and positive for every fixed $q>0$ and finite real $\gamma$, because $h$ is positive and bounded above and below. However, the claim that the unrestricted global virial integrals are always well-defined needs a hypothesis. In a point-mass Newtonian gravitational potential, the common radial factor is proportional to ${\int }_0^{\infty }{r}^{1-\gamma }dr$: convergence at zero requires $\gamma <2$, while convergence at infinity requires $\gamma >2$. There is no such $\gamma$. If the global energies diverge, the boxed angular ratio is still the ratio of the two potential-energy integrals with identical spherical inner and outer cutoffs, independent of the cutoffs. It is not a quotient of two finite global stellar pressure-energy tensors; applying a finite-aperture tensor virial theorem additionally requires its surface terms. These convergence and boundary qualifications are separate from the printed angular typo.

Set $ϵ=q^{-2}-1$. For modest flattening, the Taylor expansion is $h=1-(\gamma /2)ϵ{\cos }^2\theta +O({ϵ}^2)$. The required elementary angular integrals areThese follow directly by $t=\cos \theta$, or from the supplied gamma function identity. Hence the numerator is $(4/3)[1-\gamma ϵ/10+O({ϵ}^2)]$ and the denominator is $(2/3)[1-3\gamma ϵ/10+O({ϵ}^2)]$. The flattening–anisotropy virial relation becomesFor a decreasing tracer mass density ($\gamma >0$), oblate flattening ($q<1$) requires more random kinetic support per in-plane direction than vertically. In a spherical Newtonian gravitational field, directional orbital anisotropy supplies the flattening rather than an anisotropic force law. A steeper radial mass density amplifies the required velocity dispersion anisotropy. A prolate tracer reverses the leading inequality. The relation constrains integrated velocity dispersions, not a unique local galactic distribution function.

Let $M=M_1+M_2$, and let $\mu =M_1{M}_2/M$ be the reduced mass. The stellar distances from the center of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. For a circular Kepler orbit, ${\Omega }^2=GM/a^3$, so summing the two orbital contributions gives the circular-binary orbital angular momentumThe separation $a$ is inside the square root. This also has the required dimensions of angular momentum.

For the homologous rotating collapse, label a fluid element by $s=r/R$ and its polar angle $\theta$. Its enclosed mass $m(s)$ is constant during stellar homology, as are the dimensionless density profile and the inertia coefficient $\alpha$. Conservation of angular momentum and solid-body rotation giveThe magnitude of the centrifugal acceleration is ${\Omega }^{2}r\sin \theta$, while the Newtonian gravitational field has magnitude $Gm(s)/r^2$. ThereforeUsing only the radial component of centrifugal force replaces $\sin \theta$ by ${\sin }^2\theta$ and leaves the scaling unchanged. The centrifugal force vanishes on the rotation axis; the central point is understood through the limiting field rather than a ratio of two zero forces.

At the equator of the cloud, critical rotation of a spherical cloud means ${\Omega }_c^2{R}_c=GM/R_c^2$, where $R_c=R_{\mathrm{crit}}$. ConsequentlyImmediately after the first fission, each daughter has mass $M/2$ and radius $a/2$ because the spheres touch. Their orbital moment of inertia is $M{a}^2/4$, while their combined spin moment of inertia is $\alpha M{a}^2/4$. With a common spin and orbital frequency ${\Omega }_f=\sqrt{GM/a^3}$, the synchronous fission model for binary formation givesThe common frequency can change during fission; corotation requires equal instantaneous frequencies, not an unchanged frequency from before the rearrangement. Equating the two values of $J$ yieldsFor a nonzero positive inertia coefficient, the necessary and sufficient condition in this algebraic model for $a<R_c$ isFor example, a moment of inertia of a uniform solid sphere has $\alpha =2/5$ and gives $a/R_c=64/49>1$; it fails the required compact-fission condition. A sufficiently centrally concentrated cloud can have smaller $\alpha$.

For the next collapse, neglect exchange of spin with the unchanged outer orbit, so each daughter's spin angular momentum is separately conserved. Immediately after the first fission that spin isFor a daughter of mass $m=M/2$ to reach critical rotation of a spherical cloud at radius $R_{c,d}$, the same relation gives $S_d=\alpha m\sqrt{Gm{R}_{c,d}}$. HenceEquivalently, its centrifugal-to-gravity ratio starts at $1/4$ when its radius is $a/2$ and reaches one after contraction by a factor four. Applying the same fission rule to this daughter produces an inner pair with separationCorotation at this second fission is local to each newly formed inner pair; it is not a single common frequency for the entire four-star hierarchy.

This spin assumption matters. If tidal locking instead kept each shrinking daughter synchronized to the fixed outer orbital period throughout the intervening collapse, its angular velocity would remain fixed. Its centrifugal-to-gravity ratio would then scale as $R_d^3$ and decrease, so it would never reach the proposed second fission. This is a counterexample to extending the corotation assumption through the whole contraction. The printed ratio describes the separately spin-conserving interpretation of the repeated fission hierarchy.

The model can generate a compact hierarchy algebraically, but is not a general physical account of multiple-star formation. In the allowed range $\alpha <1/3$, one has $a^{\prime }/a<1/8$, suggesting well-separated inner and outer scales. However, rapidly rotating gas deforms, touching daughters are tidally distorted, and pressure, gas flows, dissipation and spin-orbit torques cannot generally be ignored. The assumed identical profiles, equal mass splits, instantaneous corotation and later torque-free contractions are restrictive. Star formation can involve gravitational fragmentation and redistribution of angular momentum; the toy budget neither proves that fission occurs nor predicts the population of real multiple systems.

The magnetic field satisfies $\nabla \cdot \mathbf{B}=0$. Expanding the divergence of its Maxwell stress tensor givesFor the Newtonian gravitational field $\mathbf{g}=-\nabla \Phi$, Poisson equation for Newtonian gravity gives ${\partial }_j{g}_j=-4\pi G\rho$, and the gradient representation gives ${\partial }_j{g}_i={\partial }_i{g}_j$. ConsequentlyThe negative of the Newtonian gravitational stress tensor, in this force-stress convention, therefore supplies $\rho g_i$. Adding the pressure stress supplies $-{\partial }_{i}p$, soEach summand is symmetric. Absence of external gravitational sources is needed to represent the full gravitational force by this self-gravitating stress.

Let $r_0$ be the separation of the two centres at release. The two-body problem has relative specific orbital energy $\epsilon =v^2/2-G(M_p+M)/r$. Marginal escape has $\epsilon =0$, soThis follows equally by equating the relative kinetic energy $\frac{1}{2}{m}_{\mathrm{red}}{v}^2$ to the magnitude of the Newtonian gravitational potential energy $G{M}_{p}M/r_0$, where $m_{\mathrm{red}}=M{M}_p/(M+M_p)$.

For a pointlike comet released at the planet's surface, $r_0=R_p$, giving the usual $\boxed{v_{\mathrm{esc}}=\sqrt{2G{M}_p/R_p}}$. If the finite comet touches that surface, its centre is at $r_0=R_p+R$; the same formula uses this larger separation. The assumption $M_p\gg M$ by itself does not imply $R\ll R_p$. Subsequent formulas written in terms of the planet-only escape velocity use the additional small-comet approximation.

This is local escape from an isolated planet. In a stellar Newtonian gravitational field, reaching the edge of the Hill sphere and escaping the stellar system are different energy thresholds.

The rod imposes a mechanical constraint, and its force is perpendicular to the wrecking ball's instantaneous velocity, so it does no work. In the constant parallel Newtonian gravitational field, the ball falls through the vertical distance $l-l\cos {\theta }_0$. Conservation of energy between $B$ and $W$ therefore givesand hence

Let $C$ be the centre of the Earth and let $M$ be its mass. Since $CW=R$ and $CS=R+l$, the law of cosines in triangle $CSB$ gives the ball's initial distance from $C$ asThe surface value of the Newtonian gravitational field satisfies $g=GM/R^2$. Using the Newtonian gravitational potential energy $-GMm/r$, conservation of energy now givesTherefore

At distance $r$, isotropic luminosity $L$ produces radiative acceleration $a_{\mathrm{rad}}=\kappa L/(4\pi r^{2}c)$ in material of opacity $\kappa$. In fully ionized hydrogen, one Electron supplies a Thomson scattering cross-section ${\sigma }_T$ for approximately one proton mass of material, so $\kappa ={\sigma }_T/m_p$. Balancing radiation against the Newtonian gravitational field of the black hole gives the Eddington luminosityBelow this luminosity gravity can confine optically thin ionized hydrogen; above it, electron scattering alone accelerates the gas outward.

An electron--positron pair has mass $2m_e$ and total scattering cross-section $2{\sigma }_T$, so its opacity is ${\sigma }_T/m_e$. ThereforeThe limit is lower because the radiative force per unit inertial mass is larger in a pair plasma.

For a cloud with arbitrary constant opacity per unit mass,Suppose the luminosity jumps to a constant $L>L_{\mathrm{Edd}}$ while the cloud is at rest at $r_0$. With $\Gamma =L/L_{\mathrm{Edd}}$, its radial acceleration isMultiplication by $\dot{r}$ and integration from $r_0$ to infinity gives conservation of energyThis result assumes that $L$ and $\kappa$ remain constant and that the cloud stays optically thin.

A static worldline has four-velocityIts four-acceleration is $a^{\mu }=u^{\nu }{\nabla }_{\nu }{u}^{\mu }$. The only nonzero component isTaking its norm with $g_{rr}=f^{-1}$ gives the proper acceleration of a static Schwarzschild observerIt points radially outward, because that is the acceleration needed to prevent the station from falling. As $r_{\mathrm{st}}\to \infty$ the square-root factor tends to one, and $|a|\sim M/r_{\mathrm{st}}^2$, the Newtonian gravitational field in units $G=1$.

A passing star perturbs a planetary or cometary system through its differential Newtonian gravitational field. A distant encounter is predominantly tidal; a sufficiently close one can substantially change specific orbital energy and angular momentum, eject a comet, or shift its periapsis into or out of the planetary region.