detcas – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

detcas – Determinante Konfigurationsinteraktion / Komplett aktiv
Space Self-Consistent-Field-Programm

BESCHREIBUNG

Das Programm detcas führt ein Complete-Active-Space-Self-Consistent-Field (CASSCF) durch
Optimierung von Molekülorbitalen über ein zweistufiges Verfahren in Verbindung mit dem
Determinantenkonfigurationsinteraktionsprogramm DETCI. Das Programm ist ziemlich einfach und
verwendet derzeit einen Newton-Raphson-Ansatz, um die Orbitale mithilfe eines einfachen Verfahrens zu aktualisieren
ungefähre orbitale Hessische. Die Konvergenz wird durch die direkte Umkehrung von Pulay beschleunigt
das iterative Unterraumverfahren (DIIS). Der Code wurde geschrieben, um eine allgemeinere Darstellung zu ermöglichen
Wellenfunktionen, die nicht unbedingt eine vollständige CI-Behandlung des aktiven Raums umfassen. In
Insbesondere kann jede von DETCI unterstützte CI-Wellenfunktion mit eingeschränktem aktivem Raum (RAS) unterstützt werden
verwendet, was MCSCF-Wellenfunktionen vom Typ RASSCF ermöglicht.

REFERENZEN

Ungefährer orbitaler Hessischer:

1. G. Chaban, MW Schmidt und MS Gordon, Theor. Chim. Acta 97, 88-95 (1997).

Eingeschränkter aktiver Raum CI:

1. Determinantenbasierte Konfigurationsinteraktionsalgorithmen für vollständige und eingeschränkte Konfigurationen
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen und HJ Aa.
Jensen, J. Chem. Physik. 89, 2185 (1988).

Eingeschränkter aktiver Raum SCF:

1. P.-A. Malmqvist, A. Rendell und BO Roos, J. Phys. Chem. 94, 5477 (1990).

DETCI-Programm:

1. CD Sherrill, Computational Algorithms for Large-Scale Full and Multi-Reference
Konfigurationsinteraktionswellenfunktionen, Doktorarbeit, University of Georgia, Athens,
GA, 1996.

2. CD Sherrill und HF Schaefer, The Configuration Interaction Method: Advances
in Highly Corlated Approaches, Adv. Quantenchemie. 34, 143-269 (1999).

DATEIEN ERFORDERLICH

input.dat – Eingabedatei
Datei78 – MO-Einelektronenintegrale (fzc-Operator)
Datei72 – MO-Zwei-Elektronen-Integrale
Datei73 – MO-Einteilchen-Dichtematrix
Datei74 – MO-Zwei-Teilchen-Dichtematrix
Datei 75 – MO Lagrange

DATEIEN AKTUALISIERT

Output.dat – Ausgabedatei
file14.dat – Aufzeichnung von Energien und Orbitalgradienten

SPEISUNG FORMAT

Die folgenden Befehlszeilenargumente sind verfügbar:

-ruhig Dies ergibt das gleiche Ergebnis wie DRUCKEN=0.

-o Fanname
Gibt den Dateinamen für die Ausgabedatei an. Der Standardwert ist „output.dat“.

Zusätzliche Eingaben für dieses Programm werden aus der Datei input.dat gelesen. Je häufiger verwendet
Schlüsselwörter sind:

WFNMehr = Schnur
Die unterstützten Wellenfunktionstypen sind CASSCF und RASSCF.

KONVERGENZ = ganze Zahl
Konvergenz auf dem Orbitalgradienten gewünscht. Konvergenz wird erreicht, wenn der RMS
der Fehler im Orbitalgradienten ist kleiner als 10**(-n). Der Standardwert ist 4 für
Energieberechnungen und 7 für Gradienten.

DOCC = integer_array
Dieser Vektor gibt die Anzahl der doppelt besetzten Orbitale in jedem Irrep an. Es gibt
kein Standard.

SOCC = integer_array
Dieser Vektor gibt die Anzahl der einfach besetzten Orbitale in jedem Irrep an. Es gibt
kein Standard.

ENERGIE_KONVERGENZ = ganze Zahl
Gewünschte Konvergenz der gesamten MCSCF-Energie. Der Standardwert ist 7.

FROZEN_DOCC = integer_array
Die Anzahl der doppelt besetzten Orbitale mit der niedrigsten Energie in jedem irreduziblen Orbital
Darstellung, die buchstäblich eingefroren wird (im MCSCF nicht aktualisiert). Der
Es wird eine Baumwollreihenfolge der irreduziblen Darstellungen verwendet. Der Standardwert ist
Nullvektor.

FROZEN_UOCC = integer_vector
Die Anzahl der unbesetzten Orbitale mit der höchsten Energie in jeder irreduziblen Darstellung
Dies wird buchstäblich eingefroren (im MCSCF nicht aktualisiert). Der Standardwert ist Null
Vektor.

RESTRICTED_DOCC = integer_array
Die Anzahl der doppelt besetzten Orbitale mit der niedrigsten Energie in jedem irreduziblen Orbital
Darstellung, die optimiert, aber im MCSCF doppelt besetzt bleibt.
Diese Orbitale folgen dem FROZEN_DOCC Orbitale. Der Standardwert ist Null
Vektor.

RESTRICTED_UOCC = integer_array
Die Anzahl der unbesetzten Orbitale mit der höchsten Energie in jeder irreduziblen Darstellung
die optimiert, aber im MCSCF unbesetzt bleiben. Diese Orbitale kommen
Vor dem FROZEN_UOCC Orbitale. Der Standardwert ist der Nullvektor.

NCASITER = ganze Zahl
Maximale Anzahl von Iterationen zur Optimierung der Orbitale. Diese Option sollte be
angegeben in DEFAULT Abschnitt of Eingangsmöglichkeiten: , weil es für die sichtbar sein muss
Steuerprogramm PSI. Der Standardwert ist 1.

IHR EIGENES LOGO = ganze Zahl
Diese Option bestimmt die Ausführlichkeit der Ausgabe. Ein Wert von 1 oder 2 gibt an
minimales Drucken, ein Wert von 3 gibt ausführliches Drucken an. Werte von 4 oder 5 sind
Wird zum Debuggen verwendet. Verwenden Sie Level 5 nicht, es sei denn, der Testfall ist sehr klein (z. B
STO H2O CISD).

Die weniger häufig verwendeten Schlüsselwörter sind:

DIIS_FREQ = ganze Zahl
Das DIIS-Extrapolationsverfahren wird jedes Mal versucht n Iterationen. Der Standard
ist 1.

DIIS_MAX_VECS = ganze Zahl
Die maximale Anzahl von Unterraumvektoren für die DIIS-Prozedur. Nach dieser Nummer
Anzahl der Vektoren erreicht ist, werden ältere Vektoren aus dem Unterraum gelöscht
notwendig. Der Standardwert ist 8.

DIIS_MIN_VECS = ganze Zahl
Die Mindestanzahl von Unterraumvektoren, bevor eine DIIS-Interpolation durchgeführt werden kann
durchgeführt. Der Standardwert ist 2. Werte unter 2 sind nicht sinnvoll.

DIIS_START = ganze Zahl
Die Iterationsnummer bei Pulays direkter Inversion des iterativen Unterraums (DIIS)
Das Verfahren zur Konvergenzbeschleunigung sollte eingeschaltet sein. Der Standardwert ist 3.
Vor dieser Iteration werden dem DIIS-Unterraum keine Vektoren hinzugefügt. Der erste
Der DIIS-Schritt wird erst unternommen DIIS_MIN_VECS Vektoren liegen im DIIS-Unterraum.

LAG_DATEI = ganze Zahl
Datei (Einheitsnummer) zum Lesen der Lagrange-Matrix. Der Standardwert ist
derzeit 75.

OPDM_FILE = ganze Zahl
Datei (Einheitsnummer) zum Auslesen der Ein-Partikel-Dichtematrix. Der Standardwert
liegt derzeit bei 73.

SCALE_GRAD = boolean
Gibt an, ob der Orbitalgradient durch das ungefähre (diagonale) Orbital skaliert werden soll
Hessisch. Der Standardwert ist TRUE.

SCALE_STEP = echt
Skalierungsfaktor für den Orbitalrotationsschritt. Der Standardwert ist 1.0.

TPDM_FILE = ganze Zahl
Datei (Einheitsnummer) zum Auslesen der Zwei-Partikel-Dichtematrix. Der Standardwert
liegt derzeit bei 74.

LEVEL_SHIFT = boolean
Gibt an, ob eine Pegelverschiebung der Hesse-Matrix zugelassen werden soll. Dies kann genutzt werden
um sicherzustellen, dass der Hesse-Wert für die Anfangsiterationen positiv definit ist. Wenn
Ist die Hesse-Kurve nicht positiv definitiv nahe der Konvergenz, dann kann dies ein Zeichen dafür sein
numerische Instabilitäten im MCSCF. Der Standardwert ist TRUE.

SHIFT = echt
Gibt an, um wie viel die diagonalen Elemente des Hessischen für die Ebene verschoben werden sollen
Verschiebung. Der Standardwert ist 0.01.

DETERM_MIN = echt
Minimal zulässiger Wert für die Determinante der Hesse-Matrix bei Ebenenverschiebung
ist eingeschaltet. Dies wird nicht verwendet, wenn LEVEL_SHIFT=FALSE. Der Standardwert ist 0.00001.

STEP_MAX = echt
Dies ist die maximal zulässige einzelne Orbitalrotation. Der Standardwert ist 0.30.

FORCE_STEP = boolean
Dadurch kann der Benutzer den berechneten Schritt überschreiben und einen Schritt in a erzwingen
bestimmte Richtung. Dies kann nützlich sein, wenn Sie versuchen, die Berechnung zu erzwingen
von einem Sattelpunkt aus. Der Standardwert ist FALSE.

FORCE_PAIR = ganze Zahl
Dies ist der Index für das unabhängige Paar, das gedreht werden soll, wenn
FORCE_STEP=TRUE. Dies wird ansonsten ignoriert. Es gibt keinen Standardwert.

FORCE_VALUE = echt
Dies ist der Orbitalrotationswert zum Drehen der Orbitale, der durch angegeben wird
FORCE_PAIR. Dies wird nur verwendet, wenn FORCE_STEP auf TRUE gesetzt ist. Der Standardwert ist
0.0

CHECK_HESSIAN = boolean
Berechnen Sie die Eigenwerte der Orbital-Hesse-Kurve. Dies eignet sich gut zur Überprüfung
Sattelpunktlösungen. Sollte ohne Pegelverschiebung verwendet werden. Die Standardeinstellung ist
FALSCH.

EIGEN_VECTORS = boolean
Möchten Sie auch die Eigenvektoren des Orbital-Hesses berechnen? Das ist
Wird ignoriert, wenn CHECK_HESSIAN=FALSE. Der Standardwert ist FALSE.

8 Mai, 1998 detcas(1)

Nutzen Sie Detcas online über die Dienste von onworks.net