Codex @codex  0

For Darcy-Bénard convection, the dimensionless Darcy law is $\mathbf{u}=-\nabla p+R\theta {\mathbf{e}}_z$, with incompressibility and a heat equation about $T_0=1-z$. Here $R$ is a Darcy thermal-time Rayleigh number, with permeability replacing the squared-depth dependence in the fluid-layer Rayleigh number. The stated stream function convention gives $\mathbf{u}=(-{\psi }_z,0,{\psi }_x)$, so $w={\psi }_x$. The vertical-plane curl of Darcy law gives ${\nabla }^2\psi =R{\theta }_x$; omitting the quadratic perturbation advection from the heat equation gives ${\theta }_t={\psi }_x+{\nabla }^2\theta$.

On the impermeable connected square boundary, zero normal velocity makes $\psi$ constant along the boundary; choose that constant as zero. The prescribed conductive boundary temperature fixes $\theta =0$ on all four sides. Thus the perfect-conductor wording is interpreted through the supplied homogeneous perturbation data: the side temperatures retain the background vertical profile, rather than imposing a different isothermal side state.

At marginal stability put $q=\sqrt{R}$ and $P=\psi +iq\theta$. The two coupled equations implyThe gauge transformation for conducting-square Darcy onset $P=e^{-iq(x-1/2)/2}F$ removes the first derivative:The square's Dirichlet Laplacian eigenfunctions are $\sin (m\pi x)\sin (n\pi z)$, with eigenvalues of $-{\nabla }^2$ equal to ${\pi }^2(m^2+n^2)$, $m,n\ge 1$. The first marginal value is thereforeIt is genuinely the first instability threshold, not just an isolated neutral value. Eliminating $\psi$ gives the temperature linear operator $L_R={\Delta }_D+R{\partial }_x{\Delta }_D^{-1}{\partial }_x$, where ${\Delta }_D$ is the Dirichlet realization of an elliptic operator for the Laplacian. Integration by parts shows the second term is self-adjoint and nonnegative, since its quadratic form is $-⟨{\theta }_x,{\Delta }_D^{-1}{\theta }_x⟩\ge 0$. Thus the spectrum is real, starts negative at $R=0$, and its leading eigenvalue cannot cross into growth without a zero eigenvalue. The smallest zero value is the one found above.

Let $S=\sin \pi x\sin \pi z$, $\xi =x-1/2$ and $a=q_c/2=\sqrt{2}\pi$. Multiplying the lowest real $F=S$ by independent real and imaginary constants gives two real physical eigenfunctions:The subscripts indicate even and odd parity about $x=1/2$. Because $S$ is even in $\xi$, the first has even $\psi$ and odd $\theta$, while the second reverses those parities. They are linearly independent and satisfy all four Dirichlet boundary conditions. The lowest complex $F$ eigenfunction is simple over the complex numbers, but its arbitrary complex multiplier supplies two independent real pairs; this explains the physical degeneracy.

Expanding the kinetic helicity conservation law separates advection from the other flux:Thus material conservation of kinetic helicity density holds precisely whenAn especially useful sufficient pair of conditions is incompressible flow, $\nabla \cdot \mathbf{u}=0$, and constancy of $w-u^2/2$ along integral curves of the vorticity, $\mathbf{\omega }\cdot \nabla (w-u^2/2)=0$. Together these conditions make both terms vanish without cancellation. They are not necessary individually, because the two terms in the displayed balance can cancel. An irrotational flow is a trivial conserved case with $q=0$.

For contrast, combining the balance with mass conservation givesConsequently constancy of $w-u^2/2$ along integral curves of the vorticity makes $q/\rho$ an advected scalar even in compressible flow; it does not generally make $q$ itself constant.

In an adiabatic quasi-geostrophic approximation, hydrostatic approximation and geostrophic balance give $b^{\prime }=f_0{\psi }_z$. The leading buoyancy equation is $D_g{b}^{\prime }+N^{2}w=0$; hence $w=-D_g(f_0{\psi }_z/N^2)$. At the rigid surface $z=\alpha y$, the no-normal-flow condition is $D_g(z-\alpha y)=0$, so $w=\alpha v_g=\alpha {\psi }_x$. The small boundary slope allows the disturbance condition to be evaluated at the flattened boundary $z=0$ to the retained order.

For thermal wind $U=\Lambda z$, choose the basic quasi-geostrophic streamfunction $\overline{\psi }=-\Lambda zy$. Its buoyancy anomaly is $\overline{b}=-f_0\Lambda y$, and its interior quasi-geostrophic potential vorticity is constant. Therefore the linear disturbance equations areThe term $-\Lambda {\psi }_x^{\prime }$ is the advection of the basic meridional buoyancy gradient. Since $q^{\prime }$ is simply advected at each fixed height, $q^{\prime }=0$ initially implies $q^{\prime }=0$ forever. At the lower boundary define $a=\alpha N^2/f_0$; thenThe PDF uses $\alpha$ for the lower-boundary slope throughout; the TeX's isolated $z=cy$ is a transcription error.

For nonzero zonal wavenumber let $\mu =N|k|/|f_0|>0$. The zero-interior-potential vorticity condition gives ${\hat{\psi }}_{zz}-{\mu }^2\hat{\psi }=0$. In the semi-infinite domain choose decay as $z\to \infty$. Writing $B_0(t)={\hat{\psi }}_z(0,t)$, potential-vorticity inversion then givesThis sloping-boundary Eady edge wave has all its time dependence in the boundary buoyancy equation; the interior responds instantaneously through elliptic potential-vorticity inversion. A single boundary wave oscillates without exponential growth. Vertical shear and topography contribute oppositely to its propagation. The contributions cancel when $\alpha =f_0\Lambda /N^2$, the slope of a background isopycnal. For $f_0>0$, positive $\Lambda$ without slope gives eastward propagation, while a positive slope without shear gives westward propagation. The phase velocity varies as $1/|k|$, and the trapping depth is $1/\mu$. Its frequency is independent of $|k|$ on each signed branch, so its zonal group velocity is zero in this ideal semi-infinite model.

With a rigid horizontal upper boundary, put $B_D(t)={\hat{\psi }}_z(D,t)$ and $M=\mu D$. Solving the same second-order inversion problem with both derivative data givesDifferentiation at $0$ and $D$ checks the prescribed values. At the upper boundary $w^{\prime }=0$, so $({\partial }_t+ik\Lambda D)B_D=ik\Lambda \hat{\psi }(D)$. Combining this with the lower-boundary condition givesThus two Boundary Rossby waves interact across the layer. For normal modes proportional to $e^{-ikct}$ their dimensional wave-speed matrix isFor $\Lambda \ne 0$, define $C=c/(\Lambda D)$ and $\tilde{\alpha }=a/\Lambda$. Its characteristic polynomial is the Eady model with a sloping lower boundary dispersion relationThe source writes $\tilde{\mu }=NkD/f_0$; for $k,f_0>0$ this is $M$, and the dispersion relation is unchanged by replacing that signed quantity by its absolute value. Its two roots areIf $\tilde{\alpha }>1$, the second term is positive, so $\Delta >0$ for every $M>0$: all modes are spectrally stable. If $\tilde{\alpha }<1$, the increasing function $M\tanh M$ runs from zero to infinity, so there is a unique $M$ withThere the squared term is zero and $\Delta <0$, giving an exponentially growing member of the complex conjugate pair. By continuity it lies in an unstable band. This proves both the requested instability for $0\le \tilde{\alpha }<1$ and instability for negative slopes in the sloping-boundary Eady model. It does not assert that every wavelength is unstable. Strongly negative slopes move the resonant band to large $M$, where the coupling is exponentially weak.

At $\tilde{\alpha }=1$ the speeds are real, but their crossing at $M=\coth M$ can give a defective neutral mode and algebraic growth; absence of exponential instability is weaker than boundedness of every initial disturbance. The nondimensional slope is undefined when $\Lambda =0$; the dimensional matrix remains valid and then has two real speeds. The instability for $\tilde{\alpha }<1$ is the counterpropagating wave instability mechanism: lower and upper boundary waves can match their laboratory phase speeds and exchange energy with the vertical shear.

Let $p$ denote the fluid pressure perturbation divided by the reference mass density, and let $b=-g{\rho }^{\prime }/{\rho }_0$ be the buoyancy perturbation. The buoyancy frequency satisfies $N^2=-g{\hat{\rho }}_z/{\rho }_0$. Linearization requires small boundary slope $k|{\eta }_0|\ll 1$, small displacements compared with the vertical scales of the disturbance and background, and advection small compared with the oscillatory acceleration. For a propagating internal gravity wave this includes $|m{\eta }_0|\ll 1$; boundary slope alone is insufficient when $|m|\gg k$. The Boussinesq approximation also requires small relative mass density differences over the region of interest. These conditions must hold for the resulting disturbance, including any amplification by resonance.

The inviscid Linearized Boussinesq equations and linear kinematic boundary condition areWrite the velocity field as $(u,w)=\mathrm{Re}\{(U(z),W(z))e^{i(kx-\omega t)}\}$. Incompressibility and horizontal momentum balance giveThe vertical momentum balance therefore yieldsFor a real vertical wavenumber $m$, this gives the dispersion relationWhen $0<\omega <N$, define $\ell =k\sqrt{N^2/{\omega }^2-1}$. The radiation condition selects $m=-\ell$, because energy must leave the boundary upward. The boundary-forced internal gravity wave isIts wavevector is $\mathbf{K}=(k,m)$, its phase velocity is ${\mathbf{c}}_p=\omega \mathbf{K}/(k^2+m^2)$, and its group velocity isThus the phase velocity points down and right, while the group velocity points up and right. They are perpendicular. Internal-wave polarization makes the particle motion an oscillation along the group velocity direction, with zero first-order mean transport. The constant-phase lines of an internal gravity wave are also parallel to the group velocity. Their inclination $\vartheta$ above the horizontal satisfies $\sin \vartheta =\omega /N$.

When $\omega >N$, set $\kappa =k\sqrt{1-N^2/{\omega }^2}$. Boundedness at infinity selects the evanescent waveThe horizontal and vertical velocity components are in quadrature: fluid particles describe small ellipses, and the response decays over ${\kappa }^{-1}$. There is no upward time-averaged energy flux, because $P$ and $W$ are in quadrature. A real vertical group velocity is not defined for this evanescent wave. The pattern travels horizontally with phase velocity $\omega /k$.

For $N=0$ the response is the evanescent wave with $\kappa =k$; the ratios $\omega /N$ should not be used. At the cutoff $\omega =N>0$, the bounded harmonic solution has $W=-i\omega {\eta }_0$, $U=0$: it neither decays nor has nonzero upward group velocity. It is the limiting cutoff response, rather than a localized radiating disturbance. A uniform nonzero buoyancy frequency in an infinitely deep Boussinesq approximation is itself a local idealization of the background mass density.

The arrows distinguish the phase velocity, group velocity, and oscillatory particle motion of the internal gravity wave. The right panel shows the decay envelope and particle ellipses of the evanescent wave.

In a steady state, mass conservation makes $\mathcal{F}$ independent of radius. Taking the accretion rate $\dot{M}>0$ for inward flow gives $\mathcal{F}=-\dot{M}$. The conservation of angular momentum equation then makes $-\dot{M}h+\mathcal{G}$ constant. The zero-torque inner boundary condition fixes this constant to $-\dot{M}{h}_{\mathrm{in}}$, soThe viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^3\bar{\nu }\Sigma {\Omega }^{\prime }=2\pi qh\bar{\nu }\Sigma$. Combining these expressions gives the steady accretion disk with arbitrary rotation lawFor $x=r/R_g$, the Paczyński-Wiita circular orbit formulas giveConsequentlyIn particular, $f(6)=0$ and $f\to 1$ far from the black hole, recovering the outer Keplerian accretion disk.

Inside the innermost stable circular orbit, the gas enters the plunging region of a black-hole accretion disk. Its rapid inward motion leaves little time for stresses to exchange angular momentum, motivating the zero-torque inner boundary condition. This is a thin-disc approximation; a strong magnetic stress could change it.

Matter supplied from very large radius has negligible specific orbital energy, while matter crossing the inner edge carries ${\epsilon }_{\mathrm{in}}=-\eta c^2$. With zero inner torque, no energy is supplied by a stress at that edge. If the heat released outside it escapes by radiative transfer, rather than being lost through inward advection, the integrated conservation of energy balance gives

Let $b$ denote buoyancy in a stably stratified Boussinesq approximation. Geostrophic balance and hydrostatic balance imply the thermal wind relation $f_0{U}_z=-{\bar{b}}_y$. Consequently the specified vertical shear requiresFor northern-hemisphere $f_0>0$ and $\Lambda >0$, the mean potential temperature decreases northward and increases upward: ${\bar{\Theta }}_y=-({\Theta }_0/g)f_0\Lambda$, ${\bar{\Theta }}_z=({\Theta }_0/g)N^2$. The corresponding mean isentropes slope upward toward the north, with $dz/dy=f_0\Lambda /N^2$. Stable potential-temperature stratification does not require the actual temperature to increase upward; the dry adiabatic lapse also enters that comparison.

In the rigid-lid constant-shear baroclinic model, perturbation quasi-geostrophic potential vorticity isThe mean interior potential-vorticity gradient vanishes on this constant-$f_0$, constant-$N$, linear-shear basic state. Linear advection of that perturbation therefore becomes $({\partial }_t+\Lambda z{\partial }_x)q^{\prime }=0$, the supplied interior equation. It retains slow balanced motion and omits friction, planetary-vorticity gradients and nonlinear disturbance advection.

Substitution of ${\psi }^{\prime }=\Psi (z)\cos (ly)e^{ik(x-ct)}$ givesA possible branch has zero interior perturbation potential vorticity, soBoth supplied applications have $f_0>0$, so this equals the source's expression. For an unstable mode the nonzero imaginary part of the phase speed prevents a real critical level $\Lambda z=c$; for real speeds the displayed solution remains a valid zero-interior-PV branch away from, and smoothly through, any critical level. Taking the real part gives the physical disturbance.

The rigid lids enforce $w^{\prime }=0$. Linearized adiabatic buoyancy conservation isUse the geostrophic conventions $v^{\prime }={\psi }_x^{\prime }$ and $b^{\prime }=f_0{\psi }_z^{\prime }$. At a lid this gives $({\partial }_t+\Lambda z{\partial }_x){\psi }_z^{\prime }-\Lambda {\psi }_x^{\prime }=0$, henceThese are the no-normal-flow condition expressed through boundary buoyancy advection, not conditions of zero horizontal velocity. Boundary temperature anomalies can interact even though this solution has no interior perturbation PV, allowing baroclinic instability to extract mean available potential energy.

The first equation expresses horizontal incompressibility. The other two express steady balance of the Coriolis force, a geopotential gradient, linear drag and a prescribed zonal acceleration $F_u(y)$; nonlinear advection has been neglected. The parameter $r>0$ is an inverse drag time, so the frictional acceleration is $-r(u,v)$. On the equatorial beta plane, $f=\beta y$.

Use the streamfunction conventionTake ${\partial }_x$ of the meridional momentum equation minus ${\partial }_y$ of the zonal one. The mixed geopotential derivatives cancel, andHorizontal incompressibility removes the first term, leavingThe right side is the vertical curl of the imposed acceleration, with units ${\mathrm{s}}^{-2}$. The left side balances planetary vorticity advection and frictional removal of relative vorticity. The streamfunction has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$, since it describes horizontal velocity, rather than depth-integrated volume flux.

To return the broad interior transport, the western boundary current has opposite sign and characteristic peakThe leading exponential profile retains the smaller interior offset, $v(0)\simeq (F_0/\beta )(1-L/{\delta }_W)\simeq -52.1\mathrm{m}{\mathrm{s}}^{-1}$. Because $\gamma L\simeq 0.593$ is not very small, the exact two-wall solution gives a more accurate linear-model peak:All three evaluations identify an intense southward western boundary current, with peak at $x=0,y=0$, of order tens of metres per second. They are predictions of the stated steady linear barotropic model. Nonlinear advection, stratification and the actual basin geometry could alter them substantially; the calculation does not establish that real Pluto winds attain these values.

photon orbit and the Newtonian limit. Write $r_g=GM/c^2$ and $e=E/(m{c}^2)$. For a timelike circular Schwarzschild geodesic,The denominator tends to zero at $r=3r_g$, so the energy per unit rest mass diverges as a massive circular orbit approaches that radius. The corresponding local orbital speed measured by a static observer satisfies $v_{\mathrm{loc}}^2/c^2=r_g/(r-2r_g)$ and tends to one. This suggests a null circular orbit, but one must not simply insert $m=0$ into a massive-particle formula. For a photon, the radial effective potential is proportional to $(1-2r_g/r)/r^2$; its derivative vanishes at $r=3r_g$ and its second derivative there is negative. Thus the photon sphere is at $3GM/c^2$, and its circular null orbits are unstable.

At $r\gg r_g$, expansion givesThe second term is the total Newtonian circular-orbit energy, kinetic plus potential; the relativistic energy also includes rest energy. Here $E$ is the conserved energy referred to a stationary observer at infinity, not the local kinetic energy.

The stable inner edge. At fixed angular momentum, write $\ell =L/(mc)$. The timelike radial geodesic equation can be writtenThe circular condition $V_{\mathrm{eff}}^{\prime }=0$ gives ${\ell }^2=r_g{r}^2/(r-3r_g)$. Substituting into the second derivative taken at fixed $\ell$ givesTimelike circular orbits exist for $r>3r_g$; they are stable against radial perturbations for $r>6r_g$ and unstable for $3r_g<r<6r_g$. ThereforeThe boundary orbit is marginally stable. Differentiating $e(r)$ also gives a minimum there, consistently with this fixed-angular-momentum stability calculation.

accretion efficiency and black-hole spin. Matter supplied slowly from far away initially brings rest energy $d{m}_0{c}^2$. A radiatively efficient, zero-inner-torque thin disk releases the binding energy down to the innermost stable circular orbit; the remaining energy $e_{\mathrm{ISCO}}d{m}_0{c}^2$ is carried into the hole during the plunge. ThusHere $\eta$ denotes the efficiency called $ϵ$ in this question. The ideal relation assumes negligible additional dissipation during plunge and neglects photon capture and advection. A hot radiatively inefficient flow need not have this efficiency.

For a Kerr black hole, an equatorial prograde disk has a smaller innermost stable circular orbit and greater binding-energy efficiency than the Schwarzschild value; a retrograde disk has the opposite trend. The formal extremal prograde limits are $r_{\mathrm{ISCO}}\to r_g$ and $\eta \to 1-1/\sqrt{3}\simeq 0.423$. These are ideal limiting values, not a universal realized efficiency. photon capture and the accretion geometry affect the achievable black-hole spin and escaping radiation.

radiation force and the limiting luminosity. For isotropic luminosity $L$, the outward radiation momentum flux at radius $r$ is $L/(4\pi r^{2}c)$. In fully ionized hydrogen, Thomson scattering acts mainly on the electron, while electrostatic coupling makes the electron and proton move together. Neglecting the electron mass, one electron–proton pair experiencesEquating them yieldsThis Eddington luminosity uses electron-scattering opacity, spherical or isotropic force balance and negligible relativistic corrections at the force-comparison radius. Composition, other opacity, anisotropy and photon trapping can change the relevant limit.

Exponential growth and doubling. Distinguish rest-mass supply ${\dot{m}}_0$ from black-hole mass growth $\dot{M}$:At the Eddington limit,Using the printed scattering cross section gives $t_E\simeq 4.50\times 10^8\mathrm{yr}$. The Salpeter time $t_S$ is an e-folding time; the Eddington-limited doubling time isfor the Schwarzschild thin-disk efficiency. If one approximates $\dot{M}\simeq {\dot{m}}_0$, the factor $1/(1-\eta )$ is dropped and the estimate is about $1.79\times 10^7\mathrm{yr}$. The retained-energy version is the mass-conserving expression.

The efficiency bound. A growth factor $10^8$ requires $N=\log (10^8)=18.4207$ e-folds. Eddington-limited accretion is the fastest growth allowed by the adopted luminosity bound at fixed efficiency, so growth within a time $\Delta t$ requiresThis efficiency upper bound from Eddington-limited growth is an upper bound: higher efficiency radiates more of each supplied rest mass, and the luminosity ceiling therefore allows less mass to enter per unit time. A duty cycle below unity or sub-Eddington luminosity tightens it. The approximate convention $\dot{M}={\dot{m}}_0$ would instead give $\eta \lesssim 0.121$; both figures refer to the same effect, with different treatment of the radiated mass deficit. Mergers, super-Eddington growth or a different seed history would change the premises.

Coherent prograde accretion spins up the hole, moves the inner edge inward and increases $\eta$. For fixed $L/L_{\mathrm{Edd}}$ the mass-growth e-folding time then grows, so the logarithmic growth curve becomes less steep; one cannot simply retain the initial Schwarzschild efficiency throughout. Episodic misaligned accretion can have a different black-hole spin history. In general,

energy available to affect the galaxy. Let $M_g$ denote galaxy mass, keeping it distinct from black-hole mass. At constant efficiency, negligible seed mass and no merger contribution, the emitted energy in building $M_{\mathrm{BH}}=0.002M_g$ isFor $v=300\mathrm{km}{\mathrm{s}}^{-1}$ and the nonspinning efficiency this is approximately $121$. The common small-efficiency approximation gives $2000\eta \simeq 114$, and a fiducial $\eta =0.1$ would give about $222$ with retained mass included. These differences do not alter the conclusion: the available accretion energy substantially exceeds the stated galactic binding-energy estimate.

Only a coupling fraction of order one percent would then be energetically sufficient to affect the gas strongly. radiation, winds or jets can heat gas, drive outflows and suppress further cooling and star formation, providing active-galactic-nucleus feedback and possible self-regulation of black-hole growth. energy availability alone does not prove that the whole galaxy is expelled: most radiation can escape, jets can be collimated, heated gas can cool, and collisionless stars couple differently from gas.

mass from motion. A resolved Keplerian orbit gives $M=4{\pi }^2{a}^3/(G{P}_{\mathrm{orb}}^2)$ when the orbiting object has negligible mass. At the Galactic centre, repeated infrared astrometry and spectroscopy measure sky motion, acceleration and line-of-sight velocity of stars around Sagittarius A*. Resolving a compact enclosed mass and fitting complete orbits avoids reliance on a single projected velocity. Source confusion, reference-frame stability, distance, inclination and possible extended mass must be included in the fit.^[ref]

For an unresolved stellar-mass binary in a circular orbit, the donor's velocity semiamplitude $K_2$ and orbital period give the binary mass functionIt provides a lower mass bound; donor mass and orbital inclination are needed for an actual mass estimate. For an eccentric orbit, multiply the observational combination $P_{\mathrm{orb}}{K}_2^3/(2\pi G)$ by $(1-e^2{)}^{3/2}$ when evaluating the binary mass function. Resolved gas or maser rotation and stellar-dynamical models also measure masses, provided pressure support, noncircular motion and line-of-sight projection are modeled. Resolving the black-hole sphere of influence improves separation of central mass from the host galaxy.

For active nuclei, reverberation mapping measures a lag $\tau$ between variable continuum and broad-line response. With a line velocity scale $\Delta v$,The factor $f_{\mathrm{geom}}$ accounts for geometry and kinematics and is not known simply from the lag. Cadence, line responsivity, inclination and winds require control. Rapid variability demonstrates compactness but is not by itself a unique mass measurement.

black-hole spin from strong-field emission. Define $a_{\ast }=cJ/(G{M}^2)$. Large-radius Keplerian motion is principally a mass probe and usually insensitive to black-hole spin. In the Kerr model, the radius and binding energy of the innermost stable circular orbit vary with $a_{\ast }$ and orbital orientation, making inner-disk radiation useful.

black-hole continuum fitting infers an inner emitting radius from the disk's thermal spectrum and flux normalization, then identifies it with the innermost stable circular orbit in a relativistic thin-disk model. Independent mass, distance and inclination, together with atmospheric spectral hardening, absorption and contamination, are needed. Spectra should be selected in a disk-dominated regime where truncation and strong advection are limited. A 2011 Cygnus X-1 analysis combined such spectra with independently constrained binary parameters and reported a near-extreme black-hole spin conditional on black-hole spin–orbit alignment and its disk model.^[ref]

black-hole reflection spectroscopy instead models a relativistically broadened iron-line and reflection continuum: orbital Doppler shifts and beaming broaden the profile, while gravitational redshift contributes to the extended red wing. Its innermost emitting radius can constrain black-hole spin if reflection extends to the stable inner edge. Ionization, abundance, absorption, coronal illumination, disk thickness and possible truncation are competing parameters; a broad line is not a direct black-hole spin reading independent of these assumptions. Timing or precession models can add information but introduce their own dynamical identification.

mass is principally inferred from gravitating dynamics; black-hole spin is inferred from strong-field effects with a specified spacetime and accretion model. Agreement between independent methods is especially valuable because disk inner radius, inclination and radiation-transfer assumptions can otherwise mimic a change of black-hole spin.

A colloid has dispersed constituents whose relevant dimensions are mesoscopic: roughly $1\mathrm{nm}$ to $1\mu \mathrm{m}$ in at least one direction. Interfaces and Brownian motion are consequently important. The size range is a useful convention, not a sharp boundary separating all physical behaviour. A sol of fine silicon dioxide particles in water illustrates solid in liquid; an oil-in-water emulsion illustrates liquid in another immiscible liquid; and a foam illustrates gas dispersed in a liquid, with thin intervening films supplying the small structural dimension. A colloid need not be a thermodynamic equilibrium state.^[ref]

Unprotected particles attract through Van der Waals forces. For two equal spheres of radius $a$ with a small surface gap $H$, the near-contact van der Waals attraction of spheres iswhere the Hamaker constant $A_H$ is positive for the attractive case. The large number of interacting constituent molecules can give a well much deeper than $k_{B}T$, even though individual molecular attractions are weak. Brownian motion brings particles together, and a sufficiently deep attractive well makes separation unlikely. This gives particle aggregation; solid particles can remain distinct within a cluster rather than merging into one continuous solid. Emulsion droplets can additionally coalesce if the separating film ruptures.

Charge stabilization of colloids supplies a repulsive barrier. Like-charged surfaces carry compensating counterions and electric double layers; bringing the layers together costs free energy. DLVO theory combines this repulsion with Van der Waals forces. A barrier large compared with $k_{B}T$ can make aggregation slow, even when sticking would lower the final free energy. Adding salt generally shortens the Debye–Hückel screening length and can lower the barrier; surface charge regulation and specific ion interactions also matter. A negatively charged particle and a positively charged particle need not repel.

Steric stabilization instead coats each particle with adsorbed or grafted polymer chains. In a good solvent, overlap and compression of these layers raise the osmotic pressure and reduce chain entropy, opposing close contact. Electrosteric stabilization combines the two mechanisms. Dense, securely bound layers matter: incomplete adsorption can allow polymer bridging, while nonadsorbing dissolved polymers can create an attractive depletion force. Poor solvent quality can collapse a protective layer. Reducing the density mismatch slows sedimentation but does not remove the interparticle attraction, so uniform suspension and resistance to aggregation are different requirements.

The left panel is an original qualitative potential with arbitrary coefficients: it illustrates a barrier and its reduction by shorter-ranged repulsion, without fitting an experimental DLVO theory curve or extrapolating the continuum attraction to molecular contact. The right panel plots the derived quasi-steady concentration profile in dimensionless variables.

In diffusion-controlled growth, attachment is rapid and delivery through the surrounding medium is the bottleneck. The growing surface depletes the neighboring solute, so a substantial concentration gradient develops. In reaction-controlled growth, arrival is fast compared with incorporation; many arrivals fail to attach, and the surface concentration stays near the reservoir value. The analogous distinction for particle aggregation is between near-unit sticking probability in diffusion-limited cluster aggregation and barrier-suppressed sticking in reaction-limited cluster aggregation. Here the subsequent calculation concerns accretion onto one growing particle, rather than a collision kernel for two mobile clusters.

The printed time-dependent equation is the diffusion equation, not the usual Poisson equation. For constant diffusion coefficient $D$, no bulk sources and negligible advection, Fick's laws give the number flux $\mathbf{j}=-D\nabla c$. Local conservation of mass, expressed here as conservation of the accreting constituent, then givesAn arbitrary time-independent exterior profile does not solve the full moving-boundary problem exactly. The required elementary result uses a quasi-steady approximation: diffusion establishes the exterior profile much faster than the particle radius changes. At each instantaneous radius $a$, set $\partial c/\partial t\simeq 0$. Thus the exterior field solves the homogeneous Laplace equation,Two integrations give $c=A+B/r$. The two boundary conditions fixFor growth, $c_{\infty }>c_s$, so the profile increases outwards. Its radial flux is negative, and the positive inward arrival rate isA material density is needed to turn an arrival rate into a radial velocity. Introduce the constant constituent number density $c_p$ inside the particle, using the same units as $c_s,c_{\infty }$. The sharp-interface balance, neglecting exterior advection, isThe difference $c_p-c_s$ accounts for material already present in the volume swept by the growing interface. In the usual dilute-feed approximation $c_s\ll c_p$, it can be replaced by $c_p$. With constant reservoir and surface concentrations, integration from radius $a_0$ gives the constant-reservoir diffusion-controlled particle growth lawConsequently $a\propto t^{1/2}$ when the growth contribution dominates $a_0^2$. For a finite seed at the literal limit $t\to 0$, the increment $a-a_0$ is instead initially proportional to $t$. Here “early time” means before appreciable reservoir depletion or interactions with other particles, while already beyond the rapid initial establishment of the diffusion field. Indeed, the ratio of diffusion-adjustment time $a^2/D$ to radial-growth time $a/\dot{a}$ is $(c_{\infty }-c_s)/(c_p-c_s)$; small supersaturation relative to the particle density makes the quasi-steady approximation self-consistent. No initial concentration field was specified, so an exact transient solution is not uniquely determined.

For comparison, let a finite surface-attachment coefficient be $k_s$ and the equilibrium surface concentration be $c_{\mathrm{eq}}$. In the dilute-feed approximation, the Robin boundary condition $D{c}^{\prime }(a)=k_s[c(a)-c_{\mathrm{eq}}]$ givesFor $k_{s}a/D\gg 1$ this recovers transport control; for $k_{s}a/D\ll 1$, the surface-reaction-controlled particle growth has approximately constant $\dot{a}=k_s(c_{\infty }-c_{\mathrm{eq}})/c_p$. Thus reaction control gives $a-a_0\propto t$, in contrast to the diffusion-controlled square-radius law. No numerical $c_p$ or $k_s$ is supplied or needed for these scaling results.

velocity $U$, length $\ell$ and diffusivity $D$ have dimensions $L/T$, $L$, and $L^2/T$. Their dimensionless Péclet number isdiffusion dominates for $\mathrm{Pe}\ll 1$, advection for $\mathrm{Pe}\gg 1$. At fixed $U,D$, the crossover length is $D/U$; neglecting advection requires $\ell \ll D/U$.

TODO: user-deferred organism-specific swimming-speed and detection-scale estimates. No biological numerical substitution is made.

For the generic equilibrium transport derivation, let an external potential $V(x)$ give force $F=-V^{\prime }$ and drift $u=F/\zeta$. Boltzmann equilibrium has $c\propto e^{-V/(k_{B}T)}$, hence $c^{\prime }=Fc/(k_{B}T)$. Substituting in the zero-current condition givesThe Stokes–Einstein relation in the final equality is for an isolated no-slip sphere in a Newtonian liquid at low Reynolds number. Zero flux alone is insufficient: the equilibrium distribution and mechanical mobility law are also needed.

TODO: user-deferred organism-specific absorbing-surface calculation and its biological flux interpretation.

TODO: user-deferred receptor-number, receptor-density and receptor-competition calculations. These continuations have not been solved.

Let the channel length be ${\ell }_c\gg R$ and assume uniform solution conductivity. Its Ohmic resistance scales as ${\ell }_c/(\sigma R^2)$, whereas access resistance of a reservoir mouth scales as $1/(\sigma R)$. The fraction of voltage in the entrance regions is therefore of order $R/{\ell }_c$. Most of the potential drop lies inside the long channel; the small entrance field region is neglected in the model, not asserted identically zero.

Write $x$ for inserted contour length and choose $qE>0$ to pull the chain inward. The charge already inside feels force $qEx$, givingConfinement removes polymer configurations. In the stiff-chain regime, the confinement deflection length obeys $R^2\sim {\lambda }^3/{\ell }_p$, hence $\lambda \sim (R^2{\ell }_p{)}^{1/3}$. A free-energy cost of order $k_{B}T$ per deflection gives a roughly constant resisting entropic force $f\sim k_{B}T/\lambda$. The insertion barrier of a charged polymer follows from $F(x)\simeq fx-qE{x}^2/2$, with an entry maximum at $x_{\ast }=f/(qE)$ when that point lies within the insertion interval. Its height is $f^2/(2qE)$. The field supplies work against this confinement entropy.

Put ${\zeta }_{\ell }=2\pi \eta /\ln (R/a)$, so $\gamma (x)={\zeta }_{\ell }x$. After capture, neglect the entry barrier, exterior-chain drag and thermal corrections to strong directed motion. ThenThe confined-polymer entry and recoil times have different length dependence; this entry graph is linear through the origin. Charge and drag both grow with inserted length, leaving a constant velocity. The same cancellation explains length-independent mobility in an ideal free-draining electrophoresis model. It does not include a capture waiting time or claim that every free-solution chain has identical hydrodynamic drag.

A measured superlinear entry time can arise when exterior-coil drag, tension propagation, wall adsorption or nonuniform friction adds length dependence absent from this model. For example $\gamma (x)\propto x^p$ with the same $qEx$ drive gives $\tau \propto L_0^p$ for $p>0$. Under the assumptions of uniform wall charge and polymer entry scaling, a homogeneous thin-double-layer electro-osmotic model adds a length-independent fluid velocity $u_{\mathrm{eo}}=-ϵ{\zeta }_{w}E/\eta$. It changes the drift coefficient, not its length exponent. This does not rule out charge heterogeneity, adsorption or near-cancellation followed by diffusion in a real channel; those are additional mechanisms.

With the field off, the confinement entropy drives outward motion. Assuming the constant force $f$ and the same local friction,The small final end region lies outside the long-confined-segment approximation but gives only a small correction. Thus ideal insertion scales as $L_0$, while recoil scales as $L_0^2$.

TODO: user-deferred DNA-specific numerical stopping-voltage estimate and its biological parameter-based inference. No biological quantities are evaluated.

For entropic-trap length separation, a reservoir/channel array can repeat length-dependent capture, confinement and relaxation delays, turning them into different mean residence times. The uniform inside-channel drift speed derived above is itself length independent, so that drift alone cannot separate lengths. Separation must retain length dependence of the barrier or of the capture/recoil dynamics and can be tuned with a forward-biased time-dependent field. For large bulk throughput, pressure-driven advection is preferable to diffusion alone; a pressure gradient can be combined with electric capture control. In the unobstructed continuum limit $Q=\pi R^4\Delta P/(8\eta {\ell }_c)$. Excessive driving erases trapping differences, so throughput and resolution compete; the question supplies no geometry or adsorption law selecting a unique optimal drive.

Define $A=2\cot \alpha /3+8\mathrm{Re}/15$. The linearization about $h=1$ givesThe Galilean transformation $\xi =x-2t$, $\tau =t$ has ${\partial }_t={\partial }_{\tau }-2{\partial }_{\xi }$ and ${\partial }_x={\partial }_{\xi }$. It removes the leading advection, leavingFor a normal mode proportional to $e^{\sigma \tau +ik\xi }$, the second and fourth derivatives contribute $-k^2$ and $k^4$. The dispersion relation is thereforeFor $A>0$ and $S>0$, the unstable wavenumbers satisfy $0<|k|<\sqrt{A/S}$. In the laboratory coordinates, the temporal exponent is $\sigma -2ik$, so this linear stability calculation predicts growth together with downslope transport at dimensionless speed $2$.