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.

A delayed feedback characteristic root solves the characteristic equation of a delay differential equation and supplies an exponential function mode. Its real part determines decay or growth of that mode.

A feedback delay is a time interval between an event and its returning influence. Multiple biochemical steps can introduce an effective lag; a fixed delay is an idealization, not a universal description of signal transduction. To isolate timing, let $x$ be a signed displacement from an equilibrium point and use normalized dimensionless variables:The local removal is $-x$ and the returning influence is negative at every delay. At $d=0$ the ordinary differential equation reduces to the decay already illustrated under signaling feedback. For $d>0$ a delay differential equation requires a history on $[-d,0]$, not just an initial condition at one instant. Negative values of $x$ are displacements, not negative chemical concentrations.

An exponential function $x=e^{\lambda t}$ gives the characteristic equation of a delay differential equationPutting $\lambda =i\omega$ and equating real and imaginary parts gives $\cos (\omega d)=-1/2$ and $\omega =2\sin (\omega d)$. The first crossing with $\omega >0$ therefore occurs atFor this particular mathematical model, the equilibrium point is asymptotically stable for $0\le d<d_c$. No nonnegative real delayed feedback characteristic root can satisfy the equation; the first loss of asymptotic stability is the indicated pair of imaginary delayed feedback characteristic roots. Differentiating the characteristic equation of a delay differential equation along a root yieldsThus the pair moves into growth as the delay increases through $d_c$. The neutral mode at $d_c$ is a periodic function, but its amplitude depends on history: this linear mathematical model does not select an attracting limit cycle.

The figure plots $x(t)=e^{at}\cos (bt)$ using a delayed feedback characteristic root $\lambda =a+ib$ and that same formula on the history interval. At $d=0.7$ the displayed mode decays; at $d=1.5$ it grows. These are exact modes up to numerical root accuracy, not step responses with identical histories. Molecular bounds and a changed return law would be needed to interpret what happens after an unstable biological response becomes large.

Put $u=1+x$. Then $\dot{x}=-x(t-\tau )-x(t)x(t-\tau )$, so the linearization is $\dot{x}=-x(t-\tau )$. A mode $e^{\lambda t}$ gives the characteristic equation of a delay differential equationFor real $\lambda \ge 0$, the left side is nonnegative and the right side negative. The two real graphs therefore never meet there, as drawn in the figure. There is no nonoscillatory exponential instability. If stability is lost, the unstable eigenvalues must be a nonreal conjugate pair, producing oscillatory growth; absence of a positive real root alone does not imply stability.

In this delayed carbon-dioxide ventilation feedback model, $p$ is a representative carbon-dioxide partial pressure, $K$ an effective storage capacitance, $M>0$ metabolic production, and $pV$ removal by pulmonary ventilation. The delay represents transport and sensing between the gas-exchange compartment and the ventilatory response. The Hill equation response rises sigmoidally with pressure, reaches half its maximum $\alpha$ at $p=\theta$, and saturates at $\alpha$; $n>1$ permits a steep threshold-like response.

For positive parameters, $pV(p)$ increases strictly from zero to infinity because its derivative is $V+p{V}^{\prime }>0$. Its graph therefore intersects the horizontal line $M$ exactly once: $\boxed{M=p^{\ast }V(p^{\ast }),p^{\ast }>0}$. The figure shows this argument and the first imaginary-root crossing.

For a small perturbation $z=p-p^{\ast }$, linearization givesThe sigmoidal delayed ventilation threshold follows from the characteristic equation of a delay differential equation, $\lambda +A+B{e}^{-\lambda \tau }=0$. At an imaginary root $\lambda =i\omega$, separate real and imaginary parts:For $B>A>0$, the first positive crossing hasAt zero delay the scalar equation is stable. No root can cross the imaginary axis before this first crossing, and zero itself is not a root. To check the crossing direction, implicit differentiation givesAll subsequent imaginary crossings have the same direction. Thus the equilibrium is linearly stable below ${\tau }_c$ and unstable above it. This establishes the oscillatory linear threshold; the full nonlinear Hopf bifurcation additionally requires the usual nondegeneracy conditions. If $B\le A$, no finite imaginary-root crossing occurs and the equilibrium is stable for every finite delay. For the stated response, $B/A=n{\theta }^n/({\theta }^n+(p^{\ast }{)}^n)$, showing explicitly when the gain can exceed instantaneous damping.

If $A\ll B$, then ${\omega }_c\approx B$ and ${\omega }_c{\tau }_c\approx \pi /2$. Consequently the small-amplitude period near onset is $\boxed{2\pi /{\omega }_c\approx 4{\tau }_c}$. The comparison is to the critical delay near the first instability, rather than an arbitrary delay far beyond onset. A longer circulatory transit time can move the physiological system across this delayed-feedback threshold, providing a mechanism for periodic breathing associated with impaired cardiac output.

The delayed positive term models recruitment depending on the population a time $T$ earlier; the quadratic negative term represents instantaneous density-dependent losses. Thus this delay differential equation delays recruitment rather than the strength of crowding acting on the present population.

With $x=1+y$, the linearized equation is $\dot{y}=\alpha y(t-T)-2\alpha y(t)$. Its characteristic equation of a delay differential equation is $\lambda +2\alpha =\alpha e^{-\lambda T}$. A root with nonnegative real part would satisfy $|\lambda +2\alpha |\ge 2\alpha$, whereas $|\alpha e^{-\lambda T}|\le \alpha$. This is impossible. Consequently the equilibrium is locally exponentially asymptotically stable for every $\alpha >0$ and $T>0$. There is no upper bound on $\alpha T$ in this model.

Put $x=1+y$ and retain linear terms. The delay differential equation is $\dot{y}(t)=-\alpha y(t-T)$, with characteristic equation of a delay differential equation $\lambda +\alpha e^{-\lambda T}=0$. Write $z=\lambda T$ and $q=\alpha T$. If a root has $\mathrm{Re}z\ge 0$, then $|z|=q{e}^{-\mathrm{Re}z}\le q$. Writing $z=a+ib$ gives $a=-q{e}^{-a}\cos b$. For $q<\pi /2$, $|b|<\pi /2$ and the right side is negative, contradicting $a\ge 0$. All characteristic roots therefore lie in the left half-plane.

An imaginary root $z=ib$ requires $\cos b=0$ and $b=q\sin b$. The first crossing is $q=\pi /2$, $z=\pm i\pi /2$. Implicit differentiation gives $dz/dq=-e^{-z}/(1+z)$, whose real part at the positive crossing is $(\pi /2)/(1+{\pi }^2/4)>0$. The later imaginary crossings occur at $q=\pi /2+2\pi k$ and also enter the right half-plane. Thus the equilibrium is locally exponentially asymptotically stable whenand unstable when $\alpha T>\pi /2$. At equality the linearized equation has neutral oscillatory modes; it is the Hopf bifurcation threshold. If stability is meant for the nonlinear equation itself, the borderline can also be resolved. In dimensionless time $s=t/T$, let $q_0=\pi /2$. A small center oscillation has the form $y=R\cos \theta +R^2(\cos 2\theta /5+\sin 2\theta /10)+O(R^3)$. Substitution at quadratic order gives these second-harmonic coefficients. At cubic order the fundamental forcing is $q_0{R}^3(\cos \theta /20+3\sin \theta /20)$. Projection onto the center mode divides the complex forcing by the characteristic derivative $1+i{q}_0$, yielding the Hopf normal form amplitude equationIts cubic coefficient is negative, so the equilibrium at equality is locally asymptotically stable, with algebraic rather than exponential decay. Thus the usual strictly decaying linear condition is $\alpha T<\pi /2$, while nonlinear local asymptotic stability includes $\alpha T=\pi /2$. The purely imaginary linear modes alone would not justify that borderline conclusion.

Put $x(t)=1+u(t)$. Retaining terms linear in $u$ givesThe exponential ansatz $u=e^{\lambda t}$ gives the characteristic equation of a delay differential equationAt $T=0$ its root is $\lambda =-\alpha$, and the remaining delay roots lie far in the left half-plane, so the equilibrium is stable for sufficiently small $T$.

A stability change can first occur when $\lambda =i\omega$. Separating real and imaginary parts givesFor the first positive crossing, $\omega T=\pi /3$ and $\omega =\sqrt{3}\alpha$. ThusAll positive-frequency crossings have $\omega T=\pi /3+2\pi k$, so there is no earlier imaginary-axis crossing. Implicit differentiation ofgives, with $q=\alpha T_c$,so the conjugate pair crosses into the right half-plane as $T$ increases. Therefore the stability threshold for delayed quadratic crowding is

For $K\dot{p}=M-pV(p(t-\tau ))$ and a positive steady state, linearization gives $\dot{z}=-Az-Bz(t-\tau )$, where $A=V(p^{\ast })/K$ and $B=p^{\ast }{V}^{\prime }(p^{\ast })/K$. Splitting the characteristic equation of a delay differential equation at an imaginary root yields ${\omega }^2=B^2-A^2$ and the displayed first critical delay. The root crosses into the unstable half-plane with real derivative ${\omega }^2/[(1+A{\tau }_c{)}^2+{\omega }^2{\tau }_c^2]>0$. If $A\ll B$, the period near onset is approximately $4{\tau }_c$. For a Hill equation response, $B/A=n{\theta }^n/({\theta }^n+(p^{\ast }{)}^n)$ identifies whether delayed gain can exceed instantaneous damping.

A longer circulation time increases the delay between gas exchange and its sensed feedback. The American Thoracic Society's documented case relates a change in breathing pattern during one night to a longer estimated circulation time; the simple scalar model isolates a possible feedback mechanism rather than representing every determinant of that observation.^[ref]