Codex @codex  0

For $0\le k<1$, the first burn reduces the total mass from $M+m_0$ to $M+(1-k)m_0$. The rocket equation therefore givesDetachment exerts negligible impulse, so the surviving second stage keeps velocity $v_1$. Its initial and final masses are $(1-k)(M+m_0)$ and $(1-k)M$, respectively. A second application of the rocket equation givesFor the single-stage comparison from rest, $v_{\mathrm{single}}=u\log (1+m_0/M)$. HenceEquality holds at $k=0$; for positive fuel and exhaust speed it is strict when $0<k<1$. This two-stage rocket with proportional fuel and body masses gains speed by discarding inert body mass before the last burn.

As $k\to 1^{-}$ the expression tends to $2u\log (1+m_0/M)$, but the surviving second-stage mass also tends to zero. At exactly $k=1$ there is no second stage, so its final speed is undefined and that limiting formula cannot be assigned to it. The first stage alone reaches the single-stage burnout speed before detachment. If $m_0=0$, neither design gains speed.

Use $v=|u_r|>0$ for the radial speed and distinguish the central mass $M$ from the Mach number $\mathcal{M}=v/c$. In each smooth nonzero-flow region, mass conservation, the adiabatic equation of state, and the Bernoulli equation giveHere $K>0$ is constant on a smooth isentropic flow branch, $c$ is the adiabatic sound speed, and $\mathcal{E}$ is the Bernoulli function. Either sign of $u_r$ is possible. The positive-radius nondimensionalization in the question requires $\mathcal{E}>0$; flows with $\mathcal{E}\le 0$ do not have these positive $r_s,x,\lambda$ and require another scaling.

Indeed, the supplied adiabatic process equations give $D(p{\rho }^{-3/2})/Dt=0$, so $K$ is spatially constant along a steady nonzero radial fluid flow. The adiabatic sound speed is $c^2=(\partial p/\partial \rho {)}_s$, and the specific enthalpy is $h=c^2/(\gamma -1)=2c^2$. The radial Euler momentum equation, $vdv/dr=-{\rho }^{-1}dp/dr-GM/r^2$, integrates to $\mathcal{E}$ because $dh=dp/\rho$ on this isentropic flow branch.

Put $D=[j(3K/2{)}^2{]}^{2/5}$. Eliminating the mass density from the continuity equation givesWith $y={\mathcal{M}}^{2/5}$, the Bernoulli equation becomesTherefore the physical constants of the spherical polytropic flow with adiabatic exponent three halves are

For $t>0$, $F^{\prime }(t)=4(t^5-1)/t^2$, so $F$ decreases from infinity to its unique minimum $F(1)=5$, then increases to infinity. This completely determines the algebraic curves. For $0<\lambda <1$, there are two disconnected branches at every $x>0$, one subsonic flow with $y<1$ and one supersonic flow with $y>1$. For $\lambda =1$, the line $y=x$ and a decreasing branch meet at $(1,1)$. Their slopes are $+1$ and $-1$, since the leading expansion is $F(1+\eta )=5+10{\eta }^2+O({\eta }^3)$. For $\lambda >1$, solutions exist only where $F(x)\ge 5\lambda$: if $x_{-}<1<x_{+}$ solve $F(x_{\pm })=5\lambda$, there is a forbidden interval $x_{-}<x<x_{+}$. At either endpoint the two branches turn at $y=1$ and have an infinite slope.

A smooth sonic point needs $F^{\prime }(x)=F^{\prime }(y)=0$; hence a regular sonic transition requires $\lambda =1$ and $r=r_s$. The critical speed of a polytropic flow is then $c_s^2=GM/(2r_s)=2\mathcal{E}$, and the allowed mass flux is $j=(2\mathcal{E}{)}^{5/2}{r}_s^2/(3K/2{)}^2$. Here the subscript on $c_s$ denotes evaluation at the sonic point, rather than the isothermal sound speed convention used in the next solution.

The simple transonic branch is $y=x$. Direct substitution givesThus the radial velocity has constant magnitude, while the mass density and pressure scale as $r^{-2}$ and $r^{-3}$. For an outward fluid flow, this is a subsonic flow becoming a supersonic flow as its adiabatic sound speed decreases.

On the other transonic branch, $y$ decreases as $x$ increases. At large $r$, the terms $4/y$ and $x^4$ dominate, giving $y\sim 4x^{-4}$. At small $r$, the terms $y^4$ and $4/x$ dominate, giving $y\sim 4^{1/4}{x}^{-1/4}$. The resulting asymptotic expansions areFor inward fluid flow, these are the nearly static reservoir at infinity and the inner free fall region of Bondi accretion. The specific enthalpy is smaller than $GM/r$ in the inner limit, which explains free fall. Reversing the radial velocity gives a decelerating outflow with the same profiles. The formal small-$r$ limit describes the point-mass model: if the spherical gravitating body has a finite surface, the exterior solution stops there.

A stationary, nonradiative normal shock wave conserves the mass flux and the sum $v^2/2+h$, while the Newtonian gravitational potential is continuous across its thin layer. Hence $\mathcal{E}$ and $r_s$ are unchanged. The entropy production in a perfect-gas shock increases $K=p/{\rho }^{3/2}$, soMore explicitly, the Rankine-Hugoniot condition for $\gamma =3/2$ gives, with upstream $m={\mathcal{M}}_1^2>1$,Thus the shock wave jumps vertically, at the same $x$, from the supersonic flow to a subsonic flow on a larger-$\lambda$ curve. The figure illustrates an outward flow on $y=x$ with a shock wave at $x=1.4$; the downstream subsonic flow has $\lambda \simeq 1.224$ and continues to larger $x$. The shock wave radius is an additional boundary-data choice, rather than being determined by the smooth equations alone.

The stellar relaxation time is the time on which accumulated discrete gravitational encounters change a typical star's velocity by an amount comparable to its original velocity. It concerns two-body relaxation, rather than the much faster orbital evolution in a smooth gravitational potential. A collisionless stellar system requires this time to greatly exceed the duration being studied.

In the straight-line impulse approximation, an encounter at relative speed $u$ has transverse acceleration $Gmb/(b^2+u^2{t}^2{)}^{3/2}$. ThusA small-angle gravitational encounter needs $\Delta v_{\perp }\ll u$. The transition to order-one deflections is therefore $b\sim Gm/u^2$, ignoring factors such as the equal-mass relative deflection. For the spherical estimate $u\sim v$,This is the lower cutoff of the weak-scattering estimate; closer encounters actually occur and require strong-scattering treatment.

The number of encounters in time $T$ with impact parameters in $[b,b+db]$ is $2\pi bdbnuT$. Independent random transverse directions make mean kicks cancel while their variances add. ConsequentlyEach logarithmic interval contributes equally, giving the Coulomb logarithm in stellar dynamics. Defining a deflection time by variance $v^2$ and setting $u=v$ gives $v^3/(8\pi G^2{m}^{2}n\log \Lambda )$. The paper instead uses an order-one normalization three quarters of this estimate:Both have the same physical scaling. The specified characteristic speed, the word “comparable”, the velocity-distribution average and strong-encounter cutoff do not fix that numerical coefficient uniquely. The crude equal-speed impulse calculation must not be claimed to determine $3/(32\pi )$ exactly.

For a self-gravitating, approximately virialized system, the virial theorem gives $v^2\sim GNm/R$. With $n=3N/(4\pi R^3)$, $\Lambda \sim N$ and the stellar crossing time $t_{\mathrm{dyn}}=R/v$, substitution into the paper's convention givesAn externally dominated gravitational potential or a different structural constant changes this substitution. For an N-body simulation lasting $T=K{t}_{\mathrm{dyn}}$, a tolerable fractional velocity-squared diffusion $\epsilon \ll 1$ requires $T/t_{\mathrm{rel}}\lesssim \epsilon$, henceThere is no unique smallest $N$ without $K$, the error tolerance and the force prescription. For scale, equality at $t_{\mathrm{rel}}/t_{\mathrm{dyn}}=100$ occurs near $N=7.1\times 10^3$; $N=10^5$ gives about $1.1\times 10^3$ crossing times. A calculation lasting 100 crossing times therefore needs substantially more than the first threshold to have negligible relaxation. Gravitational softening can raise the effective cutoff and reduce artificial scattering, but it also sets the spatial force resolution.

Put $u=\lambda v$ and let $n$ now mean number per area. During time $T$, the two strips of impact parameters between $b$ and $b+db$, on opposite sides of the trajectory, have combined area $2uTdb$. ThusThe same straight-line impulse approximation gives a kick $2Gm/(bu)$. Summing its squared size over these encounters givesassuming $b_{\max }\gg b_{\min }$. The relevant change is now the small random speed $u$, not the full circular speed $v$. Setting $D_v{t}_{\mathrm{rel}}\sim u^2$ yieldsThe razor-thin disk relaxation integral is dominated by nearby encounters, rather than equally by logarithmic intervals. Because kicks cease to be weak near the cutoff, its numerical coefficient is an estimate.

Here the weak-deflection cutoff must use the relative speed: $b_{\min }\sim Gm/u^2=Gm/({\lambda }^2{v}^2)$, which is larger than the spherical cutoff. Thus $t_{\mathrm{rel}}\sim u/(Gmn)$. Using $n\sim N/R^2$ and $v^2\sim GNm/R$ from the virial theorem givesFor example $n=N/(\pi R^2)$ gives the order-one coefficient $\pi /8$ in this last substitution; the paper's approximate equality suppresses such structural factors. The cancellation of $N$ is the central result.

An unsoftened point-particle system confined to an infinitesimally thin plane is therefore not collisionless over orbital times for $\lambda \sim 0.1$ in this encounter model. At fixed total mass, increasing $N$ decreases $m$ and the cutoff together; the enhanced importance of close planar encounters cancels the usual benefit of increasing particle number. This is a singular zero-thickness limit. A fixed nonzero gravitational softening length or physical thickness restores a different $N$ dependence, so the conclusion is not a prohibition on collisionless numerical models of disks.

For finite-thickness disk relaxation, assume a well-separated range $b_{\min }\ll z_0\ll R$ and approximately virialized self-gravity. The three factors have distinct origins. Replacing the characteristic encounter speed by $u=\lambda v$ shortens the diffusion time by ${\lambda }^3$: the variance-production rate has one inverse power of speed, and the target random velocity squared has two powers. Concentrating the same number of stars into volume of order $R^2{z}_0$, rather than $R^3$, raises the number volume density by $R/z_0$ and shortens the time by $z_0/R$. Finally, only separations below the thickness sample a three-dimensional encounter geometry. The logarithmic upper cutoff becomes $z_0$, and therefore the inverse-logarithm time acquires $\log (R/b_{\min })/\log (z_0/b_{\min })$. More distant encounters have the planar geometry and provide a convergent, non-logarithmic correction. A smaller Coulomb logarithm in stellar dynamics actually partly offsets the first two shortenings.

Thus, using the same cutoff for this comparison,For actual point particles the disk cutoff is $Gm/(\lambda v{)}^2$, whereas the spherical one is $Gm/v^2$; one must use the appropriate cutoff in each logarithm rather than silently identifying them. For a softened calculation it is instead of order the larger of the weak-deflection scale and the softening length. The formula is not valid when $z_0$ approaches that cutoff.

Write $h=z_0/R$ and $q=K/\epsilon$, with simulation duration $K{t}_{\mathrm{dyn}}$ and allowed fractional diffusion $\epsilon$. Combining the paper's spherical normalization and its comparison formula gives the leading criterionThe spherical comparison logarithm cancels when its numerical convention and a common comparison cutoff are used; order-one shape factors remain approximate. In an unsoftened disk with $v^2\sim GNm/R$, $z_0/b_{\min }\sim h{\lambda }^{2}N$. Hence the implicit requirement is $N/\log (h{\lambda }^{2}N)\gtrsim 8q/({\lambda }^{3}h)$, in the regime $h{\lambda }^{2}N\gg 1$. For $\lambda =h=0.1$, the stellar relaxation time is of order $10^{-4}$ of the spherical value apart from logarithms; $N\sim 10^8$ gives a time of order 100 crossing times in this estimate. Negligible relaxation over that duration requires an additional margin. A disk generally needs far more particles than a round galaxy for the same relaxation tolerance; no universal minimum follows without duration, thickness and softening. Collective spiral or bending responses are a separate source of evolution even in a collisionless stellar system.

Use the same specified homogeneous Newtonian background and Jeans swindle, but describe stars by a mass-normalized galactic distribution function. The Collisionless Boltzmann equation is ${\partial }_{t}f+\mathbf{v}\cdot {\nabla }_{x}f-\nabla \varphi \cdot {\nabla }_{v}f=0$. Let $f=f_0(\mathbf{v})+f_1{e}^{i(\mathbf{k}\cdot \mathbf{x}-\omega t)}$ and similarly perturb the gravitational potential. To first order,Eliminating $f_1$ and the nonzero gravitational potential amplitude gives the collisionless Jeans dispersion relationFor growing modes $\mathrm{Im}\omega >0$ there is no real-velocity pole. Real and damped frequencies require a causal contour prescription or analytic continuation; one must not simply integrate through a pole without specifying it.

Take $b>0$ in the Cauchy velocity distribution. Its normalization is $\int f_0{d}^{3}v={\rho }_0$. Orient $\mathbf{k}$ along the $u=v_z$ axis and integrate over the perpendicular velocities:At the marginal mode $\omega =0$, $F^{\prime }(u)/u=-2{\rho }_{0}b/[\pi (u^2+b^2{)}^2]$, with its finite limiting value at zero. The provided Gamma function integral at power two givesConsequentlyTo check that this really separates growing modes, put $c=\omega /k$ in the upper half-plane. Residue integration, or integration by parts followed by the Cauchy resolvent, gives $\int F^{\prime }(u)/(u-c)du={\rho }_0/(c+ib{)}^2$. The upper-half-plane solution has $\omega =i(\sqrt{4\pi G{\rho }_0}-kb)$, which grows exactly for $0<k<k_J$. At the threshold it approaches the marginal mode continuously. Although $b$ plays the threshold role of a velocity scale, this distribution's second moment diverges: the radial integrand for $\int v^2{f}_0{d}^{3}v$ tends to a nonzero constant at large speed. It is therefore incorrect to identify $b^2$ with a finite velocity dispersion or to infer this result by substituting an rms speed into a fluid formula.

Let $q=\sigma /L$, with $L>0$, and let $Z=X^1+i{X}^2$. Direct differentiation gives $\ddot{Z}=Z^{\prime \prime }=-Z/L^2$, while the time coordinate also satisfies the wave equation. Using the Minkowski metric,Thus both Virasoro constraints hold, and the induced worldsheet metric is $g_{\mu \nu }={\cos }^{2}q\mathrm{diag}(-1,1)$. The configuration solves the Nambu–Goto action equations in its interior.

Place the fixed end at $\sigma =0$. The first free endpoint has $Z^{\prime }=e^{it/L}\cos q=0$ at $q=\pi /2$. For the usual single unfolded rotating Nambu–Goto string,This parameter interval need not equal the earlier interval $[0,\pi ]$, since that interval was a coordinate convention. The proper length is measured on a constant-$X^0$ slice; the local motion is perpendicular to the string, so there is no longitudinal Lorentz contraction.

All points have the same phase $t/L$, hence the segment rotates rigidly with angular velocity $\Omega =1/L$. At radius $r=L\sin q$, its speed is $v=r/L$. Therefore the free tip moves at the speed of light. This is consistent with its Neumann boundary condition: $X^{\prime }=0$ and the Virasoro constraints imply ${\dot{X}}^2=0$ at the free tip. The null endpoint is a boundary limit, not a nondegenerate interior point.

In this conformal gauge, $P_m=T{\dot{X}}_m$, so the energy density per $\sigma$ is $-P_0=T$. Equivalently, the energy density per proper length is $T/\sqrt{1-r^2/L^2}$. Thusand the rotational kinetic energy isThe angular momentum about the fixed point isTherefore the Regge trajectory isThe endpoint conditions alone also admit folded extensions with parameter length $(n+\frac{1}{2})\pi L$, $n\ge 0$. Their proper length, energy and angular momentum are $(2n+1)$ times the values above, and $\beta =1/[(2n+1)\pi T]$. The displayed answer uses the implicit unfolded-segment convention.

Write $E=|\mathbf{E}|$ and $\hat{\mathbf{E}}=\mathbf{E}/E$. The constant electric force gives momentumCombining $\mathbf{p}=\gamma m\dot{\mathbf{x}}$ with the relativistic energy-momentum relation givesThe dimensional analysis from part (a) confirms that $qEt$ has the dimensions of momentum, so $qEt/(mc)$ is dimensionless and the expression has the dimensions of velocity.

The speed increases monotonically from zero and approaches the speed of light asymptotically:Integrating in the direction of the field givesAt small times, the Taylor expansion of the square root gives the nonrelativistic limitwhich is motion from rest with the constant acceleration $q\mathbf{E}/m$.

Writing $x_i=(x_i,y_i,z_i)$, rotation about the $z$-axis gives particle $i$ the speed $\omega \sqrt{x_i^2+y_i^2}$. Its kinetic energy is therefore $m_i{\omega }^2(x_i^2+y_i^2)/2$, soThis is the defining mass sum for the moment of inertia about the $z$-axis.

For the uniform cuboid, the mass density is $\rho =M/(8abc)$ and

Comparing the correction in part c with the nonrelativistic kinetic energy $m_{\nu }⟨v^2⟩/2$ gives the characteristic late-time speedThus the late-time relic-neutrino speed redshifts as $a^{-1}$.

For a stationary collisionless stellar system confined to a plane, velocity isotropy makes the galactic distribution function depend on speed through $u=v^2/2$. The Collisionless Boltzmann equation implies ${\nabla }_{x}f=f_u{\nabla }_x\Phi$, hence $f=F(E)$ with specific orbital energy $E=u+\Phi$. This conclusion requires stationarity as well as isotropy.

In an unsoftened plane with number surface density $n_2$ and relative speed $u$, an encounter strip has area $2udtdb$. Thus the velocity diffusion rate is $8G^2{m}^2{n}_2/u\int db/b^2$, dominated by the small-impact cutoff. Setting this cutoff to $Gm/u^2$ and requiring a change of order $u^2$ gives $t_{\mathrm{rel}}\sim u/(Gm{n}_2)$. With $u=\lambda v$, $n_2\sim N/R^2$ and $v^2\sim GNm/R$, this is of order $\lambda R/v$, independent of $N$ at fixed total mass. The point-particle, zero-thickness limit is singular; fixed softening or nonzero thickness changes it. Numerical coefficients are not controlled by extending weak kicks down to order-one deflections.

For a mass $m$ passing at relative speed $u$ and impact parameter $b$, integration of $Gmb/(b^2+u^2{t}^2{)}^{3/2}$ gives the transverse kick $2Gm/(bu)$. Weak deflection requires $b\gg Gm/u^2$, up to the equal-mass order-one factor. This scale supplies a small-impact cutoff for the weak-scattering diffusion approximation, not a physical exclusion of stronger encounters.

The stellar relaxation time is the time on which random gravitational encounters change a star's velocity by an amount comparable to its characteristic random velocity. A weak straight-line encounter at relative speed $u$ and impact parameter $b$ gives a kick $2Gm/(bu)$. Three-dimensional encounter counting gives $d⟨|\Delta \mathbf{v}{|}^2⟩/dt=8\pi G^2{m}^{2}n\log (b_{\max }/b_{\min })/u$. Definitions of relaxation, velocity averaging and system structure alter order-one coefficients; the robust spherical scaling is $t_{\mathrm{rel}}\sim v^3/(G^2{m}^{2}n\log \Lambda )\sim (N/\log N)t_{\mathrm{dyn}}$. Collisionless simulation requirements depend on the duration and permitted fractional relaxation, not a universal minimum particle number.

Two equal uniform rods of length $2l$ and mass $M$, one pivoted to the floor and the other endpoint sliding on it, have arch coordinates $B=(2l\cos \theta ,2l\sin \theta )$ and $C=(4l\cos \theta ,0)$. The kinetic energy decomposition about the center of mass gives $T=\frac{4}{3}M{l}^2(1+3{\sin }^2\theta ){\dot{\theta }}^2$, while gravitational potential energy is $2Mgl\sin \theta$. Conservation of energy from rest at angle $\alpha$ gives the displayed equation and hinge impact speed $\sqrt{6gl\sin \alpha }$, provided contact is maintained. The endpoint's normal reaction becomes negative during this constrained fall if $\sin \alpha >35/54$; a unilateral floor then permits detachment and the constrained formula no longer describes the entire motion.