使用cosmors溶剂化计算TFSI-阴离子，结构优化中断

vantasyoscar

看了一下输出文件，scf好像没有什么问题，

1. Cycle  10  Pass 1  IDiag  1:
   
2. RMSU=  9.84D-09    CP:  1.00D+00  1.06D+00  6.83D-01  7.14D-01  6.47D-01
   
3.                     CP:  7.92D-01  6.78D-01  8.25D-01  6.75D-01
   
4. E= -1830.22233552897     Delta-E=        0.000000000085 Rises=F Damp=F
   
5. DIIS: error= 8.54D-08 at cycle  10 NSaved=  10.
   
6. NSaved=10 IEnMin= 7 EnMin= -1830.22233554301     IErMin=10 ErrMin= 8.54D-08
   
7. ErrMax= 8.54D-08  0.00D+00 EMaxC= 1.00D-01 BMatC= 3.26D-12 BMatP= 3.11D-11

复制代码

但是结构优化时还是会提示l502错误，

1. End of input file reading file name.
   
2. Error termination via Lnk1e in /opt/software/g16_c01/g16/l502.exe at Tue Dec 24 11:39:26 2024.
   
3. Job cpu time:       0 days  0 hours  8 minutes  3.8 seconds.
   
4. Elapsed time:       0 days  0 hours  0 minutes 16.0 seconds.
   
5. File lengths (MBytes):  RWF=     87 Int=      0 D2E=      0 Chk=      9 Scr=      2
   

复制代码

不知道为什么，初始结构已经用cp2k预先优化了一遍。

Uus/pMeC6H4-/キ

先用CP2K优化再拿ASE转换成Gaussian输入文件是多此一举，直接在GaussView里面画出结构再自行编辑输入文件添加关键词多方便……

（编辑：删除了对B3P86和COSMO-RS的无知评价，还是要多学习）

vantasyoscar

Uus/pMeC6H4-/キ 发表于 2024-12-24 12:10
先用CP2K优化再拿ASE转换成Gaussian输入文件是多此一举，直接在GaussView里面画出结构再手动编辑输入文件多 ...

想要生成.cosmo文件，用于COSMO-RS软件的后续计算，我看COSMO-RS推荐的泛函是B3P86，不过之前用过B3LYP，也是加了opt关键字就报错，只计算单点就可以。

Uus/pMeC6H4-/キ

vantasyoscar 发表于 2024-12-24 12:13
想要生成.cosmo文件，用于COSMO-RS软件的后续计算，我看COSMO-RS推荐的泛函是B3P86，不过之前用过B3LYP， ...

论坛搜COSMO-RS也能看到类似的问题，可以试试先用scrf=cosmo做opt，再用scrf=cosmors做单点并得到.cosmo文件，分两步也方便确认几何优化收敛合适。