计算化学公社

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

[CASTEP/Dmol3/MS] CASTEP进行几何优化任务时晶胞内分子结构严重变形,请各位大佬指点

[复制链接 Copy URL]

47

帖子

0

威望

758

eV
积分
805

Level 4 (黑子)

本帖最后由 小庄师弟 于 2022-11-30 14:33 编辑

在使用CASTEP进行几何优化时发现晶胞内的分子结构严重变形,请问可能会是什么原因导致的呢?晶胞结构是从CCDC晶体库中下载的,分子式为C7H15Cl2N2O2P。
几何优化前收敛性测试没问题,使用GGA-PBE泛函,精度全部使用Ultra-fine,勾选DFT-D,优化晶胞fixed basis quality,Compressibility选择Soft,使用OTFG模守恒赝势。
优化过程没有完成我就停止了任务,分子结构变形太严重,结果如下图所示。


微信截图_20220607154542.png (89.5 KB, 下载次数 Times of downloads: 52)

优化前

优化前

微信截图_20220607154605.png (98.07 KB, 下载次数 Times of downloads: 38)

优化后

优化后

76

帖子

0

威望

1164

eV
积分
1240

Level 4 (黑子)

2#
发表于 Post on 2022-8-14 10:38:57 | 只看该作者 Only view this author
您好!请问你的问题解决了吗

47

帖子

0

威望

758

eV
积分
805

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-8-17 09:34:51 | 只看该作者 Only view this author
石头 发表于 2022-8-14 10:38
您好!请问你的问题解决了吗

没有,调整参数,换方法都不管用。您也遇到相同的问题了吗?
回复不及时,这几天仅限LV4及以上用户发言。

76

帖子

0

威望

1164

eV
积分
1240

Level 4 (黑子)

4#
发表于 Post on 2022-8-17 10:20:53 | 只看该作者 Only view this author
小庄师弟 发表于 2022-8-17 09:34
没有,调整参数,换方法都不管用。您也遇到相同的问题了吗?
回复不及时,这几天仅限LV4及以上用户发言 ...

是的,我也在优化一个有机分子,连苯环都变形了,但是目前仍未解决,在想是不是castep不擅长这种柔性小分子

47

帖子

0

威望

758

eV
积分
805

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-8-17 10:55:52 | 只看该作者 Only view this author
石头 发表于 2022-8-17 10:20
是的,我也在优化一个有机分子,连苯环都变形了,但是目前仍未解决,在想是不是castep不擅长这种柔性小分 ...

我之前算过一个包含9个C组成的碳链结构,优化完成后看着没问题,C链和初始结构相差不大,我想是不是和加没加H有关,我算的C链结构C上有H或者甲基,初始结构和优化完成的结构都类似锯齿状。

76

帖子

0

威望

1164

eV
积分
1240

Level 4 (黑子)

6#
发表于 Post on 2022-8-17 11:29:26 | 只看该作者 Only view this author
小庄师弟 发表于 2022-8-17 10:55
我之前算过一个包含9个C组成的碳链结构,优化完成后看着没问题,C链和初始结构相差不大,我想是不是和加 ...

难道是因为我没加氢,我试试

47

帖子

0

威望

758

eV
积分
805

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-8-17 11:43:53 | 只看该作者 Only view this author
石头 发表于 2022-8-17 11:29
难道是因为我没加氢,我试试

自己建模可以试试吧。我是在CCDC上下载的晶体文件,然后直接算,加H之后感觉更不对了。

76

帖子

0

威望

1164

eV
积分
1240

Level 4 (黑子)

8#
发表于 Post on 2022-8-17 14:39:08 | 只看该作者 Only view this author
小庄师弟 发表于 2022-8-17 11:43
自己建模可以试试吧。我是在CCDC上下载的晶体文件,然后直接算,加H之后感觉更不对了。

我就是在CCDC上下载的晶体

47

帖子

0

威望

758

eV
积分
805

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-8-17 15:41:55 | 只看该作者 Only view this author
石头 发表于 2022-8-17 14:39
我就是在CCDC上下载的晶体

那可能还是得想想别的办法。我算过体系很大的,也算过链比它长的,还是第一次出这问题。

76

帖子

0

威望

1164

eV
积分
1240

Level 4 (黑子)

10#
发表于 Post on 2022-8-18 09:58:20 | 只看该作者 Only view this author
小庄师弟 发表于 2022-8-17 15:41
那可能还是得想想别的办法。我算过体系很大的,也算过链比它长的,还是第一次出这问题。

你这也是有机小分子,为什么晶体盒子里放了好几个分子,不应该放一个吗?

47

帖子

0

威望

758

eV
积分
805

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2022-8-18 13:16:50 | 只看该作者 Only view this author
石头 发表于 2022-8-18 09:58
你这也是有机小分子,为什么晶体盒子里放了好几个分子,不应该放一个吗?

晶胞中分子个数和空间群有关吧,晶体文件中有相关信息。具体晶胞中分子个数怎么确定我也不清楚,我用的都是前人建好的晶胞。

76

帖子

0

威望

1164

eV
积分
1240

Level 4 (黑子)

12#
发表于 Post on 2022-8-18 16:08:04 | 只看该作者 Only view this author
小庄师弟 发表于 2022-8-18 13:16
晶胞中分子个数和空间群有关吧,晶体文件中有相关信息。具体晶胞中分子个数怎么确定我也不清楚,我用的都 ...

哦哦,我只放了一个分子

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 15:44 , Processed in 0.250516 second(s), 23 queries , Gzip On.

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