Oscilaciones químicas y modelo Brusselator

COMPLEXA · Ciencias de la Complejidad y Teoría del Caos
Módulo 18. Organización fuera del equilibrio y redes de reacción

176.1. Un modelo oscilador químico ideal

El Brusselator describe cómo una red de reacciones alimentada desde el exterior mantiene variaciones periódicas de concentración. Con una cinética sencilla permite estudiar estabilidad, realimentación, bifurcaciones y ciclos límite; también permite distinguir las oscilaciones amortiguadas alrededor de un estado estacionario de las oscilaciones autosostenidas bajo condiciones externas constantes. Esta distinción organiza el capítulo y delimita lo que puede deducirse con cada herramienta matemática.

Las variables dinámicas $x(t)$ e $y(t)$ son concentraciones adimensionales de los intermediarios $X$ e $Y$; $t$ es tiempo adimensional. Los parámetros positivos $a$ y $b$ son niveles efectivos de alimentación de las especies externas $A$ y $B$. Mantenerlos constantes supone reservorios de composición regulada o de variación despreciable en la escala temporal considerada.

El modelo es determinista y espacialmente homogéneo; cada concentración tiene un único valor en el reactor y quedan fuera las fluctuaciones moleculares y los gradientes espaciales. La reducción es apropiada para estudiar el mecanismo temporal ideal si la mezcla es suficientemente rápida y las poblaciones moleculares admiten una representación continua; su validez experimental debe comprobarse independientemente. La adimensionalización absorbe escalas de concentración, tiempo y constantes cinéticas, por lo que los coeficientes unitarios de las ecuaciones posteriores reflejan una elección de unidades y de parametrización. Recuperar tiempos y concentraciones físicos exige especificar las escalas originales (un período adimensional solo se expresa en segundos con la constante temporal correspondiente).

Dos variables bastan para producir una órbita periódica en un sistema autónomo continuo. Una ecuación escalar autónoma con unicidad de soluciones no puede hacerlo, porque una trayectoria no vuelve repetidamente al mismo valor con velocidades de signos distintos. En el plano, la interacción entre acumulación y consumo permite recorridos cerrados; la no linealidad autocatalítica hace que la intensidad de la respuesta dependa del estado instantáneo.

La construcción de Prigogine y Lefever (1968) es un marco ideal para estudiar inestabilidades disipativas, que requieren un sostenimiento externo. Una concentración estacionaria puede coexistir con flujos químicos permanentes y una concentración oscilante con alimentación invariable; estacionariedad y equilibrio termodinámico son propiedades distintas. Se investigan aquí estados dinámicos de un sistema abierto, sin desarrollar una contabilidad termodinámica completa de sus reservorios. El capítulo deriva los balances, identifica el único punto fijo positivo y determina cuándo pierde estabilidad; después examina el mecanismo oscilatorio, las restricciones sobre promedios temporales y un protocolo de integración reproducible. Se distinguen las conclusiones locales de las globales y las predicciones exactas de las aproximaciones válidas solo cerca del umbral oscilatorio.

176.2. Esquema de reacciones y balances

El Brusselator combina alimentación, conversión, autocatálisis y eliminación. Las concentraciones de $A$ y $B$ se mantienen fijas mediante un quimiostato; $D$ y $E$ son productos que no retroalimentan la cinética. Solo $X$ e $Y$ requieren balances dinámicos. El esquema ideal consta de cuatro reacciones irreversibles

La primera reacción produce una unidad de $X$; la segunda consume $X$ y genera $Y$ y $D$; la tercera gana netamente una unidad de $X$ y pierde una de $Y$; la cuarta elimina una unidad de $X$. En las unidades adimensionales, las velocidades efectivas son $v_1=a$, $v_2=bx$, $v_3=x^2y$ y $v_4=x$, todas no negativas en el cuadrante físico. Los parámetros $a$ y $b$ miden, respectivamente, la alimentación constante de $A$ y el efecto de la concentración fija de $B$.

El paso autocatalítico describe una amplificación efectiva de $X$ en presencia de $Y$, sin presuponer como mecanismo elemental universal una colisión simultánea de tres moléculas. Redes microscópicas con intermediarios pueden producir dependencias efectivas semejantes bajo hipótesis de reducción apropiadas.

Si la prima indica la derivada respecto del tiempo adimensional, multiplicar cada velocidad por su cambio estequiométrico y sumar las contribuciones da

La conversión entre intermediarios aparece con signos opuestos, por lo que se cancela al sumar los balances de $X$ e $Y$. La alimentación y la eliminación modifican su suma; la dinámica resultante acopla esos flujos con una velocidad autocatalítica proporcional a $x^2y$, que aumenta con $x$ y depende de la disponibilidad de $y$.

176.3. Derivación de las ecuaciones y positividad

Al sustituir las cuatro velocidades en los balances se obtiene el sistema diferencial del Brusselator. La alimentación da un término constante positivo; la conversión y la eliminación generan dos pérdidas lineales de $X$; el paso autocatalítico añade la contribución cúbica. Para $Y$, la producción lineal y el consumo autocatalítico son las únicas contribuciones retenidas. Con $a>0$ y $b>0$, las ecuaciones son

El campo vectorial es polinómico y, por tanto, localmente Lipschitz respecto de las concentraciones. Cada condición inicial finita determina una solución local única. Así, dos soluciones distintas no pasan por el mismo estado en el mismo instante. Sin embargo, la regularidad local de un campo polinómico no garantiza que todas sus soluciones permanezcan finitas para cualquier tiempo.

El dominio físico es el cuadrante cerrado de concentraciones no negativas. Sobre su frontera vertical, donde $x=0$, la primera ecuación da $x'=a>0$, de modo que el campo apunta hacia el interior sea cual sea el valor de $y$. Sobre la frontera horizontal, donde $y=0$, la segunda ecuación da $y'=bx\geq0$. Ninguna frontera permite, pues, salir hacia concentraciones negativas mientras la solución existe.

Estas condiciones, junto con la unicidad local, establecen la invariancia positiva del cuadrante. Una trayectoria iniciada con concentraciones no negativas permanece en él durante su intervalo de existencia. Si ambas concentraciones empiezan estrictamente positivas, ninguna alcanza una frontera ni la atraviesa hacia valores negativos. En el origen, la alimentación pone inmediatamente en marcha la producción de $X$; después, la conversión alimenta también la concentración de $Y$.

La invariancia positiva debe distinguirse de la acotación. Para examinar esta última conviene introducir la concentración total retenida $s=x+y$. La cancelación de los términos de conversión conduce a $s'=a-x$ y, como $x\geq0$, a una cota superior lineal en el tiempo

La cota controla las soluciones físicas en cualquier intervalo temporal finito e impide así una explosión en tiempo finito. Si el extremo del intervalo de existencia fuera finito, las concentraciones seguirían dentro de un conjunto compacto y la solución podría prolongarse. Por consiguiente, las soluciones físicas existen para todo tiempo futuro. El argumento usa la positividad y la estructura particular de los balances; los coeficientes polinómicos solos no habrían bastado.

La estimación anterior crece con $t$ y no da una cota uniforme cuando el tiempo tiende a infinito; establecer una región absorbente o un confinamiento asintótico exige un razonamiento adicional. Existencia global, permanencia en el dominio físico y acotación uniforme son propiedades diferentes, aunque puedan intervenir juntas en una demostración global.

176.4. Nulclinas y estado estacionario

Las nulclinas son los conjuntos donde se anula una componente del campo vectorial; permiten leer la dirección instantánea sin integrar las ecuaciones. La condición $y'=0$ equilibra la producción y el consumo autocatalítico de $Y$; la condición $x'=0$, la alimentación, las pérdidas lineales y la producción de $X$. En el interior físico, $x>0$, ambas condiciones definen las curvas

La nulclina de $Y$ es una hipérbola positiva decreciente. Por debajo de ella, $bx-x^2y>0$ y la concentración de $Y$ aumenta; por encima domina el consumo y disminuye. Aunque $x=0$ también anula $y'$, allí la alimentación mantiene $x'>0$, así que esa componente no contiene estados estacionarios.

La nulclina de $X$ corta el eje horizontal en $x=a/(b+1)$ y su rama físicamente relevante satisface $x\geq a/(b+1)$; para valores menores, $y$ es negativa. Como el coeficiente de $y$ en $x'$ es $x^2>0$, por encima de la curva $x'>0$ y por debajo $x'<0$; así quedan fijadas las componentes horizontales del movimiento.

En un estado estacionario, sumar los balances da $0=a-x$, luego $x_*=a$; al sustituir en $y'=0$, resulta $y_*=b/a$. La positividad de los parámetros garantiza coordenadas finitas y positivas, de modo que la intersección queda en el cuadrante físico, en el punto

La deducción prueba también que el estado estacionario físico es único, porque todo candidato debe tener primero $x=a$ y luego $y=b/a$. Unicidad y estabilidad son independientes; la ubicación de la intersección no indica si las trayectorias vecinas se acercan o se alejan, cuestión que requiere estudiar el campo alrededor del punto fijo.

Las nulclinas no son, en general, trayectorias; sobre la de $X$ el movimiento es vertical si $y'\neq0$ y sobre la de $Y$ es horizontal si $x'\neq0$. Las trayectorias suelen atravesarlas; confundirlas con órbitas distorsiona la lectura de los máximos y mínimos. Los extremos de $x(t)$ e $y(t)$ se sitúan en sus respectivas nulclinas, salvo degeneraciones, y pueden ocurrir en instantes distintos. El desfase se debe al recorrido por regiones de crecimiento y decrecimiento.

176.5. Linealización y bifurcación de Hopf

La estabilidad local del estado estacionario se obtiene expandiendo el campo alrededor de su posición. Sean $u=x-a$ y $v=y-b/a$ las desviaciones de las concentraciones. Los términos de primer orden definen la evolución lineal de perturbaciones pequeñas; los términos cuadráticos y cúbicos describen las correcciones no lineales. Al evaluar las derivadas parciales en el punto fijo se obtiene la matriz jacobiana

En un sistema bidimensional, la traza y el determinante bastan para localizar los autovalores respecto del eje imaginario. Aquí la traza es $\tau=b-1-a^2$ y el determinante es $\Delta=a^2$, siempre positivo. Los autovalores son las raíces de $\lambda^2-\tau\lambda+\Delta=0$. El determinante positivo excluye un punto silla e impide también que un autovalor se anule mientras $a>0$.

Cuando $\tau<0$, ambos autovalores tienen parte real negativa y el punto fijo es localmente asintóticamente estable; según el signo del discriminante $\tau^2-4a^2$, la aproximación puede ser nodal u oscilatoria. Si $\tau>0$, las partes reales son positivas y el punto fijo es inestable. Esta clasificación describe perturbaciones suficientemente pequeñas; el comportamiento de trayectorias alejadas requiere información que la aproximación lineal descarta.

Al aumentar $b$ con $a$ fijo, la traza cambia de signo en un valor crítico. Allí el polinomio característico se reduce a $\lambda^2+a^2=0$, de modo que los autovalores forman un par imaginario puro. La frecuencia angular lineal es $a$ en las unidades de tiempo elegidas. El umbral y el período asociado a esa frecuencia se expresan mediante

Cerca del umbral, los autovalores siguen siendo complejos y su parte real vale $\tau/2$. Su derivada respecto de $b$ es, por tanto, $1/2$, que no se anula. El cruce transversal del eje imaginario satisface así una de las condiciones de la bifurcación de Hopf. La condición restante concierne a los términos no lineales y decide si emerge una rama de órbitas periódicas y cuál es su estabilidad.

La igualdad $\tau=0$ identifica el candidato a Hopf, pero unos autovalores imaginarios puros no implican que el sistema no lineal sea un centro. En el propio umbral, la aproximación lineal deja indeterminada la evolución radial de amplitudes pequeñas. Los términos de orden superior pueden producir una deriva lenta que se acumula a tiempos largos, aunque el giro domine la apariencia inmediata.

El período $2\pi/a$ es el límite del período de las oscilaciones nacientes al aproximarse a una Hopf no degenerada. A distancia finita del umbral, la frecuencia puede recibir correcciones dependientes de la amplitud y del parámetro. Queda pendiente la relación entre crecimiento lineal y saturación no lineal, que explica por qué la pérdida de estabilidad puede generar oscilaciones sostenidas sin que la linealización determine su tamaño.

176.6. Casos numéricos didácticos y alcance del análisis lineal

La elección didáctica $a=1$ sitúa el umbral en $b_{\mathrm H}=2$ y permite comparar tres regímenes mediante cálculos algebraicos explícitos. El punto fijo es entonces $(1,b)$, la traza vale $b-2$ y el determinante es uno. Los ejemplos solo evalúan el jacobiano y su polinomio característico, sin integrar ni medir amplitudes.

Para $b=1{,}5$, los autovalores son $-0{,}25\pm0{,}9682458i$. El estado estacionario es un foco estable, en el que las perturbaciones pequeñas giran con una envolvente que, en la aproximación lineal, decrece exponencialmente. La parte imaginaria es la frecuencia angular de esa relajación local. En este régimen, una señal con varios máximos sucesivos puede ser un transitorio amortiguado, de modo que las oscilaciones visibles no prueban que haya un autooscilador.

Para $b=2$, los autovalores son $\pm i$. La dinámica lineal conserva la amplitud de sus oscilaciones, pero en el campo completo la deriva radial de las perturbaciones pequeñas la deciden los términos no lineales. Como la tasa lineal de amortiguación se anula, cerca del umbral cabe esperar transitorios lentos, que deben considerarse en simulaciones cortas.

Para $b=2{,}2$, los autovalores son $0{,}1\pm0{,}9949874i$. El punto fijo es un foco inestable y las perturbaciones pequeñas crecen mientras giran a su alrededor. El crecimiento exponencial solo describe la salida inicial del entorno lineal; al aumentar la desviación, las contribuciones no lineales ganan peso y la solución deja de seguir esa extrapolación. La tasa positiva tampoco da el valor de los máximos y mínimos finales de las concentraciones.

En el Brusselator estándar, la bifurcación de Hopf es supercrítica para $a>0$. Demostrarlo exige un cálculo no lineal, expresable mediante el primer coeficiente de Lyapunov; la propiedad puede verificarse con el cálculo de Ferčec y colaboradores (2021) y se sitúa en el marco clásico de Prigogine y Lefever (1968). Con $\mu=b-(1+a^2)$, una coordenada radial local adecuada admite una forma normal cuyo coeficiente cúbico estabilizador es positivo

En la convención de Lyapunov empleada por Ferčec y colaboradores, su ecuación (42), especializada a las tasas $(k_1,k_2,k_3,k_4)=(a,1+a^2,1,1)$, da

Numerador y denominador son positivos antes del signo menos para todo $a>0$. Esta cuenta da el control no lineal que faltaba en el criterio de traza nula. El coeficiente $g_1$ y el coeficiente radial $c(a)$ de la forma normal usan normalizaciones diferentes (el valor de $c(a)$ depende de cómo se normalice $r$); la información que comparten aquí es el signo estabilizador.

Para $\mu>0$ suficientemente pequeño aparece una órbita periódica estable, con radio local proporcional a $\mu^{1/2}$. Esa proporcionalidad es asintótica y se refiere a coordenadas transformadas. Traducirla a amplitudes de $x$ o de $y$ exige conocer la transformación y los coeficientes correspondientes, o integrar el sistema con precisión controlada.

La supercriticidad establece la existencia y la estabilidad de una rama próxima al umbral. Una bifurcación de Hopf genérica no asegura la unicidad global del ciclo, la atracción desde cualquier condición inicial ni el comportamiento para parámetros arbitrariamente alejados. Los tres casos ilustran una transición local; la cuenca, el período lejos de la rama naciente y la forma completa de la órbita requieren un análisis separado.

176.7. Mecanismo temporal y retrato de fase

La oscilación puede interpretarse como un proceso recurrente de autoamplificación, agotamiento y recuperación. La alimentación introduce $X$, mientras la conversión dependiente de $b$ transfiere parte de esa concentración a $Y$. Cuando ambas especies alcanzan una combinación adecuada, la velocidad $x^2y$ hace que el aumento de $x$ refuerce su propia producción, hasta que el consumo de $Y$ y las pérdidas limitan ese crecimiento.

La realimentación positiva depende de una reserva dinámica. El término autocatalítico aumenta $x$ y reduce $y$ con igual contribución absoluta en sus balances. Por ello, una fase de crecimiento rápido de $X$ puede consumir el intermediario necesario para sostenerla. Cuando $y$ disminuye suficientemente, la amplificación pierde intensidad. La alimentación y la conversión permiten después reconstruir las condiciones de una nueva fase de crecimiento, con un desfase entre las dos concentraciones.

Esta descripción identifica funciones dentro del mecanismo, pero no impone tres etapas discretas ni velocidades claramente separadas para todos los parámetros. Cerca de Hopf, la oscilación tiene amplitud pequeña y una forma aproximadamente sinusoidal. Más lejos del umbral, la no linealidad puede deformar la señal y redistribuir el tiempo de permanencia a lo largo del recorrido. La existencia de tramos rápidos y lentos debe verificarse para el régimen concreto considerado.

Una trayectoria temporal muestra $x(t)$ e $y(t)$ frente al mismo reloj y permite reconocer desfases, tiempos entre máximos y variaciones de amplitud durante el transitorio. El retrato de fase muestra, en cambio, los pares $(x(t),y(t))$. Una oscilación periódica corresponde a un recorrido cerrado; la velocidad con que se recorre esa curva sigue siendo información temporal adicional.

Las nulclinas ayudan a interpretar los extremos de las señales sin dibujar una trayectoria supuesta. Al cruzar $x'=0$, la concentración $x$ cambia de crecimiento a decrecimiento o viceversa, según el sentido del cruce. Lo mismo sucede con $y$ sobre $y'=0$. Las regiones delimitadas por ambas curvas codifican los signos simultáneos de las velocidades; esa lectura relaciona la forma de las señales temporales con la geometría del campo vectorial.

El sentido local de giro también puede deducirse del jacobiano. Justo a la derecha del punto fijo, una perturbación con $u>0$ y $v=0$ tiene $v'\simeq-bu<0$. En la representación usual, con $x$ horizontal e $y$ vertical, el movimiento oscilatorio local gira en sentido horario. Esta conclusión orienta la lectura del retrato cercano al equilibrio y evita atribuir un orden temporal incompatible con las ecuaciones.

El período debe estimarse a partir de recurrencias temporales equivalentes y compararse con $2\pi/a$ solo en el régimen próximo a Hopf. Fuera de ese régimen, la circulación por regiones con velocidades diferentes introduce correcciones que el jacobiano no contiene. El mecanismo cualitativo explica por qué pueden alternarse acumulación y consumo, mientras la trayectoria concreta determina cuánto dura cada parte del ciclo y qué amplitud alcanzan las concentraciones.

176.8. Predicciones sobre el ciclo límite

Las ecuaciones imponen restricciones exactas a cualquier solución periódica, aunque se desconozca su forma. Supóngase una órbita de período $T>0$ y defínase el promedio temporal de una función $q(t)$ mediante $\langle q\rangle=T^{-1}\int_0^Tq(t)\,dt$. El origen temporal puede elegirse en cualquier punto de la órbita. La periodicidad garantiza que la integral de la derivada de cada concentración se anula sobre un período completo.

La suma de los balances elimina los términos internos de conversión. Al integrar la identidad resultante entre dos instantes separados por un período, el cambio neto de $x+y$ es cero, de modo que la entrada acumulada iguala la eliminación acumulada. Se obtiene así una predicción independiente de la amplitud, de la forma de la órbita y de su distancia al umbral de Hopf

La igualdad afirma que la concentración media de $X$ coincide con su coordenada estacionaria, aunque sus valores instantáneos varíen. Esta coincidencia procede de un balance lineal y no implica simetría de la señal alrededor de su media; un máximo estrecho y elevado puede compensarse con un intervalo prolongado de valores menores. Tampoco determina los extremos, el desfase ni el período, porque muchas funciones periódicas diferentes comparten el mismo promedio temporal.

Integrado sobre un período, el balance de $Y$ indica que su producción total iguala el consumo autocatalítico total. Con el resultado anterior se obtiene una restricción sobre el promedio de un producto no lineal de concentraciones; ese promedio, fijado por los parámetros, es una velocidad media en las unidades adimensionales del modelo

En general, promediar un producto no equivale a multiplicar los promedios. La correlación temporal entre $x^2$ e $y$ influye en la velocidad media autocatalítica. Puede escribirse $\langle x^2y\rangle=\langle x^2\rangle\langle y\rangle+\operatorname{Cov}(x^2,y)$, donde la covarianza se define con el mismo promedio temporal. Además, $\langle x^2\rangle=a^2+\operatorname{Var}(x)$. Una sustitución ingenua de cada variable por su valor medio omite ambas correcciones.

Por esta razón, las identidades no fijan $\langle y\rangle=b/a$. Esa expresión corresponde exactamente al estado estacionario, donde las varianzas y covarianzas temporales se anulan. En una órbita no trivial, el valor medio de $Y$ depende conjuntamente de la variabilidad de $X$ y de su correlación con $Y$. Que una media dinámica coincida con una coordenada estacionaria debe justificarse por separado para cada variable.

Con estas restricciones puede comprobarse analíticamente una integración numérica. Una vez eliminado el transitorio y seleccionado un período completo, los promedios calculados deben aproximar $a$ y $ba$ dentro del error de integración y cuadratura. Las discrepancias pueden revelar ventanas incompletas, transitorios persistentes o resolución insuficiente. El cumplimiento de ambas identidades es una comprobación necesaria de consistencia, aunque no certifica todos los detalles de la órbita calculada.

176.9. Integración numérica y controles de convergencia

Para contrastar los balances del ciclo se integraron las ecuaciones con $a=1$, $b=2{,}01$ y $b=2{,}2$, desde $(x_0,y_0)=(1{,}1,b)$. Se usó DOP853 mediante solve_ivp de SciPy 1.17.0, con paso máximo $0{,}25$. Cada caso se calculó con tolerancias relativa/absoluta $10^{-9}/10^{-11}$ y $10^{-11}/10^{-13}$. Los horizontes fueron 4000 y 1000 unidades de tiempo, respectivamente, para dar margen a la relajación más lenta cerca del umbral. Estos resultados son simulaciones del modelo, sin datos experimentales.

Se registraron cruces ascendentes de $x=1$ y se analizó el último período completo. Los promedios se calcularon por cuadratura sobre la solución interpolada; los extremos se localizaron mediante optimización local dentro de ese período. Las cifras de la tabla permanecieron iguales con ambas tolerancias.

El primer período está próximo al límite $2\pi\approx6{,}283185$; el segundo muestra la corrección a distancia finita del umbral. Los extremos no son simétricos respecto de la media, igual a uno, especialmente para $b=2{,}2$. Los balances medios se cumplen con discrepancias menores que $10^{-9}$ en las evaluaciones realizadas. Esta concordancia comprueba la integración y la selección del período en los dos casos, aunque no demuestra unicidad global ni atracción desde todas las condiciones iniciales. El código y los resultados de la comprobación acompañan al conjunto de capítulos.

Un integrador adaptativo ajusta el paso temporal según una estimación del error local. La tolerancia absoluta controla errores cuando una concentración es pequeña; la relativa pondera el error respecto de su magnitud. Ambas deben declararse, junto con el algoritmo y la implementación. Para confiar en el resultado hay que repetir el cálculo con tolerancias más estrictas y comprobar que los observables relevantes se mantienen estables. Un algoritmo que termina sin advertencias todavía puede producir resultados insuficientemente convergidos.

El horizonte de integración debe superar las escalas de relajación del régimen estudiado. Cerca de Hopf, la parte real de los autovalores se aproxima a cero y los transitorios pueden prolongarse; una ventana aparentemente estacionaria puede conservar una deriva lenta de amplitud. Conviene separar explícitamente un intervalo descartado como transitorio de otro destinado al análisis y comprobar que las estadísticas y los tiempos de retorno permanecen estables al ampliar ambos.

Para explorar la cuenca más allá de los dos cálculos presentados, las condiciones iniciales deben variar dentro del cuadrante físico. Las coincidencias entre trayectorias tardías son evidencia numérica de atracción, con el alcance limitado del muestreo realizado. Una inicialización exactamente en el punto fijo puede permanecer allí incluso cuando este es inestable, salvo perturbaciones de redondeo. Por ello, la ausencia de movimiento desde esa condición particular no prueba su estabilidad.

La positividad analítica es un control adicional. Los métodos generales no siempre preservan concentraciones no negativas bajo cualquier paso temporal. Si aparecen valores negativos, deben revisarse las tolerancias, la resolución y la adecuación del método. Recortar arbitrariamente cada valor negativo a cero modifica las ecuaciones efectivamente integradas y puede deformar los balances. Cuando sea necesario, puede usarse un método que preserve la positividad; hay que documentar su formulación y volver a verificar la convergencia.

Los períodos pueden estimarse mediante cruces de una sección de Poincaré transversal a la órbita. Una opción es registrar los cruces de $x=a$ en una dirección especificada, siempre que $x'\neq0$ en los eventos considerados. Si $t_n$ denota el instante del cruce número $n$, los intervalos $T_n=t_{n+1}-t_n$ aproximan períodos sucesivos después del transitorio. La localización de eventos requiere interpolación controlada; tomar simplemente el punto almacenado más cercano introduce un error adicional.

La evaluación final debe combinar recurrencia geométrica, estabilidad de los períodos y comprobación de $\langle x\rangle=a$ y $\langle x^2y\rangle=ba$. Los máximos y mínimos solo deben comunicarse cuando se hayan calculado con resolución suficiente y resulten estables al refinar el procedimiento. Conservar parámetros, versión del software, tolerancias, criterio de descarte y definición de eventos permite distinguir una observación reproducible de una impresión visual dependiente del muestreo.

La elección de la variable observada afecta al diagnóstico de una oscilación. Cerca del punto de Hopf, las dos concentraciones oscilan con un desfase determinado por el autovector complejo del jacobiano. Una combinación lineal de ambas puede reducir mucho la amplitud visible, aunque la trayectoria en el plano de estados tenga una oscilación clara. Por eso conviene registrar al menos dos observables cuando se pretende reconstruir el ciclo y distinguir su fase de una simple fluctuación de amplitud.

La escala temporal de la medición también debe compararse con el período predicho. Un registro con pocos puntos por vuelta puede distorsionar los máximos y producir una frecuencia aparente diferente. En una observación corta, el acercamiento en espiral a un equilibrio puede confundirse con una órbita sostenida. Variar la duración y la frecuencia del muestreo permite comprobar si la amplitud converge a un valor persistente o sigue disminuyendo; así se evita atribuir al mecanismo químico un efecto generado por la observación.

176.10. Alcance y límites del modelo

El Brusselator muestra que una alimentación constante y una cinética no lineal pueden sostener oscilaciones internas. Su persistencia exige recursos y condiciones externas mantenidas; en un reactor real, fijar los parámetros supone reponer reactivos y gestionar productos o intercambios con el entorno. «Autónomo» significa que los coeficientes no dependen explícitamente del tiempo; el sistema no está aislado ni puede oscilar indefinidamente sin consumir energía libre.

Si se agotan los reservorios, sus concentraciones pasan a ser variables dinámicas y cambian las ecuaciones. Durante un intervalo limitado puede haber oscilaciones similares a las del sistema alimentado, con amplitud, frecuencia o duración variables según la disponibilidad de reactivos. Una larga sucesión de máximos aún puede ser transitoria. Para establecer el autosostenimiento hay que identificar el régimen externo que lo mantiene y su escala temporal de validez.

La aproximación determinista omite las fluctuaciones de poblaciones moleculares finitas. El ruido puede perturbar la fase de un ciclo estable y ensanchar la distribución de tiempos de retorno; cerca de un foco estable también puede excitar respuestas oscilatorias amortiguadas. Un pico espectral en una señal irregular no prueba, pues, que exista un ciclo límite determinista; su interpretación depende de la estabilidad subyacente y de la intensidad y estructura de las fluctuaciones.

Un autooscilador tiene una dinámica periódica atractora bajo condiciones externas constantes. Un sistema excitable permanece cerca de un estado estable y, tras una perturbación adecuada, puede realizar una gran excursión antes de recuperarse. Para distinguirlos hay que examinar atractores, umbrales y trayectorias de retorno, porque la amplificación química o los pulsos aislados no bastan para clasificar el régimen.

Los patrones espaciales plantean otro problema. La difusión convierte las concentraciones en campos dependientes de posición y tiempo y permite estudiar perturbaciones de distintas longitudes de onda. La bifurcación de Hopf es una inestabilidad de la dinámica temporal; la inestabilidad de Turing implica selección espacial por el acoplamiento entre reacción y transporte, con condiciones distintas. Ambos mecanismos forman parte del estudio de los sistemas fuera del equilibrio (Cross y Hohenberg, 1993).

El modelo homogéneo permite localizar una transición oscilatoria, interpretar su mecanismo y deducir balances medios exactos. Aplicarlo a una química específica exige justificar la reducción cinética, las escalas, la mezcla y el mantenimiento de los reservorios. El análisis lineal describe el entorno del punto fijo; la forma normal, una vecindad de la bifurcación; y la simulación da evidencia sujeta a convergencia. Las afirmaciones globales requieren argumentos adicionales acordes con su alcance.

Referencias

Prigogine, I. y Lefever, R. (1968). Symmetry breaking instabilities in dissipative systems. II. The Journal of Chemical Physics, 48, 1695–1700. DOI.

Cross, M. C. y Hohenberg, P. C. (1993). Pattern formation outside of equilibrium. Reviews of Modern Physics, 65, 851–1112. DOI.

Ferčec, B., Dukarić, M., Aybar, O. O. y Aybar, I. K. (2021). Supercritical Hopf bifurcations in two biochemical reaction systems. MATCH Communications in Mathematical and in Computer Chemistry, 85, 525–544. Texto completo.