24小时热门版块排行榜    

查看: 1336  |  回复: 16
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

scjiang

银虫 (正式写手)

[交流] g03w B.01版本中,叶立德型有机化合物计算优势构象基组选择

各位前辈好,我是刚刚接触g03w没多久的菜鸟,看了很多手册对很多问题还是不太明白,希望得到您们的帮助。

我需要计算一个N-N型叶立德化合物(含有最大原子为氯)的N-N键长,现在遇到了问题。我先描述一下我操作吧。

1, 先在Chem3D Ultra 10.0里面使用Calculations/MM2/Minimiz Energy初步优化了一下结构,点击叶立德的正电N原子和负电N原子发现显示的一个是charge:+1.00;一个是charge:-1.00,于是便用另存为选项保存为Gaussian Input file(gjf格式)。

2, 用UltraEdit10打开gjf格式的文件,发现Chem3D使用的是笛卡尔直角坐标系表示的分子几何构形,里面没有电荷信息。但是我还是把它拿来运行g3w了。我在运行说明行和其后添加的内容是:
#HF/6-31G Opt

[No title]

0  1
C 0   -3.590182    1.533814   -0.029880
C 0   -3.484144    0.123764   -0.034077
.
.
.
N 0    0.224491   -0.263694   -0.001520
.
.
N 0    0.642892   -1.677703   -0.011997
.
.
.
H 0    5.824028   -0.951299    0.044691
H 0    4.228762    0.961336    0.055007

3,运行winows单机版本的g03w后每次都跳出错误信息:
Severe Error Message # 2070

The processing of the last link ended abnormally.
All processing has been aborted.
打开输出文件.out发现最后内容是:
Rotational constants (GHZ):      1.4189680      0.2965498      0.2452883
Standard basis: 6-31G (6D, 7F)
The combination of multiplicity 1 and   115 electrons is impossible.
Error termination via Lnk1e in F:\G03W\l301.exe at Mon Apr 30 13:19:33 2007.
Job cpu time:  0 days  0 hours  0 minutes  2.0 seconds.
File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=      1 Scr=      1

4, 查找帮助手册知道了l301是用来生成基组信息的,我想应该是基组选择出错了,于是把它改成最小的基组STO-3G结果也是一样,改成6-311G和6-31G**也是如此!

以上就是我的操作了。我现在想问的有以下几条:

1,叶立德型化合物(含有最大原子为氯)在HF方法里应该用什么基组计算最优构象?

2,描述叶立德型化合物(含有最大原子为氯)的分子构象是用笛卡尔直角坐标还是用z矩阵?我直接用Chem3D 10.0转换成的输入文件是否正确?如错我应该怎么改正?

3,Charge关键字指定在计算中包含背景电荷的分布。我是不是可以在笛卡尔直角坐标里面含电荷的原子用charge关键字标示出来?标示出来又应该怎么标示?我实在没法理解03中文手册上的那个例子!!

4, 叶立德型化合物的净电荷为0,我输入的电荷与自旋多重度有没有错?

5,我是否能将得到的已经优化的输出文件转化成Chem3D那种点击原子就能看见键长键角信息的格式?最后得到的Optimized Parameters应该怎样去转化成可视化的信息?

实在是太谢谢各位走过路过的前辈帮忙了!这个问题困扰了我好长时间了都还是没有搞定!郁闷啊

我的email :wcumsjy@hotmail.com

附件就是我要计算的结构!

[ Last edited by zzgyb on 2008-3-28 at 09:19 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scjiang

银虫 (正式写手)

能够正常运行了!
再次谢谢niliu前辈!
11楼2007-04-30 16:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

niliu

铁杆木虫 (著名写手)

★ ★ ★
bloveocean(金币+3):解答非常完美,多加点。
从出错提示来看,你输入的体系有115个电子,应该是输入的分子结构有问题。不妨把整个输入文件贴出来看看。
2楼2007-04-30 15:16:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scjiang

银虫 (正式写手)

#HF/6-31G opt

[No Title]

  0  1
C 0   -3.590182    1.533814   -0.029880
C 0   -3.484144    0.123764   -0.034077
C 0   -2.242338   -0.479991   -0.025158
C 0   -1.068742    0.315800   -0.010300
C 0   -1.167222    1.720186   -0.002884
C 0   -2.454718    2.319490   -0.014380
N 0    0.224491   -0.263694   -0.001520
C 0    1.342083    0.524811    0.016961
C 0    1.267951    1.951947    0.026512
C 0    0.027363    2.529237    0.015573
N 0    0.642892   -1.677703   -0.011997
C 0    2.012969   -1.685058   -0.000172
C 0    2.474260   -0.354707    0.019852
C 0    2.920361   -2.774118   -0.008044
C 0    4.271132   -2.485702    0.007472
C 0    4.738433   -1.143828    0.031464
C 0    3.860933   -0.075820    0.037349
H 0   -4.585705    2.006994   -0.038547
H 0   -4.397562   -0.493036   -0.044867
H 0   -2.168710   -1.579086   -0.028775
H 0   -2.556594    3.416962   -0.010999
H 0    2.181956    2.564764    0.040250
H 0   -0.052314    3.628602    0.021069
H 0    0.011174   -2.510237   -0.028795
H 0    2.554575   -3.811990   -0.026889
H 0    5.001695   -3.311438    0.001083
H 0    5.824028   -0.951299    0.044691
H 0    4.228762    0.961336    0.055007
3楼2007-04-30 15:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scjiang

银虫 (正式写手)

上面就是所有的输入文件!太感谢niliu帮忙了!
4楼2007-04-30 15:24:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见