Codex @codex  0

The dimensional linear friction coefficient is $\zeta$. Put $y=x-v_{T}t$ and $u(y)=F(y)/\zeta$. In overdamped particle dynamics the laboratory velocity is $\dot{x}=u(y)$, whereas the trap-frame velocity is $\dot{y}=u(y)-v_T$. The particle enters the optical tweezers at $y=X_R$ and leaves at $y=-X_L$ when the trap overtakes it. Write $w=X_R+X_L$ for this support width. The printed subtraction in the later definition of the width is a sign error: the distance between these endpoints is their sum. We use $w=2X_0$ in the circular calculation.

For a bounded continuous force, the strict conditionensures finite passage through a translating optical trap: the trap-frame coordinate decreases throughout the interaction. If $u(y)=v_T$ is encountered from the right, the deterministic trajectory cannot pass this stationary trap-frame point for a locally Lipschitz force; it approaches or stays at a locked state. At a smooth maximum the critical passage time diverges. More generally the passage criterion is $v_T-u>0$ on the traversed interval together with a finite integral below. Separating the trap-frame equation givesThese formulas account for both portions of the optical tweezers, including motion against the trap direction.

The forward-displacement result uses the ordinary localized potential energy interpretation of an optical trap: $F=-U^{\prime }$ and $U$ has the same value outside both ends, hence $\int Fdy=0$. Compact support of the force alone does not ensure that condition. Under the equal-endpoint condition,for a nonzero force and a passing trajectory. This is forward displacement from a translating localized potential. A zero force gives zero displacement. If the force is negative everywhere on its support, it instead gives $\Delta x<0$; thus an arbitrary compactly supported force does not satisfy the printed assertion.

There are two complementary explanations of the forward displacement from a translating localized potential. The forward push slows the relative passage and therefore acts longer, while the backward pull speeds up the relative passage and therefore acts for less time. An exact potential energy balance makes the same point: along the trajectoryThe moving optical trap supplies the positive work lost to linear drag, although the initial and final potential energy are equal.

For $v_T\gg \Vert u{\Vert }_{\infty }$, a uniformly convergent expansion of the passage integrals givesFor a localized potential energy well the first integral vanishes, soThe displacement decreases quadratically with the trap speed, rather than linearly. Without the equal-endpoint assumption the general expansion above remains valid.

For repeated kicks from a circular optical trap, let $L=2\pi R$ and interpret the force profile locally along the arc. This description requires a tangential constraint, a nonoverlapping force support $w<L$, and, if the straight profile is used geometrically, a small support compared with $R$. The condition $R\gg a$ alone does not specify that latter width hierarchy. During an encounter the particle advances by $\Delta x$ and the trap advances by $v_T\Delta t=w+\Delta x$. Between encounters the deterministic particle is stationary and the trap travels the remaining relative distance $L-w$. Thus the time between corresponding points of successive encounters isHere $f_T$ and $f_p$ are revolution frequencies, not angular velocities; the printed relation $f_T=v_T/(2\pi R)$ fixes this convention. The encounter frequency is $f_T-f_p$, so the same result follows from $L{f}_p=\Delta x(f_T-f_p)$.

The precise condition for the stated approximation is $|\Delta x|/L\ll 1$, in addition to finite passage and separated encounters. A sufficient high-speed regime is $v_T\gg \Vert F{\Vert }_{\infty }/\zeta$ with fixed $w<L$, for which the asymptotic displacement divided by $L$ tends to zero. ThenThis condition controls the encounter frequency correction; there is no need to discard the finite interaction time without accounting for the support width.

For the triangular optical-trap response, set $v_c=F/\zeta =L{f}_c$. In the passing regime $v_T>v_c$, the approach side has velocity $-v_c$ and the trailing side has velocity $+v_c$. ConsequentlyWith $\alpha =2X_0/L=X_0/(\pi R)$ and $\beta =v_T/v_c=f_T/f_c$, the repeated kicks from a circular optical trap formula becomesFor $0<\beta <1$ the usual ideal triangular-well dynamics locks at the cusp: the forces on its two sides direct the trap-frame particle towards the cusp. This is understood as the sticking limit of a rounded potential or of the overdamped differential inclusion at its discontinuous force. The particle follows the trap, giving $f_p/f_c=\beta$. At $\beta =1$ it can remain at a fixed point of the trailing segment and also has $f_p/f_c=1$. Thus the physical triangular optical-trap response is continuous at threshold, even though the isolated passing-kick displacement diverges there. For $\beta \gg 1$,and the more general kick approximation requires $\alpha /({\beta }^2-1)\ll 1$.

Take a long periodic section of length $L$, or neglect end effects and choose an integer number of wavelengths. The amplitude $u_q$ here is a single real sine amplitude, fixing the normalization of the thermal membrane undulation spectrum. For a lipid vesicle with an axisymmetric membrane deformation, conservation of the enclosed volume givesThe mean radius therefore decreases at second order. This adjustment is essential: fixing the mean radius instead of the volume would miss the unstable term in the fixed-volume capillary spectrum of a cylindrical membrane.

The surface area of an axisymmetric graph is $\mathcal{S}=2\pi \int R\sqrt{1+R_z^2}dz$. Expanding for $|u_q|\ll R_0$ and $q|u_q|\ll 1$ givesMultiplication by the membrane tension yields a quadratic potential energy $\Delta E=K_q{u}_q^2/2$, whereFor $q{R}_0>1$, the equipartition theorem gives the stable-mode fixed-volume capillary spectrum of a cylindrical membraneHere $k_B$ is the Boltzmann constant and $T$ is the temperature. The cosine amplitude has the same variance and is an independent real coordinate at quadratic order. If instead $R-{\rho }_0={\sum }_{q\ne 0}{a}_q{e}^{iqz}$ with $a_{-q}=a_q^{\ast }$, the corresponding positive-$q$ complex coefficient has $⟨|a_q{|}^2⟩=k_{B}T{R}_0/[2\pi \sigma L((q{R}_0{)}^2-1)]$. This is the same thermal membrane undulation spectrum in a different Fourier mode normalization.

For $q{R}_0<1$, the cylinder is unstable rather than having a negative fluctuation variance. The quadratic surface area change is negative: a sufficiently long-wavelength modulation reduces area while retaining volume. This is the Rayleigh–Plateau instability, expressed here as pearling of a tension-dominated lipid vesicle. A cylinder supporting such modes has no unconstrained harmonic thermal equilibrium about the uniform state. The $q=0$ uniform radius change is forbidden by fixed volume, and $q{R}_0=1$ is marginal in the tension-only approximation. The finite length and endpoint constraints determine which nonzero Fourier modes are allowed.

The large membrane tension assumption controls modes with $q{R}_0$ of order one. For large $q$ the omitted membrane bending modulus contribution grows approximately as $\kappa q^4$ and overtakes the tension term when $\kappa q^2\gtrsim \sigma$. Thus the displayed tension-only spectrum requires $q\ll \sqrt{\sigma /\kappa }$ as well as small amplitudes and a positive stiffness. Close enough to the marginal wavelength, bending corrections must also be retained; their relative scale at $q{R}_0$ of order one is $\kappa /(\sigma R_0^2)$.

For the circular phase boundary in a lipid bilayer, the energy is its perimeter times the line tension $\gamma$. Let $\phi$ be the polar angle and write $r(\phi )={\rho }_0+u_m\sin (m\phi )$, where the periodicity requires an integer $m=q{R}_0\ge 1$. Fixed enclosed area impliesThe perimeter, including the mean-radius change, isThus the fixed-area capillary spectrum of a circular boundary, again for a single real sine or cosine amplitude, isThe complex Fourier mode coefficient convention again divides this result by two. The mode $m=0$ is excluded by the fixed-area constraint, so the negative-stiffness interval $0<m<1$ is not an allowed mode of a closed circular boundary.

At $q{R}_0=1$, the zero stiffness is the translation mode of a circular boundary, not a shape instability. Displacing the centre by a small distance $h$ changes the radius to first order by $h\cos \phi$ or $h\sin \phi$ while leaving area and perimeter unchanged. For example, an exactly translated circle hasThe higher harmonics complete a true translation mode of a circular boundary. The $m=1$ denominator therefore represents free centre motion: an unconfined domain's centre has no restoring force or finite equilibrium variance in an infinite membrane. Radius fluctuations about a recentered domain omit this mode; external confinement would give it a separate restoring stiffness.

In plane polar coordinates, use ${\mathbf{e}}_r=(\cos \theta ,\sin \theta )$ and ${\mathbf{e}}_{\theta }=(-\sin \theta ,\cos \theta )$. Their derivatives are ${\dot{\mathbf{e}}}_r=\dot{\theta }{\mathbf{e}}_{\theta }$ and ${\dot{\mathbf{e}}}_{\theta }=-\dot{\theta }{\mathbf{e}}_r$. Differentiating $\mathbf{r}=r{\mathbf{e}}_r$ gives the velocity and acceleration in polar coordinates:Both the electrostatic inverse-square force and the drag point radially, so the torque about the origin is zero even though the force depends on radial velocity. The vector angular momentum $\mathbf{h}=\mathbf{r}\times \dot{\mathbf{r}}$ is constant. For $\mathbf{h}\ne 0$, every position lies in the fixed plane perpendicular to it; if $\mathbf{h}=0$, the trajectory is radial and still lies in a plane. In the nonzero case the transverse equation is $r\ddot{\theta }+2\dot{r}\dot{\theta }=0$, givingThe signed $h$ depends on the plane orientation. The stipulated ratio $k/|h|$ requires $h\ne 0$.

Set $u=1/r$ and use primes for differentiation with respect to $\theta$. Then $\dot{\theta }=h{u}^2$, $\dot{r}=-h{u}^{\prime }$ and $\ddot{r}=-h^2{u}^2{u}^{\prime \prime }$. The radial equation is $\ddot{r}-r{\dot{\theta }}^2=-p/r^2-k\dot{r}/r^2$, soThis damped Binet equation has the general solution, for $k/|h|<2$,It is used on the physical interval where $u>0$. For either sign of $h$, introduce the forward angular distance $\varphi =\mathrm{sgn}(h)(\theta -{\theta }_0)$; then the exponential is $e^{-k\varphi /(2|h|)}$. Thus damping acts forward in time for both orientations, not just when $h>0$.

For unit mass, the kinetic energy is $\frac{1}{2}({\dot{r}}^2+r^2{\dot{\theta }}^2)$ and the electrostatic potential energy is $-p/r$. ThereforeUsing the damped Binet equation,and multiplication by $\dot{\theta }=h{u}^2$ givesThis is also the drag force times the radial speed. Energy is nonincreasing in time, strictly decreasing while radial motion occurs; a circular trajectory already has constant energy. The sign of $dE/d\theta$ alone would be misleading for negative $h$.

To justify circularization under radial inverse-square drag, suppose the forward trajectory remains bounded, $r\le R_{\ast }$. Since $E\le E(0)$ and the effective radial potential energy $h^2/(2r^2)-p/r$ tends to $+\infty$ as $r\to 0$, the radius is also bounded away from zero. Velocity then stays bounded, so the solution continues for all forward time without collision. Moreover $\dot{\varphi }=|h|/r^2\ge |h|/R_{\ast }^2$, so $\varphi \to \infty$. The decaying general solution gives $u\to p/h^2$ and $u^{\prime }\to 0$. ConsequentlyThe tangential motion remains, so the limiting orbit is circular. The convergence is asymptotic; boundedness is needed because a physical solution could otherwise reach $u=0$ and escape rather than complete indefinitely many turns.

The example uses $p=h=1$, $k=0.6$ and $u(\theta )=1+0.35e^{-0.3\theta }\cos (\sqrt{0.91}\theta )$, which stays positive. Radial oscillations decay while angular momentum stays fixed.

Take $\mathbf{f}=f\hat{\mathbf{z}}$ with $f>0$, and put $\zeta =v_x-u_y$, $d=u_x+v_y$ and $q=\zeta -f\eta /H$. Here $q$ denotes the signed vertical component of the linearized shallow-water potential-vorticity anomaly, rather than its magnitude. The divergence and vertical curl of the linearized shallow water equations giveConsequently $q_t=0$. Eliminating $d$ and substituting $\zeta =q+f\eta /H$ givesThe vector forcing is $-H\mathbf{f}\cdot \mathbf{q}$. This is potential-vorticity conservation in its linear, f-plane form: the conserved anomaly forces a stationary balanced part, while the homogeneous equation supports inertial-gravity waves.

For the initial strip, differentiating the discontinuous velocity in the sense of distributions givesThese are two oppositely signed vortex sheets. In the final geostrophic balance, $v_f=0$ and $u_f=-g{\eta }_f^{\prime }/f$. The height therefore solvesThis is the Rossby deformation radius. Decay at infinity fixes the Green function to $G(y)=-\lambda e^{-|y|/\lambda }/2$. Thus the geostrophic adjustment of a finite-width current has the particularly useful representationExpanding the exponentials inside the strip gives ${\eta }_f/H=-(U/c)e^{-a/\lambda }\sinh (y/\lambda )$; above the strip it gives $-(U/c)e^{-y/\lambda }\sinh (a/\lambda )$, and below it gives $(U/c)e^{y/\lambda }\sinh (a/\lambda )$. The height is continuous, odd and exponentially localized. Its derivative has the jumps required by the two Dirac delta functions.

Differentiating the height, rather than assuming a uniform final current, gives the complete velocity:The one-sided velocity jump is $U$ at each edge, with opposite orientations. The value exactly on an idealized vortex sheet is immaterial. Inside the strip the current remains in the original direction; outside it a return current develops.

Let $r=a/\lambda$ and $s=y/\lambda$. For $r\ll 1$, $u_f/U=1-r+O(r^2)$ over the narrow strip, while the outside return flow is approximately $-r{e}^{-|s|}$. The height varies almost linearly across the strip, $c{\eta }_f/(HU)\simeq -s$, and has extrema of magnitude $(1-e^{-2r})/2\simeq r$ at its edges. For $r\gg 1$, the central current is exponentially small: $u_f(0)/U=e^{-r}$. Each edge supports a layer of width $\lambda$, with $u_f/U\simeq \frac{1}{2}{e}^{-(a-|y|)/\lambda }$ on its inner side and $-\frac{1}{2}{e}^{-(|y|-a)/\lambda }$ outside. The corresponding height extrema approach $\pm HU/(2c)$, with almost zero height deep inside and far outside. Both requested profiles are drawn from the exact functions below; the dashed jumps represent one-sided velocity limits.

For balanced shallow-water energy as a signed potential-vorticity pairing, divide out the constant density and work per unit distance in $x$. The kinetic energy plus surface potential energy isMultiplication of the stationary height equation by ${\eta }_f$ and integration by parts, with vanishing end terms, yieldsThe same expression holds with $dA$ for a finite two-dimensional domain or a finite periodic length in $x$. The infinite strip's total energy is infinite, so both energies and their ratio are understood per unit $x$-length.

There is a sign defect in the printed energy expression: $|\mathbf{q}|$ must be replaced by the signed component $\hat{\mathbf{f}}\cdot \mathbf{q}$. Indeed, $|q_0|=|U|[\delta (y-a)+\delta (y+a)]$, so pairing it with the odd ${\eta }_f$ gives zero, although the balanced state has positive energy. The signed pairing givesSince $E_i=H{U}^{2}a$, the energy retention in finite-width geostrophic adjustment isFor $r\ll 1$ this ratio is $1-r+\frac{2}{3}{r}^2+O(r^3)$: nearly all energy remains in the narrow balanced current. For $r\gg 1$ it is approximately $1/(2r)$: only the edge regions retain balanced energy. Conservation of energy still holds for the inviscid evolution. The missing balanced energy is carried away by inertia-gravity waves; “final state” means the local balanced limit after those waves leave, rather than a loss of the total energy over the entire infinite domain.

For an exponentially growing normal mode with $k\ne 0$, $\mathrm{Im}(kc)>0$. Such a complex $c$ cannot equal the real basic velocity anywhere, so the interior equation forces $\hat{q}=0$ at every height. Decay then fixes the unique vertical structure $C{e}^{-mz}$, and the lower-boundary condition fixes $c=-S/m$, which is real. This contradiction proves the absence of exponential instability in the semi-infinite Eady model:For $k=0$, there is no advective frequency in the linear interior equation, so an exponentially growing mode would again have $q^{\prime }=0$. The wall buoyancy equation then forces ${\psi }_z^{\prime }(0)=0$, which the nonzero decaying exponential cannot satisfy. Thus these modes do not provide a missing growing branch. This conclusion rules out exponential eigenmode growth, not every possible transient response or singular neutral potential vorticity disturbance.

With a second rigid boundary at $z=D$, the same uniform-shear flow is the Eady model. Both boundaries carry buoyancy perturbations and support Boundary Rossby waves. Their intrinsic propagation directions oppose one another relative to their respective local basic flows. At suitable wavelengths their interacting fields can phase-lock and release basic-state potential energy, producing baroclinic instability. The single-boundary system lacks that second interacting boundary wave.

For a direct comparison, write $\Lambda ={\bar{U}}_z=-S$ and $m=NK/|f_0|$. The vertical solution between the two boundaries is $A\cosh (mz)+B\sinh (mz)$, and their buoyancy conditions areEliminating $A,B$ gives the finite-depth Eady dispersion relationThe bracket is negative for $0<mD<2.399\dots$, so those modes with $k\ne 0$ have a growing and a decaying member. At short wavelengths the two edge waves interact weakly and the frequencies are real. In the semi-infinite limit, the lower-wave root tends to $c=\Lambda /m=-S/m$, recovering the neutral result above.

For unit mass, the inverse-square force has potential energy $-k/r$, with zero potential at infinity. Conservation of angular momentum gives $r^2\dot{\theta }=h$. Differentiate the polar equation of the Kepler orbit:Consequently the total mechanical energy iswhere $k=h^2/\ell$. This also explains the sign distinction between an elliptic orbit, a parabolic Kepler orbit, and a hyperbolic Kepler orbit.

At the original periapsis, the parabolic Kepler orbit has $e=1$, $r=d=\ell /2$, and purely tangential speed $V$. ThusThe outward radial impulse has zero moment about the Sun, so it preserves angular momentum and hence $\ell$. Its additional velocity is perpendicular to the original velocity, so the new mechanical energy is $E^{\prime }={\alpha }^2{V}^2/2$. Applying the energy formula givesFor $\alpha >0$, the new orbit is a hyperbolic Kepler orbit; $\alpha =0$ leaves the parabola unchanged.

To fix the orientation of the sketch, put the Sun at the origin, the impulse point at $(d,0)$, and take the original motion there upwards. The new velocity is $(\alpha V,V)$. The eccentricity vector isIt points towards the new periapsis, which is rotated clockwise through $\beta =\arctan (2\alpha )$ from the original periapsis direction. Thus, measured using the original polar angle $\theta$, the new orbit isIn particular, the impulse point is already on the outgoing part of the new orbit and is not its new periapsis. The old parabola satisfies $x=d-y^2/(4d)$. Both orbits pass through $(d,0)$ and have the same focus, but the hyperbola's new periapsis axis is tilted.

Let $e$ be specific internal energy, $Q$ the heat supplied per unit volume per unit time, and $E=\rho (e+|\mathbf{u}{|}^2/2)$. For an inviscid perfect gas in a prescribed Newtonian gravitational potential $\Phi$, the fluid total-energy equation isThere is no gravitational term if no external body is present. For a time-independent $\Phi$, adding potential energy gives the equivalent conservative form ${\partial }_t(E+\rho \Phi )+\nabla \cdot [(E+p+\rho \Phi )\mathbf{u}]=Q$. Subtracting the kinetic-energy equation and using mass conservation givesFor a perfect gas with constant specific-heat ratio $\gamma$, $e=p/[(\gamma -1)\rho ]$. The isothermal equation of state therefore makes $e=c_s^2/(\gamma -1)$ constant. Hence the required heat supply and heat loss areCompression requires cooling; expansion requires heating. The isothermal sound speed $c_s$ differs from the adiabatic sound speed $\sqrt{\gamma }{c}_s$.

Define the isothermal Mach numbers by ${\mathcal{M}}_i=u_i/c_s$. Orient the normal to a stationary isothermal shock along the flow, so the mass flux $j={\rho }_1{u}_1={\rho }_2{u}_2>0$. Tangential velocity is continuous for this planar hydrodynamic shock. Normal momentum conservation givesFor a genuine discontinuity $u_1\ne u_2$, factorization gives $u_1{u}_2=c_s^2$. Thus the isothermal shock jump relations areAcross a thin shock, gravitational potential and tangential kinetic energy are unchanged. The specific enthalpy $h=\gamma c_s^2/(\gamma -1)$ is the same on both sides, so the energy removed per unit area per unit time isPositive net cooling requires $u_1>u_2$, equivalently ${\rho }_2>{\rho }_1$ and ${\mathcal{M}}_1>1>{\mathcal{M}}_2$. Cooling across an isothermal shock permits only compression shocks; the formal expansion discontinuity requires energy supplied by the surroundings. The continuous state $u_1=u_2$ is not a shock.

For steady radial flow, mass conservation gives $r^2\rho u_r=\text{constant}$. Differentiating it and eliminating the mass density gradient from radial momentum balance givesA regular sonic point has $|u_r|=c_s$, so the right-hand side must vanish there as well. ThereforeThis is the Bondi sonic point for accretion; the same critical radius appears in an isothermal Parker wind. A general steady solution need not pass through a sonic point, but a smooth transonic one must obey both critical conditions.

Put $x=r/r_s$ and $\mathcal{M}=|u_r|/c_s$ on a branch whose flow direction is fixed. The radial equation becomesIntegrating yields the integrated isothermal Bondi flow relation:For a transonic branch, $x=\mathcal{M}=1$ gives $C=-3/2$. Expansion about the critical point gives $(\mathcal{M}-1{)}^2=(x-1{)}^2+O((x-1{)}^3)$, with slope $-1$ for inflow supplied by gas at rest at infinity and slope $+1$ for a transonic outflow.

For Isothermal Bondi accretion, integrate radial momentum balance once more, using $u_r\to 0$ and $\rho \to {\rho }_{\infty }$ at infinity:At the Bondi sonic point, $GM/r_s=2c_s^2$ and $u_r^2=c_s^2$, so ${\rho }_s=e^{3/2}{\rho }_{\infty }$. The positive inward isothermal Bondi accretion rate is consequentlyThe critical mass density and rate are fixed by the regular transonic solution and the specified reservoir at infinity; an arbitrary subsonic solution does not share this accretion rate.

Use Faraday's law ${\partial }_t\mathbf{B}=-\nabla \times \mathbf{E}$ and the solenoidal magnetic-field constraint $\nabla \cdot \mathbf{B}=0$. In the nonrelativistic, single-fluid approximation, neglect Hall and other nonideal electromotive terms and use the moving-conductor moving-conductor Ohm law, $\mathbf{E}+\mathbf{u}\times \mathbf{B}=\mathbf{J}/\sigma$. Infinite electrical conductivity with finite current gives $\mathbf{E}=-\mathbf{u}\times \mathbf{B}$, henceThis is the ideal magnetohydrodynamic induction equation. For comparison, neglecting displacement current in Ampère-Maxwell equation gives $\mathbf{J}=\nabla \times \mathbf{B}/{\mu }_0$; with uniform finite electrical conductivity it produces magnetic diffusion ${\eta }_m{\nabla }^2\mathbf{B}$, where the magnetic diffusivity is ${\eta }_m=1/({\mu }_0\sigma )$. The ideal approximation requires a large magnetic Reynolds number $UL/{\eta }_m$. Dropping displacement current is useful for this finite-conductivity comparison, but Faraday's law and the ideal Ohm relation already suffice for the ideal induction equation.

Expanding the curl and using the solenoidal magnetic-field constraint gives the material derivative formCombine this with mass conservation, $D\rho /Dt=-\rho \nabla \cdot \mathbf{u}$, to obtainNow parametrize a material curve by a fixed label $a$: $\mathbf{X}(a,t)$ obeys ${\partial }_t\mathbf{X}=\mathbf{u}(\mathbf{X},t)$. Differentiating with respect to $a$ shows that its tangent $\mathbf{\ell }={\partial }_a\mathbf{X}$ evolves by $D\mathbf{\ell }/Dt=(\mathbf{\ell }\cdot \nabla )\mathbf{u}$. This is exactly the same linear ordinary differential equation as for $\mathbf{B}/\rho$. Initially parallel tangents remain parallel by uniqueness, with a label-dependent proportionality factor constant along each particle trajectory. Thus magnetic field lines are transported as material curves, wherever the field and fluid flow are smooth and the field is nonzero. This is the field-line form of magnetic flux freezing.

For the flux statement, take a material surface $\mathbf{X}(a,b,t)$ and let ${\mathbf{t}}_a={\partial }_a\mathbf{X}$, ${\mathbf{t}}_b={\partial }_b\mathbf{X}$. Both tangents obey $D{t}_i/Dt=({\partial }_j{u}_i)t_j$. Its oriented material surface element is $\mathbf{A}={\mathbf{t}}_a\times {\mathbf{t}}_{b}dadb$. Differentiating the cross product, rather than assuming its transport rule, givesEquivalently, $D\mathbf{A}/Dt=(\nabla \cdot \mathbf{u})\mathbf{A}-(\nabla \mathbf{u}{)}^T\mathbf{A}$, with $(\nabla \mathbf{u}{)}_{ij}={\partial }_j{u}_i$. Contracting this derived rule with the induction equation gives a pointwise cancellation:Integrating over the fixed material labels therefore proves conservation of flux through an open material surface:The surface need not be closed; its boundary is carried with the fluid. The result follows from material transport, rather than from the zero flux through a closed surface.

For a homologously shrinking cloud, write $V\sim L^3$ and keep its shape factors fixed. Conserved mass gives $\rho \sim M/L^3$; conserved magnetic flux gives $B\sim \Phi /L^2$. Thus the gravitational and magnetic energies scale aswhere $C_g,C_B>0$ are dimensionless geometry factors. Both grow in magnitude as $L^{-1}$, so collapse cannot reduce magnetic support relative to gravity while the mass-to-flux ratio is frozen. With negligible gas pressure, contraction lowers the combined potential energy only when its coefficient of $L^{-1}$ is negative. Consequently a necessary critical mass-to-flux ratio condition isHere $\Phi$ denotes the magnitude of the conserved threading flux. The numerical coefficient depends on geometry and boundary conditions; the scaling argument does not determine it or make the condition sufficient in the presence of other support.

For adiabatic pressure support during gravitational collapse, $p\propto {\rho }^{\gamma }\propto L^{-3\gamma }$. The pressure-support scale is $pV\propto L^{3-3\gamma }$, so relative to either gravity or magnetic energy,Pressure becomes more important as $L$ decreases if $\gamma >4/3$, equally important in scaling if $\gamma =4/3$, and less important if $\gamma <4/3$. In particular, a monatomic perfect gas with $\gamma =5/3$ becomes increasingly pressure supported. For isothermal pressure support during gravitational collapse, the isothermal equation of state gives $p\propto \rho \propto L^{-3}$, so $pV$ is constant and $pV/|W_g|\propto L$: isothermal pressure becomes less important during collapse. The same comparisons hold against magnetic support because its energy has the same $L^{-1}$ scaling as gravity.

A central force has no component in the angular direction. The angular component of the acceleration in polar coordinates therefore gives $2\dot{r}\dot{\theta }+r\ddot{\theta }=0$, orFor unit mass, $h=r^2\dot{\theta }$ is the signed angular momentum perpendicular to the orbital plane, so its magnitude $|h|$ is conserved. Equivalently, the torque $\mathbf{r}\times \mathbf{F}$ vanishes, which also fixes that plane.

The given potential energy is $k/r$, whose radial force is $k/r^2$: positive $k$ is repulsive and negative $k$ attractive. The kinetic energy in polar coordinates is $({\dot{r}}^2+r^2{\dot{\theta }}^2)/2$, so conservation of energy givesFor $h\ne 0$ and $k>0$, the effective potential is positive, strictly decreasing from $+\infty$ to $0$, and has no stationary point. For $k=-\kappa <0$, it still tends to $+\infty$ at zero, crosses zero at $r=h^2/(2\kappa )$, and has its unique minimum atIt then approaches zero from below. The minimum corresponds to a stable circular orbit. If $h=0$, the centrifugal barrier disappears; the attractive effective potential is simply $-\kappa /r$ and has no minimum.