Autocatálisis: amplificación y límites de recursos

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

175.1. Significado formal de la autocatálisis

La autocatálisis ocurre cuando un producto cataliza su propia formación. En el esquema mínimo, una unidad preexistente de $X$ transforma el sustrato $S$ y genera una unidad adicional; la amplificación, que depende de una semilla y de materia disponible, sigue la reacción

En la reacción, las letras denotan especies; en las ecuaciones cinéticas, concentraciones en milimoles por litro. El tiempo se expresa en horas y $k$ caracteriza la velocidad bajo condiciones experimentales especificadas. El cambio neto consume una unidad de $S$ y produce una de $X$; la unidad inicial es una de las dos unidades finales, que son equivalentes.

La presencia del producto debe cambiar la velocidad de su formación. La retroalimentación positiva por calentamiento, transporte o regulación externa no es autocatálisis si el producto no cataliza la transformación. Una curva acelerada tampoco la demuestra, porque puede deberse a inducción, a la activación de un precursor o a variaciones del medio.

El esquema puede representar una reacción elemental o resumir asociación, conversión y liberación. La ley efectiva puede incluir saturación, inhibición o intermediarios; la estequiometría global no determina sus exponentes de concentración. La teoría de redes distingue estructura estequiométrica, hipótesis cinéticas y dinámica (Feinberg, 1987).

La amplificación requiere recursos materiales y condiciones energéticas compatibles. En un recipiente cerrado, el sustrato disminuye; en un reactor alimentado, la entrada puede compensar el consumo, aunque la salida introduce pérdidas. Las fronteras cambian los balances y los estados estacionarios sostenibles.

La verificación compara velocidades con cantidades iniciales distintas de producto, controlando sustrato, temperatura y composición del medio. También identifica intermediarios, descarta impurezas catalíticas y separa una reacción basal independiente de la semilla. El modelo mínimo excluye esa producción basal, de modo que una concentración inicial nula permanece nula. La autocatálisis amplifica producto existente; la primera unidad exige un proceso adicional identificado.

175.2. Ley de acción de masas y logística sin pérdidas

Considérese un sistema cerrado, homogéneo y de volumen fijo, donde la transformación autocatalítica mínima es el único proceso. La mezcla uniforme permite describir cada especie mediante su concentración. Bajo acción de masas, la velocidad por unidad de volumen es proporcional al producto de las concentraciones de sustrato y autocatalizador. Con temperatura constante y sin transporte limitante, reacción inversa, producción basal ni degradación del producto, los balances son

Con $X(t)$ y $S(t)$ en unidades de concentración y $t$ en unidades de tiempo, la consistencia dimensional exige que $k$ tenga unidades de concentración inversa por tiempo inverso; con milimoles por litro y horas, litros por milimol y hora. La estequiometría unitaria hace coincidir el consumo de sustrato con la formación neta de producto.

Al sumar las ecuaciones, la derivada de la suma se anula y $M=S(t)+X(t)=S_0+X_0$ permanece constante; el subíndice cero indica concentraciones iniciales. Esta conservación del inventario convertible presupone equivalencia estequiométrica unitaria; en una química real con masas moleculares o rendimientos distintos se requiere una combinación ponderada. Al sustituir $S=M-X$, el sistema se reduce a una sola variable dinámica, que obedece

La ecuación logística se deriva así del balance material, con $r$ en tiempo inverso y $M$ en concentración. El límite superior corresponde al consumo de una reserva finita. Aunque la forma aparece en otros contextos, aquí cada término tiene interpretación química; la contribución cuadrática negativa refleja la disminución del sustrato disponible al aumentar el producto.

El intervalo físico $0\leq X\leq M$ es invariante, porque en los extremos la derivada es nula y en el interior es positiva para $k>0$. Los equilibrios son $X=0$ y $X=M$; en ellos, la derivada $f'(X)=k(M-2X)$ vale $kM$ y $-kM$, respectivamente. En esta clase de conservación, una semilla positiva crece y la conversión completa atrae las trayectorias interiores.

Para $X>0$, la tasa específica instantánea es $X^{-1}dX/dt=kS$, por lo que la amplificación por unidad de producto disminuye al agotarse el recurso. La reducción logística deja de ser exacta con muerte, entrada, salida o rendimiento distinto. Cada modificación exige revisar los balances de especies y determinar si queda una cantidad conservada que permita otra reducción escalar.

175.3. Solución analítica y fases

Para una semilla inicial positiva y menor que el inventario total, la logística admite una solución explícita. La separación de variables integra el cociente entre producto y sustrato restante; la constante de integración queda determinada por $X_0$. Con $M$ en unidades de concentración, $k$ en concentración inversa por tiempo inverso y $t$ en tiempo, el exponente es adimensional. La solución es

Equivalentemente, $\ln[X/(M-X)]=kMt+\ln[X_0/(M-X_0)]$. Los argumentos de los logaritmos son cocientes de concentraciones y, por tanto, adimensionales. El logaritmo de la razón producto/sustrato aumenta linealmente con el tiempo; esta forma permite estimar parámetros si la conservación postulada es válida. Un ajuste lineal satisfactorio comprueba el modelo, pero no identifica de manera única el mecanismo molecular.

Cuando el producto es una fracción pequeña del inventario, el sustrato permanece próximo a su valor inicial. Entonces $X(t)\approx X_0e^{kS_0t}$; si la semilla es muy pequeña respecto de $M$, $S_0\approx M$ y la tasa exponencial es $kM$. La aproximación pierde precisión cuando la conversión deja de ser despreciable.

El crecimiento absoluto se acelera inicialmente solo si la trayectoria comienza por debajo de la mitad del inventario, ya que $d^2X/dt^2=k(M-2X)\,dX/dt$. La inflexión ocurre en $X=M/2$, donde la velocidad alcanza $kM^2/4$, en unidades de concentración por tiempo. Después disminuye por la escasez de sustrato. Si la condición inicial supera esa mitad, hay desaceleración desde el comienzo.

En tiempos largos, la concentración se aproxima asintóticamente a $M$, mientras $M-X(t)$ decae proporcionalmente a $e^{-kMt}$. La conversión completa no se alcanza en tiempo finito en esta descripción continua; definir un tiempo de finalización exige fijar una tolerancia o fracción de conversión. Los tiempos para alcanzar fracciones especificadas son finitos y se obtienen invirtiendo la solución.

La condición $X_0=0$ se trata por separado porque la fórmula divide por la semilla; la ecuación diferencial da $X(t)=0$ para todo tiempo. Si $X_0=M$, la solución permanece en conversión completa. Semillas positivas cada vez menores implican un retraso creciente antes de alcanzar una concentración dada, aunque comparten el límite final determinista. La probabilidad de establecimiento requiere considerar explícitamente los eventos discretos cuando existen pérdidas.

175.4. Caso calculable didáctico

Este ejemplo es analítico y didáctico; sus valores no proceden de observación empírica. Se toma un inventario total $M=1\ \mathrm{mmol/L}$, una constante $k=1\ \mathrm{L/(mmol\cdot h)}$ y una semilla $X_0=0{,}01\ \mathrm{mmol/L}$. El sustrato inicial es $S_0=0{,}99\ \mathrm{mmol/L}$.

La tasa específica inicial es $kS_0=0{,}99\ \mathrm{h^{-1}}$, mientras que la velocidad inicial de producción es $kS_0X_0=0{,}0099\ \mathrm{mmol/(L\cdot h)}$. Son magnitudes distintas, pues la primera divide la producción por la cantidad presente y la segunda mide el aumento de concentración por hora. El parámetro logístico $kM=1\ \mathrm{h^{-1}}$ aproxima la tasa específica inicial; la diferencia se debe a la fracción de inventario ya convertida en producto. Con estos valores, la solución es

Despejando la exponencial, el tiempo para alcanzar una fracción $q=X/M$, con $0<q<1$, es $t_q=(kM)^{-1}\ln[(M/X_0-1)q/(1-q)]$. Si la fracción objetivo supera la inicial, el logaritmo es positivo y el tiempo pertenece al futuro de la trayectoria; si coinciden, el tiempo es cero. La expresión es dimensionalmente coherente, porque el logaritmo actúa sobre un número adimensional y el factor exterior tiene dimensión temporal. Para la mitad y el noventa por ciento de conversión resulta

A mitad de conversión, producto y sustrato tienen concentraciones de $0{,}5\ \mathrm{mmol/L}$ y la velocidad absoluta alcanza su máximo, $0{,}25\ \mathrm{mmol/(L\cdot h)}$. Al noventa por ciento quedan $0{,}1\ \mathrm{mmol/L}$ de sustrato; la tasa específica desciende a $0{,}1\ \mathrm{h^{-1}}$ y la velocidad absoluta a $0{,}09\ \mathrm{mmol/(L\cdot h)}$, aunque hay mucho más autocatalizador que al inicio.

El intervalo entre ambos objetivos es aproximadamente $2{,}197\ \mathrm{h}$, equivalente a $\ln 9/(kM)$. En una trayectoria logística con iguales parámetros, este intervalo no depende de la semilla, que solo desplaza los tiempos de llegada mediante el término logarítmico inicial. Con las demás hipótesis cinéticas fijas, aumentar $k$ comprime la escala temporal y modificar $M$ cambia tanto la concentración límite como el parámetro temporal.

La extrapolación exponencial sirve solo mientras el consumo es pequeño. Extenderla a conversiones altas sobreestima la producción y puede superar el inventario disponible; la solución exacta evita esto mediante el balance de sustrato. La conversión completa requiere tiempos progresivamente mayores, de modo que una reacción terminada depende de la tolerancia adoptada. Convertir estos resultados en una predicción experimental exige contrastar mezcla, mecanismo, pérdidas y constancia de parámetros.

175.5. Pérdidas $X\to D$ y conservación

Supóngase que el producto se transforma irreversiblemente en una especie inactiva $D$, con constante de pérdida $\delta>0$ (tiempo inverso). Aquí $D(t)$ denota la concentración degradada; más adelante, la letra designará la dilución del reactor continuo. En un recipiente cerrado de volumen fijo y con equivalencia estequiométrica unitaria entre sustrato, producto activo y degradado, los balances son

Las dos primeras derivadas suman $-\delta X$, de modo que el inventario activo $S+X$ disminuye mientras haya producto. La cantidad conservada incluye la especie inactiva, $S+X+D=C$, con $C=S_0+X_0+D_0$ constante y expresada como concentración. Esto supone que el material degradado permanece dentro del sistema; si sale del recipiente, también se contabiliza lo acumulado en la corriente de salida.

Sumar una pérdida lineal a la logística y mantener $S=M-X$ es inconsistente con este mecanismo, porque esa sustitución requiere que el inventario activo sea constante y la degradación lo reduce. Hay que mantener el sistema para sustrato y producto o usar una reducción compatible con el balance real. Una logística con pérdida podría representar otra química, como un reciclaje explícito que devuelve al sustrato utilizable cada unidad perdida.

Para $X>0$, la tasa específica es $kS-\delta$; el producto crece si $S>\delta/k$, queda instantáneamente estacionario en el umbral y decrece por debajo. Como el sustrato no aumenta, el umbral puede cruzarse; una semilla inicialmente creciente puede entonces alcanzar un máximo y luego decaer. El umbral $\delta/k$ tiene unidades de concentración y expresa el balance entre producción autocatalítica y pérdida.

En el plano de concentraciones, para valores positivos, $dX/dS=-1+\delta/(kS)$. Al integrar se obtiene $X-X_0=-(S-S_0)+(\delta/k)\ln(S/S_0)$. No es una conservación lineal, pero relaciona estados de una misma trayectoria. Si $S_0>0$, puede quedar sustrato residual cuando desaparece la actividad; disponer de materia no garantiza que haya catalizador suficiente para transformarla.

Cuando el sustrato se agota, $dX/dt\approx-\delta X$ y la actividad decrece exponencialmente. En el sistema cerrado con pérdida positiva, el modelo continuo predice que desaparece asintóticamente. A escala molecular, la última unidad puede degradarse en tiempo finito y causar extinción exacta; esta extinción probabilística difiere del acercamiento continuo a cero. Eliminado el último autocatalizador, el mecanismo no contiene ninguna vía para reiniciar la producción.

175.6. Reactor continuo (quimiostato)

Un reactor continuo perfectamente mezclado recibe sustrato a concentración $S_{\mathrm{in}}$ y expulsa una corriente con la composición interna. El volumen permanece constante al coincidir ambos caudales. La dilución $D=F/V$ (caudal volumétrico/volumen) tiene unidades de tiempo inverso. Si la alimentación no contiene autocatalizador y no hay muerte adicional al arrastre, los balances son

Las concentraciones pueden expresarse en milimoles por litro; $k$ tiene unidades de concentración inversa por tiempo inverso. Para $T=S+X$, la suma de balances da $\dot T=D(S_{\mathrm{in}}-T)$ y $T(t)=S_{\mathrm{in}}+[T(0)-S_{\mathrm{in}}]e^{-Dt}$. Así, el total activo relaja a $S_{\mathrm{in}}$, pero no se conserva en general; $T=S_{\mathrm{in}}$ es invariante si se cumple inicialmente.

El equilibrio de lavado es $(S,X)=(S_{\mathrm{in}},0)$. Una perturbación de producto crece a tasa específica aproximada $kS_{\mathrm{in}}-D$; hay invasión si $kS_{\mathrm{in}}>D$ y decaimiento si $kS_{\mathrm{in}}<D$. En la igualdad, el punto no es hiperbólico y deben considerarse términos no lineales. La comparación requiere expresar producción y retirada en la misma unidad temporal. El equilibrio activo es

El estado activo existe si $kS_{\mathrm{in}}>D$. En él, el balance entre reproducción y lavado fija el sustrato residual y la alimentación determina el producto sostenido. El Jacobiano general tiene entradas $J_{11}=-D-kX$, $J_{12}=-kS$, $J_{21}=kX$ y $J_{22}=kS-D$, todas de dimensión tiempo inverso. En el estado activo resulta

Para $D>0$ y $X^*>0$, la traza es negativa y el determinante positivo, por lo que el estado es localmente asintóticamente estable. Sus autovalores, $-D$ y $-kX^*$, son reales. En el lavado son $-D$ y $kS_{\mathrm{in}}-D$, de acuerdo con el criterio de invasión y con el cambio de estabilidad al superar el umbral.

La entrada material continua puede sostener actividad estacionaria; una dilución demasiado rápida arrastra el producto antes de compensar su salida. Los resultados suponen mezcla perfecta, cinética bilineal y ausencia de pérdidas adicionales. Saturación, retardos, inhibición o especies suplementarias pueden cambiar la estabilidad. El modelo mínimo relaja hacia el estado activo sin oscilar, aunque combine autocatálisis y alimentación continua.

175.7. Ejemplo numérico con muerte

Considérese un ejemplo didáctico de reactor continuo con $S_{\mathrm{in}}=1\ \mathrm{mmol/L}$, $k=1\ \mathrm{L/(mmol\cdot h)}$ y $D=0{,}2\ \mathrm{h^{-1}}$. Sin muerte adicional, el equilibrio activo es $S^*=0{,}2\ \mathrm{mmol/L}$ y $X^*=0{,}8\ \mathrm{mmol/L}$. La tasa de invasión desde el lavado vale $0{,}8\ \mathrm{h^{-1}}$. Estos valores proceden del modelo ideal y no describen una población o reacción medida experimentalmente.

Introdúzcase ahora una muerte o inactivación de primer orden con constante $\delta>0$, expresada en horas inversas, cuyo producto no se recicla a sustrato utilizable. El balance del sustrato conserva su forma anterior y el de producto pasa a ser $\dot X=kSX-(D+\delta)X$. Si el material inactivo se conserva en el reactor hasta salir, puede añadirse una tercera ecuación, aunque no interviene en la cinética activa supuesta. El equilibrio activo pasa a ser

La primera igualdad procede de exigir una concentración activa positiva con derivada nula. La segunda se obtiene sustituyéndola en el balance estacionario del sustrato. El factor $D/(D+\delta)$ aparece porque una parte del recurso transformado termina en material inactivo, de modo que el producto activo ya no cubre toda la diferencia entre alimentación y sustrato residual. La condición de existencia es $kS_{\mathrm{in}}>D+\delta$, que en este ejemplo equivale a $\delta<0{,}8\ \mathrm{h^{-1}}$.

El total activo satisface $\dot T=D(S_{\mathrm{in}}-T)-\delta X$, donde $T=S+X$. En el equilibrio resulta $T^*=S_{\mathrm{in}}-(\delta/D)X^*$, inferior a la alimentación cuando hay muerte y actividad positiva. Si se incluye una concentración inactiva $Z$ con balance $\dot Z=\delta X-DZ$, el total ampliado $S+X+Z$ sí relaja a $S_{\mathrm{in}}$. En particular, $Z^*=(\delta/D)X^*$ explica exactamente la fracción ausente del inventario activo.

La estabilidad local tampoco se deduce del carácter autocatalítico. Con $q=kX^*>0$, el Jacobiano activo tiene traza $-(D+q)$ y determinante $(D+\delta)q$, que garantizan estabilidad asintótica. El discriminante de su polinomio característico es $(D+q)^2-4(D+\delta)q$; puede ser negativo y producir aproximación amortiguada con autovalores complejos. Esa posibilidad no implica una oscilación autosostenida; un ciclo persistente exigiría propiedades dinámicas adicionales que habría que demostrar en las ecuaciones correspondientes.

De hecho, para concentraciones estrictamente positivas, el multiplicador de Dulac $1/X$ da una divergencia $-D/X-k<0$ en este sistema bilineal, con lo que el criterio excluye órbitas periódicas interiores. La muerte puede reducir la actividad estacionaria, elevar el sustrato residual y modificar la forma del transitorio, pero no crea automáticamente un oscilador. Una calibración experimental tendría que identificar por separado dilución e inactivación, pues sumar ambas pérdidas en la ecuación activa no las vuelve equivalentes en los balances materiales del reactor.

175.8. Autocatálisis cooperativa y umbral

La transformación cooperativa formal $2X+S\to3X$ tiene velocidad de acción de masas $kSX^2$; cada evento produce una unidad neta e incorpora dos unidades preexistentes. Aquí $k$ tiene unidades de concentración inversa al cuadrado por tiempo inverso. Una ley cuadrática experimental requiere examinar el mecanismo, pues asociaciones rápidas y otras aproximaciones pueden producir exponentes efectivos semejantes.

Con sustrato constante (mantenido por un depósito externo) y pérdida lineal, la dinámica es $\dot X=kSX^2-\delta X$, donde $\delta$ tiene unidades de tiempo inverso. Los equilibrios son cero y el umbral $X_c=\delta/(kS)$; por debajo domina la pérdida y por encima la producción. El origen es estable y el umbral inestable; la semilla determina así el resultado, a diferencia de la autocatálisis bilineal cerrada sin pérdidas. Si se omite el agotamiento del sustrato, el crecimiento superumbral es ilimitado e incluso diverge en tiempo finito.

Para representar un inventario finito con $S=M-X$ hay que justificar una suma activa constante. El reciclaje explícito $X\to S$ a velocidad $\delta X$ devuelve al recurso cada unidad perdida y conduce a

En este modelo reciclado, $M$ es una concentración conservada y la ecuación escalar es exacta. Si la pérdida produce una especie inerte, esta reducción no respeta el balance. Los equilibrios positivos cumplen $kX(M-X)=\delta$; la parábola alcanza su máximo a mitad del inventario. Hay dos raíces positivas distintas si $kM^2>4\delta$, una raíz doble en la igualdad y ninguna raíz positiva si la pérdida supera el valor crítico. Las raíces son $X_\pm=(M\pm(M^2-4\delta/k)^{1/2})/2$.

El origen es localmente estable, pues allí la derivada del campo vale $-\delta$. En una raíz positiva, la derivada se reduce a $kX(M-2X)$; $X_-$ es inestable y $X_+$, estable. Las condiciones iniciales entre cero y la raíz inferior decaen al origen; si la semilla la supera y permanece en el intervalo físico, converge al equilibrio superior.

La biestabilidad permite estados inactivo y activo con iguales parámetros y distinta historia inicial. En la raíz doble aparece una bifurcación de silla-nodo cuya clasificación requiere el término no lineal; cerca del umbral puede haber transitorios lentos. Estas conclusiones corresponden al mecanismo reciclado; atribuirlas a una degradación irreversible sin cambiar el balance confunde modelos distintos. En poblaciones moleculares pequeñas, las fluctuaciones también pueden cruzar el umbral o eliminar por completo la semilla.

Si se toman el total y la constante cooperativa iguales a uno (en las unidades apropiadas) y la pérdida igual a $0{,}16$ por unidad de tiempo, los equilibrios positivos están en $0{,}2$ y $0{,}8$. Una preparación por debajo del primero decae hacia cero, mientras que una ligeramente superior crece hacia el estado activo alto. La igualdad exacta con el umbral deja la trayectoria determinista en ese punto inestable.

La cercanía al umbral puede producir un retraso prolongado antes de observar crecimiento apreciable. Una perturbación pequeña respecto de la raíz inestable crece inicialmente con la tasa lineal local, así que el tiempo necesario para alcanzar una separación macroscópica aumenta al disminuir esa perturbación. Una latencia larga no demuestra que el mecanismo catalítico se haya activado de forma tardía; puede surgir de una preparación cercana a la separatriz desde el primer instante.

175.9. Estocasticidad y semilla pequeña

Cuando la semilla contiene pocas moléculas, la concentración continua oculta la discreción de los eventos. Una descripción estocástica asigna una propensión a cada reacción y representa mediante números enteros la cantidad de unidades de sustrato y producto. La propensión es la probabilidad por unidad de tiempo de que ocurra un evento, condicionada al estado actual. Esta formulación conserva las variaciones aleatorias entre realizaciones y la posibilidad de perder todas las unidades activas, algo que una trayectoria diferencial positiva no representa exactamente.

Sean $n_S$ y $n_X$ los números de moléculas y $\Omega$ el factor que convierte concentración en número para el volumen y las unidades elegidas. Así, $S=n_S/\Omega$ y $X=n_X/\Omega$. Para la reacción bilineal, una propensión compatible con la ley macroscópica es $a_{\mathrm{b}}=(k/\Omega)n_Sn_X$. La pérdida tiene propensión $a_{\mathrm{d}}=\delta n_X$. Ambas tienen unidades de tiempo inverso; cada evento modifica los conteos según su vector estequiométrico.

El algoritmo de simulación exacta de Gillespie genera el tiempo hasta el siguiente evento y selecciona cuál ocurre mediante esas propensiones (Gillespie, 1977). Si $a_0$ es la suma de las propensiones, el tiempo de espera tiene distribución exponencial de parámetro $a_0$ y cada canal se elige con probabilidad $a_j/a_0$. Después se actualizan los conteos y se recalculan las propensiones. La exactitud se refiere al proceso de Markov especificado; su adecuación física sigue dependiendo de mezcla, mecanismo y validez de las propensiones.

El estado $n_X=0$ es absorbente cuando ninguna reacción independiente puede generar producto. En él se anulan todas las propensiones autocatalíticas, aunque quede sustrato. Si además existen pérdidas, una semilla pequeña puede desaparecer antes de amplificarse. Sin pérdidas, una semilla positiva en el modelo bilineal cerrado no se extingue por fluctuaciones, porque ningún evento reduce su conteo. La absorción en cero y la posibilidad efectiva de alcanzarla son propiedades distintas que deben comprobarse en los canales incluidos.

La discreción también cambia la finalización del proceso sin pérdidas. Con un total finito $K=n_S+n_X$, semilla $n_0\geq1$ y solo la reacción bilineal, el conteo activo aumenta de uno en uno hasta $K$. En el estado $n$, la propensión es $(k/\Omega)n(K-n)$. La suma de las esperas medias de esos pasos da el tiempo medio de conversión completa

Es finito para cada población finita. Con $K=100$, $n_0=1$, $\Omega=100\ \mathrm{L/mmol}$ y $k=1\ \mathrm{L/(mmol\,h)}$, vale aproximadamente $10{,}355$ horas. El inventario y la concentración inicial coinciden con los del ejemplo logístico, pero el resultado es la media de un tiempo aleatorio hasta el último evento. La ecuación continua se aproxima al agotamiento sin alcanzarlo en tiempo finito; las dos descripciones difieren porque la pregunta depende de resolver la última molécula.

En el modelo cooperativo reciclado con población finita y pérdida positiva hay otra diferencia. Desde cualquier estado con actividad existe una secuencia de pérdidas que conduce a cero, mientras que cero no tiene salida. El espacio de estados es finito y no hay otra clase cerrada, por lo que la extinción eventual tiene probabilidad uno. El estado activo estable de las ecuaciones deterministas describe entonces una persistencia que puede durar tiempos muy largos, pero que acaba en cada realización finita. El tiempo de observación y el tamaño del sistema determinan qué aproximación es útil. Una producción basal que recree catalizador cambiaría esta conclusión al eliminar el carácter absorbente del cero.

Ese cambio puede comprobarse experimentalmente introduciendo una producción basal controlada $S\to X$, con propensión $\eta n_S$ y $\eta>0$. En el inventario reciclado finito, el estado sin catalizador puede volver a activarse porque todavía contiene sustrato. Con pérdida positiva, todos los conteos entre cero y el total quedan comunicados por nacimientos basales y pérdidas; el proceso es irreducible y tiene una distribución estacionaria única. El tiempo de espera para salir de cero tiene media $1/(\eta K)$. Al reducir $\eta$, las reactivaciones se espacian cada vez más y, en una ventana corta, el sistema puede parecer absorbente aunque matemáticamente haya dejado de serlo. La producción basal debe medirse por separado, porque atribuir a la autocatálisis la primera aparición de actividad confundiría la iniciación con la amplificación posterior.

La media de producto mezcla trayectorias extinguidas y supervivientes. Si $N(t)=n_X(t)$, entonces $\mathbb E[N(t)]=\Pr[N(t)>0]\,\mathbb E[N(t)\mid N(t)>0]$. Una media creciente puede coexistir con una probabilidad apreciable de extinción si unas pocas realizaciones supervivientes alcanzan tamaños grandes. Además, $\mathbb E[n_Sn_X]$ no coincide en general con el producto de las medias. La ecuación determinista para concentraciones medias necesita una aproximación de cierre que desprecie correlaciones relevantes, especialmente cerca de fronteras y umbrales.

Con sustrato aproximadamente constante, la etapa inicial bilineal con muerte se aproxima a un proceso de nacimiento y muerte por individuo, con tasas $\lambda=kS$ y $\mu=\delta$. Si $\lambda>\mu$, la probabilidad de extinción eventual desde $n_0$ unidades es $(\mu/\lambda)^{n_0}$; si $\lambda\leq\mu$, es uno. El resultado se obtiene condicionando al primer evento de un linaje. Si $q$ es su probabilidad de extinción, entonces $q=\mu/(\lambda+\mu)+[\lambda/(\lambda+\mu)]q^2$, porque una muerte extingue el linaje y un nacimiento deja dos descendencias independientes. La probabilidad buscada es la menor raíz de esa ecuación en $[0,1]$. Para $n_0$ linajes iniciales independientes, la probabilidad es $q^{n_0}$. La independencia es una propiedad de esta aproximación con recurso constante; la competencia por una reserva compartida la rompe.

Este resultado requiere un depósito sostenido y no debe extrapolarse al recipiente finito. La descripción determinista es más fiable con conteos grandes, fluctuaciones relativas pequeñas y distancias suficientes de los estados absorbentes y de las separatrices dinámicas.

175.10. Redes autocatalíticas, recursos compartidos y límites conceptuales

En una red autocatalítica, unas especies facilitan la formación de otras que sostienen su regeneración. La existencia del circuito no prueba amplificación material; hay que identificar reacciones, catalizadores y balances. La estructura delimita posibilidades; velocidades y alimentación determinan cuáles se realizan, por lo que se distinguen estequiometría y dinámica (Feinberg, 1987).

Un recurso común genera competencia, porque el crecimiento de una vía reduce lo disponible para las otras. Puede favorecerse la que mantiene actividad con menor sustrato residual, aunque inhibición, costes de mantenimiento y rendimientos cambien el resultado. En alimentación cruzada, el producto secundario de una vía alimenta otra; la persistencia exige flujos compatibles y compensar las pérdidas.

El coste catalítico incluye formación, renovación y conservación funcional. Una vía rápida puede consumir recursos auxiliares o generar inhibidores; si afectan la dinámica, hay que ampliar el modelo. Cada especie requiere un balance y cada reacción respeta restricciones materiales. En sistemas abiertos, entradas y salidas fijan inventarios sostenibles; una concentración estacionaria positiva puede ocultar recambio considerable.

La termodinámica restringe la dirección y el mantenimiento de los flujos. Los catalizadores cambian velocidades y dejan inalteradas las diferencias de energía libre y el equilibrio. La actividad persistente lejos del equilibrio requiere recursos y disipación compatibles. Los modelos de organización química disipativa destacan condiciones abiertas y no linealidades en comportamientos colectivos (Prigogine y Lefever, 1968); las oscilaciones deben demostrarse para cada red.

La autocatálisis puede formar parte de sistemas vivos; para hablar de vida, inteligencia o intención hacen falta criterios y niveles organizativos adicionales. Aun una amplificación robusta puede carecer de compartimentación, herencia estable o regulación. La prueba experimental es que el producto acelere su vía con controles apropiados que descarten cambios térmicos, transporte, contaminación y producción basal independiente.

Para contrastar el mecanismo conviene variar separadamente la cantidad inicial de catalizador y el recurso total. El modelo cooperativo predice un umbral que cambia con el total y con la tasa de pérdida. El modelo logístico sin pérdida, en cambio, amplifica cualquier cantidad inicial estrictamente positiva dentro del intervalo físico. Comparar tiempos de crecimiento y destinos finales en esas preparaciones permite distinguir ambos mecanismos con más información que una sola curva creciente.

Las mediciones deben incluir el sustrato o un balance que permita reconstruirlo. Un descenso de la velocidad puede deberse al agotamiento del recurso, a la degradación del catalizador o a una limitación de transporte. Cada explicación produce respuestas distintas al añadir sustrato, modificar la dilución o cambiar la agitación. El modelo se vuelve contrastable cuando esas intervenciones corresponden a parámetros y términos definidos y cuando las cantidades conservadas se verifican con mediciones independientes.

La relación entre volumen y semilla debe declararse al comparar experimentos. Mantener fija la concentración inicial mientras aumenta el volumen incrementa el número de moléculas activas; mantener fijo ese número hace disminuir la concentración. Son preparaciones diferentes y pueden quedar a lados opuestos de un umbral determinista. En el primer caso suelen reducirse las fluctuaciones relativas cuando existe un límite macroscópico adecuado. En el segundo, la discreción inicial continúa siendo decisiva. Un protocolo reproducible informa concentraciones, cantidades y volumen, además del método de introducción de la semilla y del tiempo transcurrido desde su preparación.

El análisis cuantitativo requiere mecanismo verificable, unidades consistentes, balances cerrados y fronteras explícitas. La estabilidad determinista describe perturbaciones continuas; las simulaciones de eventos examinan establecimiento y extinción con semilla pequeña (Gillespie, 1977). La interpretación indica especies medidas, flujos omitidos y régimen de estimación de constantes; solo tras contrastar alternativas se atribuye aceleración a autocatálisis.

Referencias

Feinberg, M. (1987). Chemical reaction network structure and the stability of complex isothermal reactors—I. The deficiency zero and deficiency one theorems. Chemical Engineering Science, 42(10), 2229–2268. DOI.

Gillespie, D. T. (1977). Exact stochastic simulation of coupled chemical reactions. The Journal of Physical Chemistry, 81(25), 2340–2361. DOI.

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