24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiemin8761

铁杆木虫 (著名写手)

应助: 13 (小学生)
金币: 8776
散金: 164
红花: 2
帖子: 2212
在线: 508.9小时
虫号: 769782
注册: 2009-05-14
性别: GG
专业: 化学动力学

[求助] 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

回复此楼

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

过渡态计算

» 猜你喜欢

华南师范大学先进光电子研究院电子纸团队诚聘2026级博士研究生已经有0人回复
QE利用声子计算的一系列lambda和Tc如何选已经有8人回复
物理化学论文润色/翻译怎么收费? 已经有196人回复
理论计算合作已经有4人回复
小木虫的论文辅导靠谱吗？有没有用过的同学？已经有1人回复
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有已经有0人回复
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后已经有81人回复
书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复
七嗪类物质合成求助已经有0人回复
东方理-中科大联合博士生招聘已经有0人回复

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

1楼 2012-11-12 11:22:38

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

xiemin8761

铁杆木虫 (著名写手)

应助: 13 (小学生)
金币: 8776
散金: 164
红花: 2
帖子: 2212
在线: 508.9小时
虫号: 769782
注册: 2009-05-14
性别: GG
专业: 化学动力学

引用回帖:

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的回帖

查看全部 19 个回答

xiemin8761

铁杆木虫 (著名写手)

应助: 13 (小学生)
金币: 8776
散金: 164
红花: 2
帖子: 2212
在线: 508.9小时
虫号: 769782
注册: 2009-05-14
性别: GG
专业: 化学动力学

还有时候出错显示
。。。
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的回帖