计算化学公社

标题: 求助：自旋极化单重态的结构算高精度单点能 [打印本页]

作者
Author: 15whyhq 时间: 2020-10-19 21:34
标题: 求助：自旋极化单重态的结构算高精度单点能
本帖最后由 15whyhq 于 2020-10-20 11:28 编辑

geom=read读取已优化好的自旋极化单重态的稳定的波函数，进行高精度单点能计算，计算完检测次高精度单点计算的波函数，发现波函数不稳定；
且电子能量与不读波函数直接保存优化的结构计算的电子能量一样，算溶解自由能也是这样，这是为什么？自旋极化单重态的结构这样算单点能合理吗？
优化结构的计算水平：ub3lyp/lal2tz(f)金属元素 6-31g(d,p)非金属元素 em=gd3bj
高精度单点计算水平：ub3lyp/def2tzvp em=gd3bj

作者
Author: sobereva 时间: 2020-10-19 23:36
没有geom=gusee关键词，那叫guess=read
应当说明具体有什么金属。主族和过渡金属的情况差异极大
既然高级别发现不稳定，用stable=opt再进一步得到稳定波函数

作者
Author: 15whyhq 时间: 2020-10-20 11:23
本帖最后由 15whyhq 于 2020-10-20 11:27 编辑

sobereva 发表于 2020-10-19 23:36
没有geom=gusee关键词，那叫guess=read
应当说明具体有什么金属。主族和过渡金属的情况差异极大
既然高级 ...

是过渡金属。读优化结构的稳定波函数，在高水平下做stable=opt吗？关键词写ub3lyp/def2tzvp em=gd3bj stable=opt geom=check guess=read 对吗

作者
Author: sobereva 时间: 2020-10-21 11:41

15whyhq 发表于 2020-10-20 11:23
是过渡金属。读优化结构的稳定波函数，在高水平下做stable=opt吗？关键词写ub3lyp/def2tzvp em=gd3bj sta ...

是

作者
Author: 15whyhq 时间: 2020-10-22 09:42

sobereva 发表于 2020-10-21 11:41
是

关键词写ub3lyp/def2tzvp em=gd3bj stable=opt geom=check guess(read,mix）
会比上面的节省时间吗？

作者
Author: sobereva 时间: 2020-10-24 09:57

15whyhq 发表于 2020-10-22 09:42
关键词写ub3lyp/def2tzvp em=gd3bj stable=opt geom=check guess(read,mix）
会比上面的节省时间吗？

意义不明
如果你读入的波函数已经是对称破缺波函数，还加个mix纯粹是捣乱

作者
Author: 15whyhq 时间: 2020-11-2 22:00
算溶剂化自由能的时候，基于优化的自旋极化单重态结构，算uM052X/6-31G* stable=opt geom=check guess=read 或uM052X/6-31G* scrf=(smd,solvent=ethanol) stable=opt geom=check guess=read ；
一直scf不收敛，g16算，SCF=vshift=400；SCF=conver=6；guess=huckel；SCF=QC；SCF=Fermi；单独试了这些关键词都不行，都是l502;
是什么原因啊？还能怎么处理呢？

作者
Author: wzkchem5 时间: 2020-11-3 08:10

15whyhq 发表于 2020-11-2 22:00
算溶剂化自由能的时候，基于优化的自旋极化单重态结构，算uM052X/6-31G* stable=opt geom=check guess=read ...

贴出完整的输出文件，或者最后至少5步scf迭代的情况（必须是输入文件加了#p的），我们看看

作者
Author: 15whyhq 时间: 2020-11-4 21:45
算uM052X/6-31G* stable=opt geom=check guess=read 或uM052X/6-31G* scrf=(smd,solvent=ethanol) stable=opt geom=check guess=read
稳定波函数的能量比不稳定波函数的能量低，这是因为什么？

作者
Author: zjxitcc 时间: 2020-11-4 21:49

15whyhq 发表于 2020-11-4 21:45
算uM052X/6-31G* stable=opt geom=check guess=read 或uM052X/6-31G* scrf=(smd,solvent=ethanol) stable=o ...

建议直接晒数据（坐标、能量，最好有输出文件）。仅凭这几句话信息量太少。

作者
Author: youxx189 时间: 2023-3-17 10:29

15whyhq 发表于 2020-11-2 22:00
算溶剂化自由能的时候，基于优化的自旋极化单重态结构，算uM052X/6-31G* stable=opt geom=check guess=read ...

你好请问解决了吗目前遇到同样的问题想请教一下

作者
Author: wzkchem5 时间: 2023-3-17 14:13

youxx189 发表于 2023-3-17 03:29
你好请问解决了吗目前遇到同样的问题想请教一下

楼上已经说了，这种问题必须贴出自己的数据，别人才能回答

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