24小时热门版块排行榜    

查看: 1657  |  回复: 18

xiemin8761

铁杆木虫 (著名写手)

[求助] 2D-scan求助

根据一篇文章,我尝试用2D-scan扫描一个分子,计算的都失败了不知什么原因,请各位高手指教。
我的输入文件是:
# opt(maxcyc=100,addredundant) upbe1pbe/6-31+g* scf=(qc,nosymm) geom=connectivity

Title Card Required

0 1
C              
H                  1            B1
H                  1            B2    2            A1
H                  1            B3    2            A2    3            D1    0
C                  1            B4    2            A3    3            D2    0
C                  5            B5    1            A4    2            D3    0
H                  6            B6    5            A5    1            D4    0
H                  6            B7    5            A6    1            D5    0
H                  6            B8    5            A7    1            D6    0
O                  5            B9    1            A8    6            D7    0
O                 10           B10    5            A9    1            D8    0
H                 11           B11   10           A10    5            D9    0
H                 11           B12   10           A11    5           D10    0

   B1             1.09041428
   B2             1.09834559
   B3             1.09588262
   B4             1.50403959
   B5             1.48882226
   B6             1.08944121
   B7             1.10029571
   B8             1.09377481
   B9             1.22943504
   B10            2.25904576
   B11            0.97274001
   B12            0.97271156
   A1           109.92873208
   A2           111.12632155
   A3           112.20520707
   A4           120.67941448
   A5           111.38556342
   A6           106.61508768
   A7           109.59825862
   A8           115.77988206
   A9           112.26541245
   A10           97.56795174
   A11           99.56438207
   D1          -119.61623552
   D2           118.83803220
   D3           169.57863237
   D4           161.78293646
   D5           -79.13473791
   D6            36.72209100
   D7          -179.31990089
   D8          -177.29025062
   D9           129.69444934
   D10         -123.23665237

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 10 2.0
6 7 1.0 8 1.0 9 1.0
7
8
9
10
11 12 1.0 13 1.0
12
13

7 11 2.33673 s 10 -0.153
7 10 2.5872 s 10 -0.178

输出文件显示的错误为:

Iteration100 RMS(Cart)=  0.00002839 RMS(Int)=  0.03650170
New curvilinear step not converged.
RedCar failed.
Error termination via Lnk1e in /usr/local/intel/g03/l103.exe at Thu Nov  8 21:54:38 2012.
Job cpu time:  0 days  3 hours  1 minutes 54.0 seconds.
File lengths (MBytes):  RWF=     29 Int=      0 D2E=      0 Chk=     19 Scr=      1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

还有时候出错显示
。。。
Iteration 96 RMS(Cart)=  0.00018544 RMS(Int)=  0.47065122
Iteration 97 RMS(Cart)=  0.00016406 RMS(Int)=  0.47055968
Iteration 98 RMS(Cart)=  0.00014510 RMS(Int)=  0.47047863
Iteration 99 RMS(Cart)=  0.00012829 RMS(Int)=  0.47040689
Iteration100 RMS(Cart)=  0.00011339 RMS(Int)=  0.47034342
New curvilinear step not converged.
RedCar failed in NewRed.
Error termination via Lnk1e in /usr/local/intel/g03/l101.exe at Mon Nov 12 12:59:51 2012.
Job cpu time:  0 days  0 hours  0 minutes  3.3 seconds.
File lengths (MBytes):  RWF=      8 Int=      0 D2E=      0 Chk=      1 Scr=      1
2楼2012-11-12 12:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你先看一下这个是否为优化好的体系
3楼2012-11-12 19:43:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by lishijunzong at 2012-11-12 19:43:50
你先看一下这个是否为优化好的体系

已经是优化好了的,而且我用高斯09计算他竟然可以计算完成,但是不scan奇怪了。
4楼2012-11-13 08:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖


xiemin8761: 金币+1 2012-11-14 12:18:49
你重新提交下,103错误时gaussian的一个bug
5楼2012-11-13 11:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

引用回帖:
5楼: Originally posted by lishijunzong at 2012-11-13 11:12:59
你重新提交下,103错误时gaussian的一个bug

不明白,重新提交啥意思?
6楼2012-11-13 12:12:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖


xiemin8761: 金币+1, 有帮助 2012-11-14 12:19:50
引用回帖:
6楼: Originally posted by xiemin8761 at 2012-11-13 12:12:07
不明白,重新提交啥意思?...

之前的作业是系统错误,你从新提交作业就行
7楼2012-11-13 18:40:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xy00083789

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
xiemin8761: 金币+1, 有帮助 2012-11-14 12:19:07
xiemin8761: 金币+3, ★★★很有帮助, 我再看看能不能算。谢谢 2012-11-14 16:33:22
楼主,原因在于扫描坐标与内坐标不对应。
是这样:你在设置scan的时候
7 11 2.33673 s 10 -0.153
7 10 2.5872 s 10 -0.178
是设置的H(7) 和 两个O(10,11)的键长扫描。但是你的内坐标定义里面 两个O分别和 5,10成键,所以无法识别被扫描的坐标,造成错误。
解决方法是,在gjf的atomlist里面先构造好你需要的键长,之后再在modredundant里面,定义对该键长进行扫描。
欢迎交流zhongyue.yang@yahoo.cn
8楼2012-11-13 18:58:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

引用回帖:
8楼: Originally posted by xy00083789 at 2012-11-13 18:58:56
楼主,原因在于扫描坐标与内坐标不对应。
是这样:你在设置scan的时候
7 11 2.33673 s 10 -0.153
7 10 2.5872 s 10 -0.178
是设置的H(7) 和 两个O(10,11)的键长扫描。但是你的内坐标定义里面 两个O分别和 5 ...

能不能具体一点,怎么构造新的键长呢?奇怪的是有一次我计算2D-scan竟然成功了,后来做的2D-scan都成功不了,不知道啥原因。谢谢您的应助
9楼2012-11-14 12:15:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

引用回帖:
7楼: Originally posted by lishijunzong at 2012-11-13 18:40:15
之前的作业是系统错误,你从新提交作业就行...

是算一次错误了之后再提交计算就能算出来?那之前的要删除吗?之前的要保留chk文件吗?能不能再详细一点啊谢谢!
10楼2012-11-14 12:18:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiemin8761 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 06:40 by gy1nBQXYQJqL
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 05:13 by gy1nBQXYQJqL
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +4 gouxfjh 2026-08-28 4/200 2026-08-29 04:46 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见