Emily Krucker-Velasquez,James W. Swan,Zachary Sherman
Abstract:We develop a computational method for modeling electrostatic interactions of arbitrarily-shaped, polarizable objects on colloidal length scales, including colloids/nanoparticles, polymers, and surfactants, dispersed in explicit ion electrolytes and nonionic solvents. Our method computes the nonuniform polarization charge distribution induced in a colloidal particle by both externally applied electric fields and local electric fields arising from other charged objects in the dispersion. This leads to expressions for electrostatic energies, forces, and torques that enable efficient molecular dynamics and Brownian dynamics simulations of colloidal dispersions in electrolytes, which can be harnessed to accurately predict structural and transport properties. We describe an implementation in which colloidal particles are modeled as rigid composites of small spherical beads that tessellate the surface of the particle. The electrostatics calculations are accelerated using a spectrally-accurate particle-mesh-Ewald technique implemented on a graphics processing unit (GPU) and regularized such that the electrostatic calculations are well-defined even for overlapping bodies. We demonstrate the effectiveness of this approach through a series of calculations, including the induced dipole moments and forces for one, two, and lattices of spherical colloids in an electric field; the induced dipole moment and torque for anisotropic particles in an electric field; the equilibrium distribution of ions in the double layer around charged colloids; the dynamics of charged colloids; and ions in the double layer around a polarizable colloid exposed to an electric field.
What problem does this paper attempt to address?
The problem that this paper aims to solve is to develop a computational method for simulating the electrostatics interactions of polarized colloidal particles of arbitrary shapes dispersed in explicit - ion electrolytes and non - ionic solvents. These particles include colloids / nanoparticles, polymers and surfactants. Specifically, the method proposed in the paper can calculate the non - uniform polarization charge distribution inside the colloidal particles induced by the externally applied electric field and the local electric field generated by other charged objects in the dispersion system. This leads to expressions for electrostatics energy, force and torque, so that molecular dynamics and Brownian dynamics simulations of colloidal dispersion systems in electrolytes can be carried out efficiently, and then the structural and transport properties can be accurately predicted.
### Main Problems and Solutions
**Problems**:
1. **Multibody electrostatics interactions**: In the colloidal dispersion system, the electrostatics interactions between multiple charged objects are complex and difficult to calculate accurately.
2. **Polarization effects**: Colloidal particles will produce non - uniform surface charge distributions under the action of external electric fields and local electric fields, which have an important impact on the interactions between particles.
3. **Computational efficiency**: Traditional computational methods have high computational costs when dealing with a large number of particles and complex geometries, and it is difficult to achieve efficient dynamic simulations.
**Solutions**:
1. **Immersed Boundary Method**: The paper proposes a new computational method, using the Immersed Boundary Method to simulate polarized colloidal particles of arbitrary shapes. This method simplifies the calculation by modeling the colloidal particles as rigid complexes composed of small spherical beads, which are inlaid on the particle surface.
2. **Spectrally - accurate Particle - Mesh Ewald technique**: In order to accelerate electrostatics calculations, the paper adopts the spectrally - accurate Particle - Mesh Ewald technique based on Graphics Processing Unit (GPU). This method not only improves the calculation speed, but also ensures the accuracy of the calculation, and can be well - defined even in the case of overlapping bodies.
3. **Verification and application**: The paper verifies the effectiveness of this method through a series of calculations, including the induced dipole moments and forces of single, double and lattice - arranged spherical colloids in an electric field, the induced dipole moments and torques of anisotropic particles in an electric field, the equilibrium distribution of ions in the electric double layer, and the kinetic behaviors of charged colloids and ions in the electric double layer.
### Formulas and Methods
1. **Electric potential equation**:
\[
\nabla\cdot\mathbf{E}=\frac{\rho}{\epsilon}, \quad \mathbf{E}=-\nabla\psi
\]
where \(\mathbf{E}\) is the electric field, \(\rho\) is the free charge density, and \(\epsilon\) is the dielectric constant.
2. **Laplace equation**:
\[
\nabla^{2}\psi = 0
\]
In the form of boundary conditions:
\[
\psi_{p}=\psi_{f}, \quad (\epsilon_{f}\mathbf{E}_{f}-\epsilon_{p}\mathbf{E}_{p})\cdot\hat{n}=\frac{q_{\alpha}}{4\pi a^{2}}
\]
where \(\psi_{p}\) and \(\psi_{f}\) are the electric potentials inside and outside the particle respectively, \(q_{\alpha}\) is the net charge of ion \(\alpha\), \(\hat{n}\) is the outer normal vector, and \(\mathbf{E}_{f}\) and \(\mathbf{E}_{p}\) are the electric fields outside and inside the particle respectively.
3. **Green's function**:
\[
G(\mathbf{x}-\mathbf{y})=\frac{1}{V}\sum_{\mathbf{k}\neq0}\frac{e^{i\mathbf{k}\cdot(\mathbf{x}-\mathbf{y})}}{k^{2}}\sum_{\ell p}\frac{4\pi(2\ell + 1)i^{-\ell}}{k^{2}}j