24小时热门版块排行榜    

查看: 1212  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

onlyyou2002

铁虫 (初入文坛)

[求助] Gaussian计算竟然连initialization这一步都过不去。有人遇到过吗?

做一个optimization的计算,结果刚读入输入文件计算就不动了, 连“initialization pass" 的信息都没显示。也没有显示任何出错信息。下面是输入和输出文件。有哪位遇到过吗?谢谢

输入文件:

%chk=x.chk
%mem=40GB
%proc=12
# opt=(maxcycle=200) MaxDisk=300GB freq=noraman b3lyp/6-31g(d,p)

test

0 1
C                    -1   -4.94288   3.45349  -5.99646
C                    -5.09602   2.65634  -4.69502
C                    -4.0512    1.54721  -4.58157
.........................
H                     4.97687   2.91147   1.87889
H                     3.24659   3.29036   1.99794
H                     2.75241   2.10036  -1.18513
H                     2.71944  -3.6743   -2.52586
H                     4.40356  -6.08022  -2.19961
H                     4.58666  -4.66745  -0.14075
H                     5.29371  -2.5623    0.8192
H                     6.75951  -3.57167   0.84371
H                     6.50185  -2.42007  -0.48227

输出文件:
-----------------------------------------------------------------
#p opt=(maxcycle=200) MaxDisk=300GB freq=noraman b3lyp/6-31g(d,p)
-----------------------------------------------------------------
1/6=200,14=-1,18=20,19=15,26=3,38=1/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7//1,2,3,16;
1/6=200,14=-1,18=20,19=15/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/6=200,14=-1,18=20,19=15/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
Leave Link    1 at Thu Jan  5 18:49:44 2012, MaxMem=10737418240 cpu:       0.4
(Enter /share/apps/g09/l101.exe)
----
test
----
Charge =  0 Multiplicity = 1
Symbolic Z-Matrix:
C                    -1   -4.94288   3.45349  -5.99646
C                    -5.09602   2.65634  -4.69502
C                    -4.0512    1.54721  -4.58157
.........................
H                     4.97687   2.91147   1.87889
H                     3.24659   3.29036   1.99794
H                     2.75241   2.10036  -1.18513
H                     2.71944  -3.6743   -2.52586
H                     4.40356  -6.08022  -2.19961
H                     4.58666  -4.66745  -0.14075
H                     5.29371  -2.5623    0.8192
H                     6.75951  -3.57167   0.84371
H                     6.50185  -2.42007  -0.48227
回复此楼
人生如梦
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyichem

新虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你采用opt=modredundant 然后把你要固定的原子序数写在坐标后面。具体例子是这样的。
opt=modredundant

坐标

X 1 F
X 2 F
18楼2012-01-08 09:44:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

pwzhou

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你仔细看看第一个C原子的坐标,有4个数。
2楼2012-01-06 10:42:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onlyyou2002

铁虫 (初入文坛)

引用回帖:
: Originally posted by pwzhou at 2012-01-06 10:42:18:
你仔细看看第一个C原子的坐标,有4个数。

第一个C原子坐标有四个数,是因为我优化的时候把它的坐标固定了。
我也尝试把它的坐标不固定,情况还是一样。
人生如梦
3楼2012-01-06 10:50:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zln1987

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你把# opt=(maxcycle=200) MaxDisk=300GB freq=noraman b3lyp/6-31g(d,p) 的# 后加上p,看看能不能出来较为详细的提示。
4楼2012-01-06 11:04:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见