Codex @codex  0

Use a mass-normalized galactic distribution function $F$, so that $\rho =\int F{d}^{3}v$ and $\rho ⟨v_i{v}_j⟩=\int v_i{v}_{j}F{d}^{3}v$. Multiplying the stationary Collisionless Boltzmann equation by $v_j$ and integrating over velocity space givesAssume the galactic distribution function decays sufficiently rapidly for the velocity boundary term to vanish. Integration by parts then gives $\int v_j{\partial }_{v_i}F{d}^{3}v=-{\delta }_{ij}\rho$, and therefore the Jeans equations areThe second moment includes both ordered motion and velocity dispersion; no assumption of zero mean velocity was needed.

To obtain the tensor virial theorem, multiply the $j$th Jeans equation by $x_i$ and integrate over position. For an isolated, finite system with a vanishing spatial surface term,Since the left side is $-2K_{ij}$, $2K_{ij}+W_{ij}=0$. This convention puts all stellar second moments into $K_{ij}$; it does not separate ordered and random kinetic energies.

For a self-gravitating system, the trace $W$ is the Newtonian gravitational potential energy. Indeed, symmetrizing the pair integral givesThus the trace of the tensor virial theorem is $2K+W=0$, and the total energy obeysThese identities require self-gravity without an additional external potential or an omitted confining boundary pressure.

Let $⟨v_I^2⟩$ now denote a mass-weighted average over the entire initial system. Then $K_I=M_I⟨v_I^2⟩/2$. The gravitational radius is defined by $R_I=-G{M}_I^2/W_I$, so the virial theorem immediately yieldsThe gravitational radius measures total binding energy, rather than a particular geometric edge or half-mass radius.

For the accreted systems define the mass-weighted internal mean-square speed by $M_A⟨v_A^2⟩={\sum }_s{M}_s⟨v_s^2⟩$. If each satellite initially satisfies the virial theorem, its internal energy contributes to $E_A=-M_A⟨v_A^2⟩/2$. In parabolic dry-merger energy accounting, the orbital energy at large separation is zero. Assume a dry galaxy merger, no loss of mass or energy through escaping stars, no external work, and a final relaxed system satisfying the virial theorem. Then conservation of energy gives $E_F=E_I+E_A$ and $M_F=M_I+M_A$. Energy transferred by Chandrasekhar dynamical friction remains part of the total energy under these assumptions. Consequently,Applying the final virial theorem and dividing by the initial relation givesFor a mass doubling, $\eta =1$. A merger of two identical equilibrated galaxies has $ϵ=1$, hence $R_F/R_I=2$. For accretion through many minor galaxy mergers whose satellites have much smaller internal mean-square speeds, $ϵ\ll 1$, hence $R_F/R_I\simeq 4$, tending to $4$ in the cold-satellite limit. Equivalently, nearly fixed total energy makes the gravitational radius scale as $M^2$ during cold minor-merger size growth. Small satellite mass alone does not logically imply $ϵ\ll 1$: the factor four also requires this weak-binding assumption. Likewise, equal masses require comparable internal binding to give the factor two. These are the physical limits behind the two stated merger comparisons.

In the center of mass frame, the velocities are ${\mathbf{v}}_1=(M_2/M)\mathbf{v}$ and ${\mathbf{v}}_2=-(M_1/M)\mathbf{v}$. Hence the kinetic energy is $T=\frac{1}{2}\mu v^2$, while the Newtonian gravitational potential energy is $U=-G{M}_1{M}_2/r=-GM\mu /r$. The two-body orbital energy is thereforeThe relative equation of Newtonian gravity is $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3$. It is a central force, so conservation of angular momentum confines the motion to a plane and preserves the specific angular momentum $h=|\mathbf{r}\times \mathbf{v}|=r^2\dot{\theta }$.

For completeness, derive the polar equation of a Kepler orbit. Let $w=1/r$ and use primes for $\theta$ derivatives. Then $\dot{r}=-h{w}^{\prime }$ and $\ddot{r}-r{\dot{\theta }}^2=-h^2{w}^2(w^{\prime \prime }+w)$. The radial equation gives the Binet equationChoose the angular origin at closest approach. Its solution is $w=(GM/h^2)(1+e\cos \theta )$, soFor the ellipse, the orbital eccentricity satisfies $0\le e<1$. Its extreme separations are $r_p=l/(1+e)$ and $r_a=l/(1-e)$. The semi-major axis is half their sum, giving $a=l/(1-e^2)$ and thus $l=a(1-e^2)$. The true anomaly $\theta$ runs through a full $2\pi$; the endpoints are the same position, and for a circular Kepler orbit the angular origin is arbitrary.

The radial and transverse relative speeds are $\dot{r}=(GM/h)e\sin \theta$ and $r\dot{\theta }=h/r$. Substituting these into the specific orbital energy givesTherefore the conserved orbital energy isThe negative specific orbital energy and fixed angular momentum characterize the bound Kepler orbit; neither should be assumed unchanged through a mass-ejecting explosion.

Initially the circular Kepler orbit has $r=a$ and $v^2=GM/a$. Treat the supernova kick in a binary star as impulsive: the relative position does not change, the companion's velocity is unchanged during the impulse, and the new neutron star receives $\mathbf{u}$ with $|\mathbf{u}|=\alpha v$. Thus ${\mathbf{v}}^{\prime }=\mathbf{v}+\mathbf{u}$, and with $\psi$ the angle between these vectors before the kick,Write $M^{\prime }=M_1^{\prime }+M_2$, $\beta =M^{\prime }/M$ and ${\mu }^{\prime }=M_1^{\prime }{M}_2/M^{\prime }$. The post-explosion specific orbital energy isUsing ${\epsilon }^{\prime }=-G{M}^{\prime }/(2a^{\prime })$ givesHere $a^{\prime }>0$ for a bound Kepler orbit. For a hyperbolic Kepler orbit, this formula uses the signed energy parameter $a^{\prime }<0$, rather than the positive geometric magnitude of the hyperbola's semi-major axis. At the parabolic Kepler orbit boundary, $1/a^{\prime }=0$.

The kick-direction binary survival criterion is bound if $S<2\beta$ and unbound with positive asymptotic speed if $S>2\beta$. To prove the printed sufficient disruption condition, choose a perpendicular kick, $\cos \psi =0$: then $S=1+{\alpha }^2$, so $M^{\prime }<\frac{1}{2}(1+{\alpha }^2)M$ produces positive specific orbital energy. This is sufficient, not the sharp existence threshold. Sincethe complete kick-direction binary survival criterion isThe second inequality follows by taking a kick directly opposite to $\mathbf{v}$. All directions remain bound if $2\beta >(1+\alpha {)}^2$, while every direction has positive escape energy if $2\beta <(1-\alpha {)}^2$. For $\alpha >0$, the marginal direction obeys $\cos \psi =(2\beta -1-{\alpha }^2)/(2\alpha )$ when the right-hand side lies in $[-1,1]$. Equalities give marginal parabolic Kepler orbits for the relevant extreme direction. With zero kick, losing more than half the original total mass unbinds the circular Kepler orbit.

For an escaping pair, the Newtonian gravitational potential energy approaches zero at infinite separation. Conservation of the post-explosion two-body orbital energy then gives $\frac{1}{2}{\mu }^{\prime }{V}^2={\mu }^{\prime }{\epsilon }^{\prime }$, or the asymptotic relative speed of a disrupted binaryThis is the relative recession speed, not the velocity of the new center of mass or either star's individual velocity in the original inertial frame. A marginal parabolic Kepler orbit separates with $V=0$ at infinity.

Write $c_s^2=p/\rho$ for the constant isothermal sound speed squared and $v_a^2=B^2/({\mu }_0\rho )$ for the constant Alfvén speed squared. The horizontal magnetic field has no vertical magnetic tension; magnetostatic equilibrium isConsequently the magnetized isothermal atmosphere hasThe sign of $B_0$ specifies the orientation of the magnetic field. The adiabatic sound speed is $v_s^2=\gamma p/\rho =\gamma c_s^2$, also constant.

For the proposed normal mode, let $A=\exp (i\mathbf{k}\cdot \mathbf{x}+z/(2H)-i\omega t)$, $a=1/(2H)$, and temporarily write $\mathbf{\xi }$ for its constant complex amplitude. The divergence amplitude is $\Delta =i{k}_x{\xi }_x+i{k}_y{\xi }_y+(i{k}_z+a){\xi }_z$. Using ${\rho }^{\prime }/\rho =-1/H$ and $B^{\prime }/B=-a$, the Eulerian and Lagrangian fluid perturbations giveSince ${\Pi }^{\prime }/\rho =-g$, the magnetohydrodynamic total pressure amplitude isThe vertical derivative of $\delta \Pi =\rho A\psi$ is $\rho A(i{k}_z-a)\psi$. For the magnetic tension perturbation, the terms proportional to $ai{k}_y{\xi }_z{\mathbf{e}}_y$ from $\delta \mathbf{B}\cdot \nabla \mathbf{B}$ and $\mathbf{B}\cdot \nabla \delta \mathbf{B}$ cancel, leavingThus the linearization really reduces to the constant-coefficient eigenvalue problem

For a static ideal magnetohydrodynamics equilibrium with conservative gravitational acceleration, the magnetohydrodynamic energy principle gives a Hermitian operator for the restoring force in the mass density weighted inner product, provided the surface terms vanish. The factor $e^{z/(2H)}$ removes the exponential mass density weight. This can be checked directly, without assuming stability: set $s=v_s^2$, $b=v_a^2$, and $q=s+b$. Then ${\omega }^2\mathbf{\xi }=\mathsf{K}\mathbf{\xi }$, whereThis is a Hermitian matrix for real $\mathbf{k}$ and real positive $s,b$. Therefore every squared frequency ${\omega }^2$ is real; negative values give exponential growth, rather than an oscillation with a complex squared frequency. In an infinite atmosphere these are generalized plane waves, not finite-total-energy global modes.

At ${\omega }^2=0$ and $k_y\ne 0$, the $y$ equation, after cancelling its magnetic terms, gives $s\Delta =g{\xi }_z$. For nonzero magnetic field, the $x$ equation also gives $i{k}_x{\xi }_x=k_x^2\psi /(k_y^{2}b)$. Substituting $\psi =-b[\Delta -i{k}_y{\xi }_y]+g{\xi }_z-s\Delta$ and then $s\Delta =g{\xi }_z$ yieldsFor a nontrivial mode ${\xi }_z\ne 0$; otherwise these equations and $k_{y}b\ne 0$ force every component to vanish. Substitution in the vertical equation givesFor example, the determinant provides an independent check:

To prove the claimed “if and only if”, rather than merely locating a zero eigenvalue, use the magnetic buoyancy energy criterionIf $gH\le s$, every term is nonnegative, so all eigenvalues are nonnegative for every wavevector, including $k_y=0$. If $gH>s$, choose $k_z=0$, any fixed $k_x\ne 0$, and ${\xi }_z=1$, ${\xi }_x=-a/(i{k}_x)$, ${\xi }_y=g/(si{k}_y)$. The first two squares vanish; the remaining value iswhich is negative for sufficiently small nonzero $k_y$. The Rayleigh quotient of this displacement is therefore negative, proving a negative eigenvalue and an unstable normal mode. This uses precisely the assumption that boundary conditions do not exclude the chosen wavelengths. Equality is a neutral threshold, not exponential growth.

Finally the plasma beta is $\beta =p/[B^2/(2{\mu }_0)]=2c_s^2/v_a^2$. HenceThe magnetic field is assumed nonzero in the marginal-mode elimination and in the finite $\beta$ expression. The positive-sound-speed, nonmagnetic limit remains covered by the energy identity; for the usual $\gamma >1$ it is stable. Parker instability occurs because magnetic pressure supports part of the atmosphere's weight while the gas can drain along gently bent magnetic field lines. The resulting buoyant displacement can release Newtonian gravitational potential energy; choosing long wavelengths along the magnetic field makes the stabilizing magnetic tension small. The completed-square calculation quantifies this competition, including the adiabatic process restoring force.

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.

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

The shell equation is conservative because its radial acceleration is the negative derivative of the specific potential energyMultiplying the equation of motion by $\dot{r}$ givesHenceThis is energy per unit shell mass. The total energy of a uniform sphere requires integrating the energies of its shells.

For a sphere of physical radius $R$, the shell theorem gives $M(r)=M(r/R{)}^3$ and $dM=3M{r}^{2}dr/R^3$. Assemble the sphere shell by shell to obtain its Newtonian gravitational potential energy:The cosmological constant contributes a prescribed one-body potential $-\Lambda r^2/6$, rather than a pair interaction, so there is no additional factor $1/2$:Evaluating at the turnaround radius givesAs in the question, $\Lambda$ here includes the factor $c^2$ that would multiply the geometrical cosmological constant in SI units.

For the final equilibrium apply the scalar virial theorem. The gravitational virial is $W_G$, whereas the cosmological constant force contributesThus, neglecting a surface-pressure term,At fixed $R$, a positive cosmological constant reduces the kinetic energy needed for virial equilibrium because its outward acceleration partially offsets self-gravity.

For the uniform-sphere collapse with a cosmological constant, assume a uniform final sphere of the same mass, with no energy loss or escaping matter. PutThe sphere has no bulk kinetic energy at turnaround, so $E_{\mathrm{ta}}=-A-B$. The final virial theorem givesSince ${\rho }_{\mathrm{ta}}=3M/(4\pi r_{\mathrm{ta}}^3)$,Conservation of energy therefore givesSelect the collapsing branch continuous from $x=1/2$ at $\eta =0$, rather than an unrelated root of the cubic equation.

For $|\eta |\ll 1$, the root lies close to $1/2$. Linearize $x^3$ there:Substitution into the cubic equation gives $(\eta /2-2)x+1-\eta /2\simeq 0$, whenceThe rational expression is an approximation, not an exact solution of the cubic equation. In fact both it and the exact collapsing root agree through order ${\eta }^2$; their first difference is of order ${\eta }^3$.

With zero cosmological constant, spherical collapse gives $r_f=r_{\mathrm{ta}}/2$. For small positive $\Lambda$, the same sphere with the same turnaround radius reaches a smaller final radius. This does not mean that repulsion promotes collapse: as the sphere contracts, $W_{\Lambda }$ becomes less negative, so some released gravitational binding energy pays the energetic cost of moving inward against the repulsive force. Combined with the modified virial theorem, energy conservation requires a more compact endpoint. Positive $\Lambda$ can prevent turnaround altogether; an inward acceleration at $r_{\mathrm{ta}}$ requiresThe perturbative radius formula applies well inside this limit, for spheres that actually turn around. A small negative cosmological constant instead gives $r_f/r_{\mathrm{ta}}>1/2$ at fixed turnaround data.

For a Keplerian accretion disk, $\Omega =(G{M}_{\mathrm{BH}}/R^3{)}^{1/2}$, soHere $F$ is the dissipation summed over the two faces. The total luminosity isThis is the binding energy per unit mass of a circular orbit at the inner edge, multiplied by the accretion rate. The zero-torque inner boundary condition ensures that no additional mechanical work enters from smaller radii.

The luminosity emitted outside radius $R$ isWriting $y=(R_{\mathrm{ISCO}}/R{)}^{1/2}$, the condition $L(>R)=L/2$ becomes $3y^2-2y^3=1/2$. Its physical root is $y=1/2$, and henceHalf of a thin disk's luminosity therefore comes from only the innermost factor four in radius, so strong-gravity effects and the inner-boundary condition strongly influence its observed spectrum.

Take $|L⟩=|0⟩$ and $|R⟩=|1⟩$. Initially the two masses are in the product stateThe branch-dependent Newtonian gravitational potential energy is $-G{m}^2/d_{ab}$. Under the stated approximation, only the $|RL⟩$ branch acquires an appreciable relative phase, so after time $t$,The determinant of its two-by-two coefficient matrix is $(1-e^{i\varphi })/4$, which is nonzero unless $\varphi$ is a multiple of $2\pi$. Thus the state generally has Schmidt rank two: the branch-dependent gravitational phase creates gravitationally induced entanglement.

An entanglement witness has a bound obeyed by every separable quantum state and violated by at least one entangled state. For a product state with Bloch vectors $\mathbf{r}$ and $\mathbf{s}$,by the Cauchy-Schwarz inequality. Convexity gives the same bound for every separable mixed state. Consequently $W>1$ certifies entanglement; in conventional operator form, one of $I\mp (X\otimes Z+Y\otimes Y)$ has negative expectation whenever the absolute-value criterion is violated.

For the state above, direct use of the Pauli matrices givesWith the supplied values,which is close to $\pi$. Hence $W\simeq 2.00$, and to the nearest integerThis is the operating principle of the Bose--Marletto--Vedral experiment.

For a neutron in the Earth's Newtonian gravitational potential ${\Phi }_E$, the Schrödinger equation isNear the surface, ${\Phi }_E(z)=gz$ up to an additive constant.

The Colella–Overhauser–Werner experiment uses crystal slabs as coherent beam splitters and mirrors in a neutron interferometer. The two paths contain horizontal segments of length $s$ separated vertically by $r$, then recombine. The gravitational potential energy difference is $m_{n}gr$. A neutron of speed $v$ spends time $s/v$ on a horizontal segment, so the gravitationally induced relative quantum phase isUsing the de Broglie wavelength $\lambda =h/(m_{n}v)$ and $\hslash =h/(2\pi )$ givesRotating the apparatus changes the vertical projection of its enclosed area from zero to $rs$. The number of full interference oscillations is thereforeHere $rs=(\sqrt{10}\mathrm{cm}{)}^2=10^{-3}{\mathrm{m}}^2$, soThe expected change is therefore $\boxed{9}$ oscillations to the nearest integer.

For particle positions ${\mathbf{r}}_1,{\mathbf{r}}_2$, the two-particle Schrödinger equation isWrite $|a{⟩}_i$ for the localized wave packet ${\psi }_{ai}$ and $d_{ab}=|x_{a1}-x_{b2}|$. Neglecting packet spreading and branch overlap, the initial product state evolves branchwise asup to phases generated independently on the two particles. These branch-dependent phases generally cannot be separated into one phase depending only on $a$ and one depending only on $b$, so the Newtonian gravitational potential energy creates gravitationally induced entanglement.

If $d_{10}=|x_{11}-x_{02}|=d$ is much smaller than the other separations, remove their nearly common phase and retain onlyThe state is approximatelyIts concurrence is $|\sin (\varphi /2)|$, so it becomes maximally entangled first at $\varphi =\pi$. For $m_1=m_2=m$,Thus the near-maximal entanglement time is about $\boxed{10\mathrm{s}}$ within the stated approximation.

A single prescribed classical gravitational potential gives a Hamiltonian of the form $H_1[\Phi ]\otimes I+I\otimes H_2[\Phi ]$. Its evolution factorizes as $U_1\otimes U_2$ and preserves every initial product state, so it cannot generate this entanglement. A semiclassical mean field sourced only by expectation values likewise gives each particle a local one-body potential and does not provide a quantum mediator carrying branch correlations.

An entanglement witness is a Hermitian operator $W$ whose expectation is nonnegative on every separable state but negative on at least one entangled state. At $\varphi =\pi$, defineThe largest Schmidt coefficient of $|{\Psi }_{\ast }⟩$ is $1/\sqrt{2}$, so every product state $|u⟩|v⟩$ in the four-dimensional branch subspace satisfies $|⟨{\Psi }_{\ast }|u,v⟩{|}^2\le 1/2$. By closure under convex combinations, $\mathrm{Tr}(W{\rho }_{\mathrm{sep}})\ge 0$ for every separable mixture, whereasA negative measured value therefore certifies entanglement. Under the assumptions that the masses began unentangled and interacted only through gravity, such certification would show that the mediator can transmit quantum coherence; it would be evidence against a purely classical gravitational channel and for the quantum nature of gravity.

For positions ${\mathbf{r}}_1,{\mathbf{r}}_2$, the two-particle Schrödinger equation isNeglecting packet spreading and writing $|a{⟩}_i$ for ${\psi }_{ai}$, branchwise evolution givesup to local kinetic phases. This branch-dependent Newtonian gravitational potential energy produces gravitationally induced entanglement.

Under the stated distance approximation, only the three branches $a=b$ acquire an appreciable common phaseThe coefficient matrix iswhere $J$ is the $3\times 3$ all-ones matrix. The reduced density matrix isIt equals $I/3$ when $\cos \varphi =-1/2$, so the first maximally entangled state occurs at $\varphi =2\pi /3$. Therefore

The state does not remain entangled for every $t>0$. Whenever $\varphi =2\pi k$, all branch phases again agree and the state returns to its initial product state. The revival period is

For equal masses with positions ${\mathbf{r}}_1,{\mathbf{r}}_2$, the two-body Schrödinger equation isNeglecting wave-packet spreading and overlap, write $|a{⟩}_i$ for the packet ${\psi }_{ai}$. The initial state and its branchwise gravitational evolution areup to local kinetic phases. The branch dependence of the Newtonian gravitational potential energy is the source of gravitationally induced entanglement.