Abstract:Adsorption of dimers is modelled using random sequential adsorption algorithm. The interaction between molecules is given by screened electrostatic potential. The paper focuses on the properties of adsorbed monolayers as well as the dependence of adsorption kinetics on interaction range. We designate random maximal coverage ratios, density autocorrelations and orientational ordering inside layers. Moreover the detailed analysis of adsorption kinetics is presented including a discussion of Feder's law validity and new numerical method for modelling diffusion driven adsorption. Results of numerical simulations are compared with experimental data obtained previously for insulin dimers.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to study the behavior of dimers during the surface adsorption process, especially the influence when the non - negligible electrostatic repulsion exists. Specifically, the author simulated the adsorption process of charged dimers through the Random Sequential Adsorption (RSA) algorithm and explored the following aspects:
1. **Maximum Random Coverage Ratio**: Study the influence of different electrostatic interaction ranges \(L_e\) on the maximum random coverage ratio \(\theta_{\text{max}}\) of the dimer monolayer. As \(L_e\) increases, the maximum random coverage ratio will decrease because electrostatic repulsion increases the distance between molecules, resulting in a decrease in coverage efficiency.
2. **Density Autocorrelation Function**: Analyze the spatial distribution characteristics between molecules in the adsorption layer, that is, describe the probability distribution of the distance between molecules by calculating the autocorrelation function \(G(r)\). The results show that as \(L_e\) increases, the average distance between molecules also increases, and the distribution becomes wider.
3. **Orientation Order**: Study the orientation order of dimers in the adsorption layer. Although there is no obvious orientation order globally, a slight enhancement of orientation order can be observed in the local range, especially near the boundary.
4. **Adsorption Kinetics**: Explore the variation law of the coverage ratio with time during the adsorption process, verify the validity of Feder's law, and propose a new numerical method to simulate the diffusion - driven adsorption process. The research shows that although the electrostatic interaction changes the maximum coverage ratio, the adsorption process still follows Feder's law.
5. **Comparison with Experimental Data**: Compare the simulation results with the insulin dimer adsorption data obtained in previous experiments to verify the rationality of the model. The results show that the simulation and experimental data are similar in shape, but there are differences in the maximum coverage ratio, which may be due to different experimental conditions.
### Summary of Main Formulas
- **Yukawa Potential**:
\[
U_{\text{el}}=\begin{cases}
\frac{U_0}{r}\exp\left[-\frac{(r - a)}{L_e}\right]&\text{for }r > a\\
\infty&\text{for }r < a
\end{cases}
\]
where \(r\) is the distance between the centers of two monomers, \(U_0\) characterizes the electrostatic properties of monomers and solvents, and \(L_e\) is the Debye screening length.
- **Debye Screening Length**:
\[
L_e = \sqrt{\frac{\epsilon_0\epsilon_rk_BT}{2e^2I}}
\]
where \(\epsilon_0\) is the vacuum permittivity, \(\epsilon_r\) is the dielectric constant of the solvent, \(k_B\) is the Boltzmann constant, \(T\) is the temperature, \(e\) is the elementary charge, and \(I\) is the ionic strength of the electrolyte solution.
- **Feder's Law**:
\[
\theta_{\text{max}}-\theta(t)=At^{-1/d}
\]
where \(t\) is the dimensionless time, \(d\) is the dimension of the collector (for the two - dimensional case \(d = 2\)), and \(A\) is the proportionality coefficient.
- **Effective Molecular Area**:
\[
S_{\text{eff}}(L_e)=\int d^2r\left\{\begin{array}{ll}
1-\exp\left[-\frac{U_{\text{el}}(L_e,r)}{k_BT}\right]&\text{if }\exp\left[-\frac{U_{\text{el}}(L_e,r)}{k_BT}\right]
\end{array}\right.
\]