| 查看: 1394 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 过渡态计算 |
» 猜你喜欢
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士等
已经有20人回复
-
物理化学论文润色/翻译怎么收费?
已经有55人回复
-
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士
已经有23人回复
-
德国Karlsruhe Institute of Technology招收电化学储能及联合培养CSC博士
已经有1人回复
-
固态锂金属电池迁移数计算求助
已经有1人回复
-
荷兰奈梅亨大学赵文博与Galimberti课题组招收理论计算CSC博士 2000 欧元/月+房补
已经有0人回复
-
26储能博士申请自荐
已经有22人回复
-
cif求助下载
已经有1人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有2人回复
-
上海交通大学--宁波东方理工大学电池方向博士招生
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
xy00083789
金虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 932.7
- 散金: 4
- 帖子: 82
- 在线: 124.6小时
- 虫号: 1157021
- 注册: 2010-11-27
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★
感谢参与，应助指数 +1
xiemin8761: 金币+1, ★有帮助 2012-11-14 12:19:07
xiemin8761: 金币+3, ★★★很有帮助, 我再看看能不能算。谢谢 2012-11-14 16:33:22
感谢参与，应助指数 +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里面，定义对该键长进行扫描。 |
8楼2012-11-13 18:58:56
xiemin8761
铁杆木虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 8666
- 散金: 188
- 红花: 2
- 帖子: 2213
- 在线: 509.4小时
- 虫号: 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
lishijunzong
铁杆木虫 (正式写手)
- 应助: 214 (大学生)
- 金币: 5034
- 红花: 10
- 帖子: 574
- 在线: 195.1小时
- 虫号: 1121020
- 注册: 2010-10-13
- 性别: GG
- 专业: 理论和计算化学
3楼2012-11-12 19:43:50
xiemin8761
铁杆木虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 8666
- 散金: 188
- 红花: 2
- 帖子: 2213
- 在线: 509.4小时
- 虫号: 769782
- 注册: 2009-05-14
- 性别: GG
- 专业: 化学动力学
4楼2012-11-13 08:13:48
10