计算化学公社

标题: 结构优化和振动分析问题求助 [打印本页]

作者
Author: avazaza 时间: 2019-5-22 20:53
标题: 结构优化和振动分析问题求助
请教各位老师，我在对一个由71个原子组成的的有机刚性大分子进行结构优化和振动分析，根据sob老师帖子中的建议选用PBE0泛函进行计算，开始用1.# opt freq pbepbe/6-31g(d,p) geom=connectivity iop(3/76=1000002500,3/77=0750007500)，发现结构优化不收敛；
然后2.在1的基础上加了int=ultrafine，但是优化收敛了但是结果中振动出现虚频；
3.于是用了更严的收敛限并定义水作溶剂（前两次忘记考虑溶剂效应了），#P opt=tight freq pbepbe/6-31g(d,p) scrf=(solvent=water) geom=connectivity iop(3/76=1000002500,3/77=0750007500) int=ultrafine，然后输出文件出现了错误，好像连结构优化也没有收敛。结果如下（附件太大无法上传）：
Minotr:  Closed shell wavefunction.
      IDoAtm=11111111111111111111111111111111111111111111111111
      IDoAtm=111111111111111111111
NEqPCM:  Using equilibrium solvation (IEInf=0, Eps=  78.3553, EpsInf= 1.7778)
      Direct CPHF calculation.
      Differentiating once with respect to electric field.
            with respect to dipole field.
      Differentiating once with respect to nuclear coordinates.
      Requested convergence is 1.0D-08 RMS, and 1.0D-07 maximum.
      Secondary convergence is 1.0D-12 RMS, and 1.0D-12 maximum.
      NewPWx=T KeepS1=F KeepF1=F KeepIn=T MapXYZ=F SortEE=F KeepMc=T.
      MDV=    33553718 using IRadAn=    1.
Generate precomputed XC quadrature information.
      Solving linear equations simultaneously, MaxMat=    80.
AlAXAO:  NMat= 216 NPMax= 16 NPMax1= 0 MaxMat= 80.
33526473 words are not enough for AlAXAO.
Error termination via Lnk1e in /gdata/lun7/software/g09/l1002.exe at Tue May 21 23:45:03 2019.

想要知道我选用的泛函和基组是否适合我的有机体系，请问各位老师有没有其他的建议，我是否应该在前面的基础上继续尝试一些解决收敛问题的方法，还有第三次计算出现的问题该如何解决呢。

自学一个月的计算小白恳请各位老师指点迷津，万分感谢！

作者
Author: sobereva 时间: 2019-5-23 00:35
你用PBE0就写PBE1PBE就完了，干嘛写一堆IOp

你传的文件里报错提示是SCF不收敛，跟几何优化不收敛完全是两码事

你贴的报错提示是振动分析步骤，都做到这了不可能优化没收敛。报错8成是给的内存上限太小，用%mem适当加大

此体系给氢加极化根本没必要，去了还能节约时间

上传较大文本型文件前记得先压缩

作者
Author: avazaza 时间: 2019-5-23 19:57

sobereva 发表于 2019-5-23 00:35
你用PBE0就写PBE1PBE就完了，干嘛写一堆IOp

你传的文件里报错提示是SCF不收敛，跟几何优化不收敛完全是 ...

非常感谢sob老师,我按照您的建议算了一遍结果是Normal termination 并且没有虚频了。
但是为什么查看电荷分布时无法相关数据呢

作者
Author: sobereva 时间: 2019-5-24 04:16

avazaza 发表于 2019-5-23 19:57
非常感谢sob老师,我按照您的建议算了一遍结果是Normal termination 并且没有虚频了。
但是为什么查看电 ...

gview的bug。用gv打开fch文件，或者改用gv6，或者自行看输出文件

作者
Author: avazaza 时间: 2019-5-24 08:52

sobereva 发表于 2019-5-24 04:16
gview的bug。用gv打开fch文件，或者改用gv6，或者自行看输出文件

明白了，非常感谢老师

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