← SciSim / Neuroscience
· · SciSim ·

Cable Theory — Space Constant λ & Time Constant τ

🧠 Tier: Standard Undergraduate · Dendrites · Myelination · Rall
Version:
§ 01
Interactive Simulation — Cable Theory & Dendritic Voltage Spread
λ space const
0.71
mm
τ_m time const
20.0
ms
V at λ
36.8
%V₀
Electrotonic L
1.41
Attenuation
-12.3
dB
d diameter
1.0
μm
Sim Time
0
ms
Playback
Preset
Cable Parameters
d diameter (μm)1.0
R_m (Ω·cm²)20000
R_a (Ω·cm)100
C_m (μF/cm²)1.0
Length (mm)2.0
V₀ injection (mV)20
T sim (ms)100
Overlays
§ 02
The Idea, Step by Step
▸ From a fading whisper to the cable equation
START — everyday picture
Pass a whisper down a long line of people. Each person hears it a little fainter, and by the far end it is almost gone. A dendrite — the thin branch that carries a synapse's signal toward the cell body — behaves the same way. A synapse delivers a small voltage "bump," and that bump leaks away through the membrane as it travels. Signals from nearby synapses arrive loud at the soma; signals from far-away synapses arrive faint. The dendrite is a leaky cable, not a perfect wire.
BUILD — the space constant and one number
How far can a signal travel before it fades? That distance is the space constant $\lambda$ (lambda): the length over which a steady voltage shrinks to about 37% (that is $1/e$) of where it started. The simplest rule for a long cable with a steady push at one end is $V(x) = V_0\,e^{-x/\lambda}$. Worked number: if $\lambda = 0.5$ mm and a synapse sits $0.5$ mm away (so $x=\lambda$), the soma feels $e^{-1}\approx 37\%$ of the original bump. Move the synapse to $1$ mm ($x=2\lambda$) and only $e^{-2}\approx 14\%$ survives. The membrane also takes time to charge: the time constant $\tau_m \approx 20$ ms sets how quickly the voltage rises and settles.
DEEPEN — the full cable equation and the sliders
The space constant comes from the cable's geometry: $\lambda = \sqrt{R_m\,d/(4R_a)}$, so a fatter dendrite reaches farther, $\lambda \propto \sqrt{d}$. Combining space and time gives the cable equation $\lambda^2\,\partial^2 V/\partial x^2 - \tau_m\,\partial V/\partial t = V - V_{\text{rest}}$, with $\tau_m = R_m C_m$. Measuring a cable in "$\lambda$ units" gives its electrotonic length $L = \ell/\lambda$: at $L=1$ the far end sees 37%, at $L=2$ about 14%, at $L=3$ only 5%. In the sim, the d diameter slider sets $\lambda$ through $\sqrt{d}$; R_m and R_a push $\lambda$ up or down; C_m controls charging speed (and $\tau_m$); Length sets $L$ and therefore how badly distal signals attenuate.
TRY THIS in the sim above
1) Drag d diameter from 1 μm to 5 μm and watch $\lambda$ grow by $\sqrt{5}\approx 2.2\times$ — the signal reaches much farther. 2) Crank R_m up toward 100,000 Ω·cm² (the effect of myelin) and see $\lambda$ jump so the bump arrives at the soma nearly intact. 3) Increase Length past about $3\lambda$ and notice the far end voltage nearly vanishes ($L>3 \Rightarrow <5\%$) — that is why distal synapses are intrinsically weaker.
§ 03
Equation Derivation
▸ The Cable Equation (Lord Kelvin 1855, Wilfrid Rall 1959)
$$\boxed{\lambda^2 \frac{\partial^2 V}{\partial x^2} - \tau_m \frac{\partial V}{\partial t} = V - V_{\text{rest}}}$$ $$\lambda = \sqrt{\frac{R_m\,d}{4\,R_a}} \;\text{(space const, mm)}, \qquad \tau_m = R_m C_m \;\text{(time const, ms)}$$

Steady-state solution (semi-infinite cable, step input at x=0):

$$V(x) = V_0\,e^{-x/\lambda} \qquad \text{(exponential decay in space)}$$

Electrotonic length L = ℓ/λ. Voltage attenuation = e^{−L}. In dB: 20 log₁₀(e^{−L}) = −8.686 L dB.

SymbolMeaningUnitTypical Value
λSpace constant — 1/e decay lengthmm0.1–2 mm (dendrites)
τ_m = R_m C_mMembrane time constantms5–30 ms
R_mSpecific membrane resistanceΩ·cm²5,000–100,000
R_aSpecific axial resistanceΩ·cm50–200
dDendrite/axon diameterμm0.1–20 μm
L = ℓ/λElectrotonic length (dimensionless)0.1–3
▸ Derivation Steps
STEP 1 — Cable Equivalent Circuit
A dendrite segment δx has: axial resistance r_a δx = 4R_a/(πd²) × δx (resists axial current); membrane capacitance c_m δx = C_m πd δx; membrane resistance r_m δx = R_m/(πd) per δx. KVL along axial direction + KCL at node → cable PDE as δx → 0. The PDE has λ = √(r_m/r_a) and τ_m = r_m c_m as the two natural length and time scales.
STEP 2 — λ ∝ √d (Diameter Dependence)
λ = √(R_m d/(4R_a)) ∝ √d. Thicker cables: lower axial resistance (more cross-section area) but also more membrane area (more leak). The net effect: λ ∝ √d. Doubling d → λ×√2. Giant squid axon (d~1 mm): λ~10 mm. Thin cortical dendrite (d~0.1 μm): λ~0.1 mm. This is why large axons transmit APs faster and further.
STEP 3 — Rall's 3/2 Power Branching Rule
A dendritic branch point is impedance-matched (no reflection) when d_parent^(3/2) = Σ d_children^(3/2). Under this condition, the full branching tree can be collapsed to a single equivalent cylinder of the same electrotonic length. This is the foundation of all compartmental neuron models (NEURON, GENESIS).
STEP 4 — Myelination: R_m↑, C_m↓
Myelin wraps axon ~150 times: R_m_eff ×100, C_m_eff ÷100. λ_myelinated ≈ 10× λ_unmyelinated. Saltatory conduction: AP jumps node-to-node (1 mm spacing). Conduction velocity: myelinated v = 6d (m/s), unmyelinated v = k√d m/s. MS destroys myelin → λ decreases → conduction block → neurological symptoms.
STEP 5 — Finite-Difference Simulation
Discretise cable into N compartments: C_m δx (dV_i/dt) = (V_{i-1}−2V_i+V_{i+1})/(r_a δx) − V_i/r_m δx + I_inj. This is an N×N tridiagonal ODE system. Stability (explicit Euler): dt < δx² τ_m/(2λ²). Crank-Nicolson is unconditionally stable and preferred for production simulations.

▸ Worked Example — λ and Electrotonic Length

d = 1 μm = 10⁻⁴ cm, R_m = 20,000 Ω·cm², R_a = 100 Ω·cm, C_m = 1 μF/cm², length ℓ = 0.5 mm.

$$\lambda = \sqrt{\frac{20000 \times 10^{-4}}{4 \times 100}} = \sqrt{\frac{2}{400}} = 0.0707 \text{ cm} = 0.707 \text{ mm}$$ $$\tau_m = R_m C_m = 20000 \times 10^{-6} \text{ s} = 20 \text{ ms}$$ $$L = \ell/\lambda = 0.5/0.707 = 0.707, \quad V(\ell)/V_0 = e^{-0.707} = 49\%, \quad \text{Attenuation} = -6.1 \text{ dB}$$
▸ References
[Rall59]Rall — Branching dendritic trees and motoneuron membrane resistivity, Exp. Neurology 1:491, 1959
[JW95]Johnston & Wu — Foundations of Cellular Neurophysiology, MIT Press, 1995. Ch. 3–5
[Koc99]Koch — Biophysics of Computation, Oxford, 1999. Ch. 2–3
[LH05]London & Häusser — Dendritic computation, Annu. Rev. Neurosci. 28:503, 2005
§ 04
Frequently Asked Questions
🔬 SimulationWhat does the Cable V(x,t) tab show — how does voltage spread along a dendrite?
The Cable tab shows a 2D heatmap: x-axis = position along dendrite (mm), y-axis = time (ms), colour = voltage (warm=depolarised, cool=rest). When current is injected at x=0, voltage rises locally and spreads with exponential decay in space (e^{−x/λ}) and exponential decay in time (e^{−t/τ_m}). Unlike a simple RC (uniform decay everywhere), the cable shows a spatiotemporal wavefront. A synapse at x=2λ delivers e^{−2} ≈ 14% of its voltage to the soma. The Lambda tab shows V(x) at steady state as the exponential decay curve — drag the diameter slider and watch λ change.
Key takeaway: Cable shows 2D spatiotemporal decay — both distance and time reduce voltage. A synapse at 2λ from soma delivers only 14% voltage. Distal synapses are intrinsically weaker.
🧠 ConceptualWhy does dendrite diameter affect λ — why does λ ∝ √d?
λ = √(R_m d/(4R_a)) ∝ √d. Thicker dendrite: lower axial resistance per unit length r_a ∝ 1/d² (larger cross-section carries more current) but higher membrane resistance per unit length r_m ∝ 1/d (more perimeter = more leak). Net: r_m/r_a ∝ (1/d)/(1/d²) = d, so λ ∝ √d. Doubling d increases λ by √2 = 1.41×. The squid giant axon (d~1 mm) has λ~10 mm (signals travel far); a thin cortical dendrite (d~0.1 μm) has λ~0.1 mm (signals decay quickly). Myelination effectively increases d by adding insulating layers.
Key takeaway: λ ∝ √d because axial conductance (∝d²) beats membrane leak (∝d) — the net effect is √d. Thicker or myelinated = larger λ = better signal transmission.
🌍 AppliedHow does cable theory explain multiple sclerosis symptoms?
MS destroys myelin, which normally increases R_m×100 and decreases C_m×100, giving λ_myelinated ≈ 10× λ_unmyelinated and conduction velocity 50–70 m/s vs 0.5–2 m/s. When myelin is lost: λ decreases dramatically, C_m increases, action potential conduction fails (conduction block). MS symptoms — visual loss (optic nerve), weakness (corticospinal tract), sensory loss — correspond to demyelination in specific white matter tracts. The cable equation directly predicts the minimum diameter for successful AP conduction after demyelination. Treatment: 4-aminopyridine (K⁺ channel blocker) prolongs AP duration to compensate for reduced λ.
Key takeaway: Myelination: R_m×100, C_m÷100, λ×10, v×10. MS destroys myelin → λ decreases → conduction block → symptoms. Cable theory predicts all these quantitatively.
💡 Non-ObviousWhy do distal dendritic synapses have the same somatic EPSP as proximal ones?
Passive cable theory predicts strong location dependence — a distal synapse at L=2 delivers only 14% voltage to soma. Yet Magee & Cook (2000) showed somatic EPSPs are similar regardless of synaptic location. The resolution: (1) Synaptic strength scaling — distal synapses are physically larger (more AMPA receptors) to compensate for attenuation (homeostatic normalisation); (2) Active dendritic conductances — NMDA receptors, Ca²⁺, h-current amplify distal EPSPs locally; (3) Dendritic spikes — strong distal inputs trigger local Na⁺/Ca²⁺ spikes that propagate to soma with less attenuation. Passive cable theory is incomplete — active dendrites change the computation fundamentally.
Key takeaway: Passive cable predicts strong location dependence; biology compensates via synaptic scaling + active conductances. Dendrites are active computational units, not passive cables.
📐 ComputationalWhat is the von Neumann stability criterion for cable simulation?
For explicit Euler on the cable equation (diffusion coefficient D = λ²/τ_m): dt < δx²/(2D) = δx² τ_m/(2λ²). For λ=0.5 mm, τ_m=20 ms, δx=0.1 mm: dt_max = 0.01×20/0.5 = 0.4 ms — safe with dt=0.01 ms. Violating stability: spatial oscillations appear and grow exponentially. Crank-Nicolson (semi-implicit) is unconditionally stable for any dt and is preferred for production cable simulations. Rule: always check if δx < λ (otherwise the spatial grid cannot resolve the decay profile).
Key takeaway: dt < δx²/(2D) for stable explicit cable simulation. Crank-Nicolson is unconditionally stable. Check δx < λ for adequate spatial resolution.
🎓 DeepWhat is Rall's 3/2 power branching rule, and what does it mean for dendritic computation?
Rall's rule: a branch point is impedance-matched when d_parent^(3/2) = Σ d_children^(3/2). Under this condition, no voltage reflection occurs at the branch point, and the entire dendritic tree can be 'unfolded' into a single equivalent cylinder — enormously simplifying analysis. Real dendritic trees often violate the 3/2 rule, creating impedance mismatches that reflect voltage waves back toward the soma. These mismatches allow individual dendritic branches to implement logical AND operations: two inputs must both be active to trigger a local dendritic spike (threshold nonlinearity). This means a single neuron with active dendrites has computational power equivalent to a multi-layer network.
Key takeaway: d_p^(3/2) = Σd_i^(3/2) → impedance-matched branching, equivalent cylinder. Real dendrites violate this, enabling nonlinear synaptic integration and dendritic AND logic.
🧠 ConceptualWhat is electrotonic length L, and what does 'electrotonically compact' mean?
L = ℓ/λ is the physical length normalised by the space constant. L=1: 37% voltage at distal end; L=2: 14%; L=3: 5%. Electrotonically compact (L < 0.3): voltage changes propagate nearly undistorted from one end to the other — any synapse effectively drives the soma. Electrotonically extended (L > 2): distal synapses are electrically isolated. Most cortical pyramidal neurons have L ≈ 1–3 (partial attenuation). Interneurons are typically compact (L < 1) — their small size means any synaptic input reaches the soma effectively. This is why interneurons integrate inputs globally while pyramidal neurons can perform compartmentalised dendritic computations.
Key takeaway: L = ℓ/λ: L < 0.3 = compact (full signal transmission), L > 2 = extended (distal isolation). Pyramidal neurons: L ≈ 1–3. Interneurons: L < 1. Determines whether dendritic computation is local or global.
§ 04 Best Resources
§ 05
Misconceptions & Common Errors
Sub-block A — Conceptual Misconceptions
"Dendrites are passive input cables — they only transmit signals to the soma."
Dendrites are active computational units containing voltage-gated Na⁺, Ca²⁺, K⁺, and h-current channels. Dendritic spikes can be generated in the apical dendrite independent of somatic firing. NMDA receptors provide coincidence detection within individual synaptic spines. Back-propagating APs interact with EPSPs for STDP. Single L5 pyramidal dendrites implement logical AND operations (two simultaneous inputs required to trigger local dendritic spike). A single neuron with active dendrites has computational power equivalent to a multi-layer ANN (Koch & Segev, 1998).
📖 London & Häusser (2005) Annu. Rev. Neurosci. 28:503; Koch — Biophysics of Computation Ch. 11
"At x = λ, the signal is negligible — λ is a hard cutoff distance."
λ is the 1/e (37%) decay length, not a cutoff. At x=λ: V=37% V₀; x=2λ: 14%; x=3λ: 5%; x=5λ: 0.7%. The signal decays asymptotically to zero, never reaching zero exactly. For functional purposes, x > 3λ is approximately negligible, but x=λ to 2λ can still drive significant somatic voltage changes, especially with active dendritic amplification. The key question is whether the attenuated EPSP exceeds the somatic threshold.
📖 Johnston & Wu — Foundations of Cellular Neurophysiology, Ch. 3
"τ_m is the time for the membrane to fully charge — after τ_m ms, V = 100% of steady state."
τ_m is the time to reach 1−1/e ≈ 63% of steady state, not 100%. V(t) = V_∞(1−e^{−t/τ_m}): at t=τ_m → 63%; at 2τ_m → 86%; at 3τ_m → 95%; at 5τ_m → 99.3%. Similarly, decay to rest: at t=τ_m → 37% remaining (not 37% decayed). This 63%/37% convention applies to ALL RC circuits. In patch-clamp experiments, τ_m is measured by fitting an exponential to voltage relaxation after a current step — the time constant is the slope of the log-linear plot.
📖 Johnston & Wu — Foundations of Cellular Neurophysiology, Ch. 4
Sub-block B — Numerical Errors
Computing λ with diameter in μm instead of cm: lambda = sqrt(Rm * d_um / (4*Ra)) gives result in wrong units.
Use consistent CGS units: d in cm (1 μm = 10⁻⁴ cm), R_m in Ω·cm², R_a in Ω·cm. Then λ = √(R_m × d_cm / (4R_a)) gives λ in cm. Multiply by 10 to get mm. Example: d = 1 μm = 10⁻⁴ cm, R_m = 20,000, R_a = 100 → λ = √(20000 × 10⁻⁴ / 400) = √0.005 = 0.0707 cm = 0.707 mm. Using d in μm gives λ in √(Ω·μm·cm) — nonsense units. Always convert to consistent CGS before computing cable parameters.
🔍 Why: Cable parameters use CGS (cm, Ω·cm) but diameters are quoted in μm — unit conversion is easy to forget, giving λ values 100× too small.
Using dt = 0.1 ms and δx = 0.5 mm for cable simulation, violating stability criterion for small λ values.
Check: dt < δx² τ_m/(2λ²). For λ = 0.2 mm (thin dendrite), τ_m = 20 ms, δx = 0.5 mm: dt_max = 0.25×20/(2×0.04) = 62.5 ms — stable. But for δx = 0.5 mm and λ = 0.2 mm: δx/λ = 2.5 — the spatial grid cannot resolve the exponential decay (need δx < λ/3 for accuracy). Rule: δx should be at most λ/3 for spatial resolution. Stability and resolution are both needed: dt for time, δx for space. Use Crank-Nicolson to decouple stability from dt choice.
🔍 Why: Students satisfy the stability criterion but ignore spatial resolution — a stable simulation can still give wrong answers if the spatial grid is too coarse relative to λ.
§ 05 References