求助：用ORCA计算SOC时报错

hebrewsnabla

emission 发表于 2024-7-26 19:05
老师，我重新用高斯使用5D 7F优化波函数，然后用MOKIT获得了gbw文件，和计算SOC的inp文件放在同一个文件 ...

你的输出看起来orca没有读取轨道

emission

hebrewsnabla 发表于 2024-7-26 19:39
你的输出看起来orca没有读取轨道

老师，我本来是把inp和同名gbw文件放在一起运行的，没有加多余的关键词。后面我翻了翻论坛的贴子，试着加上moread关键词和bgw文件的路径，但是看输出文件还是没有读取，请问老师知道这是什么问题吗？
我都输入文件是
! CAM-B3LYP D3 def2-TZVP def2/J def2-TZVP/C RIJCOSX miniprint tightSCF
!moread
%moinp "/home/tang/data_disk1/WSY/orca/inp/DBTCO-S1.gbw"
%maxcore 3000
%pal nprocs  20 end
%tddft
nroots 10
dosoc true
tda false
printlevel 3
end
* xyz 0 1


同时附上了输出文件

hebrewsnabla

emission 发表于 2024-7-27 11:59
老师，我本来是把inp和同名gbw文件放在一起运行的，没有加多余的关键词。后面我翻了翻论坛的贴子，试着加 ...

可以试试去掉miniprint

emission

hebrewsnabla 发表于 2024-7-27 15:17
可以试试去掉miniprint

谢谢老师，我按您说的去掉miniprint，这次输出结果里显示读取gbw文件了！
可是还是报错了，和之前的报错一样Negative eigenvalue found in transformed Hamiltonian，请问老师知道这是什么原因吗？应该怎么去调整呀？
附上了输出文件

hebrewsnabla

emission 发表于 2024-7-29 10:14
谢谢老师，我按您说的去掉miniprint，这次输出结果里显示读取gbw文件了！
可是还是报错了， ...

用fch传给orca的话，高斯任务的基组要和orca任务的基组相同。使用fch2mkl的话，要使用fch2mkl创建的inp文件（可以在此基础上添加td 辅助基的关键词），里面已经写好了基组可以保证基组是相同的。

wzkchem5

emission 发表于 2024-7-29 03:14
谢谢老师，我按您说的去掉miniprint，这次输出结果里显示读取gbw文件了！
可是还是报错了， ...

老生常谈的问题了，说明基态波函数不稳定。
最好的解决方法是改用NEVPT2算SOC，次好是检查泛函是否不适合、如果不适合的话更换泛函，再次好是打开TDA

emission

hebrewsnabla 发表于 2024-7-29 15:07
用fch传给orca的话，高斯任务的基组要和orca任务的基组相同。使用fch2mkl的话，要使用fch2mkl创建的inp文 ...

谢谢老师，我确实没用fch2mkl的inp文件，我试一下

emission

wzkchem5 发表于 2024-7-29 16:12
老生常谈的问题了，说明基态波函数不稳定。
最好的解决方法是改用NEVPT2算SOC，次好是检查泛函是否不适 ...

谢谢老师，不太会用NEVPT2，打算后面用不同的泛函试一下

zjxitcc

fch2mkl小程序的使用十分简单，先写一个Gaussian输入文件

1. %chk=DBTCO.chk
   
2. %mem=180GB
   
3. %nprocshared=64
   
4. #p UCAM-B3LYP/def2TZVP nosymm int=nobasistransform guess=mix stable=opt
   
5. 
   
6. title
   
7. 
   
8. 0 1
   
9. C     -1.49891114    0.07026687   -0.50088109
   
10. ...
    
11. H     -9.56460524    2.69839127    0.42449638
    

复制代码

提交Gaussian计算，结果为E(UCAM-B3LYP) = -5657.24889771, <S**2>= 0.8324，基态为开壳层单重态，有自旋污染。接着运行

1. formchk DBTCO.chk DBTCO.fch
   
2. fch2mkl DBTCO.fch -dft 'CAM-B3LYP D3'
   
3. orca_2mkl DBTCO_o -gbw
   

复制代码

此举会产生DBTCO_o.inp和DBTCO_o.gbw文件，前者含体系坐标，基组数据和基态UCAM-B3LYP关键词，后者含有轨道信息。我们打开DBTCO_o.inp文件加点激发态关键词，前13行内容如下

1. %pal nprocs 16 end
   
2. %maxcore 10800
   
3. ! UKS CAM-B3LYP D3 def2/J def2-TZVP/C TightSCF defgrid3
   
4. %tddft
   
5. nroots 5
   
6. tda false
   
7. printlevel 3
   
8. end
   
9. %scf
   
10. Thresh 1e-12
    
11. Tcut 1e-14
    
12. CNVDamp False
    
13. end

复制代码

这里不用写基组def2-TZVP，因为输入文件底下的基组数据已经体现了。提交ORCA计算，可以观察到SCF 5圈收敛，能量变化微小，立即进入TDDFT计算步骤

1. -------------------------------------------------------------------------------------------
   
2.                                       ORCA LEAN-SCF
   
3.                               memory conserving SCF solver
   
4. -------------------------------------------------------------------------------------------
   
5. 
   
6. ----------------------------------------D-I-I-S--------------------------------------------
   
7. Iteration    Energy (Eh)           Delta-E    RMSDP     MaxDP     DIISErr   Damp  Time(sec)
   
8. -------------------------------------------------------------------------------------------
   
9.                ***  Starting incremental Fock matrix formation  ***
   
10.                               *** Initializing SOSCF ***
    
11. ---------------------------------------S-O-S-C-F--------------------------------------
    
12. Iteration    Energy (Eh)           Delta-E    RMSDP     MaxDP     MaxGrad    Time(sec)
    
13. --------------------------------------------------------------------------------------
    
14.     1    -5657.2503433988385950     0.00e+00  5.96e-05  3.60e-03  1.13e-05    82.4
    
15.                *** Restarting incremental Fock matrix formation ***
    
16.     2    -5657.2503533339295245    -9.94e-06  1.24e-05  8.87e-04  2.74e-05    79.8
    
17.     3    -5657.2503537240827427    -3.90e-07  1.85e-06  1.61e-04  5.62e-06    51.3
    
18.     4    -5657.2503537342090567    -1.01e-08  4.67e-07  3.19e-05  7.71e-06    42.5
    
19.     5    -5657.2503537442062225    -1.00e-08  3.87e-07  2.59e-05  3.19e-06    37.4
    
20.                  **** Energy Check signals convergence ****
    
21. 
    
22.                *****************************************************
    
23.                *                     SUCCESS                       *
    
24.                *           SCF CONVERGED AFTER   5 CYCLES          *
    
25.                *****************************************************

复制代码

但这只能解决负激发能的问题，让你可以合理计算此级别下的激发态能量。由于各个电子态有自旋污染，无法计算旋轨耦合，所以我没写doSOC；这种情况下最理想的方式是在CASSCF-NEVPT2下计算SOC，但此体系非常大，而且有很多孤对电子和pi键，这会使活性空间比较大，同时基组又不能低于TZVP，我认为NEVPT2计算多个电子态且计算SOC不是很现实。比较经济的方案是用MRSF-TDDFT计算SOC；再次一点的方案是换个泛函使基态波函数稳定，这样就用不到UKS了。