oeprop - Online w chmurze

PROGRAM:

IMIĘ

oeprop - Program Własności Jednego Elektronu

OPIS

Program oeprop oblicza wartości oczekiwane operatorów właściwości jednego elektronu za pomocą a
jednocząstkowa macierz gęstości obliczona z wektora własnego w PSIF_CHKPT lub wczytana z
plik zewnętrzny. Obecnie jest w stanie przeprowadzić analizę populacji Mulliken,
obliczanie elektrycznych momentów wielobiegunowych przez oktopole, właściwości elektrostatyczne w atomie
centra (potencjał elektrostatyczny, pole elektryczne, gradient pola elektrycznego, elektron i
gęstość spinu, udział anizotropii dipolarnej w stałych sprzężenia nadsubtelnego),
gęstość elektronów i spinu, gradient gęstości elektronów i spinów, Laplacian elektronu i
gęstości spinu, potencjał elektrostatyczny na dowolnym dwuwymiarowym (płaskim)
siatka prostokątna i wartości orbitali molekularnych w dowolnym trójwymiarowym obszarze
prostokątna siatka. Różne możliwości obejmują obliczenia relatywistyczne
jednoelektronowe poprawki pierwszego rzędu na energię (prędkość masowa i wyrażenia Darwina),
konstrukcja naturalnych orbitali molekularnych na podstawie gęstości jednej cząsteczki odczytanej z zewnątrz
plik (NO można zapisać do PSIF_CHKPT) i obliczenie zasięgu przestrzennego – oczekiwanie
wartości operatorów X^2, Y^2, Z^2, R^2 - całkowitej gęstości elektronowej i indywidualnej
MO (jeśli CZYTAJ_OPDM = fałszywe) lub naturalne (jeśli CZYTAJ_OPDM = prawda) orbitale (MPMAX musi być ustawiony na
wartość większa niż 1 do obliczania tych jednostek). Należy stosować zakresy przestrzenne
ostrożnie, gdyż zależą one od punktu odniesienia.

LITERATURA

Analiza populacji Mullikena

1. Elektroniczna analiza populacji na funkcjach fali molekularnej LCAO-MO. RS
Mulliken, J. Chem. Fiz. 23, 1833 (1955), tamże. 23, 1841 (1955), tamże. 36, 3428
(1962).

Relacje powtarzania dla całek jednoelektronowych po kartezjańskich funkcjach Gaussa.

1. Efektywne rekurencyjne obliczanie całek molekularnych po kartezjańskim Gaussa
Funkcje. S. Obara i A. Saika, J. Phys. Chem. 84, 3963 (1986).

Podstawowe stałe fizyczne i współczynniki przeliczeniowe.

1. CRC Podręcznik chemii i fizyki. Pod redakcją DR Lide’a. Wydanie 73
(1992-1993).

AKTA WYMAGANE

input.dat - Plik wejściowy
PSIF_CHKPT – plik punktu kontrolnego

AKTA AKTUALIZACJA

wyjście.dat
dipmom.dat - Momenty dipolowe
esp.dat - Potencjał elektrostatyczny na siatce 2D
edens.dat - Gęstość elektronów na siatce 2D
edgrad.dat - Gradient gęstości elektronów na siatce 2D
edlapl.dat - Laplacian gęstości elektronowej na siatce 2D
sdens.dat - Gęstość spinu na siatce 2D
sdgrad.dat - Gradient gęstości wirowania na siatce 2D
sdlapl.dat - Laplacian gęstości spinu na siatce 2D
mo.dat - Wartości orbitalu molekularnego/gęstości na siatce 3D
mo.pov - plik wejściowy MegaPov do renderowania obrazu mo.dat
mo.cube — Orbital(y) molekularny na siatce 3D w formacie Gaussa94 Cube
dens.cube - Gęstość elektronów/spinów na siatce 3D w formacie Gaussa94 Cube

WEJŚCIE FORMAT

Większość słów kluczowych nie jest konieczna w przypadku rutynowych zadań. Następujące słowa kluczowe to
ważny:

WFN = boolean
Rodzaj funkcji falowej. To słowo kluczowe jest „makro”, które pozwala użytkownikowi ustawić większość
niezbędne słowa kluczowe. Rozpoznawane są następujące wartości:

WFN = SCF - odpowiednik CZYTAJ_OPDM = fałszywy;

WFN = DETCI - odpowiednik CZYTAJ_OPDM = prawda, PLIK_OPDM = 40, OPDM_BASIS = AO,
OPDM_FORMAT = TRÓJKĄT;

WFN = CCSD - odpowiednik EAD_OPDM = prawda, PLIK_OPDM = 79, OPDM_BASIS = AO,
OPDM_FORMAT = TRÓJKĄT;

WFN = QVCCD - odpowiednik CZYTAJ_OPDM = prawda, PLIK_OPDM = 76, OPDM_BASIS = SO
OPDM_FORMAT = TRÓJKĄT;

CZYTAJ_OPDM = boolean
Ta flaga określa, czy macierz gęstości jednocząstkowej ma być odczytywana z dysku.
Wartość domyślna to fałsz.

PLIK_OPDM = liczba całkowita
Określa numer pliku macierzy gęstości pojedynczych cząstek. Wartość domyślna to 40 (plik główny). Do
zapewniają kompatybilność wsteczną z wcześniejszymi pakietami właściwości PSI (właściwy,
cipro, ccprop) zakładany jest specjalny format pliku gęstości, gdy PLIK_OPDM = 40
(właściwości obliczeniowe z gęstości CI - cipro tryb zgodności) i PLIK_OPDM =
79 (obliczanie właściwości na podstawie gęstości CC - ccprop tryb zgodności). Jak na razie
w ogólnym przypadku onepdm musi być zapisane na samym początku pliku. w
przyszły PSI będzie miał standardowy plik onepdm.

OPDM_BASIS = ciąg
Ta opcja może nie istnieć w przyszłości. Od 1 lutego 1998 roku norma dot
format pliku onepdm nie został ustawiony. To słowo kluczowe powinno mieć wartość „SO”
(odczyt w matrycy onepdm w oparciu o SO) lub „AO” (w oparciu o AO). Wartość domyślna to „SO”.

OPDM_FORMAT = ciąg
Ta opcja może nie istnieć w przyszłości. To słowo kluczowe powinno mieć wartość albo
„TRIANG” (odczytywane w macierzy onepdm w postaci dolnego trójkąta) „KWADRAT” (w kwadracie
formularz). Wartość domyślna to „TRIANG”

ASYMM_OPDM = boolean
Ta flaga określa, czy macierz gęstości jednej cząstki musi być symetryzowana.
Musi mieć wartość true, jeśli ma zostać odczytany ogólny, niesymetryczny onepdm (na przykład z a
program połączonych klastrów). To słowo kluczowe służy wyłącznie do tworzenia kodu. Istniejący PSI
Obecnie używane kody CC dają symetryczny onepdm, dlatego nie ma potrzeby ich stosowania
to słowo kluczowe. Wartość domyślna to fałsz.

ROOT = liczba całkowita
Określa, dla którego elementu głównego przeprowadzić analizę stanu wzbudzonego. Właściwy
z dysku zostanie odczytana jedna macierz gęstości cząstek. Obecnie realizowany dla
Funkcje falowe DETCI i DETCAS.

MPMAX = liczba całkowita
Ta liczba całkowita z zakresu od 1 do 3 określa najwyższy elektryczny moment wielobiegunowy
obliczone.

MPMAX = 1 - obliczany będzie tylko elektryczny moment dipolowy (domyślnie);

MPMAX = 2 - zostaną obliczone elektryczne momenty dipolowe i kwadrupolowe; MPMAX = 3 -
zostaną obliczone elektryczne momenty dipolowe, kwadrupolowe i oktopolowe.

MP_REF liczba całkowita
Parametr ten określa punkt odniesienia dla elektrycznych momentów wielobiegunowych
obliczenie.

MP_REF = 0 (domyślnie) lub 1 - środek masy;

MP_REF = 2 - początek układu współrzędnych przestrzennych;

MP_REF = 3 - środek ładunku elektronicznego;

MP_REF = 4 - środek ładunku jądrowego;

MP_REF = 5 - środek ładunku netto.

UWAGA : Według elektrodynamiki klasycznej moment elektryczny 2^(n+1)-biegunowy
jest niezależny od punktu odniesienia tylko wtedy, gdy moment elektryczny 2^(n)-biegunowy wynosi
zanik. Oznacza to, że moment dipolowy będzie zależał od punktu odniesienia, jeśli
całkowity ładunek systemu jest różny od zera. Przez analogię elektryczny moment kwadrupolowy
będzie zależeć od punktu odniesienia, jeśli system posiada niezerowy dipol elektryczny
chwila itp.

MP_REF_XYZ = real_wektor
Wektor ten określa współrzędne punktu odniesienia. Jeśli to słowo kluczowe to
obecny na wejściu MP_REF słowo kluczowe zostanie pominięte.

NUC_ESP = boolean
Ta flaga określa, czy właściwości elektrostatyczne będą obliczane w jądrach.
Aktualna lista obejmuje potencjał elektrostatyczny, pole elektryczne, pole elektryczne
gradient, gęstość elektronów i spinu oraz udział anizotropowy w nadsubtelnej warstwie
stałe sprzężenia (dwa ostatnie wymagają ustawienia SPIN_PROP do prawdy). Domyślnie jest
prawdą.

GRID = liczba całkowita
Określa typ właściwości, który ma być oceniany na siatce.

GRID = 0 (domyślnie) - nic nie oblicza;

GRID = 1 - potencjał elektrostatyczny na siatce dwuwymiarowej;

GRID = 2 - gęstość elektronów (gęstość spinu, jeśli SPIN_PROP ma wartość true) na dwu-
siatka wymiarowa;

GRID = 3 - gradient gęstości elektronów (gradient gęstości spinu jeśli SPIN_PROP jest ustawione na
true) na dwuwymiarowej siatce;

GRID = 4 - Laplacian gęstości elektronów (Laplacian gęstości spinu if
SPIN_PROP jest ustawiona na true) na siatce dwuwymiarowej. Zgodnie z konwencją
używane w terenie, tak naprawdę wykreśla się Laplacianów
znak ujemny.

GRID = 5 - wartości orbitali molekularnych na siatce trójwymiarowej.

GRID = 6 - wartości gęstości elektronowej (gradient gęstości spinu jeśli SPIN_PROP is
ustawione na true) na trójwymiarowej siatce.

GRID_FORMAT = ciąg
Określa, w jakim formacie zostanie utworzony wynik siatki. Obecnie, PLOTTV
(domyślnie dla siatek 2-wymiarowych), MEGAPOWPLUS (dostępne dla siatek 3-w) i
GAUSSKUB(domyślnie dla siatek 3-D) są obsługiwane.

MO_TO_PLOT = wektor
Określa indeksy orbitali molekularnych, które mają być obliczone na siatce 3-D. Indeksy
można określić jako:

unsigned integer - indeks w porządku Pitzera (uporządkowany według irreps, a nie
wartości własne). Zakres od 1 do liczby MO.

liczba całkowita ze znakiem - indeks względem poziomu Fermiego. +1 oznacza LUMO, +2 oznacza sekundę
najniższy orbital wirtualny, -1 oznacza HOMO itp.

Wszystkie indeksy muszą być albo niepodpisane, albo podpisane, nie można ich mieszać i dopasowywać ani ty
przyniesie nieprzewidywalne rezultaty. Domyślnie obliczane są HOMO i LUMO.

GRID_ORIGIN = real_wektor
Określa początek siatki. Prostokątna siatka, która otacza całość
cząsteczka zostanie obliczona automatycznie, jeśli GRID_ORIGIN tam jednak brakuje
nie jest wartością domyślną dla siatek 2-wymiarowych.

GRID_JEDNOSTKA_X = real_wektor
Ten wektor określa kierunek pierwszej (x) strony siatki. Tak nie jest
muszą mieć długość jednostkową. Nie ma wartości domyślnej dla siatek 2-wymiarowych.

GRID_JEDNOSTKA_Y = real_wektor
To samo dla drugiej strony (y). Nie musi mieć jednostkowej długości ani nawet
prostopadły do GRID_JEDNOSTKA_X. Nie ma ustawień domyślnych dla siatek 2-wymiarowych.

GRID_XY0 = real_2d_vector
Określa współrzędne lewego dolnego rogu prostokąta siatki w 2D
układ współrzędnych zdefiniowany przez GRID_ORIGIN, GRID_UNIT_X i GRID_UNIT_Y. Jest
brak domyślnego.

GRID_XY1 = real_2d_vector
Określa współrzędne prawego górnego rogu prostokąta siatki w 2D
układ współrzędnych zdefiniowany przez GRID_ORIGIN, GRID_UNIT_X i GRID_UNIT_Y. Jest
brak domyślnego.

GRID_XYZ0 = real_3d_vector
Określa współrzędne lewego dolnego rogu ramki siatki w widoku 3D
układ współrzędnych zdefiniowany przez GRID_ORIGIN, GRID_UNIT_X, GRID_UNIT_Y i krzyż
produkt dwóch ostatnich. Nie ma wartości domyślnej.

GRID_XYZ1 = real_3d_vector
Określa współrzędne najbliższego prawego górnego rogu ramki siatki w 3D
układ współrzędnych zdefiniowany przez GRID_ORIGIN, GRID_UNIT_X, GRID_UNIT_Y i krzyż
produkt dwóch ostatnich. Nie ma wartości domyślnej.

NIC = liczba całkowita
Liczba punktów siatki wzdłuż kierunku x. Parametr ten musi być większy niż
1. Wartość domyślna to 20.

NIE = liczba całkowita
Taki sam jak NIC w kierunku Y. Wartość domyślna to 20.

NIZ = liczba całkowita
Taki sam jak NIC dla kierunku Z. Wartość domyślna to 20.

GRID_ZMIN = Podwójna
Dolna granica wyświetlanych wartości z dla wykresów konturowych gęstości elektronowej i jej
Laplazjan. Wartość domyślna to 0.0

GRID_ZMAX = Podwójna
Górna granica wyświetlanych wartości z dla wykresów konturowych gęstości elektronowej i jej
Laplazjan. Wartość domyślna to 3.0

EDGRAD_LOGSCALE = liczba całkowita
Steruje skalowaniem logarytmicznym wytworzonego wykresu gradientu gęstości elektronów. Obroty
skalowanie wyłączone jeżeli ustawione na zero, w innym przypadku im wyższa wartość tym mocniejsza
pole gradientu zostanie przeskalowane. Zalecana wartość (domyślna) to 5.

SPIN_PROP = boolean
Flaga do obliczania właściwości spinu (analiza populacji Mulliken alfa i beta
gęstości, gęstości spinu i wkłady anizotropowe w sprzężenie nadsubtelne
stałe w centrach atomowych). Wartość domyślna to fałsz.

WYDRUKOWAĆ = liczba całkowita
To najważniejsze słowo kluczowe – od niego zależy ilość drukowanych informacji.
Obecnie używane są następujące wartości:

WYDRUKOWAĆ = 0 - tryb cichy - drukuje tylko istotne wyniki - "kompaktowe" wyniki
Analiza populacji Mullikena, elektryczne momenty wielobiegunowe i elektrostatyka
nieruchomości;

WYDRUKOWAĆ = 1 (domyślnie) - wszystkie powyższe plus lista zadań do wykonania i lista
parametrów obliczeniowych;

WYDRUKOWAĆ = 2 - wszystkie powyższe plus macierz populacji Mulliken AO i elektroniczna oraz
jądrowe składniki elektrycznego momentu dipolowego;

WYDRUKOWAĆ = 3 - wszystko powyższe plus macierz gęstości w bazie AO i momencie dipolowym
całki w bazie AO (i SO);

WYDRUKOWAĆ = 4 - wszystkie powyższe plus informacje o zestawie podstawowym, orbitale naturalne w ujęciu
orbitali symetrii, macierz nakładania się;

WYDRUKOWAĆ >= 5 - wszystkie powyższe plus wektory współczynników sprzężenia, zawód
wektor i zmodyfikowany wektor Z w bazie MO.

DRUKUJ NR = boolean
Jeśli WRTNOS = TRUE i ta opcja jest również TRUE, wydrukowane zostaną orbitale naturalne
do wyjścia, zanim zostaną zapisane w pliku punktu kontrolnego.

WRTNOS = boolean
Jeśli TRUE, naturalne orbitale zostaną zapisane w pliku punktu kontrolnego.

GRID WYDAJNOŚĆ ROLNICZE KONSPIRATORSTWO

Obecnie oeprop generuje dane wyjściowe w postaci dwuwymiarowych siatek gotowych do wykreślenia za pomocą a
program PLOTTV wersja 1.3.2. Program został napisany przez Kenny’ego Toha
(ktoh@td2cad.intel.com), programista dla Działu Technologii CAD, Intel Corp,
Święta Clara. Jest to pakiet freeware, który można pobrać z Internetu.

Siatki trójwymiarowe są wyprowadzane w formacie odpowiednim do wykreślenia za pomocą programu MegaPow
wersja 0.5. Ten darmowy program jest poprawioną wersją POV-Ray. Jest rozwijany przez A
liczby osób i można je pobrać z Internetu (przejdź do
http://nathan.kopp.com/patched.htm aby dowiedzieć się więcej). Aby renderować MO lub gęstość
image, edytuj (jeśli to konieczne) plik poleceń ruch pow stworzone przez oeprop i wykonać megapowplus
+Imo.pow Aby uzyskać więcej opcji, uruchom megapowplus -h

30 marca 2001 r. oeprop(1)

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