计算化学公社

标题: 镧系配合物计算激发能波函数不稳定 [打印本页]

作者
Author: 声聚变 时间: 2025-10-8 15:24
标题: 镧系配合物计算激发能波函数不稳定
各位老师，我正对镧系配合物进行激发能计算，打算计算50个激发态来绘制CTS图谱，引用了用PBE0计算单点能（为了绘制lumo,homo轨道）得到的的稳定的波函数作为初猜

以下是我的输入文件，坐标及赝势定义省略（单点能计算过5重态和7重态，7重态能量更低）
%cpu=64-127
%mem=100GB
%oldchk=OR1-stableopt1.chk
%chk=OR1-excit.chk
#p pbe1pbe/genecp em=gd3bj int=fine IOp(9/40=4) TD(full, nstates=50) stable=opt guess=read

opted

0 7

计算基态单点能部分收敛很快，迭代1轮就已收敛，后续程序显示波函数不稳定，尝试自动优化后显示报错 failed to find a lower point in StbLin

请问应该如何调整使得波函数稳定

作者
Author: zjxitcc 时间: 2025-10-8 15:46
本帖最后由 zjxitcc 于 2025-10-8 15:53 编辑

您的计算有误，首要之急不是“如何调整使得波函数稳定”。TD与stable=opt不可以写在一个任务中。要么先检验基态波函数是否稳定、不写TD；要么先不检验波函数稳定性、直接做TD计算，您只能选一种，然后展示计算是否有问题，如果没问题则无需发帖；如果有问题可以进一步展示。

以“先不检验波函数稳定性、直接做TD计算”为例，您不需要写那么多关键词，从#p这一行开始，写
#p UPBE1PBE chkbasis em=gd3bj int=fine IOp(9/40=4) TD(full,nstates=50) guess=read geom=allcheck
[空行]
就可以了。标题行、电荷、自旋多重度、坐标、自定义基组、赝势全都不需要写，写多了有时会犯小错误（例如忘记复制赝势），都从chk文件里读取比较方便。理论上int=fine也不需要写，但可能您已经做了一批这样的计算并且想保持精度一致，那还是保留。

作者
Author: 声聚变 时间: 2025-10-8 15:49

zjxitcc 发表于 2025-10-8 15:46
您的计算有误，首要之急不是“如何调整使得波函数稳定”。TD与stable=opt不可以写在一个任务中。要么先检验 ...

感谢解答，我引用的是用PBE0计算得到的已经稳定的波函数，是否无需再进行检验波函数稳定性

作者
Author: zjxitcc 时间: 2025-10-8 15:58
本帖最后由 zjxitcc 于 2025-10-8 15:59 编辑

声聚变发表于 2025-10-8 15:49
感谢解答，我引用的是用PBE0计算得到的已经稳定的波函数，是否无需再进行检验波函数稳定性

guess=read是读取chk文件中的波函数，称为读取更合适，而不是引用。
如果读取前后使用的是一样的基组和赝势，则无需再检验波函数稳定性（已经检验过了，读进来的是稳定波函数），可直接进行TD计算。

作者
Author: 声聚变 时间: 2025-10-8 15:58

zjxitcc 发表于 2025-10-8 15:46
您的计算有误，首要之急不是“如何调整使得波函数稳定”。TD与stable=opt不可以写在一个任务中。要么先检验 ...

感谢解答

作者
Author: 声聚变 时间: 2025-10-8 16:18

zjxitcc 发表于 2025-10-8 15:58
guess=read是读取chk文件中的波函数，称为读取更合适，而不是引用。
如果读取前后使用的是一样的基组和 ...

感谢解答

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