Source author record

Zecheng Gan

Zecheng Gan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2026arXiv

An O(N) quasi-Ewald splitting method for nanoconfined electrostatics

Simulating the dynamics of charged particles in quasi-two-dimensional (quasi-2D) nanoconfined systems presents a significant computational challenge due to the long-range nature of electrostatic interactions and the geometric anisotropy. To address this, we introduce a novel quasi-Ewald splitting strategy tailored for particle-based simulations in such geometry. Our splitting strategy seamlessly integrates a collection of advanced numerical techniques, including optimal quadrature rules [L. N. Trefethen, SIAM Rev. 64(1)(2022), pp.132-150], fast pairwise kernel summation methods [S. Jiang and L. Greengard, Commun. Comput. Phys. 31(1)(2022), pp.1-26], and the random batch method with importance sampling in k-space [S. Jin, L. Li, Z. Xu et al., SIAM J. Sci. Comput. 43(4)(2021), pp.B937-B960]. The resulting algorithm achieves an O(N) overall computational complexity, where N denotes the total number of confined particles. Simulations of several prototype systems validate the accuracy and efficiency of our method. Furthermore, we present numerical observations specifically related to nanoconfined charged many-body systems, highlighting phenomena such as dielectric boundary effects, anisotropic diffusion, and the structure of the electrical double layer (EDL) under conditions of charge asymmetry.

preprint2020arXiv

Harmonic surface mapping algorithm for electrostatic potentials in an atomistic/continuum hybrid model for electrolyte solutions

Simulating charged many-body systems has been a computational demanding task due to the long-range nature of electrostatic interaction. For the multi-scale model of electrolytes which combines the strengths of atomistic/continuum electrolyte representations, a harmonic surface mapping algorithm is developed for fast and accurate evaluation of the electrostatic reaction potentials. Our method reformulates the reaction potential into a sum of image charges for the near-field, and a charge density on an auxiliary spherical surface for the far-field, which can be further discretized into point charges. Fast multipole method is used to accelerate the pairwise Coulomb summation. The accuracy and efficiency of our algorithm, as well as the choice of relevant numerical parameters are demonstrated in detail. As a concrete example, for charges close to the dielectric interface, our method can improve the accuracy by two orders of magnitudes compared to the Kirkwood series expansion method.

preprint2013arXiv

Efficient Implementation of the Barnes-Hut Octree Algorithm for Monte Carlo Simulations of Charged Systems

Computer simulation with Monte Carlo is an important tool to investigate the function and equilibrium properties of many systems with biological and soft matter materials solvable in solvents. The appropriate treatment of long-range electrostatic interaction is essential for these charged systems, but remains a challenging problem for large-scale simulations. We have developed an efficient Barnes-Hut treecode algorithm for electrostatic evaluation in Monte Carlo simulations of Coulomb many-body systems. The algorithm is based on a divide-and-conquer strategy and fast update of the octree data structure in each trial move through a local adjustment procedure. We test the accuracy of the tree algorithm, and use it to perform computer simulations of electric double layer near a spherical interface. It has been shown that the computational cost of the Monte Carlo method with treecode acceleration scales as $\log N$ in each move. For a typical system with ten thousand particles, by using the new algorithm, the speed has been improved by two orders of magnitude from the direct summation.

preprint2012arXiv

Effects of image charges, interfacial charge discreteness, and surface roughness on the zeta potential of spherical electric double layers

We investigate the effects of image charges, interfacial charge discreteness, and surface roughness on spherical electric double layers in electrolyte solutions with divalent counter-ions in the setting of the primitive model. By using Monte Carlo simulations and the image charge method, the zeta potential profile and the integrated charge distribution function are computed for varying surface charge strengths and salt concentrations. Systematic comparisons were carried out between three distinct models for interfacial charges: 1) SURF1 with uniform surface charges, 2) SURF2 with discrete point charges on the interface, and 3) SURF3 with discrete interfacial charges and finite excluded volume. By comparing the integrated charge distribution function (ICDF) and potential profile, we argue that the potential at the distance of one ion diameter from the macroion surface is a suitable location to define the zeta potential. In SURF2 model, we find that image charge effects strongly enhance charge inversion for monovalent interfacial charges, and strongly suppress charge inversion for multivalent interfacial charges. For SURF3, the image charge effect becomes much smaller. Finally, with image charges in action, we find that excluded volumes (in SURF3) suppress charge inversion for monovalent interfacial charges and enhance charge inversion for multivalent interfacial charges. Overall, our results demonstrate that all these aspects, i.e., image charges, interfacial charge discreteness, their excluding volumes have significant impacts on the zeta potential, and thus the structure of electric double layers.