计算化学公社

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

[Material Studio] 求助,MS使用Forcite模块对COF材料进行结构优化,不知道为什么会发生扭曲变形

[复制链接 Copy URL]

6

帖子

0

威望

29

eV
积分
35

Level 2 能力者

各位大神,我使用了MS建立了一种COF材料的晶胞,然后使用 Frocite 模块进行结构的优化,发现得不到论文中一样的2D平面结构,材料中的分子会发生扭曲变形。我看了一下之前相关的贴子,上面讨论说是跟力场有关,但我尝试了模块里面的各种力场,发现结果都差不多。我又尝试了一下精度更高的Dmol3 模块,但都是算了十几分钟后就提示SCF不收敛,不知道有什么比较好的优化方法呢?或者说是MS本身不太适合做该类材料的结构优化?


Screen Shot 2021-06-25 at 7.39.09 PM.png (229.2 KB, 下载次数 Times of downloads: 68)

优化之后的结构

优化之后的结构

Screen Shot 2021-06-25 at 7.38.42 PM.png (231.74 KB, 下载次数 Times of downloads: 65)

优化之后的结构

优化之后的结构

Screen Shot 2021-06-25 at 7.38.20 PM.png (260.96 KB, 下载次数 Times of downloads: 58)

优化前的结构

优化前的结构

39

帖子

0

威望

1984

eV
积分
2023

Level 5 (御坂)

2#
发表于 Post on 2021-6-26 16:27:34 | 只看该作者 Only view this author
属于建模时操作的问题。画好后优化的第一步不要优化晶胞,然后找出对称性后再优化晶胞

6

帖子

0

威望

29

eV
积分
35

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-6-26 19:44:49 | 只看该作者 Only view this author
swordshine 发表于 2021-6-26 16:27
属于建模时操作的问题。画好后优化的第一步不要优化晶胞,然后找出对称性后再优化晶胞

第一步没有优化晶胞呢,都是先直接优化结构而已,然后优化完想找对称性,要么就依旧是p1,要么有其他的对称性,但跟论文中的不一样。我尝试着把论文中的晶胞拿过来,复制了他的单体进来优化,发现虽然扭曲的没有那么厉害了,但还是有一些部分发生了偏转,得不到完全的平面。

6

帖子

0

威望

29

eV
积分
35

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2021-6-26 19:46:37 | 只看该作者 Only view this author

541

帖子

3

威望

6620

eV
积分
7221

Level 6 (一方通行)

5#
发表于 Post on 2021-6-26 20:05:06 | 只看该作者 Only view this author
H呢?搞笑呢?
结构画对了 成键搞对了  你用什么力场都可以。UFF都能搞。
弯折太严重 先别搞优化 先把晶胞参数调整到差不多 再去clean。看着顺眼后 再去优化。
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

39

帖子

0

威望

1984

eV
积分
2023

Level 5 (御坂)

6#
发表于 Post on 2021-6-26 21:24:33 | 只看该作者 Only view this author
本帖最后由 swordshine 于 2021-6-26 21:25 编辑

COF建模第一步可以先不画氢,加了氢原子会非常复杂,不利于找到体系里对称性不对的地方,搞定对称性后再补上氢就可以了。你可以看看这个分子的原文献。如果有一些有机化学的知识,可以知道砜的氧是平面外的,如果是常规的eclipsed堆积,那么会有很大的位阻,要避免这种位阻,显然需要做处理,即再补一层该结构,把中间的含砜的环翻个身,让两层中砜的氧在不同侧。搞定后再补上氢进一步优化。

541

帖子

3

威望

6620

eV
积分
7221

Level 6 (一方通行)

7#
发表于 Post on 2021-6-26 22:20:05 | 只看该作者 Only view this author
不同意楼上老哥关于不带H的操作。我认为 带H clean 和力场下优化才不会出错,如果是为了对称性考虑,那也应该搞完优化后 再去删H 找对称才合适。
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

6

帖子

0

威望

29

eV
积分
35

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2021-6-26 22:43:13 | 只看该作者 Only view this author
ABetaCarw 发表于 2021-6-26 20:05
H呢?搞笑呢?
结构画对了 成键搞对了  你用什么力场都可以。UFF都能搞。
弯折太严重 先别搞优化 先把晶 ...

氢都是先不加的吧,我感觉已经调整得蛮好的了,才去进行结构优化的,但还是找不到论文一样的对称性

6

帖子

0

威望

29

eV
积分
35

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2021-6-26 22:53:44 | 只看该作者 Only view this author
swordshine 发表于 2021-6-26 21:24
COF建模第一步可以先不画氢,加了氢原子会非常复杂,不利于找到体系里对称性不对的地方,搞定对称性后再补 ...

嗯嗯,论文中是提供了这个COF的晶胞文件的,他里面的确实是跟你说的一样,是加了一层的。我按照你说的,也多加了一层呢,但是一优化就是心碎的声音,不知道怎么搞了。我尝试过从他的晶胞文件中复制一层结构放到一个新的空晶胞里面进行结构优化,然后就是得不到原来的P6/M对称性了。传了一个文章的cif文件一起,有兴趣的大神也可以试一下。

Screen Shot 2021-06-26 at 10.36.07 PM.png (231.41 KB, 下载次数 Times of downloads: 67)

优化前

优化前

Screen Shot 2021-06-26 at 10.41.46 PM.png (357.91 KB, 下载次数 Times of downloads: 60)

优化后

优化后

FS-COF_AA_publ.cif

2.7 KB, 下载次数 Times of downloads: 4

论文晶胞数据

541

帖子

3

威望

6620

eV
积分
7221

Level 6 (一方通行)

10#
发表于 Post on 2021-6-26 23:15:05 | 只看该作者 Only view this author
量化计算新手村 发表于 2021-6-26 22:43
氢都是先不加的吧,我感觉已经调整得蛮好的了,才去进行结构优化的,但还是找不到论文一样的对称性

不加H  你进行力场优化或者Dmol3计算 纯属是在浪费计算资源
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

39

帖子

0

威望

1984

eV
积分
2023

Level 5 (御坂)

11#
发表于 Post on 2021-6-27 08:18:54 | 只看该作者 Only view this author
ABetaCarw 发表于 2021-6-26 22:20
不同意楼上老哥关于不带H的操作。我认为 带H clean 和力场下优化才不会出错,如果是为了对称性考虑,那也应 ...

不知你有没有实际做过COF的建模,你可以试试一开始就带上H建个模,简单结构还好,复杂的根本找不到哪里出了问题。至于建模时的初步模拟,一个forcite下的geometry optimization只要几秒钟,何来资源一说

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
ABetaCarw + 1 怼之前先审题

查看全部评分 View all ratings

513

帖子

0

威望

4394

eV
积分
4907

Level 6 (一方通行)

truffle

12#
发表于 Post on 2021-6-27 15:37:22 | 只看该作者 Only view this author
别的程序不知道,对于MS来说如果不加氢优化COF,必须手动分配力场和电荷才可能合理。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
ABetaCarw + 5 我很赞同

查看全部评分 View all ratings

No problem is insoluble in all conceivable circumstances.

192

帖子

0

威望

475

eV
积分
667

Level 4 (黑子)

13#
发表于 Post on 2021-7-30 17:44:14 | 只看该作者 Only view this author
ABetaCarw 发表于 2021-6-26 22:20
不同意楼上老哥关于不带H的操作。我认为 带H clean 和力场下优化才不会出错,如果是为了对称性考虑,那也应 ...

您好,想请教一下,我有个COF,find symmetry的时候给出来的是P1,感觉不对,我预计是P3,请问这是为什么呢,我的结构看起来没问题了

5

帖子

0

威望

93

eV
积分
98

Level 2 能力者

14#
发表于 Post on 2022-11-17 02:54:21 | 只看该作者 Only view this author
以上的回复都没答到key point上,这种情况不管加不加H,都可能变形啊。因为没有对称性限制,一开始至少要用P3空间群限制就可以了。

39

帖子

0

威望

3203

eV
积分
3242

Level 5 (御坂)

15#
发表于 Post on 2022-11-17 03:25:03 | 只看该作者 Only view this author
做了N多次结构模型的表示问题就在于没加H

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 13:08 , Processed in 0.184067 second(s), 27 queries , Gzip On.

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