quadratic integrate and fire

{{more citations needed|date=March 2023}}

The quadratic integrate and fire (QIF) model is a biological neuron model that describes action potentials in neurons. In contrast to physiologically accurate but computationally expensive neuron models like the Hodgkin–Huxley model, the QIF model seeks only to produce action potential-like patterns by ignoring the dynamics of transmembrane currents and ion channels. Thus, the QIF model is computationally efficient and has found ubiquitous use in computational neuroscience.{{Citation |last=Fourcaud-Trocmé |first=Nicolas |title=Integrate and Fire Models, Deterministic |date=2013 |url=https://doi.org/10.1007/978-1-4614-7320-6_148-1 |encyclopedia=Encyclopedia of Computational Neuroscience |pages=1–9 |editor-last=Jaeger |editor-first=Dieter |place=New York, NY |publisher=Springer |language=en |doi=10.1007/978-1-4614-7320-6_148-1 |isbn=978-1-4614-7320-6 |access-date=2023-03-01 |editor2-last=Jung |editor2-first=Ranu}}

An idealized model of neural spiking is given by the autonomous differential equation,

: $\backslash frac\{dx\}\{dt\}\; =\; x^2\; +\; I$

where $x$ represents the membrane voltage and $I\backslash geq\; 0$ represents an input current. A solution to this differential equation is the function,{{cite book |last=Ermentrout |first1=Bard |last2=Terman |first2=David |title=Mathematical Foundations of Neuroscience |edition=1 |publisher=Springer |date=July 1, 2010 |isbn=978-0-387-87708-2 |doi=10.1007/978-0-387-87708-2 }}

: $x(t)\; =\; \backslash sqrt\{I\}\; \backslash tan(\backslash sqrt\{I\}\; t\; +\; c\_0),$

where $c\_0$ is an arbitrary shift dependent on the initial condition $x(0)$ (specifically by the formula $c\_0\; =\; \backslash arctan(x(0)/\backslash sqrt\{I\})$). This solution "blows up" in finite time, namely at $t=n\backslash pi/(2\backslash sqrt\{I\})-c\_0$ for all $n\; \backslash in\backslash mathbb\{N\}$, which resembles the rhythmic action potentials generated by neurons stimulated by some input current. Thus a "spike" is said to have occurred when the solution reaches positive infinity. Just after this point in time, the solution resets to negative infinity by definition.

When implementing this model in a numerical simulation, a threshold crossing value ($V\_t$) and a reset value ($V\_r$) is assigned, so that when the solution rises above the threshold, $x(t)\; \backslash geq\; V\_t$, the solution is immediately reset to $V\_r$.

The above equation is directly related to an alternative form of the QIF model,

: $\backslash frac\{dv\}\{dt\}\; =\; -\backslash frac\{v(1-v)\}\{\backslash tau\_m\}\; +\; I$,

where $\backslash tau\_m$ is the membrane time constant.