cscf – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

cscf – löst die Hartree-Fock-Gleichungen

BESCHREIBUNG

Das Programm cscf führt das iterative Verfahren zur Lösung der Hartree-Fock-Gleichungen durch.

Dieses Programm ist auf die D2h-Symmetrie und ihre Untergruppen sowie die Orbitalbesetzungen beschränkt
müssen ganze Zahlen sein. Damit lassen sich bestimmte reine Drehimpulszustände ableiten
Eine teilweise Besetzung entarteter Orbitale kann mit den vorliegenden Codes nicht erreicht werden. Für
Beispiel: Der 2PIu-Zustand (Dublett PI u) von linearem ONO, abgeleitet von der niedrigsten Energie
Die lineare (pi u)1-Konfiguration kann nur als 2B2u (Dublett B2u) oder 2B3u berechnet werden
(Dublett B 3u)-Komponente des 2PIu (Dublett PI u)-Zustands und die daraus resultierende räumliche
Die Wellenfunktion weist keine PI-Symmetrie auf. In gewissem Sinne ist dies jedoch wünschenswert,
da die Energie eine kontinuierliche Funktion des Biegewinkels ist. Berechnung der Energie
von gebogenen Konfigurationen wie 2B2u (Dublett B 2u) oder 2B3u (Dublett B 3u) und einem reinen 2PIu
(Dublett PI u)-Zustand bei linearen Geometrien führt zu einer ausgeprägten Diskontinuität.

Zumeist entstehen Triplettzustände durch Doppelbesetzung eines doppelt entarteten Zustands
Orbital, wie der 3A2-Zustand (Triplett A 2), der sich aus (e')2 oder (e")2 ergibt
Konfigurationen in D3h-Symmetrie oder der 3SIGMAg-Zustand (Triplett-SIGMA g) eines (pi g)2 oder (pi
u)2-Konfiguration in Dinfh-Symmetrie (D unendlich h) hat die richtige räumliche Symmetrie.
Die aus denselben elektronischen Konfigurationen resultierenden Singulettzustände sind inhärent
Multikonfigurationen und werden daher nicht gut durch Einzelkonfigurationen repräsentiert
Wellenfunktionen.

REFERENZEN

PK-Datei-Methode:

1. RC Raffenetti, Chem. Physik. Lette. 20 (1973) 335.

Molekulare Symmetrie und HF-Berechnungen mit geschlossener Schale:

1. M. Dupuis und HFKing, Int. J. Quant. Chem. 11 (1977) 613.

DIIS für geschlossene Shell:

1. P. Pulay, Chem. Physik. Lette. 73 (1980) 393.

2. P. Pulay, J. Comp. Chem. 3 (1982) 556.

Kopplungskoeffizienten (Alpha und Beta) für offene Schale:

1. CCJ Roothaan, Rev. Mod. Physik. 32 (1960) 179.

Dämpfung:

1. DR Hartree, „The Calculation of Atomic Structures“ (Wiley: New York) 1957.

2. MC Zerner und M. Hehenberger, Chem. Physik. Lette. 62 (1979) 550.

Pegelverschiebung:

1. VR Saunders und IH Hillier, Int. J. Quant. Chem. 7 (1973) 699.

KONVERGIEREND CSCF

Für schwierige Fälle mit offener Schale empfiehlt sich die Verwendung einer entsprechenden geschlossenen Schale
Zuerst muss die Berechnung durchgeführt werden (ein zusätzliches Elektron hinzufügen oder entfernen) und dann dieser SCF-Vektor
als Schätzwert für die gewünschte Open-Shell-Wellenfunktion verwendet werden. In TCSCF-Fällen ist dies immer der Fall
Es ist ratsam, zuerst eine geschlossene Shell (oder vielleicht das entsprechende Triplet) SCF auszuführen und dann zu verwenden
Dies ist eine Schätzung für die TCSCF.

Für Open-Shell-Systeme wird ein Pegelverschiebungswert von 0.5 bis 3.0 empfohlen. Beginnen Sie mit a
hohen Wert (2.0 - 3.0) für die erste SCF-Berechnung und reduzieren Sie ihn dann (auf 0.5 - 1.0) für
Nachfolgende Läufe, die einen konvergenten SCF-Vektor als Ausgangspunkt verwenden.

Es ist äußerst wichtig zu beachten, dass diese Version des Codes nicht mehr unterstützt wird
OFFENER TYP. Eins sollen - neu Schlüsselwörter REFERENCE und MULTP zu angeben tippe of SCF
Installation!

SPEISUNG FORMAT

Das cscf Das Programm durchsucht den Standard-Schlüsselwortpfad (zuerst SCF und dann DEFAULT)
für die folgenden Schlüsselwörter:

LABEL = Schnur
Dies ist eine Zeichenfolge, die in die Ausgabe einbezogen werden soll. Diese Zeichenfolge wird nicht verwendet
durch das Programm. Es gibt keinen Standardwert.

WFNMehr = Schnur
Dies ist die Art von Wellenfunktion, die letztendlich gewünscht wird. Die Standardeinstellung ist SCF.

OFFENER TYP is nicht länger unterstützt

REFERENCE = Schnur
Dies gibt die Art der SCF-Berechnung an, die durchgeführt werden soll. Es kann eines davon sein RHF
(für ein Singlet mit geschlossener Schale), ROHF (für eine eingeschränkte Open-Shell-Berechnung), UHF
(für eine uneingeschränkte Berechnung mit offener Schale), TWOCON (für eine Zwei-Konfiguration
Singulett), oder SPECIAL. Wenn SPECIAL gegeben ist, dann Alpha- und Beta-Kopplung
Koeffizienten müssen mit angegeben werden ALPHA und BETA Schlüsselwörter. Die Standardeinstellung ist RHF.

MULTP= ganze Zahl
Gibt die Multiplizität des Moleküls an. Die Standardeinstellung ist Singlet.

LADUNG= ganze Zahl
Gibt die Ladung des Moleküls an. Standardmäßig ist 0.

DOCC = integer_vector
Dies gibt die Anzahl der doppelt besetzten Orbitale in jedem Irreduziblen an
Darstellung. Es gibt keinen Standardwert. Wenn dies nicht der Fall ist, wird CSCF versuchen, dies zu tun
Erraten Sie die Berufe mithilfe des Kern-Hamiltonian.

SOCC = integer_vector
Dies gibt die Anzahl der einfach besetzten Orbitale in jedem Irreduziblen an
Darstellung. Es gibt keinen Standardwert.

DERTYPE = Schnur
Dies gibt die Reihenfolge der Ableitung an, die letztendlich durchgeführt werden soll. Es ist benutzt
von dem scf Programm, um festzustellen, ob bestimmte Dateien geschrieben werden sollen, und das ist auch der Fall
Wird verwendet, um die Standardkonvergenz der Wellenfunktion zu bestimmen. Die Standardeinstellung ist
ZUERST.

MAXITER = ganze Zahl
Dies ergibt die maximale Anzahl von Iterationen. Der Standardwert ist 40.

KONVERGENZ = ganze Zahl
Dies gibt an, wie stark die Wellenfunktion konvergiert. Konvergenz ist
bestimmt durch Vergleich der RMS-Änderung in der Dichtematrix („Delta P“) mit der
gegebener Wert. Das Konvergenzkriterium beträgt 10**(-ganze Zahl). Der Standardwert ist 7, wenn
beide DERTYPE = NONE und WFNMehr = SCF gegeben sind und 10 sonst.

LEVELSHIFT = echt
Dies gibt die Pegelverschiebung an. Der Standardwert ist 1.

DIREKT = boolean
Gibt an, ob die SCF-Berechnung mit einer integralen direkten Technik durchgeführt werden soll. Der
Standard ist falsch.

PRINT_MOS = boolean
Gibt an, ob die Molekülorbitale gedruckt werden sollen oder nicht. Der Standardwert ist falsch.

Es gibt auch eine große Anzahl weniger häufig verwendeter Eingabeparameter. Wenn Sie nicht
Verstehen Sie, was die folgenden Optionen bedeuten, und stellen Sie dann sicher, dass sie nicht in Ihrem angezeigt werden
Eingang. Die Standardeinstellungen funktionieren in den allermeisten Fällen. Diese sind angegeben
mit folgenden Schlüsselwörtern:

DELETE_INTS = boolean
Integraldateien werden gelöscht, wenn WFNMehr = SCF und DERTYPE = ZUERST or DERTYPE = NONE.
Wenn Sie Integraldateien behalten möchten, legen Sie fest DELETE_INTS = falsch. Die Standardeinstellung ist
wahr.

NACHBESTELLUNG = Schnur
Der Parameter steuert die Neuordnung der Molekülorbitale. Wenn eingestellt auf VOR dann ist die
Die geschätzten Orbitale aus der Checkpoint-Datei werden neu angeordnet. Wenn eingestellt auf AFTER, konvergierte
Orbitale werden neu angeordnet, bevor sie in die Prüfpunktdatei geschrieben werden. In beiden
Häuser MOORDER Der Parameter muss angegeben werden, um die Neuordnungskarte anzugeben. Die Standardeinstellung ist
Orbitale nicht neu anordnen.

MOORDER = integer_vector
Dies gibt einen Molekülorbital-Neuordnungsvektor an. Es wird nur verwendet, wenn
NACHBESTELLUNG eingestellt ist. Dieser Vektor ordnet jedes Orbital seinem neuen Index zu, z MOORDER = (0
2 1) Gibt an, dass nach der Neuordnung die Orbitale 1 und 2 vertauscht werden. Der Rang von
Dieser Vektor entspricht der Anzahl der MOs. Die Indizes sind in der Pitzer-Reihenfolge angeordnet
(geordnet nach Symmetrie, dann nach Energie innerhalb jedes Symmetrieblocks), Basis-0. CSCF
wird wahrscheinlich scheitern, wenn der gegebene MOORDER Orbitale aus verschiedenen Irreps mischt. Dort
ist kein Standard.

ALPHA = real_vector
If OFFENER TYP = SPECIAL, dann gibt dieser Parameter die Alpha-Kopplungskoeffizienten an.
Die Anzahl der Elemente in diesem Vektor beträgt MM(MM+1)/2, wobei MM die Anzahl der Elemente ist
irreduzible Darstellungen mit einfach besetzten Molekülorbitalen. Dort
ist kein Standard.

BETA = real_vector
If OFFENER TYP = SPECIAL, dann gibt dieser Parameter die Beta-Kopplungskoeffizienten an.
Die Anzahl der Elemente in diesem Vektor beträgt MM(MM+1)/2, wobei MM die Anzahl der Elemente ist
irreduzible Darstellungen mit einfach besetzten Molekülorbitalen. Dort
ist kein Standard.

GUESS = Schnur
Diese Option bestimmt die Art der anfänglichen Schätzung des Eigenvektors, die CSCF verwenden wird.
Die einzig gültigen Optionen sind derzeit: (1) GUESS = CORE, was dazu führt, dass es verwendet wird
Kern-Hamilton-Eigenvektor, um die Berechnung zu starten; (2) GUESS = AUTO welche
führt zu einem Versuch, den MO-Vektor in der Prüfpunktdatei zu verwenden oder darauf zurückzugreifen
Kernschätzung, ob in dieser Datei kein Eigenvektor vorhanden ist. Der Standardwert, wenn AUTO.

IPRINT = ganze Zahl
Dies ist eine Druckoption. Der Standardwert ist 0.

MO_OUT = boolean
Druckt die Orbitale mit Symmetrie und Besetzungen am Ende aus
Berechnung. Der Standardwert ist wahr.

DREHEN = boolean
Wenn dies falsch ist, werden die Molekülorbitale nicht gedreht. Nur die Rotation
wirkt sich auf die virtuellen Orbitale für Systeme mit offener Schale aus. Dieser Parameter muss wahr sein
für korrelierte Gradienten und muss für zweite und höhere Ableitungen falsch sein.
Der Standardwert ist false, wenn WFNMehr = SCF und wahr sonst.

CHECK_ROT = boolean
Überprüfen Sie die oben beschriebene Molekülorbitalrotation, um sicherzustellen, dass keine Säulen entstehen
Die SCF-Eigenvektormatrix wird durch die Drehung vertauscht. Hat keine Auswirkung, wenn DREHEN =
falsch. Der Standardwert ist wahr.

CHECK_MO_ORTHOGONALITÄT = boolean
Überprüfen Sie, ob die Molekülorbitale orthonormal sind. Nur zum Debuggen nützlich. Der
Standard ist falsch.

DIIS = boolean
Dies bestimmt, ob diis verwendet wird. Der Standardwert ist wahr.

DISSTART = ganze Zahl
Dies ergibt die erste Iteration, für die DIIS verwendet wird. Der Standardwert ist 0.

NDIIS = ganze Zahl
Dies gibt die Anzahl der Fehlermatrizen an, die in der diis-Prozedur verwendet werden sollen. Der Standard
ist 6 für geschlossene Shell, 4 für offene Shell und 3 für tcscf.

DIISDAMP = echt
Dies ergibt den Dämpfungsfaktor für das diis-Verfahren. Der Standardwert ist 0.0 für
geschlossene Shell, 0.02 für offene Shell und 0.01 für tcscf.

ERHÖHUNG = echt
Dies wird in tcscf verwendet, um zu bestimmen, wie oft die ci-Koeffizienten neu berechnet werden.
Eine kleine Zahl (~0.25) führt dazu, dass sie fast bei jedem scf neu berechnet werden
Wiederholung. Der Standardwert ist 0.25.

DYN_ACC = boolean
Bei der Durchführung einer direkten SCF gibt dies an, ob die dynamische Integralgenauigkeit abschneidet
verwendet wird. Der Standardwert ist „true“ (dynamische Grenzwerte verwenden). Erste Iterationen sind
mit auf sechs Stellen genauen Integralen durchgeführt. Nachdem die Dichte konvergiert ist
Wenn 10^-5 oder 30 Iterationen abgeschlossen sind, wird die volle Integralgenauigkeit verwendet. Wenn scf
Durch diese Einstellung treten Konvergenzprobleme auf, die dazu führen, dass dynamische Abschaltungen deaktiviert werden
Variable auf false könnte helfen.

ORTHOG_ONLY = boolean
Manchmal ist dies bei CASSCF oder anderen Nicht-HF/KS-Systemen zur Orbitaloptimierung der Fall
nützlich, um MOs aus anderen Geometrien für die aktuelle Geometrie neu zu orthogonalisieren
Sie können als erste Schätzung für die neuen MOs verwendet werden. Dies kann durchgeführt werden durch
Ausführen von CSCF mit ORTHOG_ONLY = was immer dies auch sein sollte.. Nachdem die Orbitale orthogonalisiert sind, wird die
Das Programm wird beendet, ohne eine SCF-Berechnung durchzuführen. Dieses Schlüsselwort wird sein
Wird ignoriert, wenn in der Prüfpunktdatei keine vorherigen Orbitale vorhanden sind. Standardmäßig ist was immer dies auch sein sollte.
if WFNMehr = DETCAS.

30 Mai, 1991 cscf(1)

Verwenden Sie cscf online über die Dienste von onworks.net