cscf - Online in de cloud

PROGRAMMA:

NAAM

cscf - lost de Hartree-Fock-vergelijkingen op

PRODUCTBESCHRIJVING

Het programma cscf voert de iteratieve procedure uit om de Hartree-Fock-vergelijkingen op te lossen.

Dit programma is beperkt tot de D2h-symmetrie en zijn subgroepen en de orbitale beroepen
moeten gehele getallen zijn. Zo zijn bepaalde zuivere impulsmomenttoestanden afgeleid
gedeeltelijke bezetting van gedegenereerde orbitalen kan niet worden verkregen met de huidige codes. Voor
bijvoorbeeld de 2PIu-toestand (doublet PI u) van lineaire ONO afgeleid van de laagste energie
lineaire (pi u)1-configuratie kan alleen worden berekend als de 2B2u (doublet B2u) of 2B3u
(doublet B 3u) component van de 2PIu (doublet PI u) toestand, en de resulterende ruimtelijke
golffunctie zal geen PI-symmetrie hebben. In zekere zin is dit echter wenselijk,
aangezien de energie een continue functie zal zijn van de buighoek. Berekening van de energie
van gebogen configuraties als 2B2u (doublet B 2u) of 2B3u (doublet B 3u) en een pure 2PIu uitvoeren
(doublet PI u) toestand bij lineaire geometrieën resulteert in een uitgesproken discontinuïteit.

Voor het grootste deel zijn tripletstaten het gevolg van de dubbele bezetting van een dubbel gedegenereerd land
orbitaal, zoals de 3A2 (triplet A 2) toestand die voortkomt uit de (e')2 of (e")2
configuraties in D3h-symmetrie, of de 3SIGMAg-status (triplet SIGMA g) van een (pi g)2 of (pi
u)2-configuratie in Dinfh (D oneindig h)-symmetrie, zal de juiste ruimtelijke symmetrie hebben.
De singlettoestanden die het gevolg zijn van dezelfde elektronische configuraties zijn inherent
multiconfiguratie en worden als zodanig niet goed vertegenwoordigd door een enkele configuratie
golffuncties.

REFERENTIES

PK-bestand methode:

1. RC Raffenetti, Chem. Fys. Let. 20(1973)335.

Moleculaire symmetrie en HF-berekeningen met gesloten schaal:

1. M. Dupuis en HFKing, Int. J. Kwant. Chem. 11(1977)613.

DIIS voor gesloten schaal:

1. P. Pulay, Chem. Fys. Let. 73 (1980) 393.

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

Koppelingscoëfficiënten (alfa en bèta) voor open schaal:

1. CCJ Roothaan, ds. Mod. Fys. 32 (1960) 179.

Demping:

1. DR Hartree, "De berekening van atomaire structuren" (Wiley: New York) 1957.

2. MC Zerner en M. Hehenberger, Chem. Fys. Let. 62 (1979) 550.

Niveauverschuiving:

1. VR Saunders en IH Hillier, Int. J. Kwant. Chem. 7(1973)699.

CONVERGEREND CSCF

Voor moeilijke gevallen met open schaal wordt een geschikte gesloten schaal aanbevolen
eerst de berekening worden uitgevoerd (een extra elektron toevoegen of verwijderen) en vervolgens deze SCF-vector
worden gebruikt als een schatting voor de gewenste open-shell-golffunctie. Voor TCSCF-gevallen is dit altijd het geval
verstandig om eerst een gesloten shell (of misschien de juiste triplet) SCF te gebruiken en dan te gebruiken
dit is een gok voor de TCSCF.

Voor open-shell-systemen wordt een niveauverschuivingswaarde van 0.5 tot 3.0 aanbevolen. Begin met een
hoge waarde (2.0 - 3.0) voor de eerste SCF-berekening en verlaag deze vervolgens (naar 0.5 - 1.0) voor
daaropvolgende runs die een geconvergeerde SCF-vector als startpunt gebruiken.

Het is uiterst belangrijk op te merken dat deze versie van de code niet langer wordt ondersteund
OPEN TYPE. One Dan moet je . de nieuwe trefwoorden REFERENTIE en MEERDERE naar specificeren de type dan: of SCF
nodig.

INVOER FORMAT

De cscf programma zoekt via het standaard trefwoordpad (eerst SCF en DEFAULT)
voor de volgende trefwoorden:

LABEL = snaar
Dit is een tekenreeks die in de uitvoer moet worden opgenomen. Deze tekenreeks wordt niet gebruikt
door het programma. Er is geen standaard.

WFNMeer = snaar
Dit is het type golffunctie dat uiteindelijk gewenst is. De standaardwaarde is SCF.

OPEN TYPE is geen langer ondersteund

REFERENTIE = snaar
Dit specificeert het type SCF-berekening dat u wilt uitvoeren. Het kan een van zijn RHF
(voor een singlet met gesloten schaal), ROHF (voor een beperkte open shell-berekening), UHF
(voor een onbeperkte open shell-berekening), TWEECON (voor een configuratie met twee
singlet), of .. Indien . wordt gegeven, daarna alfa- en bètakoppeling
coëfficiënten moeten worden gegeven met de ALPHA en BETA trefwoorden. De standaardwaarde is RHF.

MULTP= geheel getal
Specificeert de veelheid van het molecuul. Standaard is singlet.

OPLADEN= geheel getal
Specificeert de lading van het molecuul. Standaard is 0.

DOCC = geheel getal_vector
Dit geeft het aantal dubbel bezette orbitalen in elke irreducibele
vertegenwoordiging. Er is geen standaard. Als dit niet wordt gedaan, zal CSCF dit proberen
raad de beroepen met behulp van de kern Hamiltoniaan.

socc = geheel getal_vector
Dit geeft het aantal afzonderlijk bezette orbitalen in elke irreducibele
vertegenwoordiging. Er is geen standaard.

DERTYPE = snaar
Dit specificeert de volgorde van de afgeleide die uiteindelijk moet worden uitgevoerd. Het is gebruikt
Door de scf programma om te bepalen of bepaalde bestanden moeten worden geschreven en dat gebeurt ook
gebruikt om de standaardconvergentie van de golffunctie te bepalen. De standaardwaarde is
EERST.

MAXITER = geheel getal
Dit geeft het maximale aantal iteraties. De standaardwaarde is 40.

CONVERGENTIE = geheel getal
Dit specificeert hoe strak de golffunctie zal worden geconvergeerd. Convergentie wel
bepaald door de RMS-verandering in de dichtheidsmatrix ("delta P") te vergelijken met de
gegeven waarde. Het convergentiecriterium is 10**(-geheel getal). De standaardwaarde is 7 als
zowel DERTYPE = GEEN en WFNMeer = SCF worden gegeven en anders 10.

NIVEAUVERSCHUIVING = vast
Dit specificeert de niveauverschuiving. De standaardwaarde is 1.

DIRECT = boolean
Specificeert of de SCF-berekening moet worden uitgevoerd met een integrale directe techniek. De
standaard is onwaar.

AFDRUKKEN_MOS = boolean
Specificeert of de moleculaire orbitalen wel of niet moeten worden afgedrukt. De standaardwaarde is onwaar.

Er zijn ook een groot aantal minder vaak gebruikte invoerparameters. Als je niet
Als u begrijpt wat de volgende opties betekenen, zorg er dan voor dat ze niet in uw
invoer. De standaardinstellingen zullen in de overgrote meerderheid van de gevallen werken. Deze zijn gespecificeerd
met de volgende trefwoorden:

VERWIJDEREN_INTS = boolean
Integraalbestanden worden gewist als WFNMeer = SCF en DERTYPE = EERST or DERTYPE = GEEN.
Als u integralenbestanden wilt behouden, stelt u dit in VERWIJDEREN_INTS = onwaar. De standaardwaarde is
waar.

NABESTELLEN = snaar
De parameter regelt de herschikking van moleculaire orbitalen. Indien ingesteld op VOOR dan de
denk dat de orbitalen uit het controlepuntbestand opnieuw worden gerangschikt. Indien ingesteld op NA, geconvergeerd
orbitalen worden opnieuw geordend voordat ze naar het controlepuntbestand worden geschreven. In beide
geval MOORDER parameter moet worden opgegeven om de herschikkingskaart op te geven. De standaardwaarde is
orbitalen niet opnieuw ordenen.

MOORDER = geheel getal_vector
Dit specificeert een moleculaire orbitale herordeningsvector. Het zal alleen worden gebruikt als
NABESTELLEN is ingesteld. Deze vector brengt elke orbitaal in kaart naar zijn nieuwe index, bijv MOORDER = (0
2 1) specificeert dat na het herschikken van orbitalen 1 en 2 zullen worden verwisseld. De rang van
deze vector is hetzelfde als het aantal MO's. De indexen staan ​​in Pitzer-volgorde
(geordend op symmetrie en vervolgens op energie binnen elk symmetrieblok), basis-0. CSCF
zal waarschijnlijk mislukken als de gegeven MOORDER orbitalen van verschillende irreps combineert. Daar
is geen standaard.

ALPHA = echte_vector
If OPEN TYPE = ., dan geeft deze parameter de alfakoppelingscoëfficiënten.
Het aantal elementen in deze vector is MM(MM+1)/2, waarbij MM het aantal is
onherleidbare representaties die afzonderlijk bezette moleculaire orbitalen bevatten. Daar
is geen standaard.

BETA = echte_vector
If OPEN TYPE = ., dan geeft deze parameter de bètakoppelingscoëfficiënten.
Het aantal elementen in deze vector is MM(MM+1)/2, waarbij MM het aantal is
onherleidbare representaties die afzonderlijk bezette moleculaire orbitalen bevatten. Daar
is geen standaard.

GUESS = snaar
Deze optie bepaalt het type initiële schatting van de eigenvector die CSCF zal gebruiken.
De enige geldige optie op dit moment zijn: (1) GUESS = CORE, waardoor het wordt gebruikt
kern Hamiltoniaanse eigenvector om de berekening te starten; (2) GUESS = AUTO welke
resulteert in een poging om de MO-vector in het controlepuntbestand te gebruiken, of neemt zijn toevlucht tot
kerngok als er geen eigenvector in dat bestand staat. De standaard als AUTO.

IDRUK = geheel getal
Dit is een afdrukoptie. De standaardwaarde is 0.

MO_OUT = boolean
Drukt de orbitalen af ​​met symmetrie en bezettingen aan het einde van de
berekening. Standaard is waar.

DRAAIEN = boolean
Als dit niet waar is, zullen de moleculaire orbitalen niet worden geroteerd. Alleen de rotatie
beïnvloedt de virtuele orbitalen voor open schaalsystemen. Deze parameter moet waar zijn
voor gecorreleerde gradiënten en moet onwaar zijn voor tweede en hogere afgeleiden.
De standaardwaarde is false if WFNMeer = SCF en anders waar.

CHECK_ROT = boolean
Controleer de hierboven beschreven moleculaire orbitale rotatie om er zeker van te zijn dat er geen kolommen zijn
de SCF-eigenvectormatrix worden door de rotatie verwisseld. Heeft geen effect als DRAAIEN =
vals. De standaardwaarde is waar.

CHECK_MO_ORTHOGONALITEIT = boolean
Controleer of de moleculaire orbitalen orthonormaal zijn. Alleen nuttig voor het debuggen. De
standaard is onwaar.

DIS = boolean
Dit bepaalt of diis wordt gebruikt. De standaardwaarde is waar.

DISSTART = geheel getal
Dit geeft de eerste iteratie waarvoor DIIS zal worden gebruikt. De standaardwaarde is 0.

NDIIS = geheel getal
Dit geeft het aantal foutmatrices dat moet worden gebruikt in de diis-procedure. De standaard
is 6 voor gesloten schaal, 4 voor open schaal en 3 voor tcscf.

DIISDAMP = vast
Dit geeft de dempingsfactor voor de diis-procedure. De standaardwaarde is 0.0 voor
gesloten schaal, 0.02 voor open schaal en 0.01 voor tcscf.

INCR = vast
Dit wordt in tcscf gebruikt om te bepalen hoe vaak de ci-coëfficiënten opnieuw worden berekend.
Een klein aantal (~0.25) zorgt ervoor dat ze bijna elke scf opnieuw worden berekend
iteratie. De standaardwaarde is 0.25.

DYN_ACC = boolean
Bij het uitvoeren van directe scf specificeert dit of er sprake is van dynamische integrale nauwkeurigheidsafsnijdingen
zal gebruikt worden. Standaard is waar (gebruik dynamische grenswaarden). De eerste iteraties zijn
uitgevoerd met integralen die tot op zes cijfers nauwkeurig zijn. Nadat de dichtheid is geconvergeerd naar
10^-5 of 30 iteraties zijn voltooid, volledige integrale nauwkeurigheid wordt gebruikt. Als scf
Er worden convergentieproblemen ondervonden die dynamische grenswaarden uitschakelen door dit in te stellen
variabel naar false kan helpen.

ORTHOG_ONLY = boolean
Soms is dit in CASSCF of andere niet-HF/KS-schema's voor orbitale optimalisatie wel het geval
handig om MO's uit andere geometrieën opnieuw te orthogonaliseren voor de huidige geometrie
ze kunnen worden gebruikt als een eerste schatting voor de nieuwe MO's. Dit kan worden uitgevoerd door
CSCF uitvoeren met ORTHOG_ONLY = waar. Nadat de orbitalen orthogonaal zijn gemaakt, wordt de
het programma wordt afgesloten zonder een SCF-berekening uit te voeren. Dit trefwoord zal zijn
genegeerd als er geen eerdere orbitalen in het controlepuntbestand aanwezig zijn. Standaard ingesteld op waar
if WFNMeer = DETCAS.

30 mei, 1991 cscf(1)

Gebruik cscf online met behulp van onworks.net-services