计算化学公社

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

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

   关闭 [复制链接 Copy URL]

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

2311#
发表于 Post on 2016-10-15 23:29:53 | 只看该作者 Only view this author
sobereva 发表于 2016-10-15 20:54
优化过程用不着def2这么大的
算单点的时候可以都用def2,优化时候Cl可以用6-311G**,Br、I用lanl08(d)

这样的话,三个元素基组不同会不会对反应有较大影响?lanl08(d)这个d是加了极化函数的意思吧?

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

2312#
 楼主 Author| 发表于 Post on 2016-10-16 00:00:07 | 只看该作者 Only view this author
小书童 发表于 2016-10-15 23:29
这样的话,三个元素基组不同会不会对反应有较大影响?lanl08(d)这个d是加了极化函数的意思吧?[/backcolo ...

对于优化没什么影响,本身轻重原子很难统一,也没必要统一,级别差不多就行了。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
小书童 + 1 谢谢

查看全部评分 View all ratings

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

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

2313#
发表于 Post on 2016-10-16 10:25:16 | 只看该作者 Only view this author
sob老师,请教您一个问题
我用B2PLYP双杂化泛函计算一个含铀络合物的单点能(已优化),末尾如下
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Largest valence mixing into a core orbital is  5.49D-01
Largest core mixing into a valence orbital is  4.33D-01
Excessive mixing of frozen core and valence orbitals.
Error termination via Lnk1e in /db/home/renym/g09/l801.exe at Sun Oct 16 00:11:44 2016.
Job cpu time:       0 days  1 hours  7 minutes 44.6 seconds.
File lengths (MBytes):  RWF=    277 Int=      0 D2E=      0 Chk=      5 Scr=      2

请问这是怎么回事?如果把B2PLYP换成M062X其他不变,计算就会成功。

谢谢老师

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

2314#
 楼主 Author| 发表于 Post on 2016-10-16 17:48:52 | 只看该作者 Only view this author
lastzealot 发表于 2016-10-16 10:25
sob老师,请教您一个问题
我用B2PLYP双杂化泛函计算一个含铀络合物的单点能(已优化),末尾如下
Defaul ...


甭用冻核近似
B2PLYP里写上full
北京科音自然科学研究中心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

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

2315#
发表于 Post on 2016-10-16 19:00:08 | 只看该作者 Only view this author
sobereva 发表于 2016-10-16 00:00
对于优化没什么影响,本身轻重原子很难统一,也没必要统一,级别差不多就行了。

多谢老师指点,明白怎么选择

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

2316#
发表于 Post on 2016-10-16 22:16:09 | 只看该作者 Only view this author
老师,为什么经常说虚频很小的过渡态不是真正的过渡态呢?我找到了一个-40的过渡态,趋势是我要找的,IRC也能连接到反应物和产物,但是在走IRC的时候我设定了stepsize=15,那么这个过渡态对吗?

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

2317#
 楼主 Author| 发表于 Post on 2016-10-16 23:33:13 | 只看该作者 Only view this author
小书童 发表于 2016-10-16 22:16
老师,为什么经常说虚频很小的过渡态不是真正的过渡态呢?我找到了一个-40的过渡态,趋势是我要找的,IRC也 ...


那是因为虚频非常小的往往是数值噪音所致,诸如对应体系整体运动。或者对应非常柔的分子内振动模式的过渡态,这肯定不对应于化学反应。

但也不排除极为个别情况,反应势垒就是特别浅,虚频特别小。对于这类情况,优化时候注意用tight、ultrafine,保证数值精度。并且也建议多尝试一些方法,以确认确实有这个过渡态。
北京科音自然科学研究中心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

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

2318#
发表于 Post on 2016-10-17 08:29:17 | 只看该作者 Only view this author
sobereva 发表于 2016-10-16 23:33
那是因为虚频非常小的往往是数值噪音所致,诸如对应体系整体运动。或者对应非常柔的分子内振动模式的过 ...

谨记老师的教诲,我在检验一下,是否对。我觉得对于虚频小的过渡态走IRC似乎也是个问题,经常2、3步就停止了。

4

帖子

0

威望

19

eV
积分
23

Level 1 能力者

2319#
发表于 Post on 2016-10-17 09:53:47 | 只看该作者 Only view this author
优化结构时,总是输出报错(l9999),把结果重新计算还是报错,请问该怎么解决?

3

帖子

0

威望

27

eV
积分
30

Level 2 能力者

2320#
发表于 Post on 2016-10-17 11:13:46 | 只看该作者 Only view this author
sob老师,请教您一个问题:
如果想计算O2和团簇相互作用形成的体系的总能量,但是由于O2是三线态,而那个团簇是单线态(并且带电荷),整个体系的线态和电荷应该如何设置?

128

帖子

0

威望

1486

eV
积分
1614

Level 5 (御坂)

2321#
发表于 Post on 2016-10-17 15:43:07 | 只看该作者 Only view this author
老师,请教您一下,我在做大气化学关于Criegee中间体的有关计算,最简单的是CH2OO这样子的结构,文献表明它的最稳定状态为开壳层单重态,这应该就要考虑对称破缺了吧,于是我用GUASSIAN做了结构优化用UM06-2X/6-311G(d,p),关键词加了guess=mix。计算完成后按照您在个问题中答复写的“<S**2>,不是闭壳层单重态的话,数值肯定>0。或者看末尾的Mulliken自旋布居输出,得到自旋极化单重态的话数值肯定不会都是0”看了一下,发现我的结果中这两个值都为0.是不是我理解错了?谢谢您

286

帖子

0

威望

4688

eV
积分
4974

Level 6 (一方通行)

2322#
发表于 Post on 2016-10-17 15:46:53 | 只看该作者 Only view this author
Gaussian09是不是目前只支持B2PLYP和mPW2PLYP这两种双杂化方法?700多个基函数的过渡金属配合物的单点能如果CCSD算不动是不是可以用这两种方法?(PBE0优化的基态结构)

720

帖子

0

威望

1252

eV
积分
1972

Level 5 (御坂)

2323#
发表于 Post on 2016-10-17 18:31:07 | 只看该作者 Only view this author
sobereva 发表于 2016-10-16 17:48
甭用冻核近似
B2PLYP里写上full

您的意思是不是说泛函那里写成如下形式
#p b2plyp(full)/genecp ...

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

2324#
 楼主 Author| 发表于 Post on 2016-10-18 04:31:51 | 只看该作者 Only view this author
lastzealot 发表于 2016-10-17 18:31
您的意思是不是说泛函那里写成如下形式
#p b2plyp(full)/genecp ...

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

管理员

公社社长

2325#
 楼主 Author| 发表于 Post on 2016-10-18 04:32:13 | 只看该作者 Only view this author
dreamyeye 发表于 2016-10-17 15:46
Gaussian09是不是目前只支持B2PLYP和mPW2PLYP这两种双杂化方法?700多个基函数的过渡金属配合物的单点能如 ...

Gaussian中非内置的理论方法和泛函的用法
http://sobereva.com/344

可以
北京科音自然科学研究中心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-25 18:03 , Processed in 0.216043 second(s), 28 queries , Gzip On.

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