|
|
|
CI: Configuration Interaction, Konfigurationswechselwirkung, Konfigurationsmischung
Theoretischer Hintergrund
Beispiel H2:
STO-6G, S0
!
! H2-Molekuel, CI
!
$CONTRL SCFTYP=RHF MULT=1 RUNTYP=ENERGY COORD=ZMT
CITYP=GUGA
$END
$SYSTEM TIMLIM=10 MEMORY=1000000 $END
$BASIS GBASIS=STO NGAUSS=6 $END
$CIDRT GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=1 $END
$GUGDIA NSTATE=2 $END
$GUGDM IROOT=1 NFLGDM(1)=1 $END
$GUGDM2 WSTATE(1)=1.0,0.0 $END
$GUESS GUESS=HUCKEL $END
$DATA
H2-Molekuel, CI
DNH 4
H
H 1 R
R 0.74
$END
Die folgenden Befehlsgruppen sind für eine MCSCF/CI-Rechnung notwendig (rosa markiert) bzw. empfehlenswert. Falls nur eine CI-Rechnung durchgeführt werden soll (keine MCSCF), werden die dunkelgrün markierten Befehle und Schlüsselwörter weggelassen. Bei einer CI-Rechnung (CITYP=GUGA oder ) kann RUNTYP=OPTIMIZE nur für SCFTYP=RHF angewendet werden (auch nicht kombiniert mit MCSCF). Man muss dann auf RUNTYP=TRUDGE ausweichen, was aber COORD=HINT voraussetzt.
$CONTRL SCFTYP=MCSCF MULT=2 RUNTYP=ENERGY CITYP=GUGA $END
$SYSTEM TIMLIM=1000 MEMORY=5000000 $END
$SCF NCONV=6 $END
$MCSCF MAXIT=500 MICIT=2 $END
$DET NCORE=1 NACT=10 NELS=9 NSTATE=2 $END
$DRT GROUP=C1 FORS=.TRUE. NMCC=1 NDOC=4 NALP=1 NVAL=5 $END
$CIDRT GROUP=C1 FORS=.TRUE. NFZC=1 NDOC=4 NALP=1 NVAL=5 $END
$GUGDIA NSTATE=2 $END
$GUGDM IROOT=2 NFLGDM(2)=1 $END
$GUESS GUESS=MOREAD NORB=11 $END
$DRT: GROUP=C1: Punktgruppe. MCSCF bzw. CI kann nur mit folgenden Punktgruppen arbeiten: C1, C2, CI, CS, C2V, C2H, D2, D2H, C4V, D4, D4H.
Wenn die Punktgruppe nicht angegeben wird, wird die aus der $DATA-Eingabe verwendet.
FORS=.TRUE.: Full Optimized Reaction Space ist identisch mit CAS: Complete Active Space und bedeutet, dass alle Kombinationen der ausgewählten (aktiven) Elektronen und Orbitale berücksichtigt werden. Das funktioniert in der Regel nur bei der Verwendung einer kleinen Zahl von Orbitalen und hängt auch von der Größe des Basissatzes ab (bei kleinem Basissatz kann man CAS größer machen).Will man die Zahl der aktiven Orbitale größer machen, muss man die Anzahl der aktiven Elektronen mit IEXCIT=N angeben und FORS=.FALSE. setzen bzw. weglassen (FORSE=.FALSE. ist default).
$DRT GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=3 $END
$CIDRT GROUP=C1 IEXCIT=2 NFZC=0 NDOC=1 NVAL=3 $END
NCORE=1: Zahl der eingefrorenen Orbitale
NACT=10: Zahl der aktiven Orbitale
NELS=1: Zahl der aktiven Elektronen
NFZC=1: Zahl der eingefrorenen Orbitale
NDOC=4: Zahl der doppelt-besetzten Orbitale
NALP=1: Zahl der besetzten alpha-Orbitale
NVAL=5: Zahl der unbesetzten Orbitale
NSTATE=2: Zahl der Elektronenzustände (hier: Grundzustand und erster angeregter Zustand.
IROOT=2: Nummer des Elektronenzustandes dessen Eigenschaften berechnet werden.
NFLGDM(2)=1: Die Zahl in der Klammer muss mit der von IROOT identisch sein. =1 bedeutet, dass die Dichtematrix des Zustandes 2 ausgegeben wird.
Input-Files für MCSCF/CI-Rechnungen
Output-Files