Die Vorgehensweise bei der Erstellung des Input-Files (d. h. einer Liste von Befehlen) stammt noch aus der Zeit der Lochkarteneingabe. Das wird deutlich, wenn man sich den Output-File anschaut. Dort wird das Einlesen des Input-File entsprechend dokumentiert.
Der Text des folgenden Input-Files ist zur besseren Unterscheidung von eingefügten Erläuterungen dunkelgrün dargestellt.
Eine Zeile, die in der ersten Spalte ein Ausrufungszeichen (oder irgendein anderes Zeichen) enthält, wird als Befehlszeile ignoriert.
! H2-Molekuel, excited for
PES
!
! Fuer COORD=UNIQUE muessen
die Koordinaten (x,y,z)
! der Symmetrie-Operation
entsprechen, d. h. die
! Koordinaten muessen die
gleichen Symmetrieeigenschaften
! wie das Molekuel
besitzen. Fuer lineare Molekuele entsprechend
! der verwendeten AOs.
!
! Lit.:R=0.741A, De=4.75 eV,
D0=36118 cm-1
!
Jede Befehlsreihe ("Lochkarte")
beginnt mit einem Leerzeichen und einem Befehl, welchem ein $-Zeichen
vorangestellt ist und wird mit dem Befehl $END beendet. Groß-und
Kleinschreibung wird nicht unterschieden. Dazwischen stehen Schlüsselwörter,
die die zulösende Aufgabe spezifizieren.
$CONTRL SCFTYP=RHF
MULT=1 RUNTYP=ENERGY COORD=UNIQUE CITYP=GUGA
$END
Die obere Zeile stellt die Befehlsreihe $CONTRL dar. Sie enthält Angaben zur eigentlichen Aufgabe: SCFTYP=RHF: eine RHF-Rechnung ist durchzuführen. MULT=1: Die Spinmultiplizität des zu berechnenden Systems beträgt 1 (Singulettzustände). RUNTYP=ENERGY: Die Energie des Systems bei der vorgegeben Molekülgeometrie ist zu berechnen (single-point calculation). COORD=UNIQUE: das Koordinatensystem der Geometrieeingabe ist vom Typ UNIQUE. CITYP=GUGA: Die RHF-Rechnung ist um eine CI-Rechnung vom Typ GUGA zu erweitern.
--$CONTRL EXETYP=CHECK $END
Die obere Zeile stellt eine andere $CONTRL-Befehlsreihe dar, die aber durch die vorangestellten Bindestriche unwirksam ("auskommentiert") wurde. Sie soll den Nutzer nur anzeigen, dass man gegebenenfalls auch das Schlüsselwort EXETYP=CHECK verwenden könnte.
$SYSTEM TIMLIM=1000 MEMORY=1000000 $END
Die $SYSTEM-Befehlsreihe enthält die Schlüsselwörter TIMLIM=1000: Die Zeitbegrenzung für den Job beträgt 1000 min. MEMORY=1000000: bereitgestellte Speicherkapazität = 1000000 Words.
--$SCF DIRSCF=.TRUE. $END
-- $BASIS GBASIS=DH $END
$GUESS GUESS=HUCKEL $END
$CIDRT GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=30 $END
$GUGDIA NSTATE=5 $END
NSTATE=5: Anzahl der berechneten Elektronen-Zustände (1: Grundzustand, 5: 4.angeregter Zustand)
$GUGDM IROOT=5 NFLGDM(5)=1 $END
IROOT=5: derjenige Zustand, dessen Eigenschaften berechnet werden sollen. Die
Energie und die Eigenfunktionen des Grundzustand werden immer mitgerechnet, auch
die Energien und CI der anderen Zustände. Andere Details (Density matrix,
energy components) aber nur für den Zustand, der durch IROOT festgelegt
ist.
NFLGDM(5)=1: die 5 in Klammern muss identisch IROOT sein; =1 bedeutet, dass die
Density- und natural orbital-Matrix ausgegeben wird.
$DATA
H2-Molekuel, cc-pVQZ basis set
DNH 4
H 1.0 0.00000 0.00000 0.37000
S 3
1 82.64 0.2006D-02
2 12.41 0.15343D-01
3 2.824 0.75579D-01
S 1
1 0.7977 1.0
S 1
1 0.2581 1.0
S 1
1 0.08989 1.0
P 1
1 2.292 1.0
P 1
1 0.838 1.0
P 1
1 0.292 1.0
D 1
1 2.062 1.0
D 1
1 0.662 1.0
F 1
1 1.397 1.0
$END
Wichtig ist hier die Leerzeile vor der $END-Zeile
Seitenanfang | Input-File allgemein | Arbeiten mit GAMESS |