Abstract:The Poisson-Boltzmann equation is a nonlinear elliptic equation that describes how electronic potential changes in the screening layer when a cell is immersed in an ionic solution. It is the equilibrium state formed by many charged particles interacting with each other through Coulomb forces. We propose a direct simulation approach for the dynamics of the charged particles in the solution governed by the Poisson-Nernst-Planck equation and the Poisson-Boltzmann equation can be solved automatically when the particles reach the equilibrium. The interacting $N$ particle system is simulated in the truncated external domain. Directly simulating the $N$ particle system has an $\mathcal{O}(N^2)$ computational cost in each time step. By making use of the Random Batch Method (RBM) proposed by Jin et al [J. Comput. Phys., 400(1): 108877, 2020], the cost per time step is reduced remarkably to $\mathcal{O}(N)$. Another benefit of RBM is that the hard sphere potential is no longer needed in the simulation. An error analysis is provided for the truncation and numerical experiments are performed to validate the method. This particle method is preferable in two aspects: it is simple and effective in high dimensions, especially in non-symmetric cases or with complicated membrane or molecular surface. Meanwhile, direct simulation may be preferred as people may be interested in the dynamics of the physical process.

Mean field error estimate of the random batch method for large interacting particle system

Mean field error estimate of the random batch method for large interacting particle system

On the Mean Field Limit of the Random Batch Method for Interacting Particle Systems

On Random Batch Methods (RBM) for interacting particle systems driven by Lévy processes

Ergodicity and long-time behavior of the Random Batch Method for interacting particle systems

Convergence of Random Batch Method with replacement for interacting particle systems

Uniform Error Estimates for the Random Batch Method to the First‐order Consensus Models with Antisymmetric Interaction Kernels

Reduced variance random batch methods for nonlocal PDEs

Convergence of Random Batch Method for Interacting Particles with Disparate Species and Weights

Uniform-in-time Error Estimate of the Random Batch Method for the Cucker–Smale Model

Random Batch Methods for classical and quantum interacting particle systems and statistical samplings

Convergence of Random Batch Method for interacting particles with disparate species and weights

On the Random Batch Method for second order interacting particle systems

Random Batch Method with Momentum Correction

On the mean-field limit of the Cucker-Smale model with Random Batch Method

Some remarks on the effect of the Random Batch Method on phase transition

A random batch Ewald method for charged particles in the isothermal-isobaric ensemble

Convergence analysis of an explicit method and its random batch approximation for the McKean-Vlasov equations with non-globally Lipschitz conditions

A Direct Simulation Approach for the Poisson-Boltzmann Equation Using the Random Batch Method.

A Random-Batch Monte Carlo Method for Many-Body Systems with Singular Kernels.

A collision-oriented interacting particle system for Landau-type equations and the molecular chaos