2025-11-13T03:01:10.895552

Simulating neuronal dynamics in fractional adaptive exponential integrate-and-fire models

Fikl, Jhinga, Kaslik et al.
We introduce an efficient discretization of a novel fractional-order adaptive exponential (FrAdEx) integrate-and-fire model, which is used to study the fractional-order dynamics of neuronal activities. The discretization is based on extension of L1-type methods that can accurately handle the exponential growth and the spiking mechanism of the model. This new method is implicit and uses adaptive time stepping to robustly handle the stiff system that arises due to the exponential term. The implicit nonlinear system can be solved exactly, without the need for iterative methods, making the scheme efficient while maintaining accuracy. We present a complete error model for the numerical scheme that can be extended to other integrate-and-fire models with minor changes. To show the feasibility of our approach, the numerical method has been rigorously validated and used to investigate several different spiking oscillations of the model. We observed that the fractional-order model is capable of predicting biophysical activities, which are interpreted through phase diagrams describing the transition from one firing type to another. This simple model shows significant promise, as it has sufficient expressive dynamics to reproduce several features qualitatively from a biophysical dynamical perspective.
academic

Simulating neuronal dynamics in fractional adaptive exponential integrate-and-fire models

基本信息

  • 论文ID: 2404.13554
  • 标题: Simulating neuronal dynamics in fractional adaptive exponential integrate-and-fire models
  • 作者: Alexandru Fikl, Aman Jhinga, Eva Kaslik, Argha Mondal
  • 分类: physics.bio-ph
  • 发表时间: 2024年5月30日 (arXiv v2)
  • 论文链接: https://arxiv.org/abs/2404.13554

摘要

本文提出了一种新型分数阶自适应指数(FrAdEx)积分发放模型的高效离散化方法，用于研究神经元活动的分数阶动力学。该离散化基于L1型方法的扩展，能够准确处理模型的指数增长和脉冲机制。新方法采用隐式格式并使用自适应时间步长来稳健处理指数项引起的刚性系统。隐式非线性系统可以精确求解，无需迭代方法，使该方案在保持精度的同时提高了效率。本文提供了完整的数值格式误差模型，可以稍作修改后扩展到其他积分发放模型。数值方法经过严格验证，并用于研究模型的多种脉冲振荡。研究发现，分数阶模型能够预测生物物理活动，通过相图描述从一种发放类型到另一种的转变。这个简单模型显示出巨大潜力，具有足够的表达动力学来从生物物理动力学角度定性再现多种特征。

研究背景与动机

要解决的问题

本研究旨在解决神经元动力学建模中的三个核心问题：

  1. 记忆效应建模：传统整数阶导数模型无法捕捉神经元的非局部、历史依赖行为和长期依赖关系
  2. 数值求解挑战：分数阶积分发放模型涉及弱奇异核的积分算子、指数增长的刚性系统以及状态依赖的脉冲条件
  3. 理论空白：文献中尚无关于状态依赖脉冲的分数阶脉冲系统的理论和数值分析

问题重要性

  1. 生物学真实性：自适应指数(AdEx)积分发放模型是简化神经元模型与复杂真实神经元动力学之间的重要桥梁
  2. 多时间尺度动力学：神经元电压轨迹通常表现出跨越多个时间尺度的动力学，这些动力学遵循幂律行为，可以通过分数阶微分方程有效建模
  3. 计算效率：与Hodgkin-Huxley型详细电导模型相比，积分发放模型更适合直接数学分析和数值实现

现有方法局限性

  1. 经典模型不足：传统的漏积分发放(LIF)模型基于整数阶导数，无法模拟记忆效应和连续整合过去神经元状态
  2. 分数阶方法缺失：虽然已有分数阶LIF模型，但针对更复杂的AdEx模型的分数阶扩展及其数值方法尚未充分研究
  3. 数值方法挑战：高阶方法对于需要完全非均匀网格的分数阶微分方程尚未充分发展

研究动机

作者受到以下因素驱动：

  1. 最近有研究使用局部分形导数扩展AdEx模型，但缺乏严格的数学分析
  2. 分数阶Caputo导数提供了更严格的数学框架
  3. 需要开发能够处理指数增长、状态依赖脉冲和分数阶记忆效应的高效数值方法

核心贡献

  1. 提出FrAdEx模型：首次引入分数阶自适应指数积分发放模型，将Caputo分数导数引入AdEx模型，能够捕捉神经元的记忆效应和历史依赖行为
  2. 开发高效数值方法：设计了基于L1型方法的隐式自适应时间步长格式，通过Lambert W函数实现精确求解（无需迭代）
  3. 完整误差分析：提供了针对状态依赖脉冲分数阶系统的完整误差模型（截断误差和全局误差），可扩展到其他积分发放模型
  4. 数值验证与应用：在多个基准模型（PIF、LIF、FrAdEx）上验证了方法的一阶收敛性，并展示了模型能够再现多种神经元发放模式（颤动、快速脉冲、强直性脉冲等）
  5. 开源实现：基于Python的pycaputo库提供开源实现，促进可复现性

方法详解

任务定义

研究的核心任务是数值求解以下分数阶脉冲微分方程系统：

FrAdEx模型（非量纲化形式）：

\frac{d^{\alpha_1} V}{dt^{\alpha_1}} = I - (V - E_L) + \exp(V) - w \\ \tau_w \frac{d^{\alpha_2} w}{dt^{\alpha_2}} = a(V - E_L) - w \end{cases}$$ **重置条件**： $$\text{if } V > V_{peak} \text{ then } \begin{cases} V \leftarrow V_r \\ w \leftarrow w + b \end{cases}$$ 其中： - $V(t)$：膜电位 - $w(t)$：自适应变量 - $0 < \alpha_i < 1$：分数阶 - $I(t)$：外部电流 - 参数：$(g_L, E_L, \Delta_T, V_T, \tau_w, a, b, V_r, V_{peak})$ **挑战**： 1. Caputo分数导数定义为弱奇异核的积分算子 2. 指数项导致的刚性系统 3. 状态依赖的脉冲时刻$\{\tau_m\}$未知 ### 模型架构 #### 1. Caputo分数导数的分段定义 对于分段绝对连续函数$y \in PAC([0,T];\mathbb{R}^2)$，分段Caputo导数定义为： $${}^{PC}D^{\alpha}_{0+}[y](t) = \frac{1}{\Gamma(1-\alpha)} \left[\sum_{j=0}^{m-1} \int_{\tau_j}^{\tau_{j+1}} \frac{y'(s)}{(t-s)^{\alpha}} ds + \int_{\tau_m}^{t} \frac{y'(s)}{(t-s)^{\alpha}} ds\right]$$ 其中$t \in (\tau_m, \tau_{m+1}]$，$\tau_m$为脉冲时刻。 #### 2. L1型离散化 在每个区间$[t_n, t_{n+1}]$上使用线性插值： $$y(s) \approx \frac{t_{n+1}-s}{t_{n+1}-t_n}y_n^+ + \frac{s-t_n}{t_{n+1}-t_n}y_{n+1}^-$$ 得到离散格式： $$\sum_{k=0}^{n} d_{n+1,k} \odot \frac{y_{k+1}^- - y_k^+}{\Delta t_k} = f(t_{n+1}, y_{n+1}^-)$$ 权重系数： $$d_{n+1,k} = \frac{(t_{n+1}-t_k)^{1-\alpha} - (t_{n+1}-t_{k+1})^{1-\alpha}}{\Gamma(2-\alpha)}$$ #### 3. Lambert W函数精确求解 将隐式方程改写为： $$\hat{V}^- + c_2 = c_3 \exp(\hat{V}^-)$$ 其中$c_2, c_3$为已知系数。解为： $$\hat{V}^- = -c_2 - W[-c_3 \exp(-c_2)]$$ $$\hat{w}^- = c_0 \hat{V}^- + c_1$$ 其中$W(\cdot)$为Lambert W函数，可以无需迭代直接计算。 #### 4. 自适应时间步长 基于误差指示器： $$\chi_{n+1} = \|\Gamma(1+\alpha)\| \frac{(t_{n+1}-t_n)^{\alpha}}{t_{n+1}^{\alpha} - t_n^{\alpha}} \frac{\|y_{n+1}-y_n\|}{\|y_n\|}$$ 归一化到$[0,1]$区间： $$\chi_{n+1} = \frac{\hat{\chi}_{n+1} - \chi_{min}}{\chi_{max} - \chi_{min}}$$ **自适应策略**： - 若$0 < \chi_{n+1} < 1$：保持步长$\Delta t_{n+1} = \theta \Delta t_n$ - 若$\chi_{n+1} < 0$：增大步长$\Delta t_{n+1} = \rho \Delta t_n$ - 若$\chi_{n+1} > 1$：减小步长$\Delta t_{n+1} = \sigma \Delta t_n$ #### 5. 脉冲时刻估计 当$\hat{V}_{n+1}^-$为复数时（Lambert W函数参数超出实数域），通过求解： $$c_3(\Delta t^*) \exp(-c_2(\Delta t^*) + 1) = 1$$ 得到最大允许时间步长$\Delta t_{Lambert}$，保证解的实数性。 ### 技术创新点 1. **Lambert W精确求解**： - 避免了Newton-Raphson迭代 - 直接计算解析解 - 提高计算效率 2. **自适应时间步长与脉冲处理的结合**： - 在指数增长区域自动细化网格 - 通过Lambert W函数约束保证数值稳定性 - 一阶精度估计脉冲时刻 3. **分段L1方法的误差分析**： - 首次针对状态依赖脉冲的分数阶系统 - 证明了全局误差为$O(\Delta t_{max})$ - 跳跃项在误差分析中精确抵消 4. **通用性设计**： - 方法适用于一般积分发放模型 - 最小化对重置条件的假设 - 易于扩展到其他模型 ## 实验设置 ### 数据集 本文不涉及传统意义的数据集，而是通过数值模拟验证： 1. **PIF模型**：具有已知解析解，用于收敛性验证 - 参数：$C=100$ pFms$^{\alpha-1}$, $I=160$ pA, $V_{peak}=0$ mV, $V_r=-48$ mV - 测试$\alpha \in \{0.5, 0.75, 0.95\}$ 2. **LIF模型**：用于验证自适应算法 - 附加参数：$g_L=3$ nS, $E_L=-50$ mV - $\alpha = 0.85$ 3. **FrAdEx模型**：完整模型验证 - 多组参数配置，模拟不同神经元发放模式 - $\alpha \in [0.9, 0.999]$ ### 评价指标 1. **相对$\ell^2$误差**： $$E(x, x_{ref}) = \frac{\|x - x_{ref}\|_2}{\|x_{ref}\|_2}$$ 2. **脉冲时刻误差**：比较数值脉冲时刻$\{t_n\}$与精确/参考脉冲时刻$\{\tau_m\}$ 3. **收敛阶**：通过对数-对数图验证$O(\Delta t_{max})$收敛性 4. **计算效率**： - 迭代次数与计算时间的关系 - 自适应vs固定步长的效率对比 ### 对比方法 1. **固定步长L1方法**：作为基线对比自适应方法的效率 2. **精确解**（PIF模型）：验证数值方法的正确性 3. **自收敛**（FrAdEx模型）：使用极细网格作为参考解 ### 实现细节 1. **编程环境**：Python + numpy + scipy 2. **开源库**：pycaputo（作者开发的分数阶微积分库） 3. **自适应参数**： - 安全因子：$\theta = 1.0$ - 减小因子：$\sigma = 0.5$ - 增大因子：$\rho \in [1.5, 2.0]$ - 最小步长：$\Delta t_{min} = 10^{-5}$ - 初始步长：$\Delta t_0 = 10^{-2}$ 4. **误差界限**： - $\chi_{min} = \{2^{-k} | k=0,...,7\}$ - $\chi_{max} = \{2^{1-k} | k=0,...,7\}$ ## 实验结果 ### 主要结果 #### 1. PIF模型收敛性（图3） - **测试配置**：$\alpha \in \{0.5, 0.75, 0.95\}$，$\Delta t \in \{10^{-2}, 5\times10^{-3}, 10^{-3}, 5\times10^{-4}\}$ - **结果**：所有分数阶均达到**一阶收敛** $O(\Delta t_{max})$ - **验证**：与精确解析解对比，误差从$10^{-3}$降至$10^{-5}$ - **脉冲数量**：6个脉冲时刻，所有脉冲位置均准确估计 #### 2. LIF模型自适应性（图4） - **配置**：$\alpha=0.85$，$\chi_{max} \in \{2^2, 2^{-2}, 2^{-6}\}$ - **观察**： - 在平滑区域，时间步长呈阶梯式增长（从$10^{-5}$增至$10^{-1}$） - 接近脉冲时，步长急剧减小以捕捉指数增长 - 更严格的误差界限（$\chi_{max}=2^{-6}$）导致更小的平均步长 - **效率**：自适应方法能够在保证精度的同时减少计算量 #### 3. FrAdEx模型自收敛（图5） - **配置**：$\alpha=0.9$，演化至$T=50$（无量纲），5个脉冲 - **收敛性**： - 全局误差表现出清晰的一阶收敛 - 每个脉冲时刻的误差随$\Delta t_{max}$线性减小 - 后续脉冲的误差不会累积恶化 - **误差范围**：从$10^{-1}$（粗网格）降至$10^{-3}$（细网格） #### 4. 计算效率对比（图6） - **渐近复杂度**： - 自适应和固定步长方法均为$O(N^2)$（分数阶记忆项的固有复杂度） - 实验验证了理论预测的二次缩放 - **效率优势**： - 达到$10^{-3}$相对误差： - 自适应方法：约0.1秒 - 固定步长：约1秒（**10倍差距**） - 自适应方法在相同精度下显著减少计算时间 ### 神经元发放模式再现（第8节） #### 实验组1（图7）：参数集1 - **$\alpha=0.999$**：**颤动（Chattering）** - 表现为密集间隔的脉冲簇 - 相图显示紧密的极限环 - **$\alpha=0.98$**：**快速脉冲伴宽脉冲后电位（Fast spiking with broad SAP）** - 脉冲后出现小曲率的后电位 - 相图显示更宽松的轨迹 - **$\alpha=0.93$**：**强直性脉冲伴尖锐SAP（Tonic spiking with sharp SAP）** - 膜电位在快速下冲后单调增加 - 无适应性，规律发放 #### 实验组2（图8）：参数集2 - **$\alpha=0.999$**：**强直性发放伴宽SAP** - 规律的动作电位放电 - **$\alpha=0.98$**：**规律脉冲伴尖锐SAP** - **$\alpha=0.93$**：**脉冲频率适应（Spike frequency adaptation）** - 初始脉冲间隔短，随后间隔逐渐增大 - 相图显示螺旋收敛模式 #### 实验组3（图9）：参数集3 - **$\alpha=0.999$**：**内在簇发（Intrinsic bursting）** - 初始簇发后转为规律单脉冲 - **$\alpha=0.95$**：**强直性脉冲伴尖锐SAP** - **$\alpha=0.9$**：**规律脉冲伴频率适应** ### 关键发现 1. **分数阶对发放模式的影响**： - $\alpha \to 1$：更接近整数阶行为，表现出复杂的簇发和颤动 - $\alpha$减小：发放模式趋于规律化，表现出更强的记忆效应 - 分数阶作为控制参数，可调节神经元的发放类型 2. **数值方法的鲁棒性**： - 在所有参数配置下均保持稳定 - 准确捕捉相图中的转变 - 能够处理从快速脉冲到慢适应的多时间尺度动力学 3. **生物物理意义**： - FrAdEx模型能够定性再现多种已知的神经元发放模式 - 相图清晰展示了不同发放类型之间的转变 - 简单模型具有足够的表达能力 ## 相关工作 ### 积分发放模型的发展 1. **经典模型**： - Lapicque (1907)：漏积分发放(LIF)模型 - Izhikevich (2003)：二次积分发放模型 - Fourcaud-Trocmé等 (2003)：指数积分发放模型 2. **AdEx模型**： - Brette & Gerstner (2005)：首次提出AdEx模型 - Naud等 (2008)：系统研究AdEx的发放模式 - Touboul & Brette (2008)：AdEx的动力学和分岔分析 ### 分数阶神经元模型 1. **分数阶LIF**： - Teka等 (2014)：首次提出分数阶LIF模型，研究脉冲时序适应 - Weinberg & Santamaria (2017)：历史依赖的神经元活动 2. **分形导数扩展**： - Souza等 (2024)：使用局部分形导数扩展AdEx模型 - 本文区别：使用Caputo分数导数，提供更严格的数学框架 ### 分数阶微分方程数值方法 1. **L1方法**： - Li & Zeng (2015)：经典L1方法专著 - Li & Cai (2019)：分数阶积分和导数的理论与数值逼近 2. **非均匀网格方法**： - Li等 (2017)：高阶数值方法 - Yang & Zeng (2023)：修正的L1方法 3. **自适应方法**： - Jannelli (2020)：求解分数阶微分方程的自适应过程 - 本文扩展：结合状态依赖脉冲的自适应方法 ### 脉冲微分方程 1. **固定时刻脉冲**： - Wang等 (2016)：脉冲分数阶微分方程综述 - 已有大量研究成果 2. **状态依赖脉冲**： - Lakshmikantham等 (1994)：变时刻脉冲微分方程的比较原理 - **本文贡献**：首次针对状态依赖脉冲的分数阶系统提供理论和数值分析 ### 本文优势 1. **数学严格性**：使用Caputo导数而非局部分形导数 2. **数值效率**：Lambert W精确求解，无需迭代 3. **完整误差分析**：首次针对状态依赖脉冲的分数阶系统 4. **通用性**：方法可扩展到其他积分发放模型 ## 结论与讨论 ### 主要结论 1. **模型能力**： - FrAdEx模型成功再现了多种神经元发放模式（颤动、强直性脉冲、频率适应、簇发等） - 分数阶$\alpha$可作为控制参数调节发放行为 - 较小的$\alpha$增强记忆效应，导致更规律的发放模式 2. **数值方法**： - 提出的L1型方法在所有测试案例中达到一阶收敛 - Lambert W求解策略显著提高计算效率 - 自适应时间步长成功处理多时间尺度动力学 3. **理论贡献**： - 首次为状态依赖脉冲的分数阶系统提供完整误差分析 - 证明了全局误差为$O(\Delta t_{max})$ - 方法框架可扩展到其他积分发放模型 ### 局限性 1. **收敛阶限制**： - 当前方法为一阶精度 - 对于分数阶微分方程，扩展到高阶方法需要显著努力 - 指数增长模型的高阶脉冲时刻估计尚不清楚 2. **脉冲累积误差**： - 理论分析假设脉冲数量$m$不太大 - 当$m \max(y(\tau_j^+) - y(\tau_j^-)) = O(\Delta t_{max}^{-1})$时估计可能失效 - 长时间演化可能需要额外考虑 3. **参数选择**： - 自适应算法的参数$(\chi_{min}, \chi_{max}, \theta, \sigma, \rho)$需要经验调整 - 缺乏自动选择这些参数的鲁棒方法 4. **模型验证**： - 定性再现发放模式，但缺乏与实验数据的定量对比 - 需要更多生物物理验证 5. **计算复杂度**： - $O(N^2)$复杂度对长时间模拟仍然昂贵 - 虽然Fourier方法可改进至$O(N\log N)$，但不清楚如何应用于非均匀不连续系统 ### 未来方向 1. **高阶方法**： - 开发二阶或更高阶的分数阶积分发放模型数值方法 - 探索在指数增长区域的高阶脉冲时刻估计 2. **神经元网络**： - 扩展到耦合FrAdEx神经元网络 - 研究网络中的同步和集体行为 - 处理耦合情况下Lambert W求解的挑战 3. **理论分析**： - 完整的动力学和分岔分析 - 不同参数区域的发放模式分类 - 分数阶对神经元动力学影响的理论刻画 4. **实验验证**： - 与真实神经元记录数据的定量对比 - 参数拟合和模型选择 - 确定生物学上合理的$\alpha$值范围 5. **算法改进**： - 开发自适应参数的自动选择策略 - 探索快速算法（如快速Fourier变换）在非均匀网格上的应用 - 并行化实现 ## 深度评价 ### 优点 1. **创新性强**： - 首次系统研究分数阶AdEx模型的数值方法 - Lambert W精确求解策略新颖且高效 - 首次针对状态依赖脉冲的分数阶系统提供完整误差分析 2. **理论严谨**： - 完整的数学框架（定理6.1和6.2） - 严格的误差估计和收敛性证明 - 清晰的假设和适用条件 3. **实验充分**： - 三个模型（PIF、LIF、FrAdEx）的系统验证 - 收敛性、自适应性、效率的全面测试 - 多组参数配置的神经元发放模式展示 4. **实用价值**： - 开源实现（pycaputo库） - 方法可扩展到其他模型 - 详细的算法伪代码（Algorithm 1） 5. **写作清晰**： - 结构合理，逻辑清晰 - 数学表述准确 - 图表丰富且信息量大 ### 不足 1. **高阶方法缺失**： - 仅实现一阶方法，对于某些应用可能精度不足 - 未探索高阶方法的可行性 2. **生物学验证不足**： - 仅定性再现发放模式 - 缺乏与实验数据的定量对比 - 未讨论$\alpha$的生物学意义 3. **参数敏感性分析不足**： - 自适应算法参数的选择缺乏系统指导 - 未充分探索参数空间 4. **计算效率**： - $O(N^2)$复杂度对大规模网络模拟仍然挑战 - 未探索快速算法的可能性 5. **理论分析局限**： - 脉冲累积误差的条件较强 - 未提供动力学和分岔的理论分析 ### 影响力 1. **对领域的贡献**： - 为分数阶神经元建模提供了重要工具 - 推动状态依赖脉冲分数阶系统的研究 - 为其他积分发放模型提供方法论框架 2. **实用价值**： - 开源代码促进可复现性和应用 - 方法可直接应用于神经科学研究 - 为计算神经科学提供新的建模工具 3. **可复现性**： - 详细的算法描述 - 开源实现（pycaputo） - 明确的参数设置 4. **潜在应用**： - 神经元网络模拟 - 脑疾病建模（如癫痫） - 神经形态计算 ### 适用场景 1. **理想场景**： - 需要建模神经元记忆效应的研究 - 多时间尺度神经元动力学研究 - 小到中等规模的神经元网络模拟 2. **不适用场景**： - 需要极高精度的应用（当前一阶方法） - 大规模网络实时模拟（$O(N^2)$复杂度） - 需要与实验数据精确拟合的定量研究 3. **潜在扩展**： - 其他类型的积分发放模型（如Izhikevich模型） - 随机分数阶神经元模型 - 时变参数系统 ## 参考文献（关键文献） 1. **Brette & Gerstner (2005)**: Adaptive exponential integrate-and-fire model - AdEx模型的原始论文 2. **Teka et al. (2014)**: Fractional leaky integrate-and-fire model - 分数阶LIF模型 3. **Li & Zeng (2015)**: Numerical methods for fractional calculus - 分数阶微积分数值方法经典教材 4. **Jannelli (2020)**: Adaptive procedure for fractional differential equations - 自适应时间步长方法 5. **Fečkan et al. (2012)**: Impulsive fractional differential equations - 脉冲分数阶微分方程理论基础 --- **总体评价**：这是一篇高质量的计算神经科学论文，在分数阶神经元建模和数值方法方面做出了重要贡献。理论严谨、实验充分、实用性强。主要局限在于当前仅实现一阶方法和缺乏生物学定量验证。该工作为分数阶积分发放模型的研究奠定了重要基础，具有较高的学术价值和应用潜力。