计算化学公社

标题: Gaussian能够同时限制性优化加上过渡态扫描吗 [打印本页]

作者
Author:
liurui    时间: 2022-6-20 08:44
标题: Gaussian能够同时限制性优化加上过渡态扫描吗
想要限制两原子之间距离进行过渡态搜索,需要定义片段fragment=3来指定重态和电荷,但是报错,在fragment这里

作者
Author:
zjxitcc    时间: 2022-6-20 09:17
任务类型太多,需要分开做,采用多步任务(--Link1--)。第一步先构造片段组合波函数初猜,第二步读取轨道,将初猜的SCF收敛(最好还用上stable=opt,检验波函数稳定性并确保波函数稳定),第三步读取轨道,做限制性优化(可以用opt=ts找过渡态)。

你当前展示内容没有“限制两原子之间距离”的关键词。基组有点辣眼睛,看《谈谈量子化学中基组的选择》http://sobereva.com/336
作者
Author:
liurui    时间: 2022-6-20 09:44
zjxitcc 发表于 2022-6-20 09:17
任务类型太多,需要分开做,采用多步任务(--Link1--)。第一步先构造片段组合波函数初猜,第二步读取轨道 ...

老师,分了三步做,是这个意思吗?如下:
#p ub3lyp/6-31g(d) guess=fragment=2 empiricaldispersion=gd3 scf=(sleazy,xqc)
#p ub3lyp 6-311+g(2df,2pd) scf=xqc guess=read scf=nosymm geom=allcheck empiricaldispersion=gd3 stable=opt
#p opt=(modredundant,calcfc,ts,noeigentest) freq=raman ub3lyp 6-311+g(2df,2pd) empiricaldispersion=gd3 iop(1/8=5) scf=xqc guess=read scf=nosymm geom=allcheck

B 1 2 S 15 -0.030000



作者
Author:
zjxitcc    时间: 2022-6-20 09:47
liurui 发表于 2022-6-20 09:44
老师,分了三步做,是这个意思吗?如下:
#p ub3lyp/6-31g(d) guess=fragment=2 empiricaldispersion=gd ...

(1)不要加多余的关键词(像scf=nosymm这种),scf=sleazy尤其可怕。
(2)你在1L里不是写的fragment=3么,怎么到了这里变成2块了。
(3)你这限制性优化,还找过渡态,却去算拉曼光谱freq=raman,这样算出的光谱没有意义。


作者
Author:
liurui    时间: 2022-6-20 09:55
zjxitcc 发表于 2022-6-20 09:47
(1)不要加多余的关键词(像scf=nosymm这种),scf=sleazy尤其可怕。
(2)你在1L里不是写的fragment=3 ...

第一步:
#p ub3lyp/6-31g(d) guess=fragment=3 empiricaldispersion=gd3

Required

0 1 0 2 0 -2 0 1
第二步:
#p ub3lyp 6-311+g(2df,2pd) scf=xqc guess=read geom=allcheck empiricaldispersion=gd3 stable=opt
第三步:
#p opt=(modredundant,calcfc,ts,noeigentest) ub3lyp 6-311+g(2df,2pd) empiricaldispersion=gd3 iop(1/8=5) scf=xqc guess=read scf=nosymm geom=allcheck

B 1 2 S 15 -0.030000
老师,这样可以吗,谢谢老师
作者
Author:
zjxitcc    时间: 2022-6-20 10:36
liurui 发表于 2022-6-20 09:55
第一步:
#p ub3lyp/6-31g(d) guess=fragment=3 empiricaldispersion=gd3

第三步多了一个scf=nosymm
作者
Author:
liurui    时间: 2022-6-20 10:41
zjxitcc 发表于 2022-6-20 10:36
第三步多了一个scf=nosymm

收到,谢谢老师




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3