oeprop - Online în cloud

PROGRAM:

NUME

oeprop - Programul de proprietate cu un electron

DESCRIERE

Programul oeprop calculează valorile așteptate ale operatorilor de proprietate cu un electron folosind a
matrice de densitate a unei particule calculată dintr-un vector propriu în PSIF_CHKPT sau citită dintr-un
fișier extern. În prezent, este capabil să efectueze analiza populației Mulliken,
calcularea momentelor electrice multipolare prin octopol, proprietăți electrostatice la atom
centrii (potențial electrostatic, câmp electric, gradient de câmp electric, electroni și
densitatea spinului, contribuția anizotropă dipolară la constantele de cuplare hiperfine),
densitatea de electroni și spin, gradient de densitate de electroni și spin, Laplacian de electron și
densități de spin, potențial electrostatic peste un arbitrar bidimensional (planar)
grilă dreptunghiulară și valorile orbitalilor moleculari peste un tridimensional arbitrar
grilă dreptunghiulară. Capacitățile diverse includ calculul relativistului
corecții de ordinul întâi cu un electron la energie (masă-viteză și termeni Darwin),
construirea orbitalilor moleculari naturali din densitatea unei particule citite dintr-un exterior
fișier (NO pot fi scrise în PSIF_CHKPT) și calculul întinderilor spațiale - așteptare
valori ale operatorilor X^2, Y^2, Z^2 și R^2 - ale densității electronilor totale și ale individuale
MO (dacă READ_OPDM = fals) sau natural (dacă READ_OPDM = adevărat) orbitali (MPMAX trebuie setat la
o valoare mai mare decât 1 pentru calcularea acestor entităţi). Ar trebui folosite întinderi spațiale
cu prudență, deoarece depind de punctul de referință.

REFERINȚE

Analiza populației Mulliken

1. Analiza electronică a populației pe funcțiile de undă moleculară LCAO-MO. RS
Mulliken, J. Chem. Fiz. 23, 1833 (1955), ibid. 23, 1841 (1955), ibid. 36, 3428
(1962).

Relații de recurență pentru integralele cu un electron asupra funcțiilor carteziene gaussiene.

1. Calcul recursiv eficient al integralelor moleculare peste gaussian cartezian
funcții. S. Obara și A. Saika, J. Phys. Chim. 84, 3963 (1986).

Constante fizice fundamentale și factori de conversie.

1. CRC Manual de chimie și fizică. Editat de DR Lide. ediția a 73-a
(1992-1993).

DOSARE NECESARE

input.dat - Fișier de intrare
PSIF_CHKPT - Fișier punct de control

DOSARE ACTUALIZAT

output.dat
dipmom.dat - Momente dipol
esp.dat - Potențial electrostatic pe o rețea 2D
edens.dat - Densitatea electronilor pe o grilă 2D
edgrad.dat - Gradient de densitate de electroni pe o grilă 2D
edlapl.dat - Laplacian al densității electronilor pe o grilă 2D
sdens.dat - Densitatea de rotire pe o grilă 2D
sdgrad.dat - Gradient de densitate de rotire pe o grilă 2D
sdlapl.dat - Laplacian al densității de spin pe o grilă 2D
mo.dat - Valorile orbitale/densității moleculare pe o grilă 3D
mo.pov - Fișier de intrare MegaPov pentru redarea unei imagini a mo.dat
mo.cube - Orbiti moleculari pe o grilă 3D în format Gaussian94 Cube
dens.cube - Densitate(i) de electroni/spin pe o grilă 3D în format Gaussian94 Cube

INTRARE FORMAT

Majoritatea cuvintelor cheie nu sunt necesare pentru sarcinile de rutină. Următoarele cuvinte cheie sunt
valabil:

WFN = boolean
Tipul funcției de undă. Acest cuvânt cheie este o „macro” care permite utilizatorului să seteze majoritatea
cuvintele cheie necesare. Sunt recunoscute următoarele valori:

WFN = SCF - echivalent cu READ_OPDM = fals;

WFN = DETCI - echivalent cu READ_OPDM = Adevărat, OPDM_FILE = 40, OPDM_BASIS = AO,
OPDM_FORMAT = TRIANG;

WFN = CCSD - echivalent cu EAD_OPDM = Adevărat, OPDM_FILE = 79, OPDM_BASIS = AO,
OPDM_FORMAT = TRIANG;

WFN = QVCCD - echivalent cu READ_OPDM = Adevărat, OPDM_FILE = 76, OPDM_BASIS = ASA DE,
OPDM_FORMAT = TRIANG;

READ_OPDM = boolean
Acest flag specifică dacă matricea de densitate a unei particule trebuie citită de pe disc.
Implicit este fals.

OPDM_FILE = întreg
Specifică numărul fișierului matricei cu densitatea unei particule. Implicit este 40 (fișier principal). La
să ofere compatibilitate inversă cu pachetele anterioare de proprietăți PSI (adecvat,
ciprop, ccprop) format special al fișierului de densitate se presupune când OPDM_FILE = 40
(proprietăți de calcul din densitatea CI - ciprop modul de compatibilitate) și OPDM_FILE =
79 (proprietăți de calcul din densitatea CC - ccprop Mod de compatibilitate). De acum,
în cazul generic onepdm trebuie să fie scris chiar la începutul fișierului. În
viitorul PSI va avea un fișier standard onepdm.

OPDM_BASIS = şir
Este posibil ca această opțiune să nu mai existe în viitor. Începând cu 1 februarie 1998, un standard pentru
formatul de fișier onepdm nu a fost setat. Acest cuvânt cheie ar trebui să fie setat la „SO”
(citește în matrice onepdm pe bază SO) sau „AO” (în bază AO). Implicit este „SO”.

OPDM_FORMAT = şir
Este posibil ca această opțiune să nu mai existe în viitor. Acest cuvânt cheie ar trebui să fie setat la oricare
„TRIANG” (citește în matrice onepdm în formă triunghiulară inferioară) din „PATRAT” (în pătrat
formă). Implicit este „TRIANG”

ASYMM_OPDM = boolean
Acest indicator specifică dacă matricea de densitate a unei particule trebuie să fie simetrizată.
Trebuie să fie setat la adevărat dacă se citește onepdm generic nesimetric (de exemplu, dintr-un
program cu clustere cuplate). Acest cuvânt cheie este doar pentru dezvoltarea codului. PSI existent
Codurile CC utilizate acum produc onepdm simetric, prin urmare nu este nevoie de utilizare
acest cuvânt cheie. Implicit este fals.

ROOT = întreg
Aceasta specifică pentru ce rădăcină să se facă analiza stării excitate. Apropiatul
o matrice de densitate a particulelor va fi citită de pe disc. În prezent implementat pentru
Funcțiile de undă DETCI și DETCAS.

MPMAX = întreg
Acest număr întreg între 1 și 3 specifică cel mai mare moment electric multipolar care urmează să fie
calculat.

MPMAX = 1 - se va calcula doar momentul dipolului electric (implicit);

MPMAX = 2 - se vor calcula momentele dipolului și cvadrupolului electric; MPMAX = 3 -
vor fi calculate momentele electrice de dipol, cvadrupol și octopol.

MP_REF întreg
Acest parametru specifică punctul de referință pentru momentele multipolare electrice
calcul.

MP_REF = 0 (implicit) sau 1 - centrul de masă;

MP_REF = 2 - originea sistemului de coordonate spațiale;

MP_REF = 3 - centrul de încărcare electronică;

MP_REF = 4 - centrul sarcinii nucleare;

MP_REF = 5 - centrul de încărcare netă.

ATENȚIE : Conform electrodinamicii clasice, momentul electric de 2^(n+1)-poli
este independent de punctul de referință numai dacă momentul electric 2^(n)-pol este
dispărând. Înseamnă că momentul dipol va depinde de punctul de referință dacă
sarcina totală a sistemului este diferită de zero. Prin analogie, moment cvadrupol electric
va depinde de punctul de referință dacă sistemul posedă un dipol electric diferit de zero
moment etc.

MP_REF_XYZ = vector_real
Acest vector specifică coordonatele punctului de referință. Dacă acest cuvânt cheie este
prezente în intrare MP_REF cuvântul cheie va fi ignorat.

NUC_ESP = boolean
Acest indicator specifică dacă proprietățile electrostatice vor fi calculate la nuclee.
Lista curentă include potențialul electrostatic, câmpul electric, câmpul electric
gradient, densitate de electroni și spin și contribuție anizotropă la hiperfin
constante de cuplare (ultimele două necesită setare SPIN_PROP la adevărat). Implicit este
Adevărat.

GRID = întreg
Specifică tipul de proprietate care urmează să fie evaluată pe o grilă.

GRID = 0 (implicit) - nu calculează nimic;

GRID = 1 - potenţial electrostatic pe o reţea bidimensională;

GRID = 2 - densitatea electronică (densitatea spin dacă SPIN_PROP este setat la adevărat) pe două
grilă dimensională;

GRID = 3 - gradient de densitate electronică (gradient de densitate de spin dacă SPIN_PROP este setat la
adevărat) pe o grilă bidimensională;

GRID = 4 - Laplacian al densității electronilor (Laplacian al densității de spin dacă
SPIN_PROP este setat la adevărat) pe o grilă bidimensională. Conform conventiei
folosit în câmp, ceea ce se complotează de fapt sunt laplacienii luați cu
semn negativ.

GRID = 5 - valorile orbitalilor moleculari pe o grilă tridimensională.

GRID = 6 - valorile densității electronilor (gradient de densitate de spin dacă SPIN_PROP is
setat la adevărat) pe o grilă tridimensională.

GRID_FORMAT = şir
Specifică în ce format va fi produsă ieșirea grilei. În prezent, PLOTMTV
(implicit pentru grile 2-d), MEGAPOVPLUS (disponibil pentru grile 3-d) și
GAUSSCUBE(implicit pentru grile 3-d) sunt acceptate.

MO_TO_PLOT = vector
Specifică indici ai orbitalilor moleculari care urmează să fie calculați pe grila 3-d. Indici
poate fi specificat ca:

întreg fără semn - index în ordinea Pitzer (ordonat conform irreps, nu
valori proprii). Variază de la 1 la numărul de MO.

întreg cu semn - indice față de nivelul Fermi. +1 înseamnă LUMO, +2 înseamnă secundă
cel mai jos orbital virtual, -1 înseamnă HOMO etc.

Toți indicii trebuie să fie fie nesemnați, fie semnați, nu puteți amesteca și potriviți, sau dvs
va obține rezultate imprevizibile. Implicit este să se calculeze HOMO și LUMO.

GRID_ORIGIN = vector_real
Specifică originea grilei. O cutie grilă dreptunghiulară care învăluie întregul
molecula va fi calculată automat dacă GRID_ORIGIN lipsește, totuși, acolo
nu este implicit pentru grilele 2-d.

GRID_UNIT_X = vector_real
Acest vector specifică direcția primei (x) laturi a rețelei. Nu este
trebuie să aibă unitatea de lungime. Nu există o prestație implicită pentru grilele 2-d.

GRID_UNIT_Y = vector_real
Același lucru pentru a doua parte (y). Nu trebuie să fie de unitate de lungime sau chiar
ortogonal la GRID_UNIT_X. Nu există o prestabilită pentru grilele 2-d.

GRID_XY0 = real_2d_vector
Specifică coordonatele colțului din stânga jos al dreptunghiului grilei în 2D
sistem de coordonate definit de GRID_ORIGIN, GRID_UNIT_X și GRID_UNIT_Y. Există
fără implicit.

GRID_XY1 = real_2d_vector
Specifică coordonatele colțului din dreapta sus al dreptunghiului grilei în 2D
sistem de coordonate definit de GRID_ORIGIN, GRID_UNIT_X și GRID_UNIT_Y. Există
fără implicit.

GRID_XYZ0 = real_3d_vector
Specifică coordonatele colțului din stânga jos al casetei grilei în 3D
sistem de coordonate definit de GRID_ORIGIN, GRID_UNIT_X, GRID_UNIT_Y și încrucișarea
produs al ultimelor două. Nu există implicit.

GRID_XYZ1 = real_3d_vector
Specifică coordonatele colțului din dreapta sus al casetei grilei în 3D
sistem de coordonate definit de GRID_ORIGIN, GRID_UNIT_X, GRID_UNIT_Y și încrucișarea
produs al ultimelor două. Nu există implicit.

NIX = întreg
Numărul de puncte ale grilei de-a lungul direcției x. Acest parametru trebuie să fie mai mare decât
1. Implicit este 20.

NIY = întreg
La fel ca NIX pentru direcția y. Implicit este 20.

NIZ = întreg
La fel ca NIX pentru direcția z. Implicit este 20.

GRID_ZMIN = dubla
Limita inferioară a valorilor z afișate pentru graficele de contur ale densității electronilor și ale acesteia
Laplacian. Implicit este 0.0

GRID_ZMAX = dubla
Limita superioară a valorilor z afișate pentru graficele de contur ale densității electronilor și ale acesteia
Laplacian. Implicit este 3.0

EDGRAD_LOGSCALE = întreg
Controlează scalarea logaritmică a graficului de gradient de densitate de electroni produs. Se întoarce
scalarea dacă este setată la zero, în caz contrar, valoarea mai mare - cu atât este mai puternică
câmpul gradient va fi scalat. Valoarea recomandată (implicit) este 5.

SPIN_PROP = boolean
Semnal pentru calcularea proprietăților spin (analiza populației Mulliken alfa și beta
densități, densități de spin și contribuții anizotrope la cuplarea hiperfină
constante la centrii atomici). Implicit este fals.

PRINT = întreg
Acesta este cel mai important cuvânt cheie - determină cantitatea de informații tipărite.
În prezent sunt utilizate următoarele valori:

PRINT = 0 - modul silențios - imprimați numai rezultatele esențiale - rezultate „compacte” ale
Analiza populației Mulliken, momente electrice multipolare și electrostatică
proprietăți;

PRINT = 1 (implicit) - toate cele de mai sus plus lista de sarcini de efectuat și lista
a parametrilor de calcul;

PRINT = 2 - toate cele de mai sus plus matricea populației Mulliken AO și electronice și
componentele nucleare ale momentului dipolar electric;

PRINT = 3 - toate cele de mai sus plus matrice de densitate în baza AO și moment dipol
integrale în baza AO (și SO);

PRINT = 4 - toate cele de mai sus plus informații despre setul de bază, orbitali naturali în termeni
a orbitalilor de simetrie, matricea de suprapunere;

PRINT >= 5 - toate cele de mai sus plus vectori coeficient de cuplare, o ocupație
vector și un vector Z modificat pe baza MO.

PRINT_NR = boolean
Dacă WRTNOS = TRUE și această opțiune este, de asemenea, TRUE, orbitalii naturali vor fi imprimați
la ieșire înainte de a fi scrise în fișierul punct de control.

WRTNOS = boolean
Dacă este TRUE, orbitalii naturali vor fi scrieți în fișierul punctului de control.

GRID REZULTATE AND PLOTARE

În prezent, oeprop produce rețele bidimensionale gata pentru reprezentare cu a
program PLOTMTV versiunea 1.3.2. Programul este scris de Kenny Toh
([e-mail protejat]), dezvoltator de software pentru Departamentul de tehnologie CAD, Intel Corp,
Sfanta Clara. Este un pachet gratuit și poate fi descărcat de pe Internet.

Grilele tridimensionale sunt scoase în format potrivit pentru trasarea cu un program MegaPov
versiunea 0.5. Acest program gratuit este o versiune corectată a POV-Ray. Este dezvoltat de a
numărul de persoane și pot fi descărcate de pe Internet (accesați
http://nathan.kopp.com/patched.htm pentru a afla mai multe informații). Pentru a reda un MO sau o densitate
imagine, editați (dacă este necesar) fișierul de comandă mo.pov creat de oeprop , și executați megapovplus
+Imo.pov Pentru mai multe opțiuni, rulați megapovplus -h

Martie 30, 2001 oeprop(1)

Utilizați oeprop online folosind serviciile onworks.net