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

147.1. Una onda que regresa al lugar de donde partió

Imaginemos una lámina de material excitable en la que se propaga un pulso. Cada punto se activa brevemente, atraviesa un período de recuperación y vuelve a quedar disponible para una nueva excitación. Si el impulso llega a un borde recto y el medio es homogéneo, es razonable esperar una onda que avance aproximadamente como una línea. La situación cambia cuando una parte del frente se interrumpe. El extremo libre queda rodeado por regiones con estados de recuperación diferentes y puede comenzar a curvarse. En determinadas condiciones, ese extremo gira alrededor de una zona que permanece sin activarse por completo. La onda ya no necesita una fuente que emita impulsos desde fuera, porque la propia geometría permite que la excitación vuelva a encontrar medio recuperado.

Arthur T. Winfree describió en 1972 ondas espirales de actividad química en reactivos de la familia de Belousov–Zhabotinsky. Su resultado mostró que una organización espacial periódica podía mantenerse mediante la propagación local de la actividad. No hacía falta imaginar una estructura sólida que trazara la espiral ni un reloj central que coordinara todos los puntos. La forma visible surgía de la interacción entre excitación, difusión y recuperación. Una espiral fotografiada no revela su mecanismo, pues podría ser estática, producto del transporte o un efecto óptico.

En una espiral de excitación, distintos puntos del medio pueden recorrer ciclos similares de activación con desfases espaciales. Cerca de su centro se encuentra una región singular desde el punto de vista de la fase. Lejos de ella, los frentes suelen aproximarse a trenes de ondas que avanzan hacia el exterior. La frecuencia, la separación de frentes y la trayectoria del centro dependen de la recuperación, el transporte y las fronteras. La espiral puede ser estable, desplazarse lentamente, oscilar de manera compleja o fragmentarse.

La persistencia necesita una explicación precisa. En un medio puramente difusivo, una perturbación localizada termina alisándose. En un medio excitable, una perturbación que supere un umbral puede regenerarse conforme avanza. La energía o los reactivos que sostienen cada activación provienen de las condiciones del sistema; la geometría reorganiza su uso. Una espiral persistente no es un movimiento perpetuo ni viola la disipación. En una preparación química abierta, por ejemplo, el mantenimiento prolongado puede requerir renovación de reactivos; en una preparación cerrada, el agotamiento termina alterando o suprimiendo las ondas.

Este capítulo analiza por qué ciertos frentes rotos generan espirales, cómo se identifican sus centros y qué significa hablar de un defecto espacial. También distingue entre una espiral que rota de manera regular, una espiral cuyo centro deriva y una colección de defectos que nacen y desaparecen en una dinámica irregular. Son fenómenos que pueden coexistir sin ser equivalentes. Hay que distinguirlos porque una misma imagen instantánea puede corresponder a evoluciones radicalmente distintas.

147.2. Del pulso unidimensional al frente curvo

El punto de partida es una versión bidimensional del modelo de activación y recuperación introducido en el capítulo anterior. Podemos representar la activación por u(x,y,t)u(x,y,t) y la variable lenta de recuperación por v(x,y,t)v(x,y,t). Una familia amplia de sistemas adopta la forma

∂tu=Du∇2u+f(u,v),∂tv=Dv∇2v+εg(u,v),0<ε≪1.\partial_{t}u = D_{u}\nabla^{2}u + f(u,v),\quad\quad\partial_{t}v = D_{v}\nabla^{2}v + \varepsilon g(u,v),\quad 0 < \varepsilon \ll 1.

Aquí, DuD_{u} y DvD_{v} son coeficientes de difusión no negativos, mientras que ε\varepsilon expresa una separación entre escalas de tiempo. La existencia de un régimen excitable depende de las funciones ff y gg, de sus parámetros y del estado de referencia. Es posible construir medios excitables sin esta separación extrema; la ecuación aísla una arquitectura dinámica frecuente. En una dimensión, la solución relevante puede ser un pulso que avanza con velocidad aproximadamente constante. En dos dimensiones debemos considerar además la orientación y curvatura locales de su frente.

Si la mitad superior de una onda plana se bloquea por una barrera temporal y la mitad inferior continúa propagándose, aparece un extremo libre. El avance normal del frente y la recuperación espacialmente desigual pueden hacer que ese extremo se arquee. Detrás del frente queda una región transitoriamente refractaria que el extremo no puede atravesar, de modo que la frontera entre regiones recuperadas y refractarias condiciona su trayectoria. Con todo, la formación de una espiral requiere que el medio, las dimensiones y la perturbación permitan el giro sostenido.

La rotación tiene una asimetría elemental, ya que bajo las mismas leyes locales una espiral puede girar en sentido horario y otra en sentido antihorario. En un medio homogéneo e isótropo, los dos sentidos suelen estar relacionados por reflexión. La historia de la preparación y la perturbación que inicia la actividad pueden seleccionar uno. Si el medio tiene anisotropías o una estructura quiral, la equivalencia por reflexión puede romperse. No debemos atribuir a una ley química una preferencia de giro que quizá proviene de la geometría experimental.

Conviene distinguir la velocidad del frente de la velocidad de desplazamiento de su extremo. Una onda casi plana tiene una normal bien definida y puede desplazarse con una velocidad cercana a la del pulso unidimensional. Cerca del extremo libre, la curvatura y la recuperación modifican el avance; incluso es posible que un punto concreto del frente avance lentamente mientras la figura visible rota a una frecuencia apreciable. El centro geométrico de una espiral no es, necesariamente, un punto material que gira con el patrón.

Si se observa la preparación durante un intervalo corto, una espiral puede parecer una curva estática. Para confirmar que transmite excitación hace falta seguir la variable de activación en el tiempo, identificar el frente y comprobar la recuperación entre dos pasos consecutivos. Una secuencia temporal de alta resolución permite distinguir una onda activa de la simple rotación aparente de una textura advectada. Volveremos a necesitar esta precaución al abordar el transporte dirigido y la inestabilidad espaciotemporal.

147.3. Fase espacial y singularidad en el centro de giro

La noción de fase ayuda a describir un medio donde muchas regiones atraviesan ciclos semejantes con distintos tiempos de activación. Para un oscilador estable, una fase θ\theta aumenta aproximadamente 2π2\pi por cada ciclo. En un medio excitable que permanece en reposo entre estímulos, no existe necesariamente una fase autónoma igual de bien definida. Aun así, durante una actividad periódica sostenida puede construirse una fase operacional a partir de dos observables, por ejemplo activación y recuperación, si sus trayectorias locales rodean de manera consistente un punto de referencia en el plano de estados.

Una definición ilustrativa consiste en introducir z=(u−u∗)+i(v−v∗)z = \left( u - u_{*} \right) + i\left( v - v_{*} \right) y escribir

θ(x,y,t)=atan2(v−v∗, u−u∗).\theta(x,y,t) = atan2\left( v - v_{*},\, u - u_{*} \right).

El par (u∗,v∗)\left( u_{*},v_{*} \right) debe seleccionarse de manera justificada. Esta fase deja de estar definida allí donde ambos componentes se anulan simultáneamente y puede ser poco informativa cuando una trayectoria no completa una vuelta alrededor del punto escogido. Tampoco tiene que coincidir con una fase asintótica construida desde la dinámica del oscilador. Su utilidad depende de la estabilidad del procedimiento frente al ruido, la discretización y la transformación de las variables.

Supongamos que trazamos una curva cerrada Γ\Gamma que rodea el centro de una espiral sin pasar por puntos donde la fase sea singular. Al avanzar alrededor de la curva, puede ocurrir que la fase aumente una vuelta completa. El número de enrollamiento, para una fase continua a lo largo de la curva, se escribe

q=12π∮Γ∇θ⋅dℓ.q = \frac{1}{2\pi}\oint_{\Gamma}^{}\nabla\theta \cdot d\mathbf{\ell}.

En condiciones apropiadas, qq es un entero. Un giro positivo de fase puede corresponder a q=+1q = + 1 y el giro inverso a q=−1q = - 1, según la convención elegida. El resultado depende de cómo cambia el estado alrededor de la singularidad; la espiral puede estar lejos de dibujar un círculo perfecto. Una perturbación suave que no atraviese una singularidad ni modifique los límites de la región conserva este índice. Por eso una clasificación topológica puede sobrevivir a deformaciones que alteran considerablemente la forma visible.

El ejemplo matemático más simple es el campo complejo z(x,y)=x+iyz(x,y) = x + iy. Sobre una circunferencia centrada en el origen, su argumento recorre 2π2\pi; en el centro, la fase no está definida. Si usamos z(x,y)=x−iyz(x,y) = x - iy, el recorrido cambia de signo. Estos dos campos carecen de dinámica temporal y no modelan ninguna reacción excitable; únicamente ilustran por qué un defecto puede localizarse mediante la fase y por qué su orientación tiene un signo.

En datos experimentales, la fase suele calcularse a partir de registros temporalmente filtrados o de reconstrucciones del espacio de estados. El centro detectado depende entonces de la calidad del registro y del método de interpolación. Cuando la señal se vuelve débil, una fase numérica forzada puede producir singularidades espurias. Un análisis responsable comprueba que los defectos persisten al variar de manera razonable el umbral de amplitud, la ventana de filtrado y la resolución espacial.

Un contorno antihorario rodea el cero de un campo complejo; la fase desenvuelta aumenta una vuelta para x más iy y disminuye una para x menos iy.
Figura 147.1. En el contorno circular, z = x + iy acumula un cambio de fase de +2π y z = x − iy uno de −2π. Sus índices son +1 y −1. En ambos campos el origen tiene amplitud cero y fase indefinida. El sentido del contorno se fija como antihorario. Las construcciones son estáticas: el índice no representa por sí solo una frecuencia ni una dinámica de espiral excitable. Ampliar figura ↗

147.4. Creación, desplazamiento y aniquilación de defectos topológicos

El término defecto puede resultar engañoso, ya que en este contexto no indica necesariamente un fallo del sistema. Un defecto es aquí un lugar donde se interrumpe la posibilidad de prolongar suavemente una descripción ordenada, como una fase bien definida o una orientación espacial continua. Las dislocaciones de patrones de franjas, los núcleos de vórtices y las singularidades de fase pertenecen a familias emparentadas, aunque sus variables físicas y sus grupos de simetría difieren. Compartir vocabulario topológico no implica tener el mismo mecanismo de generación.

En una región plana con fase continua salvo en defectos aislados, el índice medido en una curva cerrada suficientemente alejada de ellos equivale a la suma de las cargas topológicas encerradas. En consecuencia, los defectos pueden aparecer o desaparecer en pares de signos opuestos sin cambiar la carga neta interior, siempre que la frontera permanezca libre de singularidades y se mantengan las hipótesis que permiten definir la fase. Si un defecto atraviesa la frontera, la carga observada dentro de la región sí puede cambiar. En un medio excitable la construcción de fase y la propia aparición de singularidades exigen especial cuidado, por lo que las reglas no deben aplicarse sin comprobar el dominio en que tienen sentido.

Imaginemos una lámina con dos espirales de sentidos de giro contrarios. Si sus núcleos se acercan y la dinámica permite su encuentro, las singularidades pueden aniquilarse. La desaparición topológica no implica que toda actividad de la lámina cese de inmediato, ya que pueden seguir propagándose frentes residuales hasta alcanzar fronteras o zonas refractarias. En dirección inversa, una onda fragmentada puede dejar dos extremos que evolucionen como defectos de signo opuesto. El patrón resultante depende del balance algebraico de cargas y también de la historia de recuperación y del acoplamiento.

La conservación de la carga neta tampoco fija las posiciones de los defectos. Un defecto puede moverse por acción de una pequeña heterogeneidad, por interacción con otro núcleo o por perturbaciones temporales. En ciertos regímenes, muchas singularidades se crean y destruyen de forma irregular. El observador puede reconocer una textura aproximadamente repetitiva, pero una comparación entre instantes revela que sus unidades estructurales cambian continuamente. Este tipo de dinámica permite pasar de una espiral aislada a lo que en algunos modelos se denomina turbulencia de defectos.

En patrones de franjas hay otra distinción útil. Una dislocación puede identificarse cuando una franja termina dentro de una región sin llegar a una frontera. Su carga describe el desajuste en el número de períodos espaciales al rodear el defecto. Una espiral excitable, en cambio, se identifica por la organización temporal de la fase alrededor de un centro rotatorio. Ambas pueden aparecer como interrupciones de un patrón espacial, aunque no hay una traducción directa y universal entre sus cargas.

La topología da un criterio que resiste las deformaciones, aunque es incompleto para explicar la dinámica. Dos sistemas pueden presentar defectos con idéntico número de enrollamiento y comportamientos muy distintos. El índice no informa de la frecuencia de rotación, de las tasas de recuperación ni de la fuerza que desplaza un núcleo. Para llegar a esas magnitudes necesitamos regresar a las leyes locales y estudiar cómo interactúan con la geometría.

Dos ceros de un campo complejo tienen índices opuestos; los contornos individuales acumulan más y menos una vuelta y el contorno exterior acumula cero.
Figura 147.2. El producto z = (x + 1 + iy)(x − 1 − iy) tiene ceros en (−1, 0) y (1, 0), con índices +1 y −1. La fase a lo largo de un contorno exterior acumula cero vueltas netas, aunque encierre ambos defectos. Se verifica la aditividad con contornos que no atraviesan los ceros. Es una construcción matemática estática: no simula creación, desplazamiento ni aniquilación. Ampliar figura ↗

147.5. Curvatura del frente y nacimiento de un núcleo

Un frente curvo no avanza exactamente igual que uno plano. En algunos modelos de reacción–difusión con una escala de espesor pequeña frente al radio de curvatura, la velocidad normal puede aproximarse por una relación eikonal de la forma

cn≃c0−Defκ,c_{n} \simeq c_{0} - D_{ef}\kappa,

con c0c_{0} como velocidad de un frente plano, κ\kappa la curvatura con una convención de signo explícita y DefD_{ef} un coeficiente efectivo que depende del modelo. Para una curva convexa hacia la región todavía no activada, tomamos aquí κ>0\kappa > 0. Esta fórmula no es una ley universal para cualquier frente y su coeficiente efectivo no puede sustituirse sin justificación por la difusión de una especie. Su validez requiere hipótesis asintóticas y puede romperse cuando la curvatura es grande, la recuperación cambia bruscamente o el medio es anisótropo.

La relación contiene una intuición geométrica. Un frente muy curvado puede propagarse más lentamente porque la difusión redistribuye su contribución entre una región espacial más amplia o porque el perfil local se aparta del pulso plano de referencia. Si aproximamos la punta por un arco circular de radio RR, tenemos κ≈1/R\kappa \approx 1/R. En la aproximación eikonal, el avance normal tendería a anularse cerca de R≈Def/c0R \approx D_{ef}/c_{0}. Esta expresión no predice con exactitud el radio de núcleo de una espiral real, que depende también de la recuperación, la cinética local, la propagación tangencial y la estructura completa de la solución giratoria.

Un cálculo sencillo muestra el alcance del argumento y sus límites. En unidades adimensionalizadas, supongamos c0=2c_{0} = 2, Def=1D_{ef} = 1 y una curvatura de 0,40,4. La fórmula estimaría cn≈1,6c_{n} \approx 1,6. Para una curvatura de 11, daría cn≈1c_{n} \approx 1. Al extrapolarla hasta curvaturas mucho mayores acabaríamos obteniendo velocidades negativas. Esa extrapolación carece de fundamento si el régimen de curvatura débil ya dejó de cumplirse. El ejemplo solo ilustra la sensibilidad cualitativa del avance.

En una espiral sostenida, el extremo gira alrededor de una zona donde las propiedades de la onda son distintas de las del campo lejano. Las capas de activación y recuperación forman un circuito espacial que permite la rotación. Si el tiempo de recuperación aumenta, un frente que llega demasiado pronto encuentra medio refractario y puede bloquearse. Si el tiempo disminuye, quizá cambien la frecuencia y el tamaño del núcleo. La frecuencia de rotación no puede deducirse solamente de c0c_{0}, porque resulta del problema completo, incluidas la condición de frontera y la selección no lineal de la solución.

La relación entre curvatura y propagación también explica por qué una onda puede fracasar al intentar atravesar un estrechamiento. Un impulso que sale de un canal angosto a una región amplia debe expandirse con un radio inicial relativamente pequeño. Esa expansión puede exigir una curvatura incompatible con su capacidad local de propagación. La posibilidad de bloqueo depende de la geometría y de la excitabilidad, además del umbral específico del modelo. La misma preparación podría transmitir una onda en una dirección y bloquearla en la otra, con lo que se crearían condiciones propicias para la reentrada aunque no exista necesariamente una espiral estable.

147.6. Ritmo, longitud de onda y movimiento del centro

Cuando una espiral gira de manera aproximadamente estacionaria, su frecuencia de rotación determina el intervalo con que una región suficientemente alejada del núcleo recibe nuevos frentes. Si ese período es TT, la frecuencia es f=1/Tf = 1/T y la velocidad angular media vale ω=2π/T\omega = 2\pi/T. Un observador fijo ve pasar ondas sucesivas, mientras que el patrón en conjunto puede representarse aproximadamente en un marco que rota con velocidad angular ω\omega. Esa representación deja de ser exacta si el núcleo deriva, el medio presenta heterogeneidades significativas o la espiral cambia de forma durante cada vuelta.

En la región donde los frentes son casi planos y avanzan con velocidad cc, la distancia entre dos crestas consecutivas puede estimarse como λ≈cT\lambda \approx cT. Si, en unidades de un modelo, c=3c = 3 unidades de longitud por unidad de tiempo y T=8T = 8, la separación estimada es de unas 24 unidades de longitud. El cálculo exige que velocidad y período se refieran al mismo régimen espacial y que el frente observado sea realmente parte de un tren periódico. Cerca del núcleo, la velocidad normal cambia y la simple relación puede dejar de describir la distancia local entre brazos.

La recuperación impone otra condición. Un punto que se activó debe disponer de tiempo suficiente para recobrar excitabilidad antes de recibir el frente siguiente. En una preparación con recuperación lenta, aumentar la frecuencia puede reducir la amplitud de los impulsos, retrasar la conducción o impedir que el siguiente frente se propague. Por esta razón la longitud de onda de la excitación, entendida como extensión espacial del impulso más su estela refractaria relevante, no coincide necesariamente con la distancia geométrica entre dos máximos consecutivos de una señal. Ambas pueden usarse como medidas, pero requieren definiciones distintas.

Algunas espirales mantienen un núcleo prácticamente fijo; otras tienen meandro, es decir, una trayectoria periódica o cuasiperiódica del centro de rotación en torno a una posición media. Una espiral que meandra puede generar una trayectoria del extremo semejante a una flor de varios lóbulos sin que su actividad local se haya vuelto caótica. El centro puede presentar también una deriva secular superpuesta al meandro, de modo que, al observar la trayectoria durante mucho tiempo, la posición promedio cambia.

La frecuencia de rotación, la velocidad de deriva y las oscilaciones del núcleo son observables diferentes. Ajustar un único círculo a una secuencia de imágenes puede ocultar esa diferencia. Para medirlos conviene seguir la fase y el centro a través de varias rotaciones completas, informar el método de localización y cuantificar su sensibilidad al intervalo de muestreo. Un período aparente especialmente regular obtenido de cinco fotogramas no autoriza a concluir que existe una solución rotatoria estable.

147.7. Heterogeneidad, anisotropía y fronteras que organizan la onda

Ningún medio experimental es perfectamente homogéneo. Cambian la conductividad, la disponibilidad de reactivos, la geometría de los obstáculos y la velocidad local de recuperación. Una espiral que rota en un sistema ideal puede quedar anclada a una región heterogénea o desplazarse lentamente por su influencia. En ciertos casos un obstáculo no excitable tiene un contorno alrededor del cual una onda puede circular; en otros, el obstáculo bloquea la reentrada al reducir el espacio disponible o al alterar la trayectoria de recuperación. Cuál de los dos efectos domina se decide examinando cada sistema, porque los obstáculos no siempre favorecen las espirales.

La anisotropía exige particular atención. Si la difusión es más intensa en una dirección que en otra, la difusividad es un tensor D\mathbf{D} y el término espacial puede escribirse como ∇⋅(D∇u)\nabla \cdot \left( \mathbf{D}\nabla u \right). Cuando D\mathbf{D} es simétrico y definido positivo, el frente suele propagarse con velocidades diferentes según su orientación. Una espiral puede entonces presentar brazos elípticos o un núcleo deformado. Comparar esas formas con una espiral circular en coordenadas físicas puede sugerir diferencias dinámicas que en parte provienen de la métrica espacial del transporte.

En cuanto a las fronteras, una pared impermeable para una sustancia suele modelarse mediante un flujo normal nulo, mientras que una frontera absorbente exige otra relación. En un dominio finito las ondas emitidas por el núcleo alcanzan las paredes, pueden desaparecer o interactuar con ellas de acuerdo con el modelo. La dinámica del núcleo puede cambiar mucho antes de que el frente lo toque, porque una perturbación localizada se transmite a través del campo espacial. Es inadecuado estudiar una espiral en una caja pequeña y extrapolar sin control su frecuencia a un medio infinito.

Hay una teoría perturbativa útil para analizar la deriva de espirales casi estacionarias. Las llamadas funciones de respuesta describen cómo perturbaciones de diferentes posiciones se proyectan sobre modos asociados al movimiento y a la fase del patrón. Su concentración cerca del núcleo, que se observa en varios modelos, ayuda a explicar por qué pequeñas heterogeneidades locales pueden afectar con eficacia la trayectoria de la espiral. El trabajo de Biktasheva y colaboradores ha mostrado comparaciones cuantitativas entre estas predicciones y simulaciones. La aproximación, sin embargo, depende de que exista un patrón de referencia y de que las perturbaciones sean suficientemente pequeñas.

En un experimento podemos distinguir anclaje, deriva y meandro siguiendo la trayectoria del núcleo antes y después de modificar una heterogeneidad definida. También hay que mantener controlados los estímulos externos y las condiciones de frontera. Una espiral que parece quedar fija cerca de un obstáculo puede estar simplemente confinada por la geometría; demostrar anclaje dinámico exige mostrar una relación reproducible entre el obstáculo y la respuesta del sistema.

147.8. De la espiral bidimensional a las ondas de enrollamiento tridimensionales

En dos dimensiones el núcleo se identifica como un punto singular de fase. Si agregamos una dimensión espacial, ese conjunto puede extenderse y formar una línea de filamento. Las superficies de activación se enrollan alrededor del filamento y generan una onda de tipo scroll. El filamento puede ser recto, curvo o cerrado en un anillo. Una sección bidimensional de ese campo puede mostrar una espiral, aunque la configuración completa sea tridimensional y cambie con el tiempo de maneras que el corte no revela.

La longitud y la curvatura del filamento pueden evolucionar; en determinados regímenes un anillo se contrae y en otros puede expandirse. Se suele describir parte de este comportamiento mediante el concepto de tensión efectiva del filamento, pero el signo y la magnitud de esa tensión dependen de los parámetros y del modelo, de modo que la contracción no es una regla general. También pueden aparecer deformaciones helicoidales, interacciones con fronteras y procesos de reconexión que cambian la topología de las estructuras activas.

La geometría tridimensional introduce una dificultad de observación. Un mapa superficial puede mostrar uno o varios centros de giro que no coinciden con todos los filamentos presentes en el volumen. Dos cortes cercanos pueden contener singularidades en posiciones distintas porque atraviesan un filamento curvo. Una imagen superficial que parece mostrar el nacimiento de un par de defectos podría reflejar el movimiento de un filamento con respecto al plano registrado. Para decidir entre esas explicaciones se necesita información volumétrica o una inferencia sustentada por un modelo tridimensional verificable.

En materiales excitables biológicos, una onda rotatoria bidimensional es útil para estudiar posibilidades de reentrada, pero un tejido real puede incorporar capas, fibras, anisotropía y discontinuidades. Una onda rotatoria observada en la superficie no caracteriza automáticamente toda la actividad subyacente. Los estudios experimentales de espirales en músculo cardíaco aislado muestran que la dinámica de reentrada es físicamente plausible. De ellos no se sigue que todas las arritmias respondan a una única organización geométrica ni que puedan ignorarse otros mecanismos clínicos.

Las soluciones de filamento también muestran la diferencia entre persistencia y estabilidad. Una onda tridimensional puede seguir generando activación durante muchas vueltas mientras su filamento sufre deformaciones crecientes. El patrón persiste temporalmente, pero quizá se aproxime a un cambio de régimen o a una ruptura. Si quisiéramos comprobar estabilidad, estudiaríamos la respuesta a perturbaciones definidas y su evolución a largo plazo, además de registrar si la actividad continúa.

147.9. Qué tienen en común una reacción química, el corazón y un patrón de convección

Los experimentos con la reacción de Belousov–Zhabotinsky dieron ejemplos visuales de frentes químicos que se reorganizan en espirales. En esos sistemas, la cinética reactiva y el transporte sostienen ondas espaciales de actividad. La literatura sobre medios excitables desarrolló modelos genéricos que preservan algunos rasgos de esa propagación sin reproducir todos los pasos químicos. La comparación es útil porque separa las propiedades que dependen de la arquitectura dinámica de las propias de un reactivo particular.

En electrofisiología cardíaca se han observado ondas espirales de activación en preparaciones tisulares. Los trabajos de los grupos de Davidenko y de Pertsov incluyen demostraciones experimentales y análisis computacionales de la actividad reentrante. Los potenciales de membrana, las corrientes iónicas, la anatomía y la conducción anisótropa tienen una especificidad que no puede comprimirse en dos variables sin pérdida de información. Con un modelo reducido de excitabilidad pueden comprenderse principios de reentrada; la predicción clínica exige modelos y mediciones a la escala correspondiente.

Hay otro uso de la palabra «espiral» en la física de patrones. En sistemas convectivos, como algunos experimentos de Rayleigh–Bénard, pueden aparecer franjas que se enrollan alrededor de defectos y cambian de forma de modo irregular. Morris y colaboradores describieron experimentalmente un régimen denominado caos de defectos espirales. Esas espirales provienen de una dinámica hidrodinámica distinta y no necesitan corresponder a impulsos de activación con un período refractario. La semejanza geométrica permite estudiar preguntas compartidas sobre simetría, defectos y correlación espacial, pero no autoriza a usar indistintamente los mecanismos de un tejido excitable y un fluido convectivo.

Un patrón de Turing es un tercer contraste. En el escenario clásico de inestabilidad estacionaria, las perturbaciones espaciales crecen y pueden saturarse en una organización de manchas o franjas que permanece fija en el tiempo. En una espiral excitable hay movimiento continuo de activación alrededor de un núcleo. Sin embargo, los sistemas reales pueden combinar oscilaciones temporales, patrones estacionarios, ondas y defectos. Los nombres designan mecanismos o regímenes identificables, aunque compatibles entre sí.

De esta comparación se desprende una regla útil de lectura científica. Ante una imagen de bandas curvas, conviene preguntar qué variable está representada, qué ocurre en cada punto a través del tiempo y si la forma se mantiene cuando se modifica una perturbación. Los criterios morfológicos orientan la investigación; su valor explicativo aumenta cuando se relacionan con observables dinámicos y ecuaciones justificadas.

147.10. Cómo estudiar una espiral sin confundir visualización con explicación

Supongamos que un investigador observa una estructura que parece girar sobre una lámina experimental. El primer paso es precisar la magnitud registrada (concentración, intensidad fluorescente, potencial eléctrico o desplazamiento material). A continuación debe comprobar que la sucesión de imágenes permite seguir la activación local y que el intervalo de muestreo resuelve suficientemente el período de rotación. Si se toman imágenes demasiado separadas, un patrón rápido puede parecer inmóvil, girar en sentido contrario o mostrar una frecuencia inferior a la real por aliasing temporal.

Después puede construir mapas de fase. Si usa una transformación de dos variables, debe indicar cómo seleccionó el origen del plano de estados y dónde la fase deja de estar bien definida. Para identificar defectos, el procedimiento debería calcular el cambio de fase a lo largo de pequeños contornos cerrados, con un desenvolvimiento coherente que evite saltos artificiales de 2π2\pi. Conviene repetir la detección al cambiar la resolución y la suavización de la señal. Los defectos auténticos deben mantener una localización compatible dentro de la incertidumbre espacial esperable, aunque los datos ruidosos nunca permiten una comprobación topológica perfecta.

Una vez identificado el centro, las magnitudes de interés son la frecuencia, la velocidad del frente en regiones alejadas, el radio aproximado de la zona no activada, la trayectoria del núcleo y el tiempo de recuperación. Cada una requiere un método de medida y, cuando corresponde, un intervalo de incertidumbre. Si una espiral se desplaza, se deben separar la deriva lenta y el meandro periódico. Si aparecen varios núcleos, interesa seguir las parejas que se crean o aniquilan y registrar los encuentros con fronteras. Un conteo instantáneo de defectos pierde la mayor parte de esa información temporal.

La confrontación con un modelo puede avanzar por niveles. Un sistema reducido de reacción–difusión permite preguntar si la excitabilidad y el acoplamiento local bastan para producir la organización observada. Un modelo más detallado puede investigar cómo cambian los resultados al introducir heterogeneidad, anisotropía o geometría tridimensional. En ambos casos es necesario comparar observables definidos antes de ajustar parámetros; reproducir una imagen de espiral es una prueba débil si su frecuencia, su respuesta a una perturbación y su estabilidad resultan incompatibles con los datos.

Hay que reconocer también lo que no se ha demostrado. El hallazgo de una singularidad de fase es una descripción geométrica robusta cuando la fase está bien definida; no identifica de manera unívoca el mecanismo biológico o químico. La permanencia de un patrón durante varias vueltas tampoco prueba estabilidad indefinida. La desaparición tras cambiar una frontera puede deberse a una alteración real del mecanismo o a un artefacto de la preparación. El análisis gana precisión cuando las conclusiones se restringen a las escalas y observables efectivamente estudiados.

Las espirales muestran que, en los sistemas espacialmente extendidos, una ley local sencilla puede producir una organización que mantiene su identidad mientras cambia continuamente el estado de sus componentes. La fase y la topología describen aspectos de esa identidad; la propagación y la recuperación explican su mantenimiento. Para comprender regímenes con muchos defectos móviles, advectados o mutuamente inestables, tendremos que pasar ahora a la dinámica espaciotemporal y a los procesos que transportan las perturbaciones a través del dominio.

Referencias