Estructuras disipativas y formación de orden

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

177.1. Definición de estructura disipativa

Una estructura disipativa es una organización espacial o temporal mantenida por procesos irreversibles en un sistema que intercambia energía y, según el caso, materia con su entorno. Los flujos sostienen el desequilibrio que alimenta la organización; la materia se transforma o se transporta mientras la energía disponible se degrada por disipación. Su descripción exige identificar el sistema, sus fronteras, los intercambios con el entorno y las variables que revelan la organización. Una regularidad observable es insuficiente para establecer ese carácter; el criterio decisivo es que el patrón dependa dinámicamente de los flujos que lo mantienen durante un intervalo prolongado.

El impulso puede proceder de una diferencia de temperatura, de potencial químico o de otra magnitud que produzca transporte. Las fuentes externas mantienen esos desequilibrios mientras los procesos irreversibles consumen energía disponible. Un estado organizado puede conservar propiedades macroscópicas aproximadamente constantes aunque sus componentes materiales se renueven, continúen los intercambios y la producción de entropía sea positiva.

Conviene distinguir patrón, orden y autoorganización. Un patrón es una regularidad reconocible en los campos observados, como una modulación de concentración o una oscilación. El orden expresa correlaciones o restricciones cuya medida depende de las variables y de la escala elegidas. La autoorganización es la aparición de una organización mediante interacciones internas bajo condiciones externas determinadas. Ninguno de estos términos incluye en su definición la adaptación, la finalidad ni una función biológica.

Un cristal en equilibrio es un caso de contraste. Su periodicidad espacial puede persistir cuando se lo aísla y alcanza el equilibrio bajo condiciones apropiadas, ya que su red se conserva sin un flujo continuo. Asimismo, una configuración transitoria surgida durante una relajación puede mostrar formas regulares y después desaparecer. Clasificar estos casos requiere estudiar su mantenimiento, su estabilidad y su relación con las restricciones termodinámicas, además de describir su geometría visible.

Un criterio operativo es desconectar el impulso y observar la evolución. Si la organización depende del suministro, al retirarlo perderá su mecanismo de mantenimiento y tenderá a relajarse, aunque puedan quedar huellas materiales o transitorios largos. Esta prueba debe interpretarse junto con los tiempos característicos y las posibles barreras cinéticas. Una persistencia breve tras la desconexión no demuestra independencia del flujo; una desaparición rápida, a su vez, necesita una explicación dinámica.

El análisis de estas estructuras enlaza tres niveles complementarios. Los balances de intercambio establecen la compatibilidad termodinámica; el análisis de estabilidad del estado base identifica las perturbaciones que pueden crecer; las no linealidades determinan su eventual saturación. Los ejemplos de convección y reacción-difusión permiten separar estos niveles y comprobar qué conclusiones autoriza cada uno, sin atribuir al aspecto ordenado propiedades que aún no se han demostrado.

177.2. Equilibrio y balances de entropía

Para un sistema delimitado por fronteras, el balance de entropía distingue la producción interna del intercambio con el entorno. Esta separación permite describir sistemas abiertos sin confundir el cambio de una propiedad almacenada con los procesos irreversibles que ocurren en su interior. La elección de las fronteras debe ser explícita, porque un mismo proceso puede aparecer como intercambio o como producción interna al ampliar o reducir el sistema.

Aquí, $S$ es la entropía del sistema, $t$ el tiempo, $\Pi$ la tasa de producción interna y $\Phi$ el flujo neto de entropía hacia su interior. La entropía se mide en julios por kelvin y ambas tasas en julios por kelvin y segundo. Con esta convención, un valor negativo de $\Phi$ indica exportación neta. El intercambio puede incluir contribuciones del calor y de la materia que atraviesa las fronteras.

En un estado estacionario, la entropía almacenada permanece constante y el balance exige que la exportación compense la producción interna, de modo que puede haber producción positiva aunque la derivada temporal se anule. En equilibrio, en cambio, desaparecen las fuerzas termodinámicas macroscópicas que sostienen los flujos irreversibles. La igualdad entre entradas y salidas no basta para identificar un equilibrio, porque una corriente térmica estacionaria puede atravesar un sistema cuyas propiedades medias ya no cambian con el tiempo.

El orden espacial tampoco garantiza una entropía termodinámica baja. La medida visual de regularidad, la entropía de una distribución de observaciones y la entropía termodinámica pertenecen a descripciones diferentes, aunque puedan relacionarse mediante hipótesis adicionales. Para comparar dos estados hay que especificar las restricciones, las variables conservadas y los intercambios relevantes. Una fotografía del patrón no contiene la información necesaria para calcular la producción de entropía ni el balance completo.

Cerca del equilibrio, la termodinámica irreversible lineal relaciona flujos y fuerzas mediante coeficientes fenomenológicos. Las relaciones de reciprocidad de Onsager, válidas solo en el régimen de respuesta lineal próximo al equilibrio, requieren hipótesis de reversibilidad microscópica y el tratamiento adecuado de las simetrías temporales. Extenderlas directamente a una estructura impulsada lejos del equilibrio puede eliminar precisamente los mecanismos no lineales responsables de su aparición y estabilidad.

No existe un principio universal de máxima o mínima producción de entropía que seleccione todas las estructuras disipativas. Los resultados extremales disponibles tienen dominios de validez y restricciones particulares. La selección de una configuración requiere las ecuaciones dinámicas, los mecanismos de transporte, la geometría y las condiciones de frontera. El balance de entropía es una condición de consistencia y una herramienta de contraste, pero no determina qué perturbaciones crecerán ni qué organización terminará siendo estable.

177.3. Inestabilidad del estado base

El estudio de la formación de orden comienza con un estado base que satisface las ecuaciones y las condiciones de frontera, como una distribución homogénea de concentraciones o un perfil térmico conductor. Aunque sea estacionario, este estado puede perder estabilidad al cambiar la intensidad del impulso. La primera pregunta es si las perturbaciones pequeñas decaen o se amplifican cuando los parámetros externos se mantienen fijos durante el análisis de cada estado.

El campo de variables se escribe como el estado base más una perturbación. Al sustituir esta descomposición en las ecuaciones y conservar los términos de primer orden se obtiene un problema lineal, cuyos coeficientes dependen del estado base y de los parámetros de control. Esta aproximación describe la etapa inicial del crecimiento mientras la perturbación sea lo bastante pequeña para que los términos descartados contribuyan poco a su evolución.

En esta expresión, $\delta\mathbf{u}$ es la perturbación del vector de campos, $\mathbf{r}$ la posición, $\varepsilon$ una amplitud inicial pequeña y $\mathbf{v}_{q}$ un modo espacial etiquetado por $q$. El autovalor $\sigma$ tiene parte real $\sigma_r$ y parte imaginaria $\sigma_i$; la primera determina el crecimiento o decaimiento exponencial y la segunda fija la frecuencia angular de una oscilación. Los campos físicos se recuperan tomando combinaciones reales de los modos complejos.

Un parámetro de control mide la intensidad de una condición externa, como una diferencia térmica o una alimentación química. El umbral lineal se alcanza cuando la mayor parte real del espectro cruza cero. La geometría y las fronteras modifican los modos permitidos y, por ello, el valor crítico. Un umbral calculado para un dominio infinito no puede trasladarse sin revisión a un recipiente finito, donde las restricciones espaciales pueden excluir la perturbación más favorable.

En una inestabilidad estacionaria, un autovalor crítico real atraviesa cero; en una oscilatoria, un par complejo conjugado cruza el eje imaginario con frecuencia distinta de cero. Esta diferencia determina si la primera desviación observable tiende a formar una modulación persistente o una variación temporal. Sin embargo, la clasificación lineal no establece la estabilidad del estado emergente, que exige un análisis no lineal, ni excluye bifurcaciones posteriores al aumentar el impulso.

Las fluctuaciones siembran los modos inestables. Pueden proceder de ruido térmico, heterogeneidades del medio o pequeñas perturbaciones de preparación. Su amplificación se alimenta del desequilibrio que mantienen los flujos externos, que son el suministro energético del patrón. Cuando el crecimiento alcanza amplitudes apreciables, las interacciones no lineales pasan a ser esenciales y el análisis debe extenderse para determinar la saturación, la competencia entre modos o la evolución hacia estados temporalmente más complejos.

177.4. Convección de Rayleigh-Bénard

La convección de Rayleigh-Bénard aparece en una capa de fluido calentada desde abajo cuando la flotación supera los efectos estabilizadores de la viscosidad y de la difusión térmica. El estado base conductor transporta calor sin movimiento macroscópico del fluido. Una perturbación ascendente puede llevar material más caliente hacia regiones superiores, donde la diferencia de densidad favorece su ascenso. La intensidad de esta realimentación depende tanto del calentamiento como de las propiedades del fluido.

El número de Rayleigh $\mathrm{Ra}$ es adimensional. Aquí, $g$ es la aceleración gravitatoria, en metros por segundo cuadrado; $\alpha$, el coeficiente de expansión térmica, en inversos de kelvin; y $\Delta T$, la diferencia de temperatura entre la frontera inferior y la superior, en kelvin. La altura de la capa es $H$, en metros. La viscosidad cinemática $\nu$ y la difusividad térmica $\kappa$ se expresan en metros cuadrados por segundo.

Como el número de Rayleigh depende del cubo de la altura, a igualdad de propiedades y de diferencia térmica una capa más profunda da a la flotación una influencia relativa mayor. El cociente compara mecanismos que actúan simultáneamente, por lo que su interpretación exige conservar las hipótesis del modelo. Si las propiedades cambian considerablemente con la temperatura o aparecen otros efectos físicos relevantes, la reducción a este único parámetro puede ser insuficiente.

El valor crítico aproximado de 1707,76 corresponde a una capa horizontal lateralmente infinita, con fronteras rígidas isotérmicas, calentada desde abajo y descrita mediante la aproximación de Boussinesq. Esta aproximación retiene la variación de densidad en la fuerza de flotación y trata las demás propiedades como aproximadamente constantes. El valor citado es propio de ese problema; otras condiciones mecánicas, térmicas o geométricas producen espectros distintos y pueden modificar el umbral de aparición.

Por debajo del umbral, la conducción tiene disipación positiva y carece de rollos convectivos estables surgidos de esa inestabilidad lineal. Este régimen demuestra que puede haber disipación sin patrón macroscópico. Al superar el umbral, ciertas perturbaciones crecen y pueden organizarse en rollos que transportan calor mediante circulación. La estructura espacial se mantiene mientras el gradiente térmico suministra la energía necesaria para compensar la disipación del movimiento y del transporte térmico.

La predicción del inicio no fija por completo la configuración observada. Las paredes laterales, la relación entre las dimensiones horizontales y la altura, la preparación inicial y las interacciones no lineales intervienen en la orientación y la estabilidad de los rollos. El transporte de calor y la dinámica del campo de velocidad son observables complementarios para contrastar el modelo. Si se retira el calentamiento diferencial, el gradiente se relaja y la circulación pierde la fuente que sostenía su organización.

177.5. Reacción-difusión

Los sistemas de reacción-difusión combinan transformaciones locales y transporte espacial. En una descripción de dos componentes, las variables son concentraciones o campos de actividad cuyo significado físico debe especificarse en cada aplicación. La difusión tiende a suavizar los gradientes de cada componente, pero su acoplamiento con una cinética multivariable puede desestabilizar un estado homogéneo. El mecanismo de Turing muestra cómo una interacción local estable puede generar organización espacial cuando los componentes se transportan con velocidades diferentes (Turing, 1952).

Aquí, $x$ e $y$ son los campos dependientes de posición y tiempo; $f$ y $g$ son sus tasas locales de reacción. Los coeficientes positivos $D_x$ y $D_y$ son difusividades y $\nabla^2$ es el operador laplaciano espacial. En variables dimensionales, cada tasa tiene unidades de concentración por tiempo y cada difusividad de longitud al cuadrado por tiempo. El modelo supone transporte difusivo ordinario, coeficientes constantes y ausencia de advección explícita.

En el estado homogéneo estacionario se anulan simultáneamente ambas tasas de reacción. El jacobiano $J$ reúne las derivadas parciales de esas tasas respecto de las dos variables, evaluadas en dicho estado. La cinética sin dependencia espacial es linealmente estable cuando la traza de $J$ es negativa y su determinante positivo. Estas condiciones garantizan partes reales negativas de los dos autovalores, pero todavía no incorporan el efecto de las perturbaciones espaciales.

El número de onda $q$ caracteriza un modo propio del laplaciano con autovalor $-q^2$; $I$ es la identidad y $\sigma$ su tasa de crecimiento. Las raíces de la ecuación característica definen la relación de dispersión; el crecimiento depende de la parte real de cada raíz. En un dominio finito, solo los números de onda compatibles con la geometría y las condiciones de frontera pertenecen al espectro físicamente admisible.

Con cinética estable y difusión diagonal positiva, la traza modal permanece negativa. Si el determinante modal se vuelve negativo para algún número de onda no nulo, los autovalores reales tienen signos opuestos y aparece una dirección creciente. Esta condición identifica la inestabilidad estacionaria de Turing en el sistema de dos componentes considerado. Con difusividades iguales, ambos autovalores se desplazarían hacia partes reales menores y no se produciría este mecanismo de desestabilización.

Las condiciones pueden escribirse directamente con las derivadas de la cinética. Si $J$ tiene entradas $f_x,f_y,g_x,g_y$, evaluadas en el estado base, entonces

Es una parábola convexa en $q^2$. Además de $\operatorname{tr}J<0$ y $\det J>0$, para que tome valores negativos en algún $q^2>0$ se necesita

La primera desigualdad sitúa el mínimo en un número de onda físico; la segunda lo vuelve negativo. Difusividades distintas son necesarias en este modelo, aunque no suficientes, porque también deben cumplirse las condiciones sobre la reacción. Las desigualdades corresponden a dos especies con difusión diagonal positiva y coeficientes constantes. Otros transportes, como la difusión cruzada, o más componentes requieren su propio problema espectral.

La banda inestable solo indica escalas candidatas. Su realización exige que el dominio admita modos dentro de ella y que el suministro mantenga la cinética; después intervienen la saturación y la competencia no lineal.

177.6. El caso del Brusselator

El Brusselator es un modelo idealizado de cinética química no lineal que permite verificar explícitamente una inestabilidad inducida por difusión, porque su estado homogéneo, su jacobiano y su banda de números de onda inestables pueden calcularse directamente. El ejemplo pertenece al estudio de inestabilidades y ruptura de simetría en sistemas abiertos desarrollado por Prigogine y Lefever (1968). Los valores que se usan son parámetros del modelo y no tienen carácter empírico.

Las variables $x$ e $y$ son concentraciones adimensionales; las funciones $f$ y $g$ son las tasas de reacción en el tiempo reescalado. Los parámetros positivos $a$ y $b$ controlan la alimentación idealizada. El estado homogéneo estacionario tiene concentraciones $x_*=a$ e $y_*=b/a$, donde el asterisco marca los valores de ese estado. La alimentación externa hace posible un régimen estacionario con actividad química persistente.

Se toman los parámetros adimensionales $a=1$, $b=1{,}8$, $D_x=1$ y $D_y=10$, donde las dos últimas cantidades son las difusividades reescaladas. El jacobiano evaluado en el estado homogéneo tiene primera fila $(0{,}8,1)$ y segunda fila $(-1{,}8,-1)$. Su traza es $-0{,}2$ y su determinante es $1$. Por tanto, una perturbación homogénea pequeña decae, ya que la cinética local satisface las condiciones de estabilidad del sistema de dos variables.

La función $h$ es el determinante modal, $q$ el número de onda adimensional y $D$ la matriz diagonal de difusividades. El polinomio es negativo cuando el cuadrado del número de onda está estrictamente entre 0,2 y 0,5. La traza modal vale $-0{,}2-11q^2$ y permanece negativa. Dentro de esa banda, el determinante negativo implica un autovalor positivo y otro negativo, de modo que ciertas perturbaciones heterogéneas se amplifican.

El resultado separa dos pruebas sobre el mismo estado homogéneo, que examinan sectores diferentes del espectro. La primera evalúa la estabilidad de la reacción aislada frente a perturbaciones uniformes; la segunda incorpora transporte y variación espacial. Los extremos de la banda corresponden a crecimiento lineal nulo y requieren cuidado al interpretar tiempos de observación finitos. Pertenecer a la banda continua tampoco garantiza que exista un modo admisible con ese número de onda.

Todo el cálculo es didáctico y puede verificarse con derivación y álgebra elemental. Para estudiar un recipiente concreto todavía deben fijarse la longitud, las fronteras y las perturbaciones iniciales. La amplitud final del patrón exige resolver las ecuaciones no lineales.

177.7. Dominio unidimensional con condiciones de Neumann

La banda inestable calculada para números de onda continuos debe confrontarse con el espectro de un dominio finito. Se considera un intervalo unidimensional de longitud $L$, con difusión ordinaria y condiciones homogéneas de Neumann para ambos componentes. Estas condiciones anulan la derivada normal de las concentraciones en los extremos y, bajo el transporte difusivo supuesto, equivalen a que no haya flujo material a través de esas fronteras. La coordenada espacial se denota por $s$ para no confundirla con la concentración $x$.

Aquí, $\phi_n$ es una función propia espacial del laplaciano, $q_n$ su número de onda y $n$ un entero no negativo que enumera los modos. El modo con índice cero es uniforme y conserva la estabilidad de la cinética local. Los índices positivos describen variaciones espaciales compatibles con las fronteras. Las longitudes se expresan en las unidades reescaladas usadas para las difusividades del ejemplo, por lo que la comparación con la banda es dimensionalmente consistente.

Para $L=5$, el cuadrado del número de onda del primer modo no uniforme es aproximadamente 0,394784. Este valor cae dentro de la banda abierta entre 0,2 y 0,5 obtenida para el Brusselator, por lo que ese modo es linealmente inestable. El modo con $n=2$, de valor aproximado 1,579, queda fuera, igual que todos los índices superiores, cuyos cuadrados son todavía mayores; en este dominio solo el primer modo espacial pertenece a la banda inestable.

Para $L=3$, el primer valor no uniforme es aproximadamente 1,0966, mayor que el extremo superior de la banda. Como la sucesión crece con el cuadrado del índice, ningún modo espacial permitido cae en el intervalo inestable. El modo uniforme sigue siendo estable por las propiedades del jacobiano. Así, modificar solo la longitud puede suprimir la inestabilidad lineal espacial sin cambiar las tasas de reacción, las difusividades ni la estabilidad de la cinética homogénea.

Para cada índice positivo, la condición exacta de entrada en la banda es

El primer modo es inestable para $4{,}443<L<7{,}025$, aproximadamente; el segundo lo es para $8{,}886<L<14{,}050$. Por tanto, aumentar el tamaño no produce necesariamente una secuencia monótona de estabilidad en este intervalo. Con $L=8$, el primer modo queda por debajo de la banda y el segundo por encima, de modo que vuelve a faltar un modo creciente; con $L=10$, el segundo entra en ella. La banda continua y la lista discreta de modos cumplen funciones distintas en la predicción.

En $L=5$, el autovalor positivo da además un tiempo de crecimiento. Su valor es aproximadamente $0{,}0447$ en tiempo inverso adimensional, por lo que la amplitud necesita unas cien unidades temporales para multiplicarse por cien en la aproximación lineal. Una observación breve podría no detectar el patrón naciente aun cuando el estado base sea inestable. Esta estimación describe la amplificación inicial y pierde validez al entrar en el régimen de saturación.

La longitud filtra el espectro, ya que al variar el tamaño los valores discretos se desplazan respecto de la banda continua y pueden entrar o salir de ella. La forma cosenoidal muestra también que las fronteras restringen la fase espacial de las perturbaciones. Este razonamiento presupone condiciones de Neumann exactas y coeficientes uniformes; los intercambios en los extremos o las heterogeneidades materiales darían un problema de valores propios diferente y exigirían repetir el cálculo.

La existencia de un modo creciente hace posible, en la teoría lineal, la formación de un patrón, pero no determina su amplitud final, su perfil saturado ni su estabilidad frente a perturbaciones posteriores. Las no linealidades pueden introducir armónicos y desplazar la concentración media. Incluso con un único modo inicialmente inestable, la solución final necesita un análisis adicional.

177.8. Saturación no lineal y ecuaciones de amplitud

La teoría lineal predice crecimiento exponencial y pierde validez cuando la perturbación alcanza una magnitud apreciable. Una ecuación de amplitud describe, cerca de ciertos umbrales, la evolución lenta de la intensidad de un modo crítico. Su forma depende de las simetrías, de la naturaleza de la bifurcación y de las interacciones que sobreviven a la reducción. El modelo escalar real más sencillo permite estudiar cómo un crecimiento inicial puede quedar limitado por una realimentación no lineal.

La variable real A es una amplitud con signo; μ mide la distancia reescalada al umbral y g cuantifica la saturación cúbica. El tiempo t pertenece a la escala de la descripción reducida. Si la amplitud es dimensional, μ tiene unidades de tiempo inverso y g de tiempo inverso por amplitud al cuadrado. La simetría de cambio de signo significa que ambos signos corresponden a configuraciones equivalentes bajo una transformación apropiada del sistema.

Para valores negativos de μ, la amplitud nula es estable. Cuando el parámetro se hace positivo, esa solución pierde estabilidad y aparecen dos estados estacionarios no nulos, de amplitudes A = ±√(μ/g). La derivada del campo dinámico respecto de la amplitud vale −2μ en ambas ramas, por lo que las pequeñas desviaciones decaen. La dependencia de raíz cuadrada describe el crecimiento continuo del orden dentro del régimen de validez del modelo reducido.

Esta ecuación ilustra una bifurcación supercrítica con simetría de inversión, pero no se ha derivado exactamente del Brusselator presentado antes. Obtener una reducción específica requeriría identificar el modo crítico, proyectar las ecuaciones, calcular los coeficientes no lineales con las fronteras correspondientes y comprobar qué simetrías conserva el dominio. Sin esos pasos, no está justificado asignar un valor al coeficiente cúbico ni interpretar la amplitud como una concentración directamente medible.

Cuando g es negativo, el término cúbico refuerza el crecimiento y se necesita una contribución estabilizadora de orden superior. El coeficiente positivo c introduce una saturación de quinto orden; sus unidades aseguran que todos los términos sean tasas de cambio de amplitud. Esta forma puede producir ramas de amplitud finita antes de que el estado nulo pierda estabilidad, rasgo propio de una bifurcación subcrítica. El detalle de esas ramas depende de los coeficientes y del intervalo de parámetros en el que la expansión sea válida.

Si coexisten estados estables separados por una rama inestable, un barrido suficientemente lento puede mostrar histéresis. El ruido, las perturbaciones finitas y la velocidad de variación del control alteran las transiciones observadas. La ecuación de amplitud debe contrastarse dentro de su régimen aproximado, especialmente lejos del umbral.

177.9. Selección de patrón, defectos y ruido

Cuando varios modos son inestables, la relación de dispersión establece candidatos que compiten durante la evolución. El de mayor crecimiento puede dominar al principio si las semillas tienen magnitudes comparables, pero esa ventaja temprana no garantiza la selección final. Las interacciones no lineales, las cantidades conservadas y las fronteras modifican la competencia. La formación de patrones requiere estudiar tanto el acceso a una configuración como su estabilidad una vez establecida, distinción central en la teoría de sistemas extendidos (Cross y Hohenberg, 1993).

Un patrón puede caracterizarse por su amplitud, longitud de onda, orientación y fase. En dominios extensos, estas propiedades pueden variar lentamente entre regiones y producir configuraciones localmente ordenadas sin coherencia global. Las ecuaciones de amplitud espacial incorporan esas modulaciones cuando hay una separación adecuada de escalas; sin ella, una descripción con pocos modos puede omitir procesos importantes. La resolución necesaria depende del observable que se desea explicar y del intervalo temporal en el que se evalúa la organización.

Los defectos son regiones donde la descripción periódica local pierde regularidad o donde la fase deja de estar bien definida. Pueden separar dominios de distinta orientación, mediar cambios del número de bandas o participar en dinámicas persistentes. Su presencia no implica necesariamente desorden completo ni un fallo del modelo. Las tasas de creación, desplazamiento y aniquilación de defectos informan sobre la relajación y sobre la posibilidad de alcanzar un estado estructurado estable en el tamaño de sistema considerado.

Las fronteras pueden fijar orientaciones, favorecer fases y absorber o generar defectos. En un recipiente pequeño, estas restricciones pueden afectar a casi todo el patrón; en uno grande, su influencia puede competir con procesos internos. La historia de preparación también influye cuando hay varios atractores, ya que dos realizaciones con los mismos parámetros finales pueden terminar en configuraciones diferentes si recorren trayectorias distintas en el espacio de estados. Esa dependencia exige distinguir entre selección reproducible y multiestabilidad al comparar teoría y observación.

El ruido puede sembrar modos, desplazar fases y facilitar transiciones entre configuraciones. Sus efectos dependen de la forma en que entra en las ecuaciones y de su intensidad, correlación temporal y estructura espacial. Cerca de un umbral, la relajación lenta puede aumentar la sensibilidad a perturbaciones débiles. Un patrón fluctuante puede conservar correlaciones estadísticas estables aunque su realización instantánea cambie; por ello, la estabilidad de un campo estacionario y la estacionariedad de sus estadísticas son propiedades diferentes que conviene evaluar por separado.

También deben distinguirse los estados espacialmente estructurados y estacionarios de los estados oscilatorios o de variación irregular. La estabilidad de cada solución se determina perturbándola y examinando su evolución, sin inferirla de la estabilidad del estado homogéneo original. En todos estos regímenes, el suministro de energía y materia sostiene el desequilibrio disponible. La selección concreta depende de la dinámica y de las restricciones, mientras que el mantenimiento disipativo exige verificar los flujos y la respuesta del sistema cuando se interrumpe el impulso.

177.10. Pruebas y criterios de contraste

Identificar una estructura disipativa requiere evidencias dinámicas y termodinámicas compatibles, que expliquen la forma visible a partir de un estado base, un mecanismo de inestabilidad y un suministro mantenido. El contraste especifica las magnitudes observadas, las condiciones externas, el tiempo de registro y las hipótesis del modelo.

Variar el parámetro de control permite localizar cambios de estabilidad y comparar el umbral observado con la predicción; variar el dominio modifica el espectro espacial y puede admitir o excluir un modo creciente, como en el ejemplo del Brusselator. La comparación debe considerar la incertidumbre, los transitorios y la proximidad al umbral. Desconectar el impulso permite comprobar si el patrón depende del suministro; conviene registrar la relajación de campos y gradientes, pues las reservas internas o los tiempos de transporte pueden retrasar su desaparición. Los balances de energía, materia y entropía deben usar las mismas fronteras y convenciones.

Los espectros espaciales y temporales permiten contrastar escalas, tasas de crecimiento y frecuencias; la fase temprana informa sobre la relación de dispersión y el régimen tardío, sobre la saturación y la estabilidad no lineal. Deben separarse el ruido instrumental y las fluctuaciones del sistema. Un máximo espectral aislado no demuestra un mecanismo de Turing ni una bifurcación concreta.

El engrosamiento, o coarsening, aumenta progresivamente la escala característica y puede prolongarse. Una instantánea puede parecer estabilizada aunque cambien la longitud de onda, la densidad de defectos o las correlaciones. Identificar un atractor requiere comprobar el comportamiento asintótico y la estabilidad frente a las condiciones iniciales, dentro de los límites de duración y resolución de la observación.

De la organización espacial no puede extrapolarse una organización funcional, que exige criterios adicionales de desempeño, de acoplamiento con el entorno o de consecuencias causales. Los ejemplos establecen relaciones verificables entre impulso, transporte, inestabilidad y saturación, dentro de las ecuaciones, los parámetros y las fronteras examinados.

Referencias

Turing, A. M. (1952). The chemical basis of morphogenesis. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 237, 37–72. DOI.

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

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