Lu Guiwu, Li Chunxi, Wang Wenchuan, Wang Zihao. Electrostatic Energy Calculation and Parameter Optimization in Computer Molecular Simulation[J]. Chinese Journal of Chemical Physics , 2004, 17(5): 547-553. doi: 10.1088/1674-0068/17/5/547-553
Citation: Lu Guiwu, Li Chunxi, Wang Wenchuan, Wang Zihao. Electrostatic Energy Calculation and Parameter Optimization in Computer Molecular Simulation[J]. Chinese Journal of Chemical Physics , 2004, 17(5): 547-553. doi: 10.1088/1674-0068/17/5/547-553

Electrostatic Energy Calculation and Parameter Optimization in Computer Molecular Simulation

doi: 10.1088/1674-0068/17/5/547-553
Funds:  Project supported by the National Natural Science Foundation of China (20076001, 20236010) and the Ph.D. Programs Foundation of Ministry of Education of China (20020302039).
  • Received Date: 2003-08-25
  • Molecular simulation is a powerful tool in studying properties of complex fluids composed of charged particles such as electrolyte solutions, room temperature ionic liquids and colloid solutions, where the long-range interactions play a determinative role. Several methods have been available for treating the long-range interactions between charged particles and point dipole. These include the Ewald sum (ES), the reaction field and particle-particle particle-mesh methods. Among these approaches, the ES is most commonly used. However, several ES versions have been occurred in the literatures of molecular simulation, and some times it is difficult for one to choose the right formula to use in molecular simulation study. The coulombic interaction energy between charged particles is divided into the sum of real space, reciprocal space and self energy, and theircal culation equations are obtained respectively using electrostatics theory and Fourier transformation method. The Ewald sum formulate have been derived and the clear physical picture involved has been depicted. A Monte Carlo computer simulation for electrostatic interaction energy of charged hard sphere system has been conducted at varying conditions, and a good agreement with MSA is obtained. On this basis the effects on the simulation accuracy and efficiency of real space cut distance r(cut), convergence parameter 1 and reciprocal maximum vector K(max) have been analyzed. From the theoretical and computer simulation presented here, the optimization parameters for rcut, 1 and K(max) are obtained as 0.5 L (i.e. half box length), 5.8 and 3~5, respectively.
  • 加载中
  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Article Metrics

Article views(1487) PDF downloads(2003) Cited by()

Proportional views
Related

Electrostatic Energy Calculation and Parameter Optimization in Computer Molecular Simulation

doi: 10.1088/1674-0068/17/5/547-553
Funds:  Project supported by the National Natural Science Foundation of China (20076001, 20236010) and the Ph.D. Programs Foundation of Ministry of Education of China (20020302039).

Abstract: Molecular simulation is a powerful tool in studying properties of complex fluids composed of charged particles such as electrolyte solutions, room temperature ionic liquids and colloid solutions, where the long-range interactions play a determinative role. Several methods have been available for treating the long-range interactions between charged particles and point dipole. These include the Ewald sum (ES), the reaction field and particle-particle particle-mesh methods. Among these approaches, the ES is most commonly used. However, several ES versions have been occurred in the literatures of molecular simulation, and some times it is difficult for one to choose the right formula to use in molecular simulation study. The coulombic interaction energy between charged particles is divided into the sum of real space, reciprocal space and self energy, and theircal culation equations are obtained respectively using electrostatics theory and Fourier transformation method. The Ewald sum formulate have been derived and the clear physical picture involved has been depicted. A Monte Carlo computer simulation for electrostatic interaction energy of charged hard sphere system has been conducted at varying conditions, and a good agreement with MSA is obtained. On this basis the effects on the simulation accuracy and efficiency of real space cut distance r(cut), convergence parameter 1 and reciprocal maximum vector K(max) have been analyzed. From the theoretical and computer simulation presented here, the optimization parameters for rcut, 1 and K(max) are obtained as 0.5 L (i.e. half box length), 5.8 and 3~5, respectively.

Lu Guiwu, Li Chunxi, Wang Wenchuan, Wang Zihao. Electrostatic Energy Calculation and Parameter Optimization in Computer Molecular Simulation[J]. Chinese Journal of Chemical Physics , 2004, 17(5): 547-553. doi: 10.1088/1674-0068/17/5/547-553
Citation: Lu Guiwu, Li Chunxi, Wang Wenchuan, Wang Zihao. Electrostatic Energy Calculation and Parameter Optimization in Computer Molecular Simulation[J]. Chinese Journal of Chemical Physics , 2004, 17(5): 547-553. doi: 10.1088/1674-0068/17/5/547-553

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return