| 查看: 1041 | 回复: 25 | |||
| 当前主题已经存档。 | |||
njut007银虫 (小有名气)
|
[交流]
求助：MOPAC优化阳离子程序问题
|
||
|
用MOPAC优化一个阳离子，自己编了输入的程序，但是运行老是出错 错误情况如下： 原部分编程如下（其中第二个N为+1价）： C 0.000000 0 0.000000 0 0.000000 0 0 0 0 N 1.319960 0 0.000000 0 0.000000 0 1 0 0 C 1.319053 1 110.870396 1 0.000000 0 2 1 0 C 1.405313 0 105.679277 0 -0.017126 0 3 2 1 N 1.319383 0 105.774905 0 -0.006257 0 1 2 3 C 1.441586 0 124.512243 0 -179.994714 0 2 1 5 C 1.442559 1 124.638443 1 179.988212 1 5 1 2 H 1.070000 1 126.558156 1 179.961790 1 1 2 3 H 1.070000 0 127.160361 0 179.982874 0 3 2 1 H 1.070000 0 127.112547 0 179.993743 0 4 3 2 H 1.070000 0 109.471221 0 0.000000 0 6 2 1 H 1.070000 1 109.471221 1 -120.000000 1 6 2 1 H 1.070000 1 109.471221 1 120.000000 1 6 2 1 但是运行结果提示读入出错： The fault detected is in the last line printed GEOMETRY IS FAULTY. GEOMETRY READ IN IS： C 1.405313 0 105.679277 0 -0.017126 0 3 2 1 N 1.319383 0 105.774905 0 -0.006257 0 1 C 1.441586 0 124.512243 0 -179.9947 1 4 0 2 N 1.442559 1 1.000000 1 1.000000 1 1 5 0 O 1.070000 1 1.000000 1 1.000000 1 1 1 0 The fault detected is in the last line printed GEOMETRY IS FAULTY. GEOMETRY READ IN IS： H 1.070000 0 127.160361 0 179.982874 0 3 2 1 H 1.070000 0 127.112547 0 179.993743 0 4 H 1.070000 0 109.471221 0 0.000000 0 6 Mg 1.070000 1 1.000000 1 1.000000 -1 1 6 0 Al 1.070000 1 1.000000 1 1.000000 1 1 6 0 不知道是什么原因，会不会是跟我优化的是阳离子有关，系统默认 N 价键不符合？？望指点修改 |
回复此楼» 猜你喜欢
分子筛催化剂
已经有1人回复
-
青岛大学国家优青团队急招一名26年秋季入学博士（报名截止时间为6月17号下午1点）
已经有8人回复
-
物理化学论文润色/翻译怎么收费?
已经有178人回复
-
福建理工大学资源化学专业正式获批并启动招生
已经有32人回复
-
Crystals（IF=2.9）- Special Issue invatation- Functional Materials
已经有0人回复
-
分子筛催化方向 一篇cej共一第一 可以申请上海什么学校呢
已经有3人回复
-
求助共晶分子筛的电镜图分析，有偿
已经有1人回复
-
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有5人回复
-
关于合成MIL-125的求助
已经有0人回复
snoopyzhao
至尊木虫 (职业作家)
- QC强帖: 1
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2007-11-03 11:38:08
alwens
铁杆木虫 (正式写手)
老木虫
- 应助: 0 (幼儿园)
- 贵宾: 0.45
- 金币: 5208.8
- 散金: 50
- 帖子: 994
- 在线: 807.1小时
- 虫号: 8486
- 注册: 2003-04-27
- 性别: GG
- 专业: 药物设计与药物信息
3楼2007-11-03 14:26:28
njut007
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 144.5
- 帖子: 143
- 在线: 16.3小时
- 虫号: 436568
- 注册: 2007-08-25
- 性别: MM
- 专业: 有机分子功能材料化学
4楼2007-11-03 16:08:50
njut007
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 144.5
- 帖子: 143
- 在线: 16.3小时
- 虫号: 436568
- 注册: 2007-08-25
- 性别: MM
- 专业: 有机分子功能材料化学
5楼2007-11-03 16:20:48
njut007
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 144.5
- 帖子: 143
- 在线: 16.3小时
- 虫号: 436568
- 注册: 2007-08-25
- 性别: MM
- 专业: 有机分子功能材料化学
6楼2007-11-03 16:24:22
alwens
铁杆木虫 (正式写手)
老木虫
- 应助: 0 (幼儿园)
- 贵宾: 0.45
- 金币: 5208.8
- 散金: 50
- 帖子: 994
- 在线: 807.1小时
- 虫号: 8486
- 注册: 2003-04-27
- 性别: GG
- 专业: 药物设计与药物信息
7楼2007-11-03 17:40:19
snoopyzhao
至尊木虫 (职业作家)
- QC强帖: 1
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
|
AM1 POLAR DIPOLE ENPART PRECISE SYMMETRY BONDS CHARGE=1 [No Title] C 0.000000 0 0.000000 0 0.000000 0 0 0 0 N 1.319960 0 0.000000 0 0.000000 0 1 0 0 C 1.319053 1 110.870396 1 0.000000 0 2 1 0 C 1.405313 0 105.679277 0 -0.017126 0 3 2 1 N 1.319383 0 105.774905 0 -0.006257 0 1 2 3 C 1.441586 0 124.512243 0 -179.994714 0 2 1 5 C 1.442559 1 124.638443 1 179.988212 1 5 1 2 H 1.070000 1 126.558156 1 179.961790 1 1 2 3 H 1.070000 0 127.160361 0 179.982874 0 3 2 1 H 1.070000 0 127.112547 0 179.993743 0 4 3 2 H 1.070000 0 109.471221 0 0.000000 0 6 2 1 H 1.070000 1 109.471221 1 -120.000000 1 6 2 1 H 1.070000 1 109.471221 1 120.000000 1 6 2 1 ------ 这是你的输入文件，注意最后面的那个空行，严格意义上来说是必须的。但似乎也不绝对。 只是关键词部分似乎有一些问题。比如 chemoffice 2004 中的 mopac 2000 提示，POLAR 不能用在开壳的RHF 计算中。但 MOPAC 2007 便无此提示。 不过在 2007 中, DIPOLE 成了不可识别的关键词，呵呵 |
8楼2007-11-03 17:41:04
alwens
铁杆木虫 (正式写手)
老木虫
- 应助: 0 (幼儿园)
- 贵宾: 0.45
- 金币: 5208.8
- 散金: 50
- 帖子: 994
- 在线: 807.1小时
- 虫号: 8486
- 注册: 2003-04-27
- 性别: GG
- 专业: 药物设计与药物信息
9楼2007-11-03 17:59:52
snoopyzhao
至尊木虫 (职业作家)
- QC强帖: 1
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
10楼2007-11-03 18:05:06