计算化学公社

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

[新手求助] Gaussian 计算电子密度文件基组设置方面的疑问

[复制链接 Copy URL]

4

帖子

0

威望

23

eV
积分
27

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 mtlcct 于 2022-5-25 10:07 编辑

,这两个gaussian输入文件在计算chk文件后得到的电子密度文件上有区别的原因是什么,我算他们之间数值上差了好几个数量级?哪个是比较正确的算法呢?

255

帖子

1

威望

2682

eV
积分
2957

Level 5 (御坂)

2#
发表于 Post on 2022-5-25 11:15:53 | 只看该作者 Only view this author
我以为这是什么没用过的软件的输入文件,结果竟是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

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
mtlcct + 2 谢谢分享
snljty2 + 2

查看全部评分 View all ratings

我唯一知道的就是我一无所知,但我是化学小迷弟

4

帖子

0

威望

23

eV
积分
27

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-5-25 13:22:04 | 只看该作者 Only view this author
ChemG 发表于 2022-5-25 11:15
我以为这是什么没用过的软件的输入文件,结果竟是Gaussian。
如果只是想得到电子密度或者电子密度差,输入 ...

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

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

4#
发表于 Post on 2022-5-25 16:29:02 | 只看该作者 Only view this author
mtlcct 发表于 2022-5-25 06:22
谢谢哈,我会仔细看看的。
现在主要就是想问问就是这两个高斯输入文件造成最后电子密度cube文件数值不一 ...

先把不必要的关键词删掉,看看能否重复出这个问题,然后把删过关键词的输入文件贴上来,我们再看。
问”为什么A和B两种做法结果不一样“这类问题,必须先自行把A和B的关键词尽量简化,便于回答者看出原因,这是在网上问问题的常识。你写这么多关键词在上面,就好比给两张各画了100个人的画叫我们找不同,白白浪费回答者的时间。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

127

帖子

0

威望

2769

eV
积分
2896

Level 5 (御坂)

5#
发表于 Post on 2022-5-25 19:04:11 | 只看该作者 Only view this author
mtlcct 发表于 2022-5-25 13:22
谢谢哈,我会仔细看看的。
现在主要就是想问问就是这两个高斯输入文件造成最后电子密度cube文件数值不一 ...

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

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

4

帖子

0

威望

23

eV
积分
27

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2022-5-25 19:31:58 | 只看该作者 Only view this author
抱歉,下次我会改正
,像图中这样理解对吗?所以是这两个文件所用基组不同,而赝势是一样的?

655

帖子

1

威望

5327

eV
积分
6002

Level 6 (一方通行)

7#
发表于 Post on 2022-5-25 19:47:16 | 只看该作者 Only view this author
mtlcct 发表于 2022-5-25 19:31
抱歉,下次我会改正
,像图中这样理解对吗?所以是这两个文件所用基组不同,而赝势是一样的?

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

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
mtlcct + 3 我很赞同

查看全部评分 View all ratings

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

4

帖子

0

威望

23

eV
积分
27

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2022-5-25 19:58:37 | 只看该作者 Only view this author
北大-陶豫 发表于 2022-5-25 19:47
不是的。
你的第一个箭头是指定溴的赝势基组,第二个箭头是指定溴的赝势。
赝势方法是把一部分内层电子 ...

谢谢指导,我去仔细看看博文!

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 00:00 , Processed in 1.473085 second(s), 32 queries , Gzip On.

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