A fast spectral method for electrostatics in doubly periodic slit channels

Area of Science:

• Computational physics
• Electrochemistry
• Materials science

Background:

• Calculating electrostatic interactions in dielectric media is computationally challenging.
• Existing methods struggle with spectral accuracy in systems with dielectric jumps.
• Understanding confined electrolytes requires accurate modeling of polarization effects.

Purpose of the Study:

• To develop a fast and spectrally accurate method for computing electrostatic energy and forces in doubly periodic slabs with dielectric jumps.
• To address the limitations of existing spectral Ewald methods in handling dielectric discontinuities.
• To enable accurate simulations of confined electrolytes by incorporating polarization effects.

Main Methods:

• Ewald splitting combined with a localized near-field correction for spectral accuracy.
• Recasting the problem as a two-point boundary value problem solved with a Chebyshev method.
• Utilizing a harmonic correction to preserve spectral accuracy with dielectric jumps.
• Implementation on graphical processing units (GPUs) and integration with Brownian dynamics.

Main Results:

• The developed method achieves spectral accuracy for electrostatic calculations in dielectric slabs.
• The approach effectively handles dielectric jumps using a harmonic correction.
• Simulations of binary electrolytes in confined dielectric boundaries reveal strong charge depletion near interfaces.
• The results highlight the importance of polarization effects in confined electrolytes.

Conclusions:

• The new computational method provides a fast and accurate way to study electrostatic interactions in complex dielectric environments.
• Accurate modeling of polarization effects is essential for understanding the behavior of confined electrolytes.
• The GPU-accelerated solver combined with Brownian dynamics offers a powerful tool for theoretical and computational research in this field.