Abstract:We investigate the phenomenon of counterion condensation in a solution of highly charged rigid polyelectrolytes within the cell model. A method is proposed which -- based on the charge distribution function -- identifies both the fraction of condensed ions and the radial extension of the condensed layer. Within salt-free Poisson-Boltzmann (PB) theory it reproduces the well known fraction 1-1/xi of condensed ions for a Manning parameter xi>1. Furthermore, it predicts a weak salt dependence of this fraction and a breakdown of the concept of counterion condensation in the high salt limit. We complement our theoretical investigations with molecular dynamics simulations of a cell-like model, which constantly yield a stronger condensation than predicted by PB theory. While the agreement between theory and simulation is excellent in the monovalent, weakly charged case, it deteriorates with increasing electrostatic interaction strength and, in particular, increasing valence. For instance, at a high concentration of divalent salt and large xi our computer simulations predict charge oscillations, which mean-field theory is unable to reproduce.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the research on the counter - ion condensation phenomenon in charged rigid polyelectrolyte solutions. Specifically, the authors explored the condensation of counter - ions around highly charged rigid polyelectrolytes in the cell model. The research mainly focused on the following aspects:
1. **Proportion of counter - ion condensation**: The authors proposed a method based on the charge distribution function to identify the proportion of counter - ions in the condensed layer and its radial expansion. In the framework of the salt - free Poisson - Boltzmann (PB) theory, this method can reproduce the known proportion of condensed counter - ions \(1-\frac{1}{\xi}\), where \(\xi\) is the Manning parameter, representing the ratio of the line charge density to the Bjerrum length.
2. **Effect of salt**: The research further explored the effect of salt on the counter - ion condensation phenomenon. The authors found that as the salt concentration increases, the concept of counter - ion condensation gradually fails, especially at high salt concentrations, the structure of the condensed layer is destroyed.
3. **Comparison between theory and simulation**: In order to verify the theoretical predictions, the authors carried out molecular dynamics simulations. The simulation results show that in the case of monovalent and weakly charged, the theory and simulation results match well; however, as the strength of the electrostatic interaction increases, especially in the presence of high - valent ions, the difference between the theoretical prediction and the simulation results increases significantly.
4. **New condensation criteria**: The authors proposed a new method to define and identify the condensed layer, that is, by analyzing the inflection point in the logarithmic radial distance graph of the counter - ion distribution function \(P(r)\) to determine the boundary of the condensed layer and the proportion of condensed counter - ions. This method is not only applicable to the PB theory under salt - free conditions, but can also be extended to high - salt concentration situations, thus more accurately quantifying the degree of counter - ion condensation.
In summary, the main objective of this paper is to gain an in - depth understanding of the counter - ion condensation behavior in charged rigid polyelectrolyte solutions through theoretical analysis and computer simulations, and to explore the effect of salt on this process.