| 查看: 4596 | 回复: 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学习 |
» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
DNA甲基化位点定量试剂盒（qPCR版）-适合特定基因位点5mC定量检测
+1/83
-
锌离子混合电容器
+1/81
-
时间的眼神
+1/79
-
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+1/77
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
+1/68
-
真诚找对象
+1/58
-
江西理工大学 稀土学院(发光材料与器件研究所) 招收2026届 材料类博士研究生 2名
+2/52
-
考核制博士自荐
+1/36
-
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1-2名
+1/28
-
SCI，计算机相关可以写
+1/19
-
SCI，计算机相关可以写
+1/17
-
西北工业大学2026年博士研究生“翱翔快响”招生
+1/16
-
SCI，计算机相关可以写
+1/14
-
2026博士招生-上海大学先进耐火材料全国重点实验室-招收冶金工程博士研究生-1-2名
+1/11
-
英国斯旺西大学(Swansea University)招半导体材料和器件方向CSC博士生
+1/5
-
2026 博士自荐-机器人机构学方向
+1/4
-
液晶拓扑光子学博士招生（电子科技大学）
+1/3
-
[招募] 上海交通大学环境健康课题组科研实习生（环境健康与生物学大数据方向）
+1/1
-
中科院（广州）研究团队招聘无机非金属材料方向特别研究助理或博后
+1/1
-
福建师范大学海峡柔性电子学院-2026年9月入学博士-考核制
+1/1
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): 谢谢参与