Stokes Law at Molecular Length Scales: Effects of Intermolecular Interactions and Linear Response Theory

Subhajit Acharya,Biman Bagchi

2024-02-26

Abstract:The celebrated Stokes Law (SL) of hydrodynamics predicts that the velocity of a particle pulled through a liquid by an external force, Fex, is directly proportional to the force and inversely proportional to the friction {\zeta} acted by the medium on the particle. We investigate the range of validity of Stokes Law at molecular length scales by employing computer simulations to calculate friction by pulling a tagged particle with a constant force. We thus calculate friction for two model interaction potentials, Lennard-Jones, and soft sphere, for several particle sizes, ranging from radius (a) smaller than the solvent particles to three times larger. We next obtain friction from diffusion (D) by using Einstein's relation between diffusion and friction {\zeta} in an unperturbed liquid. We find, to our surprise, a quantitative agreement between the two at a small-to-intermediate pulling force regime for all the sizes studied. The Law does break down at a large pulling force beyond a threshold value. Importantly, the range of validity of Stokes' scheme to obtain friction increases substantially if we turn off the attractive part of the interaction potential. Additionally, we calculate the viscosity ({\eta}) of the unperturbed liquid and find a good agreement with the Stokes-Einstein relation {\zeta}=C{\eta}a for the viscosity dependence with a value of C close to 5 {\pi}, that is intermediate between the slip and stick boundary condition.

Soft Condensed Matter,Statistical Mechanics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to explore the validity and applicable range of Stokes' law at the molecular scale. Specifically, the author studied through computer simulations under different conditions (such as different inter - particle interaction potentials, particle sizes and thermodynamic states), when a constant external force is applied to a tagged particle, when the linear relationship of Stokes' law will fail.

### Main problems include:

1. **When does the linear relationship fail?**
- At what value of the external force does the linear relationship (i.e., the velocity is proportional to the external force) start to fail? Is there a threshold value, beyond which the linear relationship no longer holds?

2. **Sensitivity of the threshold force to the external force intensity:**
- How is the threshold force affected by the interaction potential between the tagged particle and the solvent molecules? In particular, will the threshold force change after the attractive part is removed?

3. **Influence of the thermodynamic state:**
- How does the thermodynamic state of the system (such as temperature and density) affect the threshold force? Does the linear region shrink as the density of the system increases?

4. **Understanding of the microscopic mechanism:**
- Can the earlier failure of the linear relationship with the enhancement of the intermolecular interaction be explained by the microscopic mechanism?

5. **Influence of particle size:**
- How does the size of the tagged particle affect the effective range of the linear relationship?

### Research methods:

To answer these questions, the author carried out "theoretical experiments" using molecular dynamics simulations to obtain the drift velocity of the tagged particle under the action of an external force. Two different interaction potentials were considered in the study: Lennard - Jones potential and soft - sphere potential, and simulations were carried out in two different thermodynamic states. In addition, the influence of tagged particles of different sizes on the results was also studied.

### Key findings:

- **Threshold of the linear relationship**: It was found that for the Lennard - Jones system, the threshold force is lower, while for the soft - sphere system, the threshold force is higher.
- **Influence of the thermodynamic state**: At a higher density, the linear relationship fails earlier, which may be due to the cage effect in high - density liquids.
- **Microscopic mechanism**: The presence of the attractive part makes the velocity vectors of the solvent molecules more easily aligned in the stretching direction, resulting in an earlier non - linear response.
- **Consistency of the two methods**: The friction coefficients calculated by Stokes' law and Einstein's relation are highly consistent numerically, which verifies the validity of the two methods.

In conclusion, through detailed simulation studies, this paper explores the applicability and limitations of Stokes' law at the molecular scale and reveals the key factors affecting its validity.