计算化学公社

标题: orca计算磷光发射能为负值，读取稳定的波函数后仍为负值 [打印本页]
作者
Author: wxthhh    时间: 2024-5-10 20:54
标题: orca计算磷光发射能为负值，读取稳定的波函数后仍为负值
本帖最后由 wxthhh 于 2024-5-10 21:48 编辑

计算团簇的磷光发射能时发现为负值，然后用Gaussian关键词stable=opt读取了稳定的波函数，使用orca读取后发现发射能仍为负值，求指导
作者
Author: zjxitcc    时间: 2024-5-10 20:59
（1）你的gjf文件里没加隐式溶剂模型；（2）不需要写IOP(3/32=2)，也永远不推荐写这个。
作者
Author: wxthhh    时间: 2024-5-10 21:33
zjxitcc 发表于 2024-5-10 20:59
（1）你的gjf文件里没加隐式溶剂模型；（2）不需要写IOP(3/32=2)，也永远不推荐写这个。

抱歉老师，我那个图片上传错了，我没加cpcm的也算了，还是负的，图片重新上传了
作者
Author: zjxitcc    时间: 2024-5-10 21:35
wxthhh 发表于 2024-5-10 21:33
抱歉老师，我那个图片上传错了，我没加cpcm的也算了，还是负的，图片重新上传了

如果你的体系没什么保密性，可以直接上传gjf/inp/fch文件，方便查看。
作者
Author: wxthhh    时间: 2024-5-10 21:48
zjxitcc 发表于 2024-5-10 21:35
如果你的体系没什么保密性，可以直接上传gjf/inp/fch文件，方便查看。

老师，已经上传了，谢谢老师这个点还在帮我看文件
作者
Author: zjxitcc    时间: 2024-5-11 15:51
本帖最后由 zjxitcc 于 2024-5-11 16:18 编辑

你的轨道传得有问题。传进去后ORCA的SCF花了20圈还没收敛，自动切换成TRAH做轨道优化，不幸优化回了不稳定的UDFT解，所以激发能（又）出现了负值。正确合理的传轨道，SCF应该是立即收敛（或者能量一开始就接近收敛值，小晃几圈收敛）。这里演示一下合理的传轨道会怎样：

你已经有了ag19-t1-stable.fch文件，里面是稳定的UDFT解，运行
  1. fch2mkl ag19-t1-stable.fch -dft "B3LYP/G"
  2. orca_2mkl ag19-t1-stable_o -gbw
复制代码
此举会产生ag19-t1-stable_o.inp/mkl/gbw三个文件，DFT关键词、基组信息、辅助基全都写好了，不需要修改/增添基组。我们打开inp文件增大内存和核数，适应自己的机器，添加你的TDDFT关键词
  1. %tddft
  2. nroots 5
  3. tprint 1e-8
  4. tda true
  5. end
复制代码
注意稳定的UDFT解有自旋污染，基态和激发态都不是自旋纯态，关键词irootmult triplet没有用。提交ORCA任务
  1. orca ag19-t1-stable_o.inp >ag19-t1-stable_o.out 2>&1 &
复制代码
在输出文件中可以看到能量一开始就极其接近收敛的值，用了7圈收敛
  1. --------------
  2. SCF ITERATIONS
  3. --------------
  4. ITER       Energy         Delta-E        Max-DP      RMS-DP      [F,P]
  5.                ***  Starting incremental Fock matrix formation  ***
  6.   0 -11258.9602531861   0.000000000000 0.00130735  0.00000822  0.0002203
  7.   1 -11258.9602654113  -0.000012225246 0.00029976  0.00000205  0.0000322
  8.                                ***Turning on DIIS***
  9.   2 -11258.9602660536  -0.000000642345 0.00025841  0.00000103  0.0000151
  10.                *** Restarting incremental Fock matrix formation ***
  11.                                    *** Resetting DIIS ***
  12.   3 -11258.9602660375   0.000000016158 0.00024493  0.00000081  0.0000315
  13.   4 -11258.9602658057   0.000000231830 0.00016764  0.00000060  0.0000631
  14.   5 -11258.9602661024  -0.000000296728 0.00004803  0.00000010  0.0000068
  15.   6 -11258.9602661255  -0.000000023087 0.00002416  0.00000005  0.0000071
  16.   7 -11258.9602661110   0.000000014508 0.00002240  0.00000003  0.0000033
  17.                             ***DIIS convergence achieved***

  19.                *****************************************************
  20.                *                     SUCCESS                       *
  21.                *           SCF CONVERGED AFTER   8 CYCLES          *
  22.                *****************************************************
复制代码
下方几十行处有<S**2>=0.985360，与Gaussian结果完全一致。很快进入TDDFT计算步骤
  1. STATE  1:  E=   0.045546 au      1.239 eV     9996.3 cm**-1 <S**2> =   0.905203
  2. STATE  2:  E=   0.046583 au      1.268 eV    10223.7 cm**-1 <S**2> =   1.036779
  3. ...
复制代码
没有负激发能。fch2mkl是MOKIT的一个小程序，用于Gaussian -> ORCA传轨道。
作者
Author: wxthhh    时间: 2024-5-11 16:21
zjxitcc 发表于 2024-5-11 15:51
你的轨道传得有问题。传进去后ORCA的SCF花了20圈还没收敛，自动切换成TRAH做轨道优化，不幸优化回了不稳定 ...

好的，我去下载一下这个传轨道的小程序算一下看看，多谢老师帮助



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3