Complete local expansion of the availability function in random sequential adsorption of aligned squares at low density: Termination at fourth order

Resumen Técnico: Expansión Local Completa de la Función de Disponibilidad en la Adsorción Secuencial Aleatoria de Cuadrados Alineados

Planteamiento del Problema
El artículo aborda la adsorción secuencial aleatoria (RSA, por sus siglas en inglés) de cuadrados alineados en una superficie bidimensional. Específicamente, se centra en el régimen de baja cobertura para derivar una expansión exacta y de forma cerrada de la función de disponibilidad, $\alpha(q)$. Esta función representa la fracción de posiciones en el plano donde se puede insertar un cuadrado adicional sin que se solape con los cuadrados depositados previamente, dada una cobertura superficial fraccional $q = \rho L^2$ (donde $\rho$ es la densidad numérica y $L$ es la longitud del lado). Mientras que el problema de la RSA unidimensional (el problema de estacionamiento de Rényi) es exactamente resoluble, los sistemas de mayor dimensión exhiben complejas correlaciones geométricas que dificultan las soluciones analíticas. Trabajos previos han proporcionado estimaciones numéricas de los coeficientes en la expansión de baja densidad de $\alpha(q)$ para diversas formas, pero han faltado expresiones de forma cerrada explícitas derivadas directamente de configuraciones de solapamiento geométrico más allá del orden principal.

Metodología
Los autores emplean un principio de inclusión-exclusión sistemático aplicado a la geometría de las zonas de exclusión.

1. Marco Geométrico: Cada cuadrado depositado de lado $L$ crea una zona de exclusión de lado $2L$ (área $4L^2$) donde el centro (o la esquina de referencia) de un nuevo cuadrado no puede colocarse. La función de disponibilidad se define como $1$ menos la fracción del plano cubierta por la unión de estas zonas de exclusión.
2. Lógica de la Expansión de Clústeres: La expansión de $\alpha(q)$ se construye corrigiendo los errores de sobreconteo inherentes a la suma de zonas de exclusión independientes:
   • Primer Orden: Resta el área de las zonas de exclusión individuales.
   • Segundo Orden: Suma de nuevo el área donde dos zonas de exclusión se solapan (corrigiendo la doble resta).
   • Tercer Orden: Resta de nuevo el área donde tres zonas de exclusión se solapan (corrigiendo la sobre-suma).
   • Cuarto Orden: Suma de nuevo el área donde cuatro zonas de exclusión se solapan.
3. Condición de Terminación: Una observación geométrica central impulsa la metodología: para cuadrados alineados, no más de cuatro cuadrados mutuamente no solapados pueden solapar simultáneamente la región de exclusión de un cuadrado de prueba. En consecuencia, la jerarquía de inclusión-exclusión para solapamientos locales termina exactamente en el cuarto orden.
4. Cálculo: Los coeficientes se derivan enumerando e integrando explícitamente todas las configuraciones de solapamiento de $n$-cuerpos ($n=1$ a $4$) geométricamente admisibles. Las integrales se factorizan en componentes unidimensionales debido a la naturaleza de los cuadrados alineados con los ejes. Se aplican factores de simetría combinatoria ($1/n!$) para dar cuenta de los clústeres ordenados frente a los no ordenados.
5. Validación: Los coeficientes analíticos se confirman mediante simulaciones numéricas en una red discreta, donde los coeficientes corresponden a conteos normalizados de clústeres de bloqueadores compatibles.

Contribuciones Clave y Resultados
La contribución principal es la derivación del polinomio completo y exacto de baja densidad para la función de disponibilidad de cuadrados alineados hasta su orden máximo geométricamente permitido. La expansión resultante es:

$$\alpha(q) = 1 - 4q + \frac{7}{2}q^2 - \frac{8}{9}q^3 + \frac{17}{288}q^4$$

Los coeficientes específicos son:

• $a_1 = -4$: Representa el área de exclusión de un solo cuadrado ($4L^2$).
• $a_2 = 7/2$: La primera corrección colectiva, que representa el solapamiento promedio de dos zonas de exclusión.
• $a_3 = -8/9$: La corrección para solapamientos triples.
• $a_4 = 17/288$: La corrección final para solapamientos cuádruples.

El artículo demuestra que este polinomio no es un truncamiento arbitrario sino la jerarquía de solapamiento local completa. Los autores calculan la raíz de este polinomio cúrtico, encontrando un punto de saturación de $q^* \approx 0.3455$.

Significado y Reivindicaciones
El artículo reclama significancia en las siguientes áreas:

1. Exactitud de la Expansión Local: Proporciona las primeras expresiones explícitas de forma cerrada para los coeficientes de la función de disponibilidad hasta el cuarto orden para cuadrados alineados, derivadas directamente de núcleos de solapamiento geométrico en lugar de ajustes numéricos.
2. Definición de la Física Local vs. No Local: El trabajo establece un límite claro entre las restricciones geométricas locales y los fenómenos de congestión (jamming) globales. El polinomio cúrtico describe con precisión el sistema hasta el punto en que las correlaciones locales se agotan ($q \approx 0.3455$). Sin embargo, el límite de congestión real para la RSA de cuadrados alineados es conocido por ser aproximadamente $0.5620$.
3. Limitaciones de los Modelos Locales: La discrepancia entre la raíz cúrtica ($0.3455$) y la verdadera densidad de congestión ($0.5620$) demuestra que el acercamiento a la congestión está controlado por correlaciones colectivas no locales generadas por la historia de deposición irreversible, las cuales no pueden ser capturadas por ninguna expansión local de orden finito.
4. Referencia Analítica: Los coeficientes derivados sirven como referencias analíticas rigurosas para probar simulaciones numéricas de RSA y teorías aproximadas de adsorción irreversible.

Los autores enfatizan que, si bien la expansión local es exacta hasta el cuarto orden, el fracaso de este polinomio para predecir el verdadero límite de congestión resalta la complejidad intrínseca de los sistemas de RSA de mayor dimensión, donde las correlaciones globales dependientes de la historia dominan el proceso de saturación.