计算化学公社

标题: 寻找纤维二糖降解反应过渡态求助 [打印本页]

作者
Author: yolo-blossom 时间: 2025-8-30 08:40
标题: 寻找纤维二糖降解反应过渡态求助
各位老师好，目前我在计算一个纤维二糖降解反应，现在在计算一个过渡态，当前使用的关键词是#p B3LYP/6-31G* SCRF=(SMD, Solvent=Water) opt=(ts,calcfc,noeigen) freq
计算一段时间以后出现的问题如下：

Error: segmentation violation
rax 0000000000000000, rbx ffffffffffffffff, rcx ffffffffffffffff
rdx 000000000005deee, rsp 00007ffec8c64db8, rbp 00007ffec8c65330
rsi 000000000000000b, rdi 000000000005deee, r8  0000000000000020
r9  0000000000000401, r10 00007ffec8c641a0, r11 0000000000000206
r12 00007ffec8c65378, r13 0000000000000000, r14 0000000000000000
r15 00000000000003e6
  /lib64/libpthread.so.0(+0xf5d0) [0x2aef4365a5d0]
  /lib64/libc.so.6(kill+0x7) [0x2aef43b9f4d7]
  /data/apps/g16//g16/l502.exe() [0x4afeb9]
  /data/apps/g16//g16/l502.exe() [0x4d1e0a]
  /data/apps/g16//g16/l502.exe() [0x54e80c]
  /data/apps/g16//g16/l502.exe() [0x52d7bf]
  /data/apps/g16//g16/l502.exe() [0x4529db]
  /data/apps/g16//g16/l502.exe() [0x450c2b]
  /data/apps/g16//g16/l502.exe() [0x450b34]
  /lib64/libc.so.6(__libc_start_main+0xf5) [0x2aef43b8b3d5]
  /data/apps/g16//g16/l502.exe() [0x450a29]

CoulSu:  requested number of processors reduced to:  35 ShMem 1 Linda.
CoulSu:  requested number of processors reduced to:  32 ShMem 1 Linda.
Error on total polarization charges =  0.06066
SCF Done:  E(RB3LYP) =  -2361.37071007    A.U. after  129 cycles
         NFock=128  Conv=0.10D-03    -V/T= 2.0088
KE= 2.340884618138D+03 PE=-1.642248426308D+04 EE= 6.272118328963D+03
SMD-CDS (non-electrostatic) energy    (kcal/mol) =    11.64
(included in total energy above)
Convergence failure -- run terminated.
Error termination via Lnk1e in /data/apps/g16//g16/l502.exe at Fri Aug 29 17:01:20 2025.
Job cpu time:    2 days  9 hours 20 minutes 36.2 seconds.
Elapsed time:    0 days  1 hours 22 minutes  3.4 seconds.
File lengths (MBytes):  RWF= 5008 Int=    0 D2E=    0 Chk=    59 Scr=    1
以上这种情况的出现，请问我需要如何优化我的初始体系或者优化我所使用的关键词呢，还请各位老师多多指教

作者
Author: 〇66 时间: 2025-8-30 09:36
可以上传一下压缩的输出文件。目前我目测大概率是因为SCF不收敛导致的，解决办法可以看卢老师的这个帖子http://sobereva.com/61。后续如果还有报错原因可以按照这个帖子进行排查http://bbs.keinsci.com/thread-4829-1-1.html

作者
Author: 2429783189 时间: 2025-8-30 22:22
不建议用SMD溶剂，难收敛，还容易出额外虚频，建议改IEFPCM。

作者
Author: yolo-blossom 时间: 2025-9-1 08:00

〇66 发表于 2025-8-30 09:36
可以上传一下压缩的输出文件。目前我目测大概率是因为SCF不收敛导致的，解决办法可以看卢老师的这个帖子htt ...

好的谢谢您的指导

作者
Author: yolo-blossom 时间: 2025-9-1 08:01

2429783189 发表于 2025-8-30 22:22
不建议用SMD溶剂，难收敛，还容易出额外虚频，建议改IEFPCM。

好的我再更换为IEFPCM试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)