计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › x2top生成拓扑文件和cm5电荷计算疑问

x2top生成拓扑文件和cm5电荷计算疑问

查看数: 10980 | 评论数: 6 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

鱼叉君

发布时间: 2019-2-26 13:59

正文摘要:

大家好，我想要模拟图中周期性材料和水的混合，设想的情况是这层材料在盒子中间，上下分布水分子，run一段时间后观察水的分布。目前是有这个材料的pdb文件，我按照Sob老师生成石墨烯拓扑结构的教程修改oplsaa力场， ...

回复 Reply

Aristotler 发表于 Post on 2020-11-2 09:34:51

sobereva 发表于 2020-10-31 09:30
用RESP电荷

另外，不为整数也不意味着就不合理，例如0.8*OPLS-2009IL离子液体力场将离子的原子电荷都 ...

好的，谢谢老师的指导

sobereva 发表于 Post on 2020-10-31 09:30:43

Aristotler 发表于 2020-10-29 13:24
老师，请教一个问题：比如我想算SO42- 用OPLS-AA力场，原子电荷是建议用CM5*1.2, 但是我本身SO42- 有两个 ...

用RESP电荷

另外，不为整数也不意味着就不合理，例如0.8*OPLS-2009IL离子液体力场将离子的原子电荷都乘了0.8，导致离子的净电荷绝对值也为0.8。实际上此力场里直接就有和硫酸根很类似的阴离子，自行稍微改一下弄成硫酸根也可以。

Aristotler 发表于 Post on 2020-10-29 13:24:56

sobereva 发表于 2019-3-8 11:56
不能太烂（不能低于B3LYP/6-31G*），但也绝对不需要太高（最高用到B3LYP/def2-TZVP就足矣），具体适合什 ...

老师，请教一个问题：比如我想算SO42- 用OPLS-AA力场，原子电荷是建议用CM5*1.2, 但是我本身SO42- 有两个负电荷，乘以1.2后就是-2.4e了，那怎么办？

sobereva 发表于 Post on 2019-3-8 11:56:12

苏玖染发表于 2019-3-8 09:23
请问下sob老师，计算cm5电荷一般用什么级别的方法和基组比较准确？

不能太烂（不能低于B3LYP/6-31G*），但也绝对不需要太高（最高用到B3LYP/def2-TZVP就足矣），具体适合什么级别看你的体系大小和计算能力

苏玖染 发表于 Post on 2019-3-8 09:23:15

sobereva 发表于 2019-2-26 16:20
1 最终到底用了什么原子电荷，是拓扑文件里[atoms]字段里体现的

2 输入文件显然不能用pdb，根本都不含有 ...

请问下sob老师，计算cm5电荷一般用什么级别的方法和基组比较准确？

sobereva 发表于 Post on 2019-2-26 16:20:21

1 最终到底用了什么原子电荷，是拓扑文件里[atoms]字段里体现的

2 输入文件显然不能用pdb，根本都不含有波函数信息，依赖于电子密度的CM5电荷显然不可能算出来。得用起码包含GTF信息的输入文件，诸如fch、wfn等，仔细看：
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）

3 是

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register