Réaction chimique oscillante - [Apprendre en ligne]

Chimie

Réaction chimique oscillante

Réaction de Belousov Zhabotinski. Oregonateur

Modèle de la réaction chimique oscillante de Belousov Zhabotinski.

Article mis en ligne le 7 février 2017

par bernard.vuilleumier

Belousov découvrit en 1958 que l’oxydation de l’acide citrique par le bromate de potassium catalysé par des ions céreux-cériques pouvait donner lieu à une réaction chimique oscillante et Zhabotinski poursuivit cette étude [1]. La réaction de Belousov Zhabotinski s’étudie d’habitude à 25 °C et consiste en un mélange comportant du bromate de potassium, de l’acide malonique ou bromomalonique et du sulfate cérique dissous dans l’acide sulfurique. Elle peut donner lieu à une grande variété de phénomènes allant des oscillations d’une période de quelques dizaines de secondes à l’apparition de phénomènes ondulatoires. Modélisation de la réaction.

Idealized Belousov-Zhabotinsky Reaction from the Wolfram Demonstrations Project by Luca Zammataro

Le modèle

Notations

Nous proposons un modèle permettant de simuler l’évolution des concentrations des trois substances principales, concentrations que nous désignerons de la manière suivante :

X =[HBrO₂] Y =[Br^-] Z =2[Ce⁴⁺]

Écrivons encore :

A=B=[BrO₃^-]

Schéma réactionnel

Le mécanisme de Noyes [2] peut alors s’exprimer par les étapes suivantes :

A+Y$\rightarrow$X X+Y$\rightarrow$P B+X$\rightarrow$2X+Z 2X$\rightarrow$Q Z$\rightarrow$Y

Ce schéma réactionnel, appelé « Oregonateur », comporte un mécanisme de catalyse croisée : Y produit X, X produit Z qui en retour produit Y.

Équations différentielles correspondant au schéma réactionnel

En examinant où chaque substance apparaît ou disparaît dans le schéma réactionnel, en admettant une variation proportionnelle aux concentrations en présence et en posant A=B, on établit les équations différentielles :

$\frac{dX}{dt}=AY+AX-XY-X^2$ $\frac{dY}{dt}=Z-AY-XY$ $\frac{dZ}{dt}=AX-Z$

Diagramme Stella traduisant les équations différentielles

Les équations différentielles se traduisent aisément en un diagramme Stella qui met ici en évidence le mécanisme de catalyse croisée :

**Oregonateur**

Le diagramme *Stella* met bien en évidence le mécanisme de catalyse croisée : Y produit X, X produit Z qui produit Y.

Constantes cinétiques et paramètre d’ajustement

Pour simuler une réaction oscillante, nous introduisons les constantes cinétiques $k_i$ et un paramètre d’ajustement f dans les équations [3] :

$\frac{dX}{dt}=k_1AY-k_2XY+k_5AX-2k_6X^2$ $\frac{dY}{dt}=-k_1AY-k_2XY+fk_7Z$ $\frac{dZ}{dt}=k_5AX-k_7Z$

En résolvant ces équations et en reportant les concentrations X, Y et Z en fonction du temps, on fait alors clairement apparaître des oscillations :

**Oscillation des concentrations en fonction du temps**

Avec un choix judicieux des concentrations initiales X₀, Y₀, Z₀, des constantes cinétiques k_i et du paramètre d’ajustement f, on fait apparaître des oscillations.

Concentrations initiales, constantes cinétiques et paramètre f

INIT X = 0.25
flux_X = k1*A*Y+k5*A*X-k2*X*Y-2*k6*X^2

INIT Y = 2.1
flux_Y = f*k7*Z-k1*A*Y-k2*X*Y

INIT Z = 0.83
fluy_Z = k5*A*X-k7*Z

A = 1
f = 1
k1 = 0.05
k2 = 0.8
k5 = 1.58
k6 = 0.2
k7 = 0.4

Utiliser un simulateur de la réaction

Conclusion

Ce type de réaction illustre le fait que le comportement thermodynamique d’un système peut être très différent loin de l’équilibre et à l’opposé de celui prévu par le théorème de production minimale d’entropie établi par Prigogine [4]. Le non-équilibre peut être une source d’ordre.

Voir aussi (from Wolfram Demonstrations Project)
– Complex Dynamics in Chemical Self-Replication
– Idealized Belousov Zhabotinsky Reaction

Wolfram Demonstrations Project : mode d’emploi