Abstract:This paper provides an error estimate for the u-series decomposition of the Coulomb interaction in molecular dynamics simulations. We show that the number of truncated Gaussians $M$ in the u-series and the base of interpolation nodes $b$ in the bilateral serial approximation are two key parameters for the algorithm accuracy, and that the errors converge as $\mathcal{O}(b^{-M})$ for the energy and $\mathcal{O}(b^{-3M})$ for the force. Error bounds due to numerical quadrature and cutoff in both the electrostatic energy and forces are obtained. Closed-form formulae are also provided, which are useful in the parameter setup for simulations under a given accuracy. The results are verified by analyzing the errors of two practical systems.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the error estimation problem when using the u - series method for Coulomb interaction decomposition in molecular dynamics simulations. Specifically, the paper focuses on the truncation error of the u - series method when dealing with near - field interactions and how these errors affect the computational accuracy of the electrostatic potential energy and force. Through analysis, the authors hope to provide a rigorous error estimate to guide the parameter selection in actual simulations and ensure that molecular dynamics simulations can be carried out efficiently under given precision requirements. ### Background of the Paper Molecular dynamics (MD) simulations are widely used numerical tools in multi - body particle systems and are applied in fields such as biophysics, chemical physics, and soft materials. In MD simulations, the interactions between particles usually follow Newton's equations of motion, and the dynamic and thermodynamic properties of the system are captured through the ensemble average of particle configurations. However, due to the existence of long - range Coulomb interactions, MD simulations face huge computational bottlenecks, especially the low scalability in large - scale parallel computing. ### The u - series Method The u - series method is a recently developed algorithm for electrostatic calculations. It represents the Coulomb potential through the bilateral series approximation of Gaussian functions (sum - of - Gaussians, SOG). Compared with the traditional Ewald decomposition method, the u - series method has the following advantages: 1. **Smoothness**: The SOG decomposition can maintain the smoothness of the truncation point at any order, improving the accuracy under the same computational cost. 2. **Communication Efficiency**: The separability of SOG helps to reduce the communication rounds in FFT usage. 3. **Balance between Computation and Communication**: The u - series method achieves a good balance between communication latency, communication bandwidth, and computational cost, and is suitable for high - precision parallel computing. ### Research Objectives Although the u - series method performs well and is accurate in practical applications, previous studies lack a rigorous estimate of the convergence rate caused by near - field truncation. This leads to difficulties in selecting optimal parameters in systems of different scales. Therefore, this paper aims to fill this gap and provide a rigorous error estimate for electrostatic energy and force. ### Main Results By estimating the near - field error of each Fourier mode, the authors prove that the decay rates of the total energy error and the force error are $ O(b^{-M}) $ and $ O(b^{-3M}) $ respectively, where $ b $ is the basis of the interpolation nodes and $ M $ is the number of Gaussian functions. In addition, closed - form formulas for the approximation error are also provided, which are very useful for parameter optimization in actual MD simulations. ### Conclusions The error estimates provided in this paper not only support the theoretical basis of the u - series method but also provide clear guidance for parameter selection in practical applications. This is of great significance for improving the efficiency and accuracy of molecular dynamics simulations.

Error estimate of the u-series method for molecular dynamics simulations

The HI Content and Extent of Low Surface Brightness Galaxies - Could LSB Galaxies be Responsible for Damped Ly-alpha Absorption?

A posteriori error estimates for Schr{ö}dinger operators discretized with linear combinations of atomic orbitals

Symmetries of the Kac-Peterson Modular Matrices of Affine Algebras

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Approaching the complete-basis limit with a truncated many-body expansion

Error Interference in Quantum Simulation

Realizing Effective Cubic-Scaling Coulomb Hole Plus Screened Exchange Approximation in Periodic Systems Via Interpolative Separable Density Fitting with a Plane-Wave Basis Set.

Multi-Reference UCCSD Variational Quantum Algorithm for Molecular Ground State Energies

Error Estimate of a Quasi-Monte Carlo Time-Splitting Pseudospectral Method for Nonlinear Schrödinger Equation with Random Potentials

The need for accuracy and smoothness in numerical simulations

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

A note on the accuracy of adaptive Gauss–Hermite quadrature

Practical error bounds for properties in plane-wave electronic structure calculations

Range-separated double-hybrid density-functional theory applied to periodic systems

Ensembles Are Required to Handle Aleatoric and Parametric Uncertainty in Molecular Dynamics Simulation

Coupling the molecular motion and collision processes in numerical simulations

Fully guaranteed and computable error bounds on the energy for periodic Kohn-Sham equations with convex density functionals

Covariant Bethe-Salpeter Approximation in Models of Strongly Correlated Electron Systems

Tight distance-dependent estimators for screening two-center and three-center short-range Coulomb integrals over Gaussian basis functions

A posteriori error estimates for fully coupled McKean-Vlasov forward-backward SDEs