A fast spectral method for electrostatics in doubly periodic slit channels

Broken symmetries in quasi-2D charged systems via negative dielectric confinement

We report spontaneous symmetry breaking (SSB) phenomena in symmetrically charged binary particle systems under planar nanoconfinement with negative dielectric constants. The SSB is triggered solely via the dielectric confinement effect, without any external fields. The mechanism of SSB is found to be caused by the strong polarization field enhanced by nanoconfinement, giving rise to charge/field oscillations in the transverse directions. Interestingly, dielectric contrast can even determine the degree of SSB in transverse and longitudinal dimensions, forming charge-separated interfacial liquids and clusters on square lattices. Furthermore, we analytically show that the formed lattice constant is determined by the dielectric mismatch and the length scale of confinement, which is validated via molecular dynamics simulations. The novel broken symmetry mechanism may provide new insights into the study of quasi-2D systems and the design of future nanodevices.

Computing hydrodynamic interactions in confined doubly periodic geometries in linear time

We develop a linearly scaling variant of the force coupling method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly periodic geometry with either a single bottom wall or two walls (slit channel) in the aperiodic direction. Our spectrally accurate Stokes solver uses the fast Fourier transform in the periodic xy plane and Chebyshev polynomials in the aperiodic z direction normal to the wall(s). We decompose the problem into two problems. The first is a doubly periodic subproblem in the presence of particles (source terms) with free-space boundary conditions in the z direction, which we solve by borrowing ideas from a recent method for rapid evaluation of electrostatic interactions in doubly periodic geometries [Maxian et al., J. Chem. Phys. 154, 204107 (2021)]. The second is a correction subproblem to impose the boundary conditions on the wall(s). Instead of the traditional Gaussian kernel, we use the exponential of a semicircle kernel to model the source terms (body force) due to the presence of particles and provide optimum values for the kernel parameters that ensure a given hydrodynamic radius with at least two digits of accuracy and rotational and translational invariance. The computation time of our solver, which is implemented in graphical processing units, scales linearly with the number of particles, and allows computations with about a million particles in less than a second for a sedimented layer of colloidal microrollers. We find that in a slit channel, a driven dense suspension of microrollers maintains the same two-layer structure as above a single wall, but moves at a substantially lower collective speed due to increased confinement.

Numerical simulations of confined Brownian-yet-non-Gaussian motion

Brownian motion is a central scientific paradigm. Recently, due to increasing efforts and interests towards miniaturization and small-scale physics or biology, the effects of confinement on such a motion have become a key topic of investigation. Essentially, when confined near a wall, a particle moves much slower than in the bulk due to friction at the boundaries. The mobility is therefore locally hindered and space-dependent, which in turn leads to the apparition of so-called multiplicative noises, and associated non-Gaussianities which remain difficult to resolve at all times. Here, we exploit simple, optimized and efficient numerical simulations to address Brownian motion in confinement in a broadrange and quantitative way. To do so, we integrate the overdamped Langevin equation governing the thermal dynamics of a negatively-buoyant single spherical colloid within a viscous fluid confined by two rigid walls, including surface charges. From the produced large set of long random trajectories, we perform a complete statistical analysis and extract all the key quantities, such as the probability distributions in displacements and their main moments. In particular, we propose a novel method to compute high-order cumulants by reducing convergence problems, and employ it to efficiently characterize the inherent non-Gaussianity of the confined process.

Frequency and field-dependent response of confined electrolytes from Brownian dynamics simulations

Using Brownian dynamics simulations, we investigate the effects of confinement, adsorption on surfaces, and ion-ion interactions on the response of confined electrolyte solutions to oscillating electric fields in the direction perpendicular to the confining walls. Nonequilibrium simulations allows to characterize the transitions between linear and nonlinear regimes when varying the magnitude and frequency of the applied field, but the linear response, characterized by the frequency-dependent conductivity, is more efficiently predicted from the equilibrium current fluctuations. To that end, we (rederive and) use the Green-Kubo relation appropriate for overdamped dynamics, which differs from the standard one for Newtonian or underdamped Langevin dynamics. This expression highlights the contributions of the underlying Brownian fluctuations and of the interactions of the particles between them and with external potentials. Although already known in the literature, this relation has rarely been used to date, beyond the static limit to determine the effective diffusion coefficient or the DC conductivity. The frequency-dependent conductivity always decays from a bulk-like behavior at high frequency to a vanishing conductivity at low frequency due to the confinement of the charge carriers by the walls. We discuss the characteristic features of the crossover between the two regimes, most importantly how the crossover frequency depends on the confining distance and the salt concentration, and the fact that adsorption on the walls may lead to significant changes both at high and low frequencies. Conversely, our results illustrate the possibility to obtain information on diffusion between walls, charge relaxation, and adsorption by analyzing the frequency-dependent conductivity.

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a "dry diffusion" approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

Superscalability of the random batch Ewald method

Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular dynamics simulations, which are major tools for problems at the nano-/micro-scale. Various algorithms are developed to calculate the pairwise Coulomb interactions to a linear scale, but poor scalability limits the size of simulated systems. Here, we use an efficient molecular dynamics algorithm with the random batch Ewald method on all-atom systems where the complete Fourier components in the Coulomb interaction are replaced by randomly selected mini-batches. By simulating the N-body systems up to 10⁸ particles using 10 000 central processing unit cores, we show that this algorithm furnishes O(N) complexity, almost perfect scalability, and an order of magnitude faster computational speed when compared to the existing state-of-the-art algorithms. Further examinations of our algorithm on distinct systems, including pure water, a micro-phase separated electrolyte, and a protein solution, demonstrate that the spatiotemporal information on all time and length scales investigated and thermodynamic quantities derived from our algorithm are in perfect agreement with those obtained from the existing algorithms. Therefore, our algorithm provides a promising solution on scalability of computing the Coulomb interaction. It is particularly useful and cost-effective to simulate ultra-large systems, which is either impossible or very costly to conduct using existing algorithms, and thus will be beneficial to a broad range of problems at nano-/micro-scales.

Conditions for electroneutrality breakdown in nanopores

It has recently been suggested that a breakdown of electroneutrality occurs in highly confined nanopores that are encompassed by a dielectric material. This work elucidates the conditions for this breakdown. We show that the breakdown within the pore results from the response of the electric field within the dielectric. Namely, we show that this response is highly sensitive to the boundary condition at the dielectric edge. The standard Neumann boundary condition of no-flux predicts that the breakdown does not occur. However, a Dirichlet boundary condition for a zero-potential predicts a breakdown. Within this latter scenario, the breakdown exhibits a dependence on the thickness of the dielectric material. Specifically, infinite thickness dielectrics do not exhibit a breakdown, while dielectrics of finite thickness do exhibit a breakdown. Numerical simulations confirm theoretical predictions. The breakdown outcomes are discussed with regard to single pore systems and multiple pore systems.