计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2341#
 楼主 Author| 发表于 Post on 2023-10-24 01:53:22 | 只看该作者 Only view this author
吴双 发表于 2023-10-23 11:47
想跟老师们请教一下,高斯优化结构的时候出现scf不收敛时,可以将log保存成新的gjf加上关键词重新优化么, ...

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

威望

5万

eV
积分
124671

管理员

公社社长

2342#
 楼主 Author| 发表于 Post on 2023-10-24 01:56:14 | 只看该作者 Only view this author
化小白 发表于 2023-10-23 11:59
老师您好,我在b3lyp/6-31G(d)水平上优化找到一个过渡态,走IRC出来的反应物和产物也确实是想要的,只优化 ...

如果是有机反应,用M06-2X或wB97XD结合6-311G**去优化过渡态,如果都优化不出来(最好用之前B3LYP/6-31G*找到的过渡态结构当初猜),大概率说明这个过渡态不存在、之前B3LYP/6-31G*优化出来的是虚假的。
北京科音自然科学研究中心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

威望

5万

eV
积分
124671

管理员

公社社长

2343#
 楼主 Author| 发表于 Post on 2023-10-24 03:52:17 | 只看该作者 Only view this author
ghkjhug 发表于 2023-10-23 22:20
在计算路易斯酸对过渡态的作用中,发现结构优化能量下降十分剧烈,结构严重不合理。
考虑过是初猜不合理导 ...

你定义基组时漏了硼的
北京科音自然科学研究中心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

31

帖子

0

威望

244

eV
积分
275

Level 3 能力者

2344#
发表于 Post on 2023-10-24 07:46:29 | 只看该作者 Only view this author
sobereva 发表于 2023-10-24 03:52
你定义基组时漏了硼的

好的,这个例子恰好忘了定义硼,不过之前也出现过定义了也没用的情况,我重新发一下

125

帖子

0

威望

2838

eV
积分
2963

Level 5 (御坂)

2345#
发表于 Post on 2023-10-24 09:26:43 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-23 21:56
有的键本来就不能用只切断这一根键的方法算出BDE,而必须用例如MTA等方法(http://sobereva.com/522)近 ...

好的老师,对于您提到的第二种做法,我这么理解是否正确:将体系作为三重态进行结构优化和频率分析,再对优化好的结构进行自旋极化单重态单点能计算,之后再计算BDE或其他热力学量。
如果是这样的话,三重态结构与自旋极化单重态结构的差异可以忽略不计吗?

26

帖子

0

威望

109

eV
积分
135

Level 2 能力者

2346#
发表于 Post on 2023-10-24 10:38:29 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-22 16:25
1. 注意前提条件,仅当不是非加弥散函数不可的时候,更高的基组才比加弥散更优
2. 用CASPT2或SF-TDDFT

好的,谢谢老师

166

帖子

0

威望

640

eV
积分
806

Level 4 (黑子)

2347#
发表于 Post on 2023-10-24 15:35:33 | 只看该作者 Only view this author
sobereva 发表于 2023-10-24 01:56
如果是有机反应,用M06-2X或wB97XD结合6-311G**去优化过渡态,如果都优化不出来(最好用之前B3LYP/6-31G* ...

老师,是有机反应,但是有过渡金属,所以一直没敢用M06-2X去优化

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

2348#
发表于 Post on 2023-10-24 16:21:27 | 只看该作者 Only view this author
MercuryLamp 发表于 2023-10-24 02:26
好的老师,对于您提到的第二种做法,我这么理解是否正确:将体系作为三重态进行结构优化和频率分析,再对 ...

单点能也用三重态计算。这样算出来的是断键+自旋翻转的总能量变化,但是只要键完全断了,自旋翻转对能量的变化应当可以忽略不计,所以这样的结果虽然不如断键直接优化结构那样严格,但也是有参考价值的。
如果结构优化用三重态,单点能用自旋极化单重态,得到的结果或许数值上是合理的,但没有明确物理意义
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
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

109

eV
积分
135

Level 2 能力者

2349#
发表于 Post on 2023-10-24 16:25:10 | 只看该作者 Only view this author
sobereva 发表于 2023-10-22 02:44
不收敛和不能用毫无直接关系
不收敛又不是calcall、recalc造成的,不用的时候照样没收敛
优化极小点甭 ...

老师请问,激发态势能面交叉,通过调整步长可能会到达不同的势能面,那是否意味着可以通过这种方式对S1的激发态结构进行优化。只要出现两种不同的优化能量就可以。

125

帖子

0

威望

2838

eV
积分
2963

Level 5 (御坂)

2350#
发表于 Post on 2023-10-24 16:38:40 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-24 16:21
单点能也用三重态计算。这样算出来的是断键+自旋翻转的总能量变化,但是只要键完全断了,自旋翻转对能量 ...

好的,谢谢老师

26

帖子

0

威望

109

eV
积分
135

Level 2 能力者

2351#
发表于 Post on 2023-10-24 16:44:41 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-22 16:25
1. 注意前提条件,仅当不是非加弥散函数不可的时候,更高的基组才比加弥散更优
2. 用CASPT2或SF-TDDFT

麻烦请问用CASPT2或SF-TDDFT计算,用什么软件合适,我目前用高斯,但是高斯好像不能做CASPT2或SF-TDDFT计算。

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

2352#
发表于 Post on 2023-10-24 20:06:04 | 只看该作者 Only view this author
木子叶 发表于 2023-10-24 09:44
麻烦请问用CASPT2或SF-TDDFT计算,用什么软件合适,我目前用高斯,但是高斯好像不能做CASPT2或SF-TDDFT计 ...

用Molpro、Q-Chem等
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
本团队长期招收研究生,有意者可私信联系

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

2353#
发表于 Post on 2023-10-24 20:07:55 | 只看该作者 Only view this author
木子叶 发表于 2023-10-24 09:25
老师请问,激发态势能面交叉,通过调整步长可能会到达不同的势能面,那是否意味着可以通过这种方式对S1的 ...

你优化的目的是什么?
如果优化的目的就是优化S1,那”出现两种不同的优化能量“显然既不是充分条件也不是必要条件。即使优化出一个能量,这个能量也可能是正确的S1能量;即便优化出两个能量,也不保证其中一定有一个是S1,也可能一个是S2另一个是S3。
如果优化的目的不仅仅是优化S1,那你把优化目的说清楚
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
本团队长期招收研究生,有意者可私信联系

166

帖子

0

威望

640

eV
积分
806

Level 4 (黑子)

2354#
发表于 Post on 2023-10-24 20:56:45 | 只看该作者 Only view this author
sobereva 发表于 2023-10-24 01:56
如果是有机反应,用M06-2X或wB97XD结合6-311G**去优化过渡态,如果都优化不出来(最好用之前B3LYP/6-31G* ...

老师,我用b3lyp/6-31+G(d)优化也可以找出来这个过渡态,算出来的是正常的,过渡态吉布斯自由能-反应物吉布斯自由能=1.1kcal/mol,差值和b3lyp/6-31G(d)是一样的

26

帖子

0

威望

109

eV
积分
135

Level 2 能力者

2355#
发表于 Post on 2023-10-24 21:01:12 | 只看该作者 Only view this author

好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 07:45 , Processed in 2.014926 second(s), 21 queries , Gzip On.

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