Related papers: MMM2D: A fast and accurate summation method for el…
Ewald summation and physically equivalent methods such as particle-mesh Ewald, kubic-harmonic expansions, or Lekner sums are commonly used to calculate long-range electrostatic interactions in computer simulations of polar and charged…
We derive analytic solutions for the potential and field in a one-dimensional system of masses or charges with periodic boundary conditions, in other words Ewald sums for one dimension. We also provide a set of tools for exploring the…
The constant potential molecular dynamics simulation method proposed by Siepmann and Sprik and reformulated later by Reed (SR-CPM) has been widely employed to investigate the metallic electrolyte/electrode interfaces, especially for…
Quantum simulations of photoexcited low-dimensional systems are pivotal for understanding how to functionalize and integrate novel two-dimensional (2D) materials in next-generation optoelectronic devices. First principles predictions are…
A direct summation method for the Madelung constant calculation is presented where a crystal lattice is constructed from linear arrays of charges or axial multipoles. An array is designed to have vanishing low order electric moments such…
The approximate computation of all gravitational forces between $N$ interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than $\mathcal{O}(N)$ operations. FMM groups…
We carry out the Ewald summation for the Rotne-Prager-Yamakawa mobility tensor, the Oseen mobility tensor and further variations of both, relevant for the hydrodynamic interactions in colloidal suspensions, where all interacting particles…
A new fully discrete linearized $H^1$-conforming Lagrange finite element method is proposed for solving the two-dimensional magneto-hydrodynamics equations based on a magnetic potential formulation. The proposed method yields numerical…
QM/MM hybrid methods employ accurate quantum (QM) models only in regions of interest (defects) and switch to computationally cheaper interatomic potential (MM) models to describe the crystalline bulk. We develop two QM/MM hybrid methods for…
In this paper, we mainly focus on the rigorous convergence analysis of two fully decoupled, unconditionally energy-stable methods for the diffuse interface two-phase magnetohydrodynamics (MHD) model. The two methods consist of the…
The very-low temperature thermal effective mass m* of paramagnetic and ferromagnetic electrons in a uniform electron fluid in two dimensions is studied. Analytical and numerical evaluations are used to meaningfully define an m*, even in the…
We present a new ab initio method for calculating effective onsite Coulomb interactions of itinerant and strongly correlated electron systems. The method is based on constrained local density functional theory formulated in terms of…
Many-body interactions arise naturally in the perturbative treatment of classical and quantum many-body systems and play a crucial role in the description of condensed matter systems. In the case of three-body interactions, the…
The lack of long-range electrostatics is a key limitation of modern machine learning interatomic potentials (MLIPs), hindering reliable applications to interfaces, charge-transfer reactions, polar and ionic materials, and biomolecules. In…
Classical Ewald methods for Coulomb and Stokes interactions rely on ``kernel-splitting," using decompositions based on Gaussians to divide the resulting potential into a near field and a far field component. Here, we show that a more…
A self-consistent theory of bulk electrolytes incorporating electrostatic and hard-core interactions on an equal level is applied to the two-dimensional Coulomb liquid with finite ion size. The ionic pair distributions, the structure…
We present new efficient (O(N log N)) methods for computing three quantities crucial to electronic structure calculations: the ionic potential, the electron-ion contribution to the Born-Oppenheimer forces, and the electron-ion contribution…
The electron-electron pair distribution functions (PDF) of the 2-D electron fluid (2DEF) in the quantum regime (at T=0) are calculated using a classical-map-hyper-netted-chain (CHNC) scheme and compared with currently available Quantum…
Ewald summation is an efficient method for computing the periodic sums that appear when considering the Green's functions of Stokes flow together with periodic boundary conditions. We show how Ewald summation, and accompanying truncation…
Multiscale finite element methods for 2D/1D problems have been studied in this work to demonstrate their excellent ability to solve real-world problems. These methods are much more efficient than conventional 3D finite element methods and…