计算化学公社

标题: 求助singlet diradical结构优化与过渡态优化关键词 [打印本页]

作者
Author: chemzhh 时间: 2019-3-12 15:04
标题: 求助singlet diradical结构优化与过渡态优化关键词
请问对于singlet diradical的结构优化以及过渡态优化，写如下关键词是否正确，其中guess=mix是否要作为关键词写入，谢谢
以纯有机体系为例
结构优化
# opt freq 5d UM062X/6-311+G(d,p) guess=mix

0 1

过渡态优化
# opt=(calcfc,ts,noeigen) freq 5d UM062X/6-311+G(d,p) guess=mix

0 1

作者
Author: zjxitcc 时间: 2019-3-12 15:44
本帖最后由 zjxitcc 于 2019-3-12 15:52 编辑

除非特殊目的5d不需要写，自己不懂啥意思也不要写；结构优化用6-311G(d), 6-311G(d,p), def2-TZVP等等足矣，不需要加弥散。建议：
(1) 结构优化
#p opt=calcfc UM062X/def2TZVP guess=mix nosymm
(2) 检验波函数稳定性，判断是否双自由基
#p UM062X/def2TZVP guess=read nosymm stable=opt
(3) 若波函数稳定，做频率分析。若不稳定，上一步会自动优化至稳定波函数，然后你要手动回到（1）。
#p UM062X/def2TZVP guess=read nosymm freq

加#p输出内容多些，万一出错，方便查错。calcfc有利于加快结构优化收敛，体系不大时建议使用。
一些双自由基特征不是很强的分子，用UM062X，在（2）时往往会发现波函数是稳定的，且<S**2>=0，即是闭壳层，观察不到双自由基特征，说明实际上用RM062X足矣。

最后，如果你算出闭壳层，但自己还是觉得是singlet diradical不死心，可以试试更换初始波函数，看sob老师的博文
《谈谈片段组合波函数与自旋极化单重态》 http://sobereva.com/82

作者
Author: chemzhh 时间: 2019-3-12 15:52
谢谢你的详尽回复。
“若不稳定，上一步会自动优化至稳定波函数，然后你要手动回到（1）。”这句话的意思是如果波函数不稳定就需要用（2）的结构进行（1）的结构优化？
另外singlet diradical过渡态的优化是否与结构优化一致，有无特殊的地方？

作者
Author: zjxitcc 时间: 2019-3-12 15:56

chemzhh 发表于 2019-3-12 15:52
谢谢你的详尽回复。
“若不稳定，上一步会自动优化至稳定波函数，然后你要手动回到（1）。”这句话的意思 ...

是的。因为结构虽然是（2）的，但是波函数变了，力也变了，因此未必是极小点，需要重新做结构优化。

TS优化没有什么特殊的地方。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)