oeprop - En ligne dans le Cloud

PROGRAMME:

Nom

oeprop - Programme de propriété à un électron

DESCRIPTION

Le programme oeprop calcule les valeurs attendues des opérateurs de propriété à un électron en utilisant un
matrice de densité à une particule calculée à partir d'un vecteur propre dans PSIF_CHKPT ou lue à partir d'un
fichier externe. Il est actuellement capable d'effectuer des analyses de population Mulliken,
calcul des moments multipolaires électriques à travers l'octopole, propriétés électrostatiques atomiques
centres (potentiel électrostatique, champ électrique, gradient de champ électrique, électron et
densité de spin, contribution anisotrope dipolaire aux constantes de couplage hyperfin),
densité d'électrons et de spins, gradient de densité d'électrons et de spins, laplacien d'électrons et
densités de spin, potentiel électrostatique sur un bidimensionnel arbitraire (planaire)
grille rectangulaire et valeurs d'orbitales moléculaires sur une échelle tridimensionnelle arbitraire
grille rectangulaire. Les capacités diverses incluent le calcul de la relation relativiste
corrections d'un électron de premier ordre à l'énergie (termes masse-vitesse et Darwin),
construction d'orbitales moléculaires naturelles à partir d'une densité de particule lue à partir d'un
fichier (les NO peuvent être écrits dans PSIF_CHKPT) et calcul des étendues spatiales - attente
valeurs des opérateurs X^2, Y^2, Z^2 et R^2 - de la densité électronique totale et de l'individu
MO (si LIRE_OPDM = faux) ou naturel (si LIRE_OPDM = vrai) orbitales (MPMAX doit être réglé sur
une valeur supérieure à 1 pour le calcul de ces entités). Les étendues spatiales doivent être utilisées
avec prudence, car ils dépendent du point de référence.

Références

Analyse de la population de Mulliken

1. Analyse électronique de la population sur les fonctions d'onde moléculaire LCAO-MO. RS
Mulliken, J. Chem. Phys. 23, 1833 (1955), idem. 23, 1841 (1955), idem. 36, 3428
(1962).

Relations de récurrence pour les intégrales à un électron sur les fonctions cartésiennes gaussiennes.

1. Calcul récursif efficace des intégrales moléculaires sur la gaussienne cartésienne
les fonctions. S. Obara et A. Saika, J. Phys. Chem. 84, 3963 (1986).

Constantes physiques fondamentales et facteurs de conversion.

1. CRC Manuel de chimie et de physique. Edité par DR Lide. 73e édition
(1992-1993).

DES DOSSIERS REQUIS

input.dat - Fichier d'entrée
PSIF_CHKPT - Fichier de point de contrôle

DES DOSSIERS ACTUALISÉ

sortie.dat
dipmom.dat - Moments dipolaires
esp.dat - Potentiel électrostatique sur une grille 2D
edens.dat - Densité électronique sur une grille 2D
edgrad.dat - Gradient de densité électronique sur une grille 2D
edlapl.dat - Laplacien de la densité électronique sur une grille 2D
sdens.dat - Densité de spin sur une grille 2D
sdgrad.dat - Gradient de densité de rotation sur une grille 2D
sdlapl.dat - Laplacien de la densité de spin sur une grille 2D
mo.dat - Valeurs orbitales moléculaires/densité sur une grille 3D
mo.pov - Fichier d'entrée MegaPov pour le rendu d'une image de mo.dat
mo.cube - Orbitale(s) moléculaire(s) sur une grille 3D au format Cube Gaussien94
dens.cube - Densité(s) d'électrons/spins sur une grille 3D au format Cube Gaussien94

CONTRIBUTION Format

La plupart des mots-clés ne sont pas nécessaires pour les tâches de routine. Les mots-clés suivants sont
valide:

WFNPlus = booléen
Type de la fonction d'onde. Ce mot-clé est une "macro" qui permet à l'utilisateur de définir la plupart des
les mots-clés nécessaires. Les valeurs suivantes sont reconnues :

WFNPlus = SCF - équivalent à LIRE_OPDM = faux;

WFNPlus = DETCI - équivalent à LIRE_OPDM = vrai, OPDM_FILE = 40 OPDM_BASE = AO,
OPDM_FORMAT = TRIANG ;

WFNPlus = CCDS - équivalent à EAD_OPDM = vrai, OPDM_FILE = 79 OPDM_BASE = AO,
OPDM_FORMAT = TRIANG ;

WFNPlus = QVCCD - équivalent à LIRE_OPDM = vrai, OPDM_FILE = 76 OPDM_BASE = SO,
OPDM_FORMAT = TRIANG ;

LIRE_OPDM = booléen
Cet indicateur spécifie si la matrice de densité à une particule doit être lue à partir du disque.
La valeur par défaut est false.

OPDM_FILE = entier
Spécifie le numéro de fichier de la matrice de densité à une particule. La valeur par défaut est 40 (fichier maître). À
fournir une compatibilité descendante avec les packages de propriétés PSI antérieurs (correct,
ciprop, ccprop) le format spécial du fichier de densité est supposé lorsque OPDM_FILE = 40
(calcul des propriétés à partir de la densité CI - ciprop mode de compatibilité) et OPDM_FILE =
79 (calcul des propriétés à partir de la densité CC - ccprop le mode de compatibilité). A partir de maintenant,
dans le cas générique, onepdm doit être écrit au tout début du fichier. Dans le
le futur PSI aura un fichier onepdm standard.

OPDM_BASE = un magnifique
Cette option peut ne plus exister à l'avenir. Depuis le 1er février 1998, une norme pour
le format de fichier onepdm n'a pas été défini. Ce mot-clé doit être défini sur "SO"
(lu dans la matrice onepdm en base SO) ou "AO" (en base AO). La valeur par défaut est « SO ».

OPDM_FORMAT = un magnifique
Cette option peut ne plus exister à l'avenir. Ce mot-clé doit être défini sur
"TRIANG" (lu dans la matrice onepdm sous forme triangulaire inférieure) de "SQUARE" (en carré
former). La valeur par défaut est "TRIANG"

ASYMM_OPDM = booléen
Ce drapeau spécifie si la matrice de densité à une particule doit être symétrisée.
Doit être défini sur true si un onepdm générique non symétrique doit être lu (par exemple, à partir d'un
programme de cluster couplé). Ce mot-clé est uniquement destiné au développement de code. PSI existant
Les codes CC actuellement utilisés produisent un onepdm symétrique, il n'est donc pas nécessaire d'utiliser
ce mot-clé. La valeur par défaut est false.

TRAITEMENT = entier
Ceci spécifie pour quelle racine effectuer l'analyse de l'état excité. Approprié
une matrice de densité de particules sera lue à partir du disque. Actuellement mis en œuvre pour
Fonctions d'onde DETCI et DETCAS.

MPMAX = entier
Cet entier compris entre 1 et 3 spécifie le moment multipolaire électrique le plus élevé à être
calculé.

MPMAX = 1 - seul le moment dipolaire électrique sera calculé (par défaut) ;

MPMAX = 2 - les moments dipolaires et quadripolaires électriques seront calculés ; MPMAX = 3 -
Les moments électriques dipolaires, quadripolaires et octopolaires seront calculés.

MP_REF entier
Ce paramètre spécifie le point de référence pour les moments multipolaires électriques
calcul.

MP_REF = 0 (par défaut) ou 1 - le centre de gravité ;

MP_REF = 2 - l'origine du système de coordonnées spatiales;

MP_REF = 3 - le centre de charge électronique ;

MP_REF = 4 - le centre de charge nucléaire ;

MP_REF = 5 - le centre de charge net.

ATTENTION : D'après l'électrodynamique classique, le moment électrique à 2^(n+1) pôles
n'est indépendant du point de référence que si le moment électrique 2^(n)-pôle est
disparition. Cela signifie que le moment dipolaire dépendra du point de référence si
la charge totale du système est non nulle. Par analogie, moment quadripolaire électrique
dépendra du point de référence si le système possède un dipôle électrique non nul
instant, etc

MP_REF_XYZ = vecteur_réel
Ce vecteur spécifie les coordonnées du point de référence. Si ce mot-clé est
présent dans l'entrée MP_REF le mot-clé sera ignoré.

NUC_ESP = booléen
Ce drapeau spécifie si les propriétés électrostatiques seront calculées au niveau des noyaux.
La liste actuelle comprend le potentiel électrostatique, le champ électrique, le champ électrique
gradient, densité d'électrons et de spin, et contribution anisotrope à l'hyperfin
constantes de couplage (les deux dernières nécessitent un réglage SPIN_PROP à vrai). La valeur par défaut est
vrai.

GRID = entier
Spécifie le type de propriété à évaluer sur une grille.

GRID = 0 (par défaut) - ne rien calculer ;

GRID = 1 - potentiel électrostatique sur une grille à deux dimensions;

GRID = 2 - densité électronique (densité de spin si SPIN_PROP est défini sur vrai) sur deux
grille dimensionnelle;

GRID = 3 - gradient de densité électronique (gradient de densité de spin si SPIN_PROP est fixé à
true) sur une grille à deux dimensions ;

GRID = 4 - Laplacien de la densité électronique (Laplacien de la densité de spin si
SPIN_PROP est défini sur vrai) sur une grille à deux dimensions. Selon la convention
utilisé sur le terrain, ce qui est réellement tracé, ce sont les Laplaciens pris avec
signe négatif.

GRID = 5 - valeurs des orbitales moléculaires sur une grille tridimensionnelle.

GRID = 6 - valeurs de la densité électronique (gradient de densité de spin si SPIN_PROP is
mis à vrai) sur une grille en trois dimensions.

GRID_FORMAT = un magnifique
Spécifie dans quel format la sortie de la grille sera produite. Actuellement, PLOTTMTV
(par défaut pour les grilles 2D), MÉGAPOVPLUS (disponible pour les grilles 3D), et
GAUSSCUBE(par défaut pour les grilles 3D) sont pris en charge.

MO_TO_PLOT = vecteur
Spécifie les indices des orbitales moléculaires à calculer sur la grille 3D. Indices
peut être spécifié comme :

entier non signé - index dans l'ordre de Pitzer (ordonné selon irreps, pas
valeurs propres). Va de 1 au nombre de MO.

entier signé - indice par rapport au niveau de Fermi. +1 signifie LUMO, +2 signifie deuxième
orbitale virtuelle la plus basse, -1 signifie HOMO, etc.

Tous les indices doivent être soit non signés, soit signés, vous ne pouvez pas mélanger et assortir, ou vous
obtiendra des résultats imprévisibles. La valeur par défaut est de calculer HOMO et LUMO.

GRID_ORIGINE = vecteur_réel
Spécifie l'origine de la grille. Un quadrillage rectangulaire qui enveloppe l'ensemble
molécule sera calculée automatiquement si GRID_ORIGINE manque, cependant, il
n'est pas défini par défaut pour les grilles 2D.

GRID_UNIT_X = vecteur_réel
Ce vecteur spécifie la direction du premier (x) côté de la grille. Ce n'est pas
doivent être de longueur unitaire. Il n'y a pas de valeur par défaut pour les grilles 2D.

GRID_UNIT_Y = vecteur_réel
De même pour le deuxième (y) côté. Il n'a pas besoin d'être de longueur unitaire ou même
orthogonal à GRID_UNIT_X. Il n'y a pas de valeur par défaut pour les grilles 2D.

GRILLE_XY0 = real_2d_vecteur
Spécifie les coordonnées du coin inférieur gauche du rectangle de la grille dans le 2D
système de coordonnées défini par GRID_ORIGIN, GRID_UNIT_X et GRID_UNIT_Y. Il y a
aucun défaut.

GRILLE_XY1 = real_2d_vecteur
Spécifie les coordonnées du coin supérieur droit du rectangle de la grille en 2D
système de coordonnées défini par GRID_ORIGIN, GRID_UNIT_X et GRID_UNIT_Y. Il y a
aucun défaut.

GRID_XYZ0 = real_3d_vecteur
Spécifie les coordonnées du coin inférieur gauche de la boîte de grille dans la 3D
système de coordonnées défini par GRID_ORIGIN, GRID_UNIT_X, GRID_UNIT_Y et la croix
produit des deux derniers. Il n'y a pas de défaut.

GRID_XYZ1 = real_3d_vecteur
Spécifie les coordonnées du coin supérieur droit proche de la boîte de grille dans la 3D
système de coordonnées défini par GRID_ORIGIN, GRID_UNIT_X, GRID_UNIT_Y et la croix
produit des deux derniers. Il n'y a pas de défaut.

NIX = entier
Le nombre de points de grille le long de la direction x. Ce paramètre doit être supérieur à
1. La valeur par défaut est 20.

NI = entier
Le même que NIX pour la direction y. La valeur par défaut est 20.

CHAÎNE DE CARACTÈRES = entier
Le même que NIX pour la direction z. La valeur par défaut est 20.

GRID_ZMIN = double
Limite inférieure des valeurs z affichées pour les tracés de contour de la densité électronique et de ses
Laplacien. La valeur par défaut est 0.0

GRID_ZMAX = double
Limite supérieure des valeurs z affichées pour les tracés de contour de la densité électronique et de ses
Laplacien. La valeur par défaut est 3.0

EDGRAD_LOGSCALE = entier
Contrôle la mise à l'échelle logarithmique du tracé du gradient de densité électronique produit. Se tourne
la mise à l'échelle s'il est défini sur zéro, sinon la valeur la plus élevée - plus le
le champ de gradient sera mis à l'échelle. La valeur recommandée (par défaut) est 5.

SPIN_PROP = booléen
Indicateur pour le calcul des propriétés de spin (analyse de population Mulliken d'alpha et bêta
densités, densités de spin et contributions anisotropes au couplage hyperfin
constantes aux centres atomiques). La valeur par défaut est false.

Print = entier
C'est le mot-clé le plus important - il détermine la quantité d'informations imprimées.
Les valeurs suivantes sont actuellement utilisées :

Print = 0 - mode silencieux - imprimer uniquement les résultats essentiels - résultats "compacts" de
Analyse de population de Mulliken, moments multipolaires électriques et électrostatique
Propriétés;

Print = 1 (par défaut) - tout ce qui précède plus la liste des tâches à effectuer et la liste
des paramètres de calcul ;

Print = 2 - tout ce qui précède plus la matrice de population Mulliken AO et électronique et
composants nucléaires du moment dipolaire électrique;

Print = 3 - tout ce qui précède plus la matrice de densité en base AO et moment dipolaire
intégrales en base AO (et SO);

Print = 4 - tous les éléments ci-dessus plus les informations de base, les orbitales naturelles en termes
des orbitales de symétrie, matrice de chevauchement ;

Print >= 5 - tout ce qui précède plus les vecteurs de coefficients de couplage, une occupation
vecteur, et un vecteur Z modifié en base MO.

PRINT_NOS = booléen
Si WRTNOS = TRUE et que cette option est également TRUE, les orbitales naturelles seront imprimées
à sortir avant qu'ils ne soient écrits dans le fichier de point de contrôle.

WRTNOS = booléen
Si VRAI, les orbitales naturelles seront écrites dans le fichier de point de contrôle.

GRID SORTIE ET TRAÇAGE

À l’heure actuelle, oeprop produit une sortie de grilles bidimensionnelles prêtes à être tracées avec un
Programme PLOTTMTV version 1.3.2. Le programme est écrit par Kenny Toh
([email protected]), développeur de logiciels pour le département Technology CAD, Intel Corp,
Sainte Claire. Il s'agit d'un logiciel gratuit et peut être téléchargé sur Internet.

Les grilles tridimensionnelles sont sorties dans un format adapté au traçage avec un programme MégaPov
version 0.5. Ce programme gratuit est une version corrigée de POV-Ray. Il est développé par un
nombre de personnes et peut être téléchargé sur Internet (allez à
http://nathan.kopp.com/patched.htm pour en savoir plus). Pour rendre un MO ou une densité
image, éditer (si nécessaire) le fichier de commandes mo.pov créé par oeprop , et exécuter mégapovplus
+Moi.pov Pour plus d'options, exécutez mégapovplus -h

30 mars 2001 oeprop(1)

Utiliser oeprop en ligne en utilisant les services onworks.net