請用此 Handle URI 來引用此文件：
http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/29233完整後設資料紀錄
| DC 欄位 | 值 | 語言 |
|---|---|---|
| dc.contributor.advisor | 趙聖德 | |
| dc.contributor.author | "Pai, Ching-Cheng" | en |
| dc.contributor.author | 白景政 | zh_TW |
| dc.date.accessioned | 2021-06-13T01:03:08Z | - |
| dc.date.available | 2010-07-30 | |
| dc.date.copyright | 2007-07-30 | |
| dc.date.issued | 2007 | |
| dc.date.submitted | 2007-07-23 | |
| dc.identifier.citation | [1] O. S. Tyagi, H. S. Bisht, and A. K. Chatterjee, J. Phys. Chem. B 108, 3010 (2004).
[2] P. Hobza and R. Zahradnik, Intermolecular Complexes : the role of van der Waals systems in physical chemistry and in the biodisciplines (Elsevier, New York, 1988). [3] G. A. Jeffrey and W. Saenger, Hydrogen Bonding in Biological Structures (Springer-Verlag, Berlin, 1991). [4] A. J. Stone, The Theory of Intermolecular Forces (Oxford University Press, Oxford, 1996). [5] H. Margenau, Rev. Mod. Phys. 11, 1 (1939). [6] A. van der Avoird, P. E. S. Wormer, and R. Moszynski, Chem. Rev. 94, 1931 (1994). [7] A. K. Rappe and E. R. Bernstein, J. Phys. Chem. A 104, 6117 (2000). [8] G. Chalasinski and M. M. Szczesniak, Chem. Rev. 100, 4227 (2000). [9] R. J. Wheatley, A. S. Tulegenov, and E. Bichoutskaia, Int. Rev. Phys. Chem. 23, 151 (2004). [10] Y. Zhao and D. G. Truhlar, J. Chem. Theory Comput. 1, 415 (2005). [11] S. Grimme, J. Comp. Chem. 25, 1463 (2004). [12] R. A. Friesner, Proc. Natl. Acad. Sci. 102, 6648 (2005). [13] D. Frenkel and B. Smit, Understanding Molecular Simulations (Academic Press, New York, 2002). [14] M. Marques and E. Gross, Annual Rev. Phys. Chem. 55, 427 (2004) [15] C. E. Dykstra, G. Frenking, K. S. Kim, and G. E. Scuseria, Theory and Applications of Computational Chemistry : The first forty years (Elsevier, Amsterdam, 2005). [16] A. Ruzsinszky, J. P. Perdew, and G. I. Csonka, J. Phys. Chem. A 109, 11015 (2005). [17] M. M. Szczesniak, G. Chalasinski, S. M. Cybulski, and S. Scheiner, J. Chem. Phys. 93, 4243 (1990). [18] R. L. Rowley and T. Pakkanen, J. Chem. Phys. 110, 3368 (1999). [19] D. E. Williams and D. J. Craycroft, J. Phys. Chem. 91, 6365 (1987). [20] S. Tsuzuki and K. Tanabe, J. Phys. Chem. 95, 2272 (1991) [21] J. J. Novoa, M. –H. Whangho, and J. M. Williams, J. Chem. Phys. 94, 4835 (1991) [22] D. H. Gay, H. Dai, and D. R. Beck, J. Chem. Phys. 95, 9106 (1991). [23] S. Tsuzuki, T. Uchimaru, K. Tanabe, and S. Kuwajima, J. Phys. Chem. 98, 1830 (1994). [24] S. Tsuzuki and H. P. Luthi, J. Chem. Phys. 114, 3949 (2001). [25] J. Nagy, D. F. Weaver, and V. H. Smith, Mol. Phys. 85, 1179 (1995). [26] E. Fraschini and A. J. Stone, J. Comput. Chem. 19, 847 (1998). [27] S. Tsuzuki, T. Uchimaru, and K. Tanabe, Chem. Phys. Lett. 287, 202 (1998). [28] S. Tsuzuki, K. Honda and T. Uchimaru, M. Mikami, J. Chem. Phys. 124, 114304 (2006). [29] B. J. Alder and T. E. Wainwright, Phys. Rev. 127, 359 (1962). [30] C. H. Mak,Phys. Rev. E, 73, 065104(R) (2006). [31] K. Takahashi, Y. Hayashi, and K. Tachibana, Jpn. J. Appl.Phys., Part 1, 38, 4561 (1999). [32] S.C. Davey, J. Budai, J.W. Goodby, R. Pindak, and D.E. Moncton, Phys. Rev. Lett. 53, 2129(1984). [33] F. Gallet, G. Deville, A. Vald`es, and F.I.B. Williams, Phys. Rev. Lett. 49, 212 (1982). [34] C.J. Guo, D.B. Mast, R. Mehrotra, Y.-Z. Ruan, M.A. Stan, and A.J. Dahm, Phys. Rev. Lett.51, 1461 (1983). [35] J.P. McTague, J. Als-Nielsen, J. Bohr, and M. Nielsen, Phys. Rev. B 25, 7765 (1982). [36] M. Nielsen, J. Als-Nielsen, J. Bohr, J.P. McTague, D.E. Moncton, and P.W. Stephens, Phys.Rev. B 35, 1419 (1987). [37] N. Greiser, G.A. Held, R. Frahm, R.L. Greene, P.M. Horn, and R.M. Suter , Phys. Rev. Lett.59, 1706 (1987). [38] R.E. Kusner, J.A. Mann, J. Kerins, and A.J. Dahm, Phys. Rev. Lett. 73, 3113 (1994). [39] K. Zahn, R. Lenke, and G. Maret, Phys. Rev. Lett. 82, 2721 (1999). [40] A. Szabo, N. S. Ostlund, Modern Quantum Chemistry. Introduction to Advanced Electronic Structure Theory (Dover, New York, 1996). [41] R. Krishnan, J. S. Binkley, R. Seeger, J. A. Pople, J. Chem. Phys. 72, 650 (1980) [42] T. H. Dunning, Jr., J. Chem. Phys. 90, 1007 (1989). [43] S. F. Boys and F. Bernardi, Mol. Phys. 19, 553 (1970). [44] T. Helgaker, W. Klopper, H. Koch, J. Noga, J. Chem. Phys. 106, 9639 (1997). [45] D. Feller, J. Chem. Phys. 96, 6104 (1992). [46] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. flannery, Numerical Recipe in C (Cambridge University Press, Cambridge, 1996). [47] Gaussian 03, Revision C.02, Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, Jr., J. A.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; and Pople, J. A.; Gaussian, Inc., Wallingford CT, 2004. [48] J. Q. Broughton, G. H. Gilmer, and J. D. Weeks, Phys. Rev. B 25, 4651 (1982). [49] Rahman, A., Phys. Rev. 136(1964)A 405. [50] A. H.-T. Li and S. D. Chao, J. Chem. Phys. 125, 094312 (2006). [51] M. A. Spackman, J. Chem. Phys. 94, 1295 (1991). [52] Q. Wu and W. Yang, J. Chem. Phys. 116, 515 (2002). [53] A. D. Becke and E. R. Johnson, J. Chem. Phys. 122, 154104 (2005). [54] A. D. Becke, Phys. Rev. A 38, 3098 (1988). [55] N. C. Handy and A. J. Cohen, Mol. Phys. 99, 403 (2001). [56] C. Adamo and V. Barone, J. Chem. Phys. 108, 664 (1998). [57] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). [58] K. Burke, J. P. Perdew, and Y. Wang, in Electronic Density Functional Theory: Recent Progress and New Directions, Ed. J. F. Dobson, G. Vignale, and M. P. Das (Plenum, 1998). [59] J. Tao, J. P. Perdew, V. N. Staroverov, and G. E. Scuseria, Phys. Rev. Lett. 91, 146401 (2003). [60] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). [61] A. D. Boese and N. C. Handy, J. Chem. Phys. 114, 5497; see also the supp. material: EPAPS Document No. E-JCPA6-114-301111 (2001). [62] J. C. Slater, Quantum Theory of Molecular and Solids. Vol. 4: The Self-Consistent Field for Molecular and Solids (McGraw-Hill, New York, 1974). [63] J. P. Perdew, Phys. Rev. B 33, 8822 (1986). [64] S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys . 58, 1200 (1980). [65] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). [66] C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 (1988). [67] A. D. Boese, N. L. Doltsinis, N. C. Handy, and M. Sprik, J. Chem. Phys. 112, 1670 (2000). [68] A. H.-T. Li and S. D. Chao, Phys. Rev. A 73, 016701 (2006). [69] Jeremy Q. Broughton, George H. Gilmer, and John D. Weeks, Phys. Rev. B 25, (7): 4651-4669 (1982). [70] B. J. Alder, T. E. Wainwriget, Phys. Rev. 127, 359 (1962). [71] C. H. Mak, Phys. Rev. E 73, 065104(R) (2006). [72] G. Strauβ, A. Bassen, H. Zweier, H. Bertagnolli, K. Todheide, A. K. Soper, J. Turner, Phys. Rev. E 53, (1996) 3505. [73] H. Stassen, J. Mol. Struct. 464, 107 (1999). [74] S. Murad, D.J. Evan, K.E. Gubbins, W.B. Streett, D.J. Tildesley, Mol. Phys. 37 725 (1979) [75] Setzmann, U.; Wagner, W. J. Phys. Chem. Ref. Data, 20,1061 (1991). [76] J. I. Petz, J. Chem. Phys., 43, 2238 (1965). [77] R. Dawson, F. Khoury, and R. Kobyashi, A. I. Ch. E. Tl, 16, 725 (1970). [78] X. Wu, M. C. Vargas, S. Nayak, V. Lotrich and G. Scoles, J. Chem. Phys. 115, 8748 (2001). [79] W. Kohn, Y. Meir, and D. E. Makarov, Phys. Rev. Lett. 80, 4153 (1998). [80] U. Zimmerli, M. Parrinello and P. Koumoutsakos, J. Chem. Phys. 120, 2693 (2004). [81] A. D. Becke and E. R. Johnson, J. Chem. Phys. 123, 154101 (2005). [82] E. R. Johnson and A. D. Becke, J. Chem. Phys. 124, 174104 (2006). [83] X. Xu and W. A. Goddard III, Proc. Natl. Acad. Sci. 101, 2673 (2004). [84] E. J. Meijer and M. Sprik, J. Chem. Phys. 105, 8684 (1996). [85] Q. Wu and W. Yang, J. Chem. Phys. 116, 515 (2002). [86] M. Lein, J. F. Dobson, and E. K. U. Gross, J. Comp. Chem. 20, 12 (1999). [87] A. J. Misquitta, B. Jeziorski, and K. Szalewicz, Phys. Rev. Lett. 91, 033201 (2003). [88] J. F. Dobson and J. Wang, Phys. Rev. Lett. 82, 2123 (1999). [89] M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist. Phy. Rev. Lett. 92, 246401 (2004). [90] V. F. Lotrich, R. J. Bartlett, and I. Grabowski, Chem. Phys. Lett. 405, 43 (2005). [91] M. Kamiya, T. Tsuneda, and K. Hirao, J. Chem. Phys. 117, 6010 (2002). [92] T. Sato, T. Tsuneda, and K. Hirao, Mol. Phys. 103, 1151 (2005). [93] J. G. Ángyán, I. C. Gerber, A. Savin, and J. Toulouse, Phys. Rev. A. 72, 012510 (2005). [94] T. A. Wesolowski, O. Parisel, Y. Ellinger, and J. Weber, J. Phys. Chem. A 101, 7818 (1997). [95] Y. Zhang, W. Pan, and W. Yang, J. Chem. Phys. 107, 7921 (1997). [96] G. E. Scuseria and V. N. Staroverov, Chap. 24 in C. E. Dykstra, G. Frenking, K. S. Kim, and G. E. Scuseria (Ed.), Theory and Applications of Computational Chemistry : The first forty years (Elsevier, Amsterdam, 2005). [97] D. J. Lacks and R. G. Gordon, Phys. Rev. A 47, 4681 (1993). [98] E. H. Lieb and S. Oxford, Int. J. Quam. Chem. 19, 427 (1981). [99] J. Tao and J. P Perdew, J. Chem. Phys. 122, 114102 (2005). [100] C. J. Pickard and R. J. Needs, Phys. Rev. Lett. 97, 045504 (2006). | |
| dc.identifier.uri | http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/29233 | - |
| dc.description.abstract | 我們使用HF、MP2和DFT來計算甲烷與矽烷雙體作用力，HF會產生純排斥的位能曲線，且基底對HF影響並不大。在MP2的計算中，基底會影響排斥指數項(Repulsion exponent)、平衡鍵長(Equilibrium bond length)、束縛能(Binding energy)與漸進行為(Asymptotic behavior)。比較加了BSSE與未加BSSE的結果，發現加了BSSE的計算其漸進行為會隨著基底增大而逐漸收斂至實驗值，而未加BSSE的計算則不規則的產生位能曲線，顯示BSSE對MP2方法來模擬漸近行為的重要性。我們使用大範圍的exchange-correlation functionals配對的DFT計算，並比較哪種配對較接近MP2結果。
我們接下來將量子化學計算出的位能曲線利用Lennard-Jones potentials fitting出鍵距與束縛能的參數，帶進牛頓方程式，模擬甲烷液體的平衡性質與動態性質並與實驗比較，發現模擬結果與實驗比較已經相當精確，這可說明量子力學對分子間勢能計算相當具精確性。接著利用三種位能模型去模擬甲烷實驗三相圖的氣化曲線與溶化曲線，幫助我們了解氣化曲線與溶化曲線的平衡性質與動態性質。 | zh_TW |
| dc.description.abstract | We have calculated the intermolecular interaction potentials of the methane dimer at the minimum-energy D3d conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order Møller-Plesset (MP2) perturbation theory and the density functional theory (DFT). The HF calculations yield unbound potentials largely due to the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the destined C6 value from molecular polarizability calculations. The slow convergence could indicate the inefficacy of using the MP2 calculations with Gaussian-type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, and compare with MP2 results.
Next, we use a curve fitting method to fit the results of quantum chemistry calculation. We can get the force parameters of bond lengths and binding energies and then input them to solve Newton’s equations. To simulate the equilibrium properties and dynamics properties of methane liquid, we perform NVT ensemble molecular dynamics simulations. Our results well reproduce the experimental data. It demonstrates that quantum chemistry calculated intermolecular interaction is very good. Beside we simulate the equilibrium properties and dynamics properties along the vapor-liquid curve and solid-liquid curve in the phase diagram using three potential models. | en |
| dc.description.provenance | Made available in DSpace on 2021-06-13T01:03:08Z (GMT). No. of bitstreams: 1 ntu-96-R94543017-1.pdf: 3665120 bytes, checksum: 4f4c5f1c84e6a4185908e78df2012e32 (MD5) Previous issue date: 2007 | en |
| dc.description.tableofcontents | 目 錄
口試委員會審定書…………………………………………………………………………..i 致謝…………………………………………………………………………….....................ii 摘要……………………………………………………………………………………..…..iii 目錄………………………………………………………………………………………....vi 圖目錄……...…………………………………………………….........................................ix 表目錄…...……..……………………………………………………………………..…...xiii 第一章 緒論……………………………………………………………………….……1 第二章 基本理論介紹……………………………………………………………….…2 Part I 量子力學部份………………………………………………………..……….……....7 2.1 多電子問題的薛丁格方程式……………………………………………………….….7 2.1.1 單電子的薛丁格方程式……………………………………………………….…7 2.1.2 歐本-海默近似法(Born-Oppenheimer Approximation)…………….…………....8 2.1.3 庖立不相容原理(Pauli exclusion principle)…………………………………….11 2.1.4 Hartree-Fock Approximation……………………………………………………..12 2.2 Møller-Plesset微擾理論(Møller-Plesset perturbation theory)…………………….....17 2.2.1 基本數學概念……………………………………………………………....…...17 2.2.2 時間無關的微擾理論……………………………………………………….......18 2.2.3 Møller-Plesset微擾理論(Time Independent Perturbation Theory)……………....20 2.3密度泛函理論(Density Functional theory)………………………………………….….22 2.3.1 密度泛函理論的基本概念………………………………………………….......22 2.3.2 Kohn-Sham模型……………………………………………………………….…22 Part II 量子力學部份……………………………………………………..………....….....25 2.4 分子動力學的基本理論…………………………………….……………...…...….…25 2.4.1 牛頓方程式(Newton’s equation)……………………………………….......…...25 2.4.2 勢能函數(Potential energy function)………………………………..………......25 2.4.3 時步(Time step)………………………………………………………….…..…..28 2.4.4 運動方程式(Equation of motion)……………………………………………......29 2.4.5 週期性邊界條件(Periodic boundary condition)……………………………..….31 2.4.6 量測 (Measurement)…………………………………………………..………….....31 2.4.7 All-Pair method……………………………………………………………..….....…..33 2.4.8 積分演算法………………………………………………………………..…….......33 2.4.9 溫度微調(Temperature adjustment)……………..……………………….…….....…36 2.5 分子動力學的平衡性質與動態性質………………………………………………....37 2.5.1 徑向分布函數(Radial distribution function)….……………………..…………...…37 2.5.2 速度自相關函數(Velocity autocorrelation function)………...……………….....39 2.5.3 擴散係數(Diffusion coefficient)……………………………...……………........41 第三章 計算方法………………………………………………………………….…..44 Part I 量子力學部份………………………………………………………..………..........44 3.1 甲烷與矽烷的計算方法…………………………………………………………..…..44 Part II 分子動力學部份..…………………………………………………..…………..….44 3.2 分子動力學的計算詳細-2D…………………………………………………….….…45 3.2.1 二維系統……………………………………………………………………..….45 3.2.2 三維系統………………………………………………………………….……..46 第四章 計算結果與模擬結果…………………………………………………..…….51 Part I 量子力學部份………………………………………………………..……......…….51 4.1 甲烷與矽烷的計算結果…………………………………………………………..…..51 4.1.1 甲烷HF的計算結果…………………………………………………….….…..51 4.1.2 甲烷MP2的計算結果…………………………………………………..….…..52 4.1.3 甲烷DFT的計算結果………………………………………………….……....55 4.1.4 矽烷HF的計算結果………………………………………………….………...63 4.1.5 矽烷MP2的計算結果………………………………………………….….…...64 4.1.6 矽烷DFT的計算結果……………………………………………………….....66 Part II 分子動力學部份..…………………………………………………..…………..….74 4.2 二維溶化模擬…………………………………………………………………….…...74 4.3 剛性分子的分子動力學模擬-平衡性質……………………...………………..…….77 4.3.1 徑向分布函數與實驗比較…………………………………………….…….….77 4.3.2徑向分布函數－模擬沿著氣化曲線與溶化曲線……………….……….……..79 4.4剛性分子的分子動力學模擬-動態性質…………………………………….…….….90 4.4.1速度自相關函數與擴散係數－模擬沿著氣化曲線與溶化曲線……….….…..90 4.4.2擴散係數－與實驗比較……………………………………………………..…..94 第五章 結論與未來展望………………………………………………………..…….96 5.1 量子力學結論…………………………………………………………………..……..96 5.2 分子動力學結論……………………………………………………………….…...…97 5.3 未來展望……………………………………………………………………….….…..98 圖目錄 圖1 跨尺度模擬的示意圖－在不同尺度下都有與之對應的學科，這些學科經過 連接，即完成跨尺度研究。...............................................................................................2 圖2 從位能曲線出發，了解系統的平衡性質與動態性質。…………………..………...6 圖3 薛丁格方程式解出的能量曲線圖。…………………………………….…...……….10 圖4 電子在空間某點出現的機率示意圖。………………………………….…................11 圖5 hole自旋軌道數與virtual自旋軌道數示意圖。………………….…….……..…..16 圖6 Hard sphere位能曲線，在位置rc時，有一個無限高的位能障礙。…….…..............26 圖7 Shifted Lennard-Jones Potential，為Lennard-Jones Potential加上ε。…………........27 圖8 Shifted Lennard-Jones Potential且中斷點在rc。…………………………….....…..28 圖9 週期性邊界條件示意圖－中間粗黑線的框框為真實模擬的格子，其他地方為假想的格子，目的是為了讓真實格子在每一時刻都保有相同的粒子數。………….32 圖10 All-pairs method－每ㄧ個原子都必須與模擬空間中所有的原子作用。.….........33 圖11 徑向分布函數－在ㄧ球殼中所包圍的粒子數總數，在與系統的平均密度做比值，即為在此球殼的徑向分布函數值。………………………………...……........38 圖12 甲烷的D3d構型。…………………………………………………………..…...…45 圖13 4C模型－4C模型只有氫(H)參與力的計算。……………………………….....…46 圖14 4C模型fit結果。…………………………………………………………...…..…..47 圖15 5C模型fit結果。……………………………………………………………….…..48 圖16 5C模型－每個原子間都參與力的計算，故稱5C模型。……………...….............48 圖17 H構型與I構型。…………………………………………………………..…..…...49 圖18 H+I 兩種構型的fit結果。………………………………………………….….......50 圖19 甲烷的HF計算結果－可發現除了STO-3G外，其餘幾乎都已達到收斂。…..…..52 圖20 甲烷的MP2計算結果－MP2的計算呈現了鍵長、束縛能與漸近行為，其中實心為未加BSSE修正量的曲線，其曲線未隨著基底增大而系統性收斂，反之有 加BSSE修正量則有系統性的收斂。………………………………………..…….53 圖21 甲烷遠處色散能比較－發現C6+C8項的結果比C6項精確。…………………….55 圖22 (a) 固定PBE correlaion functional(甲烷)。………………………………….…..…..56 圖22 (b) 固定PW91 correlaion functional(甲烷)。…………………………………...…...57 圖22 (c) 固定VWN correlaion functional(甲烷)。…………………………………...…...57 圖22 (d) 固定VP86 correlaion functional(甲烷)。………………………………….…......58 圖22 (e) 固定LYP correlaion functional(甲烷)。………………………………….….…...58 圖23 甲烷的漸近行為分析－除了MP2外，DFT在遠處誤差較大。……….............…...59 圖24 矽烷的HF計算結果－可發現除了STO-3G外，其餘幾乎都已達到收 斂。………………………………………………………………………………….63 圖25 矽烷的MP2計算結果－MP2的計算呈現了鍵長、束縛能與漸近行為，其中實心為未加BSSE修正量的曲線，其曲線未隨著基底增大而系統性收斂，反之有 加BSSE修正量則有系統性的收斂。………………………………………….....65 圖26 矽烷遠處色散能比較發現C6+C8項的結果比C6項精確。…………………….…66 圖27 (a) 固定PBE correlaion functional(矽烷)。……………………………………….....67 圖27 (b) 固定PBE correlaion functional(矽烷)。……………………………………..…..68 圖27 (c) 固定PBE correlaion functional(矽烷)。……………………………………….....68 圖27 (d) 固定PBE correlaion functional(矽烷)。………………………………………....69 圖27 (e) 固定PBE correlaion functional(矽烷)。……………………………………….....69 圖28 矽烷的漸近行為分析除了MP2外，DFT在遠處誤差較大。……………….……...70 圖29 First order van der Waals loop。…………………………………………...…….….75 圖30 在49000-50000△t連續曝光結果，發現有相轉移存在。……………….............…75 圖31 二維溶化－在某一時間點河流現象。…………………………………………..….76 圖32 甲烷為一個似金字塔的架構，中間是碳原子，旁邊有四個氫原子與碳連接。……………………………………………………………………………….....76 圖33 不同位能曲線所計算出的徑向分布函數 與實驗比較。……….………......78 圖34 不同位能曲線所計算出的徑向分布函數 與實驗比較。………………..…78 圖35 不同位能曲線所計算出的徑向分布函數 與實驗比較。……………...…...79 圖36 甲烷的三相實驗圖[75]。………………………………………………….………...80 圖37 (a) 模擬沿著氣化曲線的徑向分布函數 圖(4C模型)。……………...……..81 圖37 (b) 模擬沿著氣化曲線的徑向分布函數 圖(4C模型)。…………………….81 圖37 (c) 模擬沿著氣化曲線的徑向分布函數 圖(4C模型)。………………...…..82 圖37 (d) 模擬沿著氣化曲線的徑向分布函數 圖(5C模型)。……………….....…82 圖37 (e) 模擬沿著氣化曲線的徑向分布函數 圖(5C模型)。………………...…..83 圖37 (f) 模擬沿著氣化曲線的徑向分布函數 圖(5C模型)。…………………….83 圖37 (g) 模擬沿著氣化曲線的徑向分布函數 圖(H+I-4C模型)。……..........…...84 圖37 (h) 模擬沿著氣化曲線的徑向分布函數 圖(H+I-4C模型)。…………….....84 圖37 (i) 模擬沿著氣化曲線的徑向分布函數 圖(H+I-4C模型)。…………...…...85 圖38 (a) 模擬沿著溶化曲線的徑向分布函數 圖(4C模型)。………………..…...85 圖38 (b) 模擬沿著溶化曲線的徑向分布函數 圖(4C模型)。………………...…..86 圖38 (c) 模擬沿著溶化曲線的徑向分布函數 圖(4C模型)。………………….....86 圖38 (d) 模擬沿著溶化曲線的徑向分布函數 圖(5C模型)。………………..…...87 圖38 (e) 模擬沿著溶化曲線的徑向分布函數 圖(5C模型)。………………..…...87 圖38 (f) 模擬沿著溶化曲線的徑向分布函數 圖(5C模型)。……………..……...88 圖38 (g) 模擬沿著溶化曲線的徑向分布函數 圖(H+I-4C模型)。.........................88 圖38 (h) 模擬沿著溶化曲線的徑向分布函數 圖(H+I-4C模型)。…………...…..89 圖38 (i) 模擬沿著溶化曲線的徑向分布函數 圖(H+I-4C模型)。……………......89 圖39 (a) 模擬沿著氣化曲線的速度自相關函數(4C模型)。………………………...…...91 圖39 (b) 模擬沿著氣化曲線的速度自相關函數(5C模型)。…………………..……..….91 圖39 (c) 模擬沿著氣化曲線的速度自相關函數(H+I-4C模型)。…………………….….92 圖40 (a) 模擬沿著溶化曲線的速度自相關函數(4C模型)。…………………..……...….92 圖40 (b) 模擬沿著氣化曲線的速度自相關函數(5C模型)。…………………..………...93 圖40 (c) 模擬沿著溶化曲線的速度自相關函數(H+I-4C模型)。……………………......93 圖41 甲烷與矽烷12個不同構型示意圖，這12個構型幾何皆對稱。………………....104 圖42 H和I方向之構型。………………………………………………………..….….105 表目錄 Table 1 利用HF和MP2的方法計算甲烷雙體的位能模型，R0為位能為0時候的距離、Rm為平衡鍵長、Eb為束縛能、ω為振動頻率，其餘為fitting的參數。…….60 Table 2 比較6-311++G(3df,3pd)的基底計算甲烷90個exchange-correlation functional的組合，MP2計算平衡鍵長的結果為3.73 ，粗黑斜線的值為誤差在10%以內的組合。…………………………………………………………………….....61 Table3 比較6-311++G(3df,3pd)的基底計算甲烷90個exchange-correlation functionals 的組合，MP2計算束縛能的結果為-0.415 kcal/mol，粗黑斜線 的值為誤差在10%以內的組合。………………..……………………………....62 Table 4 利用HF和MP2的方法計算矽烷雙體的位能模型，R0為位能為0時候的距離、Rm為平衡鍵長、Eb為束縛能、ω為振動頻率，其餘為fitting的參數。…….....71 Table 5 比較6-311++G(3df,3pd)的基底計算矽烷108個exchange-correlation functionals 的組合，MP2計算平衡鍵長的結果為4.25 ，粗黑斜線的值為誤差在10%以內的組合。………………………………………….............................................72 Table 6 比較6-311++G(3df,3pd)的基底計算矽烷108個exchange-correlation functionals 的組合，MP2計算束縛能的結果為-0.454kcal/mol，粗黑斜線的值為誤差在10%以內的組合。……………………………………………………………......73 Table 7 隨著粒子數的增加，河流出現的時間也增加，增加二為溶化假設的論證。…………………………………………………………………….………....74 Table 8 使用三種位能模型計算沿著氣化曲線與溶化曲線的擴散係數，可以討論物質存在的相。……………………………………………………….……………....94 Table 9 使用三種位能模型計算不同溫度與密度下的擴散係數，並與實驗比較。…....95 Table 10沿著氣化曲線的熱力學資料。……………………………………………….....109 Table 11沿著溶化曲線的熱力學資料。……………………………………………..……109 參考文獻………………………………………………………………………….………..99 附錄A…………………………………………………………………………………......104 附錄B………………………………………………………………………………..…....105 附錄C…………………………………………………………………………………..…109 | |
| dc.language.iso | zh-TW | |
| dc.title | 利用量子化學計算分子間勢能進行甲烷流體之分子動力學模擬 | zh_TW |
| dc.title | Molecular Dynamics Simulations of Methane Fluids using Quantum Chemistry Calculated Intermolecular Interactions | en |
| dc.type | Thesis | |
| dc.date.schoolyear | 95-2 | |
| dc.description.degree | 碩士 | |
| dc.contributor.oralexamcommittee | 張家歐,楊照彥,陳俊杉 | |
| dc.subject.keyword | Hartree-Fock (HF) 近似法,M&oslash,ller-Plesset (MP) 微擾理論,密度泛函理論 (DFT),徑向分布函數 (RDF),速度自相關函數 (VAF),擴散係數,Gaussian 03 套裝軟體, | zh_TW |
| dc.subject.keyword | Hartree-Fock (HF) approximation,M&oslash,ller-Plesset (MP) perturbation theory,Density functional theory (DFT),Radial distribution function (RDF),Velocity autocorrelation function (VAF),Diffusion coefficient,Gaussian 03 package, | en |
| dc.relation.page | 103 | |
| dc.rights.note | 有償授權 | |
| dc.date.accepted | 2007-07-25 | |
| dc.contributor.author-college | 工學院 | zh_TW |
| dc.contributor.author-dept | 應用力學研究所 | zh_TW |
| 顯示於系所單位： | 應用力學研究所 | |
文件中的檔案：
| 檔案 | 大小 | 格式 | |
|---|---|---|---|
| ntu-96-1.pdf 未授權公開取用 | 3.58 MB | Adobe PDF |
系統中的文件，除了特別指名其著作權條款之外，均受到著作權保護，並且保留所有的權利。