计算化学公社

 找回密码 Forget password
 注册 Register

碳纳米管储氢的模拟问题

查看数: 9154 | 评论数: 12 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2020-6-1 17:49

正文摘要:

各位老师,我在做一个碳管储氢的课题,体系为1*1*12的CNT(32,0)+1000个H2分子,固定碳管,使用gromacs,opls力场。为了做出储氢的等温(300K)线,我从0.1Mpa开始,通过改变盒子大小来模拟不同的压强,0.1Mpa盒子 ...

回复 Reply

夜夜衣 发表于 Post on 2020-6-9 22:18:51
sobereva 发表于 Post on 2020-6-9 15:37:20
夜夜衣 发表于 2020-6-9 08:10
老师,再请教您下,我这个体系设置vdw-type  = Cut-off 可以吗。

夜夜衣 发表于 Post on 2020-6-9 08:10:11
sobereva 发表于 2020-6-6 10:46
没必要加色散校正。这只是定量改进色散作用能和色散作用对压力的贡献,但只是定量的,不合适的时候加了这 ...

老师,再请教您下,我这个体系设置vdw-type  = Cut-off 可以吗。
夜夜衣 发表于 Post on 2020-6-6 16:24:12
sobereva 发表于 2020-6-6 10:46
没必要加色散校正。这只是定量改进色散作用能和色散作用对压力的贡献,但只是定量的,不合适的时候加了这 ...

谢谢老师
sobereva 发表于 Post on 2020-6-6 10:46:01
夜夜衣 发表于 2020-6-5 11:45
可是老师,如果在较大压强下加不加色散呢?例如5Mpa 10MPa。这样分子就很密集,怎么界定这个事情呢,另外 ...

没必要加色散校正。这只是定量改进色散作用能和色散作用对压力的贡献,但只是定量的,不合适的时候加了这个甚至定性都是错的,审稿人也会批评。
另外,只要范德华作用的cutoff设得稍微大一点,色散校正的影响都会迅速降低。当前你的rvdw=1.2就已经不算小了,被忽视的色散作用对结果的影响微乎其微,绝对不会因此导致动力学行为有什么不合理性。
夜夜衣 发表于 Post on 2020-6-5 11:45:31
本帖最后由 夜夜衣 于 2020-6-6 00:33 编辑
sobereva 发表于 2020-6-3 20:55
用压浴肯定不行,到最后盒子会坍塌得非常小
色散校正绝对不能用,因为你的体系是非均匀的,有大量真空区 ...

可是老师,如果在较大压强下加不加色散呢?例如5Mpa 10MPa。这样分子就很密集,怎么界定这个事情呢,另外 如果模拟氨气这种带氢键的体系,需要加色散矫正吗?
夜夜衣 发表于 Post on 2020-6-5 11:12:04
sobereva 发表于 2020-6-3 20:55
用压浴肯定不行,到最后盒子会坍塌得非常小
色散校正绝对不能用,因为你的体系是非均匀的,有大量真空区 ...

谢谢老师
sobereva 发表于 Post on 2020-6-3 20:55:34
夜夜衣 发表于 2020-6-2 11:47
老师,您觉得我这个体系,用npt系综,同时加上热浴和压浴如何?加上色散矫正是否更合理一点,谢谢您!若 ...

用压浴肯定不行,到最后盒子会坍塌得非常小
色散校正绝对不能用,因为你的体系是非均匀的,有大量真空区,完全违背了色散校正的假设
夜夜衣 发表于 Post on 2020-6-2 21:11:32
sobereva 发表于 2020-6-2 08:33
具体那些文件暂时没时间看,至于energygrp_excl用了报错的那个问题,显然不能写HTR HTR,要不然氢分子之间 ...

老师,问题已经解决 是色散矫正的问题
夜夜衣 发表于 Post on 2020-6-2 11:47:06
本帖最后由 夜夜衣 于 2020-6-2 11:50 编辑
sobereva 发表于 2020-6-2 08:33
具体那些文件暂时没时间看,至于energygrp_excl用了报错的那个问题,显然不能写HTR HTR,要不然氢分子之间 ...

老师,您觉得我这个体系,用npt系综,同时加上热浴和压浴如何?加上色散矫正是否更合理一点,谢谢您!若用npt,那compressibility 该如何设置呢?
夜夜衣 发表于 Post on 2020-6-2 11:35:36
sobereva 发表于 2020-6-2 08:33
具体那些文件暂时没时间看,至于energygrp_excl用了报错的那个问题,显然不能写HTR HTR,要不然氢分子之间 ...

谢谢老师,先前我已经读过手册,这个问题已经解决了,
sobereva 发表于 Post on 2020-6-2 08:33:05
具体那些文件暂时没时间看,至于energygrp_excl用了报错的那个问题,显然不能写HTR HTR,要不然氢分子之间、氢分子内的氢原子间都没有任何相互作用了,不崩溃才怪。能写CNT CNT是因为CNT已经freeze住了,因此不描述CNT的原子间的相互作用也没关系。

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

GMT+8, 2026-2-21 11:31 , Processed in 0.864773 second(s), 25 queries , Gzip On.

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