计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 关于stable=opt计算后激发态优化的相关问题

关于stable=opt计算后激发态优化的相关问题

查看数: 6149 | 评论数: 5 | 收藏 Add to favorites 2

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

yanshen

发布时间: 2021-9-13 09:30

正文摘要:

本帖最后由 yanshen 于 2021-9-13 09:30 编辑在对胸腺嘧啶-二苯甲酮体系采用#p opt freq td=(triplets,nstates=4,root=1) b3lyp/6-31+g(d) em=gd3优化T1激发态，结果收敛，但出现激发能为负的情况；然后取优化结 ...

回复 Reply

yingyang 发表于 Post on 2021-11-4 19:37:12

你好，请问你出现的这种情况已经解决了吗？

zjxitcc 发表于 Post on 2021-9-13 23:59:37

yanshen 发表于 2021-9-13 15:59
请问那应该怎么计算呢？在stable=opt计算时加上限制性R，优化出闭壳层稳定的波函数？

当你从RDFT的stable=opt任务结果出来是UDFT波函数，就同时意味着不可能在RDFT下找到稳定的波函数（因为就是不稳定的）。

便宜经济的计算方案：采用spin flip TDDFT （免费软件ORCA,GAMESS；收费软件Q-Chem都有）。
稍昂贵、更精确的计算方案：CASSCF-NEVPT2

wzkchem5 发表于 Post on 2021-9-13 19:23:03

yanshen 发表于 2021-9-13 08:59
请问那应该怎么计算呢？在stable=opt计算时加上限制性R，优化出闭壳层稳定的波函数？

不是的，你这种体系（有三重态不稳定性的体系）估计所有闭壳层波函数都不稳定。stable=opt的重要作用之一就是在闭壳层波函数不稳定、稳定的波函数都是开壳层时，找到稳定的开壳层波函数。
所以，如果理论级别不变的话，正确的做法就是以开壳层单重态为参考态来做TDDFT。但是考虑到自旋污染等问题，这样做结果不准。以开壳层单重态为参考态做TDA会好一些，但是仍然无法解决自旋污染问题。建议如果计算条件允许的话，改用ORCA或Q-Chem等程序做SF-TDDFT计算，以T1为参考态，把S0当作激发态来算，这样可以很大程度上避免自旋污染问题。

yanshen 发表于 Post on 2021-9-13 15:59:53

本帖最后由 yanshen 于 2021-9-13 16:04 编辑

wzkchem5 发表于 2021-9-13 15:06
这说明体系是开壳层单重态，闭壳层单重态波函数不稳定。
注意计算时“区分α 、β电子”（准确来说叫做非 ...

请问那应该怎么计算呢？在stable=opt计算时加上限制性R，优化出闭壳层稳定的波函数？

wzkchem5 发表于 Post on 2021-9-13 15:06:56

这说明体系是开壳层单重态，闭壳层单重态波函数不稳定。
注意计算时“区分α 、β电子”（准确来说叫做非限制性计算）不意味着体系一定不是单重态，只意味着体系不是闭壳层

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register