计算化学公社

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

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

   关闭 [复制链接 Copy URL]

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

17071#
发表于 Post on 2022-11-11 16:47:54 | 只看该作者 Only view this author
应该是算IRC寻找过渡态
http://sobereva.com/400
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

1万

帖子

0

威望

8988

eV
积分
20736

Level 6 (一方通行)

17072#
发表于 Post on 2022-11-11 19:53:51 | 只看该作者 Only view this author
ccx 发表于 2022-11-11 05:52
我想问一下,就是这种大分子结构是那gausview自己点出来的,还是有什么网站可以下载

有若干种做法:
(1)从晶体数据库(例如CCDC)获取结构
(2)从自己实验组里测的晶体结构文件获取结构
(3)用包括gaussview在内的可视化软件画,但是当分子比较大,尤其是有重复单元,或者是生物大分子的时候,其他某些可视化软件可能用起来比gaussview更趁手
(4)从计算的文章的supporting information里获取原子坐标
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120102

管理员

公社社长

17073#
 楼主 Author| 发表于 Post on 2022-11-11 22:10:53 | 只看该作者 Only view this author
Mr.zhen 发表于 2022-11-11 16:28
如题,应当如何通过g16计算反应过程中的相对自由能谱呢?在阅读论文时发现了下图,但文中对此图的计算方法 ...

优化极小点和过渡态,然后算自由能
参考
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460http://bbs.keinsci.com/thread-12093-1-1.html
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327http://bbs.keinsci.com/thread-3345-1-1.html
文中不可能对计算级别完全没说明。有的时候正文里没说但放到补充材料里说。
北京科音自然科学研究中心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

259

帖子

0

威望

849

eV
积分
1108

Level 4 (黑子)

17074#
发表于 Post on 2022-11-12 22:40:13 | 只看该作者 Only view this author
请问寻找过渡态时,opt=(calcfc,noeigentest,ts),默认的算法是RFO,也可以选择GEDIIS,GDIIS(在opt选项里面添加GEDIIS或者GDIIS)。sob在网页http://sobereva.com/44中提及好多方法,不知道这些方法在gjf文件中的输入格式怎么写?例如: P-RFO法,AH方法,QA法,TRIM法,Dimer方法。

1万

帖子

0

威望

8988

eV
积分
20736

Level 6 (一方通行)

17075#
发表于 Post on 2022-11-13 01:07:16 | 只看该作者 Only view this author
lufenghxu 发表于 2022-11-12 15:40
请问寻找过渡态时,opt=(calcfc,noeigentest,ts),默认的算法是RFO,也可以选择GEDIIS,GDIIS(在opt选项里面 ...

http://gaussian.com/opt/
但凡是这里没提到的、但是sob老师提到过的方法,都是高斯不支持、只有其他某些程序才支持的方法
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

259

帖子

0

威望

849

eV
积分
1108

Level 4 (黑子)

17076#
发表于 Post on 2022-11-13 09:40:21 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-13 01:07
看http://gaussian.com/opt/
但凡是这里没提到的、但是sob老师提到过的方法,都是高斯不支持、只有其他 ...

感谢大佬!

93

帖子

1

威望

1399

eV
积分
1512

Level 5 (御坂)

17077#
发表于 Post on 2022-11-14 10:01:37 | 只看该作者 Only view this author
lurensan 发表于 2022-11-11 10:19
老师你好,我直接用的文献中的过渡态坐标,使用文中写的方法结构优化并不收敛。想着直接进行频率分析试着 ...

频率计算只能确认过渡态(必要不充分),不能找到。你用了superfine格点后确实消灭了一个虚频,可能人家真是在这个格点下找的,但可能由于数据精度损失或者软件版本不同,在你这还是有虚频。你可以重新用superfine优化过渡态后计算频率确认看看。

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

17078#
发表于 Post on 2022-11-15 10:17:53 | 只看该作者 Only view this author
本帖最后由 糖糖糖豆9988 于 2022-11-15 10:54 编辑

请教老师们,如何理解异面加成反应呢?定义说的是反应在π键的异侧进行,生成的新键在反应体系的相反两面。如右图显示的异面加成,上面的乙烯利用的是下、上两瓣所以是a,下面的乙烯利用的是上,上两瓣所以是s,所以就是2s+2a环加成反应?

202211151017453531..png (205.89 KB, 下载次数 Times of downloads: 40)

202211151017453531..png

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

17079#
发表于 Post on 2022-11-15 21:01:37 | 只看该作者 Only view this author
请教各位老师,文献里提到的从“核构型”角度看,同面-异面加成反应虽然是热允许的,但是不利。请问核构型角度是什么意思呢?是空间位阻吗?谢谢老师

202211152058118602..png (82.53 KB, 下载次数 Times of downloads: 34)

202211152058118602..png

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

17080#
发表于 Post on 2022-11-15 21:31:08 | 只看该作者 Only view this author
请教各位老师,文献里提到通过计算可以证明2+2环加成反应的逆反应在同面-同面反应时是热禁阻的。但是对于四元环的开环反应,我只知道去寻找协同/一步反应开环的过渡态,但是如何区分是2s+2s的同面-同面,还是2s+2a的同面-异面反应呢?谢谢老师

202211152126401391..png (167.67 KB, 下载次数 Times of downloads: 33)

202211152126401391..png

106

帖子

0

威望

426

eV
积分
532

Level 4 (黑子)

17081#
发表于 Post on 2022-11-16 10:46:06 | 只看该作者 Only view this author
有一个小问题想询问各位老师,
在文章里看到有人优化完结构进行计算,这是这个结构和我直接在PUBCHEM下载下来的有很大出入。
pubchem应该还是挺权威的,这。。。。该怎么判断

18

帖子

0

威望

131

eV
积分
149

Level 2 能力者

17082#
发表于 Post on 2022-11-16 10:48:10 | 只看该作者 Only view this author
老师们,关于发射光谱光谱,我的理解是先计算吸收光谱看从基态跃迁到哪个激发态,然后再算感兴趣的激发态的情况,请问这样理解对吗

160

帖子

0

威望

606

eV
积分
766

Level 4 (黑子)

17083#
发表于 Post on 2022-11-16 11:30:13 | 只看该作者 Only view this author
老师,我是做过渡金属有机催化的,用的是b3lyp/6-31G(d)找的过渡态,这是我找的一个过渡态的过程,这个过渡态最后找失败了,最终变成了反应物,想问一下老师,这个圈里的这一段会是过渡态吗?

EO3EL5UE0`HU@)S707}`6%X.png (36.88 KB, 下载次数 Times of downloads: 16)

EO3EL5UE0`HU@)S707}`6%X.png

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

17084#
发表于 Post on 2022-11-16 11:34:40 | 只看该作者 Only view this author
化小白 发表于 2022-11-16 11:30
老师,我是做过渡金属有机催化的,用的是b3lyp/6-31G(d)找的过渡态,这是我找的一个过渡态的过程,这个过渡 ...

不好说,你有类似体系参考吗?体系有多大?b3lyp/6-31G(d)未必是好级别。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

160

帖子

0

威望

606

eV
积分
766

Level 4 (黑子)

17085#
发表于 Post on 2022-11-16 15:05:54 | 只看该作者 Only view this author
chands 发表于 2022-11-16 11:34
不好说,你有类似体系参考吗?体系有多大?b3lyp/6-31G(d)未必是好级别。

有的老师,有篇文献中跟我的路径一模一样,就是反应物和产物不一样,过渡金属相近,但是人家的过渡态能垒就不高,我的能垒特别高,他用的是M11/6-31G(d)/lan2DZ,我用的是b3lyp/6-31G(d)/SDD,加了色散矫正

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 01:10 , Processed in 0.165279 second(s), 22 queries , Gzip On.

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