Codex @codex  0

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 $K={\rho }_0|\mathbf{u}{|}^2/2$ be the kinetic energy per unit volume. Dot the momentum equation with ${\rho }_0\mathbf{u}$. The Coriolis acceleration does no work because $\mathbf{u}\cdot ({\mathbf{e}}_z\times \mathbf{u})=0$. Incompressibility turns the advection and pressure terms into divergences, yieldingThus $\mathbf{F}=(K+P)\mathbf{u}$ is the energy flux and $S={\rho }_0{\Omega }^2(3x-2qz)u_x$ is the source in the requested kinetic-energy form.

To decide about conservation of energy, absorb the ordinary radial shearing-sheet tidal potential ${\Phi }_t=-3{\Omega }^2{x}^2/2$ into the density. Its material derivative cancels the $3x$ work term:For $q=0$ this is the usual conservative Jacobi energy in a shearing sheet, subject to vanishing boundary energy flux. For $q\ne 0$, the extra body force has nonzero curl, $\nabla \times [-2{\Omega }^{2}qz{\mathbf{e}}_x]=-2{\Omega }^{2}q{\mathbf{e}}_y$, so no scalar potential can absorb it. A potential $2{\Omega }^{2}qxz$ would necessarily introduce a vertical force $-2{\Omega }^{2}qx{\mathbf{e}}_z$, absent from the model. Therefore this vertical-shear model does not generically conserve the standard total mechanical energy: its maintained shear/irradiation background can supply or remove energy. This is the energy balance of the incompressible vertical-shear model.

For the equal-width Hazel model with stable density stratification $J\ge 0$,The Miles–Howard theorem therefore excludes exponentially growing regular normal modes for $J\ge 1/4$. The candidate neutral modes of the equal-width Hazel model instead lie onThe parabola touches the sufficient-stability threshold at its maximum. There is no contradiction: neutrality has $c_i=0$, whereas the proof in part (a) assumes $c_i>0$. Moreover, the $k=1/2$ profile has a singular critical-layer derivative and logarithmically divergent horizontal kinetic energy; it is not a regular growing mode satisfying the proof's hypotheses.

For $J<1/4$, the local gradient Richardson number condition permits instability but does not establish it. Substitution of a neutral ansatz alone also does not determine on which side of the curve unstable eigenvalues lie. It identifies the formal neutral curve; concluding a full stability boundary requires additional continuation analysis of the eigenvalues. The nondecaying $k=0$ endpoint and the smooth $k=1,J=0$ endpoint have the distinct qualifications described above.

Set $y=\tanh z$ and $a=1-k$, so $dy/dz=1-y^2$ and $\hat{w}=(1-y^2{)}^{k/2}{y}^a$. Away from $y=0$, the logarithmic derivative givesFor $c=0$, ${\overline{U}}^{\prime \prime }/\overline{U}=-2(1-y^2)$ and $N^2/{\overline{U}}^2=J(1-y^2)/y^2$. The Taylor–Goldstein equation residual divided by $\hat{w}$ consequently simplifies toThus the neutral mode of the equal-width Hazel model requiresThe qualifier “solution” needs care at the critical level of an internal gravity wave $z=0$. For $0<k<1$, $\hat{w}\sim z^{1-k}$ and ${\hat{w}}^{\prime }\sim (1-k)z^{-k}$: it is a solution separately on either side, not a classical smooth eigenfunction across the critical level of an internal gravity wave. A branch for negative $\tanh z$ is also needed. For example, a limit from $c_i>0$ assigns $(\tanh z{)}^{1-k}=|\tanh z{|}^{1-k}{e}^{-i\pi (1-k)}$ on $z<0$.

The critical-layer regularity of a neutral Hazel mode further gives local finite horizontal kinetic energy only for $0<k<1/2$, since incompressible flow gives $\hat{u}=i{\hat{w}}^{\prime }/k$ and ${\int }_0^{ϵ}{z}^{-2k}dz$ then converges. At $k=1$, $J=0$ and $\hat{w}=\mathrm{sech}z$ is smooth after the removable ${\overline{U}}^{\prime \prime }/\overline{U}$ quotient is continued. At $k=0$, $J=0$ and the formal $\tanh z$ profile does not decay at infinity, so it fails the remote endpoint condition. These qualifications prevent interpreting the whole closed parameter interval as a family of classical decaying modes.

Use $\nu =R{e}^{-1}>0$, $\kappa =P{e}^{-1}>0$, $\mathbf{U}=\overline{\mathbf{u}}+\mathbf{u}$, and $\Theta =\overline{\theta }+\theta$. Work along a sufficiently regular direct trajectory satisfying the given constraints. The zero spatial mean of the scalar field is conserved by incompressible flow, impermeable walls and zero scalar flux, so minimizing $J=(\Theta (T),\Theta (T))$ is equivalent to minimizing scalar variance, up to the fixed domain volume.

The printed functional fixes the initial state to a candidate ${\mathbf{u}}_0$; it contains no term that enforces its kinetic energy. For the optimization over that candidate, add the real Lagrange multiplier constraintEquivalently, one can restrict all control variations to the sphere in a normed vector space of fixed kinetic energy. This term changes the initial-control optimality condition, not the interior adjoint equations.

Let $\mathbf{v}=\delta \mathbf{u}$, $\sigma =\delta \theta$ and $\pi =\delta p$. Linearization of the momentum and scalar transport residuals givesBoth appearances of the perturbation velocity in the nonlinear momentum term have been differentiated. In particular, the coefficient is the gradient of the total velocity, not just the base shear.

For the negative-constraint convention of the functional, integration by parts gives the interior coefficients of $\mathbf{v},\sigma ,\pi$ asThus the adjoint equations for Boussinesq scalar mixing areThe transpose is essential: the $j$th component of $(\nabla \mathbf{U}{)}^T{\mathbf{u}}^{†}$ is ${\sum }_i{u}_i^{†}{\partial }_j{U}_i$. The coupling $-{\theta }^{†}\nabla \Theta$ transposes advection of the scalar by a velocity perturbation; $-R{i}_B{u}_y^{†}$ transposes buoyancy feedback. Dropping the latter would give a passive scalar adjoint, not the active scalar problem.

These equations are integrated backward, not forward. If $\tau =T-t$, they readwith direct coefficients evaluated at $t=T-\tau$. Both terms from the diffusion equation now have the usual forward sign in $\tau$. A direct-adjoint looping method stores or reconstructs the forward trajectory, solves these equations backward, and uses the initial adjoint as the control gradient. The endpoint and fixed-energy conditions below give necessary conditions for a local optimizer, not a global optimality theorem.

The temporal integration by parts terms areThe terminal states are free and the objective has no explicit terminal-velocity dependence. ThereforeThe factor two follows from the displayed objective without a $1/2$ prefactor. The initial scalar field perturbation is fixed, so $\delta \theta (0)=0$ and there is no independently prescribed adjoint-scalar condition at $t=0$. Its value there is obtained by backward integration. Independent variation of the initial velocity state gives ${\mathbf{u}}^{†}(0)={\mathbf{u}}_0^{†}$.

Variation of the candidate ${\mathbf{u}}_0$ in the energy-augmented functional gives the fixed-energy initial-condition optimality condition. For the usual divergence-free no-slip control space $H$, write $P_H$ for its orthogonal projection, the Leray-Helmholtz projection in the usual incompressible velocity space. This is the standard formal control condition when pressure is determined by incompressible flow and normal momentum. If the independently imposed direct pressure flux is retained as an additional constraint, the admissible control variations must also satisfy its compatibility conditions along the trajectory; they need not fill the usual space $H$. The standard projection formula alone does not prove stationarity for that more restricted problem. ThenFor the natural solenoidal adjoint space, $P_H{\mathbf{u}}^{†}(0)={\mathbf{u}}^{†}(0)$. Equivalently, its component tangent to the sphere in a normed vector space of fixed kinetic energy vanishes:Here $E_0>0$ is needed for the sphere in a normed vector space to be regular and for division by $2E_0$. The multiplier is real and can have either sign; normalizing the initial adjoint with a prescribed positive sign is not a general necessary condition for minimization. If $E_0=0$, the only feasible initial velocity is zero and this tangent formula for the sphere in a normed vector space is inapplicable. First-order stationarity alone also allows maxima or saddles; a local minimum requires the appropriate nonnegative constrained second variation.

The spatial boundary terms vanish with periodic adjoints in $x,z$, homogeneous no-slip boundary conditions ${\mathbf{u}}^{†}=0$ at $y=\pm 1$, and ${\partial }_y{\theta }^{†}=0$ there. To see this, velocity variations vanish at the wall while their normal derivatives need not, so the viscous boundary term forces the adjoint velocity to vanish. Scalar variations have zero normal derivative but free values, forcing the adjoint scalar's normal derivative to vanish. Pressure variation gives adjoint incompressible flow; no independent terminal or initial datum is assigned to $p^{†}$. Its additive time-dependent constant may be fixed by a zero spatial mean.

In particular, a homogeneous Neumann condition for the direct pressure does not by variational transposition require ${\partial }_y{p}^{†}=0$. For a smooth adjoint solution, the wall-normal adjoint momentum equation instead suppliesat the flat walls, because the direct scalar normal derivative is zero there. A zero adjoint-pressure derivative would be an additional compatibility choice, valid only when this right-hand side vanishes.

There are also literal direct-data compatibility qualifications in the printed setup. The initial total scalar field is $\overline{\theta }=-\mathrm{erf}(30y)$, whose wall derivative is $-60e^{-900}/\sqrt{\pi }\ne 0$. It is extraordinarily small but mathematically not zero. Thus exact initial scalar data and exact zero wall flux are incompatible for a classical solution smooth at $t=0$. One may use a parabolic mild solution of an abstract Cauchy problem with the boundary condition enforced for $t>0$, or replace the initial profile by an exactly compatible smooth one; these are conventions, not an equality $e^{-900}=0$. The reference scalar field is also not a stationary profile of the diffusion equation at finite $Pe$: the full $\Delta \Theta$ term must be retained.

The pressure compatibility at a no-slip wall is similarly important. Direct normal momentum gives ${\partial }_{y}p=\nu \Delta u_y-R{i}_B\theta$. The separately imposed zero pressure derivative requires this right side to vanish. It is not automatic for arbitrary fixed-energy initial velocities. For example, the divergence-free no-slip field generated by the periodic stream function $S=A(1-y^2{)}^2\sin (\alpha x)$ with $\alpha =\pi /L_x>0$ has $u_x={\partial }_{y}S$, $u_y=-{\partial }_{x}S$ and $\Delta u_y{|}_{y=\pm 1}=-8A\alpha \cos (\alpha x)\ne 0$, while $\theta (0)=0$. With $A\ne 0$ adjusted to any positive energy, no classical solution smooth up to the initial wall can satisfy that extra pressure condition. The displayed adjoints are the formal necessary equations along admissible smooth trajectories; a well-posed physical formulation normally determines pressure from incompressible flow and normal momentum, rather than imposing independent homogeneous pressure flux for every control.

The axial angular momentum per unit volume is $mn(r)r{u}_{\theta }(r)=\hslash n(r)$ for a singly quantized positive quantum vortex. Thus every particle in the prescribed shell contributes $\hslash$:This yieldsThis is $\hslash$ times the shell particle number. For a wall-truncated cloud with $R=R_0<R_r$, the corresponding result is $L_z=\pi H_0\hslash n_0[R^2-a_0^2-(R^4-a_0^4)/(2R_r^2)]$. Reversing the quantized circulation reverses $L_z$ without changing the flow kinetic energy. For winding number $q$, $L_z$ is multiplied by $q$.

For repulsive interactions $g>0$, the Thomas–Fermi approximation for a condensate neglects the density-gradient kinetic energy away from boundaries and vortex cores. The stationary Gross–Pitaevskii equation then gives $\mu =V_{\mathrm{ext}}(r)+gn(r)$ wherever the condensate is present. At the axis $\mu =g{n}_0$, henceHere $(q{)}_{+}=\max (q,0)$ and the support is also restricted to the physical bucket $r\le R_0$. The next parts' explicit outer cutoff $R_r$ assumes $R_r\le R_0$; if $R_0<R_r$ the wall truncates the parabola and the corresponding integrals must use $R_0$, not the natural trap radius. A hard-wall boundary layer or the smooth trap-edge healing region is beyond the Thomas–Fermi approximation for a condensate. For an untruncated radial profile, the total particle number is $N=\pi H_0{n}_0{R}_r^2/2$.

The four-momentum isThusMultiplication by $c$ identifies the first two terms as the rest energy $m{c}^2$ and the Newtonian kinetic energy $m|\mathbf{u}{|}^2/2$.

For the fixed-target reaction, the initial invariant isAt threshold the four final particles are at rest in their centre-of-momentum frame, giving invariant $(4mc{)}^2$. Hence $\gamma =7$ and the least laboratory speed is

The conservation of energy to use here is the total kinetic energy plus internal energy inside the shock, not just the internal energy of a uniform-pressure core. With no radiative loss and negligible initial ambient energy,Substitute the adiabatic superbubble similarity solution and define the finite positive dimensionless integralThen $\dot{E}t=4\pi {\rho }_1{\dot{R}}^2{R}^{3}I$. Since $\dot{R}=3R/(5t)$ and $R^5={\alpha }^5(\dot{E}/{\rho }_1)t^3$,This is the energy normalization of a continuously driven spherical shock. It fixes $\alpha$ once the similarity profiles are determined. If the profiles describe only the shocked ambient layer, the energy of the interior hot bubble must also be included before using this normalization. A radiatively cooled shell has a different energy budget; its familiar numerical coefficient should not be substituted into this loss-free model.

For an ideal gas with constant specific-heat ratio $\gamma$, the specific internal energy is $u=c_{V}T$ and $P=\rho (c_P-c_V)T=(\gamma -1)\rho u$. HenceThe stellar virial theorem consequently reads $\ddot{I}/2=2\mathcal{T}+3(\gamma -1)U-3P_{S}V+\Omega$. For a static nonrotating star with negligible surface pressure, $\mathcal{T}=\ddot{I}=P_S=0$, so $\Omega =-3(\gamma -1)U$.

Fully ionized nonrelativistic ions and electrons form a monatomic ideal gas, with $\gamma =5/3$, provided electron degeneracy pressure and radiation pressure are negligible. ThusBeing steady alone does not eliminate rotation or other macroscopic kinetic energy; those are additional assumptions needed for this relation.

Work in an inertial frame centred on the centre of mass, and follow a fixed material mass, with no stellar wind or accretion through its boundary. Assume nonrelativistic motion, Newtonian self-gravity, no external gravitational field, and sufficiently regular integrable fields. DefineHere $I$ is the scalar second mass moment, rather than a moment of inertia about one axis. The kinetic energy $\mathcal{T}$ measures resolved internal motions relative to the centre of mass; random microscopic motion contributes to pressure and internal energy, and must not be counted again in $\mathcal{T}$.

The paper's symmetric compressive tensor $\mathbb{P}$ is minus the tensile stress tensor, so its force per unit volume is $-{\partial }_j{P}_{ij}$. Differentiating the material integral twice and substituting the equation of motion givesIntegration by parts and the divergence theorem turn the stress contribution intoFor the requested scalar stellar virial theorem, assume isotropic pressure, $P_{ij}=P{\delta }_{ij}$, with uniform surface pressure $P_S$. Then $\mathrm{tr}\mathbb{P}=3P$ and the surface integral is $P_S\oint \mathbf{r}\cdot \mathbf{n}dS=3P_{S}V$. No spherical shape is needed for this step.

Finally, symmetrizing the pair integral for the Newtonian gravitational potential givesThereforeAnisotropic magnetic or viscous stresses retain their full trace and boundary terms; a changing integration mass adds transport terms.

The specific orbital energy released in a black-hole accretion flow reaching an inner radius $r_{\mathrm{in}}$ is of order $GM/r_{\mathrm{in}}$. For a circular Keplerian orbit, the binding energy is more precisely $GM/(2r_{\mathrm{in}})$, with the other half of the potential-energy decrease retained as orbital kinetic energy. Multiplying by the mass accretion rate gives the accretion luminosity scaleFor an inner radius of order the Schwarzschild radius, $r_{\mathrm{in}}\sim 2GM/c^2$, this becomesThis is an energy-release scaling, not unit radiative efficiency of black-hole accretion. The actual inner orbit, stresses and advected energy determine the coefficient and the emitted luminosity; the innermost stable circular orbit of a Schwarzschild black hole lies outside its horizon.

Let $\mathbf{\zeta }=\nabla \times \mathbf{u}$ and $\mathbf{\omega }=\mathbf{\zeta }+2\Omega {\mathbf{e}}_z$ be the relative vorticity and absolute vorticity. The nonlinear acceleration identity rewrites the momentum equation asTaking the curl removes the pressure, kinetic energy gradient and shearing-sheet tidal potential. ThusBoth $\nabla \cdot \mathbf{u}$ and $\nabla \cdot \mathbf{\omega }$ vanish. The vorticity equation is consequentlyFor the absolute-vorticity flux tensor $T_{ij}=u_i{\omega }_j-{\omega }_i{u}_j$, use the flux convention $(\nabla \cdot \mathbf{T}{)}_j={\partial }_i{T}_{ij}$, with the first index denoting transport direction. This givesThe scalar $\theta$ is a buoyancy displacement variable, so dimensional consistency gives it units of length; the printed potential-temperature terminology presupposes this normalization. Its curl source directly forces only the $y,z$ components of absolute vorticity. Hence radial vorticity conservation in a shearing sheet takes the local flux formIts volume integral is constant when the net boundary flux vanishes. However $D_t{\omega }_x=\mathbf{\omega }\cdot \nabla u_x$ still permits vortex stretching; the radial component is not generally a materially conserved scalar in three dimensions. Using divergence on the second index of the stated tensor would reverse the transport sign, so the tensor convention matters.

Write $H_{\pm }=H_0\pm d$, with the upper sign referring to $y>0$, and use an f-plane with Coriolis parameter $f_0>0$ for the stated directions. The formulas retain the signed $f_0$; a positive Rossby deformation radius means $R_{\pm }=\sqrt{g{H}_{\pm }}/|f_0|$. Let $\rho$ be the constant fluid density.

Far from the depth step the problem is independent of $y$. The linearized shallow water equations imply conservation of $v_x-f_0\eta /H_{\pm }$. Initially this is $-f_0{\eta }_0\mathrm{sgn}(x)/H_{\pm }$. In geostrophic balance, $u_s=0$ and $v_s=g{\eta }_{s,x}/f_0$, soBoundedness and matching of height and its derivative at $x=0$ give geostrophic adjustment of a surface-height jump:The potential energy release per unit length in $y$ is the finite difference between two individually infinite energies:For comparison, the final kinetic energy is $\rho g{\eta }_0^2{R}_{\pm }/2$, so the energy radiated in inertia-gravity waves is $\rho g{\eta }_0^2{R}_{\pm }$. Integrating the depth-weighted velocity over $x$ gives the volume transportsFor $f_0,{\eta }_0>0$, $Q_{-}$ is towards the step and $Q_{+}$ is away from it. Their difference is $4g{\eta }_{0}d/f_0$; a global steady state cannot simply join the two far-field currents without additional transport along the step.

To obtain the rotating shallow-water height equation with variable depth, set $\zeta =v_x-u_y$ and $J_H={\partial }_x(Hv)-{\partial }_y(Hu)$. Direct differentiation of the original momentum and continuity equations givesConsequentlyAt the step, $H_y=2d\delta (y)$ is a Dirac delta distribution; the original matching conditions are used rather than assuming a smooth depth there. Velocity elimination raises the time order: initially ${\eta }_t=0$ and ${\eta }_{tt}=g{\nabla }_h\cdot (H{\nabla }_h\eta )$, as well as the prescribed initial height.

For a decaying harmonic disturbance, the height equation on each constant-depth side givesSolving the two momentum equations yieldsContinuity of $H\hat{v}$ at the step gives the implicit step-trapped topographic Rossby wave dispersion relationThis derivation assumes ${\omega }^2\ne f_0^2$ and positive decay rates; it describes the low-frequency trapped branch, rather than the radiating inertia-gravity waves.

In the quasi-geostrophic approximation, replace $f_0^2-{\omega }^2$ by $f_0^2$. The resulting explicit frequency, and then its leading small-$d/H_0$ form, areFor $d>0$, the phase velocity has the sign of $f_0$. Thus propagation is eastward in the Northern Hemisphere, with the shallower side on the right. Differentiating the leading dispersion relation givesIn the long-wave limit the group velocity and phase velocity coincide:

To find the slow amplitude, integrate the low-frequency height equation across the step. This gives ${\partial }_t[H{\eta }_y]=f_0[H]{\eta }_x$, where brackets mean the upper-side minus lower-side value. For the given exponentially trapped outer profile, $[H{\eta }_y]=2H_{0}A/R_D$ and ${\eta }_x(x,0)=2{\eta }_0\delta (x)-A_x$. Hence long-wave transport along a depth step obeysThe method of characteristics, interpreted for the Dirac delta distribution source, givesFor $c_T>0$ this is $2{\eta }_0$ in $0<x<c_{T}t$ and zero outside: a surface depression is left in the wake of the wave. The sharp $x$-jump is a long-wave outer approximation. The fast adjustment region of width $O(R_D)$ near the origin, and the detailed wavefront, are unresolved. Keeping a smoothed background ${\eta }_s(x)$ instead would give $A_t+c_T{A}_x=c_T{\eta }_s^{\prime }(x)$ and $A={\eta }_s(x)-{\eta }_s(x-c_{T}t)$ at this order.

The along-step geostrophic flow is $u=-g{\eta }_y/f_0$. Integrating separately over the two sides givesIn the wake, $T_{-}=Q_{-}$ and $T_{+}=-Q_{+}$. The shallow-side transport is eastward and the deep-side transport westward, with magnitudes equal to the respective far-field currents. The net along-step transport is $T_{-}+T_{+}=-4g{\eta }_{0}d/f_0$, balancing their mismatch. These comparisons use the stated leading long-wave and small-depth-contrast approximations.

A principal body frame is an orthonormal basis ${\mathbf{e}}_a(t)$ fixed in the rotating body and aligned with its principal axes. If $R(t)$ is the rotation matrix whose columns are these vectors, differentiating $R^{T}R=I$ shows that $\dot{R}{R}^T$ is skew-symmetric. Every skew-symmetric linear map in three dimensions is cross product with a unique vector $\mathbf{\omega }$, so

For any vector $\mathbf{A}=A_a{\mathbf{e}}_a$,Applying this to the conserved angular momentumgives the Euler equations for a torque-free rigid body:

Multiplying these equations respectively by ${\omega }_i$ and summing proves conservation of kinetic energy; multiplying by $I_i{\omega }_i$ and summing proves conservation of squared angular momentum. Thus the torque-free rigid-body invariants areSolving these two linear equations for the first two squared components givesUsing $I_3{\dot{\omega }}_3=(I_1-I_2){\omega }_1{\omega }_2$ yields the quartic equation for one component of torque-free angular velocity:

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.