计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4471#
 楼主 Author| 发表于 Post on 2017-9-20 11:27:55 | 只看该作者 Only view this author
erya123 发表于 2017-9-20 09:49
gaussian09D版本有没有内置的NBO?如何查询NBO版本?

都内置了,内置的是最后一个免费版本3.1
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

13

帖子

0

威望

55

eV
积分
68

Level 2 能力者

4472#
发表于 Post on 2017-9-20 12:05:37 | 只看该作者 Only view this author

计算Gd配合物homo-lumo报错

%chk=Gd.chk
#P opt b3lyp/genecp

Gd

0 1
Gd                -3.78840000    8.05240000    6.44690000
O                 -4.18510000    8.55090000    3.89690000
H                 -3.71940000    3.19460000    2.72990000

Gd 0
MWB46
****
C H N O 0
6-31G*
****

Gd 0
MWB46


本人新手,需要计算Gd配合物的能级,我的输入文件如上,用高斯计算一直报错提示2070,请大家帮忙看看哪里有问题











6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4473#
 楼主 Author| 发表于 Post on 2017-9-20 12:08:04 | 只看该作者 Only view this author
万山红遍 发表于 2017-9-20 10:09
确实是Fe原子中4个类似于花瓣形的d轨道,它们应该是简并的,但是轨道能量看起来是不同的,不知道这跟我计 ...

你算个Mn,就会发现所有d轨道能量都相同。
Fe的5个d轨道本身并不简并,毕竟占据方式其中一个与另外的不同,感受到的势也不同,会造成能量分裂
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4474#
 楼主 Author| 发表于 Post on 2017-9-20 12:09:30 | 只看该作者 Only view this author
sobereva 发表于 2017-9-20 12:08
你算个Mn,就会发现所有d轨道能量都相同。
Fe的5个d轨道本身并不简并,毕竟占据方式其中一个与另外的不 ...

计算电离能的时候也注意看占据情况。比如应该从s轨道电离的,但实际计算的时候是从d轨道电离的(因为初猜轨道能量不合理所致),那就错了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4475#
 楼主 Author| 发表于 Post on 2017-9-20 12:16:03 | 只看该作者 Only view this author
tj2m 发表于 2017-9-20 12:05
%chk=Gd.chk
#P opt b3lyp/genecp

体系里都没有N、C,定义它们的基组肯定出错
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

13

帖子

0

威望

55

eV
积分
68

Level 2 能力者

4476#
发表于 Post on 2017-9-20 13:10:43 | 只看该作者 Only view this author
本帖最后由 tj2m 于 2017-9-20 16:28 编辑
sobereva 发表于 2017-9-20 12:16
体系里都没有N、C,定义它们的基组肯定出错

老师,体系有Gd C H O N原子,我是简化了一下,坐标太长了,只写了第一个和最后一个坐标
General basis read from cards:  (5D, 7F)
Centers:   1 48
MWB46
IA or ICore out of range in MWBBas
Error termination via Lnk1e in C:\G09W\l301.exe
Job cpu time:  0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=      9 Int=      0 D2E=      0 Chk=      1 Scr=      1

这个是错误信息

还有一个问题请教一下,对于配合物中的Gd元素,用什么基组、泛函合适呢

52

帖子

0

威望

562

eV
积分
614

Level 4 (黑子)

4477#
发表于 Post on 2017-9-20 16:00:25 | 只看该作者 Only view this author
sobereva 发表于 2017-9-19 21:10
ORCA是通过命令行使用的,不是双击执行的

谢谢老师,自己查了琢磨一上午终于可以用了

54

帖子

0

威望

354

eV
积分
408

Level 3 能力者

4478#
发表于 Post on 2017-9-20 16:06:18 | 只看该作者 Only view this author
sobereva 发表于 2017-9-20 12:09
计算电离能的时候也注意看占据情况。比如应该从s轨道电离的,但实际计算的时候是从d轨道电离的(因为初猜 ...

sob老师您好,我查看了一下铁原子的电子排布,发现4s轨道只有一个电子,而有两个d轨道是双电子,这样的话铁原子是3d7 4s1而不是3d6 4s2,为什么计算会出现这样的错误,应该如何修改才能使铁原子变成正常的电子排布啊。我估计这就是我电离能计算错误的原因。。。因为Fe→Fe+就变成了失去第一个d电子了。

7

帖子

0

威望

51

eV
积分
58

Level 2 能力者

4479#
发表于 Post on 2017-9-20 16:24:20 | 只看该作者 Only view this author
老师您好:
我用gaussian09优化一个多酸结构,输入文件中的各个有相互作用原子是连接在一起的,可是输出文件用gview打开却有一部分原来相连的原子不连接了,这表明在优化过后这些原子没有相互作用时结构是最低能量,还是可以通过哪个命令来避免这个情况发生?
                                                                                                                                                                                                   谢谢

54

帖子

0

威望

354

eV
积分
408

Level 3 能力者

4480#
发表于 Post on 2017-9-20 16:25:19 | 只看该作者 Only view this author
sob老师您好,还有个小问题,基态吸收光能到达激发态,在垂直激发的过程中,体系的自旋多重度始终是不变的是吧,比如基态是三重态,那垂直跃迁后也只能到达三重态。

49

帖子

0

威望

282

eV
积分
331

Level 3 能力者

4481#
发表于 Post on 2017-9-20 17:09:49 | 只看该作者 Only view this author

高斯中进行相互作用能计算

    在进行相互作用能计算的时候加了BSSE校正,但是计算出来的,我从gif文件中看到有好多组E总能量和E(BSSE),这个我应该选择哪一个能量呀!

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4482#
 楼主 Author| 发表于 Post on 2017-9-20 18:30:14 | 只看该作者 Only view this author
戴上墨镜看世界 发表于 2017-9-20 17:09
在进行相互作用能计算的时候加了BSSE校正,但是计算出来的,我从gif文件中看到有好多组E总能量和E(BSS ...

谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4483#
 楼主 Author| 发表于 Post on 2017-9-20 18:30:23 | 只看该作者 Only view this author
万山红遍 发表于 2017-9-20 16:25
sob老师您好,还有个小问题,基态吸收光能到达激发态,在垂直激发的过程中,体系的自旋多重度始终是不变的 ...

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4484#
 楼主 Author| 发表于 Post on 2017-9-20 18:32:15 | 只看该作者 Only view this author
夏洛特 发表于 2017-9-20 16:24
老师您好:
我用gaussian09优化一个多酸结构,输入文件中的各个有相互作用原子是连接在一起的,可是输出文 ...

没显示成键不代表没有相互作用,只不过此时键长比gview默认的判断标准长而已
只要计算级别合理,优化出什么样就是什么样。是否存在相互作用,存在什么样的相互作用,应当通过Multiwfn进行考察,看此文5.2节
Multiwfn波函数分析程序的意义、功能与用途
http://sobereva.com/184
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125259

管理员

公社社长

4485#
 楼主 Author| 发表于 Post on 2017-9-20 18:33:16 | 只看该作者 Only view this author
万山红遍 发表于 2017-9-20 16:06
sob老师您好,我查看了一下铁原子的电子排布,发现4s轨道只有一个电子,而有两个d轨道是双电子,这样的话 ...

用guess(save,only)计算一遍,这会把初猜轨道保存到chk里。观看轨道,如果一开始就是错的,通过guess=alter调整占据方式使初猜轨道和实际一致。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

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

GMT+8, 2026-2-28 23:07 , Processed in 0.292287 second(s), 22 queries , Gzip On.

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