The second box represents a reservoir exchanging carbon with the atmosphere, with equal capacity in the chosen concentration units. For a unit atmospheric pulse, , and no subsequent emissions, total is conserved while . Thus
The atmospheric burden falls from one to one half with time constant , not . The model has a permanent airborne fraction because neither box contains an external sink. The carbon-response panel of the original figure in S12 shows this explicitly.
Patch-clamp series resistance, also called patch-clamp access resistance, is the effective electrical resistance between the pipette's sensed voltage and the cell interior. It includes tip access and any remaining obstruction after break-in. For injected electric current into the biological cell, an ideal resistive access path gives
Thus through produces a difference: a substantial possible voltage clamp error. The electric current in this equation is the access-path current, excluding any bypass seal leak current. In current clamp, the same access drop contaminates the pipette voltage during injection. Bridge compensation subtracts an estimated ; an estimate that changes during recording leaves an error.
An independent passive voltage clamp calculation makes both speed and error visible. Write and , with fixed leak reversal potential , constant membrane electrical resistance and membrane capacitance . Assume an ideal command at the pipette node and omit seal leak current and pipette capacitance. Electric charge conservation gives
For a step from , the ordinary differential equation gives
With the figure's parameters and , the time constant is , the final membrane rise is , and access-path electric current falls from to . Only when does the time constant approach and the final voltage approach the command. Actual voltage clamp also depends on amplifier feedback control and the space clamp problem.
Choose the following original phenomenological mathematical model solely to calculate a paired-pulse readout. It is not experimental data, a fit to the teaching experiment, or a simulation of the Hodgkin-Huxley model. The conditioning action potential defines time zero. Test pulses have fixed width ; denotes their inward electric current amplitude. Impose an absolute neuronal refractory period . No finite test amplitude generates an event for . After that hold, stipulate a required amplitude
with , dimensionless and recovery time constant . Emit a test event exactly when , including equality. For this ideal event readout, assign events to pulse onset and assume the conditioning event is successfully evoked. There is no calculated membrane potential, spike waveform, response latency or stochastic success probability.
Figure 1.
Chosen refractory hold and required test-pulse amplitude
. Original phenomenological paired-pulse recovery rule, not measured data or a channel simulation. The grey region forbids a second event; after , the required electric current falls from toward . Horizontal lines show two chosen test amplitudes; their intersections give different first permitted test delays. The dashed baseline is approached asymptotically. Details.
The exponential decay represents a combined recovery deficit; it is neither the Hodgkin-Huxley model's nor a measured potassium electrical conductance. One such curve cannot separately identify sodium recovery and persistent potassium electrical conductance. A constant pulse-width amplitude criterion also does not establish an equivalent rule for other pulse widths or voltage trajectories.
At the hold boundary, . At ,
at ,
Thus a test fails at because of the imposed hold, fails at because recovery is insufficient for that amplitude, and succeeds at . At , a pulse would instead succeed: failure of the weaker pulse was relative refractoriness.
Figure 2.
Paired pulses and event readout at three chosen delays
. Each row is a separately reset ideal trial with a conditioning pulse assumed to evoke an action potential at time zero and a test pulse. Rectangles are imposed electric current pulses; vertical ticks are abstract events, not voltage recordings. Only the test passes the stipulated recovery rule. The and failures have different causes. Details.
For , solve using the logarithm. The earliest allowed delay is
For this is ; for it is . If , the imposed hold is the only delay restriction. If , no finite delay passes the rule: the required amplitude stays strictly above its asymptotic baseline. A finite experimental “fully recovered” time therefore also needs a tolerance and a readout definition.
RC charging circuit 2026-10-24
An ideal constant-voltage source charging an initially uncharged capacitor through a series resistor gives and . These follow from Kirchhoff's voltage law and . A smaller time constant means faster fractional settling; it does not fix the initial current or final charge unless the separate values of and are specified.
RC low-pass filter 2026-10-18
An RC low-pass filter places a resistor between input and output, with a capacitor from output to their common return. The ideal model prescribes the input electric potential difference and leaves the output unloaded. It is a first-order low-pass filter with unit eventual gain.[ref]
Figure 1. . Original schematic: the output is across the capacitor, not across the resistor. The orange dashed resistor is the optional electrical load used in RC low-pass output loading; omit it for the ideal unloaded model. The blue arrow defines positive electric current through . Details.
For electrical resistance and capacitance , Ohm's law gives the electric current through the resistor as . The capacitor relation gives . Equating them by Kirchhoff's current law yields the ordinary differential equation
The time constant has units of time. A rapid change in electric potential difference requires a large electric current into the capacitor; the finite electrical resistance limits that rate. A slowly maintained voltage instead lets the capacitor approach the input voltage. The same equation therefore predicts both smoothing of sustained oscillations and a rounded response to a sudden change.
For the original constructed example throughout this article, take
These stipulated ideal values specify an analytical model, rather than measurements of particular components.
RC low-pass laboratory 2026-10-18
The University of Colorado Boulder's Lab 3: Filters, §§7.1–7.4, specifies measuring the resistor and capacitor, then driving the assembled RC low-pass filter with a function generator. A two-channel oscilloscope reads actual input and output using high electrical impedance. Sweep the frequency, record the voltage amplitude ratio, and locate . A control compares a cable with a oscilloscope probe on a resistive voltage divider, isolating measurement capacitance before using the probe for the filter. A slow square wave then lets discharge timing estimate the time constant independently. Agreement supports the RC low-pass frequency response; source output impedance, probe loading, and instrument signal bandwidth limit it. These are prescribed measurements, with no measured dataset supplied here.[ref]
The two-channel setup can also test the predicted phase difference directly: for positive frequency , measure the output's time lag between corresponding features of the two sinusoids. For this single-stage model the phase lies between and , so choose the lag within a quarter cycle and calculate . At the cutoff frequency the predicted lag is one eighth of a cycle. This tests the timing part of the frequency response as well as the voltage ratio.
A real electrical load can draw electric current from the output. For a load resistor across the capacitor, Kirchhoff's current law instead gives
Therefore the loaded continuous-time frequency response is
Loading reduces the eventual gain and shortens the time constant. With the constructed values and , , , and the corner relative to the loaded low-frequency gain moves to . A step now approaches . A quicker rise alone would conceal the error in the final voltage.
With the ideal input voltage source set to zero, the output impedance seen by the electrical load is
This follows from Thévenin's theorem: the resistor and capacitor are then in parallel. A sufficient small-loading condition is over the frequencies of interest. An oscilloscope is itself an electrical load: its input electrical resistance and cable or probe capacitance can both change the response. Pure extra parallel capacitance increases and lowers the corner, whereas a finite parallel electrical resistance gives the different gain-and-corner changes above.
The definition of input also matters. If voltage is prescribed inside a source with series electrical resistance , the unloaded response to that internal source uses . If the actual voltage at the filter input terminal is measured, the ratio of output to that measured input still uses . Measuring both terminals avoids mistaking a source setting for the voltage delivered to the filter.
RC low-pass step response 2026-10-18
The step response describes a sudden persistent input change. Suppose the input becomes a constant at , and the initial capacitor voltage is . Finite electric current leaves that voltage continuous. Subtract from the RC low-pass filter equation: the remaining error satisfies . Separating variables gives
For an initially uncharged capacitor and a to step in the constructed example,
At one time constant, , the voltage is , about of its final value. The voltage error then continues its exponential decay. This is different from a pure delay: the output starts changing immediately, yet takes time to approach the target.
Figure 1. . Original analytical curves for and ; no experimental points. Left: a Bode plot of voltage gain and output-minus-input phase difference, with cutoff frequency . Right: the response to a to step, compared with doubling at fixed . Doubling the time constant halves the cutoff frequency and doubles the time to reach the same fraction of the final voltage. Details.
Define settling time relative to the total voltage change: require thereafter, with . The exponential decay is monotone, so
For settling, ; the to rise takes . Reducing the cutoff frequency by doubling doubles both times. This reciprocal relation is the useful connection between signal bandwidth and response speed: stronger suppression of rapid noise comes with slower tracking of wanted changes.
Reference signal 2026-10-14
A reference signal specifies the desired physical plant output as a function of time. A fixed reference signal is often called a setpoint; a changing reference signal describes a trajectory to follow. A reference signal is a goal, whereas the control input is the command chosen to achieve that goal. [ref]
Asking for a temperature of does not specify a heater command of : dimensional analysis distinguishes the units of these quantities. The reference signal and the measured output must describe the same quantity in the same coordinates before their difference can define a tracking error. An unconverted sensor readout cannot be subtracted directly from a desired temperature.
For example, let be temperature rise above a fixed ambient baseline , and write the desired temperature as . A Heaviside step function gives for , with before the change: the target rise becomes at once, although a physical controlled plant generally cannot jump at once. A ramp reference signal instead asks the temperature rise to increase at a fixed rate . If and , the first reference signal asks for a lasting one-degree rise; the second asks for a one-degree rise at and a two-degree rise at . Here degree Celsius differences are also kelvin differences. A target that grows forever is useful for analytical tests but eventually exceeds a real heater's range.
To make the distinction between a goal and a command concrete, the examples in steady-state tracking error use the first-order low-pass filter as a new controlled plant, rather than the accumulating controlled plant used in the earlier feedback demonstration:
The time constant has units of time, has units of plant output per control input, and the additive control disturbance has the units of . This is an equivalent output offset, not the rate-like in . A heater interpretation uses in degree Celsius differences and in percentage points of heater command, with in degrees per percentage point; the control input is a command, not a direct measurement of delivered power. With exact measurement, holding requires . Thus even a constant reference signal can require a persistent nonzero control input. The numerical model below is deliberately synthetic and has no fitted laboratory parameters.
Connect one physical resistor across an ideal lossless capacitor , with no applied source and an ideal readout that draws no electric current. Represent the resistor by its noisy Thévenin equivalent. The resulting RC low-pass frequency response gives the stationary capacitor voltage power spectral density
Its mean square follows by the effective noise bandwidth calculation:
This is also the equipartition theorem: the mean capacitor energy is . Raising makes the low-frequency spectrum taller and its width narrower in exactly compensating proportions. It changes the equilibration time , not the final variance.
The integral to infinity belongs to the ideal classical model. For a physical electrical circuit, assume thermal equilibrium, linear components with negligible extra dissipation, and a corner far below the frequencies where parasitic effects, frequency-dependent electrical resistance, or invalidate the model; the omitted high-frequency tail must be negligible. A finite-bandwidth instrument measures only part of this variance. Its input capacitance contributes to the total , and its finite input electrical resistance introduces both loading and its own thermal electrical noise.
For and , the predicted equilibrium variance is and the root mean square voltage is . At , and . At , and ; the equilibrium root mean square voltage is identical.
Figure 1.
Changing electrical resistance changes noise timescales but preserves
. Original calculated comparison at and . Left: the physical parallel resistor and capacitor and their noisy Thévenin equivalent. Middle: equilibrium one-sided voltage spectra for and have the same full integral. Right: mean square voltage grows toward the same equilibrium value from an exactly discharged initial capacitor. Curves describe ensemble expected values, not measured or simulated voltage traces. Details.
To obtain the growth curve, prepare exactly and connect the resistor at while maintaining its temperature. The RC low-pass filter's causal impulse response is for . Classical ideal white noise with two-sided power spectral density gives
The factor here is the two-sided time-domain noise intensity; the equilibrium spectral integral above used the one-sided . They give the same long-time limit. The variance approaches equilibrium with time constant , although a deterministic initial voltage decays with time constant . If the capacitor already starts in thermal equilibrium, its initial variance supplies the missing exponentially decaying term and the total variance stays constant.
Protein kinase activity adds modification, while protein tyrosine phosphatase activity can remove it. An original one-site model shows why a receptor phosphorylation assay cannot determine either process from a single modified fraction. Let be the fraction of a fixed receptor population carrying one selected phosphotyrosine residue. Stipulate constant effective first-order rate constants , with ATP availability absorbed into , and no receptor synthesis, degradation or trafficking:
Here is the steady state fraction and the time constant. The ordinary differential equation has the exponential function solution
The balance equation and a single relaxation step give this result; it is an illustrative model, rather than a fit or a claim about every RTK.
With stipulated rate constants , and . Doubling both rate constants leaves but changes to . The balanced phosphate turnover per site,
doubles from to . Identical steady state modification therefore permits different catalytic fluxes. Removing the addition term instead gives ; this illustrates why a perturbation that stops protein kinase activity can expose removal dynamics, subject to the model's unchanged-removal assumption. A maintained steady state with ongoing ATP consumption is not chemical equilibrium.
RL circuit 2026-10-14
An RL circuit combines electrical resistance and self-inductance . In the series model, the same electric current flows through both components. An ideal battery maintains terminal electric potential difference ; the drops in the direction of positive electric current obey Ohm's law and the inductor relation:
This is Kirchhoff's voltage law with the inductor's terminal drop, rather than its oppositely signed induced electromotive force. For constant , , and , the first-order linear differential equation has solution
The time constant measures exponential decay of the difference from the final electric current. [ref].
For a stipulated example, take electrical resistance , self-inductance , terminal electric potential difference , and initial electric current . Then
With measured numerically in seconds, this is in amperes. At one time constant,
and the terminal electric potential difference drops are
Their sum is exactly . At , the inductor takes the entire terminal electric potential difference while the electric current is still zero. At steady electric current, the inductor's terminal drop is zero and the entire drop is across the electrical resistance.
Figure 1.
Electric current and terminal electric potential difference drops during an RL circuit response
. The battery maintains across electrical resistance and self-inductance in series. The dashed vertical line marks time constant ; uses the positive terminal-drop convention. Details.
The final magnetic energy stored by the inductor is . At one time constant, the electric current has reached about of its final value, but its magnetic energy has reached only about :
This difference follows from the square of electric current in the inductor's magnetic energy.
For release, take at the start of a separate experiment with initial electric current . Remove the battery from the model while maintaining a closed path through the same electrical resistance and inductor. Then
Now : the induced electromotive force supports the continuing positive electric current as its magnitude decreases. This closed-loop model describes exponential decay, not opening the sole path and imposing an instantaneous jump to zero electric current. The inductor's magnetic energy and the power of Joule heating are
Their time constant is . The total energy dissipated by Joule heating is exactly the initial magnetic energy:
The second diagram in inductor shows this energy transfer. A finite closed path lets us predict both how long the electric current persists and how much energy becomes heat. [ref].
RL low-pass filter Created 2026-10-18 Updated 2026-10-20
An RL circuit with the output measured across its resistor is a first-order low-pass filter. For an ideal inductor, an ideal voltage source and a high-resistance output instrument, the electrical impedances and give
The gain approaches one at low angular frequency and at high angular frequency; at it is . The time-domain equation is and , so the time constant is . Instrument loading, source resistance and nonideal inductance can change the response.
Shunting inhibition 2026-10-14
Shunting inhibition is suppression of another input's effect through increased synaptic conductance: opening additional ion channels lowers effective electrical resistance, reducing the voltage response to a given added electric current. When the added reversal potential equals the resting membrane potential, there can be no voltage change at rest and still be inhibition of a simultaneous input. See 1.4 Limitations of the Leaky Integrate-and-Fire Model.
Use the passive synaptic potential model below with leak , and capacitance . A sustained inward test electric current of alone gives a steady voltage increase of , so . Effective electrical resistance is and the time constant is .
Now hold a shunting conductance open with . With no test electric current, rest stays at . With the same test electric current, total conductance becomes , so the steady voltage increase is only : . Effective electrical resistance is and the time constant is . A faster response can therefore be smaller. At the original rest the shunting electric current is zero; at it is , opposing the voltage excursion.
The factor of three reduction applies to this constant-electric current, constant-conductance example. A finite pulse, conductance-driven input or spatially extended dendrite need not have the same reduction. The interesting feature is conditional action: the shunt can leave a quiet neuron's voltage unchanged yet reshape its response when another input arrives. A model that represents all inhibition as a fixed outward electric current misses this behavior.
Signaling adaptation 2026-10-13
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 jumps from to at and remains , a Heaviside step function. Let be a fast activating branch, a slower inhibitory branch, and a readout. The activities are dimensionless variables; physical time is in seconds. Choose time constants second and seconds, with initial conditions :
Both branches read directly. Neither branch depends on , and does not change : the subtraction only defines the readout. Thus the slow opposing influence is feedforward, with no closed feedback control loop.
For , solving the two ordinary differential equations gives
Writing for the numerical time in seconds, these are , and . Initially , then rises faster than , producing a pulse. At long times and both approach , so returns to despite sustained . For this upward step , hence . 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.
Figure 1.
A sustained unit input and a transient readout in an incoherent feedforward loop toy
. Original computed curves for time constants of and seconds, zero initial conditions, and normalized dimensionless variables. The upper panel shows the sustained input and the two branches; the lower panel shows the readout , its maximum near seconds, and its approach to the zero baseline. Values and time scales are chosen for illustration, not fitted to a measured signal transduction pathway. Details.
Find the pulse maximum by differentiating the exponential functions. With numerical time in seconds, the derivative is
The stationary condition is , so the natural logarithm gives
Before this time the derivative is positive; afterwards it is negative. Together with and , this proves that it is the unique pulse maximum. As another numerical check, at seconds , and : the slow branch is still catching up.
The subtraction 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 leaves for . 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 2026-10-13
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 activate relay . If inhibits , increased eventually produces an opposing influence: negative feedback. If activates , increased 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.
Figure 1. . Original causal schematics. Arrows activate; flat-ended lines inhibit. The upper loops return from downstream relay to earlier relay . In the lower incoherent feedforward loop, input drives both a fast activating branch and a slower inhibitory branch ; their contributions meet at readout . There is no return arrow from to either branch or to . Shapes and connections illustrate the distinction, rather than a measured pathway. Details.
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 be a fixed normalized input, a normalized response, a dimensionless feedback gain, and a time constant. Suppose response removal contributes and a return proportional to the response contributes . The ordinary differential equation is
Setting the derivative to zero gives the equilibrium point . A displacement obeys , so its decaying exponential function establishes asymptotic stability. For a unit input switched on at , initial condition , and second,
The negative feedback reduces the eventual response from to , 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]
An ideal voltage source and series resistor feed parallel branches and . Thévenin's theorem gives and . Starting with zero inductor current, and . Switch-on heat in is . Disconnecting the source leaves the stored magnetic energy to dissipate through , with time constant .
A spark gap can rapidly discharge a charged capacitor through an inductor. After excitation, the stored energy runs down through Joule heating, radiation and other losses. A key held down can permit repeated discharges: each discharge starts a damped train of rapid oscillations. The spark repetition frequency, the electrical oscillation frequency within each train, and the rate of keyed symbols are different quantities. Repetition does not make each train a maintained sinusoid.
An original lumped model replaces the conducting spark gap and other losses by a constant total electrical resistance in an isolated series RLC circuit. For electric charge on the capacitor, the Kirchhoff cycle law gives
For , its characteristic roots are , with , proving the stated damped form. The electric current has the same exponential decay scale. Its wave envelope falls by per cycle; the energy envelope falls twice as fast in the exponent. A maintained sinusoid instead requires continuing input power balancing losses.
With , and , the model gives , Q factor , amplitude time constant , and damped frequency . These are chosen component values. A real spark gap switches nonlinearly, an antenna radiates, and coupling exchanges energy between circuits; one constant- free decay does not model the complete transmitter. Rapid damping also broadens the frequency spectrum, which helps explain why short trains make sharp receiver selection difficult.
Synaptic potential 2026-10-14
A synaptic potential is a postsynaptic change in membrane potential caused by activity at a synapse. Unlike the regenerative action potential, a passive synaptic potential is graded: its size depends on input strength, the initial membrane potential and the receiving cell membrane's properties. Local voltage changes can decay before reaching the soma; their timing and overlap help determine whether an action potential follows. See Ion Channels and the Electrical Properties of Membranes.
Depolarization means that the interior becomes less negative relative to the exterior; hyperpolarization means that it becomes more negative. Excitation and inhibition instead describe whether an input promotes or opposes postsynaptic action potential generation. An inhibitory postsynaptic potential can be depolarizing, and shunting inhibition can occur with no voltage change at rest. The receptor, permeating ions, ionic gradients and reversal potential matter. See Purves and colleagues: excitation, inhibition and depolarizing inhibition.
A worked passive model treats the receiving membrane as one compartment with capacitance , leak conductance , and activated synaptic conductances :
All electric currents are inward-positive. Each term vanishes at its reversal potential ; . With constant parameters, the ordinary differential equation has
Here is total conductance, is effective electrical resistance, is the steady state and the time constant. Without imposed electric current, is a conductance-weighted average of the reversal potentials. This approximation omits spatial structure, changing ionic gradients and active voltage-gated ion channels; it calculates passive voltage responses, not actual action potential waveforms. See 3.1 Synapses for the mechanism behind the terms; the inward-positive convention here reverses that source's synaptic-current sign.
Take the same illustrative leak and capacitance as under shunting inhibition, start at and set . Opening with gives and . Its initial inward electric current is . Adding with changes to and . This is shunting inhibition of the same excitatory drive.
For comparison, with alone gives : hyperpolarization. With instead, it alone gives : depolarization. Yet adding this depolarizing conductance to the excitatory drive lowers that drive's from to . At an illustrative firing threshold , its electric current is , opposing threshold approach. The combined trace initially rises faster than the excitatory drive alone, before settling below it. Thus the sign of a voltage change at rest alone does not determine excitation or inhibition; a claim about firing must specify input strength, timing and the neuron's state.
Figure 1. . Computed illustrations, not recordings. Left: inward electric current for a conductance with four reversal potentials. Upper right: constant-conductance steps at time zero in the one-compartment model with , and . Lower right: the same test electric current with and without of added conductance at . All plotted trajectories are passive; no action potential generation or reset is included. Details.
This deliberately synthetic differential equation model illustrates response direction and lag. Every plotted number is chosen for teaching; none is a biological measurement or a fitted human parameter. Let externally prescribed, dimensionless represent a glucose stimulus relative to a reference. Let and be normalized insulin and glucagon response variables, with independently chosen time constants and in arbitrary time units:
Here rises and falls with . Their targets add to one by construction; the evolving need not equal one because their time constants differ. Start with and . Set from to , from to , and otherwise. Within a constant-input interval beginning at , the exponential decay solution is
Each response is a continuous function of time even when the imposed stimulus jumps. The chosen starting values stay between zero and one, since each solution is a weighted arithmetic mean of its interval's starting value and target.
Figure 1.
Synthetic insulin and glucagon responses to imposed glucose steps
. Original illustration of the synthetic insulin-glucagon response model. All curves are constructed examples in arbitrary units. Top: imposed dimensionless stimulus. Middle: responses with unequal time constants. Bottom: illustrative liver release and tissue uptake tendencies, which do not determine the imposed stimulus. No experimental data or clinical ranges are represented. Details.
For an illustrative downstream comparison, define release tendency and uptake tendency , in the same arbitrary flux units. At the reference target, . At the high-stimulus targets, and ; at the low-stimulus targets, and . These chosen functions encode greater uptake with insulin and greater release with glucagon, moderated by insulin; they are not established enzyme kinetics. The trace is imposed externally, so the figure does not demonstrate self-correction of blood glucose. A closed model would also need a conservation of mass equation, nutrient input, finite glycogen stores, and additional tissue interactions. The qualitative restoring loops are shown separately in the insulin circuit.