计算化学公社

标题: Gaussian 计算电子密度文件基组设置方面的疑问 [打印本页]

作者
Author: mtlcct 时间: 2022-5-25 10:07
标题: Gaussian 计算电子密度文件基组设置方面的疑问
本帖最后由 mtlcct 于 2022-5-25 10:07 编辑

(, 下载次数 Times of downloads: 19) 和 (, 下载次数 Times of downloads: 23) ,这两个gaussian输入文件在计算chk文件后得到的电子密度文件上有区别的原因是什么，我算他们之间数值上差了好几个数量级？哪个是比较正确的算法呢？

作者
Author: ChemG 时间: 2022-5-25 11:15
我以为这是什么没用过的软件的输入文件，结果竟是Gaussian。
如果只是想得到电子密度或者电子密度差，输入文件不必这么复杂吧，体系不大还不如直接B3LYP/def2TZVP，Lanl2DZ太糙了
建议看看社长的相关博文
谈谈赝势基组的选用 http://sobereva.com/373
简谈量子化学计算中DFT泛函的选择 http://sobereva.com/272
谈谈量子化学中基组的选择 http://sobereva.com/336
计算化学中的一些常见不良写法和用词 http://sobereva.com/298
使用Multiwfn作电子密度差图 http://sobereva.com/113

作者
Author: mtlcct 时间: 2022-5-25 13:22

ChemG 发表于 2022-5-25 11:15
我以为这是什么没用过的软件的输入文件，结果竟是Gaussian。
如果只是想得到电子密度或者电子密度差，输入 ...

谢谢哈，我会仔细看看的。
现在主要就是想问问就是这两个高斯输入文件造成最后电子密度cube文件数值不一样的原因是什么啊？

作者
Author: wzkchem5 时间: 2022-5-25 16:29

mtlcct 发表于 2022-5-25 06:22
谢谢哈，我会仔细看看的。
现在主要就是想问问就是这两个高斯输入文件造成最后电子密度cube文件数值不一 ...

先把不必要的关键词删掉，看看能否重复出这个问题，然后把删过关键词的输入文件贴上来，我们再看。
问”为什么A和B两种做法结果不一样“这类问题，必须先自行把A和B的关键词尽量简化，便于回答者看出原因，这是在网上问问题的常识。你写这么多关键词在上面，就好比给两张各画了100个人的画叫我们找不同，白白浪费回答者的时间。

作者
Author: fatpig 时间: 2022-5-25 19:04

mtlcct 发表于 2022-5-25 13:22
谢谢哈，我会仔细看看的。
现在主要就是想问问就是这两个高斯输入文件造成最后电子密度cube文件数值不一 ...

为什么？这还用问为什么？你两次用的基组都不一样，你算啥都不一样啊...

你第一张图手动输入的 Br 的基组是 lanl2dzdp 而不是 lanl2dz ...

作者
Author: mtlcct 时间: 2022-5-25 19:31
抱歉，下次我会改正
(, 下载次数 Times of downloads: 28) ，像图中这样理解对吗？所以是这两个文件所用基组不同，而赝势是一样的？

作者
Author: 北大-陶豫 时间: 2022-5-25 19:47

mtlcct 发表于 2022-5-25 19:31
抱歉，下次我会改正
，像图中这样理解对吗？所以是这两个文件所用基组不同，而赝势是一样的？

不是的。
你的第一个箭头是指定溴的赝势基组，第二个箭头是指定溴的赝势。
赝势方法是把一部分内层电子赝掉，所以要定义赝势（如何描述被赝掉的内层电子）和赝势基组（如何描述没有被赝掉的价电子）两部分。仔细看上面提到的卢老师的博文。

作者
Author: mtlcct 时间: 2022-5-25 19:58

北大-陶豫发表于 2022-5-25 19:47
不是的。
你的第一个箭头是指定溴的赝势基组，第二个箭头是指定溴的赝势。
赝势方法是把一部分内层电子 ...

谢谢指导，我去仔细看看博文！

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