摘要
隨著日新月異的時代演變，本篇論文首先指出科技的大方向—輕薄短小，而分子動力學剛好就是ㄧ種研究小尺度的好方法，接著介紹一些微觀尺度的演進，進一步對未來的我們，在判斷工程相關的事物時，能有簡單的方向和想法，並且介紹聚丙烯材料的相關特性及應用。
接著對於分子動力學、聚丙烯、溶液性質、擴散性質，將相關研究結果以簡單的敘述做文獻回顧，再將研究的方法、勢能函式列在第三章，其中包括模擬流程圖、初始化方法、演算法、力的計算方式、週期性邊界及烷類和聚丙烯的勢能參數。
利用分子動力學模擬聚丙烯溶液的末端末端距離、迴旋半徑和動態性質—擴散係數，分別用濃度和鏈長作為操縱變因，討論模擬結果。由研究結果發現，當重量百分濃度增加時，擴散物質的體積與濃度的關係平緩，故擴散係數與濃度的關係亦平緩；當鏈長增加時，則聚丙烯所佔有的體積變大(亦可由末端末端距離和迴旋半徑看出)，因而前進時容易受阻，故擴散係數下降。
關鍵字：分子動力學、聚丙烯、擴散係數

This study pointed out the direction of science technology, small, and molecular dynamics was one of the good methods to research small scale world. We introduce the developments of small scale to help us to solve engineering problems that can make us have simple directions and ideas. This study introduced the features and applications of polypropylene simply.
There were easy words for related researches about molecular dynamics、polypropylene、solution properties、diffusion properties. In chapter 3, we showed study methods and potential functions including a simulation flow chat、initial methods、algorithm、calculation method of force、periodic boundary condition and potential parameter of alkane series and polypropylene.
Molecular dynamics simulation was used to measure end-to-end distance、radius of gyration and diffusion coefficient of polypropylene solution. The chain length and concentration were control variances to discuss simulation results. Because the change for diffusion substance was small, the diffusion coefficient changed little with weight percent concentration, too. The diffusion coefficient decreased with increasing chain length, because the raising size of polypropylene hindered substance moving.
Key Word：Molecular Dynamics、Polypropylene、Diffusion Coefficient

參考文獻
[1] James Elliott，Department of Materials Science & Metallurgy，University of Cambridge，http://www.msm.cam.ac.uk/Department/DeptInfo/StaffProfiles/ResearchFigs/Elliott.jpg
[2] 第三媒體，http://idoimg.3mt.com.cn/article/upload/200801/0801291710972031.jpg，2008
[3] Jymingyue，互動百科，http://a2.att.hudong.com/73/78/01300000231569122174780737423_s.jpg
[4] 許鴻德，□塑膠擔當未來汽車先進材料大任（上）□，http://icar.yam.com/photo_data/20071210000692/XL20071210000692_3518.jpg，10 October 2007
[5] 翔營時代電子科技有限公司，http://www.amwink.com/Product/image/6336235787792187501470512362.JPG
[6] V8rik，維基百科，http://upload.wikimedia.org/wikipedia/en/b/bf/Atacticpolymer.gif
[7] Minestrone Soup，維基百科，http://en.wikipedia.org/wiki/Image:Isotactic_polymer.png
[8] V8rik，維基百科，http://upload.wikimedia.org/wikipedia/en/6/60/Syndiotacticpolymer.gif
[9] 易豐達膠帶製品有限公司，互動百科，http://www.hudong.com/wiki/%E8%81%9A%E4%B8%99%E7%83%AF
[10] Harutun G. Karian ，HANDBOOK OF POLYPROPY1ENE AND POLYPROPYLENE COMPOSITES，RheTech. Inc，2003
[11] Chen，University of California, Los Angeles，維基百科，http://zh.wikipedia.org/w/index.php?title=%E5%88%86%E5%AD%90%E5%8A%A8%E5%8A%9B%E5%AD%A6&variant=zh-tw
[12] M. Meunier，□Diffusion coefficients of small gas molecules in amorphous cis-1,4-polybutadiene estimated by molecular dynamics simulations □，J. Chem. Phys. ，Vol. 123， pp.134906-1~ 134906-7 ，2005
[13] Mechthild Franke, Nicolaus Rehse， “Nucleation of branches in elastomeric polypropylene ”，Polymer， Vol. 49， Issue 20， pp. 4328-4331， 2008
[14] Fukuda M, Kuwajima S，J.Chem Phys，Vol.108(7)， pp.3001-3009，FEB 15 1998

[15] Stelios J. Antoniads, Christina T. Samara, and Doros N. Theodorou， Macromolecules， Vol.32, pp.8635 - 8644，1999
[16] Georgia Evangelia Logotheti and Doros N. Theodorou, □Segmental and Chain Dynamics of isotactic Polypropylene Melts□， Macromolecules，Vol. 40，pp. 2235-2245，2007
[17] Ronaldo Giro and Mar□lia J. Caldas，□Calcium depositionon poly(para-phenylene vinylene) : Molecular dynamics simulations□，PHYSICAL REVIEW B，Vol. 76， pp. 161303，2007
[18] M. Druchok, B. Hribar-Lee, H. Krienke, V. Vlachy，□A Molecular Dynamics study of short-chain polyelectrolytes in explicit water: Toward the origin of ion-specific effects□，Chemical Physics Letter，Vol. 450，pp. 281 - 285，2008
[19] Rong Zhang, Dunsheng Zheng, Yufang Pan, Sanlai Luo, Wenjuan Wu, Haoran Li，□All-atom simulation and excess properties study on intermolecular interactions of amide–water system□， Journal of Molecular Structure，Vol. 875，pp. 96–100，2008
[20] 王子瑜、曹□光，中央大學化學工程與材料工程所，□布朗運動、郎之萬方程式、與布朗動力學□，物理雙月刊，（廿七卷三期）2005年6月

[21] Mathias Putz,John G. Curro and Gary S. Grest，□Self-consistent　integral　equation　theory　for　polyolefins:　Comparison　to　molecular　dynamics　simulations　and　x-ray　scattering□，J.Chem.Phys.，Vol. 114，No.6，pp.2847 - 2860，8 February 2001
[22] Lee SH, Lee H, Pak H, et al.，Bulletin of The Korean Chemical Society，Vol:17，Issue:8，pp. 735-744 ，1996
[23] David Heine and David T. Wu, John G. Curro and Gary S. Grest，□Role of intramolecular energy on polyolefin miscibility : Isotactic polypropylene polyethylene blends□，J.Chem.Phys.，Vol. 118，No.2，pp.914 - 924，8 January 2003
[24] 葛林•費雪班（Glenn Fishbine）著╱劉世平譯：《奈米商機》（台灣培生教育出版（股）公司）
[25] Rakwoo Chang and Arun Yethiraj，□Solvent effects on the collapse dynamics of polymers□，J.Chem.Phys.，Vol. 114， No. 17， pp. 7688–7699，1 May 2001
[26] J. M. Haile，Molecular Dynamics Simulation: Elementary Methods，Wiley，1997
[27] Junfang Zhang, S K Choi, Mayela Rivero，□Molecular dynamics simulation of polyacrylamides in potassium montmorillonite clay hydrates□，JOURNAL OF PHYSICS B: ATOMIC, MOLECULAR AND OPTICAL PHYSICS，Vol. 40， pp. 545–553，19 January 2007
[28] Georgia Tsolou, Vlasis G. Mavrantzas, Doros N. Theodorou，□Detailed Atomistic Molecular Dynamics Simulation of cis-1,4-Poly(butadiene)□，Macromolecules，Vol. 38 ，No. 4， pp1478 -1492，2005
[29] R. Byron Bird, Warren E. Stewart, Edwin N. Lightfoot，Transport Phenomena，Wiley，2006，pp.513-538、pp.866
[30] Hans Christian □ttinger，□Translational diffusivity from the Zimm model□，J.Chem.Phys.，Vol. 87， pp. 3156–3165，1987
[31] C. Bennemann, J. Baschnagel, W. and K. Binder，□Molecular-dynamics simulation next term of a glassy polymer melt: previous term Rouse model next term and cage effect□，Computational and Theoretical Polymer Science，Vol. 9， pp. 217–226，15 February 2000
[32] Piyush G. Gigras and Bamin Khomami，□Adaptive configuration fields: a new multi-scale simulation next term technique for previous term reptation next term-based previous term models next term with a stochastic strain measure and local variations of life span distribution□，Journal of Non-Newtonian Fluid Mechanics，Vol. 108， pp. 99–122，30 December 2002
[33] William G. Hoover，□Nose-Hoover thermostat□，PHYSICAL REVIEW A，Vol. 31， pp. 1695-1697，1985
[34] Luis F. Rull and Juan J. Morales，□Molecular-dynamics pressure thermostat□，PHYSICAL REVIEW A，Vol. 38， pp. 4309-4312，1988
[35] S. J. V. Frankland, V. M. Harik, G. M. Odegard, D.W. Brenner, T. S. Gates，□Molecular-dynamics pressure thermostat□，Composites Scienceand Technology，Vol. 63， pp. 1655-1661，2003
[36] Thomas R. Cuthbert, Norman J. Wagner, Michael E. Paulaitis, Giovanni Murgia and Bruno D’Aguanno，□Molecular Dynamics Simulation of Penetrant Diffusionin Amorphous Polypropylene: Diffusion Mechanisms and Simulation Size Effects□，Macromolecules，Vol. 32， pp. 5017-5028，1999
[37] S. G. Lambrakos, J. P. Boris, I. Chandrasekhar, and B. Gaber，□A Vectorized Near-Neighbors Algorithm of Order N for Molecular Dynamic Simulation□，Ann. NY Acad. Sci.，Vol. 85，pp 482 ，1986
[38] S. Gupta，□Vectorization of Molecular Dynamic Simulation for Fluids of Nonspherical Molecules□，Comp. Phys. Commun.，Vol. 48， pp. 197，1988
[39] M. P. Allen, D. J. Tildesley，Computer Simulation of Liquid，Oxford，pp.228-230，1991
[40] Richard J. Sadus，Molecular Simulation of Fluids，Elsevier pp.310-313，1999，
[41] 林明獻《矽晶半導體材料技術》（全華科技圖書股份有限公司），頁2-19
[42] D. C. Rapaport，The Art of Molecular Dynamics Simulation，Cambridge，pp. 278-285，2004
[43] B. Minisini, F. Tsobnang，□Molecular dynamics study of specific interactions in grafted polypropylene organomodified clay nanocomposite □，Composites : Part A ，Vol. 36，pp. 539–544，2005
[44] 林弘凡，運用分子動力學與平行運算於奈米流場之研究，碩士論文，國立清華大學，2004
[45] 洪玉蕙，立體正規性對聚丙烯在溶液中的影響：分子動力學模擬，碩士論文，國立清華大學，2007
[46] 黃丘宏，應用分子動力學計算奈米尺度下乙醇及聚丙烯溶液物性之研究，碩士論文，國立清華大學，2006
[47] 朱敬平，材料加工輸送現象，課程講義，國立台灣科技大學高分子工程系，2008