Related papers: Random Batch Sum-of-Gaussians Method for Molecular…
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…
In this work, we develop a random batch sum-of-Gaussians (RBSOG) method for molecular dynamics simulations of charged systems in the isothermal-isobaric (NPT) ensemble. We introduce an SOG splitting of the pressure-related $1/r^3$ kernel,…
Yukawa systems have drawn widespread interest across various applications. In this paper, we introduce a novel random batch sum-of-Gaussians (RBSOG) algorithm for molecular dynamics simulations of 3D Yukawa systems with periodic boundary…
This paper introduces a random-batch molecular dynamics (RBMD) package for fast simulations of particle systems at the nano/micro scale. Different from existing packages, the RBMD uses random batch methods for nonbonded interactions of…
We develop a random batch Ewald (RBE) method for molecular dynamics simulations of particle systems with long-range Coulomb interactions, which achieves an $O(N)$ complexity in each step of simulating the $N$-body systems. The RBE method is…
The embedded atom method (EAM) is one of the most widely used many-body, short-range potentials in molecular dynamics simulations, particularly for metallic systems. To enhance the efficiency of calculating these short-range interactions,…
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…
We have extended the multilevel summation (MLS) method, originally developed to evaluate long-range Coulombic interactions in molecular dynamics (MD) simulations [Skeel et al., J. Comput. Chem., 23, 673 (2002)], to handle dispersion…
Machine-learning interatomic potentials have emerged as a revolutionary class of force-field models in molecular simulations, delivering quantum-mechanical accuracy at a fraction of the computational cost and enabling the simulation of…
The stochastic particle method based on Bhatnagar-Gross-Krook (BGK) or ellipsoidal statistical BGK (ESBGK) model approximates the pairwise collisions in the Boltzmann equation using a relaxation process. Therefore, it is more efficient to…
We present a hybrid method for the simulation of colloidal systems, that combines molecular dynamics (MD) with the Lattice-Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the…
Simulating hypersonic flow around a space vehicle is challenging because of the multiscale and nonequilibrium nature inherent in these flows. To effectively deal with such flows, a novel particle particle hybrid scheme combining the…
Immersed boundary-lattice Boltzmann method (IB-LBM) has been widely used for simulation of particle-laden flows recently. However, it was limited to small-scale simulations with no more than O(103) particles. Here, we expand IB-LBM for…
Extended Lagrangian Born-Oppenheimer molecular dynamics (XL-BOMD) [Phys. Rev. Lett. vol. 100, 123004 (2008)] is combined with Kohn-Sham density functional theory (DFT) using a DFT+U correction based on the Hubbard model. This combined…
We develop Random Batch Methods for interacting particle systems with large number of particles. These methods use small but random batches for particle interactions, thus the computational cost is reduced from $O(N^2)$ per time step to…
Modeling and direct numerical simulation of particle-laden flows have a tremendous variety of applications in science and engineering across a vast spectrum of scales from pollution dispersion in the atmosphere, to fluidization in the…
In many real-world scenarios, the underlying random fluctuations are non-Gaussian, particularly in contexts where heavy-tailed data distributions arise. A typical example of such non-Gaussian phenomena calls for L\'evy noise, which…
We propose an accurate, efficient, and low-memory sum-of-Gaussians tensor neural network (SOG-TNN) algorithm for solving the high-dimensional Schr\"{o}dinger equation. The SOG-TNN utilizes a low-rank tensor product representation of the…
In this work we present an efficient implementation of Canonical Monte Carlo simulation for Coulomb many body systems on graphics processing units (GPU). Our method takes advantage of the GPU Single Instruction, Multiple Data (SIMD)…
Precise physical descriptions of molecules can be obtained by solving the Schrodinger equation; however, these calculations are intractable and even approximations can be cumbersome. Force fields, which estimate interatomic potentials based…