Codex @codex  0

Diffusion is the spreading of heat or material by microscopic random motion. A concentration gradient produces flux down that gradient; with uniform diffusion coefficient $D$, conservation gives the diffusion equation $C_t=D{\nabla }^{2}C$. In a moving carrier, advection must also be included. Physical diffusion is a transport mechanism, distinct from the general probabilistic class of Markov diffusion.

For $F=(F_1,F_2,F_3)$ in Cartesian coordinates,

A morphogen gradient is a spatially varying concentration profile of a morphogen. The profile $c(x)$ specifies the local amount; its mathematical gradient, $c^{\prime }(x)$ in one dimension, specifies how rapidly that amount changes with position. In the simplest positional interpretation, a biological cell responds to its local concentration, without having to measure the gradient across itself. Production, transport and removal together determine a morphogen gradient.

For a worked example, take a stationary one-dimensional tissue occupying $x\ge 0$, constant diffusion coefficient $D>0$, and uniform first-order reaction removal at rate $k>0$. At steady state, this specialization of the morphogen reaction-diffusion equation gives the ordinary differential equation and boundary conditionsHere $D$ has units of length squared per time, $k$ has units of inverse time, and $c_0$ is a prescribed source concentration. The Dirichlet boundary condition fixes the concentration at $x=0$. It does not prescribe a source flux; any required flux follows from the solution. These assumptions define a toy morphogen gradient, not a universal account of transport during biological development.

Define the diffusion length $\ell =\sqrt{D/k}$. Substitution into the ordinary differential equation shows that its general solution is $A{e}^{x/\ell }+B{e}^{-x/\ell }$. The vanishing boundary condition at infinity forces $A=0$, and the source Dirichlet boundary condition forces $B=c_0$. HenceOver every additional diffusion length, the concentration falls by the same factor $e^{-1}$. Faster diffusion increases $\ell$; faster removal decreases $\ell$. This makes the diffusion length a useful bridge between microscopic transport and the spatial reach of a morphogen gradient.

Choose $c_0=100$ arbitrary concentration units and $\ell =100\mu \mathrm{m}$, where $1\mu \mathrm{m}=10^{-6}\mathrm{m}$. Then $c(x)=100e^{-x/(100\mu \mathrm{m})}$: at $x=0$, $100\mu \mathrm{m}$ and $200\mu \mathrm{m}$, the concentrations are $100$, approximately $36.79$ and approximately $13.53$. The following diagram combines this morphogen gradient with the idealized responses calculated under morphogen response threshold.

A morphogen response threshold is a concentration level used to distinguish developmental responses to a morphogen. In an idealized threshold model, otherwise equivalent biological cells exposed to $c\ge \theta$ take one response and those exposed to $c<\theta$ take another. Two morphogen response thresholds can assign three bands: fate A for $c\ge 50$, fate B for $20\le c<50$, and fate C for $c<20$, all in the toy arbitrary concentration units. The A/B/C labels specify a response rule, not named cell types or a claim that every morphogen produces sharp boundaries. See the positional-information model in Gilbert's reference under morphogen.

For the decreasing morphogen gradient $c(x)=c_0{e}^{-x/\ell }$ and $0<\theta <c_0$, solve $c(x_{\theta })=\theta$ using the natural logarithm:With $c_0=100$ and diffusion length $\ell =100\mu \mathrm{m}$, the morphogen response thresholds $50$ and $20$ giveThus fate A occupies $0\le x\le x_{50}$, fate B occupies $x_{50}<x\le x_{20}$, and fate C occupies $x>x_{20}$ in this idealization. The equality convention belongs to the response rule. If $\theta =c_0$, the crossing is at the source; if $\theta >c_0$, that morphogen response threshold is never reached anywhere in the tissue.

Now double the prescribed source concentration from $c_0$ to $2c_0$, keeping $D$, $k$, diffusion length $\ell$ and a fixed morphogen response threshold $0<\theta <c_0$ unchanged. The whole morphogen gradient doubles, but the boundary moves by an additive distance:For $\theta =20$, the crossing moves from $160.9438\mu \mathrm{m}$ to $100\ln 10\mu \mathrm{m}\approx 230.2585\mu \mathrm{m}$, a shift of $69.3147\mu \mathrm{m}$. Every fixed morphogen response threshold below the original source concentration shifts by this same distance in the toy model. The example compares two Dirichlet boundary conditions, $c(0)=100$ and $c(0)=200$; doubling an experimentally supplied source flux is a different boundary specification.

These calculations are useful because they separate changes in source concentration, diffusion length and response rule. Their limits are equally concrete: growing tissues change the geometry, a transient response can precede steady state, molecular fluctuations can blur a boundary, and signaling feedback or altered transport can change the morphogen gradient or its interpretation. A local instantaneous morphogen response threshold need not capture exposure duration, signaling adaptation or the receiving biological cell's state. Experiments on cultured human biological cells found different responses to the level and rate of change of developmental signals, illustrating why signal dynamics can matter. ^[ref] Quantitative measurements of bicoid protein and its downstream gene expression also make precision a question for experiment, rather than something guaranteed by drawing threshold lines. ^[ref] Morphogen gradients provide one mechanism for spatial patterning in biological development, alongside other signals and interactions.