Startseite

 

 

MCSCF

Multi-Configuration Self-Consistent Field, CASSCF (Complete Active Space) is an example of this approach.

Man versteht darunter eine CI-Rechnung, bei der nicht nur die Koeffizienten vor den Determinanten (Elektronenkonfigurationen) optimiert werden, sondern auch noch die Koeffizienten der AO's innerhalb eines MO's (für die sonst die Koeffizienten der RHF-Funktion genommen werden).

$CONTRL  SCFTYP=MCSCF

$DET NCORE=6 NACT=10 NELS=10 $END

NCORE: Zahl der eingefrorenen Orbitale, NACT: Zahl der aktiven Orbitale, NELS: Zahl der aktiven Elektronen.

 $DRT GROUP=D2H FORS=.TRUE. NMCC=6 NDOC=5 NALP=0   VAL=5  $END

oder

 $DRT GROUP=D2H FORS=.FALSE. IEXCIT=2 NMCC=6 NDOC=5 NALP=0 NVAL=5  $END

Zu beachten ist, dass die Schlüsselwörter in $DRT und $CIDRT zum Teil unterschiedlich sind:

$CIDRT GROUP=C1 IEXCIT=2 NFZC=11 NDOC=10 NAOS=0 NBOS=0 NVAL=10 $END
! The CIDRT group does not know D6h, therefore here C1
! IEXCIT: electron excitation level, here singles and doubles
! NFZC: number of frozen core orbitals, here all sigma MOs
! NDOC: number of doubly occupied MOs, NAOS, NBOS: number of singly occupied
! alpha and beta spins, NVAL: number of empty MOs, here 5

 

GROUP: Symmetriegruppe, FORS: Full Optimized Reaction Space, NMCC: Zahl der eingefrorenen Orbitale, NDOC: Zahl der doppelt-besetzten Orbitale, NALP=0: Zahl der alpha-Orbitale, NBET=0: Zahl der beta-Orbitale, NVAL: Zahl der unbesetzten Orbitale. 
IEXCIT=2: Zahl der aktiven Elektronen für den Fall, dass FORS=.FALSE..

Für die Berechnung elektronisch angeregter Zustände müssen die beiden folgenden Befehlsgruppen angegeben werden! Die Geometrieoptimierung erfolgt aber trotz dieser Angabe für den Grundzustand!

$GUGDIA NSTATE=2 $END
$GUGDM IROOT=2 $END

NSTATE: Anzahl der Elektronenzustände, die berechnet werden sollen
IROOT=2: Elektronenzustand, dessen Eigenschaften berechnet werden sollen und falls RUNTYP=OPTIMIZE, dessen Geometrie optimiert werden soll. 


Eine MCSCF-Rechnung erfordert fast immer Start-MO's, die mit GUESS=MOREAD eingelesen werden. Deshalb sollte man zuerst eine SCFTYP=RHF- oder UHF-Rechnung durchführen und die MO's aus dem Punchfile (name.pnc) in den Input-File übernehmen. 
Häufig ist es auch notwendig, die Optimierungskriterien zu verschärfen:

 $STATPT OPTTOL=1.0E-5 $END
 $SCF NCONV=6 $END
 $MCSCF MAXIT=500 MICIT=2 ACURCY=1.0E-06 $END

MAXIT: iterations, MICIT: microiterations within a single MCSCF iteration,! ACURCY:asymmetry in the Langrangian matrix

 

aus der Dokumentation

Beispiele: H2, O2,

  Anfang der Seite verwandte Themen: