| 查看: 906 | 回复: 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 价键不符合？？望指点修改 |
回复此楼» 猜你喜欢
国家级人才课题组招收2026年入学博士
已经有2人回复
-
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有173人回复
-
散金
已经有87人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
已经有34人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
南通大学电气与自动化学院2026年优秀博士招聘公告
已经有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
- 金币: 5248.8
- 散金: 50
- 帖子: 994
- 在线: 805.8小时
- 虫号: 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
- 金币: 5248.8
- 散金: 50
- 帖子: 994
- 在线: 805.8小时
- 虫号: 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
- 金币: 5248.8
- 散金: 50
- 帖子: 994
- 在线: 805.8小时
- 虫号: 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