detcas - Online în cloud

PROGRAM:

NUME

detcas - Interacțiunea de configurare determinantă / Activ complet
Program de câmp auto-consistent

DESCRIERE

Programul detcas efectuează un câmp auto-consistent în spațiu complet activ (CASSCF)
optimizarea orbitalilor moleculari printr-o procedură în două etape în legătură cu
program de interacțiune de configurare determinant DETCI. Programul este destul de simplu și
în prezent folosește o abordare Newton-Raphson pentru a actualiza orbitalii, folosind un simplu
Hessianul orbital aproximativ. Convergența este accelerată folosind inversarea directă a lui Pulay
procedura subspațială iterativă (DIIS). Codul a fost scris pentru a permite mai multe generale
funcții de undă care nu prezintă neapărat un tratament CI complet al spațiului activ. În
în special, orice funcție de undă CI de spațiu activ restricționat (RAS) suportată de DETCI poate fi
utilizat, permițând funcții de undă MCSCF de tip RASSCF.

REFERINȚE

Hessian orbital aproximativ:

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

Spațiu activ restricționat CI:

1. Algoritmi de interacțiune pentru configurație bazată pe determinanți pentru complete și restricționate
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen și HJ Aa.
Jensen, J. Chem. Fiz. 89, 2185 (1988).

Spațiu activ restricționat SCF:

1. P.-A. Malmqvist, A. Rendell și BO Roos, J. Phys. Chim. 94, 5477 (1990).

Programul DETCI:

1. CD Sherrill, algoritmi de calcul pentru o scară largă completă și multi-referință
Configuration Interaction Wavefunctions, teză de doctorat, Universitatea din Georgia, Atena,
GA, 1996.

2. CD Sherrill și HF Schaefer, The Configuration Interaction Method: Advances
în Abordări înalt corelate, Adv. Quantum Chem. 34, 143-269 (1999).

DOSARE NECESARE

input.dat - Fișier de intrare
file78 - MO integrale cu un electron (operator fzc)
fisier72 - MO integrale cu doi electroni
file73 - MO matrice de densitate de o particulă
file74 - MO matrice de densitate cu două particule
dosar75 - MO Lagrangian

DOSARE ACTUALIZAT

output.dat - Fișier de ieșire
fisier14.dat - Inregistrarea energiilor si a gradientilor orbitali

INTRARE FORMAT

Următoarele argumente din linia de comandă sunt disponibile:

-Liniște Acest lucru dă același rezultat ca PRINT=0.

-o fname
Oferă numele fișierului pentru fișierul de ieșire. Implicit, output.dat.

Intrarea suplimentară pentru acest program este citită din fișierul input.dat. Cel mai des folosit
cuvintele cheie sunt:

WFN = şir
Tipurile de funcții de undă acceptate sunt CASSCF și RASSCF.

CONVERGENŢĂ = întreg
Convergența dorită pe gradientul orbital. Convergența se realizează atunci când RMS
a erorii în gradientul orbital este mai mică de 10**(-n). Valoarea implicită este 4 pentru
calcule de energie și 7 pentru gradienți.

DOCC = matrice_întregi
Acest vector dă numărul de orbitali dublu ocupați în fiecare irrep. Există
fără implicit.

SOCC = matrice_întregi
Acest vector oferă numărul de orbitali ocupați individual în fiecare irrep. Există
fără implicit.

ENERGY_CONVERGENCE = întreg
Convergența dorită asupra energiei totale MCSCF. Valoarea implicită este 7.

FROZEN_DOCC = matrice_întregi
Numărul de orbitali de cea mai mică energie ocupați dublu în fiecare ireductibil
reprezentare care va fi literalmente înghețată (nu este actualizată în MCSCF). The
Se folosește ordonarea din bumbac a reprezentărilor iredicibile. Valoarea implicită este
vector zero.

FROZEN_UOCC = vector_întreg
Numărul de orbitali neocupați cu cea mai mare energie din fiecare reprezentare ireductibilă
care va fi literalmente înghețat (nu este actualizat în MCSCF). Valoarea implicită este zero
vector.

RESTRICTED_DOCC = matrice_întregi
Numărul de orbitali de cea mai mică energie ocupați dublu în fiecare ireductibil
reprezentare care va fi optimizată dar păstrată dublu ocupată în MCSCF.
Acești orbitali vin după FROZEN_DOCC orbitali. Valoarea implicită este zero
vector.

RESTRICTED_UOCC = matrice_întregi
Numărul de orbitali neocupați cu cea mai mare energie din fiecare reprezentare ireductibilă
care va fi optimizat dar păstrat neocupat în MCSCF. Acești orbitali vin
inainte de FROZEN_UOCC orbitali. Valoarea implicită este vectorul zero.

NCASITER = întreg
Număr maxim de iterații pentru optimizarea orbitalilor. Această opțiune să be
specificată in il DEFAULT secțiune of intrare, pentru că trebuie să fie vizibil pentru
program de control PSI. Implicit la 1.

PRINT = întreg
Această opțiune determină gradul de nivel al ieșirii. O valoare de 1 sau 2 specifică
tipărire minimă, o valoare de 3 specifică tipărirea verbosă. Valorile de 4 sau 5 sunt
folosit pentru depanare. Nu utilizați nivelul 5 decât dacă cazul de testare este foarte mic (de ex
STO H2O CISD).

Cuvintele cheie mai puțin utilizate sunt:

DIIS_FREQ = întreg
Procedura de extrapolare DIIS va fi încercată la fiecare n iterații. Implicit
este 1.

DIIS_MAX_VECS = întreg
Numărul maxim de vectori subspațiali pentru procedura DIIS. După acest număr
de vectori sunt atinse, vectorii mai vechi vor fi eliminati din subspațiu ca
necesar. Valoarea implicită este 8.

DIIS_MIN_VECS = întreg
Numărul minim de vectori subspațiali înainte de o interpolare DIIS poate fi
efectuat. Valoarea implicită este 2. Nu are sens să aveți valori mai mici de 2.

DIIS_START = întreg
Numărul de iterație când Pulay’s Direct Inversion of the Iterative Subspace (DIIS)
procedura de accelerare a convergenței ar trebui activată. Valoarea implicită este 3.
Înainte de această iterație, vectorii nu sunt adăugați la subspațiul DIIS. Primul
Pasul DIIS nu va fi făcut până când DIIS_MIN_VECS vectorii sunt în subspațiul DIIS.

LAG_FILE = întreg
Fișier (numărul unității) pentru citirea matricei lagrangiane. Valoarea implicită este
in prezent 75.

OPDM_FILE = întreg
Fișier (numărul unității) pentru citirea matricei de densitate a unei particule. Valoarea implicită
are in prezent 73 de ani.

SCALE_GRAD = boolean
Spune dacă să scalați gradientul orbital după orbital aproximativ (diagonal).
Hessian. Valoarea implicită este TRUE.

SCALE_STEP = real
Factor de scară pentru etapa de rotație orbitală. Implicit este 1.0.

TPDM_FILE = întreg
Fișier (numărul unității) pentru citirea matricei de densitate cu două particule. Valoarea implicită
are in prezent 74 de ani.

LEVEL_SHIFT = boolean
Spune dacă se permite schimbarea nivelului matricei hessian. Acesta poate fi folosit
pentru a se asigura că hessianul este definit pozitiv pentru iterațiile de început. Dacă
Hessianul nu este pozitiv definit aproape de convergență, atunci acesta poate fi un semn al
instabilități numerice în MCSCF. Valoarea implicită este TRUE.

SHIFT = real
Spune cât de mult ar trebui să fie deplasate elementele diagonale ale hessianului pentru nivel
schimbare. Implicit este 0.01.

DETERM_MIN = real
Valoarea minimă admisă pentru determinantul matricei hessian în cazul deplasării nivelului
este pornit. Aceasta nu este utilizată dacă LEVEL_SHIFT=FALSE. Implicit este 0.00001.

STEP_MAX = real
Aceasta este rotația orbitală unică maximă permisă. Implicit este 0.30.

FORCE_STEP = boolean
Acest lucru permite utilizatorului să depășească pasul calculat și să forțeze un pas în a
direcție anume. Acest lucru poate fi util dacă încercați să forțați calculul
dintr-un punct de şa. Valoarea implicită este FALSE.

FORCE_PAIR = întreg
Acesta este indicele pentru perechea independentă care urmează să fie rotită dacă
FORCE_STEP=ADEVĂRAT. Acest lucru este ignorat altfel. Nu există implicit.

FORCE_VALUE = real
Aceasta este valoarea rotației orbitale pentru rotirea orbitalilor specificati de
FORCE_PAIR. Acesta este utilizat numai dacă FORCE_STEP este setat la TRUE. Valoarea implicită este
0.0.

CHECK_HESSIAN = boolean
Calculați valorile proprii ale hessianului orbital. Acest lucru este bun pentru a verifica
soluții de punct de șa. Ar trebui folosit fără schimbarea nivelului. Valoarea implicită este
FALS.

EIGEN_VECTORS = boolean
Doriți să calculați și vectorii proprii ai hessianului orbital. Aceasta este
ignorat dacă CHECK_HESSIAN=FALSE. Implicit este FALSE.

8 Mai, 1998 detcas(1)

Utilizați detcas online folosind serviciile onworks.net