KA1_HM1992#
- class braincell.channel.KA1_HM1992(size, g_max=Quantity(30., 'mS / cm^2'), temp=Quantity(309.15, 'K'), q10_p=1.0, temp_ref_p=Quantity(309.15, 'K'), q10_q=1.0, temp_ref_q=Quantity(309.15, 'K'), V_sh=Quantity(0., 'mV'), name=None)#
Huguenard & McCormick 1992 A-type potassium current (IA1).
The first of the two components into which (Huguenard & McCormick, 1992) [1] splits the rapidly inactivating transient potassium current \(I_A\) of thalamic relay neurons, with \(p^4 q\) HH gating and an ohmic driving force:
\[\begin{split}\begin{aligned} p_\infty &= \frac{1}{1 + \exp(-(V' + 60) / 8.5)} \\ \tau_p &= \frac{1} {\exp((V' + 35.8) / 19.7) + \exp(-(V' + 79.7) / 12.7)} + 0.37 \\ q_\infty &= \frac{1}{1 + \exp((V' + 78) / 6)} \\ \tau_q &= \begin{cases} \left[\exp((V' + 46) / 5) + \exp(-(V' + 238) / 37.5)\right]^{-1} & V' < -63 \\ 19 & V' \geq -63 \end{cases} \end{aligned}\end{split}\]where \(V' = (V - V_{sh}) / \mathrm{mV}\) and the time constants are in milliseconds, further scaled per gate by
gate_phi().- Parameters:
size (
int|Sequence[int] |integer|Sequence[integer]) – Channel state shape.g_max (
Array|ndarray|bool|number|bool|int|float|complex|Quantity|Callable) – Maximal conductance density. Defaults to30.0 mS/cm2(see Notes).temp (
Array|ndarray|bool|number|bool|int|float|complex|Quantity) – Absolute temperature driving the Q10 factors, default 36 degrees Celsius.q10_p (
Array|ndarray|bool|number|bool|int|float|complex|Quantity|Callable) – Q10 scaling factor for the activation gate, default1.0.temp_ref_p (
Array|ndarray|bool|number|bool|int|float|complex|Quantity) – Reference temperature forq10_p, default 36 degrees Celsius.q10_q (
Array|ndarray|bool|number|bool|int|float|complex|Quantity|Callable) – Q10 scaling factor for the inactivation gate, default1.0.temp_ref_q (
Array|ndarray|bool|number|bool|int|float|complex|Quantity) – Reference temperature forq10_q, default 36 degrees Celsius.V_sh (
Array|ndarray|bool|number|bool|int|float|complex|Quantity|Callable) – Threshold shift applied to both gates’ rates, default0.0 mV.
See also
KA2_HM1992The other \(I_A\) component of the same model; it differs from this class only in
p_infandg_max(see Notes).KK2A_HM1992Slowly inactivating \(I_{K2}\) component of the same model.
KK2B_HM1992The other \(I_{K2}\) component.
Notes
The paper’s published abstract enumerates exactly four currents, which map onto BrainCell as: \(I_T\) ->
CaT_HM1992; \(I_A\) ->KA1_HM1992/KA2_HM1992; \(I_{K2}\) ->KK2A_HM1992/KK2B_HM1992; and \(I_h\) ->HCN_HM1992.One divergence between the shipped pair and the paper’s own description, recorded rather than corrected. The abstract says \(I_A\) “was modeled by assuming two components with different time constants of inactivation”. In BrainCell the two components’ inactivation is identical:
KA1_HM1992andKA2_HM1992carry the sameq_inf, the same piecewisetau_qand even the sametau_p, and differ only in the activation midpoint and slope (-60 / 8.5 here, -36 / 20 there) and in the default conductance (30 against 20 mS/cm2). The \(I_{K2}\) pair does differ in inactivation, as the paper describes for that current.The
tau_qbranches do not meet: at \(V' = -63\) the lower branch evaluates to about 23.4 ms while the class returns the constant 19 ms. The discontinuity is in the source parameterisation, not introduced here; the code applies the expression strictly below -63 and the constant at and above it.Both Q10 factors default to
1.0, so with the shipped defaults neither gate is temperature-scaled at any temperature; each gate has its ownq10/temp_refpair and is scaled independently. Theg_maxdefault is a BrainCell value: the attribution for this key was established from the paper’s current inventory, not by comparing conductance densities.References