| 查看: 4770 | 回复: 109 | ||||
[交流]
用GPUMD计算的石墨烯拉伸曲线
|
||||
|
最近在写程序，交流一下测试结果。我的目标是写一个在某些方面比lammp好用的小巧的程序。我的程序叫GPUMD，意思是graphics processing units molecular dynamics。主要特点是用GPU加速多体势，但所有的功能也有纯CPU版本。 经常看到有人讨论怎么计算stress-strain曲线。今天我用自己的程序算了一下石墨烯的拉伸曲线，感觉结果对势模型依赖比较大。文献中一般都是用AIREBO势做这个事情，但我的程序现在只有Tersoff势。不改变参数时算出的 fracture strain非常大。改了两个势参数（将Tersoff势中的两个半径参数都增加 0.2 Angstrom）之后结果好一些了。这个改变不影响杨氏模量Y的值。我计算出来的值很接近Y=1TPa。 先贴一下我的程序做这个事情的脚本(初始坐标以及相关的信息要从另一个文件读入）。看上去是不是比lammps简单？ ========================================================================= # This input computes the stress-strain relation for a pristine graphene sample potential potentials/graphene_2.tersoff # file name which contain the potential parameters velocity 300.0 # 300 K # equilibration: ensemble npt_ber 300 0.005 0.0 0.0 0.0 0.0005 # NPT (Berendsen method) time_step 1.0 # 1 fs run 10000 # 10 ps # compute the stress-strain relation: ensemble npt_ber 300 0.005 0.0 0.0 0.0 0.0005 # NPT (Berendsen method) compute_stress_strain 0.0005 1 # strain rate = 0.0005 A/fs; along the zigzag direction dump_thermo 1 # dump results every time step run 25000 # 25 ps ========================================================================= 体系大小为960个原子，模拟总时间为35ps，用我的程序的串行CPU版本耗时38秒。我感觉比lammps的串行版本快两三倍。 再说说我的算法。我要计算x方向（选为zigzag方向）的stress-strain曲线。在x,y方向都用周期边界条件，并采用NPT系综，控制温度在300开附近，y方向的压强为零附近。x方向的盒子长度则以0.0005 A/fs的速率增大（单轴拉伸）。我的体系大概是5nm x 5nm的大小。拉伸25ps相当于将盒子拉长1.25 nm，得到的最大形变为1.25/5=25%。这个时候石墨烯已经被拉断了。结果见下图：  untitled.png |
回复此楼» 收录本帖的淘帖专辑推荐
 GPUMD学习 |
» 猜你喜欢
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有14人回复
-
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了
已经有11人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有11人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有10人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
医学超声影像负责人招聘-中国科学院赣江创新研究院
+1/978
-
西湖大学2026年秋季入学物理学、光学、电子信息方向博士生有名额速来！！！
+2/254
-
985教授征女友
+1/252
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘
+2/244
-
龙凤Tai——写给恋人的第100封情书
+1/85
-
江西理工大学联合中国科学院赣江创新研究院招收2026级博士研究生
+1/82
-
广州大学“长江学者”教授团队2026年海内外高层次人才招聘（环境/化学/生物）
+1/77
-
上海交通大学大气环境科学课题组招收2026年入学博士生
+1/67
-
香港中文大学医学院 诚聘 研究助理教授 (医工结合/生物信息学方向)
+1/34
-
中国农业大学安杰课题组招聘科研助理（表现优异者可提供读博机会）
+1/31
-
同济大学脑机智能团队脑机接口方向招生招聘
+1/27
-
南京医科大学国家级高层次青年人才团队招收博士后
+1/25
-
湘潭大学“过程强化与绿色化工”创新团队补招2026年秋入学博士生
+2/22
-
中国科学技术大学环境系招生
+1/9
-
【东南大学博士后、科研助理招聘】
+1/5
-
澳科大招收2026年秋季药物递送/生物材料方向全奖博士研究生（3月5日18:00截止）
+1/4
-
【科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/3
-
大连海事大学船舶洁净能源研究中心2026年博士研究生招生启事
+1/3
-
美国密苏里大学“柔性电子”课题组诚招博士研究生
+1/2
-
太原理工大学集成电路学院院长团队招收2026年博士研究生
+1/2
2楼2016-08-12 06:05:44
12楼2016-08-12 07:59:09
21楼2016-08-12 08:36:30
30楼2016-08-12 09:12:33
简单回复
2016-08-12 06:27
回复
brucefan(金币+1): 谢谢参与
2016-08-12 06:45
回复
brucefan(金币+1): 谢谢参与
2016-08-12 07:04
回复
brucefan(金币+1): 谢谢参与
hydzp6楼
2016-08-12 07:32
回复
brucefan(金币+1): 谢谢参与
a
tzynew7楼
2016-08-12 07:38
回复
brucefan(金币+1): 谢谢参与
leic8楼
2016-08-12 07:40
回复
brucefan(金币+1): 谢谢参与
天涯20129楼
2016-08-12 07:54
回复
brucefan(金币+1): 谢谢参与
alaqi10楼
2016-08-12 07:56
回复
brucefan(金币+1): 谢谢参与
超级老快11楼
2016-08-12 07:57
回复
brucefan(金币+1): 谢谢参与
miaojiabing13楼
2016-08-12 07:59
回复
brucefan(金币+1): 谢谢参与
laom201314楼
2016-08-12 08:00
回复
brucefan(金币+1): 谢谢参与
jiechen15楼
2016-08-12 08:01
回复
brucefan(金币+1): 谢谢参与
祝福
sunyanwei16楼
2016-08-12 08:11
回复
brucefan(金币+1): 谢谢参与
张屠户17楼
2016-08-12 08:15
回复
zhchzhsh207618楼
2016-08-12 08:24
回复
brucefan(金币+1): 谢谢参与
祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福祝福
bj200898919楼
2016-08-12 08:31
回复
brucefan(金币+1): 谢谢参与
李慧玲20楼
2016-08-12 08:36
回复
brucefan(金币+1): 谢谢参与
纳米材料75322楼
2016-08-12 08:44
回复
brucefan(金币+1): 谢谢参与
youngen23楼
2016-08-12 08:48
回复
brucefan(金币+1): 谢谢参与
落花剑雨24楼
2016-08-12 08:49
回复
brucefan(金币+1): 谢谢参与
灵之红蝶25楼
2016-08-12 08:51
回复
灵之红蝶26楼
2016-08-12 08:51
回复
caledonian27楼
2016-08-12 09:00
回复
brucefan(金币+1): 谢谢参与
star_zhang28楼
2016-08-12 09:09
回复
brucefan(金币+1): 谢谢参与
祝福
springburn29楼
2016-08-12 09:12
回复
brucefan(金币+1): 谢谢参与
tornado13331楼
2016-08-12 09:12
回复
brucefan(金币+1): 谢谢参与
gengxin6032楼
2016-08-12 09:14
回复
brucefan(金币+1): 谢谢参与
jyjydd33楼
2016-08-12 09:18
回复
brucefan(金币+1): 谢谢参与
雪藏剑34楼
2016-08-12 09:19
回复
粲然一笑7935楼
2016-08-12 09:21
回复
brucefan(金币+1): 谢谢参与
343536楼
2016-08-12 09:26
回复
brucefan(金币+1): 谢谢参与
lvluohua37楼
2016-08-12 09:28
回复
brucefan(金币+1): 谢谢参与
悠悠梅梅38楼
2016-08-12 09:29
回复
brucefan(金币+1): 谢谢参与
dctx060839楼
2016-08-12 09:30
回复
brucefan(金币+1): 谢谢参与
sz_kebaoyuan40楼
2016-08-12 09:33
回复
brucefan(金币+1): 谢谢参与
各种爆41楼
2016-08-12 09:37
回复
brucefan(金币+1): 谢谢参与
hol 发自小木虫IOS客户端
yanyuxiaoxi42楼
2016-08-12 09:41
回复
brucefan(金币+1): 谢谢参与
。 发自小木虫Android客户端
wtbccq50543楼
2016-08-12 09:56
回复
brucefan(金币+1): 谢谢参与
朝都44楼
2016-08-12 10:03
回复
brucefan(金币+1): 谢谢参与
好 发自小木虫Android客户端
pingyanwe45楼
2016-08-12 10:14
回复
brucefan(金币+1): 谢谢参与
假大空46楼
2016-08-12 10:18
回复
brucefan(金币+1): 谢谢参与
顶
baiyang31347楼
2016-08-12 10:21
回复
brucefan(金币+1): 谢谢参与
liucheng88348楼
2016-08-12 10:34
回复
brucefan(金币+1): 谢谢参与
yuzhichao49楼
2016-08-12 10:36
回复
brucefan(金币+1): 谢谢参与
祝福
tan91450楼
2016-08-12 10:45
回复
brucefan(金币+1): 谢谢参与