detcas - Online in de Cloud

PROGRAMMA:

NAAM

detcas - Determinant Configuratie Interactie / Volledig Actief
Ruimte-zelf-consistent-veldprogramma

PRODUCTBESCHRIJVING

Het programma details voert een volledig actieve ruimte zelf-consistent veld (CASSCF) uit
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, waarbij gebruik wordt gemaakt van een eenvoudige
benaderende orbitale Hessiaan. Convergentie wordt versneld met behulp van Pulay's directe inversie van
de iteratieve subruimte (DIIS)-procedure. De code is geschreven om meer algemene
golffuncties die niet noodzakelijkerwijs een volledige CI-behandeling van de actieve ruimte omvatten. In
in het bijzonder kan elke beperkte actieve ruimte (RAS) CI-golffunctie die door DETCI wordt ondersteund,
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 Highly Correlated Approaches, Adv. Quantum Chem. 34, 143-269 (1999).

FILES VERPLICHT

input.dat - Invoerbestand
file78 - MO één-elektron integralen (fzc operator)
file72 - MO twee-elektronenintegralen
file73 - MO-eendeeltjesdichtheidsmatrix
file74 - MO-matrix met twee-deeltjesdichtheid
file75 - MO Lagrangiaanse

FILES UPDATE

output.dat - Uitvoerbestand
file14.dat - Registratie van energieën en orbitaalgradië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 gelezen uit het bestand input.dat. De meest gebruikte
trefwoorden zijn:

WFNMeer = snaar
De ondersteunde golffunctietypen zijn CASSCF en RASSCF.

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

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
Gewenste convergentie 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 bevroren zal worden (niet bijgewerkt in de MCSCF). De
Er wordt gebruik gemaakt van de katoenen volgorde van de onherkenbare representaties. De standaard is de
nulvector.

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

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

BEPERKT_UOCC = geheel getal_array
Het aantal onbezette orbitalen met de hoogste energie in elke onherleidbare representatie
die geoptimaliseerd zullen worden maar vrijgehouden in de MCSCF. Deze orbitalen komen
voor de BEVROREN_UOCC orbitalen. De standaard is de nulvector.

NCASITER = geheel getal
Maximaal aantal iteraties om de orbitalen te optimaliseren. Deze optie moet be
gespecificeerd in the DEFAULT sectie of invoer, omdat het zichtbaar moet zijn voor de
besturingsprogramma PSI. Standaardwaarde is 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 aantal
van vectoren worden bereikt, worden oudere vectoren uit de subruimte verwijderd
noodzakelijk. De standaardwaarde is 8.

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

DIIS_START = geheel getal
Het iteratienummer bij Pulay's directe omkering van de iteratieve deelruimte (DIIS)
De procedure voor versnelling van convergentie moet ingeschakeld zijn. De standaardwaarde is 3.
Vóór deze iteratie worden er geen vectoren toegevoegd aan de DIIS-subruimte. De eerste
DIIS-stappen zullen pas worden genomen als DIIS_MIN_VECS vectoren bevinden zich in de DIIS-deelruimte.

LAG_BESTAND = 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 van één deeltje. De standaardwaarde
bedraagt ​​momenteel 73.

SCHAAL_GRAAD = boolean
Geeft aan of de orbitaalgradiënt moet worden geschaald met de benaderende (diagonale) orbitaalgradiënt.
Hessisch. De standaardwaarde is WAAR.

SCHAALSTAP = vast
Schaalfactor voor de orbitale rotatiestap. Standaard is dit 1.0.

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

NIVEAUVERSCHUIVING = boolean
Geeft aan of niveauverschuiving van de hessische matrix moet worden toegestaan. Dit kan worden gebruikt
om ervoor te zorgen dat de hessiaan positief definitief is voor de beginiteraties. Als
de hessiaan is niet positief definitief nabij convergentie, dan kan dit een teken zijn van
Numerieke instabiliteiten in de MCSCF. De standaardwaarde is TRUE.

SHIFT = vast
Geeft aan hoeveel de diagonale elementen van het jute moeten worden verschoven voor een waterpas
verschuiven. Standaard is 0.01.

DETERM_MIN = vast
Minimaal toegestane waarde voor de determinant van de Hessische matrix bij niveauverschuiving
is ingeschakeld. Dit wordt niet gebruikt als LEVEL_SHIFT=FALSE. Standaard is dit 0.00001.

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

FORCE_STAP = boolean
Hierdoor kan de gebruiker de berekende stap overschrijven en een stap in een
bepaalde richting. Dit kan handig zijn als je 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 gedraaid als
FORCE_STEP=TRUE. Anders wordt dit genegeerd. Er is geen standaardwaarde.

FORCE_WAARDE = vast
Dit is de orbitale rotatiewaarde voor het roteren van de orbitalen zoals 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 orbitaalhessiaan. Dit is handig om te controleren op
zadelpuntoplossingen. Moeten worden gebruikt zonder niveauverschuiving. De standaard is
FALSE.

EIGEN_VECTOREN = boolean
Wil je ook de eigenvectoren van de orbitaalhessiaan berekenen? Dit is
Genegeerd indien CHECK_HESSIAN=FALSE. Standaard is FALSE.

8 mei, 1998 details(1)

Gebruik detcas online met behulp van onworks.net-services