cscf — Online w chmurze

PROGRAM:

IMIĘ

cscf - rozwiązuje równania Hartree-Focka

OPIS

Program cscf przeprowadza procedurę iteracyjną w celu rozwiązania równań Hartree-Focka.

Program ten ogranicza się do symetrii D2h i jej podgrup oraz obiektów orbitalnych
muszą być liczbami całkowitymi. Zatem pewne stany czystego momentu pędu pochodzą z
Za pomocą obecnych kodów nie można uzyskać częściowego zajęcia zdegenerowanych orbitali. Dla
na przykład stan 2PIu (dublet PI u) liniowego ONO uzyskany z najniższej energii
konfigurację liniową (pi u)1 można obliczyć tylko jako 2B2u (dublet B2u) lub 2B3u
(dublet B 3u) składowa stanu 2PIu (dublet PI u) i wynikająca z tego przestrzenność
funkcja falowa nie będzie miała symetrii PI. W pewnym sensie jest to jednak pożądane,
ponieważ energia będzie ciągłą funkcją kąta zgięcia. Obliczanie energii
konfiguracji wygiętych jako 2B2u (dublet B 2u) lub 2B3u (dublet B 3u) i wykonanie czystego 2PIu
(dublet PI u) przy geometrii liniowej powoduje wyraźną nieciągłość.

W przeważającej części stany trójkowe wynikające z podwójnego zajęcia podwójnie zdegenerowanego
orbitalny, taki jak stan 3A2 (triplet A 2) wynikający z (e')2 lub (e")2
konfiguracje w symetrii D3h lub stan 3SIGMAg (triplet SIGMA g) stanu (pi g)2 lub (pi
konfiguracja u)2 w symetrii Dinfh (D nieskończoność h), będzie miała odpowiednią symetrię przestrzenną.
Stany singletowe wynikające z tych samych konfiguracji elektronicznych są z natury
wielokonfiguracyjne i jako takie nie są dobrze reprezentowane przez pojedynczą konfigurację
funkcje falowe.

LITERATURA

Metoda pliku PK:

1. RC Raffenetti, Chem. Fiz. Łotysz. 20 (1973) 335.

Obliczenia symetrii molekularnej i HF z zamkniętą powłoką:

1. M.Dupuis i HFKing, Int. J. Quant. Chem. 11 (1977) 613.

DIIS dla zamkniętej powłoki:

1. P. Pulay, Chem. Fiz. Łotysz. 73 (1980) 393.

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

Współczynniki sprzężenia (alfa i beta) dla otwartej powłoki:

1. CCJ Roothaan, ks. Mod. Fiz. 32 (1960) 179.

Tłumienie:

1. DR Hartree, „Obliczanie struktur atomowych” (Wiley: Nowy Jork) 1957.

2. MC Zerner i M. Hehenberger, Chem. Fiz. Łotysz. 62 (1979) 550.

Zmiana poziomu:

1. VR Saunders i IH Hillier, Int. J. Quant. Chem. 7 (1973) 699.

ZBIEŻNY CSCF

W przypadku trudnych przypadków otwartych zaleca się zastosowanie odpowiedniego zamkniętego korpusu
najpierw należy przeprowadzić obliczenia (dodać lub usunąć dodatkowy elektron), a następnie ustawić wektor SCF
można wykorzystać jako przypuszczenie pożądanej funkcji falowej otwartej powłoki. W przypadku TCSCF tak jest zawsze
mądrze jest najpierw uruchomić zamkniętą powłokę (lub być może odpowiednią trójkę) SCF, a następnie użyć
to jako przypuszczenie dla TCSCF.

W przypadku systemów z otwartą powłoką zalecana jest wartość przesunięcia poziomu od 0.5 do 3.0. Zacznij od
wysoką wartość (2.0 - 3.0) dla pierwszego obliczenia SCF, a następnie zmniejsz ją (do 0.5 - 1.0) dla
kolejne przebiegi, w których punktem początkowym jest zbieżny wektor SCF.

Niezwykle ważne jest, aby pamiętać, że ta wersja kodu nie jest już obsługiwana
TYP OTWARTY. jeden musi posługiwać się dotychczasowy nowa słowa kluczowe REFERENCE i WIELE do sprecyzować dotychczasowy rodzaj of SCF
potrzebne.

WEJŚCIE FORMAT

Połączenia cscf program przeszukuje domyślną ścieżkę słów kluczowych (najpierw SCF , a następnie DEFAULT)
dla następujących słów kluczowych:

LABEL = ciąg
Jest to ciąg znaków, który należy uwzględnić w wynikach. Ten ciąg nie jest używany
przez program. Nie ma wartości domyślnej.

WFN = ciąg
Jest to ostatecznie pożądany typ funkcji falowej. Wartość domyślna to SCF.

TYP OTWARTY is Nie dłużej utrzymany

REFERENCE = ciąg
Określa typ obliczeń SCF, które chcemy wykonać. Może to być jeden z RFN
(dla singletu o zamkniętej skorupie), ROHF (dla ograniczonych obliczeń w otwartej powłoce), UHF
(dla nieograniczonych obliczeń w otwartej powłoce), DWA KON (dla konfiguracji dwóch
singlet) lub SPECIAL. Jeśli SPECIAL jest dane, a następnie sprzężenie alfa i beta
współczynniki należy podać z ALPHA i BETA słowa kluczowe. Wartość domyślna to RFN.

MULTP= liczba całkowita
Określa krotność cząsteczki. Wartość domyślna to singlet.

OPŁATA= liczba całkowita
Określa ładunek cząsteczki. Domyślnie wynosi 0.

DOCC = wektor_całkowity
Daje to liczbę podwójnie zajętych orbitali w każdym nieredukowalnym
reprezentacja. Nie ma wartości domyślnej. Jeśli nie zostanie to podane, CSCF spróbuje to zrobić
zgadnij zawody, używając rdzenia Hamiltona.

SOCC = wektor_całkowity
Daje to liczbę pojedynczo zajętych orbitali w każdym nieredukowalnym
reprezentacja. Nie ma wartości domyślnej.

DERTYP = ciąg
Określa to kolejność pochodnych, która ma zostać ostatecznie wykonana. To jest używane
przez SCF program, aby określić, czy określone pliki mają zostać zapisane i tak też jest
używany do określenia domyślnej zbieżności funkcji falowej. Wartość domyślna to
PIERWSZY.

MAXITER = liczba całkowita
Daje to maksymalną liczbę iteracji. Wartość domyślna to 40.

KONWERGENCJA = liczba całkowita
Określa, jak mocno funkcja falowa będzie zbieżna. Konwergencja jest
określa się poprzez porównanie zmiany RMS w macierzy gęstości („delta P”) z
podana wartość. Kryterium zbieżności wynosi 10**(-liczba całkowita). Wartość domyślna to 7, jeśli
obie DERTYP = BRAK i WFN = SCF są podane, a 10 w przeciwnym razie.

PRZESUNIĘCIE POZIOMU = real
Określa przesunięcie poziomu. Wartość domyślna to 1.

BEZPOŚREDNI = boolean
Określa, czy wykonać obliczenia SCF przy użyciu techniki całkowania bezpośredniego. The
wartość domyślna to fałsz.

DRUKUJ_MOS = boolean
Określa, czy drukować orbitale molekularne, czy nie. Wartość domyślna to fałsz.

Istnieje również duża liczba rzadziej używanych parametrów wejściowych. Jeśli nie
zrozumiesz, co oznaczają poniższe opcje, a następnie upewnij się, że nie pojawiają się one w pliku
wejście. Domyślne ustawienia będą działać w zdecydowanej większości przypadków. Są one określone
z następującymi słowami kluczowymi:

DELETE_INTS = boolean
Pliki całkowe zostaną usunięte, jeśli WFN = SCF i DERTYP = PIERWSZY or DERTYP = BRAK.
Jeśli chcesz zachować pliki całek, ustaw DELETE_INTS = fałsz. Wartość domyślna to
prawdą.

ZAMÓW PONOWNIE = ciąg
Parametr kontroluje zmianę kolejności orbitali molekularnych. Jeśli ustawione na PRZED a później
zgadnij, orbitale z pliku punktu kontrolnego zostaną ponownie uporządkowane. Jeśli ustawione na PO, zbiegły się
kolejność orbitali zostanie zmieniona przed zapisaniem do pliku punktu kontrolnego. W obu
walizka MOORDER Aby określić mapę zmiany kolejności, należy podać parametr. Wartość domyślna to
nie zmieniać kolejności orbitali.

MOORDER = wektor_całkowity
Określa wektor zmiany kolejności orbitalu molekularnego. Będzie używany tylko wtedy, gdy
ZAMÓW PONOWNIE jest ustawiony. Ten wektor odwzorowuje każdy orbital na jego nowy indeks, np MOORDER = (0
2 1) określa, że ​​po zmianie kolejności orbitale 1 i 2 zostaną zamienione. Ranga
wektor ten jest taki sam jak liczba MO. Indeksy są w porządku Pitzera
(uporządkowane według symetrii, następnie według energii w każdym bloku symetrii), podstawa-0. CSCF
prawdopodobnie zakończy się niepowodzeniem, jeśli dany MOORDER miesza orbitale z różnych przedstawicieli. Tam
nie jest domyślny.

ALPHA = real_wektor
If TYP OTWARTY = SPECIAL, wówczas ten parametr podaje współczynniki sprzężenia alfa.
Liczba elementów w tym wektorze to MM(MM+1)/2, gdzie MM jest liczbą
nieredukowalne reprezentacje zawierające pojedynczo zajęte orbitale molekularne. Tam
nie jest domyślny.

BETA = real_wektor
If TYP OTWARTY = SPECIAL, wówczas ten parametr daje współczynniki sprzężenia beta.
Liczba elementów w tym wektorze to MM(MM+1)/2, gdzie MM jest liczbą
nieredukowalne reprezentacje zawierające pojedynczo zajęte orbitale molekularne. Tam
nie jest domyślny.

ODGADNĄĆ = ciąg
Ta opcja określa typ początkowego odgadnięcia wektora własnego, którego użyje CSCF.
Jedyną ważną opcją w tej chwili są: (1) ODGADNĄĆ = CORE, co powoduje jego użycie
rdzeń wektora własnego Hamiltona, aby rozpocząć obliczenia; (2) ODGADNĄĆ = AUTO który
powoduje próbę użycia wektora MO w pliku punktu kontrolnego lub ucieka się do
podstawowe przypuszczenie, jeśli w tym pliku nie ma wektora własnego. Domyślne jeśli AUTO.

IDRUK = liczba całkowita
To jest opcja drukowania. Wartość domyślna to 0.

MO_OUT = boolean
Wypisuje orbitale z symetrią i zajęciami na końcu
obliczenie. Wartość domyślna to prawda.

OBRACAĆ SIĘ = boolean
Orbitale molekularne nie zostaną obrócone, jeśli jest to fałszywe. Tylko rotacja
wpływa na orbitale wirtualne dla systemów z otwartą powłoką. Ten parametr musi mieć wartość true
dla skorelowanych gradientów i musi być fałszywe dla drugiej i wyższych pochodnych.
Wartość domyślna to false, jeśli WFN = SCF i prawdą jest, że jest inaczej.

CHECK_ROT = boolean
Sprawdź opisaną powyżej rotację orbit molekularnych, aby upewnić się, że nie ma żadnych kolumn
macierz wektorów własnych SCF są zamieniane przez obrót. Nie ma żadnego efektu, jeśli OBRACAĆ SIĘ =
fałszywy. Wartość domyślna to prawda.

CHECK_MO_ORTHOGONALNOŚĆ = boolean
Sprawdź, czy orbitale molekularne są ortonormalne. Przydatne tylko do debugowania. The
wartość domyślna to fałsz.

DIIS = boolean
To określa, czy diis będzie używane. Wartość domyślna to prawda.

ROZRUCH = liczba całkowita
Daje to pierwszą iterację, w której zostanie użyty DIIS. Wartość domyślna to 0.

NDIIS = liczba całkowita
Daje to liczbę macierzy błędów do wykorzystania w procedurze diis. Domyślny
wynosi 6 dla powłoki zamkniętej, 4 dla powłoki otwartej i 3 dla tcscf.

DISDAMP = real
Daje to współczynnik tłumienia dla procedury diis. Wartość domyślna to 0.0 dla
powłoka zamknięta, 0.02 dla powłoki otwartej i 0.01 dla tcscf.

WZR = real
Jest to wykorzystywane w tcscf w celu określenia częstotliwości ponownego obliczania współczynników ci.
Mała liczba (~0.25) spowoduje, że będą one przeliczane prawie przy każdym scf
iteracja. Wartość domyślna to 0.25.

DYN_ACC = boolean
Podczas wykonywania bezpośredniego scf określa to, czy odcięcia dynamicznej dokładności całkowania
będzie użyty. Wartość domyślna to true (użyj dynamicznych odcięć). Początkowe iteracje są
wykonywane z całkami z dokładnością do sześciu cyfr. Po osiągnięciu gęstości
Zakończono 10^-5 lub 30 iteracji, używana jest pełna dokładność całkowa. Jeśli scf
Występują problemy ze zbieżnością, wyłączając dynamiczne odcięcia poprzez ustawienie tej opcji
zmienna na false może pomóc.

ORTHOG_ONLY = boolean
Czasami tak jest w CASSCF lub innych schematach optymalizacji orbity innych niż HF/KS
przydatne do reortogonalizacji MO z innych geometrii dla bieżącej geometrii, tzw
można je wykorzystać jako wstępne przypuszczenie dla nowych MO. Można to wykonać poprzez
uruchamianie CSCF z ORTHOG_ONLY = prawdziwy. Po ortogonalizacji orbitali następuje:
program zakończy działanie bez wykonywania obliczeń SCF. To słowo kluczowe będzie
ignorowane, jeśli w pliku punktu kontrolnego nie ma poprzednich orbitali. Domyślnie prawdziwy
if WFN = DETCAS.

30 maj, 1991 cscf(1)

Korzystaj z cscf online, korzystając z usług onworks.net