A Generalized Model for Disordered Random Sequential Adsorption with Charge-Dependent Deposition

Technical Summary: A Generalized Model for Disordered Random Sequential Adsorption with Charge-Dependent Deposition

Problem Statement
Random Sequential Adsorption (RSA) is a fundamental nonequilibrium model describing irreversible particle deposition, exemplified by Renyi's one-dimensional car-parking problem. In standard continuum RSA, all geometrically admissible deposition coordinates are sampled uniformly, regardless of the existing configuration. However, many physical adsorption processes—such as colloidal binding, protein adsorption, and surface reactions—are influenced by surface heterogeneity, external fields, or particle-substrate interactions that bias the arrival or attachment position.

Previous extensions of RSA often introduce disorder via a pre-existing "quenched" landscape (e.g., random sites or patterned substrates) or modify trial rates and acceptance probabilities. A distinct mechanism, where the deposition measure itself evolves dynamically based on the state of the system, has been less explored. Specifically, the authors seek to model a scenario where the deposition coordinate within an available gap is determined by a conditional density dependent on the charges bounding the gap and the endpoint charges of the incoming particle. This creates a "dynamically generated disorder" where geometry and charge states evolve together as a multitype random fragmentation process.

Methodology
The authors propose a generalized one-dimensional RSA model for charged unit segments. The core innovation is the replacement of the uniform deposition kernel with a charge-dependent beta kernel.

1. The Model:

   • Consider an available gap of length $x$ bounded by charges $\alpha, \beta \in \{+, -\}$.
   • An incoming unit segment with endpoint charges $\sigma, \tau$ is deposited at position $u$ (left endpoint).
   • The conditional probability density $P_{\alpha\beta}^{\sigma\tau}(u|x)$ is defined by a normalized Beta distribution:
     $$P_{\alpha\beta}^{\sigma\tau}(u|x) \propto u^{a_{\alpha\sigma}-1}(x-1-u)^{b_{\tau\beta}-1}$$
   • The parameters $a_{\alpha\sigma} = 1 + \lambda q_\alpha q_\sigma$ and $b_{\tau\beta} = 1 + \lambda q_\tau q_\beta$ depend on the interaction parameter $0 \le \lambda < 1$ and the charge products ($q_+ = 1, q_- = -1$).
   • Physical Interpretation: $\lambda$ acts as an effective inverse temperature ($\lambda = \varepsilon/k_B T$). At $\lambda=0$, the model recovers uniform car parking. As $\lambda \to 1$, deposition becomes increasingly localized near boundaries with opposite charges (attractive) and suppressed near like charges (repulsive). The limit $\lambda \to 1$ represents an integrability threshold rather than a zero-temperature limit.

2. Analytical Framework:

   • The authors formulate a four-component probability-generating function (PGF) equation, $G_{\alpha\beta}(z, x)$, which encodes the complete occupation statistics for a gap of length $x$ with boundary charges $\alpha, \beta$.
   • This PGF satisfies an integral recursion derived from the conditional independence of the two daughter gaps created by a deposition event.
   • From the PGF, closed recursive equations are derived for the mean occupation ($f_{\alpha\beta}$) and the variance ($s_{\alpha\beta}$).

3. Numerical and Exact Solutions:

   • Exact Finite-Shell Solutions: Analytical solutions are derived for intervals $x \in [n, n+1)$ (shells). The first non-trivial shell ($2 \le x < 3$) yields exact expressions for the distribution, mean, and variance.
   • Gauss–Jacobi Quadrature: For larger systems, the authors employ Gauss–Jacobi quadrature to handle the endpoint singularities of the beta kernels efficiently, particularly as $\lambda \to 1$.
   • Monte Carlo Simulations: A rejection-free Monte Carlo algorithm is developed to sample the final jammed state directly from the conditional kernel, serving as an independent validation of the quadrature results.

Key Results

1. Exact Recursions and Benchmarks: The paper provides exact recursions for the full occupation distribution. The analytical solutions for the first few shells serve as rigorous benchmarks, confirming that the Gauss–Jacobi and Monte Carlo methods reproduce the exact mean and variance to high precision (within $10^{-4}$).
2. Asymptotic Self-Averaging: A central theoretical result is the proof that the jammed occupation density is self-averaging.
   • While initial boundary charges affect finite-size behavior, the asymptotic mean density ($\rho$) and variance density ($\chi$) are independent of the boundary states ($\alpha, \beta$).
   • This is proven using the Perron–Frobenius theorem applied to the transition matrix of the charge states. The system "forgets" the initial boundary conditions in the bulk, as repeated deposition replaces the initial charge information with the particle distribution statistics.
3. Impact of Charge Selectivity ($\lambda$) and Composition:
   • Coverage: Charge-selective deposition generally increases the jammed density compared to uniform RSA. This is because boundary-favored events (opposite charges) place segments near the edges of gaps, leaving larger usable daughter intervals.
   • Particle Mixtures: The effect is most pronounced for mixtures dominated by opposite endpoint charges ($p_{+-}, p_{-+}$). At $\lambda = 0.99$, this mixture achieves a density of $\rho \approx 0.9451$, significantly higher than the uniform RSA limit ($\approx 0.7476$) or the equal-charge dominated case ($\approx 0.8670$).
   • Fluctuations: The variance density ($\chi$) exhibits distinct behaviors. For opposite-charge dominated mixtures, high density is accompanied by strongly reduced fluctuations ($\chi \approx 0.0243$), indicating a more reproducible jammed state. Conversely, equal-charge dominated mixtures show increased fluctuations due to the competition between suppressed and favored channels.

Significance and Claims
The paper claims to provide a tractable connection between three distinct regimes:

1. Uniform RSA: Recovered at $\lambda = 0$.
2. Dynamically Generated Disorder: The model where the deposition measure evolves with the system state, distinct from quenched disorder.
3. Singular Charged-Segment Dynamics: The model connects continuously to the deterministic, charge-selective rules studied in the authors' previous work as $\lambda \to 1$.

The authors emphasize that this framework allows for the analytical study of complete occupation statistics (not just the mean) in a continuum RSA model with state-dependent deposition. The demonstration that the jammed density is self-averaging despite the complex, history-dependent charge evolution is a key theoretical contribution. The results suggest that charge composition and interaction strength can be used to control both the efficiency (packing density) and the reproducibility (fluctuations) of disordered jammed layers.

The paper does not propose specific experimental setups or future applications beyond the theoretical modeling of adsorption processes; its contribution is strictly within the domain of statistical mechanics and nonequilibrium stochastic processes.