24小时热门版块排行榜    

查看: 1624  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 16/800 2026-08-07 13:25 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 基金中了 +11 laoda193707 2026-08-06 11/550 2026-08-07 10:54 by TT和NN
[基金申请] 听说今天filecode变了 +21 布布和一二 2026-08-06 44/2200 2026-08-07 10:11 by 布布和一二
[基金申请] 固定端突然变了,今天 +4 archvillain 2026-08-06 8/400 2026-08-07 08:06 by 带你看大海
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见