主客体分子结构优化，出现L502错误无法解决

CCm

是大环主体分子和客体分子结合的基态结构优化不成功，命令行大概是（几次加加减减命令行都是同样错误，记不清下列输出文件具体对应的命令行，但总体差不多）#T b3lyp/STO-3G opt(maxcyc=1024) scf(maxcyc=1024) em=gd3bj，看了答疑帖子但是一头雾水，SCF=NoVarAcc，SCF=NoIncFock，IOp(5/37=N这几个我都加过，依旧还是L502这个问题，输入文件见附件，输出文件因提示无法上传那么大的，所以将末尾贴在下方，请老师解惑 input.com (16.19 KB, 下载次数 Times of downloads: 10)

2024-3-26 00:04 上传 Uploaded

点击下载
Click to download
输入文件

1. <blockquote> Rotational constants (GHZ):           0.0107539           0.0040743           0.0038394

复制代码

CCm

1. Rotational constants (GHZ):           0.0107539           0.0040743           0.0038394
   
2.    827 basis functions,  2481 primitive gaussians,   827 cartesian basis functions
   
3.    519 alpha electrons      519 beta electrons
   
4.        nuclear repulsion energy     35130.2351277791 Hartrees.
   
5. NAtoms=  279 NActive=  279 NUniq=  279 SFac= 1.00D+00 NAtFMM=   60 NAOKFM=T Big=T
   
6. Nuclear repulsion after empirical dispersion term =    35129.5949089217 Hartrees.
   
7. Initial guess from the checkpoint file:  "me-beta-CD+cfp-3.chk"
   
8. B after Tr=     0.000000    0.000000   -0.000000
   
9.          Rot=    0.999740    0.022483   -0.003756    0.000311 Ang=   2.61 deg.
   
10. ExpMin= 1.69D-01 ExpMax= 1.31D+02 ExpMxC= 1.31D+02 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
    
11. Harris functional with IExCor=  402 and IRadAn=       5 diagonalized for initial guess.
    
12. HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         5 IDoV= 1 UseB2=F ITyADJ=14
    
13. ICtDFT=  3500011 ScaDFX=  1.000000  1.000000  1.000000  1.000000
    
14. FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=        2001
    
15.          NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
    
16.          wScrn=  0.000000 ICntrl=       500 IOpCl=  0 I1Cent=   200000004 NGrid=           0
    
17.          NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
    
18. Petite list used in FoFCou.
    
19. Restarting incremental Fock formation.
    
20. Restarting incremental Fock formation.
    
21. Restarting incremental Fock formation.
    
22. Restarting incremental Fock formation.
    
23. Restarting incremental Fock formation.
    
24. Restarting incremental Fock formation.
    
25. >>>>>>>>>> Convergence criterion not met.
    
26. SCF Done:  E(RB3LYP) =  -6605.54032045     A.U. after  129 cycles
    
27.             NFock=128  Conv=0.33D-05     -V/T= 2.0175
    
28. Convergence failure -- run terminated.
    
29. Error termination via Lnk1e in /home/sun/g16/l502.exe at Sun Mar 17 15:40:15 2024.
    
30. Job cpu time:       0 days 22 hours  2 minutes 23.9 seconds.
    
31. Elapsed time:       0 days  1 hours 51 minutes  5.1 seconds.
    
32. File lengths (MBytes):  RWF=    514 Int=      0 D2E=      0 Chk=    102 Scr=      1

复制代码

sobereva

CCm 发表于 2024-3-26 00:07

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

sobereva

缺乏最基本的常识，完全在做毫无意义的计算

好好把这些看了再做计算
谈谈量子化学中基组的选择（你拿STO-3G做计算没有丝毫意义，还不如像样的半经验级别方法准确）
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
解决SCF不收敛问题的方法（一个字一个字阅读！此文你要么没读过，要么看了但完全没过脑子）
http://sobereva.com/61
别乱加IOp
置顶的社员必读贴明确说了，较大文本文件必须压缩后再上传，认真阅读社员必读贴

这么昂贵的任务，写#T干嘛，当然得用#P监控

这种体系先用Gaussian自带的PM6D3或者xtb程序支持的GFN-xTB等半经验级别预优化，再用DFT结合个像样的基组搞。


另外，强烈建议仔细阅读
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

zhuwu

请问一下你的主客体分子是怎么放到一起的？

CCm

sobereva 发表于 2024-3-26 00:19
缺乏最基本的常识，完全在做毫无意义的计算

好好把这些看了再做计算

好的老师，非常抱歉，我会阅读您给的文章