Medios excitables y ondas de recuperación

Año 2 · Módulo 15: Dinámica espacial y sistemas continuos

146.1. Una perturbación que avanza y desaparece detrás de sí

Una onda de Fisher–KPP deja a su paso una región que, en el modelo ideal, permanece cerca del estado ocupado. La transición separa dos estados diferentes y progresa porque uno de ellos puede invadir al otro. Hay fenómenos espaciales que producen una imagen distinta. Un impulso atraviesa un tejido; las células que encuentra responden intensamente durante un intervalo breve y después regresan a una condición de reposo. Esta dinámica caracteriza a los medios excitables.

Una membrana nerviosa es un ejemplo histórico. Un estímulo localizado puede desencadenar un impulso que se transmite por regiones contiguas. Cada región responde al cambio eléctrico que recibe de sus vecinas; los mecanismos de recuperación reducen después la excitación y modifican temporalmente la respuesta ante un segundo estímulo. Se obtiene así una secuencia espacial organizada en la que conviven regiones todavía no activadas, regiones activas y regiones en recuperación. Modelos afines también se han aplicado a la actividad eléctrica cardíaca, a ciertas reacciones químicas y a ondas intracelulares de señalización, aunque las variables físicas concretas difieren considerablemente.

La excitabilidad presupone que el reposo puede ser estable. Pequeñas perturbaciones decaen y la actividad espontánea no invade todo el medio por la mera presencia de una fluctuación infinitesimal. Un estímulo suficientemente intenso, de duración y extensión adecuadas, puede llevar al sistema hacia una excursión grande antes de que regrese al equilibrio. El tamaño de la respuesta es entonces desproporcionado respecto de diferencias pequeñas entre estímulos cercanos al umbral. El regreso al reposo no suele ser inmediato, porque una variable lenta conserva la historia reciente de la activación.

Este comportamiento plantea preguntas que no aparecían en el frente monostable del capítulo anterior. ¿Por qué una perturbación local alcanza el umbral de sus vecinas? ¿Cómo se limita en el tiempo un impulso para que la activación no se vuelva permanente? ¿Por qué una segunda onda puede fracasar en una región que hace poco condujo la primera? ¿De qué dependen la velocidad, la anchura y la dirección de propagación? Responder exige estudiar a la vez las reacciones locales y el acoplamiento espacial.

En este capítulo usaremos un modelo de dos variables inspirado en la familia FitzHugh–Nagumo. La primera variable representa una activación rápida; la segunda, una recuperación más lenta. Su combinación produce una descripción matemática manejable de impulsos y refractariedad. El objetivo será identificar, sin pretender una fisiología completa de cada tejido, las condiciones dinámicas bajo las cuales un sistema inicialmente estable puede transportar una señal de amplitud finita y dejar detrás un medio que necesita tiempo para volver a responder.

146.2. La arquitectura mínima de activación y recuperación

Sea $u(x,t)$ una variable rápida de activación y $v(x,t)$ una variable lenta de recuperación. Podemos considerar el sistema unidimensional

con $D > 0$, $0 < a < 1$, $\varepsilon > 0$ pequeño, $\gamma > 0$ y parámetros $\beta$ e $I$ que desplazan la dinámica local. Es una de varias normalizaciones posibles dentro de la familia de modelos de FitzHugh–Nagumo. Las variables y coeficientes son adimensionales en esta presentación, salvo que posteriormente se definan escalas físicas de longitud y tiempo. El término espacial difunde la activación; hemos supuesto que la recuperación no tiene difusión propia, una simplificación que puede levantarse si el fenómeno la requiere.

La no linealidad cúbica tiene una función específica. Para $v \approx 0$ e $I = 0$, el término $f(u) = u(1 - u)(u - a)$ es negativo cuando $u$ está entre cero y $a$ y positivo cuando está entre $a$ y uno; por encima de uno vuelve a ser negativo. Produce un umbral aproximado entre el regreso hacia el reposo y el crecimiento hacia una fase activa. En el sistema completo ese umbral no es un único número universal, porque el valor de $v$, el estímulo y el acoplamiento espacial modifican la respuesta.

El término $- v$ reduce la activación cuando la recuperación se acumula. La ecuación lenta integra parte de la actividad de $u$ y hace que sus consecuencias persistan después de la fase ascendente del impulso. El parámetro $\varepsilon$ fija la separación de escalas temporales; cuando es pequeño, $u$ puede cambiar mucho mientras $v$ apenas varía. Esa asimetría es la fuente geométrica principal de los pulsos de activación y recuperación. Sin una variable que registre el pasado reciente, una descripción escalar puede generar frentes, pero tiene dificultades para representar de manera robusta un pulso aislado de tipo excitable con reposo estable a ambos lados.

Para fijar ideas, tomemos $a = 0,2$, $\gamma = 1$, $\beta = 0$ e $I = 0$. El punto $(u,v) = (0,0)$ es un equilibrio homogéneo. La matriz jacobiana de la cinética local en ese punto es

En dos dimensiones estas condiciones bastan para la estabilidad lineal local del reposo. El equilibrio puede recuperarse después de perturbaciones pequeñas sin impedir excursiones grandes provocadas por estímulos más intensos. Estabilidad lineal y excitabilidad son, por tanto, propiedades compatibles; la primera describe la respuesta infinitesimal y la segunda depende de la geometría no lineal de trayectorias que se apartan de manera considerable del equilibrio.

146.3. Umbral, excursión y retorno en el plano de fases

Sin transporte, el estado de cada punto sigue un sistema de ecuaciones diferenciales ordinarias. Las curvas donde una variable deja instantáneamente de cambiar se llaman nulclinas. Para el sistema anterior con $I$ constante, las nulclinas son

La primera es cúbica y reúne los estados en los que la activación rápida queda momentáneamente equilibrada. La segunda es la recta en la que $v$ deja de cambiar. Sus intersecciones son los equilibrios de la dinámica local. Saber dónde se cortan no explica todavía el recorrido después de un estímulo; para ello hay que considerar la dirección del campo vectorial y la diferencia entre las velocidades de las dos variables.

Cuando $\varepsilon$ es pequeño, las trayectorias muestran movimientos rápidos de $u$ mientras $v$ cambia poco, alternados con períodos en que el estado se mueve más lentamente cerca de ramas de la nulclina cúbica. Un estímulo débil puede alterar $u$ ligeramente y provocar un retorno directo hacia el equilibrio. Un estímulo suficientemente fuerte puede llevar el estado a otra región del plano, desde la cual realiza una excursión de gran amplitud antes de volver. Durante la fase activa, $v$ aumenta; ese aumento modifica el balance de la ecuación rápida y termina facilitando la caída de $u$.

La palabra umbral debe usarse con cuidado. En la aproximación muy rápida con $v$ fijado cerca de cero, el valor $a$ separa dos tendencias del término cúbico escalar. En el sistema dinámico completo puede existir una superficie umbral definida por la organización de trayectorias, cuya posición depende de $v$ y de la forma temporal del estímulo. En el medio espacial aparece además un umbral de extensión, de modo que un estímulo fuerte aplicado sobre una región demasiado estrecha puede disiparse antes de reclutar suficiente tejido vecino. Por ello, la excitabilidad no se reduce a que la variable supere cierta cifra durante un instante.

La geometría también aclara la diferencia entre un medio excitable y uno oscilatorio. En el régimen excitable, el sistema vuelve a un estado de reposo estable después de una perturbación transitoria. En ciertos regímenes de parámetros o forzamientos, ese reposo puede perder estabilidad y aparecer un ciclo límite, con lo que la actividad se repite incluso sin pulsos externos separados. Aunque la estructura de variables pueda ser parecida, las soluciones cualitativas difieren; una onda inducida por un estímulo en un tejido en reposo y una secuencia espontánea de oscilaciones sincronizadas no deberían tratarse como el mismo fenómeno.

En este modelo, la recuperación modifica el futuro sin exigir una memoria externa a las ecuaciones. Dos regiones con idéntico valor actual de $u$ pueden responder de manera diferente porque tienen distintos valores de $v$. Una acaba de recuperar la condición basal; la otra todavía conserva una señal de inhibición tras su activación. Esta propiedad permite representar refractariedad dentro de un sistema determinista de pocas dimensiones.

146.4. Cómo una célula activada recluta a sus vecinas

La difusión de $u$ conecta posiciones que, consideradas aisladamente, regresarían al reposo ante una alteración pequeña. Supongamos que una región localizada recibe un estímulo y entra en la fase activa. La variable de activación adquiere entonces una pendiente espacial pronunciada en sus bordes. El término $Du_{xx}$ redistribuye activación hacia las regiones cercanas. Si el aporte es suficiente para desplazar el estado de una región vecina más allá de su umbral dinámico, esa región desarrolla su propia respuesta intensa. La señal avanza como una sucesión de activaciones locales acopladas.

Esta explicación no requiere que una porción fija de materia viaje desde el origen hasta el extremo del tejido. En el caso eléctrico, la propagación es la transmisión de una perturbación de voltaje y corrientes locales; las moléculas y estructuras celulares que participan tienen movimientos y escalas propios. La analogía difusiva de una ecuación de reacción–difusión captura un acoplamiento espacial efectivo. Para derivarla cuantitativamente a partir de una membrana real harían falta parámetros eléctricos, geometría y leyes constitutivas específicas.

La capacidad de una región activa para estimular a otra está limitada por un problema de fuente y sumidero. La fuente debe ejercer suficiente influencia durante un intervalo adecuado; el sumidero incluye el volumen o el área todavía inactivos, sus pérdidas y su proximidad al umbral. Un estímulo intenso pero extremadamente breve puede fracasar, del mismo modo que un estímulo sostenido distribuido sobre una región pequeña. Por ello, amplitud, duración y anchura espacial contribuyen a determinar si surge un pulso que se propaga o si la perturbación decae localmente.

La comparación de escalas ayuda. Si el término de reacción rápida tiene una tasa característica $k$, una longitud de difusión aproximada es $\ell \sim \sqrt{D/k}$. Una zona activada mucho más estrecha que esa escala puede dispersarse considerablemente antes de sostener la activación de su entorno, aunque el resultado concreto dependa de la no linealidad y de las condiciones iniciales. La estimación dimensional orienta el análisis, pero no es una fórmula universal para el tamaño crítico de un estímulo. Ese tamaño puede obtenerse mediante simulación, análisis de soluciones umbral o experimentación controlada.

El estado de recuperación modifica también la conducción. Si una región tiene $v$ elevado, una perturbación que habría sido suficiente en tejido recuperado puede resultar ineficaz. La propagación se bloquea cuando el frente encuentra material que todavía no puede responder. Esta propiedad hace posibles interacciones muy diferentes de las ondas lineales de amplitud pequeña. En particular, dos pulsos que se aproximan pueden desaparecer al encontrarse, porque cada uno llega a una región cuya historia reciente impide una nueva activación.

Estos comportamientos son pruebas experimentales de que existe recuperación. Un cambio de velocidad ante un segundo impulso cercano, un bloqueo dependiente del intervalo entre estímulos o la desaparición de frentes que colisionan informan sobre la dinámica interna del medio. Sin datos adicionales, sin embargo, no demuestran que el mecanismo microscópico tenga exactamente dos variables o una reacción cúbica. El modelo da un marco conceptual que después debe contrastarse con las propiedades materiales concretas.

146.5. Del frente permanente al pulso viajero

En un medio excitable esperamos a menudo un perfil espacial localizado, en el que la activación está cerca del reposo delante del pulso, aumenta en la región activa y vuelve a valores bajos detrás. La recuperación permanece elevada durante parte de esa región posterior y decae más lentamente. Busquemos una solución viajera con $u(x,t) = U(z)$, $v(x,t) = V(z)$ y $z = x - ct$. Sin estímulo externo, las ecuaciones se convierten en

En una onda de invasión como Fisher–KPP, los dos extremos espaciales correspondían a estados diferentes. Para un pulso excitable aislado ambos extremos pueden corresponder al mismo reposo, con $(U,V)( - \infty) = (U,V)( + \infty) = \left( u_{0},v_{0} \right)$. La trayectoria espacial, considerada como sistema dinámico en la variable $z$, debe alejarse de ese equilibrio y regresar a él. Esa estructura se denomina órbita homoclínica en el sistema de ecuaciones ordinarias del perfil, si se satisfacen las condiciones dinámicas pertinentes. Su existencia no está garantizada para cualquier elección de parámetros; la velocidad se determina junto con el perfil y las condiciones de conexión.

El pulso puede interpretarse, de manera aproximada, como una interfaz delantera de excitación, una zona activa y una interfaz posterior de desactivación. La interfaz delantera encuentra material recuperado, con $V$ cercano al valor basal. La posterior encuentra una región en que $V$ se ha acumulado y modifica la reacción. Las dos interfaces son distintas y deben ajustarse mutuamente para que todo el perfil avance con la misma velocidad. Esta coordinación es una de las razones por las que el problema completo suele resolverse mediante análisis de perturbaciones singulares o métodos numéricos para problemas de contorno.

Podemos derivar una identidad útil integrando la primera ecuación del perfil sobre toda la recta. Si $U$ y $U\prime$ regresan a sus valores de reposo en ambos extremos, los términos derivados no contribuyen y queda

En un pulso completo, el crecimiento neto local de activación debe compensarse con la recuperación integrada sobre la región espacial recorrida. La identidad difiere de la obtenida en Fisher–KPP porque allí $U( - \infty) \neq U( + \infty)$. Sin otras condiciones no fija la velocidad de propagación, pero es una prueba de coherencia para un perfil calculado numéricamente. En intervalos finitos deberán incluirse los términos de borde correspondientes.

Otra precaución es distinguir la velocidad de fase de una onda viajera de la velocidad aparente de un máximo que cambia de amplitud. Durante el establecimiento de un pulso, el frente delantero puede avanzar mientras su parte posterior todavía se organiza. Medir la distancia recorrida por el máximo en ese transitorio puede dar una velocidad que no coincide con la de la solución viajera finalmente formada. La estabilidad del ancho y de la distancia entre perfiles de $U$ y $V$ ayuda a reconocer cuándo la aproximación viajera es adecuada.

146.6. Un frente rápido resoluble y sus límites

La separación de escalas permite estudiar el borde delantero cuando la variable lenta permanece aproximadamente congelada. Si el medio por delante está completamente recuperado y $V \approx 0$, la ecuación rápida se aproxima a $u_{t} = Du_{xx} + u(1 - u)(u - a).$

Es la ecuación bistable cúbica que examinamos como contraste en el capítulo anterior. Para $0 < a < 1$ admite una solución de frente decreciente $U(z) = 1/\left( 1 + e^{kz} \right)$, con $z = x - ct$. La derivada es $U\prime = - kU(1 - U)$ y la segunda derivada $U'' = k^{2}U(1 - U)(1 - 2U)$. Al sustituir ambas en $DU'' + cU\prime + U(1 - U)(U - a) = 0$ y dividir por $U(1 - U)$, obtenemos

Para que la identidad sea válida a lo largo de todo el perfil, el coeficiente de $U$ y el término constante deben anularse por separado. De ahí resultan

La expresión corresponde al frente escalar con la recuperación fija en cero. Es una aproximación para la parte delantera de un pulso excitable cuando la separación de tiempos y el estado del medio justifican esa hipótesis. En el sistema FitzHugh–Nagumo de dos variables, la variable $V$, la forma del borde posterior y la exigencia de construir un pulso completo pueden desplazar sustancialmente la velocidad o impedir la propagación sostenida.

Consideremos $D = 0,5$ y $a = 0,2$ en las unidades adimensionales de este modelo. Entonces $k = 1$, la longitud de variación característica es $1/k = 1$ y

Si un cálculo de la ecuación escalar da una velocidad cercana a $0,3$ después del transitorio y tras refinar la malla, tendremos una verificación del algoritmo en ese límite. Si una simulación completa del sistema de dos variables da otro número, no hay necesariamente un error; primero hay que comprobar cuánto cambia $V$ en el borde delantero y si el pulso realmente alcanzó un régimen estacionario en su sistema de referencia.

El signo de $1 - 2a$ tiene una lectura intuitiva. Un aumento de $a$ reduce la fuerza con que la región activada invade el reposo en la aproximación escalar. Cuando $a = 1/2$, el frente bistable aislado queda estacionario; para valores mayores se invierte su dirección según la orientación elegida. En el medio excitable completo esta observación orienta la discusión de propagación o fracaso de la excitación, aunque la frontera exacta entre ambos regímenes dependerá de más parámetros que $a$.

146.7. Recuperación lenta, refractariedad y distancia entre impulsos

Durante una fase activa prolongada podemos aproximar $u \approx 1$. Con $\beta = 0$, la ecuación de recuperación local se reduce a $\dot{v} \approx \varepsilon(1 - \gamma v)$. Si $v(0) = v_{0}$, su solución es

La aproximación describe el aumento gradual de la variable de recuperación mientras la activación permanece alta. En algún momento, la respuesta rápida deja de sostener ese estado y $u$ cae. Después, cuando $u \approx 0$, la recuperación satisface $\dot{v} \approx - \varepsilon\gamma v$ y disminuye como $v(t) = v_{post}e^{- \varepsilon\gamma t}$. La escala temporal $\tau_{rec} = 1/(\varepsilon\gamma)$ caracteriza esa relajación aproximada, aunque la duración de la refractariedad funcional depende además del umbral necesario para responder de nuevo.

Con $\varepsilon = 0,01$ y $\gamma = 1$, la escala anterior es de cien unidades de tiempo adimensional. Ese valor no equivale a un período refractario de exactamente cien unidades, porque una región puede recuperar su capacidad de conducción antes de que $v$ vuelva prácticamente a cero; esa capacidad, además, puede depender de la intensidad y la duración del siguiente estímulo. Para medir una curva de recuperación, podemos aplicar dos pulsos separados por un intervalo variable y registrar la amplitud o la posibilidad de propagación del segundo. El experimento produce una relación funcional, más informativa que un único corte temporal.

La recuperación convierte la dirección espacial en una propiedad emergente. Si el primer pulso avanza hacia la derecha, la región que deja a su izquierda permanece temporalmente menos excitable. El frente puede reclutar medio fresco a la derecha, pero encuentra material en recuperación detrás. En una geometría regular esta asimetría favorece la propagación ordenada. La ley de difusión sigue siendo simétrica; la dirección se origina en la historia espacialmente distribuida de la variable lenta.

En un tren de pulsos periódicos que se desplazan con velocidad aproximadamente constante $c$, la separación espacial entre pulsos puede aproximarse por $\Lambda \approx cT$, donde $T$ es el período con que una posición recibe activaciones. Esta relación cinemática es sencilla, pero el período permitido depende también del estado de recuperación. Si se estimula demasiado deprisa, el pulso siguiente encuentra un valor elevado de $v$ y puede reducir su velocidad, cambiar de amplitud o bloquearse. La velocidad puede depender, además, del intervalo diastólico o de recuperación, una dependencia conocida en fisiología cardíaca como restitución de conducción.

La longitud refractaria, como escala orientativa, puede estimarse mediante $L_{ref} \sim c\tau_{ref}$ si el frente se desplaza aproximadamente a velocidad constante y el tiempo de recuperación se define de manera consistente. La relación ayuda a evaluar si una trayectoria espacial es suficientemente larga para alojar un pulso sin que su frente encuentre la cola refractaria. La existencia de una onda recurrente depende además de la geometría, de las interacciones entre pulsos y de la heterogeneidad del tejido, que siguen siendo determinantes.

146.8. Curvatura, anisotropía y bloqueo de conducción

En una franja recta, un frente excitable tiene una vecindad parecida a lo largo de toda su extensión. En una superficie, la geometría puede cambiar la facilidad de propagación. Un frente muy convexo hacia la región que intenta activar distribuye su influencia entre un área receptora cada vez mayor. La pérdida de eficacia debida a la curvatura se estudia mediante relaciones de tipo eikonal. En el límite de frente delgado, curvatura pequeña y recuperación casi uniforme junto a la interfaz, una aproximación común es

donde $V_{n}$ es la velocidad normal, $c_{0}$ la velocidad de un frente plano en condiciones equivalentes y $\kappa$ la curvatura del frente con la convención geométrica adecuada. La fórmula es local y aproximada. Si el radio de curvatura se vuelve comparable al espesor del frente, o la recuperación cambia significativamente a lo largo de la interfaz, es necesario resolver el sistema espacial completo.

El efecto de la curvatura es especialmente marcado cuando una excitación atraviesa un estrechamiento y emerge hacia una región ancha. El frente pasa de reclutar un número pequeño de puntos vecinos a influir sobre una región receptora mucho mayor. Puede propagarse, ralentizarse o bloquearse según la capacidad excitadora y el estado del medio. El fenómeno tiene una interpretación dinámica de fuente y sumidero; su presencia no implica que la difusión deje de ser simétrica. La propagación depende de la relación entre el estímulo efectivo que llega a cada posición y la dinámica local que ese estímulo desencadena.

En muchos tejidos la transmisión no es isótropa. Las propiedades de acoplamiento pueden ser mayores a lo largo de una dirección estructural que en la dirección transversal. La ecuación de activación se generaliza entonces, con un tensor positivo definido $\mathbf{D}$, a

mientras la recuperación conserva su ecuación propia. La velocidad y la forma del frente dependen de su orientación respecto de los ejes principales de $\mathbf{D}$. En un medio homogéneo con tensor constante, una transformación lineal de coordenadas puede ayudar a comparar la propagación anisótropa con el caso isótropo. Si el tensor varía de una región a otra, las interfaces pueden deformarse y aparecer efectos de conducción que una única difusividad escalar oculta.

Las interacciones entre pulsos muestran otra diferencia respecto de las ondas lineales. En la ecuación de ondas ideal, dos perturbaciones pequeñas pueden superponerse y atravesarse. En un medio excitable, cada pulso altera el estado por donde pasa. Dos frentes que se aproximan pueden encontrarse con regiones recién activadas, perder la capacidad de seguir avanzando y aniquilarse en la colisión. En otros contextos, obstáculos o regiones no excitables pueden fragmentar una onda y cambiar el patrón espacial de propagación. Esos fenómenos necesitan la no linealidad y la memoria local de recuperación.

Cuando los frentes se curvan, se bloquean parcialmente o rodean obstáculos, aparecen condiciones para estructuras espaciales más complejas, incluidas ondas giratorias persistentes. Su descripción exige estudiar la organización de los extremos de frente, los defectos y la topología de las interfaces. Ese problema se abordará por separado, porque no se deduce de la mera presencia de una onda viajera unidimensional.

146.9. Verificación numérica de un pulso y sus errores habituales

Para simular un medio excitable conviene comenzar con un dominio unidimensional y un equilibrio basal numéricamente comprobado. Puede estimularse una región pequeña mediante $I(x,t)$ durante un intervalo conocido, seguido de un período sin forzamiento. El primer resultado que interesa es cualitativo. ¿La perturbación desaparece, produce una excursión local que luego se extingue o genera un pulso que sigue viajando después de retirar el estímulo? La respuesta depende de la amplitud, la duración y la extensión espacial del estímulo. Registrar solo la actividad máxima impediría distinguir estas posibilidades.

Una discretización explícita de la activación usa un factor $D\Delta t/\Delta x^{2}$ que multiplica la segunda diferencia espacial, además de la reacción local y el estímulo. En una dimensión la parte difusiva pura exige habitualmente $D\Delta t/\Delta x^{2} \leq 1/2$ para la estabilidad del método de Euler explícito. Sin embargo, la ecuación acoplada contiene una dinámica rápida no lineal y una recuperación lenta; el paso temporal debe resolver ambas y controlar el error de integración. Un método implícito o semimplícito puede mejorar la estabilidad difusiva sin garantizar automáticamente precisión de velocidad o de umbral.

La anchura de un frente puede ser bastante menor que la distancia recorrida por el pulso. Una malla incapaz de resolver esa anchura introduce errores en el flujo que llega a las posiciones vecinas y puede cambiar artificialmente la velocidad o provocar bloqueo. El problema es especialmente grave cerca del umbral de propagación, cuando una pequeña variación de las condiciones numéricas decide si la onda se sostiene. Por ello debe repetirse el cálculo reduciendo $\Delta x$, ajustando $\Delta t$ de manera controlada y comparando los resultados para varios niveles de resolución.

Conviene medir por separado la posición del frente delantero, el máximo de activación y el borde posterior. La distancia entre el frente delantero y el borde posterior aproxima la anchura espacial del pulso; el intervalo entre sus pasos por un mismo punto da su duración local. Los resultados deberían hacerse aproximadamente estacionarios si existe un pulso viajero estable en el dominio considerado. La variable de recuperación requiere mediciones propias (amplitud máxima, demora respecto de $u$ y tiempo necesario para aproximarse al estado basal).

Las condiciones de frontera pueden alterar el diagnóstico. Una frontera periódica devuelve un pulso al lugar del que salió y puede producir un tren de excitaciones que no existiría en un tejido abierto. Una condición reflectante para $u$ restringe el flujo, pero no necesariamente representa correctamente un extremo físico. Si el propósito es medir una velocidad en un medio extenso, es preferible colocar los bordes suficientemente lejos de la región observada y comprobar que la velocidad no cambia al ampliar el dominio. El registro del tipo de frontera forma parte de la reproducibilidad del experimento.

Finalmente, una auditoría útil combina pruebas con distinto grado de complejidad. Primero comprueba la estabilidad del reposo mediante el jacobiano. Después reproduce el frente escalar rápido con la velocidad exacta del apartado 146.6. Solo entonces estudia pulsos del sistema completo y su sensibilidad a estímulos y parámetros. Esta secuencia separa errores del código de fenómenos auténticos del modelo. Un gráfico vistoso de propagación no sustituye una comprobación de convergencia, particularmente cuando se estudian umbrales, bloqueos y velocidades pequeñas.

146.10. Qué puede inferirse de una onda de excitación observada

Supongamos que un tejido cultivado responde a un estímulo localizado con una banda de actividad que se aleja del punto de origen. Detrás de la banda, la señal desciende y permanece durante un tiempo menos sensible a una segunda estimulación. Esta observación es compatible con una dinámica excitable, pero deja abiertas muchas hipótesis mecanicistas. Podría haber variables eléctricas, procesos bioquímicos, acoplamientos mecánicos o varios mecanismos actuando de manera conjunta. El modelo de dos variables captura una clase de comportamientos; identificar sus variables con cantidades físicas específicas requiere evidencia adicional.

Una investigación cuantitativa comenzaría por describir las escalas espaciales y temporales de la señal medida. Hay que distinguir la resolución del sistema de observación de la duración real del evento, porque una onda rápida puede parecer amplia cuando se muestrea lentamente. Después mediríamos la velocidad del frente, la anchura de la zona activa, el tiempo de recuperación y la dependencia de la respuesta respecto de la amplitud y duración del estímulo. Repetiríamos las observaciones ante estímulos separados por intervalos variables para determinar si existe refractariedad funcional.

Estas mediciones permiten confrontar modelos. Un sistema escalar con una transición permanente puede reproducir un borde que avanza, pero no la recuperación observada detrás. Un modelo excitable puede describir activación seguida de retorno, pero aun así fallar en la relación entre velocidad e intervalo de recuperación. Un modelo con advección puede producir una banda móvil por transporte del material, mientras que uno de reacción–difusión predice activación secuencial del medio local. La discriminación requiere observables que distingan esas posibilidades, más allá de una semejanza visual entre simulación y experimento.

La historia de la neurofisiología ilustra la utilidad de separar niveles de descripción. Hodgkin y Huxley desarrollaron una representación cuantitativa de corrientes de membrana con variables de compuerta dependientes del voltaje. FitzHugh exploró cómo una dinámica con pocas variables podía reproducir de manera cualitativa estados excitables y oscilatorios; Nagumo y colaboradores aplicaron además esos principios a una línea activa de transmisión de pulsos. Las formulaciones no son intercambiables. El modelo simplificado ayuda a entender la organización dinámica; uno biofísicamente detallado permite contrastar otros aspectos de las corrientes y propiedades de conducción.

En aplicaciones biomédicas debe extremarse la distinción entre una predicción matemática y una conclusión clínica. Que un modelo produzca bloqueo o una onda recurrente bajo ciertos parámetros no demuestra que esos parámetros representen un tejido individual ni permite, sin otras pruebas, establecer un diagnóstico o una intervención. La relación entre modelo y fenómeno exige calibración, validación con observables independientes y un análisis explícito de sensibilidad. Hay que ser especialmente prudente cuando varias combinaciones de parámetros producen señales eléctricas superficialmente parecidas.

La noción central de este capítulo es la recuperación espacialmente distribuida. La señal viaja porque una región activa recluta a otra; su duración es finita porque la propia activación induce procesos que terminan inhibiéndola. Las regiones que acaba de recorrer conservan temporalmente un valor elevado de $v$, que cambia las posibilidades de propagación de ondas posteriores y permite que las interacciones entre frentes sean fuertemente no lineales. Cuando las interfaces se curvan y encuentran obstáculos, la recuperación también condiciona la supervivencia de sus extremos y la aparición de estructuras giratorias. El capítulo C147 estudiará esa organización mediante defectos, espirales y persistencia espacial.

En tejido cardíaco, el paso desde una ecuación excitable genérica hasta una explicación de la propagación eléctrica exige incorporar propiedades celulares y estructurales adicionales. La revisión de Alonso, Bär y Echebarria (2016) describe esa variedad de mecanismos y evita interpretar cada espiral simulada como equivalente a un circuito fisiopatológico demostrado.

Referencias

• Hodgkin, A. L., y Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology, 117(4), 500–544. https://doi.org/10.1113/jphysiol.1952.sp004764
• FitzHugh, R. (1961). Impulses and physiological states in theoretical models of nerve membrane. Biophysical Journal, 1(6), 445–466. https://doi.org/10.1016/S0006-3495(61)86902-6
• Nagumo, J., Arimoto, S., y Yoshizawa, S. (1962). An active pulse transmission line simulating nerve axon. Proceedings of the IRE, 50(10), 2061–2070. https://doi.org/10.1109/JRPROC.1962.288235
• Barkley, D. (1991). A model for fast computer simulation of waves in excitable media. Physica D: Nonlinear Phenomena, 49(1–2), 61–70. https://doi.org/10.1016/0167-2789(91)90194-E
• Keener, J., y Sneyd, J. (2009). Mathematical Physiology (2.ª ed., 2 vols.). Springer.
• Winfree, A. T. (1987). When Time Breaks Down: The Three-Dimensional Dynamics of Electrochemical Waves and Cardiac Arrhythmias. Princeton University Press.
• Alonso, S., Bär, M., y Echebarria, B. (2016). Nonlinear physics of electrical wave propagation in the heart: A review. Reports on Progress in Physics, 79(9), 096601. https://doi.org/10.1088/0034-4885/79/9/096601