| 查看: 1084 | 回复: 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 价键不符合？？望指点修改 |
回复此楼» 猜你喜欢
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有52人回复
-
分子动力学模拟合作需求
已经有4人回复
-
求ENDNOTE破解软件安装包
已经有7人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有0人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有0人回复
njut007
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 144.5
- 帖子: 143
- 在线: 16.3小时
- 虫号: 436568
- 注册: 2007-08-25
- 性别: MM
- 专业: 有机分子功能材料化学
|
结构之前写错了，根据高斯结果调整程序如下： AM1 POLAR DIPOLE ENPART PRECISE SYMMETRY BONDS CHARGE=1 1 C 0.000000 0 0.000000 0 0.000000 0 0 0 0 2 N 1.319960 0 0.000000 0 0.000000 0 1 0 0 3 C 1.319053 1 110.870396 1 0.000000 0 2 1 0 4 C 1.405313 1 105.679277 1 -0.017126 0 3 2 1 5 N 1.319383 1 105.774905 1 -0.006257 0 4 3 2 6 C 1.441586 1 124.512243 1 -179.994714 0 2 1 5 7 C 1.442559 1 124.638443 1 179.988212 1 5 4 3 8 H 1.070000 1 126.558156 1 179.961790 1 1 5 4 9 H 1.070000 1 127.160361 1 179.982874 1 3 2 1 10 H 1.070000 1 127.112547 1 179.993743 1 4 3 2 11 H 1.070000 1 109.471221 1 0.000000 1 6 2 1 12 H 1.070000 1 109.471221 1 -120.000000 1 6 2 1 13 H 1.070000 1 109.471221 1 120.000000 1 6 2 1 14 H 1.070000 1 109.471221 1 -179.955346 1 7 5 4 15 H 1.070000 1 109.471221 1 -59.955346 1 7 5 4 16 H 1.070000 1 109.471221 1 60.044654 1 7 5 4 运行后出错的结果和主题贴的情况是一样的 之前的程序后面的连接性出了错，少了个C原子，不好意思 要换用这组数据。 问题还是解决不了  我用的是MOPAC20002，不去掉DIPOLE也能运行出来，但是结果就是 GEOMETRY IS FAULTY. 后面的错误 如主题贴 [ Last edited by njut007 on 2007-11-3 at 20:19 ] |
14楼2007-11-03 20:08:57
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