计算化学公社

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

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

   关闭 [复制链接 Copy URL]

146

帖子

0

威望

1329

eV
积分
1475

Level 4 (黑子)

6736#
发表于 Post on 2018-5-24 22:40:45 | 只看该作者 Only view this author
本帖最后由 刘小英雄steven 于 2018-5-25 09:13 编辑
sobereva 发表于 2018-5-23 19:09
看到底是SCF部分太慢,还是计算双杂化的MP2部分太慢。如果是后者,用RIJCOSX结合诸如PBE0-D3来算。如果即 ...

Sob老师,我勉强采用了pwpb95/def2tzvp级别用orca计算单点能(关键词!RI-PWPB95 D3  gcp(dft/TZ)  def2-TZVP  def2/J def2-TZVP/C RIJK grid4  nopop ),但是得到的能量约为-2911.47828a.u,而此前我在高斯里面opt得到的能量为约为-1777.26264a.u(b3lyp/6-31g*,fe和pt均为sdd),这个能量差距是不可思议的。

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

6737#
 楼主 Author| 发表于 Post on 2018-5-25 01:39:18 | 只看该作者 Only view this author
刘小英雄steven 发表于 2018-5-24 22:40
Sob老师,我勉强采用了pwpb95/def2tzvp级别用orca计算单点能(关键词!RI-PWPB95 D3  gcp(dft/TZ)  def2-T ...

SDD对Fe是赝势基组,而def2-TZVP对前四周期都是全电子基组,这个差异会极大地影响能量
北京科音自然科学研究中心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

139

帖子

0

威望

473

eV
积分
612

Level 4 (黑子)

6738#
发表于 Post on 2018-5-25 09:01:43 | 只看该作者 Only view this author

高斯输出文件l1.exe报错

跑IRC,出现了下面情况的报错,想请教一下为什么
%chk=i1.chk
%mem=16GB
-----------------------------------------------
#p irc=(calcfc,lqa) PM7 scrf=(smd,solvent=c6h6)
-----------------------------------------------
QPErr --- A syntax error was detected in the input line.
#p irc=(calcfc,lqa) PM7 scrf=(smd,solven
                     '
Last state= "GCL"
TCursr=      3523 LCursr=        20
Error termination via Lnk1e in /opt/soft/Gauss/Gauss09.E01/g09/l1.exe at Fri May 25 08:38:42 2018.
Job cpu time:       0 days  0 hours  0 minutes  0.6 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

541

帖子

3

威望

6622

eV
积分
7223

Level 6 (一方通行)

6739#
发表于 Post on 2018-5-25 09:52:01 | 只看该作者 Only view this author
本帖最后由 ABetaCarw 于 2018-5-25 09:55 编辑

solvent=c6h6?
solvent=Benzene

http://sobereva.com/g09/k_scrf.htm 这个网页翻到最后有惊喜
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

113

帖子

0

威望

271

eV
积分
384

Level 3 能力者

6740#
发表于 Post on 2018-5-25 10:11:47 | 只看该作者 Only view this author
本帖最后由 erya123 于 2018-5-25 11:05 编辑

老师,我计算反应历程中其中一个过渡态(单重态波函数不稳定,但是三重态反应能垒高)的自旋破缺态,该反应物带一个单位正电荷,结果如下:
Mulliken charges and spin densities:
               1          2
     1  C   -0.180948  -0.004020
     2  C   -0.174825   0.026676
     3  C   -0.209284  -0.126737
     4  C   -0.179639   0.008406
     5  C   -0.243374  -0.001037
     6  C   -0.214428   0.000381
     7  H    0.196473   0.005517
     8  H    0.211601  -0.004539
    14  H    0.158710   0.000032
    15  H    0.142596  -0.000210
    16  H    0.138531  -0.000194
    17  H    0.154594   0.000093
    18  H   -0.028641  -0.000476
    19  Ni  -0.224051   0.137467
    20  Al   0.621635  -0.001258
Sum of Mulliken charges =   1.00000   0.00000
Mulliken charges and spin densities with hydrogens summed into heavy atoms:
               1          2
     1  C    0.137859  -0.003487
     2  C    0.202802  -0.001613
     3  C   -0.012811  -0.121220
     4  C    0.166578  -0.008479
     5  C    0.057932  -0.001214
     6  C    0.078697   0.000280
    19  Ni  -0.252693   0.136992
    20  Al   0.621635  -0.001258

1.这两个(红色)自旋布局有什么不一样?2.以上计算结果能用来说明单重态波函数不稳定吗?(个人理解好像不符合自旋破缺态特征,但是为什么单重态不稳定)

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

6741#
 楼主 Author| 发表于 Post on 2018-5-25 11:29:19 | 只看该作者 Only view this author
大慧慧 发表于 2018-5-25 09:01
跑IRC,出现了下面情况的报错,想请教一下为什么
%chk=i1.chk
%mem=16GB

PM7在G16才支持
solvent=benzene,溶剂写法别自己编,看手册
北京科音自然科学研究中心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
积分
125197

管理员

公社社长

6742#
 楼主 Author| 发表于 Post on 2018-5-25 11:43:09 | 只看该作者 Only view this author
erya123 发表于 2018-5-25 10:11
老师,我计算反应历程中其中一个过渡态(单重态波函数不稳定,但是三重态反应能垒高)的自旋破缺态,该反应 ...

第二次输出的是把氢原子的数值归到与之相连的重原子后的数值

没法说明。只能说明闭壳层单重态可能不稳定(需要跟闭壳层单重态对比能量)
北京科音自然科学研究中心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

146

帖子

0

威望

1329

eV
积分
1475

Level 4 (黑子)

6743#
发表于 Post on 2018-5-25 12:31:20 | 只看该作者 Only view this author
sobereva 发表于 2018-5-25 01:39
SDD对Fe是赝势基组,而def2-TZVP对前四周期都是全电子基组,这个差异会极大地影响能量

老师,我是不是可以理解为:由于我的研究对象(如吸附能,能垒)看的是相对的能量,所以绝对能量差异意义不大,我的工作依然是合理的。

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

6744#
 楼主 Author| 发表于 Post on 2018-5-25 13:31:21 | 只看该作者 Only view this author
刘小英雄steven 发表于 2018-5-25 12:31
老师,我是不是可以理解为:由于我的研究对象(如吸附能,能垒)看的是相对的能量,所以绝对能量差异意义 ...

北京科音自然科学研究中心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

95

帖子

0

威望

237

eV
积分
332

Level 3 能力者

6745#
发表于 Post on 2018-5-26 09:11:20 | 只看该作者 Only view this author
老师,激发态和过渡态有什么区别?还是一样?

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

6746#
 楼主 Author| 发表于 Post on 2018-5-26 13:29:57 | 只看该作者 Only view this author
lily19830130 发表于 2018-5-26 09:11
老师,激发态和过渡态有什么区别?还是一样?

完全是八竿子打不着的两码事
基态和激发态对应不同势能面,而过渡态是势能面上的一阶鞍点,因此基态可以有过渡态,激发态也可以有过渡态
北京科音自然科学研究中心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

38

帖子

0

威望

135

eV
积分
173

Level 3 能力者

6747#
发表于 Post on 2018-5-26 16:00:18 | 只看该作者 Only view this author

计算CCSDT单点出现502错误怎么解决吗

本帖最后由 猴猴是个胖猴 于 2018-5-26 16:08 编辑

请问老师,计算CCSDT单点出现502错误怎么解决吗?加了scf=qc出现508错误和913错误,加了scf=256,scf=conver=5也不能解决,现在不知道该怎么办。还希望老师能够赐教。


QQ图片20180526160716.png (10.88 KB, 下载次数 Times of downloads: 12)

QQ图片20180526160716.png

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

6748#
发表于 Post on 2018-5-26 16:47:01 | 只看该作者 Only view this author
高斯根本就不能算CCSDT,不要乱写
scf=256是什么鬼。。。

所有能解决SCF不收敛问题的方法都在 http://sobereva.com/61 文章里

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

6749#
 楼主 Author| 发表于 Post on 2018-5-26 16:59:38 | 只看该作者 Only view this author
猴猴是个胖猴 发表于 2018-5-26 16:00
请问老师,计算CCSDT单点出现502错误怎么解决吗?加了scf=qc出现508错误和913错误,加了scf=256,scf=conver ...


看6759L的回复,另外仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331
量子化学中的一些常见不良写法和用词
http://sobereva.com/298
北京科音自然科学研究中心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

113

帖子

0

威望

271

eV
积分
384

Level 3 能力者

6750#
发表于 Post on 2018-5-27 00:06:27 | 只看该作者 Only view this author
老师,经过计算单重态的反应路径能垒较低,反应路径中中间体,过渡态的波函数是不是都必须是稳定的?还有有其他什么情况会导致不稳定呢?两态反应。。。。

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 18:42 , Processed in 1.185255 second(s), 22 queries , Gzip On.

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