计算化学公社

标题: gamess输入坐标问题请教 [打印本页]

作者
Author: 413 时间: 2019-11-26 22:26
标题: gamess输入坐标问题请教
请教一下大家，我在用gamess做测试时，算了一个水分子，也输入了一个水分子的坐标，可是为什么log中出现了4个水？
以下为gamess的部分输出
INPUT CARD> $contrl scftyp=rhf runtyp=gradient dfttyp=b3lyp $end
INPUT CARD> $SYSTEM TIMLIM=1 $END
INPUT CARD> $BASIS  GBASIS=N31 NGAUSS=6 NDFUNC=1 NPFUNC=1 DIFFSP=.T. $END
INPUT CARD> $GUESS  GUESS=HUCKEL $END
INPUT CARD> $DATA
INPUT CARD>h2o test
INPUT CARD>Cnv  2
INPUT CARD>
INPUT CARD>O 8.0          6.00372672 -1.18719932 0.00000000
INPUT CARD>H 1.0          6.96372672 -1.18719932 0.00000000
INPUT CARD>H 1.0          5.68327213 -0.28226348 0.00000000
INPUT CARD> $END
INPUT CARD>
1000000 WORDS OF MEMORY AVAILABLE

   BASIS OPTIONS
   -------------
   GBASIS=N31       IGAUSS=    6    POLAR=POPN31
   NDFUNC=    1    NFFUNC=    0    DIFFSP=    T
   NPFUNC=    1    DIFFS=    F    BASNAM=

   RUN TITLE
   ---------
Methylene...3-B-1 state...UHF/STO-2G

THE POINT GROUP OF THE MOLECULE IS CNV
THE ORDER OF THE PRINCIPAL AXIS IS    2

ATOM    ATOMIC                   COORDINATES (BOHR)
         CHARGE       X                Y                Z
O          8.0 -11.3453984060       2.2434814076       0.0000000000
O          8.0 11.3453984060       2.2434814076       0.0000000000
O          8.0 -11.3453984060    -2.2434814076       0.0000000000
O          8.0 11.3453984060    -2.2434814076       0.0000000000
H          1.0 -13.1595353542       2.2434814076       0.0000000000
H          1.0 13.1595353542       2.2434814076       0.0000000000
H          1.0 -13.1595353542    -2.2434814076       0.0000000000
H          1.0 13.1595353542    -2.2434814076       0.0000000000
H          1.0 -10.7398270394       0.5334006335       0.0000000000
H          1.0 10.7398270394       0.5334006335       0.0000000000
H          1.0 -10.7398270394    -0.5334006335       0.0000000000
H          1.0 10.7398270394    -0.5334006335       0.0000000000

.............

                     ----------------------
                     GRADIENT OF THE ENERGY
                     ----------------------

UNITS ARE HARTREE/BOHR E'X             E'Y             E'Z
1 O             -0.007125295    -0.087805310    0.000000000
2 O             0.007125295    -0.087805310    0.000000000
3 O             -0.007125295    0.087805310    0.000000000
4 O             0.007125295    0.087805310    0.000000000
5 H             0.004337114    0.005865359    0.000000000
6 H             -0.004337114    0.005865359    0.000000000
7 H             0.004337114    -0.005865359    0.000000000
8 H             -0.004337114    -0.005865359    0.000000000
9 H             0.002792047    -0.402880189    0.000000000
10 H             -0.002792047    -0.402880189    0.000000000
11 H             0.002792047    0.402880189    0.000000000
12 H             -0.002792047    0.402880189    0.000000000

作者
Author: zjxitcc 时间: 2019-11-26 22:31

因为你开启了对称性，此时只能输入对称唯一的原子坐标，也就是一个O和一个H。要么就不用对称性。

作者
Author: 413 时间: 2019-11-26 22:33

zjxitcc 发表于 2019-11-26 22:31
因为你开启了对称性，此时只能输入对称唯一的原子坐标，也就是一个O和一个H。要么就不用对称性。

请教一下，是这个吗？Cnv 2，指定了C2v就相当于使用了对称性？

是不是我后面计算所有的分子都用C1就可以了？

作者
Author: zjxitcc 时间: 2019-11-26 22:39

413 发表于 2019-11-26 22:33
请教一下，是这个吗？Cnv 2，指定了C2v就相当于使用了对称性？

是不是我后面计算所有的分子都用C1就 ...

是的。

不用对称性写C1[空格]1.

看你想用还是不想用了。

作者
Author: 413 时间: 2019-11-26 22:40

zjxitcc 发表于 2019-11-26 22:39
是的。

不用对称性写C1[空格]1.

不想用，算的分子基本都没有对称性。
谢谢谢谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)