Gaussian计算自由基反应时过渡态自旋多重度设置疑问

JohnYao

zjxitcc 发表于 2022-9-19 16:33
如此循环下去，直至稳定为止。若过程中出现这些现象：波函数频繁不稳定、激发能为负、在多次稳定+优化的 ...

老师您好，浏览众多帖子之后，关于流程还是有些细节上的小困惑，抱歉打搅你了，想请教下老师就是在做自旋计划单重态分子的计算时，
第一步，通过#p UB3LYP 6-31G(d) em=gd3bj guess=mix nosymm stable=opt 得到稳定波函数以后，
第二步，通过#p UB3LYP 6-31G(d) em=gd3bj opt(ts,calcfc,noeigen) freq guess(read) 得到优化后过渡态结构。且验证第二步结果波函数依旧稳定
第三步，通过#p M062X def2tzvp em=gd3 guess=mix 对第二步优化完的结构计算高精度单点能后，验证发现第三步波函数不稳定。
问题一，此时重新优化波函数的话是该读取第三步的结果 使用 #p M062X def2tzvp em=gd3 guess=read nosymm stable=opt关键词，从而重新获得稳定波函数吗？
问题二，在得到M062X def2tzvp高精度下的稳定波函数后，是否需要使用 #M062X def2tzvp em=gd3 opt(ts,calcfc,noeigen) freq guess=read 在高精度下重新对结构进行优化；或者是直接可以从上一步高精度下得到的稳定波函数结果中读取单点能

zjxitcc

第二步和第三步建议也加上nosymm。
问题一：可以，这是一种做法。另一种做法是在第三步时就应该使用guess=read读取第二步的轨道，说不定结果就是稳定波函数，就不用后面这么一着了。
问题二：既然叫高精度单点，往往意味着这个计算级别你只能算得动单点，做不动结构优化+频率分析，所以当然是不用做结构优化或频率分析了。

JohnYao

zjxitcc 发表于 2023-12-2 20:49
第二步和第三步建议也加上nosymm。
问题一：可以，这是一种做法。另一种做法是在第三步时就应该使用guess= ...

收到 感谢老师的耐心解答 万分感谢您的帮助

ljh123

sobereva 发表于 2022-9-18 04:43
当前计算级别太差。本来文献里逆向势垒就几乎可忽略不计，最起码也应当用M06-2X/6-31G**来算。而且在这个势 ...

sob老师，我的反应物是Fe(IV)和苯酚的络合物，总自旋多重度为5，电荷为0。产物是苯酚酚羟基单电子转移给Fe(IV)生成Fe(III)和苯酚阳离子自由基，总自旋多重度为7，电荷为0。所以我想请问下，我可以用AIMD模拟出这个电子转移反应吗？（比如以反应物开始，自旋多重度设为5，这个模拟过程会出现自旋多重度变成7吗？）

sobereva

ljh123 发表于 2024-11-25 19:18
sob老师，我的反应物是Fe(IV)和苯酚的络合物，总自旋多重度为5，电荷为0。产物是苯酚酚羟基单电子转移给F ...

Gaussian、ORCA做AIMD，自旋多重度没法中途变
CP2K做AIMD时靠RELAX_MULTIPLICITY可以允许自旋多重度自发发生改变