Machine learning interatomic potentials(MLIPs)often neglect long-range interactions,such as electrostatic and dispersion forces.In this work,we introduce a straightforward and efficient method to account for long-rang...Machine learning interatomic potentials(MLIPs)often neglect long-range interactions,such as electrostatic and dispersion forces.In this work,we introduce a straightforward and efficient method to account for long-range interactions by learning a hidden variable from local atomic descriptors and applying an Ewald summation to this variable.We demonstrate that in systems including charged and polar molecular dimers,bulk water,and water-vapor interface,standard short-ranged MLIPs can lead to unphysical predictions even when employing message passing.The long-range models effectively eliminate these artifacts,with only about twice the computational cost of short-range MLIPs.展开更多
We extend the recent formulation of the Ewald sum for electrostatics in a two-dimensionally periodic three-dimensional multi- atom layer or two-dimensional single-atom layer system with a rectangular periodic boundary...We extend the recent formulation of the Ewald sum for electrostatics in a two-dimensionally periodic three-dimensional multi- atom layer or two-dimensional single-atom layer system with a rectangular periodic boundary condition （J Chem Theory, Comput, 2014, 10： 534-542） to that with a parallelogrammic periodic boundary condition in general. Following the discussion of an efficient implementation of the formula, we suggest a simple setup of parameters using a relatively smaller screening factor and the associated larger real space cutoff distance to reach an optimized algorithm of an order N computational cost. The connection between the previous application of the Ewald sum to ionic crystal systems and the future application to mo- lecular self-assembly or disassembly systems on solid surfaces or at liquid-liquid interfaces ate illustrated to demonstrate the applicability of the present work to simulate the self-assembly process and to produce dynamical, structural and thermody- namic properties of experimental self-assembly systems of interest.展开更多
利用简单多项式基组拟合变化平缓的长程作用部分,并增加多个限定条件以提高优化势与库仑势的局部重合度,对Natoli优化法进行改进.在不增加计算复杂度的前提下,得到平均偏差明显降低的优化势;建议优化时划分r空间与k空间的准则取k c r c...利用简单多项式基组拟合变化平缓的长程作用部分,并增加多个限定条件以提高优化势与库仑势的局部重合度,对Natoli优化法进行改进.在不增加计算复杂度的前提下,得到平均偏差明显降低的优化势;建议优化时划分r空间与k空间的准则取k c r c≥15.在稠密液氢分子模拟中,优化法获得的能量及压力均小于Natoli法与Ewald法的结果.展开更多
Ewald summation method, based on Non-Uniform FFTs (ENUF) to compute the electrostatic interactions and forces, is implemented in two different particle simulation schemes to model molecular and soft matter, in classic...Ewald summation method, based on Non-Uniform FFTs (ENUF) to compute the electrostatic interactions and forces, is implemented in two different particle simulation schemes to model molecular and soft matter, in classical all-atom Molecular Dynamics and in Dissipative Particle Dynamics for coarse-grained particles. The method combines the traditional Ewald method with a non-uniform fast Fourier transform library (NFFT), making it highly efficient. It scales linearly with the number of particles as , while being both robust and accurate. It conserves both energy and the momentum to float point accuracy. As demonstrated here, it is straight- forward to implement the method in existing computer simulation codes to treat the electrostatic interactions either between point-charges or charge distributions. It should be an attractive alternative to mesh-based Ewald methods.展开更多
文摘Machine learning interatomic potentials(MLIPs)often neglect long-range interactions,such as electrostatic and dispersion forces.In this work,we introduce a straightforward and efficient method to account for long-range interactions by learning a hidden variable from local atomic descriptors and applying an Ewald summation to this variable.We demonstrate that in systems including charged and polar molecular dimers,bulk water,and water-vapor interface,standard short-ranged MLIPs can lead to unphysical predictions even when employing message passing.The long-range models effectively eliminate these artifacts,with only about twice the computational cost of short-range MLIPs.
基金supported by the National Natural Science Foundation of China(91127015,21103063(Z.H.))
文摘We extend the recent formulation of the Ewald sum for electrostatics in a two-dimensionally periodic three-dimensional multi- atom layer or two-dimensional single-atom layer system with a rectangular periodic boundary condition （J Chem Theory, Comput, 2014, 10： 534-542） to that with a parallelogrammic periodic boundary condition in general. Following the discussion of an efficient implementation of the formula, we suggest a simple setup of parameters using a relatively smaller screening factor and the associated larger real space cutoff distance to reach an optimized algorithm of an order N computational cost. The connection between the previous application of the Ewald sum to ionic crystal systems and the future application to mo- lecular self-assembly or disassembly systems on solid surfaces or at liquid-liquid interfaces ate illustrated to demonstrate the applicability of the present work to simulate the self-assembly process and to produce dynamical, structural and thermody- namic properties of experimental self-assembly systems of interest.
文摘利用简单多项式基组拟合变化平缓的长程作用部分,并增加多个限定条件以提高优化势与库仑势的局部重合度,对Natoli优化法进行改进.在不增加计算复杂度的前提下,得到平均偏差明显降低的优化势;建议优化时划分r空间与k空间的准则取k c r c≥15.在稠密液氢分子模拟中,优化法获得的能量及压力均小于Natoli法与Ewald法的结果.
文摘Ewald summation method, based on Non-Uniform FFTs (ENUF) to compute the electrostatic interactions and forces, is implemented in two different particle simulation schemes to model molecular and soft matter, in classical all-atom Molecular Dynamics and in Dissipative Particle Dynamics for coarse-grained particles. The method combines the traditional Ewald method with a non-uniform fast Fourier transform library (NFFT), making it highly efficient. It scales linearly with the number of particles as , while being both robust and accurate. It conserves both energy and the momentum to float point accuracy. As demonstrated here, it is straight- forward to implement the method in existing computer simulation codes to treat the electrostatic interactions either between point-charges or charge distributions. It should be an attractive alternative to mesh-based Ewald methods.
