Simulación de dinámicas neuronales en modelos fraccionarios de integración y disparo exponencial adaptativo

Información Básica

• ID del Artículo: 2404.13554
• Título: Simulating neuronal dynamics in fractional adaptive exponential integrate-and-fire models
• Autores: Alexandru Fikl, Aman Jhinga, Eva Kaslik, Argha Mondal
• Clasificación: physics.bio-ph
• Fecha de Publicación: 30 de mayo de 2024 (arXiv v2)
• Enlace del Artículo: https://arxiv.org/abs/2404.13554

Resumen

Este artículo propone un método de discretización eficiente para un nuevo modelo fraccionario de integración y disparo exponencial adaptativo (FrAdEx), utilizado para estudiar la dinámica fraccionaria de la actividad neuronal. La discretización se basa en una extensión del método tipo L1, capaz de manejar con precisión el crecimiento exponencial y los mecanismos de disparo del modelo. El nuevo método utiliza un esquema implícito con pasos de tiempo adaptativos para manejar de manera robusta los sistemas rígidos causados por términos exponenciales. El sistema no lineal implícito puede resolverse exactamente sin métodos iterativos, mejorando la eficiencia mientras se mantiene la precisión. El artículo proporciona un modelo de error de formato numérico completo, extensible a otros modelos de integración y disparo con ligeras modificaciones. El método numérico se valida rigurosamente y se utiliza para estudiar múltiples oscilaciones de disparo del modelo. Se encuentra que el modelo fraccionario puede predecir actividad biofísica, describiendo transiciones entre tipos de disparo a través de diagramas de fase. Este modelo simple demuestra un potencial considerable, con suficiente expresividad dinámica para reproducir cualitativamente múltiples características desde una perspectiva de dinámica biofísica.

Antecedentes de Investigación y Motivación

Problemas a Resolver

Esta investigación aborda tres problemas fundamentales en el modelado de dinámicas neuronales:

1. Modelado de Efectos de Memoria: Los modelos tradicionales con derivadas de orden entero no pueden capturar el comportamiento no local, dependiente del historial y las dependencias a largo plazo de las neuronas
2. Desafíos de Resolución Numérica: Los modelos fraccionarios de integración y disparo implican operadores integrales con núcleos débilmente singulares, sistemas rígidos con crecimiento exponencial y condiciones de disparo dependientes del estado
3. Vacío Teórico: Falta de teoría y análisis numérico en la literatura para sistemas fraccionarios de pulsos con impactos dependientes del estado

Importancia del Problema

1. Realismo Biológico: El modelo adaptativo exponencial de integración y disparo (AdEx) es un puente importante entre modelos neuronales simplificados y la dinámica neuronal real compleja
2. Dinámicas de Múltiples Escalas Temporales: Las trayectorias de voltaje neuronal típicamente exhiben dinámicas que abarcan múltiples escalas temporales, siguiendo comportamientos de ley de potencia que pueden modelarse efectivamente mediante ecuaciones diferenciales fraccionarias
3. Eficiencia Computacional: Los modelos de integración y disparo son más adecuados para análisis matemático directo e implementación numérica en comparación con modelos de conductancia detallados tipo Hodgkin-Huxley

Limitaciones de Métodos Existentes

1. Insuficiencia de Modelos Clásicos: Los modelos tradicionales de integración y disparo con fugas (LIF) basados en derivadas de orden entero no pueden simular efectos de memoria e integración continua de estados neuronales pasados
2. Ausencia de Métodos Fraccionarios: Aunque existen modelos LIF fraccionarios, la extensión fraccionaria a modelos AdEx más complejos y sus métodos numéricos aún no se han investigado suficientemente
3. Desafíos de Métodos Numéricos: Los métodos de orden superior para ecuaciones diferenciales fraccionarias que requieren mallas completamente no uniformes aún no se han desarrollado completamente

Motivación de la Investigación

Los autores están impulsados por los siguientes factores:

1. Investigaciones recientes utilizan derivadas fractales locales para extender el modelo AdEx, pero carecen de análisis matemático riguroso
2. La derivada fraccionaria de Caputo proporciona un marco matemático más riguroso
3. Necesidad de desarrollar métodos numéricos eficientes que manejen crecimiento exponencial, pulsos dependientes del estado y efectos de memoria fraccionaria

Contribuciones Principales

1. Propuesta del Modelo FrAdEx: Introducción del primer modelo fraccionario de integración y disparo exponencial adaptativo, incorporando la derivada fraccionaria de Caputo en el modelo AdEx para capturar efectos de memoria e historial de neuronas
2. Desarrollo de Método Numérico Eficiente: Diseño de un esquema implícito con pasos de tiempo adaptativos basado en el método tipo L1, con resolución exacta mediante la función Lambert W (sin iteración)
3. Análisis de Error Completo: Proporciona un modelo de error completo para sistemas fraccionarios con pulsos dependientes del estado (error de truncamiento y error global), extensible a otros modelos de integración y disparo
4. Validación Numérica y Aplicación: Validación de convergencia de primer orden en múltiples modelos de referencia (PIF, LIF, FrAdEx), demostrando que el modelo puede reproducir múltiples patrones de disparo neuronal (temblor, disparo rápido, disparo tónico, etc.)
5. Implementación de Código Abierto: Implementación de código abierto basada en la biblioteca Python pycaputo, promoviendo la reproducibilidad

Explicación Detallada del Método

Definición de la Tarea

La tarea central es resolver numéricamente el siguiente sistema de ecuaciones diferenciales fraccionarias con pulsos:

Modelo FrAdEx (forma adimensional):