Random sequential adsorption

Random sequential adsorption (RSA) refers to a process where particles are randomly introduced in a system, and if they do not overlap any previously adsorbed particle, they adsorb and remain fixed for the rest of the process. RSA can be carried out in computer simulation, in a mathematical analysis, or in experiments. It was first studied by one-dimensional models: the attachment of pendant groups in a polymer chain by Paul Flory, and the car-parking problem by Alfréd Rényi.^[1] Other early works include those of Benjamin Widom.^[2] In two and higher dimensions many systems have been studied by computer simulation, including in 2d, disks, randomly oriented squares and rectangles, aligned squares and rectangles, various other shapes, etc.

An important result is the maximum surface coverage, called the saturation coverage or the packing fraction. On this page we list that coverage for many systems.

The blocking process has been studied in detail in terms of the random sequential adsorption (RSA) model.^[3] The simplest RSA model related to deposition of spherical particles considers irreversible adsorption of circular disks. One disk after another is placed randomly at a surface. Once a disk is placed, it sticks at the same spot, and cannot be removed. When an attempt to deposit a disk would result in an overlap with an already deposited disk, this attempt is rejected. Within this model, the surface is initially filled rapidly, but the more one approaches saturation the slower the surface is being filled. Within the RSA model, saturation is sometimes referred to as jamming. For circular disks, saturation occurs at a coverage of 0.547. When the depositing particles are polydisperse, much higher surface coverage can be reached, since the small particles will be able to deposit into the holes in between the larger deposited particles. On the other hand, rod like particles may lead to much smaller coverage, since a few misaligned rods may block a large portion of the surface.

For the one-dimensional parking-car problem, Renyi^[1] has shown that the maximum coverage is equal to

${\displaystyle \theta _{1}=\int _{0}^{\infty }\exp \left(-2\int _{0}^{x}{\frac {1-e^{-y}}{y}}dy\right)dx=0.7475979202534\ldots }$

the so-called Renyi car-parking constant.^[4]

Then followed the conjecture of Ilona Palásti,^[5] who proposed that the coverage of d-dimensional aligned squares, cubes and hypercubes is equal to θ₁^d. This conjecture led to a great deal of work arguing in favor of it, against it, and finally computer simulations in two and three dimensions showing that it was a good approximation but not exact. The accuracy of this conjecture in higher dimensions is not known.

For ${\displaystyle k}$-mers on a one-dimensional lattice, we have for the fraction of vertices covered,^[6]

${\displaystyle \theta _{k}=k\int _{0}^{\infty }\exp \left(-u-2\sum _{j=1}^{k-1}{\frac {1-e^{-ju}}{j}}\right)du=k\int _{0}^{1}\exp \left(-2\sum _{j=1}^{k-1}{\frac {1-v^{j}}{j}}\right)dv}$

When ${\displaystyle k}$ goes to infinity, this gives the Renyi result above. For k = 2, this gives the Flory ^[7] result ${\displaystyle \theta _{1}=1-e^{-2}}$.

For percolation thresholds related to random sequentially adsorbed particles, see Percolation threshold.

system	Saturated coverage ${\displaystyle \theta _{k}}$(fraction of sites filled)
dimers	${\displaystyle 1-e^{-2}=0.86466472}$[7]
trimers	${\displaystyle {\frac {3{\sqrt {\pi }}({\text{erfi}}(2)-{\text{erfi}}(1))}{2e^{4}}}\approx 0.82365296}$[6]
k = 4	${\displaystyle 0.80389348}$[6]
k = 10	${\displaystyle 0.76957741}$[6]
k = 100	${\displaystyle 0.74976335}$[6]
k = 1000	${\displaystyle 0.74781413}$[6]
k = 10000	${\displaystyle 0.74761954}$[6]
k = 100000	${\displaystyle 0.74760008}$[6]
k = ${\displaystyle \infty }$	${\displaystyle 0.74759792}$[1]

Asymptotic behavior: ${\displaystyle \theta _{k}\sim \theta _{\infty }+0.2162/k+\ldots }$ .

R = size ratio of segments. Assume equal rates of adsorption

system	Saturated coverage ${\displaystyle \theta }$(fraction of line filled)
R = 1	0.74759792[1]
R = 1.05	0.7544753(62) [9]
R = 1.1	0.7599829(63) [9]
R = 2	0.7941038(58) [9]

system	Saturated coverage ${\displaystyle \theta _{k}}$(fraction of sites filled)
dimers k = 2	0.906820(2),[10] 0.906,[11] 0.9068,[12] 0.9062,[13] 0.906,[14] 0.905(9),[15] 0.906,[11] 0.906823(2),[16]
trimers k = 3	${\displaystyle }$[6] 0.846,[11] 0.8366 [12]
k = 4	${\displaystyle }$0.8094 [13] 0.81[11]
k = 5	0.7868 [11]
k = 6	0.7703 [11]
k = 7	0.7579 [11]
k = 8	0.7479,[13] 0.747[11]
k = 9	0.7405[11]
k = 16	0.7103,[13] 0.71[11]
k = 32	0.6892,[13] 0.689,[11] 0.6893(4)[17]
k = 48	0.6809(5),[17]
k = 64	0.6755,[13] 0.678,[11] 0.6765(6)[17]
k = 96	0.6714(5)[17]
k = 128	0.6686,[13] 0.668(9),[15] 0.668[11] 0.6682(6)[17]
k = 192	0.6655(7)[17]
k = 256	0.6628[13] 0.665,[11] 0.6637(6)[17]
k = 384	0.6634(6)[17]
k = 512	0.6618,[13] 0.6628(9)[17]
k = 1024	0.6592 [13]
k = 2048	0.6596 [13]
k = 4096	0.6575[13]
k = 8192	0.6571 [13]
k = 16384	0.6561 [13]
k = ∞	0.660(2),[17] 0.583(10),[18]

Asymptotic behavior: ${\displaystyle \theta _{k}\sim \theta _{\infty }+\ldots }$ .

system	Saturated coverage ${\displaystyle \theta _{k}}$(fraction of sites filled)
dimers k = 2	0.9142(12),[19]
k = 3	0.8364(6),[19]
k = 4	0.7892(5),[19]
k = 5	0.7584(6),[19]
k = 6	0.7371(7),[19]
k = 8	0.7091(6),[19]
k = 10	0.6912(6),[19]
k = 12	0.6786(6),[19]
k = 20	0.6515(6),[19]
k = 30	0.6362(6),[19]
k = 40	0.6276(6),[19]
k = 50	0.6220(7),[19]
k = 60	0.6183(6),[19]
k = 70	0.6153(6),[19]
k = 80	0.6129(7),[19]
k = 90	0.6108(7),[19]
k = 100	0.6090(8),[19]
k = 128	0.6060(13),[19]

system	Saturated coverage ${\displaystyle \theta _{k}}$(fraction of sites filled)
Square lattice with NN exclusion	0.3641323(1),[20] 0.36413(1),[21] 0.3641330(5),[22]
Honeycomb lattice with NN exclusion	0.37913944(1),[20] 0.38(1),[2] 0.379[23]

.

system	Saturated coverage ${\displaystyle \theta _{k}}$(fraction of sites filled)
k = 2	0.74793(1),[24] 0.747943(37),[25] 0.749(1),[26]
k = 3	0.67961(1),[24] 0.681(1),[26]
k = 4	0.64793(1),[24] 0.647927(22)[25] 0.646(1),[26]
k = 5	0.62968(1)[24] 0.628(1),[26]
k = 8	0.603355(55)[25] 0.603(1),[26]
k = 10	0.59476(4)[24] 0.593(1),[26]
k = 15	0.583(1),[26]
k = 16	0.582233(39)[25]
k = 20	0.57807(5)[24] 0.578(1),[26]
k = 30	0.574(1),[26]
k = 32	0.571916(27)[25]
k = 50	0.56841(10)[24]
k = 64	0.567077(40)[25]
k = 100	0.56516(10)[24]
k = 128	0.564405(51)[25]
k = 256	0.563074(52)[25]
k = 512	0.562647(31)[25]
k = 1024	0.562346(33)[25]
k = 4096	0.562127(33)[25]
k = 16384	0.562038(33)[25]

For k = ∞, see "2d aligned squares" below. Asymptotic behavior:^[25] ${\displaystyle \theta _{k}\sim \theta _{\infty }+0.316/k+0.114/k^{2}\ldots }$ . See also ^[27]

• Adsorption
• Particle deposition
• Percolation threshold