楊慶生， 尚軍軍， 劉 夏

(北京工業(yè)大學(xué)機械工程與應(yīng)用電子技術(shù)學(xué)院， 北京 100124)

?

多層石墨烯的粗?；肿觿恿W(xué)模型與力學(xué)性能

楊慶生， 尚軍軍， 劉 夏

(北京工業(yè)大學(xué)機械工程與應(yīng)用電子技術(shù)學(xué)院， 北京 100124)

為了有效研究多層石墨烯等石墨烯集合體的力學(xué)性能，避免由于實驗難度大、全原子計算成本高等帶來的困難，基于Tersoff勢函數(shù)提出一種可用于計算大規(guī)模石墨烯的粗?；肿觿恿W(xué)方法. 該方法計算時只需輸入原子坐標(biāo)，大大減少了模型的數(shù)據(jù)處理，提高了計算效率. 通過與全原子模型計算結(jié)果比較，本文提出的粗?；肿觿恿W(xué)方法可以準(zhǔn)確預(yù)測多層石墨烯的拉伸破壞性能，對制備高性能石墨烯纖維有一定的指導(dǎo)意義.

多層石墨烯；粗?；Ｐ?；分子動力學(xué)模擬；力學(xué)性能；破壞強度

石墨烯是由單層碳原子通過sp2共價鍵結(jié)合而成的二維六邊形網(wǎng)格結(jié)構(gòu)，其特殊的結(jié)構(gòu)使其具有優(yōu)異的力學(xué)、電學(xué)、熱學(xué)、化學(xué)等特性[1-3]，是已知的最薄、最堅硬、最有韌性的新型材料，具有廣闊的應(yīng)用前景. 近年來，大規(guī)模石墨烯的成功制備使得石墨烯增強復(fù)合材料得到了很好的發(fā)展[4-5]. 石墨烯纖維[6]、石墨烯紙[7]、石墨烯氣凝膠[8]等石墨烯集合體的問世進一步展示了石墨烯在納米科技領(lǐng)域的應(yīng)用潛力. 而石墨烯的尺寸、排布方式等對石墨烯集合體[9]及石墨烯復(fù)合材料[10]的性能有很大影響，因此充分掌握石墨烯材料的宏觀變形與微觀結(jié)構(gòu)演變之間的關(guān)聯(lián)具有重要意義. 由于尺度限制，實驗方法在觀測材料的納米尺度演化規(guī)律方面具有很大的局限性，數(shù)值模擬方法成為重要的研究手段，分子動力學(xué)是主要研究方法之一. 大規(guī)模石墨烯及其多層級結(jié)構(gòu)，如多層石墨烯等，具有跨尺度特性，全原子模擬費時費力. 粗?；肿觿恿W(xué)在模擬介觀尺度結(jié)構(gòu)的同時又保留了系統(tǒng)的分子特性，是一種有效的多尺度模擬方法，被廣泛應(yīng)用于跨尺度納米材料的研究. 如Liu等[11]采用粗?；肿觿恿W(xué)方法深入研究了碳納米管(carbon nanotube，CNT)連續(xù)纖維在加捻和拉伸條件下的力學(xué)行為及微觀結(jié)構(gòu)演變機理. 分子粗?；Ｐ鸵矔艿讲牧辖Y(jié)構(gòu)的影響，Liu等[12]還根據(jù)CNT纖維的組成結(jié)構(gòu)發(fā)展了一種粗?；肿觿恿W(xué)模型研究其拉伸性能. 很多學(xué)者采用粗粒化分子動力學(xué)方法研究有關(guān)石墨烯的問題. Titov等[13]用粗粒化分子動力學(xué)方法研究了石墨烯在脂膜中的夾層結(jié)構(gòu)；Wu等[14]模擬了烷基- 聚氧化乙烯在石墨烯表面的吸附及自組裝行為；Pandey等[15]研究了肽與石墨烯片結(jié)合的穩(wěn)定性；Cranford等[16]模擬了單層和多層石墨烯的自折疊行為；Cranford等[17]研究了加捻超長多層石墨烯帶的結(jié)構(gòu)和構(gòu)象行為；Li等[18]研究了石墨烯與脂雙層膜作用的機理；Wang等[19]重點研究了石墨烯厚度、氧化石墨烯、脂質(zhì)分子表面吸附性對石墨烯與脂雙層膜相互作用的影響. 目前，石墨烯的粗?；肿觿恿W(xué)模型主要有2類[20]，一類是Martini力場模型[21]，一類是基于應(yīng)變能守恒的降階模型[16]，由于這些方法在計算大變形問題時誤差較大且不能反映石墨烯的層間剪切響應(yīng)，因此不適用于模擬多層石墨烯的力學(xué)行為. Ruiz等[20]發(fā)展了一種可模擬石墨烯多層級結(jié)構(gòu)大變形及破壞的粗?；肿觿恿W(xué)方法. 該方法的系統(tǒng)勢能由鍵能、鍵角能、二面角能及非成鍵勢能構(gòu)成，根據(jù)應(yīng)變能守恒原理對各個參數(shù)進行校驗，得到粗?；肿觿恿W(xué)函數(shù)，其中鍵能函數(shù)形式為Morse勢. 雖然該方法可以很好地模擬多層石墨烯的力學(xué)性能，但是由于其勢能包含鍵角能和二面角能，在建模時需要模型的角度參數(shù)和二面角參數(shù)，因此當(dāng)模型復(fù)雜時，建模比較困難. 本文基于Tersoff勢函數(shù)，發(fā)展了一種可以對大規(guī)模石墨烯進行數(shù)值模擬的粗?；肿觿恿W(xué)方法，模型輸入數(shù)據(jù)簡單，對多層石墨烯的拉伸破壞進行了有效預(yù)測.

1 粗粒化模型與參數(shù)

本文以Tersoff勢函數(shù)形式表示粗粒之間的成鍵作用，Tersoff勢函數(shù)的總能量函數(shù)為

(1)

Vij(r)=[aijfR(rij)+bijfA(rij)]fC(rij)

(2)

式中：fR為排斥勢；fA為吸引勢；fC為相互作用的截斷函數(shù)；aij為階段距離，當(dāng)只考慮第一近鄰距離時，aij≈1；bij為鍵序函數(shù)，用于表征系統(tǒng)的多體效應(yīng)，它考慮了i、j原子及其周圍鍵空間分布的影響[25-26].

Lernnard-Jones勢用來描述同一分子體系內(nèi)非化學(xué)鍵作用的兩原子間的勢能，表達式為

(3)

式中：r為兩原子之間的距離；ε為能量標(biāo)度參數(shù)；σ為碰撞直徑參數(shù).

對不同的n，勢能函數(shù)中參數(shù)的取值也不相同，對于本文所研究的情況，即n=2，Tersoff勢函數(shù)中發(fā)生變化的參數(shù)及變化后的取值如表1所示. Lernnard-Jones勢函數(shù)的參數(shù)取值為ε=0.022 92 eV，σ=0.346 9 nm.

表1 粗?；Ｐ偷膮?shù)

2 數(shù)值模擬

本文所有的分子動力學(xué)計算過程采用LAMMPS程序[27]執(zhí)行. 為了驗證粗粒化模型的準(zhǔn)確性，對單層、雙層和三層的石墨烯分別進行了全原子模型和粗粒化模型的拉伸模擬，并對結(jié)果進行了比較. 模擬中單片石墨烯的尺寸為11.912 nm×20.277 nm，其中x和y方向分別沿著石墨烯的zigzag和armchair方向. 該尺寸石墨烯對應(yīng)的全原子模型中原子的個數(shù)為9 408，粗?；Ｐ椭写至５膫€數(shù)為2 352. 雙層石墨烯和三層石墨烯分別由2層和3層單片石墨烯堆疊而成，石墨烯片的層間距離為0.335 nm. 系統(tǒng)經(jīng)過能量最小化后，采用NVE微正則系綜，在300 K溫度下對模型進行充分弛豫，時間步長設(shè)置為1 fs. 模型的底端固定，沿y方向拉伸，拉伸過程中溫度保持不變，直到模型發(fā)生破壞.

本文用粗?；肿觿恿W(xué)方法對多層石墨烯的拉伸性能進行了模擬，研究了石墨烯片不同排布方式對整體性能的影響. 為此選取4種排布方式的模型進行模擬，如圖2所示，每個模型層數(shù)均為11，層間距離為0.335 nm，所有模型的整體尺寸均為45.928 nm×7.374 nm×3.685 nm. 當(dāng)一層中有多片石墨烯時，相鄰2片之間的距離為0.984 nm. A模型中，每層石墨烯都是完整的，模型中粗粒個數(shù)為37 224. B模型的石墨烯片有規(guī)律的排布，1，3，…，11層的石墨烯排布完全相同，每層有2片尺寸相同的石墨烯；2，4，…，10層的石墨烯排布完全相同，每層有3片尺寸相同的石墨烯，模型中總的粗粒個數(shù)為36 360. C模型中的石墨烯片隨機分布. 為了與B模型具有可比性，C模型粗粒個數(shù)與B模型粗粒個數(shù)相同，每層均由2或3片石墨烯構(gòu)成，石墨烯片的尺寸隨機. 為了研究連續(xù)石墨烯片對多層石墨烯性能的影響，將C模型中的2、6、11層石墨烯替換為一整片完整石墨烯，其余層數(shù)石墨烯的組成方式不變，形成D模型，模型中粗粒個數(shù)為36 576，如圖2(d)所示. 模擬過程如下：系統(tǒng)先經(jīng)過能量最小化后，在NVE微正則系綜中進行充分弛豫，弛豫溫度控制在300 K，時間步長為1 fs. 系統(tǒng)得到最優(yōu)化初始構(gòu)型后沿y方向拉伸，直到發(fā)生斷裂破壞.

3 結(jié)果與討論

3.1 粗粒化模型準(zhǔn)確性驗證

由圖3可以看出，對于單層石墨烯拉伸，粗?；Ｐ团c全原子模型結(jié)果吻合得很好，單片石墨烯的彈性模量為0.94 TPa，與實驗[28]和數(shù)值模擬[29]結(jié)果一致. 單片石墨烯的破壞強度為177 GPa，比實驗測得數(shù)值(130±10)GPa[28]偏高，這可能是由于數(shù)值模擬模型為理想化材料所引起的. 圖4所示為全原子模型和粗?；Ｐ屠炱茐倪^程，二者破壞位置都是從底部固定端開始，裂紋擴展方向與軸線大致呈45°.

圖5為雙層石墨烯全原子模型和粗粒化模型的結(jié)果比較. 從勢能和應(yīng)力與拉伸應(yīng)變的變化關(guān)系可以看出，全原子模型在第1層石墨烯發(fā)生斷裂后，第2層石墨烯又被拉伸一段位移后斷裂，2層分別斷裂區(qū)分明顯. 而對于粗?；Ｐ?，2層石墨烯斷裂時間間隔很短，模型在一層發(fā)生破壞后，整體迅速破壞，破壞模式如圖6所示. 雖然2種模型在破壞過程中裂紋擴展的形式有所不同，但整體的破壞應(yīng)變和破壞強度非常相近. 由于模型發(fā)生破壞后材料已經(jīng)失效，而破壞之前2種模型的力學(xué)性能基本一致，因此可以用粗?；Ｐ蛠泶嫒幽Ｐ?

三層石墨烯的全原子模型和粗?；Ｐ屠炱茐牡慕Y(jié)果與雙層石墨烯類似，2種模型的破壞強度和破壞應(yīng)變相差很小，全原子模型不同層的破壞相隔時間較長，而粗?；Ｐ透鲗訋缀跬瑫r斷裂，如圖7所示. 具體每層的破壞模式如圖8所示.

3.2 不同排布方式的多層石墨烯的數(shù)值模擬結(jié)果

不同排布方式的多層石墨烯力學(xué)性能差異較大，因此圖2所示模型的數(shù)值結(jié)果將在不同的圖中表示. A模型的計算結(jié)果如圖9所示. 由圖9(b)可以看出，當(dāng)多層石墨烯每層都由完整石墨烯片構(gòu)成時，其破壞強度與單層石墨烯相差不大，但是破壞應(yīng)變明顯增加. 圖10所示為B和C模型的計算結(jié)果. 從圖10(b)中可以看出，B和C模型的破壞強度遠低于A模型，這是因為B、C模型中，每層都是多片石墨烯組成，主要靠層間的范德華力來承受拉力，不同層的石墨烯片發(fā)生相對滑移引起整個結(jié)構(gòu)失效，而單個石墨烯片并沒有發(fā)生破壞. 由于B模型中石墨烯片發(fā)生相對滑移時，相互作用的表面積大于C模型，因此B模型的破壞強度要大于C模型. 為了研究石墨烯片中的成鍵作用和非成鍵作用對拉伸性能的共同影響，對D模型進行了模擬，結(jié)果如圖11所示. 由結(jié)果可以看出，D模型的破壞強度比B、C模型明顯增強.

4 結(jié)論

1) 本文提出的粗粒化分子動力學(xué)方法可以準(zhǔn)確預(yù)測大規(guī)模石墨烯的拉伸性能.

2) 當(dāng)多層石墨烯中每層的石墨烯片都不連續(xù)時，材料整體的破壞強度大大降低. 反之，含有連續(xù)石墨烯片的多層石墨烯具有較高的破壞強度.

[1] NOVOSELOV K S, GEIM A K, MOROZOV S V, et al. Electric field effect in atomically thin carbon films[J]. Science, 2004, 306(5696): 666-669.

[2] GEIM A K, NOVOSELOV K S. The rise of graphene[J]. Nature Materials, 2007, 6(3): 183-191.

[3] GEIM A K. Graphene: status and prospects[J]. Science, 2009, 324(5934): 1530-1534.

[4] ZHAO X, ZHANG Q, CHEN D, et al. Enhanced mechanical properties of graphene-based poly(vinyl alcohol) composites[J]. Macromolecules, 2010, 43(5): 2357-2363.

[5] KUILLA T, BHADRA S, YAO D, et al. Recent advances in graphene based polymer composites[J]. Progress in Polymer Science, 2010, 35(11): 1350-1375.

[6] XU Z, GAO C. Graphene chiral liquid crystals and macroscopic assembled fibres[J]. Nature Communications, 2011, 2: 571.

[7] CHEN H, MüLLER M B, GILMORE K J, et al. Mechanically strong, electrically conductive, and biocompatible graphene paper[J]. Advanced Materials, 2008, 20(18): 3557-3561.

[8] ZHANG X, SUI Z, XU B, et al. Mechanically strong and highly conductive graphene aerogel and its use as electrodes for electrochemical power sources[J]. Journal of Materials Chemistry, 2011, 21(18): 6494.

[9] CONG H P, CHEN J F, YU S H. Graphene-based macroscopic assemblies and architectures: an emerging material system[J]. Chemical Society Reviews, 2014, 43(21): 7295-7325.

[10] ZHANG M, LI Y, SU Z, et al. Recent advances in the synthesis and applications of graphene-polymer nanocomposites[J]. Polymer Chemistry, 2015, 6(34): 6107-6124.

[11] LIU X, LU W, AYALA O M, et al. Microstructural evolution of carbon nanotube fibers: deformation and strength mechanism[J]. Nanoscale, 2013, 5(5): 2002.

[12] LIU X, YANG Q. Plastic deformation and failure mechanisms of collapsed-carbon nanotube fibers by coarse-grained molecular dynamic simulations[J]. International Journal of Plasticity, 2015, 64: 104-112.

[13] TITOV A V, KRL P, PEARSON R. Sandwiched graphene-membrane superstructures[J]. ACS Nano, 2010, 4(1): 229-234.

[14] WU D, YANG X. Coarse-grained molecular simulation of self-assembly for nonionic surfactants on graphene nanostructures[J]. The Journal of Physical Chemistry B, 2012, 116(39): 12048-12056.

[15] PANDEY R B, KUANG Z, FARMER B L, et al. Stability of peptide (P1 and P2) binding to a graphene sheet via an all-atom to all-residue coarse-grained approach[J]. Soft Matter, 2012, 8(35): 9101.

[16] CRANFORD S, SEN D, BUEHLER M J. Meso-origami: folding multilayer graphene sheets[J]. Applied Physics Letters, 2009, 95(12): 123121.

[17] CRANFORD S, BUEHLER M J. Twisted and coiled ultralong multilayer graphene ribbons[J]. Modelling and Simulation in Materials Science and Engineering, 2011, 19(5): 785-791

[18] LI Y, YUAN H, VON DEM BUSSCHE A, et al. Graphene microsheets enter cells through spontaneous membrane penetration at edge asperities and corner sites[J]. Proceedings of the National Academy of Sciences, 2013, 110(30): 12295-12300.

[19] WANG J, WEI Y, SHI X, et al. Cellular entry of graphene nanosheets: the role of thickness, oxidation and surface adsorption[J]. RSC Advances, 2013, 3(36): 15776-15782.

[20] RUIZ L, XIA W, MENG Z, et al. A coarse-grained model for the mechanical behavior of multi-layer graphene[J]. Carbon, 2015, 82: 103-115.

[21] MARRINK S J, RISSELADA H J, YEFIMOV S, et al. The MARTINI force field: coarse grained model for biomolecular simulations[J]. The Journal of Physical Chemistry B, 2007, 111(27): 7812-7824.

[22] MONTICELLI L, KANDASAMY S K, PERIOLE X, et al. The MARTINI coarse-grained force field: extension to proteins[J]. Journal of Chemical Theory and Computation, 2008, 4(5): 819-834.

[23] HSU D D, XIA W, ARTURO S G, et al. Systematic method for thermomechanically consistent coarse-graining: a universal model for methacrylate-based polymers[J]. Journal of Chemical Theory and Computation, 2014, 10(6): 2514-2527.

[24] XIA W, MISHRA S, KETEN S. Substrate vs. free surface: competing effects on the glass transition of polymer thin films[J]. Polymer, 2013, 54(21): 5942-5951.

[25] TERSOFF J. Modeling solid-state chemistry-interatomic potentials for multicomponent systems[J]. Physical Review B, 1989, 39(8): 5566-5568.

[26] TERSOFF J. Erratum: modeling solid-state chemistry: interatomic potentials for multicomponent systems[J]. Physical Review B, 1990, 41(5): 3248.

[27] PLIMPTON S. Fast parallel algorithms for short-range molecular-dynamics[J]. Journal of Computational Physics, 1995, 117(1): 1-19.

[28] LEE C, WEI X, KYSAR J W, et al. Measurement of the elastic properties and intrinsic strength of monolayer graphene[J]. Science, 2008, 321(5887): 385-388.

[29] JIANG J, WANG J, LI B. Young’s modulus of graphene: a molecular dynamics study[J]. Physical Review B, 2009, 80(11): 113405.

(責(zé)任編輯 呂小紅)

Coarse-grained Molecular Dynamics Model and Mechanical Properties of Multi-layer Graphene

YANG Qingsheng, SHANG Junjun, LIU Xia

(College of Mechanical Engineering and Applied Electronics Technology, Beijing University of Technology,Beijing 100124, China)

To investigate the behavior of multi-layer graphene (MLG) and other graphene assemblies effectively, a coarse-grained molecular dynamics (CG-MD) method was developed, which conquers the limitations of experimental observations and full atom simulation to some extent. The CG-MD method based on Tersoff potential was used to simulate large scale graphene. The coordinates of the model were needed in the process of computation, which greatly reduced the time of modeling. The comparison between the results of full atom model and CG-MD model was made. It is proved that the CG-MD model can predict the tensile behavior of MLG accurately and it has great significance to producing graphene fiber with high performance.

multi-layer graphene；coarse-grained model; molecular dynamics simulation；mechanical properties；failure strength

2016- 08- 08

國家自然科學(xué)基金資助項目(11472020)

楊慶生(1962—)， 男， 教授， 博士生導(dǎo)師， 主要從事新型材料與結(jié)構(gòu)的力學(xué)問題方面的研究， E-mail: qsyang@bjut.edu.cn

TB 303

A

0254-0037(2016)12-1798-07

10.11936/bjutxb2016080021