BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zhenli Xu

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

ResearcherUnclaimed source record
physics.comp-phcond-mat.softphysics.chem-phcond-mat.stat-mechmath.NANumerical AnalysisComputational Engineering, Finance, and Science

Catalog footprint

What is connected

19works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

19 published item(s)

preprint2026arXiv

Random Batch Sum-of-Gaussians Method for Molecular Dynamics of Born-Mayer-Huggins Systems

The Born-Mayer-Huggins (BMH) potential, which combines Coulomb interactions with dispersion and short-range exponential repulsion, is widely used for ionic materials such as molten salts. However, large-scale molecular dynamics simulations of BMH systems are often limited by computation, communication, and memory costs. We recently proposed the random batch sum-of-Gaussians (RBSOG) method, which accelerates Coulomb calculations by using a sum-of-Gaussians (SOG) decomposition to split the potential into short- and long-range parts and by applying importance sampling in Fourier space for the long-range part. In this work, we extend the RBSOG to BMH systems and incorporate a random batch list (RBL) scheme to further accelerate the short-range part, yielding a unified framework for efficient simulations with the BMH potential. The combination of the SOG decomposition and the RBL enables an efficient and scalable treatment of both long- and short-range interactions in BMH system, particularly the RBL well handles the medium-range exponential repulsion and dispersion by the random batch neighbor list. Error estimate is provided to show the theoretical convergence of the RBL force. We evaluate the framework on molten NaCl and mixed alkali halide with up to $5\times10^6$ atoms on $2048$ CPU cores. Compared to the Ewald-based particle-particle particle-mesh method and the RBSOG-only method, our method achieves approximately $4\sim10\times$ and $2\times$ speedups while using $1000$ cores, respectively, under the same level of structural and thermodynamic accuracy and with a reduced memory usage. These results demonstrate the attractive performance of our method in accuracy and scalability for MD simulations with long-range interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved random batch Ewald method in molecular dynamics simulations

The random batch Ewald (RBE) is an efficient and accurate method for molecular dynamics (MD) simulations of physical systems at the nano-/micro- scale. The method shows great potential to solve the computational bottleneck of long-range interactions, motivating a necessity to accelerating short-range components of the non-bonded interactions for a further speedup of MD simulations. In this work, we present an improved RBE method for the non-bonding interactions by introducing the random batch idea to constructing neighbor lists for the treatment of both the short-range part of the Ewald splitting and the Lennard-Jones potential. The efficiency of the novel neighbor list algorithm owes to the stochastic minibatch strategy which can significantly reduce the total number of neighbors. We obtan the error estimate and convergence by theoretical analysis and implement the improved RBE method in the LAMMPS package. Benchmark simulations are performed to demonstrate the accuracy and stability of the algorithm. Numerical tests on computer performance by conducting large-scaled MD simulations for systems including up to 0.1 billion water molecules, run on massive cluster with up to 50 thousand CPU cores, demonstrating the attractive features such as the high parallel scalability and memory-saving of the method in comparison to the existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Random batch sum-of-Gaussians method for molecular dynamics simulations of particle systems

We develop an accurate, highly efficient and scalable random batch sum-of-Gaussians (RBSOG) method for molecular dynamics simulations of systems with long-range interactions. The idea of the RBSOG method is based on a sum-of-Gaussians decomposition of the Coulomb kernel, and then a random batch importance sampling on the Fourier space is employed for approximating the summation of the Fourier expansion of the Gaussians with large bandwidths (the long-range components). The importance sampling significantly reduces the computational cost, resulting in a scalable algorithm by avoiding the use of communication-intensive fast Fourier transform. Theoretical analysis is present to demonstrate the unbiasedness of the approximate force, the controllability of variance and the weak convergence of the algorithm. The resulting method has $\mathcal{O}(N)$ complexity with low communication latency. Accurate simulation results on both dynamical and equilibrium properties of benchmark problems are reported to illustrate the attractive performance of the method. Simulations on parallel computing are also performed to show the high parallel efficiency. The RBSOG method can be straightforwardly extended to more general interactions with long ranged kernels, and thus is promising to construct fast algorithms of a series of molecular dynamics methods for various interacting kernels.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A kernel-independent sum-of-Gaussians method by de la Vallée-Poussin sums

Approximation of interacting kernels by sum of Gaussians (SOG) is frequently required in many applications of scientific and engineering computing in order to construct efficient algorithms for kernel summation or convolution problems. In this paper, we propose a kernel-independent SOG method by introducing the de la Vallée-Poussin sum and Chebyshev polynomials. The SOG works for general interacting kernels and the lower bound of Gaussian bandwidths is tunable and thus the Gaussians can be easily summed by fast Gaussian algorithms. The number of Gaussians can be further reduced via the model reduction based on the balanced truncation based on the square root method. Numerical results on the accuracy and model reduction efficiency show attractive performance of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

HSMA: An O(N) electrostatics package implemented in LAMMPS

We implement two recently developed fast Coulomb solvers, HSMA3D [J. Chem. Phys. 149 (8) (2018) 084111] and HSMA2D [J. Chem. Phys. 152 (13) (2020) 134109], into a new user package HSMA for molecular dynamics simulation engine LAMMPS. The HSMA package is designed for efficient and accurate modeling of electrostatic interactions in 3D and 2D periodic systems with dielectric effects at the O(N) cost. The implementation is hybrid MPI and OpenMP parallelized and compatible with existing LAMMPS functionalities. The vectorization technique following AVX512 instructions is adopted for acceleration. To establish the validity of our implementation, we have presented extensive comparisons to the widely used particle-particle particle-mesh (PPPM) algorithm in LAMMPS and other dielectric solvers. With the proper choice of algorithm parameters and parallelization setup, the package enables calculations of electrostatic interactions that outperform the standard PPPM in speed for a wide range of particle numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Linear-Scaling Selected Inversion based on Hierarchical Interpolative Factorization for Self Green's Function for Modified Poisson-Boltzmann Equation in Two Dimensions

This paper studies an efficient numerical method for solving modified Poisson-Boltzmann (MPB) equations with the self Green's function as a state equation to describe electrostatic correlations in ionic systems. Previously, the most expensive point of the MPB solver is the evaluation of Green's function. The evaluation of Green's function requires solving high-dimensional partial differential equations, which is the computational bottleneck for solving MPB equations. Numerically, the MPB solver only requires the evaluation of Green's function as the diagonal part of the inverse of the discrete elliptic differential operator of the Debye-Hückel equation. Therefore, we develop a fast algorithm by a coupling of the selected inversion and hierarchical interpolative factorization. By the interpolative factorization, our new selected inverse algorithm achieves linear scaling to compute the diagonal of the inverse of this discrete operator. The accuracy and efficiency of the proposed algorithm will be demonstrated by extensive numerical results for solving MPB equations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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 scaling but the poor scalability limits the size of simulated systems. Here, we conduct 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 100 million particles using $10$ thousand CPU 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, micro-phase-separated electrolyte and 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 breakthrough solution on scalability of computing the Coulomb interaction. It is particularly useful and cost-effective to simulate ultra-large systems, which was either impossible or very costing to conduct using existing algorithms, thus would benefit the broad community of sciences.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

L1-based reduced over collocation and hyper reduction for steady state and time-dependent nonlinear equations

The task of repeatedly solving parametrized partial differential equations (pPDEs) in, e.g. optimization or interactive applications, makes it imperative to design highly efficient and equally accurate surrogate models. The reduced basis method (RBM) presents as such an option. Enabled by a mathematically rigorous error estimator, RBM constructs a low-dimensional subspace of the parameter-induced high fidelity solution manifold from which an approximate solution is computed. It can improve efficiency by several orders of magnitudes leveraging an offline-online decomposition procedure. However, this decomposition, usually through the empirical interpolation method (EIM) when the PDE is nonlinear or its parameter dependence nonaffine, is either challenging to implement, or severely degrades online efficiency. In this paper, we augment and extend the EIM approach as a direct solver, as opposed to an assistant, for solving nonlinear pPDEs on the reduced level. The resulting method, called Reduced Over-Collocation method (ROC), is stable and capable of avoiding the efficiency degradation inherent to a traditional application of EIM. Two critical ingredients of the scheme are collocation at about twice as many locations as the dimension of the reduced solution space, and an efficient L1-norm-based error indicator for the strategic selection of the parameter values to build the reduced solution space. Together, these two ingredients render the proposed L1-ROC scheme both offline- and online-efficient. A distinctive feature is that the efficiency degradation appearing in alternative RBM approaches that utilize EIM for nonlinear and nonaffine problems is circumvented, both in the offline and online stages. Numerical tests on different families of time-dependent and steady-state nonlinear problems demonstrate the high efficiency and accuracy of L1-ROC and its superior stability performance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Understanding depletion induced like-charge attraction from self-consistent field model

The interaction force between likely charged particles/surfaces is usually repulsive due to the Coulomb interaction. However, the counterintuitive like-charge attraction in electrolytes has been frequently observed in experiments, which has been theoretically debated for a long time. It is widely known that the mean field Poisson-Boltzmann theory cannot explain or predict this anomalous feature since it ignores many-body properties. In this paper, we develop efficient algorithm and perform the force calculation between two interfaces using a set of self-consistent equations which properly takes into account the electrostatic correlation and the dielectric-boundary effects. By solving the equations and calculating the pressure with the Debye-charging process, we show that the self-consistent equations could be used to study the attraction between like-charge surfaces from weak-coupling to mediate-coupling regime, and that the attraction is due to the electrostatics-driven entropic force which is significantly enhanced by the dielectric depletion of mobile ions. A systematic investigation shows that the interaction forces can be tuned by material permittivity, ionic size and valence, and salt concentration, and that the like-charge attraction exists only for specific regime of these parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Multi-scale Monte Carlo Method for Electrolytes

Artifacts arise in the simulations of electrolytes using periodic boundary conditions (PBC). We show the origin of these artifacts are the periodic image charges and the constraint of charge neutrality inside the simulation box, both of which are unphysical from the view point of real systems. To cure these problems, we introduce a multi-scale Monte Carlo method, where ions inside a spherical cavity are simulated explicitly, whilst ions outside are treated implicitly using continuum theory. Using the method of Debye charging, we explicitly derive the effective interactions between ions inside the cavity, arising due to the fluctuations of ions outside. We find that these effective interactions consist of two types: 1) a constant cavity potential due to the asymmetry of the electrolyte, and 2) a reaction potential that depends on the positions of all ions inside. Combining the Grand Canonical Monte Carlo (GCMC) with a recently developed fast algorithm based of image charge method, we perform a multi-scale Monte Carlo simulation of symmetric electrolytes, and compare it with other simulation methods, including PBC+GCMC method, as well as large scale Monte Carlo simulation. We demonstrate that our multi-scale MC method is capable of capturing the correct physics of a large system using a small scale simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electrostatic interaction in the presence of dielectric interfaces and polarization-induced like-charge attraction

Electrostatic polarization is important in many nano-/micro-scale physical systems such as colloidal suspensions, biopolymers, and nanomaterials assembly. The calculation of polarization potential requires an efficient algorithm for solving 3D Poisson's equation. We have developed a useful image charge method to rapid evaluation of the Green's function of the Poisson's equation in the presence of spherical dielectric discontinuities. This paper presents an extensive study of this method by giving an convergence analysis and developing a coarse-graining algorithm. The use of the coarse graining could reduce the number of image charges to around a dozen, by 1-2 orders of magnitude. We use the algorithm to investigate the interaction force between likely charged spheres in different dielectric environments. We find the size and charge asymmetry leads to an attraction between like charges, in agreement with existing results. Furthermore, we study three-body interaction and find in the presence of an external interface, the interaction force depends on the curvature of the interface, and behavior a non-monotonic electrostatic force.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-consistent field model for strong electrostatic correlations and inhomogeneous dielectric media

Electrostatic correlations and variable permittivity of electrolytes are essential for exploring many chemical and physical properties of interfaces in aqueous solutions. We propose a continuum electrostatic model for the treatment of these effects in the framework of the self-consistent field theory. The model incorporates a space-or field-dependent dielectric permittivity and an excluded ion-size effect for the correlation energy. This results in a self-energy modified Poisson-Nernst-Planck or Poisson-Boltzmann equation together with state equations for the self energy and the dielectric function. We show that the ionic size is of significant importance in predicting a finite self energy for an ion in an inhomogeneous medium. Asymptotic approximation is proposed for the solution of a generalized Debye-Huckel equation, which has been shown to capture the ionic correlation and dielectric self energy. Through simulating ionic distribution surrounding a macroion, the modified self-consistent field model is shown to agree with particle-based Monte Carlo simulations. Numerical results for symmetric and asymmetric electrolytes demonstrate that the model is able to predict the charge inversion at high correlation regime in the presence of multivalent interfacial ions which is beyond the mean-field theory, and also show strong effect to double layer structure due to the space- or field-dependent dielectric permittivity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-energy-modified Poisson-Nernst-Planck equations: WKB approximation and finite-difference approaches

We propose a modified Poisson-Nernst-Planck (PNP) model to investigate charge transport in electrolytes of inhomogeneous dielectric environment. The model includes the ionic polarization due to the dielectric inhomogeneity and the ion-ion correlation. This is achieved by the self energy of test ions through solving a generalized Debye-Hückel (DH) equation. We develop numerical methods for the system composed of the PNP and DH equations. Particularly, towards the numerical challenge of solving the high-dimensional DH equation, we developed an analytical WKB approximation and a numerical approach based on the selective inversion of sparse matrices. The model and numerical methods are validated by simulating the charge diffusion in electrolytes between two electrodes, for which effects of dielectrics and correlation are investigated by comparing the results with the prediction by the classical PNP theory. We find that, at the length scale of the interface separation comparable to the Bjerrum length, the results of the modified equations are significantly different from the classical PNP predictions mostly due to the dielectric effect. It is also shown that when the ion self energy is in weak or mediate strength, the WKB approximation presents a high accuracy, compared to precise finite-difference results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulation of Electric Double Layers around Charged Colloids in Aqueous Solution of Variable Permittivity

The ion distribution around charged colloids in solution has been investigated intensely during the last decade. However, few theoretical approaches have included the influence of variation in the dielectric permittivity within the system, let alone in the surrounding solvent. In this article, we introduce two relatively new methods that can solve the Poisson equation for systems with varying permittivity. The harmonic interpolation method (HIM) approximately solves the Green's function in terms of a spherical harmonics series, and thus provides analytical ion-ion potentials for the Hamiltonian of charged systems. The Maxwell Equations Molecular Dynamics (MEMD) algorithm features a local approach to electrostatics, allowing for arbitrary local changes of the dielectric constant. We show that the results of both methods are in very good agreement. We also found that the renormalized charge of the colloid, and with it the effective far field interaction, significantly changes if the dielectric properties within the vicinity of the colloid are changed.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mellin Transform and Image Charge Method for Dielectric Sphere in an Electrolyte

We revisit the image charge method for the Green's function problem of the Poisson-Boltzmann equation for a dielectric sphere immersed in ionic solutions. Using finite Mellin transformation, we represent the reaction potential due to a source charge inside the sphere in terms of one dimensional distribution of image charges. The image charges are generically composed of a point image at the Kelvin point and a line image extending from the Kelvin point to infinity with an oscillatory line charge strength. We further develop an efficient and accurate algorithm for discretization of the line image using Padé approximation and finite fraction expansion. Finally we illustrate the power of our method by applying it in a multiscale reaction-field Monte Carlo simulation of monovalent electrolytes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solving Fluctuation-Enhanced Poisson-Boltzmann Equations

Electrostatic correlations and fluctuations in ionic systems can be described within an extended Poisson-Boltzmann theory using a Gaussian variational form. The resulting equations are challenging to solve because they require the solution of a non-linear partial differential equation for the pair correlation function. This has limited existing studies to simple approximations or to one-dimensional geometries. In this paper we show that the numerical solution of the equations is greatly simplified by the use of selective inversion of a finite difference operator which occurs in the theory. This selective inversion preserves the sparse structure of the problem and leads to substantial savings in computer effort. In one and two dimensions further simplifications are made by using a mixture of selective inversion and Fourier techniques.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

From Over-charging to Like-charge Attraction in the Weak Coupling Regime

Despite decades of intensive studies, the effective interactions between strongly charged colloids still remain elusive. Here we show that a strongly charged surface with a layer of condensed counter- ions behaves effectively as a conductor, due to the mobile nature of the condensed ions. An external source charge in its vicinity is therefore attracted towards the surface, due to the image charge effect. This mechanism leads to correlational energies for counter-ions condensed on two distinct surfaces, as well as for free ions in the bulk. Generalizing Debye-Huckel theory and image charge methods, we analytically calculate these correlation energies for the two-plates problem, at the iso-electric point, where condensed counterions precisely balance the bare surface charges. At this point, the effective interaction between two plates is always attractive at small separation and repulsive at large separation.

0 citations0 reviewsNo code linkedOpen access