La reacción de Belousov–Zhabotinsky (BZ) es una familia de reacciones químicas oscilantes descubiertas por Borís Belousov en 1951 y ampliadas por Anatol Zhabotinsky en la década de 1960 — uno de los ejemplos más espectaculares de autoorganización en desequilibrio en química. En una mezcla BZ típica — ácido malónico, bromato, ácido sulfúrico y un catalizador metálico — la disolución oscila espontáneamente entre estados oxidado y reducido, formando ondas químicas concéntricas o espirales visibles a simple vista. Estos patrones surgen porque el sistema actúa como un medio excitable: un estado de reposo, un umbral de excitación, un periodo refractario — el mismo comportamiento que muestran el músculo cardíaco y las fibras nerviosas.
Este simulador utiliza un autómata celular de tres estados (el modelo de Greenberg–Hastings) para aproximar la dinámica excitable de la reacción BZ en una cuadrícula bidimensional. Haz clic en el lienzo para plantar «chispas» asimétricas que inician centros espirales. Ajusta el umbral de excitación, el número de pasos refractarios y la frecuencia de encendido espontáneo para explorar distintos modos de onda. Cuatro esquemas de color (Fuego, Océano, Espectro, Mono) representan los estados de las células.
¿Por qué oscila la reacción BZ?
La reacción BZ oscila porque contiene un bucle de retroalimentación negativa con retardo: el bromato oxida el catalizador metálico, produciendo ácido bromoso (HBrO₂), que se amplifica a sí mismo de forma autocatalítica — una retroalimentación positiva. Cuando la concentración del catalizador alcanza su pico, se acumula ácido bromomalónico, que inhibe la autocatálisis y devuelve el sistema a su estado reducido. Esta autocatálisis inhibida por el bromuro crea un ciclo límite en el espacio de concentraciones químicas. El periodo característico en una mezcla BZ bien agitada suele ser de 30 a 120 segundos.
¿Qué es un medio excitable y en qué se diferencia de un oscilador?
Un medio excitable tiene un estado de reposo estable que requiere una perturbación umbral para dispararse — a diferencia de un verdadero oscilador, que oscila continuamente sin entrada externa. Una vez disparada, una célula excitable produce una respuesta de gran amplitud y luego entra en un periodo refractario durante el cual no puede dispararse de nuevo. El músculo cardíaco y los axones nerviosos (ecuaciones de Hodgkin–Huxley) son los medios excitables biológicos clásicos. La reacción BZ es excitable en condiciones subumbrales y oscilatoria por encima de un límite de parámetros — una bifurcación de Hopf.
¿Por qué se forman espirales en lugar de simples anillos expansivos?
Las espirales surgen de condiciones iniciales asimétricas — en este simulador, una semilla de encendido en forma de medio disco crea la asimetría necesaria. Una onda circular en expansión se aniquila cuando sus frentes de onda colisionan; pero si parte del frente queda bloqueada o llega más tarde, se forma una punta libre que se curva hacia adentro, iniciando una espiral autosostenida. Las ondas espirales son atractores: perturbarlas tiende a restaurar la espiral, lo que las hace muy robustas.
El Oregonator es un modelo simplificado de EDO con tres variables de la reacción BZ, propuesto por Field y Noyes en 1974. Reduce el complejo mecanismo FKN a tres especies clave: HBrO₂ (X), el catalizador oxidado (Y) y el ácido bromomalónico (Z). Pese a su simplicidad, este modelo reproduce cuantitativamente el periodo de oscilación y la velocidad de onda de los experimentos BZ reales.
En el autómata celular de Greenberg–Hastings usado aquí, el umbral determina cuántas vecinas excitadas necesita una célula en reposo para dispararse. Un umbral bajo (1 vecina) hace que las ondas se propaguen rápidamente, y las espirales giran rápido con longitud de onda corta. Un umbral alto (5–6) produce anillos concéntricos lentos y muy espaciados que se extinguen fácilmente. En una mezcla BZ real, el umbral de excitabilidad efectivo se regula cambiando la concentración de bromato o de ácido.
Sí. La fibrilación cardíaca (tanto auricular como ventricular) implica ondas espirales rotatorias de activación eléctrica que alteran el latido periódico normal; la desfibrilación termina estas espirales. La ameba social Dictyostelium discoideum usa ondas espirales de cAMP para coordinar la agregación durante la inanición. La depresión cortical propagada — la base neurológica del aura de la migraña — se propaga como una onda de despolarización lenta a unos 3 mm/min.
El periodo refractario determina cuánto tiempo debe esperar una célula tras dispararse antes de poder volver a excitarse. Un mayor número de pasos aumenta la longitud de onda espacial de las espirales — las células detrás del frente de onda tardan más en recuperarse, así que la siguiente cresta debe recorrer más distancia. Un periodo refractario muy largo puede suprimir por completo la formación de espirales, dejando solo ondas aisladas emitidas desde regiones marcapasos.
Sí. Los sistemas de reacción-difusión BZ se han propuesto y demostrado como sustratos de computación no convencional. Adamatzky y colaboradores han usado ondas BZ en medios de gel para resolver problemas de laberinto, implementar puertas lógicas (colisión de ondas → AND/OR) e incluso procesar imágenes. El paralelismo inherente al sistema — cada punto espacial evoluciona simultáneamente — ofrece ventajas potenciales para ciertos cálculos masivamente paralelos.
La frecuencia de encendido espontáneo controla con qué frecuencia las células en reposo se encienden sin suficientes vecinas excitadas — representando el ruido térmico o los sitios de nucleación heterogénea en una mezcla BZ real. Con frecuencia cero, la cuadrícula acabará alcanzando un punto fijo silencioso salvo que se siembren chispas manualmente. Una frecuencia alta crea un encendido caótico y denso que destruye las espirales organizadas y produce ondas concéntricas desordenadas — análogo a un experimento BZ a alta temperatura donde dominan las fluctuaciones térmicas.
El artículo de Alan Turing de 1952 «La base química de la morfogénesis» predijo que dos especies químicas difusoras — un activador y un inhibidor —, con el inhibidor difundiendo más rápido que el activador, podían romper espontáneamente la simetría espacial y generar patrones periódicos. La reacción BZ demuestra esta clase de inestabilidad de Turing en tres dimensiones. El mismo marco matemático explica los patrones del pelaje animal, el espaciado de los folículos pilosos en la piel de los mamíferos y la distancia entre dedos en las extremidades de los vertebrados.
La aniquilación de ondas al colisionar es un sello distintivo de los medios excitables y surge del periodo refractario. Cuando dos frentes se encuentran, las células en el punto de colisión ya han sido excitadas por ambos frentes que se aproximan; en lugar de generar una nueva excitación, esas células entran en estado refractario. Detrás de cada frente hay una zona refractaria que no puede propagar la excitación en dirección inversa, por lo que ninguna onda continúa más allá del punto de colisión. Esto es fundamentalmente distinto de los medios oscilatorios, donde las ondas se atraviesan entre sí.