24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 559  |  回复: 7
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

素婕求学

铁虫 (小有名气)

[交流] 【求助】有自己编程算CNT或Graphene的嘛?急于求助中。。。。。。

本人自己编程算的graphene, 用的是second generation REBO potential,在NVT系综下模拟,没有加任何边界条件。现在用不同的时间步长来算,发现时间步长大一点的,系统算到1000步就崩溃了;时间步长小一点的,系统温度能达到设定温度,从能量变化趋势上来看,总能量呈下降趋势,动能没什么变化(这应该和加了温度调控有关)。我想问下,这样的情况属于正常吗?有没有人算过相类似的东西呢?你在NVT系综下算时,系统呈现出什么状态呢?我现在都不知道该如何调这个程序呢,茫然中~~~~~~~~~
望高手指点,万分感谢。。。。。。。。。。。。。。。。。

另外,我在想一个问题,对于像carbon nanotube or graphene等有具体尺寸和结构本身就属于纳米尺度的材料,在用分子动力学模拟时,同样需要建立模拟盒,需要使用一定的边界条件,如周期性边界条件嘛? 如果需要使用周期性边界条件,该如何使用呢?好像不是简单的按照任一原子离开元胞,就有镜像的原子进来的原则吧?

希望高手赐教,这两天要给老板交结果,心急如焚啦~~~~~~~~~

[ Last edited by 素婕求学 on 2008-7-29 at 09:25 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

素婕求学

铁虫 (小有名气)

引用回帖:
Originally posted by acridine at 2008-7-29 11:12:
周期性边界条件不会破坏键的,比如简单的一维情况:
|-a-b-c-d-| (注意这里d-a是成键的).
虽然如果d原子从右侧出去会从左侧回来,变成
|-d-a-b-c-| (注意这里c-d是成键的),
但是d-a以及c-d的成键关系是没变的

嗯,我明白你的意思

可是按照周期性条件来处理,原子d 的坐标会被替换成xr=xr-box*nint(xr/box)
那本来原子在右侧,和它相邻的原子也在右侧,现在原子d的坐标被替换到左边去了,那计算距离时,它就不满足成键的距离了呀
是不是我的公式弄错了呀
哎,着急呀~
6楼2008-07-29 13:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

素婕求学

铁虫 (小有名气)

高手们都起床了嘛?
自己来顶
2楼2008-07-29 09:26:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
时间步长过大确实容易导致结构变乱,尤其是温度较高或远离平衡结构的情况,具体时间步长的量级和势函数也有关系,可以参考使用同一势函数的文献
做NVT体系,最初的驰豫阶段总能量确实应是逐渐降低的,动能不变,而势能在逐渐降低,但驰豫一段时间之后总能应该就趋于恒定,因为已经找到势能最低的结构,如果是这样的话程序就应该没太大问题了.
周期边界条件就是你说的那样吧,计算距离时注意与镜像原子的距离就可以了
3楼2008-07-29 09:49:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

素婕求学

铁虫 (小有名气)

引用回帖:
Originally posted by acridine at 2008-7-29 09:49:
时间步长过大确实容易导致结构变乱,尤其是温度较高或远离平衡结构的情况,具体时间步长的量级和势函数也有关系,可以参考使用同一势函数的文献
做NVT体系,最初的驰豫阶段总能量确实应是逐渐降低的,动能不变,而势能 ...

非常谢谢!!!!!!!!!!

不过,我还是觉得周期性边界条件在碳管或者graphene中的应用应该有特殊处理的方法,因为如果应用镜像原则来处理,从元胞左边出去的碳原子,会从右边补充进来,那碳原子间的键作用关系就会打破,六边形形状的结构就不能保持了。这个到底该怎么处理呢?
跪请大侠们指点。。。。。。。。。。。。。。。。。。
4楼2008-07-29 10:26:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
[考研] 085600材料与化工专硕329 求调剂 +18 额cc 2026-04-06 19/950 2026-04-07 19:38 by qlm5820
[考研] 生物工程求调剂 +13 喜欢还是不甘心 2026-04-05 13/650 2026-04-07 16:55 by Ecowxq666!
[考研] 265求调剂 +17 小木虫085600 2026-04-06 19/950 2026-04-07 12:04 by Sammy2
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 生物与医药086000调剂一志愿西北农林320分 +3 美美女士 2026-04-03 3/150 2026-04-05 21:55 by 学员8dgXkO
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 【求调剂】新能源材料本科,一志愿211,初试321 +6 求调剂学校, 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
信息提示
请填处理意见