An Analytic Approach To Maximize Entropy For Computing Equilibrium Densities Of K-Mers On Linear Chains

The irreversible adsorption of polyatomic (or k-mers) on linear chains is related to phenomena such as the adsorption of colloids, long molecules, and proteins on solid substrates. This process generates jammed or blocked final states. In the case of k = 2, the binomial coefficient computes the number of final states. By the canonical ensemble, the Boltzmann-Gibbs-Shannon entropy function is obtained by using Stirling's approximation, and its equilibrium density rho(eq,2) is its maximum at the thermodynamic limit with value rho(eq,2) approximate to 0.822 991 17. Moreover, since at the same energy we have several possible configurations, we obtain the state probability density. Maximizing the entropy, it converges to a Gaussian distribution N(rho(eq,2), sigma(2)(eq,2)) as L -> infinity L -> infinity.In this article, we generalize this analysis to k > 2 to maximize the entropy and to get the equilibrium densities rho(eq,k). We first develop a complete combinatorial analysis to get the generalized recurrence formula (GRF) for counting all blocked configuration states on a chain of length L with fixed k, which corresponds to a generalized truncated Fibonacci sequence. The configuration states for allocating N k-mers is related with the general binomial coefficient[GRAPHICS]. Since Stirling's approximation cannot be used for GRF, we numerically compute the state probability density and approximate rho(eq,k) and sigma(eq,k) for large k-mers with high precision for k-mers up to k = 1, 000, 000. We highlight that rho(eq,k) decreases from k = 2, horizontal ellipsis , 8 reaching a minimum at k = 9 and then increases with an asymptotic value rho(eq,infinity) = 0.9285685. We compared with jamming densities obtained by RSA and at k approximate to 16, both curves intersect and ergodicity is not broken since rho(jam,k) approximate to rho(eq,k). In the case of sigma(2)(eq,k), it grows similarly with asymptotic value sigma(2)(eq,infinity) = 0.8629597. Since the similar behavior for large values, we found the limit relationship sigma(2)(eq,k) = rho(2)(eq,k) L -> infinity for any k. Finally, as k -> infinity, we get the Gaussian distribution for the continuous blocked irreversible adsorption or equivalent to the irreversible blocked car parking problem.