| 查看: 4396 | 回复: 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学习 |
» 猜你喜欢
博士读完未来一定会好吗
已经有23人回复
-
导师想让我从独立一作变成了共一第一
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
Ei期刊青年编委招募（工程设计方向）
+1/372
-
哈尔滨工业大学深圳校区—国家级青年人才王龙龙教授团队—诚聘新能源电池方向博士后
+1/181
-
转录组数据求助
+1/83
-
中国医学科学院放射医学研究所、中国医学科学院天津医学健康研究院彭鑫课题组招生招聘
+1/79
-
广州，征女友
+1/77
-
燕山大学亚稳材料全国重点实验室2026年硕士/博士研究生招生信息
+1/67
-
广州医科大学招聘微塑料生物毒理纳米材料方向博士后2名
+1/33
-
校长团队招博士生和博士后
+1/30
-
浙江农林大学森林食物资源挖掘与利用全国重点实验室2026年博士生招生
+1/29
-
澳门科技大学2026年数学博士招生——计算物理与数学课题组: 相场与计算流体动力学
+1/11
-
博士/硕士招生
+1/9
-
招聘2026年入学博士生
+1/7
-
HBI的合成交流
+1/7
-
国家“双一流”建设高校-南京林业大学-国家级青年人才团队招聘
+1/7
-
招若干有机合成人员 (中山大学)
+1/6
-
北京理工大学（珠海）医疗机器人实验室现招收2026年博士研究生
+1/6
-
2026年-中国科学院环境分析化学方向“申请考核”博士生
+1/4
-
深圳先进院三院院士成会明团队诚聘液流电池/高分子合成等方向博后、科研助理、工程师
+1/3
-
海南大学化学院 招聘 材料与电化学方向——研究助理，博士（2026年入学）
+1/1
-
华中农业大学生命科学技术学院 “酶智能设计与农药生物合成”课题组招聘启事
+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): 谢谢参与