计算化学公社

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

[GROMACS] 求助算十二胺在锂云母层表面的吸附,能量最小化后固体层结构改变、有机分子插入内部

[复制链接 Copy URL]

25

帖子

0

威望

585

eV
积分
610

Level 4 (黑子)

想要研究十二胺(CH₃(CH₂)₁₁NH₂)在锂云母层表面的吸附行为。在模拟中,十二胺采用GAFF力场,锂云母层采用CLAYFF力场。然而,在体系能量最小化后,发现锂云母层结构失稳,且十二胺分子插入了矿物层内部。
恳请各位老师指教,出现这一现象可能的原因有哪些,以及应如何调整优化?同时附上十二胺与锂云母的拓扑文件,烦请各位老师帮忙审阅其中是否存在问题。

能量最小化前后对比.png (483.4 KB, 下载次数 Times of downloads: 0)

能量最小化前后对比.png

1.itp

34.23 KB, 下载次数 Times of downloads: 0

十二胺拓扑文件

2.itp

6.94 KB, 下载次数 Times of downloads: 1

锂云母拓扑文件

849

帖子

4

威望

1915

eV
积分
2844

Level 5 (御坂)

A Student

2#
发表于 Post on 2025-12-2 19:43:50 | 只看该作者 Only view this author
看上去有个分子跨pbc了没处理就直接扩胞,固体夹在同一分子的两个片段中间了。优化时当然会跑回去变完整分子。先确保分子完整地方在盒内。
敬仰一针见血的指责,厌倦别有用心的赞美。

25

帖子

0

威望

585

eV
积分
610

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-12-3 11:09:57 | 只看该作者 Only view this author
student0618 发表于 2025-12-2 19:43
看上去有个分子跨pbc了没处理就直接扩胞,固体夹在同一分子的两个片段中间了。优化时当然会跑回去变完整分 ...

非常感谢您的指点,这确实是问题所在。但是在优化过程中,我发现尽管对固体层添加位置限制,能量最小化之后结构仍然不能很好的维持,会出现某些原子翻转的情况,想请教您有遇到过类似的情况嘛,有没有解决的方向呢?再次感谢您的指导!

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 22:12 , Processed in 0.515346 second(s), 23 queries , Gzip On.

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