24小时热门版块排行榜    

查看: 1417  |  回复: 5

xiuhui

木虫 (正式写手)

[求助] 算频率出错。

最近,我用G03计算一系列有机分子体系,当烷基链增大时,如C21,计算频率就会出现 CPHF failed to converge in LinEq2.错误。不知道该如何解决,请高手指教。
计算方法如下。
******************************************
%Nproc=4
Will use up to    4 processors via shared memory.
%nproc=4
Will use up to    4 processors via shared memory.
%chk=C11-xianan-C11.chk
%mem=10GB
-----------------------------------------------------
# b3lyp/6-311+g(d,p) opt(tight,nolinear) freq=noraman

错误信息:
  1 vectors were produced by pass991.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass992.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass993.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass994.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass995.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass996.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass997.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass998.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        13.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass999.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
PrismC:  NFx=      2048 NFxT=        14 NFxU=        14.
     1 vectors were produced by pass***.
CPHF failed to converge in LinEq2.
回复此楼

» 收录本帖的淘帖专辑推荐

高斯专辑

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gmy1990: 金币+2 2012-04-17 22:52:58
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-04-17 20:14:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

书万里

木虫 (著名写手)

历史控+原著控

引用回帖:
2楼: Originally posted by beefly at 2012-04-17 20:14:18:
参见第九项
http://people.smu.edu/wzou/qcsc/faq/gaussian.html

谢谢beefly专家,我也学了点东西,
喂人民服雾。
3楼2012-04-18 11:13:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiuhui

木虫 (正式写手)

引用回帖:
2楼: Originally posted by beefly at 2012-04-17 20:14:18:
参见第九项
http://people.smu.edu/wzou/qcsc/faq/gaussian.html

第九项是这样的

9. 问: HF/DFT解析频率计算不收敛,该如何解决?
  答:  
    在原子数很少的情况下(通常5个原子以下),可以考虑改用数值频率(Freq=Numerical或Freq=DoubleNumer)。或者修改控制CPHF/CPKS方程收敛的IOP选项,有以下几个。在单点能输出之后,找到类似于下面的输出(如果计算加上了#P开关,位于"(Enter .../g09/l1002.exe)"的下一行):

Minotr: Closed shell wavefunction.
IDoAtm=111
Direct CPHF calculation.
Differentiating once with respect to electric field.
with respect to dipole field.
Differentiating once with respect to nuclear coordinates.
Using symmetry in CPHF.
Requested convergence is 1.0D-10 RMS, and 1.0D-09 maximum.
Secondary convergence is 1.0D-12 RMS, and 1.0D-12 maximum.

最后两行定义了解CPHF/CPKS方程的四个收敛标准。
第一个标准由IOp(10/7=N)控制:10^(-N),默认N=8
第二个标准由系统设置,等于第一个标准乘10
第三、四个标准通常相等,由IOp(10/16=N)控制:10^(-N),默认N=12。用于溶剂计算。

如果是气态分子计算,把IOp(10/7)改小一点就行了。其它可能用到的还有:
IOp(10/73):可以增大CPHF方程的迭代次数,默认1000
IOp(10/18=3):CPHF方程仅迭代一次,不检查收敛。可以用来检查定性结果,通常会带来几个到几十个波数的振动频率误差,不能用于正式计算。

另外还需要注意,IOP不能在多步任务中传递,因此如果使用上面的IOP选项,结构优化和频率计算要分开算。


我算的大分子体系,针对大分子有什么解决方案吗?
4楼2012-04-18 17:24:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiuhui

木虫 (正式写手)

引用回帖:
2楼: Originally posted by beefly at 2012-04-17 20:14:18:
参见第九项
http://people.smu.edu/wzou/qcsc/faq/gaussian.html

多谢beefly专家,但是貌似第九项对于小分子,我算的分子比较大,原子数>50的话,也能这么处理吗?多谢啦!
5楼2012-04-18 22:01:18
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiuhui: 金币+20, ★★★很有帮助 2012-04-22 22:20:43
频率计算加上开关iop(10/33=1),可以判断属于以下情况中的哪一种。

如果CPHF有收敛趋势,可以降低收敛阈值IOp(10/7),或者增加迭代次数IOp(10/73),也可以同时设置这两个参数。

如果CPHF没有收敛趋势(通常迭代是发散的)。少数情况下是数值误差造成的,用同样的结构重新做频率计算,运气好的话也许能避免这个问题。一般的解决办法是换理论方法(泛函和基组)重新优化结构和计算频率。如果原子数在20以下,也可以算半数值频率(Freq=Numerical),但这通常比解析频率耗时。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2012-04-19 00:44:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiuhui 的主题更新
信息提示
请填处理意见