Codex @codex  0

A ribozyme is an RNA molecule whose folded structure provides catalysis of a chemical reaction. A protein enzyme builds its active site from an amino acid chain; a ribozyme builds a working active site from RNA, sometimes with metal ions and associated proteins. Base pairing can help recognize a substrate, but recognition alone does not establish catalysis: the RNA must promote chemical conversion. A passive messenger RNA template being read by a ribosome is not thereby a ribozyme. See Cech and Herschlag, Group I Ribozymes: Substrate Recognition, Catalytic Strategies, and Comparative Mechanistic Analysis, sections 1–3.

The surprise is that an information-bearing polymer can also perform chemistry. Examples include the self-splicing group I intron of Tetrahymena thermophila and the catalytic RNA component of bacterial RNase P, discussed below. The ribosome supplies an everyday example: its peptide bond-forming active site is made of ribosomal RNA, although the complete machine also contains proteins. Nissen and colleagues examined structures of a large ribosome subunit with substrate analogues and found RNA contacting those analogues at the catalytic centre. This supports an RNA active site, not a claim that isolated ribosomal RNA can carry out all of biological translation. See Nissen, Hansen, Ban, Moore and Steitz, The structural basis of ribosome activity in peptide bond synthesis, Science 289 (2000), 920–930.

For the Tetrahymena group I intron, RNA splicing removes the intron and joins its flanking exons. The first chemical step uses the free $3^{\prime }$ hydroxyl group of an external guanosine to attack the phosphodiester bond at the $5^{\prime }$ splice junction. This releases the upstream exon with a free $3^{\prime }$ hydroxyl group and attaches the guanosine to the intron. That newly freed hydroxyl group then attacks the phosphodiester bond at the other splice junction, joining the exons and releasing the linear intron. Each step exchanges one backbone connection for another; it is not simply separation of base pairs. See Cech, Zaug and Grabowski, 1981: the phosphoester-transfer mechanism.

The motivating experiment began with an unexpected control. As Cech recounts, Zaug purified a ribosomal RNA precursor from Tetrahymena using a detergent and phenol extraction, then compared test tubes with and without added cell nucleus extract. Incubation in salts and nucleotides released the same amount of small intron RNA in both cases, detected by gel electrophoresis. Repetition made a missing ingredient or accidental omission less plausible. Activity also survived treatments intended to destroy proteins. Nevertheless, absence of added extract did not by itself exclude a tightly attached, resistant protein, or prove that the small RNA had been newly excised rather than released from an already processed complex. See Thomas R. Cech, 1989 Nobel lecture, pp. 655–659: the unexpected control and alternative explanations.

An independent chemical readout addressed the latter problem. Cech mixed precursor RNA with phosphorus-32-labelled guanosine bearing three phosphate groups (GTP), and used gel electrophoresis followed by exposure to X-ray film to locate the label. Label appeared in the excised intron with magnesium present; it was absent when magnesium was omitted or labelled ATP replaced the labelled guanosine reagent. This experiment is described in Cech, 1989 Nobel lecture, pp. 657–658: labelled-guanosine controls. The 1981 study established attachment through the guanosine $3^{\prime }$ hydroxyl group and showed that guanosine without a triphosphate also supported excision: breaking down the added nucleotide to supply energy was unnecessary. See Cech, Zaug and Grabowski, In vitro splicing of the ribosomal RNA precursor of Tetrahymena, Cell 27 (1981), 487–496: primary abstract. Separately, sequencing and terminal analysis found an extra $5^{\prime }$ guanosine absent from the corresponding DNA sequence. See Zaug and Cech, The intervening sequence excised from the ribosomal RNA precursor of Tetrahymena contains a 5-terminal guanosine residue not encoded by the DNA, Nucleic Acids Research 10 (1982), 2823–2838. Together, incorporation and sequence evidence support a new covalent bond, beyond a change in how RNA associates.

Kruger, Grabowski, Zaug, Sands, Gottschling and Cech's 1982 test removed dependence on a precursor purified from Tetrahymena. They cloned the intron and adjacent exons into a DNA template, made precursor RNA with purified bacterial RNA polymerase under conditions suppressing RNA splicing, and removed proteins by detergent–phenol extraction. Denaturing gel electrophoresis separated precursor species. Purified precursor was then incubated with the guanosine reagent, magnesium and salt; matched samples without this incubation supplied a precursor baseline. Their standard assay used $10\mathrm{mM}$ magnesium chloride, $0.1\mathrm{mM}$ of that guanosine reagent, $100\mathrm{mM}$ ammonium sulfate and a buffer solution at pH $7.5$, for $30$ minutes at temperature $37{}^{\circ }\mathrm{C}$, with specified variations in individual experiments. Gel electrophoresis and radioactive-film detection showed excised intron RNA; sequencing verified its $5^{\prime }$ junction. See Kruger et al., Self-Splicing RNA, Cell 31 (1982), 147–157, figures 3–4 and Experimental Procedures.

Excision alone could leave disconnected exons. To test joining, the authors paired labelled RNA with DNA spanning the joined-exon sequence, removed exposed RNA with enzymes, and separated the protected pieces by denaturing gel electrophoresis. Spliced samples yielded protected RNA long enough to span the junction. See Kruger et al., 1982, figure 6: independent evidence for exon ligation. The experimental logic matters: an intron-sized band identifies a candidate product, while sequence and junction assays test its chemical identity.

These controls support sufficiency of the prepared RNA under the assay conditions. They do not show that every RNA self-splices, or that proteins are dispensable inside a biological cell. Many group I introns use proteins to help attain an active fold. See Cech and Herschlag, Group I Ribozymes, section 6: protein facilitation. Nor does a failed assay necessarily identify an inactive nucleotide sequence: unsuitable ions, pH, or temperature can prevent a working active site from forming.

The distinction between self-processing and repeated use became especially clear with bacterial RNase P. In 1983 Guerrier-Takada, Gardiner, Marsh, Pace and Altman separated the RNA and protein components of RNase P from two bacteria. Purified RNA cleaved precursor transfer RNA under conditions containing $60\mathrm{mM}$ magnesium ions, or $10\mathrm{mM}$ magnesium ions plus $1\mathrm{mM}$ spermidine. The protein components alone had no catalytic activity. At $5$–$10\mathrm{mM}$ magnesium without spermidine, both components were needed. Thus the comparison demonstrated an RNA catalyst while also explaining why protein assistance could matter under other conditions. RNase P produces the mature $5^{\prime }$ end of transfer RNA by removing its leader; it acts on a separate substrate and can be used again. See Guerrier-Takada et al., The RNA moiety of ribonuclease P is the catalytic subunit of the enzyme, Cell 35 (1983), 849–857: primary abstract.

Self-splicing demonstrates RNA chemistry, but one precursor completing one rearrangement is not yet evidence for catalytic turnover: repeated conversion of separate substrates by a regenerated catalyst. The self-splicing precursor changes its own connectivity and incorporates external guanosine. In 1986 Zaug and Cech showed that a shortened Tetrahymena intron RNA could cleave and rejoin separate short RNA substrates in a sequence-dependent reaction. This demonstrated an RNA enzyme acting beyond its own self-processing. See Zaug and Cech, The intervening sequence RNA of Tetrahymena is an enzyme, Science 231 (1986), 470–475: primary abstract.

Read the second figure as a counting argument. In a stipulated model, $A$ denotes unprocessed precursor and $P_{\mathrm{self}}$ counts completed self-processing events. Set $A(0)=R_0$, $P_{\mathrm{self}}(0)=0$ and $dA/d\tau =-A$. The exponential function solution gives $P_{\mathrm{self}}/R_0=1-e^{-\tau }$, never exceeding one. Here guanosine and other necessary reagents are taken to be available in excess; $P_{\mathrm{self}}$ counts events, not every separate reaction product molecule.

For a second stipulated model, let an unchanged ribozyme $R$ perform $R+S⟶R+P$ on a separate substrate $S$. Assume the effective law of mass action $dS/dt=-kRS$, no reverse reaction, no ribozyme loss and no inhibited or long-lived bound states. Set $R=R_0$, $S(0)=5R_0$, $P(0)=0$ and $\tau =k{R}_{0}t$. Then $S/R_0=5e^{-\tau }$ and $P/R_0=5(1-e^{-\tau })$, with $S+P=5R_0$. At $\tau =1$, self-processing gives about $0.632$ completed events per initial RNA, whereas the regenerated ribozyme gives about $3.161$ reaction product molecules per initial ribozyme. More than one substrate conversion per initial ribozyme, with appropriate background and identity controls, is evidence for turnover. Product accumulation below that threshold alone cannot distinguish a reusable catalyst from one-shot chemistry. These models illustrate that distinction; actual enzyme kinetics can include folding delays, binding, saturation and product release.

Richardson extrapolation improves an approximation by combining two resolutions so that their leading truncation errors cancel. The useful surprise is that the size of that leading error need not be known: its scaling under refinement is enough. This lets an existing finite difference method or quadrature rule produce a more accurate answer from two runs, without first deriving a new method from scratch.

Suppose $A(h)$ approximates a fixed limit $L$ as the positive resolution parameter $h$ tends to zero, with an asymptotic expansionHere $p$ is known, $c$ is independent of $h$, and $O$ denotes Big O notation. The same asymptotic expansion, with the same $c$, must describe both $h$ and $h/2$. Nested grids make it natural to reuse samples, although nesting itself does not establish this error model. Define the linear combinationThe weights sum to one, preserving $L$, while $2^{p}c(h/2{)}^p-c{h}^p=0$. Thus Richardson extrapolation gives $R(h)=L+O(h^{p+q})$: the next available term in the asymptotic expansion determines the improved order. Justin Solomon's numerical analysis text develops this cancellation for finite differences and applies it to numerical integration.^[ref]

Geometrically, regard the approximation as a function of $t=h^p$. Its first two terms are a straight line $L+ct$. The two computed points lie at positive $t$; extending the line through them to $t=0$ gives $R(h)$. This is why the operation is extrapolation: it estimates a value outside the sampled parameter interval. The remaining curvature corresponds to the uncancelled terms of the asymptotic expansion.

Worked example: a centered derivative of sine. For the smooth function $f(x)=\sin x$, the exact derivative at $x=0$ is $L=\cos 0=1$. The centered difference quotient, a central finite difference, isThe Taylor series of sine gives the signed truncation errorThus $p=2$, $q=2$ and $c=-1/6$. Halving $h$ divides the leading truncation error by four. Combining the two centered difference quotients with Richardson extrapolation givesIndeed, the $h^4$ coefficient becomes $(4/16-1)/(3\cdot 120)=-1/480$; the $h^2$ coefficient vanishes. This is the even-power case of Richardson extrapolation discussed in Persson's authored lecture notes.^[ref]

For $h=0.4$ and $h/2=0.2$, the actual computations areThe absolute error drops from about $2.6454\times 10^{-2}$ for $A(0.4)$ to $5.3080\times 10^{-5}$ for $R(0.4)$, an improvement by roughly a factor of 498. Refining again gives $R(0.2)\approx 0.999996670632607$, whose absolute error is about $3.3294\times 10^{-6}$. The refined truncation errors decrease by nearly four for $A$ and nearly sixteen for $R$, reflecting their second and fourth orders.

The logarithmic scale on each axis turns power laws into straight lines. The dashed reference curves are $h^2/6$ and $h^4/480$, the magnitudes of the leading truncation errors. The plotted errors are computed with 70-digit decimal arithmetic, rather than clipped to a display floor; this isolates the asymptotic expansion from roundoff error.

Worked example: integrating an exponential. Approximate the definite integral of the exponential function on $[0,1]$:For numerical integration, the composite trapezoidal rule adds the areas of trapezoids between successive samples; this is the quadrature use of the endpoint average that also occurs in the ODE trapezoidal rule. On $n$ equal intervals, with $h=1/n$, defineSumming the geometric series gives $T_n=(e-1)\frac{h}{2}\frac{e^h+1}{e^h-1}$. Expanding the exponential function near zero yieldsThe leading truncation error again has $p=2$, $q=2$. On the uniform nested grids with $n$ and $2n$ intervals, Richardson extrapolation therefore givesThis linear combination is exactly composite Simpson's rule on the finer grid, whose spacing is $h/2$. On each coarse panel, an old endpoint receives weight $h/6$ and the new midpoint weight $2h/3$; shared old endpoints receive weight $h/3$. These are the usual Simpson's rule weights $1,4,2,\dots ,4,1$ multiplied by $h/6$. The equivalence is explained directly in LeMesurier's university notes.^[ref]

For two coarse panels and four fine panels,The absolute errors are about $0.0356493$, $0.00894008$ and $0.0000370135$, respectively. Reusing the fine-grid samples makes Richardson extrapolation especially attractive here: the resulting Simpson's rule needs no samples beyond those already used by $T_{2n}$.

When the improvement is justified. Richardson extrapolation requires the stated asymptotic expansion and the correct positive $p$. If the actual leading exponent is $r\ne p$, its coefficient after combination is $c(2^{p-r}-1)/(2^p-1)$, so the leading truncation error generally remains. Merely knowing $A(h)-L=O(h^p)$ does not supply a common coefficient $c$ or a higher-order remainder. A discontinuity of a function, a singular endpoint, or unrelated changes between grids can invalidate the assumed asymptotic expansion; refinement alone does not guarantee improvement for every numerical problem.

Roundoff error and errors in the sampled values survive the cancellation of truncation error. If the coarse and fine values have absolute perturbations bounded by $\epsilon$, the perturbation in $R(h)$ is bounded by $(2^p+1)\epsilon /(2^p-1)$; for $p=2$ this is $5\epsilon /3$. In floating-point arithmetic, a centered difference quotient can also amplify sample errors through division by $h$, so an indefinitely smaller $h$ need not be more accurate. Persson's notes describe this balance between truncation error and roundoff error.^[ref]

The same cancellation idea appears in step-doubling error estimation. For an order-$p$ one-step method for an ordinary differential equation, two half steps each contribute local truncation error of order $h^{p+1}$, so their accumulated leading local error has ratio $2^{-p}$ to one full step. This explains the denominator $2^p-1$ there; the appropriate exponent follows from that local error model, rather than automatically from a generic use of Richardson extrapolation.

The ideal rocket equation gives the change in signed velocity produced by ejecting rocket propellant. Assume one-dimensional Newtonian mechanics in an inertial frame, no external force, and a fixed positive rocket exhaust speed $u$ directed rearward relative to the instantaneous rocket. The rocket's instantaneous mass $M$ includes its remaining rocket propellant; burning it gives $dM<0$. For initial and final masses $M_{\mathrm{initial}}\ge M_{\mathrm{final}}>0$,Here $\ln$ is the natural logarithm, and its argument is the rocket mass ratio. This is the constant-rocket exhaust speed, zero-pressure-correction version of the ideal rocket equation. ^[ref]

For a toy rocket with $u=2km/s$, $M_{\mathrm{initial}}=1000\mathrm{kg}$ and rocket dry mass $M_{\mathrm{final}}=250\mathrm{kg}$ after all rocket propellant is spent, the rocket mass ratio is $4$ andBurning from $1000$ to $500\mathrm{kg}$ gains $2\ln 2km/s\simeq 1.3863km/s$. Burning from $500$ to $250\mathrm{kg}$ gains exactly the same velocity increment, using only $250\mathrm{kg}$ of rocket propellant instead of $500\mathrm{kg}$. Both burns halve the remaining mass. The second burn accelerates a lighter rocket; equal fractional mass losses give equal velocity gains.

The natural logarithm makes the rocket equation both useful and revealing: multiplying successive rocket mass ratios adds their velocity gains. Conversely, the required rocket mass ratio is the exponential function $M_{\mathrm{initial}}/M_{\mathrm{final}}=e^{\Delta v/u}$. Each extra $u\ln 2$ of ideal velocity gain doubles the required initial mass for fixed final mass. This explains why a larger rocket exhaust speed is valuable and why carrying ever more rocket propellant has diminishing returns. Changing the burn duration alone does not change this ideal endpoint result.

For a short proof outline, follow the rocket together with a tiny parcel about to be expelled. Before ejection, their total momentum is $Mv$. After ejection, the remaining rocket has mass $M+dM$ and velocity $v+dv$, while the parcel has positive mass $-dM$ and, to first order, inertial velocity $v-u$. Momentum conservation givesDiscarding the second-order product $dMdv$ gives $Mdv=-udM$; taking the integral gives $\Delta v=-u{\int }_{M_{\mathrm{initial}}}^{M_{\mathrm{final}}}dM/M$. Evaluating the integral yields the natural logarithm above. This is a variable-mass system calculation: the rocket alone is not a closed system, and its momentum changes as the expelled material carries away momentum. The fixed collection consisting of the rocket and all its expelled material has conserved total momentum when external forces vanish.

Fixed rocket exhaust speed means fixed relative speed at release, not fixed inertial velocity. As the rocket's velocity $v$ changes, newly emitted material has changing inertial velocity $v-u$; previously released parcels coast at their own release velocities in this ideal model. If $v>u$, the exhaust even moves forward in the inertial frame, although it is emitted rearward relative to the rocket. No surrounding medium is needed for this exchange of momentum.

Real flight adds gravity, drag force and possibly a pressure correction to rocket thrust. A varying rocket thrust due solely to a changing mass loss rate still gives the same ideal rocket equation if $u$ is constant; it changes the acceleration and burn duration. A changing rocket exhaust speed instead requires $\Delta v=-\int udM/M$, and external forces require their own contribution to the velocity change. The existing constant mass-loss propulsion with quadratic drag and constant-burn rocket with ground friction illustrate such additional forces; the forward-exhaust rocket equation reverses the ejection direction and therefore the sign of the recoil.

For $a>0$, let $g(x)=f(ax)$. Whenever the Mellin transform of $f$ is absolutely convergent,The Mellin fundamental strip is unchanged: the corresponding integrals of absolute values differ by the finite positive factor $a^{-\sigma }$. The change of variables $y=ax$ gives $x^{s-1}dx=a^{-s}{y}^{s-1}dy$, proving the rule. For $s=\sigma +i\tau$, the transform's absolute value changes by $a^{-\sigma }$, while the factor $e^{-i\tau \log a}$ changes its phase. On the line $\sigma =0$, wherever absolute convergence holds, the transform's absolute value is preserved.

For the Mellin transform of an interval indicator, replacing $1_{[1,4]}(x)$ by $1_{[1,4]}(2x)=1_{[1/2,2]}(x)$ givesAt the removable singularity $s=0$, both sides have value $\log 4$. Scaling leaves the width measured by the natural logarithm unchanged. Likewise, the exponential function $e^{-ax}$ has Mellin transform $a^{-s}\Gamma (s)$ on $\mathrm{Re}s>0$.

The sedimentation coefficient $s$ is radial drift velocity divided by local centrifugal acceleration. It has units of time and combines particle properties with those of the surrounding liquid. In a dilute, deterministic mathematical model with negligible particle inertia, an effective displaced volume $V$, particle mass $m$, local liquid density ${\rho }_f$ and linear friction coefficient $f>0$ giveThe subtraction represents buoyancy. Equal mass does not imply equal $s$: changing shape changes drag force, and changing the medium changes both buoyancy and viscosity.^[ref]

For an isolated sphere of radius $b$ and uniform density ${\rho }_p$, Stokes's law supplies $f=6\pi \eta b$, where $\eta$ is dynamic viscosity, and $m=(4\pi /3)b^3{\rho }_p$. ThereforeThis original force-balance check predicts fourfold $s$ on doubling $b$ with everything else fixed, and zero drift when ${\rho }_p={\rho }_f$. It applies to dilute isolated spheres in Stokes flow, not arbitrary folded proteins or crowded suspensions. If $s$ and $\omega$ remain constant in a uniform medium, the ordinary differential equation $dr/dt=s{\omega }^{2}r$ has the exponential function solution $r(t)=r(0)e^{s{\omega }^{2}t}$. Arrival at an outer boundary therefore depends on starting position as well as $s$; a separation has no universal size cutoff.

Take identical rigid spheres of radius $a$, particle density ${\rho }_p$ and displaced liquid density ${\rho }_f$, at uniform absolute temperature $T$. Let $z$ increase upwards from a reference layer and write $\Delta \rho ={\rho }_p-{\rho }_f>0$. The buoyancy correction gives the effective gravitational massThe actual grain mass is ${\rho }_{p}V$; the displaced liquid subtracts ${\rho }_{f}V$ from the effective gravitational mass. In a dilute, noninteracting particle suspension, the Boltzmann distribution therefore givesHere $n$ is number density and $k_B$ is the Boltzmann constant. Equivalently, the dilute osmotic pressure $\Pi =n{k}_{B}T$ supports the grains above it: $d\Pi /dz=-n{m}_{b}g$. Dividing by $n$ and integrating yields the same exponential function. No dynamic viscosity appears in this equilibrium argument. If ${\rho }_p<{\rho }_f$, the population instead increases upwards; at equal densities this gravitational model predicts uniform $n$.

A microscope counts a finite optical slab, not an infinitesimal plane. For equal viewed area $A$ and slab thickness $\delta$, an ideal uniform detection window centred at $z$ givesThus equal optical volumes preserve the exponential function count ratio exactly, even for finite $\delta$, provided the windows lie wholly inside the accessible chamber and have the same detection response. Unequal depth response or overlapping a boundary can corrupt this inference. Repeated observations estimate the mean count; closely spaced visits need not be independent because grains remain in view.

For an original chosen model, set $a=0.250\mu \mathrm{m}$, ${\rho }_p=1200\mathrm{kg}{\mathrm{m}}^{-3}$, ${\rho }_f=1000\mathrm{kg}{\mathrm{m}}^{-3}$, $T=298\mathrm{K}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$. These are illustrative inputs, not Jean Baptiste Perrin's measurements. Using $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ givesIf equal observation effort gives expected count $400$ in the reference layer, the expected counts $20$, $40$ and $60\mu \mathrm{m}$ higher are approximately $214.3$, $114.8$ and $61.49$. Fractions describe expected values over observations, not fractional grains. Using the actual grain mass instead of its buoyancy-corrected value would make the inferred Avogadro constant six times too small for these chosen densities.

A semi-log plot uses a logarithmic scale on one axis and a linear scale on the other. Which axis is logarithmic matters. Plotting $x$ against $\log y$ makes positive exponential functions appear as straight lines; plotting $\log x$ against $y$ instead makes an affine function of a logarithm appear straight. The plotting operations change positions while retaining the original tick values.^[ref]

Take the original dimensionless mathematical model $y(t)=3e^{0.4t}$ for $0\le t\le 8$. It has exponential growth, with $y(0)=3$, $y(4)\approx 14.8591$, and $y(8)\approx 73.5976$. Taking a base-ten logarithm givesThe slope in these transformed positions is $0.4/\ln 10\approx 0.173718$, although the original relative growth rate $y^{\prime }/y$ is $0.4$. Equal steps of two in $t$ multiply $y$ by $e^{0.8}\approx 2.22554$. The corresponding ordinary gains increase, so the left plot bends while the right plot is straight.

One practical use is checking a proposed exponential distribution for failure times. For rate $\lambda >0$ and $t\ge 0$, its cumulative distribution function is $F(t)=1-e^{-\lambda t}$, so the reciprocal survival function $y=1/[1-F(t)]$ satisfies ${\log }_{10}y=\lambda t/\ln 10$. Estimated probabilities can therefore be placed on a logarithmic vertical axis and compared with a straight line. Agreement supports compatibility with the chosen statistical model; a roughly straight finite sample does not establish a unique physical mechanism.^[ref]

In a chain $N_1\to N_2\to N_3$ with distinct decay constants, the Bateman equation expresses each population as a linear combination of the exponential functions $e^{-{\lambda }_{i}t}$. When two consecutive constants coalesce to $\lambda$, the repeated Laplace transform pole produces a term proportional to $t{e}^{-\lambda t}$.

Signaling adaptation is a decline of a response toward its earlier level while the initiating signal remains present. Specify which response is measured and its baseline: complete signaling adaptation means eventual return to that baseline; partial signaling adaptation leaves a lasting displacement. The continuing signal distinguishes signaling adaptation from simply switching off the stimulus. Signaling adaptation describes behavior, while signaling feedback describes a causal arrangement, so neither implies the other.

By reducing a persistent response to background ligand concentration, signaling adaptation can make a biological cell sensitive to changes rather than only to absolute levels. Resetting the response helps leave room for another change over a useful range of ligand concentrations. This is valuable in chemotaxis: swimming bacteria encounter changing chemical conditions, and their signal transduction influences whether a run continues or a tumble occurs. Signaling adaptation helps separate a change encountered during motion from a lasting background. This usefulness does not imply exact return or equal sensitivity at every possible concentration. ^{[ref] [ref]}

Here is a solvable incoherent feedforward loop toy that exhibits complete signaling adaptation without signaling feedback. The input $s(t)$ jumps from $0$ to $1$ at $t=0$ and remains $1$, a Heaviside step function. Let $x$ be a fast activating branch, $a$ a slower inhibitory branch, and $r=x-a$ a readout. The activities $s,x,a,r$ are dimensionless variables; physical time is in seconds. Choose time constants ${\tau }_x=1$ second and ${\tau }_a=5$ seconds, with initial conditions $x(0)=a(0)=0$:Both branches read $s$ directly. Neither branch depends on $r$, and $a$ does not change $x$: the subtraction only defines the readout. Thus the slow opposing influence is feedforward, with no closed feedback control loop.

For $t\ge 0$, solving the two ordinary differential equations givesWriting $t$ for the numerical time in seconds, these are $x=1-e^{-t}$, $a=1-e^{-t/5}$ and $r=e^{-t/5}-e^{-t}$. Initially $r=0$, then $x$ rises faster than $a$, producing a pulse. At long times $x$ and $a$ both approach $1$, so $r$ returns to $0$ despite sustained $s=1$. For this upward step $e^{-t/5}\ge e^{-t}$, hence $r(t)\ge 0$. The striking feature of this signaling adaptation toy is that both branches remain active while the readout becomes quiet: the continuing background and the response to its arrival are separable.

Find the pulse maximum by differentiating the exponential functions. With numerical time in seconds, the derivative isThe stationary condition is $e^{-4t/5}=1/5$, so the natural logarithm givesBefore this time the derivative is positive; afterwards it is negative. Together with $r(0)=0$ and ${\lim }_{t\to \infty }r(t)=0$, this proves that it is the unique pulse maximum. As another numerical check, at $t=10$ seconds $x\simeq 0.999955$, $a\simeq 0.864665$ and $r\simeq 0.135290$: the slow branch is still catching up.

The subtraction $r=x-a$ is a mathematical readout, not a claim that an actual chemical concentration equals this difference for arbitrary inputs. Nonnegativity was established only for the specified upward step; other input histories can give a negative signed readout. Complete signaling adaptation here also relies on equal eventual branch gains: changing the readout to $r=x-0.9a$ leaves $r_{\ast }=0.1$ for $s=1$. This toy illustrates cancellation and different time constants; it does not establish molecular implementation or robustness to arbitrary parameter changes. Ma and colleagues' 2009 study separately analyzed biochemical mathematical models with negative feedback and incoherent feedforward loops, identifying additional conditions for signaling adaptation in those models. ^[ref]

A historical connection to feedback control comes from Yi, Huang, Simon and Doyle's 2000 analysis of a bacterial chemotaxis model. They identified integral action in its regulation of receptor activity: accumulated deviations alter a slower adjustment, which in turn changes the measured activity. At a stationary state, an adjustment that integrates a nonzero error cannot remain constant, explaining how zero error can follow. Their conclusion depended on stated biochemical assumptions and stable behavior; it was not a claim that every biological cell uses that mechanism. A transient pulse alone, as the incoherent feedforward loop toy shows, cannot identify the underlying signaling feedback architecture. ^[ref]

Signaling feedback returns the effect of a downstream event to an earlier stage of signal transduction, closing a loop analogous to feedback control. Trace the causal arrows: if an earlier relay activates a later relay and the later relay changes the earlier relay's activity, the return path closes the loop. Negative feedback opposes the initiating change; positive feedback reinforces it. Signaling feedback is a possible property of signal transduction, not a requirement of every pathway.

For a schematic example, let relay $A$ activate relay $B$. If $B$ inhibits $A$, increased $A$ eventually produces an opposing influence: negative feedback. If $B$ activates $A$, increased $A$ produces a reinforcing influence: positive feedback. These signs describe effects on activity, rather than the direction in which molecules travel. A closed loop containing two inhibitions can have positive feedback, because inhibiting an inhibitor reinforces the original change. The sign belongs to the whole loop.

Negative feedback can restrain excessive activation and help preserve a usable range of responses. Positive feedback can support a decisive transition: Alberts and colleagues describe an enzyme whose product increases its own activity. Such reinforcement can contribute to switchlike behavior, but a positive feedback arrow alone does not prove bistability. A different structure, an incoherent feedforward loop, sends opposing influences forward from the same input; inhibition there is not a downstream return. ^{[ref] [ref]}

A worked mathematical model shows why negative feedback does not automatically produce complete signaling adaptation. Let $s$ be a fixed normalized input, $r$ a normalized response, $g\ge 0$ a dimensionless feedback gain, and $\tau >0$ a time constant. Suppose response removal contributes $-r$ and a return proportional to the response contributes $-gr$. The ordinary differential equation isSetting the derivative to zero gives the equilibrium point $r_{\ast }=s/(1+g)$. A displacement $\delta r=r-r_{\ast }$ obeys $\tau d(\delta r)/dt=-(1+g)\delta r$, so its decaying exponential function establishes asymptotic stability. For a unit input switched on at $t=0$, initial condition $r(0)=0$, $g=2$ and $\tau =1$ second,The negative feedback reduces the eventual response from $1$ to $1/3$, but this response rises monotonically and never returns to its zero baseline. Attenuation and complete signaling adaptation are different properties. This deliberately simple mathematical model specifies no molecular pathway; integral action provides a contrasting way to remove a lasting offset in a stable feedback control system. ^[ref]

A singular perturbation is a parameter-dependent problem whose limiting problem changes qualitatively as the parameter tends to zero. In an ordinary differential equation, a small coefficient of the highest derivative can reduce the order of the limiting differential equation, preventing it from satisfying all the original boundary conditions. Different length or time scales often require different asymptotic expansions.^[ref]

Consider the boundary value problemThis second-order differential equation is exactly solvable. Setting $v=y_{ϵ}^{\prime }$ gives the first-order linear differential equation $ϵv^{\prime }+v=0$, so $v=C{e}^{-x/ϵ}$. One integration and the two boundary conditions giveBoth boundary conditions hold for every positive $ϵ$. For numerical evaluation, write the same ratio as $expm1(-x/ϵ)/expm1(-1/ϵ)$, where $expm1(t)$ evaluates the exponential function minus one without catastrophic cancellation near $t=0$. Keep the denominator for finite $ϵ$: at $ϵ=0.2$ it is approximately $0.993262$, rather than exactly one.

Setting $ϵ=0$ gives $y_0^{\prime }=0$. The outer solution selected by the right boundary condition is therefore $y_0=1$, which fails the left boundary condition. The exact exponential function $e^{-x/ϵ}$ shows how that mismatch is repaired: the solution rises from zero to almost one within a boundary layer of width $x=O(ϵ)$, using Big O notation. For each fixed $x>0$, $y_{ϵ}(x)\to 1$, while $y_{ϵ}(0)=0$ always. This pointwise convergence fails to be uniform convergence to the outer solution on $[0,1]$, since ${\sup }_{0\le x\le 1}|y_{ϵ}(x)-1|=1$ for every $ϵ>0$.

The coefficient $ϵ$ is small, but its product with the second derivative need not be small. Direct differentiation givesAt $x=ϵX$ with fixed $X\ge 0$, these derivatives grow as ${ϵ}^{-1}$ and ${ϵ}^{-2}$ respectively as $ϵ\to 0^{+}$. Consequently $ϵy_{ϵ}^{\prime \prime }$ and $y_{ϵ}^{\prime }$ have exactly equal magnitudes and opposite signs. Even away from the boundary layer, their magnitudes remain equal in this example; both become exponentially small at each fixed positive $x$. The outer expansion neglects an exponentially small transient there, rather than establishing that one derivative term is smaller than the other everywhere.

To predict the width without knowing the exact solution, suppose an order-one change occurs over distance $\delta$. The derivative scales are $y^{\prime }=O({\delta }^{-1})$ and $y^{\prime \prime }=O({\delta }^{-2})$. The dominant balance $ϵ/{\delta }^2\sim 1/\delta$ gives $\delta \sim ϵ$. The resulting inner variable $X=x/ϵ$ makes the boundary layer visible at a fixed scale. This is the useful feature of singular perturbation: an apparently disappearing term identifies the missing structure when one measures its derivatives on the appropriate scale.

Let $\kappa >0$ be an inverse length. In three spatial dimensions, seek the Green's function that decays at infinity and satisfieswhere ${\nabla }^2$ is the Laplacian and ${\delta }^{(3)}$ the three-dimensional Dirac delta function. The point source and the differential operator have spherical symmetry. For $r>0$, the radial Laplacian givesWriting $u(r)=r{G}_{\kappa }(r)$ reduces this to $u^{\prime \prime }-{\kappa }^{2}u=0$. Decay at infinity removes the growing exponential function, leaving $G_{\kappa }(r)=C{e}^{-\kappa r}/r$.

The Dirac delta function fixes the coefficient. Integrating over a ball $B_{ϵ}$ and using the divergence theorem givesThe first term tends to $4\pi C$, while the second is $O({ϵ}^2)$. Hence $C=1/(4\pi )$ andThis solves the source equation in the sense of distributions, including at the origin. Merely verifying the equation for $r>0$ would leave the normalization undetermined. The limit $\kappa \to 0$ at fixed $r>0$ is $1/(4\pi r)$, the Green function of the Laplacian for the convention $-{\nabla }^2$.

For an attractive scalar toy interaction, absorb constants into $A>0$ and setHere $A$ has units of energy times length, and negative $F_r$ points inward. At equal $A$, comparison with the unscreened model $V_0=-A/r$ givesThus the potential energy and force have the same exponential scale but different suppression factors. At $x=1$, these ratios are approximately $0.368$ and $0.736$; at $x=3$, they are $0.0498$ and $0.199$. Finite range means rapid decay, not an abrupt cutoff.

Statistical mechanics connects microstates to macrostates through probability distributions. Its microcanonical ensemble, canonical ensemble and grand canonical ensemble specify different constraints on exchanges with the surroundings. An ensemble average describes repeated preparations under those constraints; it need not be the value in every individual preparation.

This makes the subject useful whenever tracking every microscopic detail would obscure the question being asked. To predict a heat capacity, we need how the expected value of energy changes with temperature, rather than a record of every unit's state. The partition function packages that calculation into one reusable object.

In his 1902 book, J. Willard Gibbs systematically developed the canonical ensemble. His preface contrasts following one system through time with studying a collection of systems with different microscopic conditions. Chapter IV defines the canonical ensemble through an exponential function of energy and establishes its normalization. This supplied a general strategy for predicting collective measurements by averaging over possible microscopic conditions, with the partition function organizing the calculation. ^{[ref] [ref]}

If $n\ge 1$, $f$ is real and $C^n$ on an interval, and $f^{(n)}$ is nowhere zero there, the displayed inequality holds. If the errors were equal, a best degree-at-most-$n-1$ polynomial would also be best at degree $n$. The Chebyshev alternation theorem and intermediate value theorem would give $n+1$ zeros of the error. Applying Rolle's theorem $n$ times would force a zero of $f^{(n)}$, a contradiction. The zero-error case already contradicts $f^{(n)}\ne 0$. In particular this proves strict improvement for the exponential function at every positive degree step.

On $S_h=\{z:0<\mathrm{Im}z<h\}$, let $\beta \in S_h$ and write $e^{\pi \beta /h}=A+iB$ with $B>0$. The exponential function $E(z)=e^{\pi z/h}$ is a conformal map from the strip bijectively onto the upper half-plane, sending $ih/2$ to $i$. The affine map $w\mapsto A+Bw$ is an automorphism of the half-plane sending $i$ to $E(\beta )$. Conjugate it by $E$, using the principal complex logarithm as the inverse. The displayed formula is a conformal automorphism of the strip and satisfies $g_{\beta }(ih/2)=\beta$. The branch range $0<\mathrm{Im}\mathrm{Log}w<\pi$ is what prevents an unwanted vertical period in this inverse.

For a function with derivatives of every order at $a$, its Taylor series is the power seriesIts partial sums are Taylor polynomials. At $a=0$ it is called a Maclaurin series. Forming the coefficients, proving that its partial sums form a convergent sequence, and proving that their limit equals $f(x)$ are separate tasks. A smooth function provides all the coefficients, but a real analytic function is additionally represented locally by the resulting power series.

For the exponential function, every derivative of $e^x$ at zero equals one, givingHere equality holds for every real $x$: the Taylor remainder is bounded by $e^{|x|}|x{|}^{N+1}/(N+1)!$, which tends to zero for each fixed $x$. At $x=1$, the degree-five Taylor polynomial is $163/60\approx 2.71666667$, and the absolute difference from $e$ is less than $e/720$. The factorial denominator eventually overwhelms every fixed power, allowing local derivative information to reconstruct this function throughout the real line.

The geometric series shows why the region of validity matters. For $f(x)=1/(1-x)$, the Taylor series at zero is ${\sum }_{k\ge 0}{x}^k$, with radius of convergence one. For $x\ne 1$, its finite Taylor polynomial satisfies the exact identityAt $x=1/2$ and $N=5$, the remainder is $1/32$. At $x=0.9$, the same degree-five polynomial has value $4.68559$ and remainder $5.31441$, although the Taylor series still converges to $10$ there. Convergence of the infinite expression does not guarantee accuracy from a particular small number of terms. For $|x|>1$, even the terms fail to tend to zero, so this power series cannot represent the function. The standard expansions and the conditions for convergent series are discussed in OpenStax, working with Taylor series.

A stronger caution comes from the flat functionIt is a smooth function with $h^{(n)}(0)=0$ for every $n$, so every Taylor polynomial at zero is zero and its Taylor series converges everywhere to zero. Nevertheless $h(1/2)=e^{-4}>0$. Off zero, repeated differentiation produces a polynomial in $1/x$ times $e^{-1/x^2}$; the exponential decay makes each required limit at zero vanish. This explains how all the local coefficients can miss a positive function arbitrarily close to the centre. Real analyticity supplies the additional condition that rules out this behavior.

Release of electrons from a heated material; one source of electrons in an X-ray tube.

Heating a metal gives some electrons enough energy to escape its surface. A heated diode cathode can therefore supply an electric current across a vacuum gap to a collecting diode anode. This makes temperature a control for the supply of electrons. Conversely, measuring that supply at several temperatures lets us estimate the surface work function: the energy barrier measured from the Fermi level to the nearby vacuum level. A Richardson plot turns the strong exponential function dependence into a slope.

Three questions must be kept separate: how many electrons cross the surface, how many traverse the gap, and how many the collecting diode anode records. The Richardson-Dushman law answers the first under specified surface assumptions. Space charge can obstruct the second, while incomplete collection can alter the third.

For any fixed real $n$-by-$n$ matrix $A$, as the real parameter $t$ tends to zero,The big O notation means the remainder is bounded in absolute value by $C{t}^2$ for some constant $C$ near zero. The determinant stays positive for sufficiently small $t$, so its oriented volume factor is also its ordinary volume factor. Thus the trace is the first-order fractional volume change, measured from the identity matrix.

A proof follows directly from multilinearity of the determinant. Each column of $I+tA$ is $e_j+tA{e}_j$. The constant term chooses all the standard basis columns and equals one. A term containing exactly one factor $t$ replaces just column $j$ by $A{e}_j$; expanding that vector in the standard basis leaves only its $e_j$ component, $a_{jj}$, because every other choice repeats a column. Summing gives $\mathrm{tr}A$. Terms with two or more replacements have order $t^2$ or higher. This also proves the derivative of the determinant at the identity matrix.

This coefficient has a direct visual meaning. At the point $e_j$, the linear map $x\mapsto Ax$ has component $a_{jj}$ along direction $e_j$. In a small step $x\mapsto x+tAx$, these are the first-order changes of the square's or box's sides in their own directions. Summing them gives its fractional area or volume change; the off-diagonal entries initially change shape rather than contributing to this first-order total.

For the constant-coefficient linear system of differential equations $x^{\prime }=Ax$, the exact linear flow is $x(t)=e^{tA}x(0)$, using the matrix exponential. Its volume factor isThis is the matrix exponential determinant identity. One can derive it without diagonalization: set $E(t)=e^{tA}$, whose matrix inverse is $E(-t)$, and use $E^{\prime }=AE$ in Jacobi's formula to get $V^{\prime }=\det E\mathrm{tr}(E^{-1}AE)=(\mathrm{tr}A)V$. The scalar differential equation with initial value one gives the displayed exponential function. Thus positive trace causes volume expansion, negative trace causes contraction, and zero trace preserves volume, even when individual directions behave differently.

In two dimensions, compare four generators at $t=0.4$:The first and third diagonal matrices expand or contract both directions. The second expands horizontally and contracts vertically by reciprocal factors, making a rectangle of unchanged area. For the last nilpotent matrix, $A^2=0$, so its matrix exponential is exactly $I+tA$; it sends $(x,y)$ to $(x+ty,y)$, a shear mapping. Zero trace allows substantial shape change and does not mean the flow map is stationary.

The first-order statement and the exact flow map should be distinguished. For $A=\mathrm{diag}(1,-1)$, $\det (I+tA)=1-t^2$, whereas $\det (e^{tA})=1$ for all real $t$. Zero trace eliminates the first-order change in the approximate step, and it preserves volume exactly in the matrix exponential flow. The examples use an abstract parameter $t$ and coordinate vectors; no physical units are assumed.

For motion confined to a plane but retaining potential energy $-\kappa /r$, separation in polar coordinates gives angular wavefunctions $e^{im\varphi }$ with integer $m$. The regular radial factor behaves as $r^{|m|}$ at the origin. Removing this power and the decaying exponential function leaves an Generalized Laguerre polynomial differential equation. polynomial termination requires $N=n_r+|m|+1/2$, yielding the displayed energy. At fixed $q=n_r+|m|$, the spinless quantum degeneracy is $2q+1$. This is not the logarithmic potential of electrostatics formulated entirely in two spatial dimensions.

Two-state ion-channel gating models an individual gate as a random switch between closed and open states. Write $p(t)$ for its open probability; this is the $X$ in the header equation. An individual gate's state is always closed or open, rather than partly open. A population of identical gates can have a fraction open near $p(t)$, and averaging repeated identically prepared trials estimates the same probability. A finite population's observed fraction fluctuates. See Purves and colleagues, Neuroscience, 2nd edition (2001): Ion Channels Underlying Action Potentials for single-channel steps and their population averages.

At fixed membrane potential $V$, let constant rates $\alpha ,\beta >0$ govern closed-to-open and open-to-closed transitions respectively, in ${\mathrm{s}}^{-1}$. This is a two-state continuous-time Markov chain. Over a short interval $\Delta t$, an initially closed gate opens with probability $\alpha \Delta t+o(\Delta t)$, and an initially open gate closes with probability $\beta \Delta t+o(\Delta t)$. Gain minus loss gives the deterministic ordinary differential equation for the probability:Here $p_{\infty }$ is the fixed-voltage limiting open probability, and $\tau$ is a time constant. The ratio of the rates fixes the limit, while their sum fixes the speed. With initial probability $p(0)=p_0\in [0,1]$,For a short derivation, subtract $p_{\infty }$: the deviation obeys $d(p-p_{\infty })/dt=-(p-p_{\infty })/\tau$, yielding the exponential function. After one time constant, $1/e$ of the initial deviation remains. The solution stays between $p_0$ and $p_{\infty }$, so it remains a valid probability. The rate and time constant formulations are described in Gerstner and colleagues, Neuronal Dynamics, §2.2.2: transition rates and channel kinetics.

For a worked toy condition A, choose $\alpha =\beta =50{\mathrm{s}}^{-1}$ and $p_0=0.1$. Then $p_{\infty }=0.5$, $\tau =0.01\mathrm{s}=10\mathrm{ms}$, andIn a separately prepared condition B, choose $\alpha =200{\mathrm{s}}^{-1}$, $\beta =50{\mathrm{s}}^{-1}$ and the same $p_0=0.1$. Now $p_{\infty }=0.8$, $\tau =4\mathrm{ms}$, and $p(10\mathrm{ms})=0.8-0.7e^{-2.5}\approx 0.742540$. Condition B approaches a higher limit more quickly. These are two stipulated constant-rate conditions with $V$ held fixed in each experiment, not an empirical function relating rates to voltage or a sequential switch between conditions. The initial $p_0=0.1$ need not be the limiting probability at the newly imposed voltage.

In the Hodgkin-Huxley model, the rates depend on membrane potential: ${\alpha }_X(V)$ and ${\beta }_X(V)$. Its $X_{\infty }(V)={\alpha }_X(V)/({\alpha }_X(V)+{\beta }_X(V))$ is the limit for voltage held at that $V$. During an action potential, changing $V$ changes the rates and gates feed back onto $V$ through ionic electric currents; substituting the present voltage into the constant-rate solution does not solve that coupled system. Gating describes availability of a conducting pathway. Even an open ion channel carries zero net ionic electric current at its reversal potential.

The useful bridge is from abrupt random molecular events to a smooth population prediction. A channel requiring several independent permissive gates in an ion channel has an open probability equal to their product; this can make channel availability much smaller than each gate's availability. The Hodgkin-Huxley model uses ${\bar{g}}_{\mathrm{Na}}{m}^{3}h$ and ${\bar{g}}_{\mathrm{K}}{n}^4$ for sodium and potassium electrical conductances, respectively.