本研究利用耗散粒子動力學(Dissipative Particle Dynamics, DPD)模擬，為了擴大空間與時間的尺度，進而研究在介觀尺度下水分子的靜電感應行為 (Electrostatic Induction Behavior)。過去，關於靜電感應的勢能，在分子動力學上常利用Ewald Summation，因為可以避開在分子間的距離大於截斷半徑時，其靜電勢能會發散的問題。Ewald Summation應用在介觀系統耗散粒子動力學模擬時，在分子間距離為零時，其勢能會發生發散。直到近幾年，這個問題已經被解決。因此，本研究目的分為五部分：
一、 Ewald summation系統的建立以及系統驗證。
二、 從微觀到介觀，Ewald summation修正項的驗證。
三、 建立Ewald summation在耗散粒子動力學模擬系統上。
四、 擴散系統的建立以及系統驗證。
五、 利用耗散粒子動力學，模擬介觀尺度下的水分子，了解其靜電感應行為，並且預測擴散係數。
從最後結果，可以知道隨著部分電荷增加，擴散係數會降低。固定電荷下，隨著Nm加大，擴散係數也會遞減，原因是因為質量影響大過靜電力的影響。

This research, use dissipative particle dynamics simulations,in order to expand the space and time scales,and study the water molecules in the mesoscale electrostatic induction behavior.In past,the electrostatic induction potential regular used Ewald summation,in order to avoid the molecular distance is greater than cut-off radius,the electrostatic potential energy will be divergent problems.Ewald Summation application to mesoscopic systems dissipative particle dynamics simulation,the distance between the molecules is zero,the potential divergence will occur.Until recent years,this problem has been resolved.Therefore,the purpose of this research is divided into five parts: I.Set up the Ewald summation system and verify system. II.From micro to meso,verify Ewald summation correction equation. III.Set up the Ewald summation in the dissipative particle dynamics simulation system. IV.Set up the Diffusion system and verify system. V.Use the dissipative particle dynamics simulation,understand the water molecules in the mesoscale electrostatic induction behavior and prediction the diffusion coefficient. From the result,Know increase the partial charge,the diffusion coefficient will decrease.Fixed charge,the Nm increase,the diffusion coefficient will decrease,because the mass impact is large than static electricity. Key words: Dissipative particle dynamics 、Ewald summation、Electrostatic induction、Water molecules

[1] P. J. Hoogerbrugge and J. M. V. A. Koelma, "Simulating Microscopic Hydrodynamic Phenomena with Dissipative Particle Dynamics.," Europhys. Lett., 155 (1992).
[2] P. EspanoL and P. Warren, "Statistical Mechanics of Dissipative Particle Dynamics.," Europhys. Lett. 30, 191 (1995).
[3] R. D. Groot and S. D. Stoyanov, "Mesoscopic Model for Colloidal Particles, Powders, and Granular Solids," Phys. Rev. E 78 (2008).
[4] R. D. Groot, "Electrostatic Interactions in Dissipative Particle Dynamics—Simulation of Polyelectrolytes and Anionic Surfactants," J. Chem. Phys. 118, 11265 (2003).
[5] R. D. Groot, "Mesoscopic Simulation of Polymer-Surfactant Aggregation," Langmuir 16, 7493 (2000).
[6] R. D. Groot, T. J. Madden, and D. J. Tildesley, "On the Role of Hydrodynamic Interactions in Block Copolymer Microphase Separation," J. Chem. Phys. 110, 9739 (1999).
[7] R. D. Groot and T. J. Madden, "Dynamic Simulation of Diblock Copolymer Microphase Separation," J. Chem. Phys. 108, 8713 (1998).
[8] R. D. Groot and P. B. Warren, "Dissipative Particle Dynamics: Bridging the Gap between Atomistic and Mesoscopic Simulation," J. Chem. Phys. 107, 4423 (1997).
[9] Z. Wu, Q. Cui, and A. Yethiraj, "A New Coarse-Grained Model for Water: The Importance of Electrostatic Interactions," J. Chem. Phys. 114, 10524 (2010).
[10] 邱佑宗, "耗散粒子動力學計算方法簡介及應用," 工業材料 213 (2004).
[11] G. M. Minerva, E. Mayoral, M. E. Velázquez, and J. Alejandreb, "Electrostatic Interactions in Dissipative Particle Dynamics Using the Ewald Sums," J. Chem. Phys. 125, 224107 (2006).
[12] C. Ibergay, P. Malfreyt, and D. J. Tildesley, "Electrostatic Interactions in Dissipative Particle Dynamics: Toward a Mesoscale Modeling of the Polyelectrolyte Brushes," J. Chem. Theory Comput. 5, 3245 (2009).
[13] S. Nose, "A Unified Formulation of the Constant Temperature Molecular Dynamics Methods," J. Chem. Phys. 81, 511 (1984).
[14] M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (1991)
[15] M. Kawata and U. Nagashima, "Particle Mesh Ewald Method for Three-Dimensional System with Two-Dimensional Periodicity," Chem. Phys. Lett. 340, 165 (2001).
[16] A. Bro´dka, "Optical Parameter Values of the Ewald Method with Electrostatic Layer Correction for Coulomb Interactions in Slab Geometry," J. Mol. Struct. 792-793, 56 (2006).
[17] I. C. Yeh and M. L. Berkowitza, "Ewald Summation for Systems with Slab Geometry," J. Chem. Phys. 111, 3155 (1999).
[18] C. L. Brooks, B. M. Pettitt, and M. Karplus, "Structure and Energetic Effects of Truncating Long Ranged Interactions in Ionic and Polar Fluids," J. Chem. Phys. 83, 5897 (1985).
[19] J. H. Haile, Molecular Dynamics Simulation (1991)