计算化学公社

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

[CASTEP/Dmol3/MS] castep结构优化后出现原子消失

[复制链接 Copy URL]

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

各位老师好,我在用MS castep模块进行结构优化然后能量计算时发现一个很奇怪的现象,结构优化时结构还是正常的,如图一,但把结构优化后的结构再投入能量计算时发现有部分原子消失了(如图2),或者将display style先改成in cell,再改成default,也是这种情况。想请教造成这种现象的可能原因是什么呢?

图1.png (43.71 KB, 下载次数 Times of downloads: 83)

图1

图1

图2.png (49.41 KB, 下载次数 Times of downloads: 108)

图2

图2

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2019-10-21 12:22:01 | 只看该作者 Only view this author
消失正常啊,边界上重复的原子就不显示了啊。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-10-21 21:30:45 | 只看该作者 Only view this author
卡开发发 发表于 2019-10-21 12:22
消失正常啊,边界上重复的原子就不显示了啊。

明白了,谢谢卡卡老师。再请教您一个问题,我现在在很多文献中看到用第一性原理计算过渡态,进而得到原子在固体中的扩散系数的,但没有分子的,请问用过渡态理论可以计算分子在固体中的扩散系数吗?比如氢气分子

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2019-10-21 22:44:07 | 只看该作者 Only view this author
战泽天 发表于 2019-10-21 21:30
明白了,谢谢卡卡老师。再请教您一个问题,我现在在很多文献中看到用第一性原理计算过渡态,进而得到原子 ...

原子和分子应该并没有什么显著差别,原理上可以这样做。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-10-21 22:46:00 | 只看该作者 Only view this author
卡开发发 发表于 2019-10-21 22:44
原子和分子应该并没有什么显著差别,原理上可以这样做。

好的,谢谢卡卡老师的指导

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2019-10-21 23:01:24 | 只看该作者 Only view this author
战泽天 发表于 2019-10-21 22:46
好的,谢谢卡卡老师的指导

要做扩散的话格子可能还得再大点。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-10-26 08:50:41 | 只看该作者 Only view this author
卡开发发 发表于 2019-10-21 23:01
要做扩散的话格子可能还得再大点。

抱歉,卡卡老师,这几天一直没有上网站,没看您的消息。这个格子大是根据什么判断的呢?我一直都是用之前师兄做吸附的结构,没有考虑过这个问题

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

8#
发表于 Post on 2019-10-26 11:22:51 | 只看该作者 Only view this author
战泽天 发表于 2019-10-26 08:50
抱歉,卡卡老师,这几天一直没有上网站,没看您的消息。这个格子大是根据什么判断的呢?我一直都是用之前 ...

尽量避免扩散的粒子间有相互作用即可。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2019-10-26 11:35:19 | 只看该作者 Only view this author
卡开发发 发表于 2019-10-26 11:22
尽量避免扩散的粒子间有相互作用即可。

明白了。再请教您一个关于结构优化的问题,我在fine的精度下进行结构优化,虽然并没有收敛,但得到一个相对优化的结构,我用这个结构调整参数重新优化,跟用之前完全没有优化的结构相比,可以节省优化时间吗?

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

10#
发表于 Post on 2019-10-26 14:08:41 | 只看该作者 Only view this author
战泽天 发表于 2019-10-26 11:35
明白了。再请教您一个关于结构优化的问题,我在fine的精度下进行结构优化,虽然并没有收敛,但得到一个相 ...

看调整哪些参数,有些可以节省,有些反而不行。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2019-10-26 14:27:30 | 只看该作者 Only view this author
卡开发发 发表于 2019-10-26 14:08
看调整哪些参数,有些可以节省,有些反而不行。

想要把cut-off energy提高一些,从400eV调到500eV,然后把SCF tolerance 降低,从1.0e-6调到5.0e-6,调整这两个参数的话可以节省时间吗?

3809

帖子

3

威望

1万

eV
积分
20343

Level 6 (一方通行)

围观吃瓜群众

12#
发表于 Post on 2019-10-26 16:09:30 | 只看该作者 Only view this author
战泽天 发表于 2019-10-26 14:27
想要把cut-off energy提高一些,从400eV调到500eV,然后把SCF tolerance 降低,从1.0e-6调到5.0e-6,调整 ...

可以这么做,但效果怎样我不好说。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

43

帖子

0

威望

613

eV
积分
656

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2019-10-26 20:07:06 | 只看该作者 Only view this author
卡开发发 发表于 2019-10-26 16:09
可以这么做,但效果怎样我不好说。

好的

105

帖子

0

威望

1133

eV
积分
1238

Level 4 (黑子)

14#
发表于 Post on 2021-6-8 20:34:32 | 只看该作者 Only view this author
有没有能显示的方法啊?

17

帖子

0

威望

101

eV
积分
118

Level 2 能力者

15#
发表于 Post on 2021-10-15 13:23:22 | 只看该作者 Only view this author
本帖最后由 寻找Medusa 于 2021-10-15 13:40 编辑
卡开发发 发表于 2019-10-21 12:22
消失正常啊,边界上重复的原子就不显示了啊。

卡卡老师,做吸附的时候出现这种情况还能再继续进行计算吗

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 18:21 , Processed in 0.177054 second(s), 23 queries , Gzip On.

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