🌀 Réaction de Belousov-Zhabotinsky

L'Oregonator est un modèle EDO simplifié à trois variables de la réaction de BZ, proposé par Field et Noyes en 1974, distillant le mécanisme complexe FKN (Field-Körös-Noyes) en trois espèces clés : HBrO₂ (X), le catalyseur oxydé (Y), et l'acide bromomalonique (Z). Les équations dX/dt = s(qY + X(1−X) − XY)/ε, dY/dt = (−qY − XY + fZ)/ε', dZ/dt = X − Z décrivent les boucles de rétroaction positive (X) et négative (Y, Z). Malgré sa simplicité, l'Oregonator reproduit quantitativement la période d'oscillation et la vitesse d'onde des expériences BZ réelles.

Dans l'automate cellulaire de Greenberg-Hastings utilisé ici, le seuil détermine combien de voisins de Moore excités (sur 8) une cellule au repos doit avoir pour se déclencher. Un seuil bas (1 voisin) signifie qu'un seul voisin excité déclenche le tir — les ondes se propagent rapidement, tendent à remplir tout le domaine, et les spirales tournent vite avec une courte longueur d'onde. Un seuil élevé (5-6) signifie que seules les régions densément excitées déclenchent leurs voisins, produisant des anneaux concentriques largement espacés, lents et qui s'éteignent facilement. Les milieux BZ réels ont un seuil d'excitabilité effectif ajustable en modifiant la concentration de bromate ou d'acide.

Oui — les ondes en spirale excitables apparaissent largement en biologie. La fibrillation cardiaque (atriale comme ventriculaire) implique des ondes en spirale rotatives d'activation électrique qui perturbent le battement cardiaque périodique normal ; la défibrillation met fin à ces spirales en dépolarisant simultanément tout le tissu. L'amibe sociale Dictyostelium discoideum utilise des ondes en spirale d'AMPc pour coordonner l'agrégation en période de famine. La dépression corticale envahissante — la base neurologique de l'aura migraineuse — se propage sous forme d'une onde de dépolarisation lente à environ 3 mm/min. La réaction de BZ a été utilisée comme modèle physique pour étudier tous ces phénomènes.

La période réfractaire détermine combien de temps une cellule doit attendre après son déclenchement avant de pouvoir être réexcitée. Une période réfractaire plus longue (plus de pas) augmente la longueur d'onde spatiale des spirales — les cellules derrière le front d'onde mettent plus de temps à récupérer, si bien que la prochaine crête d'onde doit parcourir une plus grande distance avant de rencontrer une région rétablie. Des périodes réfractaires très longues peuvent supprimer entièrement la formation de spirales, ne permettant qu'une émission d'ondes cibles isolées à partir de régions de stimulation. Ce paramètre correspond directement à la période réfractaire relative en électrophysiologie cardiaque, qui régit la vulnérabilité aux arythmies.

Oui — les systèmes de réaction-diffusion BZ ont été proposés et démontrés comme substrats de calcul non conventionnels. Adamatzky et ses collègues ont utilisé des ondes BZ dans des milieux gélifiés pour résoudre des problèmes de labyrinthe (l'onde trouve le chemin le plus court), implémenter des portes logiques (collision d'ondes → ET/OU), et même effectuer du traitement d'images. Le parallélisme inhérent au système — chaque emplacement spatial évolue simultanément — offre des avantages potentiels pour certains calculs massivement parallèles. Plus concrètement, le calcul par réaction-diffusion éclaire la compréhension du calcul morphogénétique en biologie du développement.

Le taux de déclenchement spontané contrôle la fréquence à laquelle les cellules au repos s'allument sans avoir suffisamment de voisins excités — représentant le bruit thermique ou des sites de nucléation hétérogènes dans un milieu BZ réel. À taux spontané nul, la grille finira par atteindre un point fixe silencieux à moins d'être réamorcée manuellement. Des taux spontanés faibles maintiennent le motif en nucléant occasionnellement de nouvelles spirales lorsque celles existantes se dégradent. Des taux spontanés élevés créent un déclenchement dense et chaotique qui détruit les spirales organisées et produit des ondes cibles désordonnées — analogue à une expérience BZ menée à haute température où les fluctuations thermiques dominent.

L'article d'Alan Turing de 1952 « The Chemical Basis of Morphogenesis » prédisait que deux espèces chimiques diffusantes — un activateur et un inhibiteur — l'inhibiteur diffusant plus vite que l'activateur, pouvaient spontanément briser la symétrie spatiale et générer des motifs périodiques (points, rayures, labyrinthes). La réaction de BZ démontre cette classe d'instabilité de Turing en trois dimensions : les ondes en spirale et les ondes cibles sont des motifs de Turing spatio-temporels. Le même cadre mathématique explique les motifs du pelage animal (taches de léopard, rayures de zèbre), l'espacement des doigts dans les membres en développement, et l'espacement des photorécepteurs dans l'œil.

L'annihilation d'onde lors d'une collision est une caractéristique des milieux excitables et découle de la période réfractaire. Lorsque deux fronts d'onde se rencontrent, les cellules au point de collision viennent d'être excitées par les deux fronts qui s'approchaient ; au lieu de générer une nouvelle excitation, ces cellules entrent dans l'état réfractaire. Derrière chaque front se trouve une zone réfractaire qui ne peut pas propager l'excitation dans la direction inverse, si bien qu'aucune des deux ondes ne continue au-delà de la collision — elles s'annihilent. C'est fondamentalement différent des milieux oscillatoires où deux ondes se traversent l'une l'autre, et des ondes linéaires où se produit une interférence constructive/destructive.