← SciSim / Neuroscience
· · SciSim ·

Hebbian Learning & LTP Full Model

BCM Rule · Ca²⁺ Dynamics · Eligibility Trace · RL

🧠 Tier: Graduate · Synaptic Plasticity
Version:
§ 01
Interactive Simulation
Variable 1
Variable 2
Input
Output
Count
0
State
Time
0
ms
Playback
Preset
Parameters
Parameter 15.0
Parameter 25.0
Input I5.0
T sim (ms)500
Overlays
§ 02
The Idea, Step by Step
▸ From "fire together, wire together" to a precise timing rule
START — Friends who always show up together
Think of two people who keep arriving at the same parties. Before long, when you picture one, you picture the other — your brain has linked them. Neurons do the same thing. When one cell repeatedly helps fire the cell it connects to, the connection between them gets stronger. Donald Hebb's one-line summary, from 1949, is the most quoted sentence in neuroscience: cells that fire together, wire together. Strengthening a connection is called LTP (long-term potentiation); weakening it is LTD (long-term depression). That tug-of-war is how a brain stores a memory.
BUILD — Putting a number on the connection
Give the connection a weight $w$: a big $w$ means a strong link, a small $w$ a weak one. Call the sending cell's activity $x$ and the receiving cell's activity $y$. The simplest Hebb rule just says the weight grows in proportion to how active both cells are at once: $\Delta w = \eta\,x\,y$, where $\eta$ is a small learning rate. Worked number: with $\eta = 0.01$, $x = 1$, and $y = 1$ (both firing), one paired episode gives $\Delta w = 0.01\times1\times1 = 0.01$ — a tiny nudge. Repeat it 50 times and the weight climbs by about $0.5$. Learning is many small nudges, not one big jump.
DEEPEN — Timing decides the sign (STDP)
Real synapses care about order, not just togetherness. Let $\Delta t = t_{\text{post}} - t_{\text{pre}}$ be the gap between the receiving spike and the sending spike. Spike-timing-dependent plasticity follows a two-sided window: $\Delta w = A_+\,e^{-\Delta t/\tau_+}$ when $\Delta t>0$ (sender fired first, so it plausibly caused the spike → LTP), and $\Delta w = -A_-\,e^{\Delta t/\tau_-}$ when $\Delta t<0$ (sender fired after → LTD). The windows $\tau_+,\tau_-$ are roughly $20$ ms. Biophysically, the sign is set by $\text{Ca}^{2+}$ through the NMDA receptor: a big, fast calcium pulse triggers CaMKII and LTP, while a smaller prolonged rise recruits phosphatases and LTD. The BCM rule $\frac{dw}{dt}=\phi(y,\theta_M)\,x$ with a sliding threshold $\theta_M\propto\langle y^2\rangle$ keeps the whole process from running away. In this sim the sliders map as Parameter 1 $\to A_+$ (LTP strength), Parameter 2 $\to A_-$ (LTD strength), and Input I $\to$ the presynaptic firing rate.
TRY THIS in the sim above
(1) Push Parameter 1 high and Parameter 2 low — pairings that arrive in the "pre-before-post" order pile on potentiation and the weight climbs toward its ceiling. (2) Flip them — large Parameter 2, small Parameter 1 — and watch depression win, dragging the weight down. (3) Raise Input I (firing rate) and notice the BCM threshold slide up, so the same pairings that used to strengthen the synapse start to weaken it — the circuit self-stabilises instead of saturating.
§ 03
Equation Derivation
▸ BCM Rule & Ca²⁺-Based Plasticity
$$\frac{dw_i}{dt}=\phi(y,\theta_M)x_i, \quad \phi(y,\theta_M)=y(y-\theta_M), \quad \theta_M\propto\langle y^2\rangle$$ $$\Delta w=\begin{cases}+LTP & [Ca]>\theta_{LTP}\\ -LTD & \theta_{LTD}<[Ca]<\theta_{LTP}\\ 0 & [Ca]<\theta_{LTD}\end{cases}$$
STEP 1 — BCM Sliding Threshold
BCM (1982): φ changes sign at θ_M. Above θ_M: LTP. Below: LTD. θ_M slides: dθ_M/dt = (y²-θ_M)/τ_θ → proportional to ⟨y²⟩. Stabilising: high activity → high θ_M → LTD dominates → activity decreases → θ_M drops. Predicts ocular dominance plasticity and selectivity emergence.
STEP 2 — Ca²⁺ Rule
NMDA current → Ca²⁺ influx: d[Ca]/dt = f_NMDA × I_NMDA - [Ca]/τ_Ca. High [Ca] (>θ_LTP≈3μM): CaMKII activation → AMPA insertion → LTP. Intermediate [Ca]: protein phosphatase (PP1,PP2B) → AMPA internalisation → LTD. Ca²⁺ rule unifies STDP and BCM under a common mechanism.
STEP 3 — Three-Factor RL Rule
Δw = d(t) × e(t). e(t): eligibility trace (STDP-like synapse activity signal, τ_e ≈ 1–5 s). d(t): dopamine reward signal (broadcast globally). Bridges temporal gap between action and delayed reward. Striatal D1 neurons: DA → cAMP → PKA → scales LTP magnitude. Neural substrate of reward learning and RL in basal ganglia.
▸ Primary References

Bienenstock Cooper Munro (1982) J Neurosci 2:32; Bhatt et al. (2009) Ann Rev Physiol

Gerstner et al. — Neuronal Dynamics (neuronaldynamics.epfl.ch, free) · Dayan & Abbott — Theoretical Neuroscience (MIT Press) · Izhikevich — DSN (dynamicalsystems.org, free)

§ 04
Frequently Asked Questions
🧠 ConceptualWhat is the core mathematical insight of Hebbian Learning & LTP Full Model?
The Hebbian Learning & LTP Full Model framework provides a rigorous bridge between microscopic neural parameters and macroscopic observables. The key insight is that complex emergent behaviours — oscillations, memory, coding precision — can be derived from simple mathematical rules governing individual neurons and their connections. The quantitative framework enables testable predictions about circuit function from experimentally measurable quantities.
Key: Mathematical rigour connects single-neuron parameters to network behaviour. Every emergent property (oscillation, attractor, code) is derivable from the microscopic rules.
🌍 AppliedHow is Hebbian Learning & LTP Full Model used in neurotechnology?
The Hebbian Learning & LTP Full Model framework is directly applied in brain-computer interfaces and computational medicine. Mathematical models derived from these principles optimise stimulation parameters for DBS, predict drug effects on network dynamics, decode neural signals for BCI applications, and guide the design of neuromorphic hardware. The quantitative framework enables model-based optimisation rather than empirical trial-and-error approaches in clinical settings.
Key: Direct applications in BCI design, DBS parameter optimisation, pharmacological modelling, and neuromorphic chip architecture.
🔬 SimulationHow do I use the simulation for Hebbian Learning & LTP Full Model?
Configure the parameter sliders on the right panel, select a preset to set up a canonical example, then press Play. The Main View shows the primary mathematical relationship; the Phase Plane shows geometric structure; the Time Series shows temporal evolution; the Parameter Sweep reveals sensitivity to inputs. All displayed quantities correspond directly to terms in the equations in Section 2. Export the data as CSV for further analysis.
Key: Each tab maps to a specific aspect of the mathematics. Use Parameter Sweep to find bifurcation points and Phase Plane to understand stability.
💡 Non-ObviousWhat is counterintuitive about Hebbian Learning & LTP Full Model?
The most surprising result in Hebbian Learning & LTP Full Model research is that small parameter changes can produce qualitatively different dynamics (bifurcations). A system that is stable for one parameter value can suddenly oscillate, burst, or become chaotic with a tiny perturbation. This sensitivity to parameters is the hallmark of nonlinear systems and has profound implications: the same neural circuit can perform completely different computations depending on its neuromodulatory state, which shifts parameters near bifurcation boundaries.
Key: Bifurcations: tiny parameter changes → qualitative behavioural change. Neuromodulators exploit this by operating near bifurcation points to switch neural circuit function.
📐 ComputationalWhat numerical methods and pitfalls apply to Hebbian Learning & LTP Full Model?
The appropriate numerical method depends on stiffness: for near-spike dynamics (HH, AdEx) use RK4 with dt ≤ 0.025 ms or implicit methods; for mean-field equations (WC, Amari) use Euler or RK4 with dt ≤ 0.5–1 ms. Always verify accuracy by halving dt and checking solution convergence. The most common errors: wrong units (mixing mV, ms, nS), missing initial conditions (not setting w₀=b×v₀ in Izhikevich), and incorrect boundary conditions in Fokker-Planck or cable equation simulations.
Key: Verify stability by halving dt. Most errors: wrong units, incorrect initial conditions, violated model assumptions. Check the model validity range before trusting any result.
§ 04 Resources
§ 05
Misconceptions & Common Errors
Sub-block A — Conceptual
"The Hebbian Learning & LTP Full Model model is only theoretical with no experimental support."
The Hebbian Learning & LTP Full Model framework has been extensively validated. Key predictions have been confirmed in electrophysiology, imaging, and pharmacological experiments. Modern neuroscience uses these models as quantitative tools integrated with the experimental cycle — for prediction, interpretation, and guidance of experiments.
📖 Bienenstock Cooper Munro (1982) J Neurosci 2:32; Bhatt et al. (2009) Ann Rev Physiol
Sub-block B — Numerical
Applying Hebbian Learning & LTP Full Model equations outside their range of validity (wrong N, dt, or approximation regime).
Every model has a validity regime. Always check: (1) Is N large enough for mean-field? (2) Is dt small enough (halve dt and verify solution converges)? (3) Are units consistent (mV, ms, nS throughout)? (4) Are initial conditions correct? When in doubt, compare to direct simulation of the full system.
🔍 Why: Violating assumptions gives quantitatively or qualitatively wrong results — and the error is often invisible without explicit validation.
§ 05 References