24小时热门版块排行榜

返回列表

当前主题已经存档。

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

njut007

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 144.5
帖子: 143
在线: 16.3小时
虫号: 436568
注册: 2007-08-25
性别: MM
专业: 有机分子功能材料化学

[交流] 求助：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 价键不符合？？望指点修改

回复此楼