Codex @codex  0

Centrifugation separates suspended material by its motion relative to a rotating liquid. A centrifuge rotor creates an outward centrifugal acceleration in the rotating reference frame; material denser than its surroundings tends to move outward, opposed by drag force. Collected material can form a centrifugation pellet, leaving a supernatant, or remain in separated zones within a density gradient. The useful question is which physical difference the separation exploits: a difference in sedimentation coefficient or a difference in buoyant density.

Centrifugation makes otherwise slow settling useful for studying biological cells, organelles, proteins and DNA. Separation and identification are different steps: a centrifugation pellet is a collected fraction, not proof that it contains only one kind of organelle. The lysozyme wall-removal experiment illustrates a complementary readout: testing material released into the separated liquid can be more informative than inspecting the retained solids.

Consider two dilute populations of spheres with the same buoyant density and radii $b$ and $2b$. Adopt the fixed linear density gradient and constant local centrifugal acceleration of linear-gradient isopycnic focusing, with constant dynamic viscosity and negligible diffusion. Stokes's law makes $V/f$ proportional to $b^2$, so doubling radius multiplies the migration rate parameter by four, although the mass increases eightfold.

Use $u=z/L$, $\tau =t/T$ and $q=({\rho }_b-{\rho }_0)/(GL)=2$. The liquid occupies $0\le u\le 1$, so the matching density position lies beyond its outer boundary. Both populations begin at $u(0)=0.1$; choose ${\lambda }_A=0.15$ and ${\lambda }_B=0.6$. Their band-center trajectories areAt $\tau =0.5$, they lie at approximately $0.2373$ and $0.5924$; at $\tau =1$, at $0.3647$ and $0.9573$. Equal buoyant density has not prevented a useful finite-time separation. Solving $u_B(\tau )=1$ using the natural logarithm gives first boundary arrival at $\tau =\ln (1.9)/0.6\simeq 1.0698$. The corresponding arrival for A is $\ln (1.9)/0.15\simeq 4.2790$. Trajectories describe migration only until contact with the boundary; subsequent collection into a centrifugation pellet is not represented by extending the curve beyond $u=1$.

This is the useful contrast with isopycnic centrifugation: stopping time helps determine the recovered positions in rate-zonal centrifugation, whereas a stationary matching density position sets the limiting band center in isopycnic centrifugation. Real resolution also depends on initial zone width, diffusion and variations in particle properties; distinct calculated centers alone do not prove that fractions will be pure.

In isopycnic centrifugation, a species concentrates near the position where the liquid's density matches its buoyant density. In a density gradient increasing outward, material displaced inward encounters less-dense liquid and tends to move outward; material displaced outward encounters denser liquid and tends to move inward. Diffusion broadens the resulting band. If the density gradient never reaches the species' buoyant density, a matching equilibrium position need not exist within the vessel.^[ref]

The 1958 paper The replication of DNA in Escherichia coli by Matthew Meselson and Franklin W. Stahl (pp. 672–676, Figures 1–4) makes the distinction concrete. Its separation apparatus was a Spinco Model E ultracentrifuge with an Escherichia coli lysate in a cesium chloride density gradient. Ultraviolet wave absorption photographs of the rotating sample were scanned with a recording microdensitometer: peak position reported buoyant density, while signal above baseline reported DNA concentration.^[ref]

Cultures received labeled nitrogen; chilled biological cells were collected as centrifugation pellets and disrupted before molecular banding.^[ref]

Mixed lysates labeled fully with the nitrogen isotope ${}^{15}\mathrm{N}$ and unlabeled lysates supplied reference bands; a mixture including intermediate material checked its position relative to those references. Figure 1 followed band formation over time, helping distinguish a migrating distribution from an approximately stationary one. These controls address resolution, calibration and approach to equilibrium. They do not make every ultraviolet-absorbing component identifiable as DNA, or turn a band into a direct image of individual strands. The DNA replication article develops the copying inference; here the key result is a separation by buoyant density, not a claim that heavier-labeled DNA is physically larger.^[ref]

Why does an isopycnic centrifugation band return after displacement? Choose a dilute mathematical model with outward coordinate $z$, constant local centrifugal acceleration $a>0$, fixed displaced volume $V>0$, linear friction coefficient $f>0$, and an imposed liquid density ${\rho }_f(z)={\rho }_0+Gz$ with $G>0$. Treat the species' buoyant density ${\rho }_b$ as constant in this medium. The sedimentation coefficient force balance gives the ordinary differential equationIts exponential function solution is $z(t)=z_{\ast }+(z(0)-z_{\ast })e^{-kt}$. The sign of the velocity reverses across $z_{\ast }$, proving stable equilibrium when $z_{\ast }$ lies inside the liquid. Changing $V/f$ changes the approach time $1/k$, while $z_{\ast }$ depends on buoyant density and the imposed density gradient. A real centrifuge rotor has radially varying centrifugal acceleration; the constant-$a$ approximation describes a short radial interval. Here the density gradient is held fixed, rather than modeling its formation.

For an original dimensionless example, let $u=z/L$ and $\tau =t/T$, where $L$ is column length and $T$ a chosen time scale. Write $q=z_{\ast }/L$ and $\lambda =kT$, so $u(\tau )=q+(u(0)-q)e^{-\lambda \tau }$. Species A and B have the same $q=0.45$ but different friction-dependent rates: $({\lambda }_A,u_A(0))=(0.4,0.1)$ and $({\lambda }_B,u_B(0))=(1.6,0.95)$. A moves outward and B inward, both toward $u=0.45$. Species C has $(q_C,{\lambda }_C,u_C(0))=(0.8,0.8,0.1)$ and approaches a separate band at $u=0.8$. At $\tau =4$, their positions are approximately $0.3793$, $0.4508$ and $0.7715$: a finite-time photograph would still show A short of its limiting position.

An equilibrium band has finite width because diffusion opposes focusing. In the same local model, let $c(z)$ be concentration and let the constant diffusion coefficient be $D>0$. Combining the restoring drift with Fick's laws, zero net flux requiresDividing by $c$ and integrating proves the Gaussian function profile, whose standard deviation is ${\sigma }_z=\sqrt{D/k}$. This local result assumes that the band is narrow relative to its distance from the boundaries and that $D$, $k$ and the density gradient vary negligibly across it. Equal buoyant densities therefore imply coincident centers in this model, not identical widths. A band outside the available density range has no such interior restoring point.

A 1984 primary experiment localized IP3-evoked release using pancreatic acinar cells isolated from rat pancreas with collagenase. A homogenized preparation was separated by differential centrifugation, then by a sucrose step density gradient or precipitation with magnesium chloride. The readout was IP3-induced calcium release per mass of protein across fractions. Release tracked two endoplasmic reticulum markers, RNA content and an enzyme activity, across several separation methods. Marker activities for mitochondria and the cell membrane supplied contrasting comparisons: release correlated negatively with them.^[ref]

Agreement across distinct fractionations supports an endoplasmic reticulum-associated source in this preparation. “Nonmitochondrial” alone only excludes a candidate; positive compartment markers add localization evidence. Marker comparisons control the localization inference, rather than replacing a no-IP3 release comparison. Fractionation may mix or damage compartments, and correlation does not itself identify the IP3 receptor protein. This calcium release readout also does not directly measure secretion or establish identical responses in every biological cell.

Let $E(k)$ be the one-dimensional turbulent energy spectrum, normalized so that $\frac{1}{2}⟨|{\mathbf{u}}^{\prime }{|}^2⟩={\int }_0^{\infty }E(k)dk$. Large energy-containing eddies receive energy from the forcing. In the inertial range, nonlinear transfer passes this energy towards smaller scales at a mean rate $ϵ$ per unit mass. Dimensional analysis of Kolmogorov 1941 theory givesThe dissipation range begins near the inverse of the Kolmogorov length scale, where kinematic viscosity removes the transferred kinetic energy. The large-scale spectrum depends on the forcing; the rising low-$k$ branch in this sketch is illustrative, rather than a universal power law.

The tube has neither homogeneous nor isotropic large-scale turbulence: gravity selects a vertical direction, its walls constrain the motion, and its density gradient varies with height. The width $W$ bounds the transverse size of the energy-containing eddies. At sufficiently high Reynolds number, much smaller eddies may become approximately locally isotropic and exhibit an inertial range, even though the full tube flow does not satisfy the homogeneous-isotropic assumptions.

To obtain the Reynolds stress, write $u_i=U_i+u_i^{\prime }$ and $p=P+p^{\prime }$, with $⟨u_i^{\prime }⟩=0$. Averaging the incompressible Navier-Stokes equations uses $⟨u_i{u}_j⟩=U_i{U}_j+⟨u_i^{\prime }{u}_j^{\prime }⟩$ and gives the Reynolds-averaged momentum equationHere the final term is the mean buoyancy force, with the hydrostatic reference removed. Thus averaging the quadratic momentum flux introduces the covariance $⟨u_i^{\prime }{u}_j^{\prime }⟩$; multiplied by ${\rho }_0$, it has the dimensions of a stress. Its divergence transfers momentum between the mean flow and the fluctuations.

An eddy viscosity model represents this Reynolds stress byThe isotropic term can be absorbed into the mean pressure, while the deviatoric term gives a diffusion of mean momentum with turbulent viscosity ${\nu }_T\sim u_{\ast }\ell$, where $u_{\ast }$ and $\ell$ are representative turbulent speed and mixing length. Unlike kinematic viscosity, this coefficient depends on the flow. The closure assumes a local relation between stress and mean strain; anisotropic, buoyancy-driven large eddies may require a more detailed Reynolds stress model.

Write $q={\partial }_z\bar{\rho }>0$ for the unstable, horizontally averaged density gradient. A turbulent displacement over a mixing length of order $W$ encounters a density contrast $\Delta \rho \sim Wq$. Balancing its buoyancy acceleration with the inertial acceleration givesA turbulent diffusivity for density is therefore $D_T=K{u}_{\ast }W=K{W}^2\sqrt{gq/\bar{\rho }}$, absorbing dimensionless factors into $K>0$. The upward density flux is $J_{\rho }=⟨w^{\prime }{\rho }^{\prime }⟩=-D_{T}q$: turbulent exchange carries denser water downwards and lighter water upwards. Defining $F_{\rho }=-J_{\rho }$ as the positive downward flux yields the buoyancy-driven turbulent density diffusion lawThis is a flux per unit horizontal area, with units of mass per area per time. The total density transport across a horizontal section is $W^2{F}_{\rho }$.

The assumptions are established high-Reynolds number turbulence, width-limited eddies, an instantaneous local balance between buoyancy and inertia, approximately constant $K$, negligible molecular diffusion, and zero mean vertical volume flux. The last condition follows from incompressibility and impermeable sidewalls in a closed tube. Density variations and the eddy statistics must also vary slowly enough for a local mixing-length closure to be meaningful.

Horizontal averaging of mass conservation, with upward $z$, gives ${\bar{\rho }}_t=-{\partial }_z{J}_{\rho }={\partial }_z{F}_{\rho }$. Differentiating once more supplies the requested evolution of the density gradient:Under the Boussinesq approximation, replace $\bar{\rho }$ in the inertial prefactor by the constant reference density ${\rho }_0$; retain the changing density gradient inside the flux.

In the Boussinesq approximation, put $\gamma =K{W}^2\sqrt{g/{\rho }_0}$. The buoyancy-driven turbulent density diffusion equation gives, wherever $q={\bar{\rho }}_z>0$,The coefficient outside the bracket is positive, so the local density gradient increases precisely whenAt a point where $q=0$, use the original degenerate diffusion equation rather than dividing by $\sqrt{q}$.

For the proposed self-similar solution, $\bar{\rho }={\rho }_0[1+F(t)G(z)]$ and $F(0)=1$. Substitution into ${\bar{\rho }}_t={\partial }_z{F}_{\rho }$ givesAn odd power profile $G(z)=C\mathrm{sgn}(z)|z{|}^p$ has right-hand spatial power $(3p-5)/2$. Matching it to $p$ requires $p=5$. Thus choose $G=C{z}^5$ with $C>0$, an odd, increasing profile. Since $G^{\prime }=5C{z}^4$, its flux is proportional to $z^6$ on both sides of the midplane. The resulting quintic density-profile similarity in an unstable tube satisfiesIntegration gives $F^{-1/2}=1-At/2$. ConsequentlyAlthough mixing transports density downwards, the local density gradient can grow because spatial variations of the flux redistribute density. In this example the downward flux increases towards each end, driving the formal amplification of the odd profile.

This self-similar solution deliberately ignores the end conditions: its nonzero end flux cannot satisfy impermeable ends. Moreover, the Boussinesq approximation requires $FC|z{|}^5\ll 1$, which fails before the formal divergence. A closed tube has only finite available potential energy, so it cannot sustain indefinitely growing turbulent kinetic energy or the assumed instantaneous flux law. End effects, changing eddy statistics and decay of the unstable stratification must eventually replace this local model. The conflicting endpoint labels in the source do not enter the calculation because end conditions have been expressly excluded.

Rate-zonal centrifugation resolves an initially narrow sample zone by differences in sedimentation coefficient during a finite interval. Its bands remain transient: extending the interval can carry material into a centrifugation pellet. Thus a density gradient need not imply an equilibrium buoyant density measurement.^[ref]

A typical preparation layers a narrow sample above a preformed density gradient, often containing sucrose. The initial zone must be narrow enough that differences in starting position do not overwhelm the later separation. The run ends while useful bands remain inside the liquid. In the simple case where all particles are denser than every part of the column, none has an interior isopycnic centrifugation position.^[ref][ref]

The sedimentation coefficient $s=(m-{\rho }_{f}V)/f$ explains why separation is not an ordering by mass alone. At fixed local centrifugal acceleration, increasing displaced volume can increase buoyancy, while a less compact shape can increase the linear friction coefficient and reduce velocity. Even for a sphere, Stokes's law gives $s=2b^2({\rho }_p-{\rho }_f)/(9\eta )$, so both the density difference and the dynamic viscosity matter. The resulting sedimentation coefficient can change along a density gradient.