计算化学公社

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

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

   关闭 [复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4921#
 楼主 Author| 发表于 Post on 2017-10-28 21:18:11 | 只看该作者 Only view this author
niubaotian85 发表于 2017-10-28 21:04
sob老师,请问:
1.在使用homo-lumo gap讨论团簇是否稳定的时候,是用优化结构得到的gap好,还是在优化基 ...

1 如果优化和单点级别一样,得到的gap也一样
2 HF成分越大gap越大。纯泛函会低估gap,M06-2X这样HF高达54%的泛函会高估gap。通常建议用B3LYP这样20% HF成份左右的算gap
3 用gap讨论稳定性很不可靠,看这些讨论
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=600
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=1137
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=1207
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=1272
北京科音自然科学研究中心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

99

帖子

0

威望

333

eV
积分
432

Level 3 能力者

4922#
发表于 Post on 2017-10-28 22:11:02 | 只看该作者 Only view this author
sobereva 发表于 2017-10-28 21:18
1 如果优化和单点级别一样,得到的gap也一样
2 HF成分越大gap越大。纯泛函会低估gap,M06-2X这样HF高达5 ...

非常感谢

13

帖子

0

威望

62

eV
积分
75

Level 2 能力者

4923#
发表于 Post on 2017-10-30 10:18:51 | 只看该作者 Only view this author
sob老师好,我跑IRC的时候输入关键词如下:
#p irc=(calcfc,LQA,maxpoints=20,stepsize=10,noeigen) b3lyp/6-31+g(d)
geom=connectivity
但是报错说:
QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.
#p irc=(calcfc,LQA,maxpoints=20,stepsize----
然后有个小角标是指着LQA前面的,的。。我已经dos2unix过了,但是还是跑不起来。。。
不知道这个要怎么改呢?

100

帖子

0

威望

414

eV
积分
514

Level 4 (黑子)

4924#
发表于 Post on 2017-10-30 11:21:07 | 只看该作者 Only view this author
sob老师,用m062x/6-311+G(2d,2p)计算优化和频率,结果显示热力学数据(Sum of electronic and thermal Enthalpies)与实验值吻合度较高,之后换了m062x/def2-tzvpd结合BSSE校正计算(校正因子+单点能),结果反而与前者更差了,现在这种情况能不能用前者的方法来计算呢?

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4925#
 楼主 Author| 发表于 Post on 2017-10-30 15:24:15 | 只看该作者 Only view this author
flumense 发表于 2017-10-30 10:18
sob老师好,我跑IRC的时候输入关键词如下:
#p irc=(calcfc,LQA,maxpoints=20,stepsize=10,noeigen) b3lyp ...

IRC里不能写noeigen,本身IRC就不做本征值检测
北京科音自然科学研究中心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
积分
125245

管理员

公社社长

4926#
 楼主 Author| 发表于 Post on 2017-10-30 15:25:40 | 只看该作者 Only view this author
不懂计算 发表于 2017-10-30 11:21
sob老师,用m062x/6-311+G(2d,2p)计算优化和频率,结果显示热力学数据(Sum of electronic and thermal E ...

如果你用的counterpoise,本身这就不适合用在优化和振动分析场合,又导致耗时极高(没有解析导数)又没有什么意义
北京科音自然科学研究中心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

100

帖子

0

威望

414

eV
积分
514

Level 4 (黑子)

4927#
发表于 Post on 2017-10-30 19:59:21 | 只看该作者 Only view this author
sobereva 发表于 2017-10-30 15:25
如果你用的counterpoise,本身这就不适合用在优化和振动分析场合,又导致耗时极高(没有解析导数)又没有 ...

老师,是用counterpoise计算单点,再读取6-311+G(2d,2p)基组下算的Thermal correction to Enthalpy,二者相加,其结果反而比低精度基组直接计算的更差了,读取的是Sum of electronic and thermal Enthalpies

53

帖子

0

威望

733

eV
积分
786

Level 4 (黑子)

4928#
发表于 Post on 2017-10-30 23:26:21 | 只看该作者 Only view this author
本帖最后由 登登 于 2017-10-30 23:31 编辑

老师,在chemdraw2010里的H显示有自己的规则吗?怎么选择性的在某个C上的H以键的形式显示出来,就像下面图中的碳上的H.


6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4929#
 楼主 Author| 发表于 Post on 2017-10-31 00:04:24 | 只看该作者 Only view this author
不懂计算 发表于 2017-10-30 19:59
老师,是用counterpoise计算单点,再读取6-311+G(2d,2p)基组下算的Thermal correction to Enthalpy,二 ...

counterpoise不是用来计算单点的,是用来计算BSSE校正后的相互作用能的
北京科音自然科学研究中心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

威望

62

eV
积分
75

Level 2 能力者

4930#
发表于 Post on 2017-10-31 09:52:57 | 只看该作者 Only view this author
sobereva 发表于 2017-10-30 15:24
IRC里不能写noeigen,本身IRC就不做本征值检测

可是我之前没写noeigen的时候,跑IRC结果报了l9999的错误:
Optimization stopped.
    -- Wrong number of negative eigenvalues: Desired=  0 Actual=***
    -- Flag reset to prevent archiving.
然后按照常见问题那帖子的说法应该是要加noeigen关键词?
这是我之前报l9999错误的关键词:
#p irc=(calcfc,maxpoints=20,stepsize=10,LQA) b3lyp/6-31+g(d)
geom=connectivity
麻烦sob老师了!

4

帖子

0

威望

35

eV
积分
39

Level 2 能力者

4931#
发表于 Post on 2017-10-31 17:02:20 | 只看该作者 Only view this author
本帖最后由 小时候头很油 于 2017-10-31 17:05 编辑

sob老师 很少见到C-C键的等键反应,这个等键反应合理吗?自己想算生成焓设计了一条
学习!

97

帖子

0

威望

337

eV
积分
434

Level 3 能力者

4932#
发表于 Post on 2017-11-1 12:18:57 | 只看该作者 Only view this author
请问老师,对于一般分子的电子总能量(单位 A.U.),如果进行比较的话,精确到小数点后多少位就可以了呢,因为高斯默认的位数很长,我不知道怎么取舍比较好。

201

帖子

2

威望

1354

eV
积分
1595

Level 5 (御坂)

4933#
发表于 Post on 2017-11-1 13:51:40 | 只看该作者 Only view this author
本帖最后由 steven 于 2017-11-1 13:52 编辑

sob老师好! 请教您两个问题:
1.我用qst2方法找过渡态时,总是出现这个错误,自己不确定是什么问题,输出文件末尾是这样的.
Old curvilinear step not converged, using linear step:
SCX= 1.12D+02 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
RedCar/ORedCr failed for GTrans.
Error termination via Lnk1e in /opt/soft/Gauss/Gauss09.E01/g09/l101.exe at Wed Nov  1 13:48:01 2017.
我确定原子对应关系都是合理的,键连关系我也弄成一致的了.

2.有一个经验性的结论是没有化学意义的怪异的结构,SCF越爱不收敛。这是什么原因?  
随便搭的一个奇怪的结构,虽然化学上是不存在的,但是此时此刻它也是有对应的波函数,只不过是SCF方式难以求得,还是说这结构压根儿就没有波函数?

谢谢老师!!

55

帖子

0

威望

778

eV
积分
833

Level 4 (黑子)

4934#
发表于 Post on 2017-11-2 11:36:42 | 只看该作者 Only view this author

求助

分子的 Ag/Bu 代表什么?

145

帖子

0

威望

3940

eV
积分
4085

Level 6 (一方通行)

4935#
发表于 Post on 2017-11-2 11:41:14 | 只看该作者 Only view this author

本版积分规则 Credits rule

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

GMT+8, 2026-2-28 18:17 , Processed in 0.236671 second(s), 28 queries , Gzip On.

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