计算化学公社

 找回密码 Forget password
 注册 Register
Views: 486|回复 Reply: 11

[新手求助] 硅宾SiH2=Si+H2反应过程过渡态的搜索求助

[复制链接 Copy URL]

50

帖子

0

威望

315

eV
积分
365

Level 3 能力者

发表于 Post on 2022-1-13 23:07:48 | 显示全部楼层 Show all |阅读模式 Reading model
最近在搜索硅宾类物质SiH2、SiHCl、SiCl2分解成Si的反应过渡态,但是目前只搜索出来了SiHCl=Si+HCl的过渡态,尝试了多种SiH2=Si+H2的可能过渡态,但是都计算出错,遂希望公社里大佬们给点意见。
这是我搜索SiHCl时的gjf文件
%rwf=TS1.rwf
%NoSave
%chk=TS1.chk
%nprocshared=64
%mem=5GB
#p opt=(calcall,ts,noeigen,gdiis) freq ub3lyp def2tzvp em=gd3bj


G1TS


0 3
Cl                 0.03251700   -1.14117600    0.00000000
Si                 0.03251700    1.39437400    0.00000000
H                 -1.00803400   -0.12124000    0.00000000



但是当我把Cl换成H之后却死活搜不出来过渡态,是我使用的方法有问题吗?还是我过渡态的初猜没有设好?
希望有大佬能给点意见,谢谢!
附件里给出了我搜索SiHCl过渡态时的gjf文件和log文件,以及我搜SiH2时的gjf文件和log文件。

SiHClTS1.gjf

326 Bytes, 下载次数 Times of downloads: 3

SiHClTS1.log

137.34 KB, 下载次数 Times of downloads: 2

SiH2TS3.gjf

318 Bytes, 下载次数 Times of downloads: 4

SiH2TS3.7z

52.34 KB, 下载次数 Times of downloads: 1

log文件

3万

帖子

99

威望

4万

eV
积分
81475

管理员

公社社长+计算化学玩家

发表于 Post on 2022-1-14 06:49:18 | 显示全部楼层 Show all
做个对应解离的柔性扫描试试
另外,SiH2基态是单重态,你当三重态算不妥
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

50

帖子

0

威望

315

eV
积分
365

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-1-14 11:33:32 | 显示全部楼层 Show all
sobereva 发表于 2022-1-14 06:49
做个对应解离的柔性扫描试试
另外,SiH2基态是单重态,你当三重态算不妥

sob老师好,算单重态的话,SiHCl搜不出来过渡态了,气相硅三重态更稳定一些,所以我当三重态算的。

77

帖子

0

威望

380

eV
积分
457

Level 3 能力者

发表于 Post on 2022-1-14 11:41:17 | 显示全部楼层 Show all
本帖最后由 dnlx 于 2022-1-14 12:06 编辑

过渡态会不会类似于甲醛分解成CO和氢气(初级班有这个例子)?如果二者相似,过渡态应该是单重态,一个Si-H键拉长,另一个H在該键上面,呈钝角三角形(补充:就是你摆的SiHCl过渡态这个形状)。

38

帖子

0

威望

893

eV
积分
931

Level 4 (黑子)

发表于 Post on 2022-1-14 13:05:17 | 显示全部楼层 Show all
这个东西有很多文献都算过了,如 Int. J. Quantum Chem. 2013, 113, 1735-1746.,你根据文献中的数据构建结构,事半功倍。

3万

帖子

99

威望

4万

eV
积分
81475

管理员

公社社长+计算化学玩家

发表于 Post on 2022-1-14 13:07:31 | 显示全部楼层 Show all
shiguang 发表于 2022-1-14 11:33
sob老师好,算单重态的话,SiHCl搜不出来过渡态了,气相硅三重态更稳定一些,所以我当三重态算的。

你研究的是解离过程,原本的SiH2是单重态,对于找过渡态,显然就算只考虑一个自旋态,那也应当是考虑单重态
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

50

帖子

0

威望

315

eV
积分
365

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-1-17 11:12:20 | 显示全部楼层 Show all
本帖最后由 shiguang 于 2022-1-17 11:20 编辑
joeson 发表于 2022-1-14 13:05
这个东西有很多文献都算过了,如 Int. J. Quantum Chem. 2013, 113, 1735-1746.,你根据文献中的数据构建结 ...

这个文献里说他用CCSD(T)/6-311+G(d,p)优化结构,貌似CCSD(T)算不了结构优化
我用DFT方法试了试,发现无论是B3LYP、M06-2X或者wB97XD全部报错,都优化不出来,是我用的泛函不对吗?
基组我都是用的def2TZVP
还有其他的能查到的文献吗?
谢谢大佬回复! 1642389239(1).jpg

TS3M062X.gjf

332 Bytes, 下载次数 Times of downloads: 1

TS3M062X.log

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

TS3wb97xd.gjf

337 Bytes, 下载次数 Times of downloads: 0

TS3wb97xd.7z

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

TS3B3LYP.gjf

337 Bytes, 下载次数 Times of downloads: 0

TS3B3LYP.log

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

3万

帖子

99

威望

4万

eV
积分
81475

管理员

公社社长+计算化学玩家

发表于 Post on 2022-1-17 11:19:45 | 显示全部楼层 Show all
shiguang 发表于 2022-1-17 11:12
这个文献里说他用CCSD(T)/6-311+G(d,p)优化结构,貌似CCSD(T)算不了结构优化
我用DFT方法试了试,发现无 ...

给CCSD(T)用6-311+G(d,p)基组是严重缺乏常识的行为
几何优化用CCSD(T)比CCSD几乎没有实际益处,在Gaussian里后者还没解析梯度而只能用数值梯度来优化

Clipboard01.png

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

50

帖子

0

威望

315

eV
积分
365

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-1-17 11:50:00 | 显示全部楼层 Show all
sobereva 发表于 2022-1-17 11:19
给CCSD(T)用6-311+G(d,p)基组是严重缺乏常识的行为
几何优化用CCSD(T)比CCSD几乎没有实际益处,在Gaussi ...

好的,sob老师,我试试

50

帖子

0

威望

315

eV
积分
365

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-1-17 13:57:03 | 显示全部楼层 Show all
joeson 发表于 2022-1-14 13:05
这个东西有很多文献都算过了,如 Int. J. Quantum Chem. 2013, 113, 1735-1746.,你根据文献中的数据构建结 ...

严重怀疑这篇文献数据瞎编的,我用CCSD(T)/6-311+G(d,p)也没有优化出来

CCSDT.gjf

400 Bytes, 下载次数 Times of downloads: 1

CCSDT.7z

165.53 KB, 下载次数 Times of downloads: 2

38

帖子

0

威望

893

eV
积分
931

Level 4 (黑子)

发表于 Post on 2022-1-19 14:54:53 | 显示全部楼层 Show all
shiguang 发表于 2022-1-17 11:12
这个文献里说他用CCSD(T)/6-311+G(d,p)优化结构,貌似CCSD(T)算不了结构优化
我用DFT方法试了试,发现无 ...

我认为当前结构具有静态相关性,一直报错可以试试用纯泛函计算,然后以结果为初猜,在稍好的级别下计算。那篇文献确实有问题,应该是用蓝拱门算的。其他文献我也是一扫而过,不可能记得太清楚。

50

帖子

0

威望

315

eV
积分
365

Level 3 能力者

 楼主 Author| 发表于 Post on 2022-1-20 11:15:47 | 显示全部楼层 Show all
joeson 发表于 2022-1-19 14:54
我认为当前结构具有静态相关性,一直报错可以试试用纯泛函计算,然后以结果为初猜,在稍好的级别下计算。 ...

好的,谢谢joeson

本版积分规则 Credits rule

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

GMT+8, 2022-5-18 00:48 , Processed in 0.402620 second(s), 26 queries .

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