Посчитать орбиты в Gaussian

igprot

Надо посчитать в Gaussian -е виды потенциальных кривых электронных орбит (основных и возбужденных) в зависимости от межядерного рассояния. Для начала хотел бы разобраться с наиболее простыми молекулами - H2 и H2+
Хотел бы знать из чего состоит входной файл и что именно нужно бать из выходного.
В Gaussian -е я полное ламье =) только начал изучать, но времени уже не хватает.

Serge_FL2

виды потенциальных кривых электронных орбит
1.Ты ставишь меня в ступор реальне этой фразой.
2. H2+ считается по сути аналитически, программа тут не нужна-в книжке формулы посмотри(сделай мо лкао из орбиталей атома водорода, которые в аналитическом виде есть). А в книжке небось еще и все формулы готовые есть.
3.Гауссиан это зло, пользуй гамесс.

terl

Надо посчитать в Gaussian -е виды потенциальных кривых электронных орбит (основных и возбужденных) в зависимости от межядерного рассояния.
Это называется электронные термы основного и возбуженных состояний или просто поверхности потенциальной энергии.

sidy_mos

Подскажите пожалуйста как считать термы основного и возбужденых состояний у молекул в gaussian-е или еще в чем нибудь? Если кто-то может подсказать как считается в гамесе, то скажите в нем.
В gaussian-е я понял как считать полную энергию и разницу энергий для заданого междуядерного расстояния. Осталось только понять как считать терм основного состояния и тогда можно будет сложить его энергию и их разницу.

terl

Могу про GAMESS рассказать.
Нужно написать текстовый входной файл. Для расчета электронного терма основного состояния H2 (в рамках ограниченного метода Хартри-Фока, в базисе 6-31G, впрочем если это водород, то на 6- пофигу) он будет такой:

$contrl
scftyp=rhf
runtyp=surface
icharg=0
mult=1
coord=unique
units=angs
$end

$scf
dirscf=.t.
$end

$guess
guess=huckel
$end

$basis
gbasis=n31
ngauss=6
$end

$surf
ivec1(1)=1,2
igrp1(1)=2
orig1=-0.4
disp1=0.05
ndisp1=86
$end

$data
h2-gs
C1
H 1.0 0.0 0.0 0.0
H 1.0 0.0 0.0 0.7
$end

Что все это значит можно посмотреть в документации к GAMESS
На выходе получится тоже текстовый файл. В случае PCGAMESS в конце будет таблица зависимости полной энергии от шага. В GAMESS(US) не помню пишется ли в конце сводная таблица, но в любом случае энергии можно вытащить скажем grep'ом по строке "FINAL ENERGY IS".
Если нужны зависимости орбитальных энергий от расстояний, то придется генерировать координаты атомов внешней утилитой и для каждой геометрии считать энергию (runtyp=energy а потом вытаскивать их из выходного файла. Ну или надо заставить гамесс выдавать вектора и орбитальные энергии на каждом шаге сканирования, но че-то не найду как.
Если нужны электронные термы возбужденных состояний - считай конфигурационное взаимодействие (cityp) или многоконфигурационное взаимодействие (mcscf). Но об этом в следующей серии. А еще лучше у спроси - он должен лучше знать.
зы
и да, пробелы в начале строк во входном файле _нужны_

strelok69

не пашет :smirk:

1 ******************************************************
* GAMESS VERSION = 6 JUN 1999 *
* FROM IOWA STATE UNIVERSITY *
* M.W.SCHMIDT, K.K.BALDRIDGE, J.A.BOATZ, S.T.ELBERT, *
* M.S.GORDON, J.H.JENSEN, S.KOSEKI, N.MATSUNAGA, *
* K.A.NGUYEN, S.J.SU, T.L.WINDUS, *
* TOGETHER WITH M.DUPUIS, J.A.MONTGOMERY *
* J.COMPUT.CHEM. 14, 1347-1363(1993) *
*******Intel x86 (Win32,Linux,OS/2,DOS) VERSION*******
* PC GAMESS version 6.4, build number 2606 *
* Compiled on Thursday, 04-03-2004, 15:45:41 *
* Intel specific optimization, bug fixes, *
* code changes, and additional functionality - *
* copyright (c) 1994, 2004 by Alex. A. Granovsky, *
* Laboratory of Chemical Cybernetics, *
* Moscow State University, Moscow, Russia. *
* Some parts of this program include code due to *
* work of Jim Kress, Peter Burger, and Robert Ponec. *
******************************************************
* PC GAMESS URLs: *
* http://classic.chem.msu.su/gran/gamess/index.html *
* http://quantum-2.chem.msu.ru/gran/gamess/index.html*
* E-mail: classic.chem.msu.su *
* This program may not be redistributed without *
* the specific, written permission of its developers.*
******************************************************


Pentium II / Win32 PC GAMESS version running under Windows NT
Running on Intel CPU: Brand ID 22, Family 6, Model 13, Stepping 8
CPU Brand String : Intel(R) Pentium(R) M processor 1.86GHz
CPU Features : CMOV, MMX, SSE, SSE2
Data cache size : L1 32 KB, L2 0 KB, L3 0 KB
Operating System successfully passed SSE support test.

Warning: this PC GAMESS version was optimized for a different CPU type!
The overall performance degradation is expected.
For better performance, use Pentium III optimized PC GAMESS version.


EXECUTION OF GAMESS BEGUN 0:58:08 LT 25-OCT-2008

ECHO OF THE FIRST FEW INPUT CARDS -
INPUT CARD> $contrl
INPUT CARD> scftyp=rhf
INPUT CARD> runtyp=surface
INPUT CARD> icharg=0
INPUT CARD> mult=1
INPUT CARD> coord=unique
INPUT CARD> units=angs
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $scf
INPUT CARD> dirscf=.t.
INPUT CARD> diis=.t.
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $guess
INPUT CARD> guess=huckel
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $basis
INPUT CARD> gbasis=n31
INPUT CARD> ngauss=6
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $surf
INPUT CARD> ivec1(1)=1,2
INPUT CARD> igrp1(1)=2
INPUT CARD> orig1=-0.4
INPUT CARD> disp1=0.05
INPUT CARD> ndisp1=86
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $data
INPUT CARD>h2-gs
INPUT CARD>C1
INPUT CARD>H 1.0 0.0 0.0 0.0
INPUT CARD>H 1.0 0.0 0.0 0.7
INPUT CARD> $end
2000000 WORDS OF MEMORY AVAILABLE

BASIS OPTIONS
-------------
GBASIS=N31 IGAUSS= 6 POLAR=NONE
NDFUNC= 0 DIFFSP= F
NPFUNC= 0 DIFFS= F


RUN TITLE
---------
h2-gs

THE POINT GROUP OF THE MOLECULE IS C1
THE ORDER OF THE PRINCIPAL AXIS IS 0

ATOM ATOMIC COORDINATES (BOHR)
CHARGE X Y Z
H 1.0 0.0000000000 0.0000000000 0.0000000000
H 1.0 0.0000000000 0.0000000000 1.3228081914

INTERNUCLEAR DISTANCES (ANGS.)
------------------------------

H H

1 H 0.0000000 0.7000000 *
2 H 0.7000000 * 0.0000000

* ... LESS THAN 3.000


ATOMIC BASIS SET
----------------
THE CONTRACTED PRIMITIVE FUNCTIONS HAVE BEEN UNNORMALIZED
THE CONTRACTED BASIS FUNCTIONS ARE NOW NORMALIZED TO UNITY

SHELL TYPE PRIM EXPONENT CONTRACTION COEFFICIENTS

H

1 S 1 18.731137 0.214935 ( 0.033495)
1 S 2 2.825394 0.364571 ( 0.234727)
1 S 3 0.640122 0.415051 ( 0.813757)

2 S 4 0.161278 0.181381 ( 1.000000)

H

3 S 5 18.731137 0.214935 ( 0.033495)
3 S 6 2.825394 0.364571 ( 0.234727)
3 S 7 0.640122 0.415051 ( 0.813757)

4 S 8 0.161278 0.181381 ( 1.000000)

TOTAL NUMBER OF SHELLS = 4
TOTAL NUMBER OF BASIS FUNCTIONS = 4
NUMBER OF ELECTRONS = 2
CHARGE OF MOLECULE = 0
STATE MULTIPLICITY = 1
NUMBER OF OCCUPIED ORBITALS (ALPHA) = 1
NUMBER OF OCCUPIED ORBITALS (BETA ) = 1
TOTAL NUMBER OF ATOMS = 2
THE NUCLEAR REPULSION ENERGY IS 0.7559674989

THIS MOLECULE IS RECOGNIZED AS BEING LINEAR.

$CONTRL OPTIONS
---------------
SCFTYP=RHF RUNTYP=SURFACE EXETYP=RUN
MPLEVL= 0 LOCAL =NONE UNITS =ANGS
MULT = 1 ICHARG= 0 MAXIT = 30
NPRINT= 7 IREST = 0 COORD =UNIQUE
ECP =NONE NORMF = 0 NORMP = 0
ITOL = 20 ICUT = 9 NZVAR = 0
NOSYM = 0 INTTYP=POPLE GEOM =INPUT
PLTORB= F MOLPLT= F RPAC = F
AIMPAC= F FRIEND= CITYP =NONE
DFTTYP=NONE

$SYSTEM OPTIONS
---------------
KDIAG = 0 MEMORY= 2000000 TIMLIM= 36000.0 SEC.
COREFL= F PTIME = F R = F
BALTYP=LOOP

----------------
PROPERTIES INPUT
----------------

MOMENTS FIELD POTENTIAL DENSITY
IEMOM = 1 IEFLD = 0 IEPOT = 0 IEDEN = 0
WHERE =COMASS WHERE =NUCLEI WHERE =NUCLEI WHERE =NUCLEI
OUTPUT=BOTH OUTPUT=BOTH OUTPUT=BOTH OUTPUT=BOTH
IEMINT= 0 IEFINT= 0 IEDINT= 0
MORB = 0

EXTRAPOLATION IN EFFECT
DIIS IN EFFECT

----------------------
INTEGRAL INPUT OPTIONS
----------------------
NOPK = 1 NORDER= 0 SCHWRZ= T

ATTENTION! AO INTEGRALS WILL BE PACKED.
THRESHOLD FOR PACKING PKTHR = 0.10000000D-01

-------------------------------
INTEGRAL TRANSFORMATION OPTIONS
-------------------------------
NWORD = 0 CUTOFF = 1.0E-09
MPTRAN = 0 DIRTRF = T
AOINTS =DUP IREST = 0

------------------------------------------
THE POINT GROUP IS C1 , NAXIS= 0, ORDER= 1
------------------------------------------

DIMENSIONS OF THE SYMMETRY SUBSPACES ARE
A = 4

..... DONE SETTING UP THE RUN .....

CPU TIME: STEP = 0.02 , TOTAL = 0.1 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.17 , TOTAL = 0.2 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 9.04%, TOTAL = 45.20%

---------------------------
POTENTIAL SURFACE MAP INPUT
---------------------------
COORD 1 LYING ALONG ATOM PAIR 1 2
HAS ORIGIN=-0.400, DISPLACEMENT= 0.050 AND 86 STEPS.
GROUP 1 CONTAINS 1 ATOMS:
2
...... END OF ONE-ELECTRON INTEGRALS ......

CPU TIME: STEP = 0.00 , TOTAL = 0.1 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.00 , TOTAL = 0.2 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 44.66%

-------------
GUESS OPTIONS
-------------
GUESS =HUCKEL NORB = 0 NORDER= 0
MIX = F PRTMO = F SYMDEN= F
TOLZ = 1.0E-08 TOLE = 1.0E-05

INITIAL GUESS ORBITALS GENERATED BY HUCKEL ROUTINE.
HUCKEL GUESS REQUIRES 1424 WORDS.

SYMMETRIES FOR INITIAL GUESS ORBITALS FOLLOW. BOTH SET(S).
1 ORBITALS ARE OCCUPIED ( 0 CORE ORBITALS).
1=A 2=A 3=A 4=A
...... END OF INITIAL ORBITAL SELECTION ......

CPU TIME: STEP = 0.00 , TOTAL = 0.1 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.01 , TOTAL = 0.2 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 43.07%
...... END OF TWO-ELECTRON INTEGRALS .....

CPU TIME: STEP = 0.00 , TOTAL = 0.1 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.01 , TOTAL = 0.2 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 40.17%

DIRECT SCF CALCULATION, SCHWRZ=T FDIFF=T

NONZERO BLOCKS
ITER EX DEM TOTAL ENERGY E CHANGE DENSITY CHANGE DIIS ERROR INTEGRALS SKIPPED
* * * INITIATING DIIS PROCEDURE * * *
1 0 0 -0.615477086 -0.615477086 0.083273977 0.265268887 55 0
2 1 0 -0.644522586 -0.029045499 0.014881237 0.038303534 55 0
3 2 0 -0.645154533 -0.000631948 0.000000000 0.000541897 55 0
DIIS EQUATIONS ARE SINGULAR, BOMBING


EXECUTION OF GAMESS TERMINATED ABNORMALLY AT 0:58:08 LT 25-OCT-2008
197280 WORDS OF DYNAMIC MEMORY USED

CPU TIME: STEP = 0.00 , TOTAL = 0.1 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.01 , TOTAL = 0.2 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 38.19%
STOP IN ABRT

terl

убери diis в $scf

sidy_mos

Спасибо за советы.
Встает ребром следующий вопрос, где можно найти программу gamess ? =)

sidy_mos

Наконец-то обзавелся Гамессом и начал экспериментировать =)
Я правильно понял, что код, который тут приведен он для основного состояния?
А можно считать еще и возбужденные?

terl

да, то был входной файл для основного состояния.
вот _пример_ расчета терма первого возбужденного состояния в рамках КВ:

$contrl
scftyp=rohf
runtyp=surface
icharg=0
mult=1
coord=unique
units=angs
cityp=guga
$end

$scf
dirscf=.t.
! diis=.t.
npunch=2
$end

$basis
gbasis=n31
ngauss=6
$end

$cidrt
fors=.t.
nfzc=0
ndoc=1
nval=3
$end

$gugdia
nstate=3
$end

$gugdm
iroot=2
NFLGDM(1)=1,1,1
$end


$surf
ivec1(1)=1,2
igrp1(1)=2
orig1=-0.4
disp1=0.1
ndisp1=40
$end

$data
h2-e1
C1
H 1.0 0.0 0.0 0.0
H 1.0 0.0 0.0 0.7
$end


зы: я правда хз че получится, возможно нужно добавить диффузных функций. вобщем это пример )

sidy_mos

Что-то он ругается... :
 1         ******************************************************
* GAMESS VERSION = 6 JUN 1999 *
* FROM IOWA STATE UNIVERSITY *
* M.W.SCHMIDT, K.K.BALDRIDGE, J.A.BOATZ, S.T.ELBERT, *
* M.S.GORDON, J.H.JENSEN, S.KOSEKI, N.MATSUNAGA, *
* K.A.NGUYEN, S.J.SU, T.L.WINDUS, *
* TOGETHER WITH M.DUPUIS, J.A.MONTGOMERY *
* J.COMPUT.CHEM. 14, 1347-1363(1993) *
*******Intel x86 (Win32,Linux,OS/2,DOS) VERSION*******
* PC GAMESS version 6.2, build number 2068 *
* Compiled on Friday, 21-09-2001, 00:37:32 *
* Intel specific optimization, bug fixes, *
* code changes, and additional functionality - *
* copyright (c) 1994, 2001 by Alex. A. Granovsky, *
* Laboratory of Chemical Cybernetics, *
* Moscow State University, Moscow, Russia. *
* Some parts of this program include code due to *
* work of Jim Kress, Peter Burger, and Robert Ponec. *
******************************************************
* PC GAMESS URLs: *
* http://classic.chem.msu.su/gran/gamess/index.html *
* http://quantum-2.chem.msu.ru/gran/gamess/index.html*
* E-mail: classic.chem.msu.su *
* This program may not be redistributed without *
* the specific, written permission of its developers.*
******************************************************


Pentium III/ Win32 PC GAMESS version running under Windows NT
Running on Intel CPU: Brand ID 22, Family 6, Model 13, Stepping 8
CPU Features : CMOV, MMX, SSE, SSE2
Data cache size : L1 16 KB, L2 0 KB
Operating system successfully passed SSE support test.


EXECUTION OF GAMESS BEGUN 18:18:22 LT 30-OCT-2008

ECHO OF THE FIRST FEW INPUT CARDS -
INPUT CARD>! H2 surface scan
INPUT CARD>!
INPUT CARD> $contrl
INPUT CARD> scftyp=rhf
INPUT CARD> runtyp=surface
INPUT CARD> icharg=0
INPUT CARD> mult=1
INPUT CARD> coord=unique
INPUT CARD> units=angs
INPUT CARD> cityp=guga
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $scf
INPUT CARD> dirscf=.t.
INPUT CARD>! diis=.t.
INPUT CARD> npunch=2
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $basis
INPUT CARD> gbasis=n31
INPUT CARD> ngauss=6
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $cidrt
INPUT CARD> fors=.t.
INPUT CARD> nfzc=0
INPUT CARD> ndoc=1
INPUT CARD> nval=3
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $gugdia
INPUT CARD> nstate=3
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $gugdm
INPUT CARD> iroot=2
INPUT CARD> NFLGDM(1)=1,1,1
INPUT CARD> $end
INPUT CARD>
INPUT CARD>
INPUT CARD> $surf
INPUT CARD> ivec1(1)=1,2
INPUT CARD> igrp1(1)=2
INPUT CARD> orig1=-0.4
INPUT CARD> disp1=0.1
INPUT CARD> ndisp1=40
INPUT CARD> $end
INPUT CARD>
INPUT CARD> $data
INPUT CARD>h2-e1
2000000 WORDS OF MEMORY AVAILABLE

BASIS OPTIONS
-------------
GBASIS=N31 IGAUSS= 6 POLAR=NONE
NDFUNC= 0 DIFFSP= F
NPFUNC= 0 DIFFS= F


RUN TITLE
---------
h2-e1

THE POINT GROUP OF THE MOLECULE IS C1
THE ORDER OF THE PRINCIPAL AXIS IS 0

ATOM ATOMIC COORDINATES (BOHR)
CHARGE X Y Z
H 1.0 0.0000000000 0.0000000000 0.0000000000
H 1.0 0.0000000000 0.0000000000 1.3228081914

INTERNUCLEAR DISTANCES (ANGS.)
------------------------------

H H

1 H 0.0000000 0.7000000 *
2 H 0.7000000 * 0.0000000

* ... LESS THAN 3.000


ATOMIC BASIS SET
----------------
THE CONTRACTED PRIMITIVE FUNCTIONS HAVE BEEN UNNORMALIZED
THE CONTRACTED BASIS FUNCTIONS ARE NOW NORMALIZED TO UNITY

SHELL TYPE PRIM EXPONENT CONTRACTION COEFFICIENTS

H

1 S 1 18.731137 0.214935 ( 0.033495)
1 S 2 2.825394 0.364571 ( 0.234727)
1 S 3 0.640122 0.415051 ( 0.813757)

2 S 4 0.161278 0.181381 ( 1.000000)

H

3 S 5 18.731137 0.214935 ( 0.033495)
3 S 6 2.825394 0.364571 ( 0.234727)
3 S 7 0.640122 0.415051 ( 0.813757)

4 S 8 0.161278 0.181381 ( 1.000000)

TOTAL NUMBER OF SHELLS = 4
TOTAL NUMBER OF BASIS FUNCTIONS = 4
NUMBER OF ELECTRONS = 2
CHARGE OF MOLECULE = 0
STATE MULTIPLICITY = 1
NUMBER OF OCCUPIED ORBITALS (ALPHA) = 1
NUMBER OF OCCUPIED ORBITALS (BETA ) = 1
TOTAL NUMBER OF ATOMS = 2
THE NUCLEAR REPULSION ENERGY IS 0.7559674989

THIS MOLECULE IS RECOGNIZED AS BEING LINEAR.

$CONTRL OPTIONS
---------------
SCFTYP=RHF RUNTYP=SURFACE EXETYP=RUN
MPLEVL= 0 LOCAL =NONE UNITS =ANGS
MULT = 1 ICHARG= 0 MAXIT = 30
NPRINT= 7 IREST = 0 COORD =UNIQUE
ECP =NONE NORMF = 0 NORMP = 0
ITOL = 20 ICUT = 9 NZVAR = 0
NOSYM = 0 INTTYP=POPLE GEOM =INPUT
PLTORB= F MOLPLT= F RPAC = F
AIMPAC= F FRIEND= CITYP =GUGA

$SYSTEM OPTIONS
---------------
KDIAG = 0 MEMORY= 2000000 TIMLIM= 36000.0 SEC.
COREFL= F PTIME = F R = F
BALTYP=NXTVAL

----------------
PROPERTIES INPUT
----------------

MOMENTS FIELD POTENTIAL DENSITY
IEMOM = 1 IEFLD = 0 IEPOT = 0 IEDEN = 0
WHERE =COMASS WHERE =NUCLEI WHERE =NUCLEI WHERE =NUCLEI
OUTPUT=BOTH OUTPUT=BOTH OUTPUT=BOTH OUTPUT=BOTH
IEMINT= 0 IEFINT= 0 IEDINT= 0
MORB = 0

EXTRAPOLATION IN EFFECT

----------------------
INTEGRAL INPUT OPTIONS
----------------------
NOPK = 1 NORDER= 0 SCHWRZ= T

ATTENTION! AO INTEGRALS WILL BE PACKED.
THRESHOLD FOR PACKING PKTHR = 0.10000000D-01

-------------------------------
INTEGRAL TRANSFORMATION OPTIONS
-------------------------------
NWORD = 0 CUTOFF = 1.0E-09
MPTRAN = 0 DIRTRF = T
AOINTS =DUP IREST = 0

------------------------------------------
THE POINT GROUP IS C1 , NAXIS= 0, ORDER= 1
------------------------------------------

DIMENSIONS OF THE SYMMETRY SUBSPACES ARE
A = 4

..... DONE SETTING UP THE RUN .....

CPU TIME: STEP = 0.02 , TOTAL = 0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.36 , TOTAL = 0.4 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 4.34%, TOTAL = 13.02%

---------------------------
POTENTIAL SURFACE MAP INPUT
---------------------------
COORD 1 LYING ALONG ATOM PAIR 1 2
HAS ORIGIN=-0.400, DISPLACEMENT= 0.100 AND 40 STEPS.
GROUP 1 CONTAINS 1 ATOMS:
2
...... END OF ONE-ELECTRON INTEGRALS ......

CPU TIME: STEP = 0.00 , TOTAL = 0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.01 , TOTAL = 0.4 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 12.67%

-------------
GUESS OPTIONS
-------------
GUESS =HUCKEL NORB = 0 NORDER= 0
MIX = F PRTMO = F SYMDEN= F
TOLZ = 1.0E-08 TOLE = 1.0E-05

INITIAL GUESS ORBITALS GENERATED BY HUCKEL ROUTINE.
HUCKEL GUESS REQUIRES 1412 WORDS.

SYMMETRIES FOR INITIAL GUESS ORBITALS FOLLOW. BOTH SET(S).
1 ORBITALS ARE OCCUPIED ( 0 CORE ORBITALS).
1=A 2=A 3=A 4=A
...... END OF INITIAL ORBITAL SELECTION ......

CPU TIME: STEP = 0.00 , TOTAL = 0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.07 , TOTAL = 0.4 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 10.65%
...... END OF TWO-ELECTRON INTEGRALS .....

CPU TIME: STEP = 0.00 , TOTAL = 0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.01 , TOTAL = 0.4 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 10.42%

DIRECT SCF CALCULATION, SCHWRZ=T FDIFF=T

NONZERO BLOCKS
ITER EX DEM TOTAL ENERGY E CHANGE DENSITY CHANGE DIIS ERROR INTEGRALS SKIPPED
1 0 0 -0.615477086 -0.615477086 0.083273977 0.000000000 55 0
2 1 0 -0.644522586 -0.029045499 0.012699048 0.000000000 55 0
3 2 0 -0.645138919 -0.000616334 0.002015003 0.000000000 55 0
4 3 0 -0.645154260 -0.000015340 0.000321835 0.000000000 55 0
5 4 0 -0.645154650 -0.000000391 0.000051458 0.000000000 55 0
6 5 0 -0.645154660 -0.000000010 0.000008229 0.000000000 55 0
7 6 0 -0.645154661 0.000000000 0.000001316 0.000000000 55 0

-----------------
DENSITY CONVERGED
-----------------
TIME TO FORM FOCK OPERATORS= 0.0 SECONDS ( 0.0 SEC/ITER)
FOCK TIME ON FIRST ITERATION= 0.0, LAST ITERATION= 0.0
TIME TO SOLVE SCF EQUATIONS= 0.0 SECONDS ( 0.0 SEC/ITER)

FINAL ENERGY IS -0.6451546607 AFTER 7 ITERATIONS
...... END OF RHF CALCULATION ......

CPU TIME: STEP = 0.00 , TOTAL = 0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.06 , TOTAL = 0.5 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 0.00%, TOTAL = 9.19%

GUGA CI RUN OPTIONS NRNFG NPFLG
-------------------------------------

-DRT- TABLE 1 -5
TRANSFORMATION 1 -5
ENERGY MATRIX 1 -5
DIAGONALIZATION 1 -5
1E-DENSITY MATRIX 1 -5
2E-DENSITY MATRIX 0 -5
LAGRANGIAN MATRIX 0 -5
POINT GROUP= IS UNKNOWN TO THE DRT CODE.


EXECUTION OF GAMESS TERMINATED ABNORMALLY AT 18:18:22 LT 30-OCT-2008
100000 WORDS OF DYNAMIC MEMORY USED

CPU TIME: STEP = 0.00 , TOTAL = 0.0 SECONDS ( 0.0 MIN)
WALL CLOCK TIME: STEP = 0.00 , TOTAL = 0.5 SECONDS ( 0.0 MIN)
CPU UTILIZATION: STEP = 100.00%, TOTAL = 9.19%
STOP IN ABRT

DanilaZ

PC GAMESS version 6.2
Дело видимо в этом. Попробуй $CIDRT:GROUP=C1

sidy_mos

Работает, спасибо! =)

sidy_mos

Народ, а информация о том, какая орбиталь занята содержится тут:
 
$cidrt
fors=.t.
nfzc=1
ndoc=2
nval=1
GROUP=C1
$end

А то я даже с почитыванием мануала что-то не совсем разберусь как с заполненными орбиталями играться и подбирать нужные конфигурации. Подскажите пожалуйста =)

terl

какая орбиталь занята содержится тут
не орбиталь, а номер возбужденного состояния задается словом IROOT в группе $GUGDM
вообще, судя по вопросам, тебе стоит прочитать про метод конфигурационного взаимодействия.
Оставить комментарий
Имя или ник:
Комментарий: