计算化学公社

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

[GROMACS] 求助:模拟丙烷在石英壁面接触角时,丙烷无法成液滴状

[复制链接 Copy URL]

17

帖子

0

威望

195

eV
积分
212

Level 3 能力者

各位老师好,我想模拟丙烷在全羟基化石英壁面下的接触角,丙烷采用TraPPE-UA力场,液滴采用常压低温(180K)密度值填充的液球,石英壁面采用clayff力场,top文件利用sobtop生成,整体壁面呈电中性,然后模拟180K甚至更低温度下丙烷的接触角,壁面freeze住了,最小化的时候液滴没有散开,正品模拟的时候,就丙烷液球就开始自由扩散了,这样就没法进行接触角的模拟,尝试过cut-off计算库仑力方式,可以维持液滴状,但不知道可不可以,采用PME计算液滴就会自由扩散,请问有大佬知道该如何解决么,烦请指点一下,感谢。

f8f443b8c03d00a43d77ac7dd9391c7.png (119.59 KB, 下载次数 Times of downloads: 6)

f8f443b8c03d00a43d77ac7dd9391c7.png

C3H8.pdb

334 Bytes, 下载次数 Times of downloads: 0

system.top

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

SiO2.top

389 Bytes, 下载次数 Times of downloads: 5

NVT.mdp

2.69 KB, 下载次数 Times of downloads: 3

em.mdp

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

C3H8.itp

725 Bytes, 下载次数 Times of downloads: 5

325

帖子

0

威望

5170

eV
积分
5495

Level 6 (一方通行)

2#
发表于 Post on 2025-4-15 20:58:25 | 只看该作者 Only view this author
确实挺有意思,itp里面丙烷的电荷都是0,这样按理说静电加和用哪种方式都没有区别呢,纯粹是vdw作用

17

帖子

0

威望

195

eV
积分
212

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-4-16 11:59:03 | 只看该作者 Only view this author
mol 发表于 2025-4-15 20:58
确实挺有意思,itp里面丙烷的电荷都是0,这样按理说静电加和用哪种方式都没有区别呢,纯粹是vdw作用

那么请问大佬,像这样采用cutoff进行计算得到的结果是不是也没问题

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 17:32 , Processed in 0.153504 second(s), 23 queries , Gzip On.

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