24小时热门版块排行榜    

查看: 1171  |  回复: 2

中科院苦行僧

木虫 (小有名气)

[求助] 计算频率分析 Tx not orthogonal to T.

用CIS优化构型正确 计算频率时失败 一直报错:

Inverted reduced A of dimension   428 with in-core refinement.
Max Tx.T =      0.000011726838 for Degree =    10
Tx not orthogonal to T.
Error termination via Lnk1e in /public1/software/gaussian09/g09/l1014.exe at Sun Sep 27 13:21:31 2015.
Job cpu time:  0 days  0 hours 55 minutes 19.7 seconds.
File lengths (MBytes):  RWF=   2202 Int=      0 D2E=      0 Chk=      5 Scr=      2

输入文件:
%nprocshared=4
Will use up to    4 processors via shared memory.
%mem=2000MW
%chk=S1-CIS-3Cl.chk
-------------------------------------------
#p freq=noraman rcis=(root=1)/6-311++g(d,p)
-------------------------------------------
1/10=4,30=1,38=1/1,3;
2/12=2,17=6,18=5,40=1/2;
3/5=4,6=6,7=1111,11=1,16=1,25=1,30=1,71=2,116=1/1,2,8,3;
4//1;
5/5=2,38=5,98=1/2;
8/6=1,10=2,107=1/1,4;
9/8=1,15=3,42=3/14;
10/5=4,20=4/2;
11/6=1,8=1,12=-1,16=1,42=3/1,2,10;
10/6=2,21=1,28=1/2,14;
6/7=2,8=2,9=2,10=2,18=1/1;
7/8=1,10=1,12=6,13=5,14=1,25=1,44=2/1,2,3,16;
1/10=4,30=1/3;
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : wrong-S1-CIS-3Cl-trans-freq.log
  • 2015-10-08 13:04:06, 52.67 K

» 猜你喜欢

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

中科院苦行僧

木虫 (小有名气)

这个问题有大神碰到过吗,跪求啊
beayardstickogquality
2楼2015-10-09 10:34:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

Sobereva在帖子《量子化学计算中帮助几何优化收敛的常用方法》中说:

RFO、GDIIS、GEDIIS这三种优化方法在走步时都需要Hessian矩阵(力常数矩阵),然而计算Hessian矩阵是很耗时的。因此,在默认情况下,Gaussian会通过价力场的方法近似估计出初始的Hessian矩阵,在每一步优化中只精确计算梯度,利用梯度对原先Hessian矩阵进行修正不断得到新的Hessian矩阵。所以,在默认情况下,优化从头到尾使用的Hessian矩阵都只是近似的(这也是为什么在优化出来的结构上做Freq往往会显示还没收敛,因为Freq用的是精确的Hessian矩阵)。当Hessian矩阵离精确值偏差较大,就会造成收敛缓慢,或者始终不收敛。为解决这个问题,可以用opt=calcfc,这会在优化的第一步使用精确计算的Hessian矩阵,但是仍可能后续优化过程中Hessian矩阵逐渐变得越来越不精确而依然收敛失败。opt=calcall则不仅在第一步,在后续的每一步中也都精确计算Hessian矩阵,这使得很多优化失败的情况都能得到解决,优化所需步数通常也会减少很多,而且能够保证最终优化结果准确(因为最终判断是否收敛时是基于精确的Hessian矩阵所得结果),但代价是每一步计算量会很大。

原帖:量子化学计算中帮助几何优化收敛的常用方法
http://www.sobereva.com/164
3楼2015-10-19 14:24:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员PzX7SI 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 03:34 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 03:30 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 03:27 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 4/200 2026-08-04 03:25 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 02:10 by Ocx3OKOzGU4q
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 8/400 2026-08-03 16:02 by yanghuping
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见