计算化学公社

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

[综合交流] 如何计算S1/S0的圆锥交叉点?

[复制链接 Copy URL]

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
大家好!
最近研究了一个有45个原子的纯有机分子,想找一下S1/S0的圆锥交叉点,对于CASSCF来说我这个分子应该是跑不动了,
请问还有其他什么方法可以计算其S1/S0的圆锥交叉点吗?谢谢指导。

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

17#
 楼主 Author| 发表于 Post on 2020-12-25 12:31:34 | 只看该作者 Only view this author
pwzhou 发表于 2020-12-25 09:30
正常100圈应该可以收敛的,你看看收敛趋势,如果100圈没收敛但是一直在持续变小,可以把100继续加大

嗯嗯,感谢回复,我研究研究,之前没怎么用过Q-chem,感谢您的指导!

368

帖子

0

威望

4691

eV
积分
5059

Level 6 (一方通行)

16#
发表于 Post on 2020-12-25 09:30:09 | 只看该作者 Only view this author
暖空 发表于 2020-12-25 08:34
好的,感谢您的回复。
如果把最大循环圈调到100最后还是没有收敛的话,请问除了调整初始构型,还有其他 ...

正常100圈应该可以收敛的,你看看收敛趋势,如果100圈没收敛但是一直在持续变小,可以把100继续加大

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2020-12-25 08:34:21 | 只看该作者 Only view this author
pwzhou 发表于 2020-12-24 21:10
MAX_CIS_CYCLES = 100
加上这个,默认循环30圈没有收敛。

好的,感谢您的回复。
如果把最大循环圈调到100最后还是没有收敛的话,请问除了调整初始构型,还有其他什么方法么?
谢谢老师~

368

帖子

0

威望

4691

eV
积分
5059

Level 6 (一方通行)

14#
发表于 Post on 2020-12-24 21:10:49 | 只看该作者 Only view this author
暖空 发表于 2020-12-24 20:07
您好!
我在尝试用Q-chem的spin-flip TDDFT计算S1/S0的圆锥交叉点,根据手册上(Q-chem manual 4.3)给的 ...

MAX_CIS_CYCLES = 100
加上这个,默认循环30圈没有收敛。

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2020-12-24 20:07:45 | 只看该作者 Only view this author
pwzhou 发表于 2020-11-26 13:23
45个原子用(12,12)没啥问题,可以算动,不过别用Gaussian,Gaussian的CASSCF活化空间超过(10,8)收敛太差。F ...

您好!
我在尝试用Q-chem的spin-flip TDDFT计算S1/S0的圆锥交叉点,根据手册上(Q-chem manual 4.3)给的算例,分别用branching-plane update method 和penalty constrained method 两种方法测试了我一个体系分子,但是提交上去没跑多久都报错了,能麻烦您帮忙看一下具体是哪里出问题了么?谢谢。
branching_ts-TRAJ21-438fs.out (17.15 KB, 下载次数 Times of downloads: 6) ts-TRAJ21-438fs.out (16.89 KB, 下载次数 Times of downloads: 3)

4289

帖子

4

威望

9546

eV
积分
13915

Level 6 (一方通行)

MOKIT开发者

12#
发表于 Post on 2020-12-3 10:55:06 | 只看该作者 Only view this author
暖空 发表于 2020-12-3 10:18
我去orca论坛上搜了一下,现在ORCA找交叉点只能用TDDFT。虽然手册上说的可以用CASSCF,但是有BUG,而且在 ...

专业多参考计算,早日投身OpenMolcas/Molpro,打开新世界。
自动做多参考态计算的程序MOKIT

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2020-12-3 10:18:40 | 只看该作者 Only view this author
biogon 发表于 2020-11-26 15:10
用orca的CASSCF应该是算的动的

我去orca论坛上搜了一下,现在ORCA找交叉点只能用TDDFT。虽然手册上说的可以用CASSCF,但是有BUG,而且在最新版的4.2.1中还没修复

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2020-12-3 10:16:03 | 只看该作者 Only view this author
niobium 发表于 2020-11-26 21:03
ORCA可以用TDDFT优化交叉点,靠不靠谱就不知道了

用ORCA的TDDFT确实优化出来了一个交叉点构型,但是得到的构型很不合理

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2020-11-30 19:23:10 | 只看该作者 Only view this author
niobium 发表于 2020-11-26 21:03
ORCA可以用TDDFT优化交叉点,靠不靠谱就不知道了

嗯嗯,我主要是想看一下交叉点的结构

538

帖子

1

威望

5765

eV
积分
6323

Level 6 (一方通行)

8#
发表于 Post on 2020-11-26 21:03:16 | 只看该作者 Only view this author
ORCA可以用TDDFT优化交叉点,靠不靠谱就不知道了

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2020-11-26 19:00:53 | 只看该作者 Only view this author
biogon 发表于 2020-11-26 15:10
用orca的CASSCF应该是算的动的

嗯嗯,我再搜点儿资料了解一下,感谢回复~

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2020-11-26 18:59:15 | 只看该作者 Only view this author
pwzhou 发表于 2020-11-26 13:23
45个原子用(12,12)没啥问题,可以算动,不过别用Gaussian,Gaussian的CASSCF活化空间超过(10,8)收敛太差。F ...

感谢回复~
之前和老板讨论的时候,老板建议用Spin-Flip计算,但是搜了点儿资料,没看太明白所以就搁置了,我再查一下,谢谢指点!

131

帖子

0

威望

927

eV
积分
1058

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-11-26 18:54:48 | 只看该作者 Only view this author
喵星大佬 发表于 2020-11-26 11:40
跑动跑不动看你的活性空间怎么选了,太大的共轭体导致过多轨道要进入活性空间的话应该是不行

不过一般情 ...

S1和S0的构型变化不太大,主要差异在于一个苯环发生了点扭曲,组内有orca和molcas,我再搜点儿资料研究研究,感谢指导

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

4#
发表于 Post on 2020-11-26 15:10:31 | 只看该作者 Only view this author
本帖最后由 biogon 于 2020-11-26 15:13 编辑

用orca的CASSCF应该是算的动的

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

GMT+8, 2026-2-22 10:26 , Processed in 0.200149 second(s), 24 queries , Gzip On.

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