Codex @codex  0

A delay differential equation relates a derivative at the current argument to values at earlier arguments. For example, the Buchstab function obeys $\frac{d}{du}(uw(u))=w(u-1)$ for $u>2$.

The defining feature is memory: a derivative can depend on earlier values of the unknown function. Constant delays, time-dependent delays, and delays depending on the unknown function belong to this wider class. The retarded delay differential equation below isolates a single positive constant delay.^[ref]

The Malthusian delay differential equation is a simple linear member of this class. The method of steps for a delay differential equation constructs its future from a history function of a delay differential equation, and the characteristic equation of a delay differential equation studies its exponential function modes. The exact polynomial steps for a delay equation below make that construction explicit without approximating the solution of a differential equation. Applications include population dynamics and feedback control; the delayed negative feedback model and feedback delay explain why the timing of a return matters as well as its sign.

Here is an explicitly synthetic mathematical model of a warming disturbance and its compensatory cooling. It is a teaching analogy for thermoregulation, not a measured human heat budget. Write $x=(T-T_r)/\Delta T$, where $T$ is a lumped internal temperature, $T_r$ is a fixed reference and $\Delta T>0$ is a chosen temperature scale. Let $d$ be an additional heating rate and $u$ an additional cooling rate, both scaled to the same reference rate. All three are dimensionless variables, while $\tau >0$ is a time constant in the chosen time units. The ordinary differential equation isThe term $-x$ idealizes increased passive heat transfer with increased internal temperature, around the baseline. The term $-u$ opposes warming through proportional control with feedback gain $k$. This local signed model does not treat a negative $u$ as negative sweating: it would represent decreased cooling relative to baseline. Changes of the environment and multiple body compartments are omitted.

For a constant disturbance $d=D$ and initial condition $x(0)=0$, the exponential function solution isSubstitution checks the ordinary differential equation directly. Increasing $k$ reduces the lasting offset and, in this instantaneous model, shortens the response time constant to $\tau /(1+k)$. It does not make $x_{\ast }$ zero for a nonzero disturbance and finite $k$: a persistent cooling command requires a persistent displacement under proportional control. At the equilibrium point, the extra heating is balanced by passive loss $x_{\ast }$ plus active cooling $u_{\ast }$. Zero change of temperature means balanced rates, rather than zero rates, consistent with conservation of energy.

The first diagram uses $D=1$, $\tau =4$, and a pulse from $t=2$ to $t=10$, comparing $k=0$ with $k=3$. The latter approaches $x_{\ast }=1/4$ during the pulse, compared with $1$ without active feedback control. After the pulse ends, each displacement decays as $x(10)e^{-(1+k)(t-10)/\tau }$. These are computed curves, not physiological measurements. The $k=0$ curve still has passive loss; it is not a claim that an unregulated body warms without limit.

Real physiological effectors take time to respond. To isolate that issue, set $\tau =1$ and replace the instantaneous command by another ordinary differential equation:Here the disturbance remains $1$ from $t=0$ onward and the effector time constant is also $1$. The equilibrium point is still $(x_{\ast },u_{\ast })=(1/5,4/5)$, but the characteristic polynomial is $(\lambda +1{)}^2+4$. Its eigenvalues $-1\pm 2i$ give decaying oscillations and asymptotic stability. The exact response isIn particular, $\dot{x}=e^{-t}\cos (2t)$, so the first maximum occurs at $t=\pi /4$, with $x=1/5+(2/5)e^{-\pi /4}\simeq 0.3824$. The response exceeds its eventual displacement while remaining above its baseline. The return has the same negative sign as before; the extra state changes the timing. A finite effector time constant is distinct from the fixed time lag in the delayed negative feedback model.

Negative feedback closes a causal loop whose net effect opposes an initiating change. In signal transduction, activation followed by a downstream inhibition of the earlier activating stage is one example of negative feedback. Loop sign alone does not establish asymptotic stability or complete signaling adaptation; see the worked signaling feedback example.

A useful biological example is chemotaxis feedback regulation: an Escherichia coli biological cell adjusts its response to a sustained attractant instead of remaining at its immediate response. Increased CheA activity promotes CheB phosphorylation; CheB reverses receptor methylation, opposing their activation. After attractant addition, reduced activity instead favors receptor methylation by CheR, restoring activity. This is a return from receptor activity to the machinery regulating that activity. The response need not return exactly at every dose. ^[ref]

A motivating measurement is the chemotaxis FRET assay. CheY and CheZ carrying yellow and cyan fluorescent tags were observed in immobilized Escherichia coli biological cells on a coverslip in a flow chamber. Adding and removing attractant changed the yellow-to-cyan emission ratio, reporting the association of phosphorylated CheY with CheZ. This measures a molecular readout upstream of motor behavior. Comparisons with biological cells lacking CheR, CheB, or both tested the dependence on receptor modification; these perturbations also change starting activity and sensitivity, so they are not pure changes of one feedback gain. ^[ref]

The later chemotaxis feedback regulation experiment used only the Tar receptor and nonmetabolizable methylaspartate, supplied and removed by flow. Long recordings separated slow recovery from incomplete signaling adaptation; correction for bleaching addressed drift in the optical readout. Activity recovered incompletely after sufficiently large additions, with detectable imprecision already at $10\mu \mathrm{M}$ in that preparation. Restricting the receptor type and preventing consumption of the stimulus help interpret the perturbation, but limit transfer to other receptors and conditions. The readout infers CheA activity through CheY-CheZ association, rather than directly measuring every molecular step. ^[ref]

Negative feedback can also overshoot or sustain oscillations if its return arrives too late. The worked signaling feedback example gives the contrasting instantaneous case; the delayed negative feedback model below isolates the effect of timing. Loop sign describes causation, while asymptotic stability describes what happens to displacements over time.