计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3668|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] Gaussian能够同时限制性优化加上过渡态扫描吗

[复制链接 Copy URL]

262

帖子

0

威望

1508

eV
积分
1770

Level 5 (御坂)

想要限制两原子之间距离进行过渡态搜索,需要定义片段fragment=3来指定重态和电荷,但是报错,在fragment这里

捕获.JPG (18.19 KB, 下载次数 Times of downloads: 13)

捕获.JPG

4289

帖子

4

威望

9538

eV
积分
13907

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2022-6-20 09:17:32 | 只看该作者 Only view this author
任务类型太多,需要分开做,采用多步任务(--Link1--)。第一步先构造片段组合波函数初猜,第二步读取轨道,将初猜的SCF收敛(最好还用上stable=opt,检验波函数稳定性并确保波函数稳定),第三步读取轨道,做限制性优化(可以用opt=ts找过渡态)。

你当前展示内容没有“限制两原子之间距离”的关键词。基组有点辣眼睛,看《谈谈量子化学中基组的选择》http://sobereva.com/336
自动做多参考态计算的程序MOKIT

262

帖子

0

威望

1508

eV
积分
1770

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2022-6-20 09:44:27 | 只看该作者 Only view this author
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


4289

帖子

4

威望

9538

eV
积分
13907

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2022-6-20 09:47:55 | 只看该作者 Only view this author
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,这样算出的光谱没有意义。

评分 Rate

参与人数
Participants 2
eV +9 收起 理由
Reason
liurui + 4 谢谢
lonemen + 5 牛!

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

262

帖子

0

威望

1508

eV
积分
1770

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2022-6-20 09:55:52 | 只看该作者 Only view this author
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
老师,这样可以吗,谢谢老师

4289

帖子

4

威望

9538

eV
积分
13907

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2022-6-20 10:36:02 | 只看该作者 Only view this author
liurui 发表于 2022-6-20 09:55
第一步:
#p ub3lyp/6-31g(d) guess=fragment=3 empiricaldispersion=gd3

第三步多了一个scf=nosymm
自动做多参考态计算的程序MOKIT

262

帖子

0

威望

1508

eV
积分
1770

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2022-6-20 10:41:07 | 只看该作者 Only view this author
zjxitcc 发表于 2022-6-20 10:36
第三步多了一个scf=nosymm

收到,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 18:49 , Processed in 0.202853 second(s), 24 queries , Gzip On.

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