计算化学公社

标题: MCM-41簇模型的ESP电荷计算求助 [打印本页]

作者
Author: lisa25496 时间: 2019-5-20 09:39
标题: MCM-41簇模型的ESP电荷计算求助
按照文献，MCM-41是根据SiO2_21A.3d.msi建超胞挖孔，然后加氢饱和等一些列操作完成的，为精确计算整个结构的电荷，文献采用计算簇模型ESP电荷的方法进然后取其平均得到整个结构的ESP电荷。
模拟采用DMOL OPT任务完成，应用GGA PBE 函数，选用DNP基组，原子核采用semicore pseudopotential描述。（均按照文献进行设置）

计算结果如下：
原子       电荷/e       文献电荷/e
Si       1.2187       0.1222
O       -0.6866    -0.06157
H       0.42911    0.0318
几乎与文献电荷差十倍，不知道问题出在什么地方。计算文件太大上传不了，只能上传我的簇模型和输出文件，请各位大佬指点一二！

作者
Author: sobereva 时间: 2019-5-20 10:21
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的标题“ESP 电荷计算”改了，以后务必注意

如置顶的新社员必读贴和论坛首页的公告栏所示，上传超过500KB的文本型文件上传前必须先压缩再上传，以节约论坛空间、加快下载速度。以后务必注意

作者
Author: sobereva 时间: 2019-5-20 10:22
对于当前的簇模型的拟合静电势电荷，文献里的电荷明显不对，数值不可能那么小，要么误解了数值含义

作者
Author: lisa25496 时间: 2019-5-20 10:33

sobereva 发表于 2019-5-20 10:22
对于当前的簇模型的拟合静电势电荷，文献里的电荷明显不对，数值不可能那么小，要么误解了数值含义

谢谢sob老师，根据我的理解，我是将最后得到的ESP电荷中的chg各个原子进行了数值平均最后得到的各个原子的电荷，不知道有没有理解错。文献值应该不可能算错，因为好几篇文献中的电荷值都是这么多。另外计算文件太大，网站上传不了，不知道能不能发给您看看？

作者
Author: sobereva 时间: 2019-5-21 01:40

lisa25496 发表于 2019-5-20 10:33
谢谢sob老师，根据我的理解，我是将最后得到的ESP电荷中的chg各个原子进行了数值平均最后得到的各个原子 ...

我没看你过说的文献，我只能很确切地告诉你，对于你当前的这种模型、对于拟合静电势电荷，你算出来的数量级才是正合理的。诸如对氢来说，每个氢的电荷数值都不小，是否取平均不会影响数量级。

我很反感dmol3，此程序详情我不清楚，给我输出文件也没用
记得以后凡是大文件，一律通过网盘贴链接的方式放在帖子里

更多关于拟合静电势电荷的信息看下文，里面有一些例子
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

作者
Author: lisa25496 时间: 2019-5-21 10:03

sobereva 发表于 2019-5-21 01:40
我没看你过说的文献，我只能很确切地告诉你，对于你当前的这种模型、对于拟合静电势电荷，你算出来的数量 ...

好的谢谢老师指点

作者
Author: RUSTYMI 时间: 2023-5-4 11:01
你好，我最近也在做MCM-41分子筛的相关计算，请问您是否可以把本帖提到的文献共享一下呢？谢谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)