detcas - Online in de cloud

PROGRAMMA:

NAAM

detcas - Interactie determinantconfiguratie / Compleet actief
Space Self-Consistent-Field-programma

PRODUCTBESCHRIJVING

Het programma details voert een zelfconsistent veld met volledige actieve ruimte uit (CASSCF)
optimalisatie van moleculaire orbitalen via een tweestapsprocedure in combinatie met de
determinant configuratie interactieprogramma DETCI. Het programma is vrij eenvoudig en
gebruikt momenteel een Newton-Raphson-benadering om de orbitalen bij te werken, met behulp van een eenvoudige
geschatte orbitale Hessiaan. De convergentie wordt versneld met behulp van Pulay's directe inversie van
de iteratieve deelruimte (DIIS) procedure. De code is geschreven om meer algemeenheid mogelijk te maken
golffuncties die niet noodzakelijkerwijs een volledige CI-behandeling van de actieve ruimte omvatten. In
in het bijzonder kan elke CI-golffunctie met beperkte actieve ruimte (RAS) die door DETCI wordt ondersteund, dat zijn
gebruikt, waardoor MCSCF-golffuncties van het RASSCF-type mogelijk zijn.

REFERENTIES

Geschatte orbitale Hessiaan:

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

Beperkte actieve ruimte-CI:

1. Op determinanten gebaseerde configuratie-interactie-algoritmen voor volledig en beperkt
Configuratie Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen en HJ Aa.
Jensen, J. Chem. Fys. 89, 2185 (1988).

Beperkte actieve ruimte SCF:

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

DETCI-programma:

1. CD Sherrill, computeralgoritmen voor grootschalige volledige en multi-referentie
Configuratie-interactiegolffuncties, proefschrift, Universiteit van Georgia, Athene,
GA, 1996.

2. CD Sherrill en HF Schaefer, De configuratie-interactiemethode: vooruitgang
in sterk gecorreleerde benaderingen, Adv. Kwantumchemie 34, 143-269 (1999).

FILES VERPLICHT

input.dat - Invoerbestand
file78 - MO-integralen van één elektron (fzc-operator)
file72 - MO twee-elektronenintegralen
file73 - MO-dichtheidsmatrix met één deeltjes
file74 - MO-dichtheidsmatrix met twee deeltjes
file75 - MO Lagrangiaans

FILES UPDATE

output.dat - Uitvoerbestand
file14.dat - Record van energieën en orbitale gradiënten

INVOER FORMAT

De volgende opdrachtregelargumenten zijn beschikbaar:

-rustig Dit geeft hetzelfde resultaat als AFDRUKKEN=0.

-o fnaam
Geeft de bestandsnaam voor het uitvoerbestand. Standaard ingesteld op output.dat.

Extra invoer voor dit programma wordt uit het bestand input.dat gelezen. Hoe vaker gebruikt
trefwoorden zijn:

WFNMeer = snaar
De ondersteunde golffunctietypen zijn CASSCF en RASSCF.

CONVERGENTIE = geheel getal
Convergentie gewenst op de orbitale gradiënt. Convergentie wordt bereikt wanneer de RMS
van de fout in de orbitale gradiënt is minder dan 10**(-n). De standaardwaarde is 4 voor
energieberekeningen en 7 voor gradiënten.

DOCC = geheel getal_array
Deze vector geeft het aantal dubbel bezette orbitalen in elke irrep. Er bestaat
geen standaard.

socc = geheel getal_array
Deze vector geeft het aantal afzonderlijk bezette orbitalen in elke irrep. Er bestaat
geen standaard.

ENERGIE_CONVERGENTIE = geheel getal
Convergentie gewenst op de totale MCSCF-energie. De standaardwaarde is 7.

BEVROREN_DOCC = geheel getal_array
Het aantal dubbel bezette orbitalen met de laagste energie in elk onherleidbaar
representatie die letterlijk zal worden bevroren (niet bijgewerkt in de MCSCF). De
Er wordt gebruik gemaakt van een katoenen volgorde van de onherleidbare representaties. De standaardinstelling is de
nul-vector.

BEVROREN_UOCC = geheel getal_vector
Het aantal onbezette orbitalen met de hoogste energie in elke onherleidbare representatie
die letterlijk zal worden bevroren (niet bijgewerkt in de MCSCF). De standaardwaarde is nul
vector.

RESTRICTED_DOCC = geheel getal_array
Het aantal dubbel bezette orbitalen met de laagste energie in elk onherleidbaar
representatie die zal worden geoptimaliseerd maar dubbel bezet zal blijven in de MCSCF.
Deze orbitalen komen na de BEVROREN_DOCC orbitalen. De standaardwaarde is nul
vector.

RESTRICTED_UOCC = geheel getal_array
Het aantal onbezette orbitalen met de hoogste energie in elke onherleidbare representatie
die zal worden geoptimaliseerd maar onbezet blijft in de MCSCF. Deze orbitalen komen
voor de BEVROREN_UOCC orbitalen. De standaardwaarde is de nulvector.

NCASITER = geheel getal
Maximaal aantal iteraties om de orbitalen te optimaliseren. Deze optie moet be
gespecificeerd in de DEFAULT sectie of invoer, omdat het zichtbaar moet zijn voor de
besturingsprogramma PSI. Standaard ingesteld op 1.

BEDRUK = geheel getal
Deze optie bepaalt de breedsprakigheid van de uitvoer. Een waarde van 1 of 2 geeft aan
minimaal afdrukken, een waarde van 3 specificeert uitgebreid afdrukken. Waarden van 4 of 5 zijn
gebruikt voor het debuggen. Gebruik niveau 5 alleen als de testcase erg klein is (bijv
STO H2OCISD).

De minder vaak gebruikte trefwoorden zijn:

DIIS_FREQ = geheel getal
De DIIS-extrapolatieprocedure zal elke keer worden geprobeerd n iteraties. De standaard
is 1.

DIIS_MAX_VECS = geheel getal
Het maximale aantal subruimtevectoren voor de DIIS-procedure. Na dit nummer
Als het aantal vectoren wordt bereikt, worden oudere vectoren uit de subruimte verwijderd als
nodig. De standaardwaarde is 8.

DIIS_MIN_VECS = geheel getal
Het minimale aantal subruimtevectoren vóór een DIIS-interpolatie kan zijn
uitgevoerd. De standaardwaarde is 2. Het heeft geen zin om waarden kleiner dan 2 te hebben.

DIIS_START = geheel getal
Het iteratienummer wanneer Pulay's directe inversie van de iteratieve subruimte (DIIS)
procedure voor het versnellen van de convergentie moet worden ingeschakeld. De standaardwaarde is 3.
Voorafgaand aan deze iteratie worden er geen vectoren toegevoegd aan de DIIS-subruimte. De eerste
De DIIS-stap wordt pas gezet DIIS_MIN_VECS vectoren bevinden zich in de DIIS-deelruimte.

LAG_FILE = geheel getal
Bestand (eenheidsnummer) voor het lezen van de lagrangiaanse matrix. De standaardwaarde is
momenteel 75.

OPDM_FILE = geheel getal
Bestand (eenheidsnummer) voor het lezen van de dichtheidsmatrix met één deeltje. De standaardwaarde
bedraagt ​​momenteel 73.

SCALE_GRAD = boolean
Geeft aan of de orbitale gradiënt moet worden geschaald met de geschatte (diagonale) orbitaal
Hessisch. De standaardwaarde is WAAR.

SCALE_STEP = vast
Schaalfactor voor de orbitale rotatiestap. Standaard is 1.0.

TPDM_FILE = geheel getal
Bestand (eenheidsnummer) voor het lezen van de tweedeeltjesdichtheidsmatrix. De standaardwaarde
bedraagt ​​momenteel 74.

LEVEL_SHIFT = boolean
Geeft aan of niveauverschuiving van de jutematrix moet worden toegestaan. Dit kan worden gebruikt
om ervoor te zorgen dat de jute positief definitief is voor de beginiteraties. Als
de jute is niet positief definitief in de buurt van convergentie, dan kan dit een teken zijn van
numerieke instabiliteiten in de MCSCF. De standaardwaarde is WAAR.

SHIFT = vast
Geeft aan hoeveel de diagonale elementen van de jute moeten worden verschoven voor niveau
verschuiven. De standaardwaarde is 0.01.

DETERM_MIN = vast
Minimaal toegestane waarde voor de determinant van de jutematrix bij niveauverschuiving
staat aan. Dit wordt niet gebruikt als LEVEL_SHIFT=FALSE. De standaardwaarde is 0.00001.

STEP_MAX = vast
Dit is de maximaal toegestane enkele orbitale rotatie. De standaardwaarde is 0.30.

FORCE_STEP = boolean
Hierdoor kan de gebruiker de berekende stap negeren en een stap in a forceren
bepaalde richting. Dit kan handig zijn als u de berekening wilt forceren
vanaf een zadelpunt. De standaardwaarde is FALSE.

FORCE_PAIR = geheel getal
Dit is de index voor het onafhankelijke paar dat moet worden geroteerd als
FORCE_STEP=WAAR. Anders wordt dit genegeerd. Er is geen standaard.

FORCE_VALUE = vast
Dit is de orbitale rotatiewaarde voor het roteren van de orbitalen gespecificeerd door
FORCE_PAIR. Dit wordt alleen gebruikt als FORCE_STEP is ingesteld op TRUE. De standaardwaarde is
0.0.

CHECK_HESSIAN = boolean
Bereken de eigenwaarden van de orbitale jute. Dit is goed om te controleren
zadelpunt oplossingen. Moet worden gebruikt zonder niveauverschuiving. De standaardwaarde is
FALSE.

EIGEN_VECTORS = boolean
Wilt u ook de eigenvectoren van de orbitale jute berekenen? Dit is
genegeerd als CHECK_HESSIAN=FALSE. De standaardinstelling is FALSE.

8 mei, 1998 details(1)

Gebruik detcas online met behulp van onworks.net-services