CN106934137A - 一种分析石墨烯组装体的粗粒化分子动力学方法 - Google Patents

一种分析石墨烯组装体的粗粒化分子动力学方法 Download PDF

Info

Publication number
CN106934137A
CN106934137A CN201710125683.7A CN201710125683A CN106934137A CN 106934137 A CN106934137 A CN 106934137A CN 201710125683 A CN201710125683 A CN 201710125683A CN 106934137 A CN106934137 A CN 106934137A
Authority
CN
China
Prior art keywords
graphene
coarse
assembly
molecular dynamics
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710125683.7A
Other languages
English (en)
Inventor
刘夏
杨庆生
尚军军
张雪娇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201710125683.7A priority Critical patent/CN106934137A/zh
Publication of CN106934137A publication Critical patent/CN106934137A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种分析石墨烯组装体的粗粒化分子动力学方法,属于计算机分子模拟技术领域。该方法包括以下步骤:a)构建单层石墨烯的粗粒化模型;b)根据组装形式构建石墨烯组装体的粗粒化模型;c)设置粗粒之间的相互作用;d)进行粗粒化模拟计算。通过本发明可以简便地构建石墨烯组装体计算模型,并在弛豫过程中准确地得到实际组装体模型。所采用的Tersoff势函数与Lennard‑Jones势函数可以很好地表征组装体内化学键的断裂与形成,准确地预测石墨烯的自组装行为。所采用的Tersoff势函数与Lennard‑Jones势函数可以准确地预测石墨烯组装体的力学行为。

Description

一种分析石墨烯组装体的粗粒化分子动力学方法
技术领域
本发明属于计算机分子模拟技术领域,具体涉及一种研究石墨烯组装体的粗粒化分子动力学方法。
背景技术
近几年,石墨烯组装体成为非常热的话题,吸引了众多的实验、理论和模拟研究,在许多领域都将引起新的技术革命。石墨烯组装体是指以单层石墨烯为基本组装单元,通过共价连接或者氢键、静电作用、π-π作用、范德华力等非共价作用,形成更高层次的有序功能体。石墨烯组装体是实现石墨烯从纳米材料转变为宏观功能材料的关键一步。通过调控石墨烯组装体的形貌和结构,不仅能够更好地实现石墨烯优异的力学、电学和热学等性能,而且还能拓展新的功能和应用。因此,研究石墨烯组装体具有重要的意义。
尽管在近期,提出了许多研究石墨烯的粗粒化分子动力学方法,但都是采用粗粒-弹簧粗粒化模型。该模型是将一个分子链简化为一串粗粒,粗粒间由弹簧连接,适用于具有极大长径比的聚合物分子链或纳米纤维。然而,在构建粗粒-弹簧模型时,除了需要确定粗粒的坐标,还需指定键连信息(两个粗粒间)、键角(三个粗粒间)信息以及二面角信息(四个粗粒间),建模过程非常复杂,不适于复杂体系的研究,且无法预测组装体的准确组装结构。因此,我们提供了一种新的研究单层石墨烯及石墨烯组装体的分子动力学方法,该方法建模过程简单,仅需确定单层石墨烯上各粗粒的初始坐标,且选用的Tersoff多体势与Lennard-Jones远程势的叠加形式可以准确地预测石墨烯组装体的结构。
发明内容
本发明要解决的技术问题是提供一种分析石墨烯组装体的粗粒化分子动力学方法。
为解决上述技术问题,本发明采用以下技术方案:
本发明公开一种分析石墨烯组装体的粗粒化分子动力学方法,该方法包括以下步骤:a)构建单层石墨烯的粗粒化模型;b)根据组装形式构建石墨烯组装体的粗粒化模型;c)设置粗粒之间的相互作用;d)进行粗粒化模拟计算。
所述石墨烯粗粒分子包含4个碳原子,石墨烯粗粒分子的初始位置为4个碳原子的质量中心,石墨烯粗粒分子的质量为一个碳原子的4倍,石墨烯粗粒分子间的键长为全原子碳-碳键的2倍。
建模时只需要石墨烯粗粒分子的坐标,不需要指定键连、键角与二面角等连接信息。
所述石墨烯粗粒分子间相互作用的力场形式为Tersoff多体势与Lennard-Jones远程势的叠加。
具体地说,该方法包含以下步骤:
1)在分子建模软件Accelrys Material Studio中建立单层石墨烯模型,石墨烯粗粒分子以六元环形式连接,石墨烯粗粒分子间距离为0.284nm。
2)构建石墨烯组装体的初始结构,并将所有石墨烯粗粒分子的坐标输出到坐标文件。
3)采用LAMMPS程序对石墨烯组装体进行计算,在计算命令中指定对势函数类型为Tersoff和LJ/cut,命令如下:
pair_style hybrid/overlay tersofflj/cut 10
pair_coeff**tersoffSiC.tersoffC
pair_coeff**lj/cut 0.022923.469
函数中各项参数在Graphene.tersoff中。格式如下:
C C C 3.01.00.0380494.3484-.57058.727510.000000157241.105951386.84.10.61.743955574.4
本发明的有益效果如下:
通过本发明所述的一种研究石墨烯组装体的粗粒化分子动力学方法,可以简便地构建石墨烯组装体计算模型,并在弛豫过程中准确地得到实际组装体模型。
所采用的Tersoff势函数与Lennard-Jones势函数可以很好地表征组装体内化学键的断裂与形成,准确地预测石墨烯的自组装行为。
所采用的Tersoff势函数与Lennard-Jones势函数可以准确地预测石墨烯组装体的力学行为。
附图说明
图1为双层石墨烯粗粒化模型,其中a)为初始结构;b)为拉伸断裂。
图2为无定型块体结构,其中a)为能量弛豫前;b)为能量弛豫后。
图3为单层石墨烯结构,其中a)为能量弛豫前;b)为能量弛豫后。
图4为石墨烯分子动力学模型的粗粒化映射;a)为全原子坐标;b)为粗粒坐标。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。
实施例1
双层石墨烯拉伸破坏的粗粒化分子动力学模拟
建立长度为20nm、宽度为10nm的双层石墨烯粗粒化模型,如图1中的a)所示,将初始结构的坐标保存在坐标文件。用LAMMPS程序对石墨烯粗粒化模型进行充分的能量弛豫,使其得到最优构型。将双层石墨烯的最底端1nm固定,并对最顶端1nm处的所有粗粒施加0.5nN的力,最终将双层石墨烯拉断,如图b)所示。
实施例2
石墨烯气凝胶的粗粒化分子动力学模型
在Accelrys Material Studio软件中建立长度与宽度均为5nm的单层石墨烯,用10片石墨烯构建密度为0.2g/cm3的无定型块体结构,如图2中的a)所示。然后,将所有原子坐标导出,并乘以2,保存为LAMMPS可以读取的坐标文件。采用LAMMPS程序,对无定型结构施加静水压力,石墨烯片自组装为气凝胶结构,如图b)所示。
实施例3
石墨烯三维搭接的粗粒化分子动力学模拟
在Accelrys Material Studio软件中建立3片长度与宽度均为5nm的单层石墨烯,并按图3中的a)摆放在一起。此时,每片石墨烯相互独立,相互之间没有共价键。然后,将所有原子坐标导出,并乘以2,保存为LAMMPS可以读取的坐标文件。采用LAMMPS程序,对系统进行能量弛豫。在能量弛豫过程中,三片石墨烯之间自动形成了共价连接,自组装为一个三维结构。
图4为石墨烯分子动力学模型的粗粒化映射;a)为全原子坐标;b)为粗粒坐标。

Claims (8)

1.一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:该方法包括以下步骤:a)构建单层石墨烯的粗粒化模型;b)根据组装形式构建石墨烯组装体的粗粒化模型;c)设置粗粒之间的相互作用;d)进行粗粒化模拟计算。
2.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:所述石墨烯粗粒分子包含4个碳原子,石墨烯粗粒分子的初始位置为4个碳原子的质量中心,石墨烯粗粒分子的质量为一个碳原子的4倍,石墨烯粗粒分子间的键长为全原子碳-碳键的2倍。
3.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:建模时只需要石墨烯粗粒分子的坐标,不需要指定键连、键角与二面角等连接信息。
4.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:所述石墨烯粗粒分子间相互作用的力场形式为Tersoff多体势与Lennard-Jones远程势的叠加。
5.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:具体地说,该方法包含以下步骤:
1)在分子建模软件Accelrys Material Studio中建立单层石墨烯模型,石墨烯粗粒分子以六元环形式连接,石墨烯粗粒分子间距离为0.284nm;
2)构建石墨烯组装体的初始结构,并将所有石墨烯粗粒分子的坐标输出到坐标文件;
3)采用LAMMPS程序对石墨烯组装体进行计算,在计算命令中指定对势函数类型为Tersoff和LJ/cut,命令如下:
pair_style hybrid/overlay tersofflj/cut 10
pair_coeff**tersoffSiC.tersoffC
pair_coeff**lj/cut 0.022923.469
函数中各项参数在Graphene.tersoff中;格式如下:
C C C 3.01.00.0380494.3484-.57058.727510.000000157241.105951386.84.10.61.743955574.4。
6.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:
双层石墨烯拉伸破坏的粗粒化分子动力学模拟
建立长度为20nm、宽度为10nm的双层石墨烯粗粒化模型,将初始结构的坐标保存在坐标文件;用LAMMPS程序对石墨烯粗粒化模型进行充分的能量弛豫,使其得到最优构型;将双层石墨烯的最底端1nm固定,并对最顶端1nm处的所有粗粒施加0.5nN的力,最终将双层石墨烯拉断。
7.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:
石墨烯气凝胶的粗粒化分子动力学模型
在Accelrys Material Studio软件中建立长度与宽度均为5nm的单层石墨烯,用10片石墨烯构建密度为0.2g/cm3的无定型块体结构;然后,将所有原子坐标导出,并乘以2,保存为LAMMPS可以读取的坐标文件;采用LAMMPS程序,对无定型结构施加静水压力,石墨烯片自组装为气凝胶结构。
8.根据权利要求1所述的一种分析石墨烯组装体的粗粒化分子动力学方法,其特征在于:
石墨烯三维搭接的粗粒化分子动力学模拟
在Accelrys Material Studio软件中建立3片长度与宽度均为5nm的单层石墨烯,摆放在一起;此时,每片石墨烯相互独立,相互之间没有共价键;然后,将所有原子坐标导出,并乘以2,保存为LAMMPS读取的坐标文件;采用LAMMPS程序,对系统进行能量弛豫;在能量弛豫过程中,三片石墨烯之间自动形成了共价连接,自组装为一个三维结构。
CN201710125683.7A 2017-03-05 2017-03-05 一种分析石墨烯组装体的粗粒化分子动力学方法 Pending CN106934137A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710125683.7A CN106934137A (zh) 2017-03-05 2017-03-05 一种分析石墨烯组装体的粗粒化分子动力学方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710125683.7A CN106934137A (zh) 2017-03-05 2017-03-05 一种分析石墨烯组装体的粗粒化分子动力学方法

Publications (1)

Publication Number Publication Date
CN106934137A true CN106934137A (zh) 2017-07-07

Family

ID=59423305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710125683.7A Pending CN106934137A (zh) 2017-03-05 2017-03-05 一种分析石墨烯组装体的粗粒化分子动力学方法

Country Status (1)

Country Link
CN (1) CN106934137A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416119A (zh) * 2018-02-12 2018-08-17 大连理工大学 一种基于分子动力学模拟测试碳纳米管力学性能的方法
CN109378041A (zh) * 2018-12-03 2019-02-22 中国石油大学(北京) 介尺度模拟石油分子粗粒化模型及其构建方法和预测方法
CN110070918A (zh) * 2019-04-02 2019-07-30 重庆邮电大学 基于分子间相互作用的粗粒化方法
CN110442887A (zh) * 2018-05-03 2019-11-12 中国科学院化学研究所 一种采用系统粗粒化分子动力学研究重复聚酰亚胺体系热力学性质的模型构建方法
CN111081323A (zh) * 2019-12-19 2020-04-28 哈尔滨工业大学 一种基于Tersoff力场的石墨烯多级粗粒化方法
CN112768007A (zh) * 2021-01-29 2021-05-07 太原理工大学 一种基于Martini力场的色素炭黑粗粒化模型建立方法
CN112992279A (zh) * 2021-03-02 2021-06-18 南京理工大学 基于分子动力学的氢化石墨烯纳米盒的模拟方法
CN113035283A (zh) * 2021-01-12 2021-06-25 南京理工大学 基于分子动力学的多层富勒烯单向压缩的模拟方法
CN114864019A (zh) * 2022-06-02 2022-08-05 江南大学 包含扭转行为的碳纳米管粗粒化势函数构建方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101489197A (zh) * 2009-01-22 2009-07-22 中国科学院计算技术研究所 一种电子设备间数据传输的方法和系统
CN101651781A (zh) * 2008-08-11 2010-02-17 奥林巴斯映像株式会社 便携终端装置和外部设备
CN101877021A (zh) * 2009-12-28 2010-11-03 中国农业大学 动力学分析、三维造型和有限元分析软件间数据自动转换的方法
CN104598679A (zh) * 2015-01-14 2015-05-06 中国石油大学(华东) 一种自组装纳米卷状结构的制备方法及其模拟验证方法
US9362919B1 (en) * 2014-12-22 2016-06-07 University Of Notre Dame Du Lac Devices for utilizing symFETs for low-power information processing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101651781A (zh) * 2008-08-11 2010-02-17 奥林巴斯映像株式会社 便携终端装置和外部设备
CN101489197A (zh) * 2009-01-22 2009-07-22 中国科学院计算技术研究所 一种电子设备间数据传输的方法和系统
CN101877021A (zh) * 2009-12-28 2010-11-03 中国农业大学 动力学分析、三维造型和有限元分析软件间数据自动转换的方法
CN101877021B (zh) * 2009-12-28 2012-06-20 中国农业大学 动力学分析、三维造型和有限元分析软件间数据自动转换的方法
US9362919B1 (en) * 2014-12-22 2016-06-07 University Of Notre Dame Du Lac Devices for utilizing symFETs for low-power information processing
CN104598679A (zh) * 2015-01-14 2015-05-06 中国石油大学(华东) 一种自组装纳米卷状结构的制备方法及其模拟验证方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘鑫: "《三维石墨烯的水热法制备及其性能研究》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *
杨庆生等: "《多层石墨烯的粗粒化分子动力学模型与力学性能》", 《北京工业大学学报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416119A (zh) * 2018-02-12 2018-08-17 大连理工大学 一种基于分子动力学模拟测试碳纳米管力学性能的方法
CN110442887B (zh) * 2018-05-03 2021-04-06 中国科学院化学研究所 一种研究重复聚酰亚胺体系热力学性质的模型构建方法
CN110442887A (zh) * 2018-05-03 2019-11-12 中国科学院化学研究所 一种采用系统粗粒化分子动力学研究重复聚酰亚胺体系热力学性质的模型构建方法
CN109378041B (zh) * 2018-12-03 2021-01-08 中国石油大学(北京) 介尺度模拟石油分子粗粒化模型及其构建方法和预测方法
CN109378041A (zh) * 2018-12-03 2019-02-22 中国石油大学(北京) 介尺度模拟石油分子粗粒化模型及其构建方法和预测方法
CN110070918A (zh) * 2019-04-02 2019-07-30 重庆邮电大学 基于分子间相互作用的粗粒化方法
CN110070918B (zh) * 2019-04-02 2022-12-27 重庆邮电大学 基于分子间相互作用的粗粒化方法
CN111081323A (zh) * 2019-12-19 2020-04-28 哈尔滨工业大学 一种基于Tersoff力场的石墨烯多级粗粒化方法
CN113035283A (zh) * 2021-01-12 2021-06-25 南京理工大学 基于分子动力学的多层富勒烯单向压缩的模拟方法
CN113035283B (zh) * 2021-01-12 2022-08-23 南京理工大学 基于分子动力学的多层富勒烯单向压缩的模拟方法
CN112768007A (zh) * 2021-01-29 2021-05-07 太原理工大学 一种基于Martini力场的色素炭黑粗粒化模型建立方法
CN112992279A (zh) * 2021-03-02 2021-06-18 南京理工大学 基于分子动力学的氢化石墨烯纳米盒的模拟方法
CN112992279B (zh) * 2021-03-02 2022-08-26 南京理工大学 基于分子动力学的氢化石墨烯纳米盒的模拟方法
CN114864019A (zh) * 2022-06-02 2022-08-05 江南大学 包含扭转行为的碳纳米管粗粒化势函数构建方法

Similar Documents

Publication Publication Date Title
CN106934137A (zh) 一种分析石墨烯组装体的粗粒化分子动力学方法
Liu et al. Molecular dynamics simulation on interfacial mechanical properties of polymer nanocomposites with wrinkled graphene
Zhang et al. Nanomechanics of carbon honeycomb cellular structures
Hu et al. Multi-scale numerical simulations on piezoresistivity of CNT/polymer nanocomposites
Moon et al. Interfacial strengthening between graphene and polymer through Stone-Thrower-Wales defects: Ab initio and molecular dynamics simulations
Hu et al. Prediction of thermal expansion properties of carbon nanotubes using molecular dynamics simulations
Chandra et al. Pullout strength of graphene and carbon nanotube/epoxy composites
Baimova et al. From flat graphene to bulk carbon nanostructures
Yang et al. Sensitivity-Tunable strain sensors based on carbon nanotube@ carbon nanocoil hybrid networks
Wang et al. Enhanced mechanical properties of carbon nanotube networks by mobile and discrete binders
Jensen et al. Simulating the effects of carbon nanotube continuity and interfacial bonding on composite strength and stiffness
Goel et al. Analysis of temperature, helicity and size effect on the mechanical properties of carbon nanotubes using molecular dynamics simulation
Kumar et al. A rod model for three dimensional deformations of single-walled carbon nanotubes
Baric et al. Compaction-induced restructuring of aggregated nanoparticle films using the discrete element method
Liao et al. Understanding the role of self-adhesion in crumpling behaviors of sheet macromolecules
Liao et al. Crumpling defective graphene sheets
Tonel et al. Interactions of graphene derivatives with glutamate-neurotransmitter: a parallel first principles-docking investigation
Luo et al. Highly stretchable conductors based on expanded graphite macroconfined in tubular rubber
Alian et al. Hybrid molecular dynamics–finite element simulations of the elastic behavior of polycrystalline graphene
Chen et al. Does thermal percolation exist in graphene-reinforced polymer composites? A molecular dynamics answer
Fernández Sánchez-Romate et al. Electrical transport mechanisms in graphene nanoplatelet doped polydimethylsiloxane and application to ultrasensitive temperature sensors
Akbarzadeh et al. H2 adsorption on Ag‐nanocluster/single‐walled carbon nanotube composites: A molecular dynamics study on the effects of nanocluster size, diameter, and chirality of nanotube
Jam et al. Characterizing elastic properties of carbon nanotube‐based composites by using an equivalent fiber
Niu et al. Mechanical and thermal properties of carbon nanotubes in carbon nanotube fibers under tension–torsion loading
Nasiri et al. Multilayer structures of graphene and Pt nanoparticles: A multiscale computational study

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20170707

RJ01 Rejection of invention patent application after publication