计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2910|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[MOPAC] mopac半经验优化石墨周围氢原子巨慢

[复制链接 Copy URL]

82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本人想优化一下石墨周围的氢原子来将其作为基底进行下一步计算,出现了如下两种情况:
一、本人更改了双层石墨晶胞的向量(1 1 0 / -1 1 0 / 0 0 1)使其变成矩形晶胞,然后取消晶胞后在边缘碳原子上添加了氢原子,固定所有的C原子利用PM7优化了H,共计1054个原子,28核只用了23min就优化好了
二、没有更改石墨晶胞向量,以其菱形晶胞进行相同的优化,共计1002个原子,同样是固定C优化H,30核跑了4个多小时都跑不完一个循环,请问这是什么原因呀?
图1、图2为我所优化的两种石墨
图3、图4为这两种的参数设置(完全一样)


82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2022-11-14 11:56:49 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-13 15:25
第二个结构的C-H键太长,被MS判断成是不成键的,所以对于石墨烯共振式的分析结果也就不可信了。
这种问 ...

嗯嗯明白了,谢谢

82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-11-14 11:54:45 | 只看该作者 Only view this author
sobereva 发表于 2022-11-13 21:01
只要氢的位置加得不离谱,石墨烯片段尺寸足够大的情况,边缘的氢是否优化都无所谓,不会对中心区域产生显著 ...

好的好的,谢谢sob老师

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

6#
发表于 Post on 2022-11-13 21:01:53 | 只看该作者 Only view this author
只要氢的位置加得不离谱,石墨烯片段尺寸足够大的情况,边缘的氢是否优化都无所谓,不会对中心区域产生显著影响
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

5#
发表于 Post on 2022-11-13 15:25:08 | 只看该作者 Only view this author
Rosefinch 发表于 2022-11-13 05:51
您好,我按照您说的检查了一下我的石墨片段,但是发现对于第一种可以很快优化好的石墨结构,利用MS重新生 ...

第二个结构的C-H键太长,被MS判断成是不成键的,所以对于石墨烯共振式的分析结果也就不可信了。
这种问题建议还是手画,确保只有单键和双键,没有1.5键。乍看起来工作量巨大,但是其实找到规律以后很快就能判断出来能不能画出没有自由基的共振式,毕竟石墨烯的结构是极其有规律的
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-11-13 12:51:32 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-13 02:02
石墨烯按某些方式截断时得到的是闭壳层片段,可以画出单双键交替的共振式;按另一些方式截断时得到开壳层片 ...

您好,我按照您说的检查了一下我的石墨片段,但是发现对于第一种可以很快优化好的石墨结构,利用MS重新生成键会出现几个带着自由基的碳原子;而对于第二个无法优化的却是完整的闭壳层结构,请问这是怎么回事呀?我附上了我的两种石墨结构文件。 type1.pdb (86.55 KB, 下载次数 Times of downloads: 1) type2.pdb (82.29 KB, 下载次数 Times of downloads: 1)



82

帖子

0

威望

1197

eV
积分
1279

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-11-13 11:45:25 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-13 02:02
石墨烯按某些方式截断时得到的是闭壳层片段,可以画出单双键交替的共振式;按另一些方式截断时得到开壳层片 ...

感谢感谢,受教了

1万

帖子

0

威望

9886

eV
积分
22130

Level 6 (一方通行)

2#
发表于 Post on 2022-11-13 02:02:06 | 只看该作者 Only view this author
本帖最后由 wzkchem5 于 2022-11-12 19:03 编辑

石墨烯按某些方式截断时得到的是闭壳层片段,可以画出单双键交替的共振式;按另一些方式截断时得到开壳层片段,不管怎么画共振式都会画出很多自由基。后者SCF收敛会困难得多,而且电子结构比较离域,不利于利用MOZYME等加速算法,导致每步SCF迭代的时间也可能变慢。
如无特殊原因,应该尽量截断成闭壳层片段,因为完整的石墨烯是闭壳层的,或者更严格地说即使有自由基也全在边缘上,对于不涉及石墨烯边缘的研究而言可以看作是闭壳层的

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
413 + 3 不明觉厉

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 00:04 , Processed in 0.681387 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list