计算化学公社

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

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

   关闭 [复制链接 Copy URL]

6万

帖子

99

威望

6万

eV
积分
126450

管理员

公社社长

15751#
 楼主 Author| 发表于 Post on 2022-6-19 10:21:42 | 只看该作者 Only view this author
xiayu1985 发表于 2022-6-18 21:58
两个中间体的能量分别是-1981.643673hartree,-1981.643696hartree.结构差异对比如图所示,主要是左下侧苯环 ...

怀疑有过渡态的话,可以用QST2做一下。就算有过渡态,势垒也可以忽略不计,即实际研究中不需要明确区分这两种构象
北京科音自然科学研究中心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

145

帖子

1

威望

1792

eV
积分
1957

Level 5 (御坂)

15752#
发表于 Post on 2022-6-19 13:17:14 | 只看该作者 Only view this author
请问,在寻找自旋极化单重态的稳定波函数时,能不能用计算速度较快的b3lyp以及玩具机组3-21g,找到后在用我的m062x/6-311g**算过渡态

1万

帖子

0

威望

1万

eV
积分
22423

Level 6 (一方通行)

15753#
发表于 Post on 2022-6-19 15:46:40 | 只看该作者 Only view this author
WOOOOWOOOO 发表于 2022-6-19 06:17
请问,在寻找自旋极化单重态的稳定波函数时,能不能用计算速度较快的b3lyp以及玩具机组3-21g,找到后在用我 ...

可以,前提是必须检查m062x/6-311g**收敛的波函数是否正确。只要最后收敛对了,第一步不管用多差的泛函和基组都没关系,换句话说这个问题只有结果是重要的,结果正确的前提下过程不重要

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
WOOOOWOOOO + 3 牛!

查看全部评分 View all ratings

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

威望

1万

eV
积分
22423

Level 6 (一方通行)

15754#
发表于 Post on 2022-6-19 15:49:58 | 只看该作者 Only view this author
sobereva 发表于 2022-6-19 03:21
怀疑有过渡态的话,可以用QST2做一下。就算有过渡态,势垒也可以忽略不计,即实际研究中不需要明确区分这 ...

如果只是为了找到最稳定的结构,那确实可以不用区分。但如果要做Boltzmann平均的话,这个影响还是比较可观的,因为如果有两个很相近的构象,相比只有一个构象而言,对配分函数的贡献要多一倍。
我最近也遇到了这种问题,我们一般是用RMSD自动判断两个结构是否一致,RMSD<0.125埃就认为是一致的。但是有的时候两个结构的RMSD会到0.11埃左右,所以我一直不放心会不会把两个本来独立的构象当成了一个来算
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
本团队长期招收研究生,有意者可私信联系

34

帖子

0

威望

217

eV
积分
251

Level 3 能力者

15755#
发表于 Post on 2022-6-19 18:03:08 | 只看该作者 Only view this author
请问各位老师用DFT算氧气和其他单线态的物质的反应时  自旋多重度应该怎么考虑   是0 1   /  0   3 都算一次看哪个能量低或者势垒低   还是应该怎么做   ?在相同的初始结构,我用0  3 算的时候反应物有虚频,而用0 1 算就没有虚频   ?这又应该怎么办呢?

1万

帖子

0

威望

1万

eV
积分
22423

Level 6 (一方通行)

15756#
发表于 Post on 2022-6-19 20:06:05 | 只看该作者 Only view this author
cmd27 发表于 2022-6-19 11:03
请问各位老师用DFT算氧气和其他单线态的物质的反应时  自旋多重度应该怎么考虑   是0 1   /  0   3 都算一 ...

用化学常识判断。
如果化学常识不能判断,那就查文献类似体系判断。
如果文献没做过类似体系,那就把所有可能的多重度都算一遍,取能量最低的,但是注意:(1)同一个基元反应的反应物、过渡态和产物的能量最低的自旋多重度可能是不同的,所以不能只挑一个分子算,反应物、过度态和产物都必须分别确定自旋多重度;(2)必须确保波函数是稳定的,尤其当计算单重态时更是如此。
此外,有虚频并不代表这个多重度不对,应该尝试通过扰动结构、把收敛限设严、加大格点等常规方法来消虚频。
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
本团队长期招收研究生,有意者可私信联系

34

帖子

0

威望

217

eV
积分
251

Level 3 能力者

15757#
发表于 Post on 2022-6-19 20:48:32 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-19 20:06
用化学常识判断。
如果化学常识不能判断,那就查文献类似体系判断。
如果文献没做过类似体系,那就把所 ...

老师  您的意思是反应物  过渡态  产物的自旋多重度可以不一样嘛?那不是自旋禁阻嘛?

11

帖子

0

威望

117

eV
积分
128

Level 2 能力者

15758#
发表于 Post on 2022-6-19 22:37:57 | 只看该作者 Only view this author

求助:关于带有自由基的分子和非中性分子的分子轨道如何用高斯画出

本帖最后由 chennanjie 于 2022-6-19 22:39 编辑

是这样的,我会用高斯软件包画中性分子的分子轨道,也就是HOMO和LUMO图。但是我想要画出带有自由基的分子的分子轨道图,却不知从何下手。


在文献中我有见过别人画出类似的图片。我的是中性分子A,发生了两次单电子的氧化反应,分别氧化为+1价的自由基阳离子A·+,和+2价的阳离子A2+

因此我想请教一下大家怎么画,我真的是没有人可以问,本人在一个普通211里面读硕士,反向是做有机液流电池活性材料的,材料计算这块可以说是一窍不通,

希望大家多多帮助,谢谢!
以下是我研究的分子四硫富瓦烯(TTF):



202206192229494651..png (14.02 KB, 下载次数 Times of downloads: 15)

TTF的得失电子反应图

TTF的得失电子反应图

1万

帖子

0

威望

1万

eV
积分
22423

Level 6 (一方通行)

15759#
发表于 Post on 2022-6-20 00:29:21 | 只看该作者 Only view this author
cmd27 发表于 2022-6-19 13:48
老师  您的意思是反应物  过渡态  产物的自旋多重度可以不一样嘛?那不是自旋禁阻嘛?

自旋禁阻说明反应慢,但不说明反应不可能发生。而且即使慢也没慢多少,比如磷光就是自旋禁阻过程,但是寿命其实也就在毫秒量级。
此外还有一种可能是这样的,比如单重态势能面1RC->1TS->1P,三重态势能面3RC->3TS->3P,其中3RC比1RC高,但是3TS比1TS低。所以反应物肯定取1RC,但是1RC和3RC有平衡,可以有很小平衡浓度的3RC存在,而3RC经过3TS反应非常快,拉动1RC和3RC的反应平衡向右进行。这时自由能垒就必须取作G(3TS)-G(1RC)。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
cmd27 + 2 赞!

查看全部评分 View all ratings

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

18

帖子

0

威望

782

eV
积分
800

Level 4 (黑子)

15760#
发表于 Post on 2022-6-20 10:03:31 | 只看该作者 Only view this author

Gaussian提交opt任务后报错,需要大家帮忙!

我想算的是弱相互作用,用到的关键词是#p opt freq m062xd3/def2tzvp scrf=(solvent=water),提交任务后报错
QPErr --- A syntax error was detected in the input line.
#p opt freq m062xd3/def2tzvp scrf=(solve
             '
Last state= "GCL"
TCursr=      3918 LCursr=        12
Error termination via Lnk1e in /softwares/gaussian/g16/l1.exe at Sun Jun 19 21:52:51 2022.
Job cpu time:       0 days  0 hours  0 minutes  0.8 seconds.
Elapsed time:       0 days  0 hours  0 minutes  0.1 seconds.
File lengths (MBytes):  RWF=      6 Int=      0 D2E=      0 Chk=      1 Scr=      1
应该是输入文件的错误,但是不知道怎么改才好,请大佬指教


jtyh-s-PE (4).log

4.41 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

6万

eV
积分
126450

管理员

公社社长

15761#
 楼主 Author| 发表于 Post on 2022-6-20 10:30:56 | 只看该作者 Only view this author
wzkchem5 发表于 2022-6-19 15:49
如果只是为了找到最稳定的结构,那确实可以不用区分。但如果要做Boltzmann平均的话,这个影响还是比较可 ...

所有情况如果都只用RMSD自动判断而不人工检查,对个别情况有点危险。由于RMSD本质上体现的是平均结构差异,倘若一个很大体系只有极小部分结构有差异,诸如C60连一个取代基时,仅取代基局部朝向有异导致不同的构象,二者间RMSD可能很小。如果再加个最大内坐标差异判断指标,会更可靠。
北京科音自然科学研究中心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

4326

帖子

4

威望

9777

eV
积分
14183

Level 6 (一方通行)

MOKIT开发者

15762#
发表于 Post on 2022-6-20 10:36:39 | 只看该作者 Only view this author
M062X em=GD3
自动做多参考态计算的程序MOKIT

145

帖子

1

威望

1792

eV
积分
1957

Level 5 (御坂)

15763#
发表于 Post on 2022-6-20 10:40:44 | 只看该作者 Only view this author
请问,我在输入文件里加个stable老是报错,是什么原因?

输入文件.png (27.45 KB, 下载次数 Times of downloads: 10)

输入文件.png

报错信息.png (38.97 KB, 下载次数 Times of downloads: 8)

报错信息.png

18

帖子

0

威望

782

eV
积分
800

Level 4 (黑子)

15764#
发表于 Post on 2022-6-20 10:44:36 | 只看该作者 Only view this author

修改后可以进行计算了,非常感谢您的帮助

6万

帖子

99

威望

6万

eV
积分
126450

管理员

公社社长

15765#
 楼主 Author| 发表于 Post on 2022-6-20 10:48:29 | 只看该作者 Only view this author
WOOOOWOOOO 发表于 2022-6-20 10:40
请问,我在输入文件里加个stable老是报错,是什么原因?

stable任务别和别的一起做,别图省事

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
WOOOOWOOOO + 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

本版积分规则 Credits rule

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

GMT+8, 2026-4-25 07:45 , Processed in 0.230007 second(s), 26 queries , Gzip On.

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