Related papers: MMM2D: A fast and accurate summation method for el…
The calculation of the MP2 correlation energy for extended systems can be viewed as a multi-dimensional integral in the thermodynamic limit, and the standard method for evaluating the MP2 energy can be viewed as a trapezoidal quadrature…
Electric double layers (EDL) with counterions only, say electrons with the elementary charge $-e$, in thermal equilibrium at the inverse temperature $\beta$ are considered. In particular, we study the effect of the surface charge modulation…
The 2D/1D multiscale finite element method (MSFEM) is an efficient way to simulate rotating machines in which each iron sheet is exposed to the same field. It allows the reduction of the three dimensional sheet to a two dimensional…
Based on the plenty method, this paper describes a numerical method for 2D non-smooth contact problems with Coulomb friction and bilateral constraints and its application to the simulation of statics and dynamics for a frictional…
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…
The field of two-dimensional (2D) materials has grown dramatically in the last two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computation applications. These…
In the first two papers in this series, we developed new shifted potential (SP), gradient shifted force (GSF), and Taylor shifted force (TSF) real-space methods for multipole interactions in condensed phase simulations. Here, we discuss the…
We consider two related problems: the first is the minimization of the "Coulomb renormalized energy" of Sandier-Serfaty, which corresponds to the total Coulomb interaction of point charges in a uniform neutralizing background (or rather…
We describe a set of techniques for performing large scale ab initio calculations using multigrid accelerations and a real-space grid as a basis. The multigrid methods provide effective convergence acceleration and preconditioning on all…
A variational approach is used to calculate free energy and conformational properties in polyelectrolytes. The true bond and Coulomb potentials are approximated by a trial isotropic harmonic energy containing monomer-monomer force constants…
We show that finite-range alternatives to the standard long-range BKS pair potential for silica might be used in molecular dynamics simulations. We study two such models that can be efficiently simulated since no Ewald summation is…
We review our recently developed electronic structure calculation methods used for the dynamics of large-scale solids or liquids with an efficient algorithm for large scale simultaneous linear equations. The electronic structure calculation…
Most current machine learning interatomic potentials (MLIPs) rely on short-range approximations, without explicit treatment of long-range electrostatics. To address this, we recently developed the Latent Ewald Summation (LES) method, which…
We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin.…
We implement a total-energy minimization scheme to allow for relaxation of atomic positions in density functional calculations for two-dimensional (2D) systems using a mixed basis set. The basis functions consist of products of 2D plane…
The formation energy of a solid surface can be extracted from slab calculations if the bulk energy per atom is known. It has been pointed out previously that the resulting surface energy will diverge with slab thickness if the bulk energy…
Recently, Gimbutas et al derived an elegant representation for the Green's functions of Stokes flow in a half-space. We present a fast summation method for sums involving these half-space Green's functions (stokeslets, stresslets and…
We develop a computational method for modeling electrostatic interactions of arbitrarily-shaped, polarizable objects on colloidal length scales, including colloids/nanoparticles, polymers, and surfactants, dispersed in explicit ion…
Null collapse is an implosive process whereby MHD waves focus their energy in the vicinity of a null point, forming a current sheet and initiating magnetic reconnection. We consider, for the first time, the case of collapsing 3D magnetic…
Fast Ewald summation efficiently evaluates Coulomb interactions and is widely used in molecular dynamics simulations. It is based on a split into a short-range and a long-range part, where evaluation of the latter is accelerated using the…