铱配合物TDDFT计算报错

王好赊

老师好！我在计算一类铱(III)配合物的激发态时遇到了问题，想请各位老师帮忙看一看。

几何优化：泛函：B3LYP-D3 ；基组：Ir用lanl08赝势基组，其它原子基组：6-31G** ，正常结束。


激发态计算：用PBE0-D3(BJ)泛函，Ir使用lanl08赝势基组，其它原子基组：6-31G**。输入代码如下：

1. %mem=800MW
   
2. %nprocshared=42
   
3. %lindaworker=ibfat01
   
4. %chk=m9-H-pbe0-td.chk
   
5. #p pbe1pbe/genecp em=GD3BJ TD(direct,nstates=5) scrf=(SMD,solvent=chloroform)
   
6. 
   
7. 1
   
8. 
   
9. 1 1
   
10. 原子坐标
    
11. 
    
12. C N H 0
    
13. 6-31G**
    
14. ****
    
15. Ir     0
    
16. S    1   1.00
    
17.       2.3500000              1.0000000
    
18. S    1   1.00
    
19.       1.5820000              1.0000000
    
20. S    1   1.00
    
21.       0.5018000              1.0000000
    
22. S    1   1.00
    
23.       0.2500000              1.0000000
    
24. S    1   1.00
    
25.       0.0598000              1.0000000
    
26. P    1   1.00
    
27.       2.7920000              1.0000000
    
28. P    1   1.00
    
29.       1.5410000              1.0000000
    
30. P    1   1.00
    
31.       0.5100000              1.0000000
    
32. P    1   1.00
    
33.       0.0980000              1.0000000
    
34. P    1   1.00
    
35.       0.0290000              1.0000000
    
36. D    1   1.00
    
37.       1.2400000              1.0000000
    
38. D    1   1.00
    
39.       0.4647000              1.0000000
    
40. D    1   1.00
    
41.       0.1529000              1.0000000
    
42. F    1   1.00
    
43.       0.9380000              1.0000000
    
44. ****
    
45. 
    
46. IR     0
    
47. IR-ECP     4     60
    
48. g potential
    
49.   5
    
50. 1    823.5880147             -0.1578014
    
51. 2    364.6613336          -1517.5270446
    
52. 2     55.7082801           -316.5306529
    
53. 2     12.0464544            -91.8880941
    
54. 2      3.5120610             -9.2241773
    
55. s-g potential
    
56.   6
    
57. 0    188.0490770              3.1578014
    
58. 1    340.4194712             26.8322577
    
59. 2    128.2373673            800.4250007
    
60. 2     33.8644961            369.4050683
    
61. 2      4.7560005            242.4171899
    
62. 2      3.9649974           -118.2173282
    
63. p-g potential
    
64.   5
    
65. 0    289.7291139              2.1578014
    
66. 1     87.4633789             61.9678610
    
67. 2     30.4363766            269.0581986
    
68. 2      4.0553412            231.1654793
    
69. 2      3.5525341           -133.6952667
    
70. d-g potential
    
71.   5
    
72. 0    136.4017106              3.1578014
    
73. 1     95.0776925             45.9349803
    
74. 2     49.2258410            359.0344668
    
75. 2     15.0874145            176.4740119
    
76. 2      4.0405764             54.5155286
    
77. f-g potential
    
78.   5
    
79. 0    127.3507908              3.9546197
    
80. 1     66.2364374             52.9773655
    
81. 2     34.4299229            274.8643383
    
82. 2     10.1995721            137.2047338
    
83. 2      2.5409702             14.8633305
    
84. 
    

复制代码

out文件报错：

1. 
   
2. Error on total polarization charges =  0.05630
   
3. SCF Done:  E(RPBE1PBE) =  -1984.06894314     A.U. after   36 cycles
   
4.             NFock= 36  Conv=0.53D-08     -V/T= 2.0493
   
5. KE= 1.890898980220D+03 PE=-1.653219335135D+04 EE= 6.668612880178D+03
   
6. SMD-CDS (non-electrostatic) energy       (kcal/mol) =     -15.60
   
7. (included in total energy above)
   
8. Leave Link  502 at Thu Apr 10 15:45:41 2025, MaxMem=   838860800 cpu:     23920.3
   
9. (Enter /home/t1000003449/g09/l801.exe)
   
10. DoSCS=F DFT=T ScalE2(SS,OS)=  1.000000  1.000000
    
11. ExpMin= 2.90D-02 ExpMax= 4.17D+03 ExpMxC= 6.27D+02 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
    
12. HarFok:  IExCor=  205 AccDes= 0.00D+00 IRadAn=         5 IDoV=-2 UseB2=F ITyADJ=14
    
13. ICtDFT= 12500011 ScaDFX=  1.000000  1.000000  1.000000  1.000000
    
14. Defaulting to unpruned grid for atomic number  77.
    
15. Defaulting to unpruned grid for atomic number  77.
    
16. Defaulting to unpruned grid for atomic number  77.
    
17. Defaulting to unpruned grid for atomic number  77.
    
18. Defaulting to unpruned grid for atomic number  77.
    
19. Defaulting to unpruned grid for atomic number  77.
    
20. Defaulting to unpruned grid for atomic number  77.
    
21. Defaulting to unpruned grid for atomic number  77.
    
22. Defaulting to unpruned grid for atomic number  77.
    
23. Defaulting to unpruned grid for atomic number  77.
    
24. Defaulting to unpruned grid for atomic number  77.
    
25. Defaulting to unpruned grid for atomic number  77.
    
26. Defaulting to unpruned grid for atomic number  77.
    
27. Defaulting to unpruned grid for atomic number  77.
    
28. Defaulting to unpruned grid for atomic number  77.
    
29. Defaulting to unpruned grid for atomic number  77.
    
30. Defaulting to unpruned grid for atomic number  77.
    
31. Defaulting to unpruned grid for atomic number  77.
    
32. Defaulting to unpruned grid for atomic number  77.
    
33. Defaulting to unpruned grid for atomic number  77.
    
34. Defaulting to unpruned grid for atomic number  77.
    
35. Defaulting to unpruned grid for atomic number  77.
    
36. Defaulting to unpruned grid for atomic number  77.
    
37. Defaulting to unpruned grid for atomic number  77.
    
38. Defaulting to unpruned grid for atomic number  77.
    
39. Defaulting to unpruned grid for atomic number  77.
    
40. Defaulting to unpruned grid for atomic number  77.
    
41. Defaulting to unpruned grid for atomic number  77.
    
42. Defaulting to unpruned grid for atomic number  77.
    
43. Defaulting to unpruned grid for atomic number  77.
    
44. Defaulting to unpruned grid for atomic number  77.
    
45. Defaulting to unpruned grid for atomic number  77.
    
46. Defaulting to unpruned grid for atomic number  77.
    
47. Defaulting to unpruned grid for atomic number  77.
    
48. Defaulting to unpruned grid for atomic number  77.
    
49. Defaulting to unpruned grid for atomic number  77.
    
50. Defaulting to unpruned grid for atomic number  77.
    
51. Defaulting to unpruned grid for atomic number  77.
    
52. Defaulting to unpruned grid for atomic number  77.
    
53. Defaulting to unpruned grid for atomic number  77.
    
54. Defaulting to unpruned grid for atomic number  77.
    
55. Defaulting to unpruned grid for atomic number  77.
    
56. Consistency failure #1 in FindCO.
    
57. Error termination via Lnk1e in /home/t1000003449/g09/l801.exe at Thu Apr 10 15:45:44 2025.
    
58. Job cpu time:       0 days  6 hours 50 minutes 37.5 seconds.
    
59. File lengths (MBytes):  RWF=    962 Int=      0 D2E=      0 Chk=     11 Scr=      2
    

复制代码

savorw

Warning!  O  atom    2 has  8 valence electrons but only   0 basis functions.
This is less than a minimal basis set!

这个警告都没看到啊。O原子都没有给基组，不出错才怪。

王好赊

savorw 发表于 2025-4-10 17:43
Warning!  O  atom    2 has  8 valence electrons but only   0 basis functions.
This is less than a  ...

我天，马虎了我，感谢感谢，谢谢老师

YoungYL

您好，我最近也在做铱配合物方面的计算，在计算几何优化过程中出现了问题，能联系您学习一下吗？

北大-陶豫

YoungYL 发表于 2025-8-5 21:13
您好，我最近也在做铱配合物方面的计算，在计算几何优化过程中出现了问题，能联系您学习一下吗？

http://sobereva.com/79

zjxitcc

YoungYL 发表于 2025-8-5 21:13
您好，我最近也在做铱配合物方面的计算，在计算几何优化过程中出现了问题，能联系您学习一下吗？

有问题直接在论坛上直接问，不要“能联系您学习一下吗”，你也不知道他了解多少、会不会讲错。即使他是专家，你也得通过导师付费向其请教（注意不是叫你付费）。量子化学计算问题一般是普遍性的，并不依赖于哪种常见元素。