detci - En ligne dans le Cloud

PROGRAMME:

Nom

detci - Programme d'interaction de configuration de déterminant

DESCRIPTION

Le programme détective diagonalise l'opérateur hamiltonien électronique non relativiste dans un
base des déterminants de Slater. L'ensemble des déterminants utilisés (espace CI) peut être choisi dans un
de diverses façons. Le programme peut gérer n'importe quel espace CI qui peut être formulé comme un
Espace actif restreint CI. Cela inclut CISD, CISDT, CISDTQ, etc., jusqu'à Full CI, comme
ainsi que des CI multiréférences dans lesquels les références sont choisies comme tous les déterminants dans
qui jusqu'à n électrons sont excités dans un espace MO actif. Cela inclut CISD[T],
CISD[TQ] et CI de second ordre (SOCI).

Références

Espace actif restreint CI :

1. Algorithmes d'interaction de configuration basés sur des déterminants pour complet et restreint
Espaces d'interaction de configuration, J. Olsen, BO Roos, P. Jorgensen et HJ Aa.
Jensen, J. Chem. Phys. 89, 2185 (1988).

2. Passage de la limite d'un milliard dans les calculs d'interaction de configuration complète (FCI),
J. Olsen, P. Jorgensen et J. Simons, Chem. Phys. Lett. 169, 463 (1990).

Sous-espaces virtuels tertiaires (RAS IV) :

1. Fonctions d'onde variationnelles compactes incorporant un triple et un quadruple limités
Excitations, CD Sherrill et HF Schaefer, J. Phys. Chem. 100, 6069-6075
(1996).

Programme DETCI :

1. CD Sherrill, algorithmes de calcul pour des références complètes et multiples à grande échelle
Configuration Interaction Wavefunctions, thèse de doctorat, Université de Géorgie, Athènes,
AG, 1996.

DES DOSSIERS REQUIS

input.dat - Fichier d'entrée
file71 - Intégrales à un électron transformées
file72 - Intégrales à deux électrons transformées

TEMPORAIRE DES DOSSIERS OCCASION

file50 - Diagonale de l'hamiltonien
file51 - vecteurs CI
file52 - Vecteurs Sigma
file53 - fichier D (vecteurs de correction)

DES DOSSIERS ACTUALISÉ

output.dat - Fichier de sortie

CONTRIBUTION Format

Les arguments de ligne de commande suivants sont disponibles :

-silencieux Cela donne le même résultat que IMPRIMER=0.

-o fnom
Donne le nom du fichier de sortie. La valeur par défaut est output.dat.

-e Cette option entraîne l'écriture de l'énergie ou des énergies totales du CI dans un fichier appelé
detci_energies.dat.

-c Plus-value
Donne une convergence plus lâche sur le vecteur CI, utile dans les calculs DETCAS. Les
valeur est un nombre réel, pas un entier comme dans CONVERGENCE. La convergence utilisée
sera le perdant de Plus-value et CONVERGENCE.

L'entrée supplémentaire pour ce programme est lue à partir du fichier Les mots-clés suivants sont
valide:

CONVERGENCE = entier
Convergence souhaitée sur le vecteur CI. La convergence est atteinte lorsque le RMS du
l'erreur dans le vecteur CI est inférieure à 10**(-n). La valeur par défaut est 4 pour les énergies et 7
pour les dégradés. Ce n'est pas le même critère de convergence du vecteur CI que celui trouvé dans
GUGACI.

DOCC = tableau_entiers
Ce vecteur donne le nombre d'orbitales doublement occupées dans chaque irrep. Il y a
aucun défaut.

SOCC = tableau_entiers
Ce vecteur donne le nombre d'orbitales occupées individuellement dans chaque irrep. Il y a
aucun défaut.

DIAG_METHOD = un magnifique
Ceci spécifie quelle méthode doit être utilisée pour diagonaliser l'hamiltonien. Les
les options valides sont: RER, pour former la matrice H entière et diagonaliser en utilisant libciomr
pour obtenir toutes les valeurs propres (nb nécessite une mémoire ÉNORME); OLSEN, pour utiliser celle d'Olsen
méthode du sous-espace inverse préconditionné (1990); MITROUCHENKOV, pour utiliser un 2x2
méthode Olsen/Davidson ; et DAVIDSON (ou SEM) pour utiliser l'extension simultanée de Liu
Méthode, qui est identique à la méthode de Davidson si l'on ne trouve qu'une seule racine.
Il existe également un mode de débogage SEM, SEMTESTL’ SEM la méthode est la plus
robuste, mais il nécessite également 2(N*M)+1 vecteurs CI sur le disque, où N est le maximum
nombre d'itérations et M est le nombre de racines.

PRÉCONDITIONNEUR = un magnifique
Ceci spécifie le type de préconditionneur à utiliser dans la diagonalisation sélectionnée
méthode. Les options valides sont : DAVIDSON qui se rapproche de la matrice hamiltonienne
par les éléments diagonaux ; H0BLOC_INV qui utilise un hamiltonien exact de
H0_BLOCKSIZE et l'inverse explicitement ; GEN_DAVIDSON qui fait une spectrale
décomposition de H0BLOCK ; ITER_INV en utilisant une approche itérative pour obtenir le
vecteur de correction de H0BLOCK. Le H0BLOCK_INV, GEN_DAVIDSON et ITER_INV
les approches sont toutes formellement équivalentes mais l'ITER_INV est moins calculatoire
cher. La valeur par défaut est DAVIDSON.

RÉFÉRENCE = un magnifique
Ceci spécifie le type de fonction de référence. C'est RHF ou ROHF. UHF et
TWOCON ne sont pas pris en charge. Pour ROHF, une multiplicité de 1 implique une coquille ouverte
maillot. Le programme fonctionnera pour les singlets open-shell, mais il n'a pas été
correctement adapté pour utiliser une référence à deux déterminants correcte dans ce cas, donc
l'exécution avec des références singulet open-shell n'est pas conseillée, sauf pour les CI complets.

MISE À JOUR = un magnifique
DAVIDSON utilise la formule standard de vecteur de mise à jour ou de correction DAVIDSON, tandis que
OLSEN utilise le vecteur de correction OLSEN. La valeur par défaut est DAVIDSON.

HD_OTF = booléen
Si VRAI les éléments diagonaux de la matrice hamiltonienne sont calculés à la volée,
sinon, un vecteur d'élément diagonal est écrit dans un fichier séparé sur le disque. Défaut
est vrai.

HD_AVE = un magnifique
HD_EXACT utilise les énergies diagonales exactes qui se traduisent par des vecteurs d'expansion qui
briser la symétrie de spin. HD_KAVE fait la moyenne des énergies diagonales sur un couplage de spin
ensemble produisant des vecteurs d'expansion purs de spin. ORB_ENER utilise la somme des orbitales
approximation de l'énergie donnant des vecteurs d'expansion de spin purs mais double généralement le
nombre d'itérations de davidson. ÉVANGÉLISTI utilise les sommes et les différences des orbitales
énergies avec l'énergie de référence SCF pour produire des vecteurs d'expansion de spin pur.
LEININGER approximation qui soustrait la contribution d'un électron de la
énergies orbitales, multiplie par 0.5 et ajoute la contribution d'un électron
dans, produisant des vecteurs d'expansion de spin pur et développé par votre serviteur et fonctionne comme
ainsi que ÉVANGÉLISTI.

FICHIER NOD = booléen
Uniquement possible si NUM_ROOTS = 1. Utilise le dernier espace vectoriel du fichier BVEC pour
écrivez scratch DVEC plutôt que d'utiliser un fichier DVEC séparé.

ENERGIE_CONVERGENCE = entier
Convergence souhaitée sur l'énergie CI. La valeur par défaut est 6 pour les énergies de point unique
et 8 pour les dégradés ou CASSCF.

EX_LVL = entier
Niveau d'excitation pour les excitations dans les orbitales virtuelles (par défaut 2, c'est-à-dire CISD).

VAL_EX_LVL = entier
Niveau d'excitation pour les références dans les orbitales de RAS II. La valeur par défaut est zéro.

FCI = booléen
Si cet indicateur est défini sur VRAI, alors le stockage des chaînes est simplifié pour un Full
CI et le calcul nécessite moins de frais généraux. Cependant, les résultats finaux devraient
être identiques à ceux lorsque FCI = FAUX. Peut provoquer des résultats imprévisibles si FCI =
VRAI mais EX_LVL n'est pas compatible avec un IC complet.

CONGELÉ_DOCC = tableau_entiers
Le nombre d'orbitales les plus énergétiques doublement occupées dans chaque irréductible
représentation à partir de laquelle il n'y aura pas d'excitations. La commande de coton de la
représentations irréductibles est utilisé. La valeur par défaut est le vecteur zéro.

CONGELÉ_UOCC = vecteur_entier
Le nombre d'orbitales inoccupées les plus énergétiques dans chaque représentation irréductible
dans lequel il n'y aura pas d'excitations. La valeur par défaut est le vecteur zéro.

FREEZE_CORE = booléen
Cette option détermine si les orbitales de noyau gelées doivent être incluses
implicitement (vrai) ou explicitement (faux). Dans le premier cas, l'énergie
les contributions des orbitales de noyau gelées sont repliées dans l'un électron
intégrales et dans "l'énergie du noyau gelé" calculée par le programme de transformation.
La valeur par défaut est true.

EXPORT_VECTOR = booléen
Ceci spécifie s'il faut stocker le(s) vecteur(s) convergé(s) à la fin de l'analyse. Les
le(s) vecteur(s) est(sont) stocké(s) dans un format transparent tel que d'autres programmes puissent l'utiliser
facilement. Le format est spécifié dans src/lib/libqt/slaterdset.h. La valeur par défaut est
faux

NUM_EXPORTATION = entier
If EXPORT_VECTOR est défini sur vrai, ceci détermine le nombre de vecteurs qui
doivent être exportés à la fin de l'exécution. La valeur par défaut est 1.

GUESS_VECTOR = un magnifique
Ceci spécifie le type de vecteur de supposition à utiliser dans l'itération CI. Actuellement
uniquement utilisé par la méthode d'itération SEM. Les valeurs acceptées sont UNITÉ pour un vecteur unitaire
deviner (NUM_ROOTS et NUM_INIT_VECS doivent tous les deux être 1); H0_BLOC utiliser des vecteurs propres
à partir de la sous-matrice H0 BLOCK (par défaut) ; DFILE à utiliser NUM_ROOTS auparavant convergé
vecteurs dans le fichier D ; et MP2 pour utiliser la fonction d'onde MP2 comme une supposition (ne fonctionne pas
à l'heure actuelle).

H0_BLOCKSIZE = entier
Ce paramètre spécifie la taille du bloc "H0" de l'hamiltonien qui est
résolu exactement. Les n déterminants avec la plus faible énergie SCF sont sélectionnés, et un
la sous-matrice de l'hamiltonien est formée à l'aide de ces déterminants. Cette sous-matrice est
utilisé pour accélérer la convergence des itérations de CI dans le OLSEN et MITROUCHENKOV
schémas d'itération, et aussi pour trouver une bonne estimation de départ pour le SEM méthode si
GUESS_VECTOR = H0_BLOC. La valeur par défaut est 40. Notez que le programme peut changer le
taille donnée pour Ms=0 cas (Ms0 = VRAI) s'il détermine que le bloc H0 comprend
un seul membre d'une paire de déterminants liés par la symétrie d'inversion temporelle. Pour
de très petites tailles de blocs, cela pourrait éventuellement éliminer tout le bloc H0 ; les
le programme devrait imprimer des avertissements si cela se produit.

H0_BLOCK_COUPLING_SIZE = entier
Paramètres qui spécifient la taille du bloc de couplage dans le
préconditionneur davidson. La valeur par défaut est 1000.

MAX_DET = entier
Définit le nombre maximum de déterminants ; si l'espace CI est plus grand que cela, le
le programme s'interrompt. Cette option existe pour s'assurer que les calculs très volumineux ne sont pas
courir par accident. Au cours de la phase de développement actuelle, la valeur par défaut est 10000, mais
il sera levé d'ici peu.

MAXITRE = entier
Nombre maximum d'itérations pour diagonaliser l'hamiltonien. La valeur par défaut est 12.

Ms0 = booléen
If VRAI, utilisez la composante Ms=0 de l'état. Par défaut à VRAI si coque fermée et
à FAUX autrement. En lien avec S paramètre.

NIMPRIMER = entier
Cette valeur spécifie le nombre de déterminants qui seront imprimés, ainsi que
leurs coefficients, dans la liste des déterminants les plus importants dans l'IC final
vecteur. La valeur par défaut est 20.

NUM_ROOTS = entier
Cette valeur donne le nombre de racines qui doivent être obtenues à partir du séculaire
équations. La valeur par défaut est un. Si plus d'une racine est requise, définissez DIAG_METHOD
à SEM (ou, pour les très petits cas, RER or SEMTEST).

NUM_INIT_VECS = entier
Le nombre de vecteurs initiaux à utiliser dans la procédure itérative CI. Par défaut à
le nombre de racines.

OPDM = booléen
If VRAI calculer la matrice de densité à une particule et faire OPDM_WRITE par défaut
VRAI. La valeur par défaut de OPDM is FAUX.

OPDM_FILE = entier
Fichier (numéro d'unité) pour l'écriture de la matrice de densité à une particule si OPDM_WRITE =
VRAI. La valeur par défaut est actuellement 73.

OPDM_WRITE = booléen
Indicateur indiquant s'il faut ou non écrire la matrice de densité à une particule sur le disque.

OPDM_PRINT = booléen
Indicateur d'impression ou non de la matrice de densité à une particule.

OPDM_DIAG = booléen
Indicateur indiquant s'il faut ou non diagonaliser la matrice de densité à une particule.

WRTNOS = booléen
Indicateur d'écriture ou non des orbitales naturelles CI dans PSIF_CHKPT.

FICHIER ORBS = entier
Fichier (numéro d'unité) pour l'écriture de diverses orbitales naturelles CI. La valeur par défaut est
76.

OPDM_AVE = booléen
Indicateur de moyenne ou non de l'OPDM sur plusieurs racines afin d'obtenir un
matrice de densité à une particule moyenne d'état. Cette matrice de densité peut être diagonalisée
pour obtenir les orbitales naturelles CI.

ORBS_ROOT = entier
Drapeau pour définir le nombre de racine pour lequel les orbitales naturelles CI sont écrites
PSIF_CHKPT. La valeur par défaut est 1 (racine la plus basse).

Print = entier
Cette option détermine la verbosité de la sortie. Une valeur de 1 ou 2 spécifie
impression minimale, une valeur de 3 spécifie une impression détaillée. Les valeurs de 4 ou 5 sont
utilisé pour le débogage. N'utilisez pas le niveau 5 sauf si le cas de test est très petit (par ex.
DSIC STO H2O).

TRAITEMENT = entier
La racine pour laquelle écrire la matrice de densité à deux particules (la densité à une particule
les matrices sont écrites pour toutes les racines). Utile pour un CASSCF ou un CI spécifique à un état
optimisation sur un état excité.

S = entier
La valeur du nombre quantique de spin S est donnée par cette option. La valeur par défaut est 0
(maillot). La seule chose à laquelle cela est réellement utilisé est de déterminer la phase de
la moitié redondante du vecteur CI lorsque la composante Ms=0 est utilisée (c'est-à-dire, Ms0 =
VRAI). Pour les cas où S n'est pas un entier, ce paramètre n'a pas besoin d'être saisi
car un tel état ne peut pas avoir une composante Ms=0.

TPDM = booléen
If VRAI calculer la matrice de densité à deux particules et faire TPDM_WRITE par défaut
VRAI. La valeur par défaut de TPDM is FAUX.

TPDM_FILE = entier
Fichier (numéro d'unité) pour l'écriture de la matrice de densité à deux particules si TPDM_WRITE =
VRAI. La valeur par défaut est actuellement 74.

TPDM_WRITE = booléen
Indicateur indiquant s'il faut ou non écrire la matrice de densité à deux particules sur le disque.

TPDM_PRINT = booléen
Indicateur d'impression ou non de la matrice de densité à deux particules. Typiquement un très
mauvaise idée sauf pour le débogage de petits cas.

Il existe également des entrées moins couramment utilisées, que les utilisateurs novices de PSI n'auront pas besoin de
utiliser.

BENDAZZOLI = booléen
Utilisez des routines pour calculer le sigma sur la base des articles de Bendazzoli et al.
Semble être plus lent et ne vaut pas la peine ; peut finir par disparaître. Fonctionne uniquement pour
IC complet et je ne me souviens pas si je pouvais voir à quel point leur schéma intelligent pourrait être
étendu au RAS en général.

CALC_SSQ = booléen
Si VRAI, calculez la valeur attendue de l'opérateur S^2 pour le CI final
fonction d'onde pour chaque racine. En principe, DETCI devrait fournir S^2 fonctions propres.
Le défaut est faux.

COLLAPSE_SIZE entier
Donne le nombre de vecteurs à retenir lorsque le sous-espace de Davidson est réduit (voir
MAXNVECT au dessous de). S'il est supérieur à un, le sous-espace réduit conserve le meilleur
estimation du vecteur CI pour les n itérations précédentes. La valeur par défaut est 1.

FIRST_TMP_UNIT = entier
Donne le numéro de fichier (unité) associé au premier fichier scratch utilisé par DETCI.
Les autres fichiers de travail sont numérotés consécutivement à partir de ce point, dans l'ordre
H(diag), C, S, D. Chacun de ces fichiers logiques occupe un certain nombre de fichiers physiques
spécifié par les paramètres d'entrée encore plus obscurs NUM_HD_TMP_UNITS,
NUM_C_TMP_UNITS, NUM_S_TMP_UNITS, NUM_D_TMP_UNITS. L'utilisateur peut également spécifier
différents points de départ pour chacun de ces ensembles en utilisant les paramètres
FIRST_HD_TMP_UNIT et ainsi de suite. Le fractionnement d'un fichier sur plusieurs unités peut aider
éviter le problème de taille d'entier dans le traitement de fichiers volumineux présents dans
DETCI et dans les bibliothèques d'E/S PSI ; mais encore une fois, je ne l'ai pas testé pour voir ce que
arrive. La première unité de chaque section est imprimée sous la rubrique FICHIERS dans
la sortie de paramètre commençant l'exécution de DETCI.

FZC = booléen
Détermine si les orbitales centrales gelées sont traitées comme vraiment gelées (c'est-à-dire,
totalement absent du calcul, FZC = VRAI) ou s'ils sont présents mais
restreint à être doublement occupé (FZC = FAUX). Dans le programme GUGA CI, c'est
la distinction entre ce qu'il appelle les orbitales FZC et COR. Généralement, le
les intégrales pour les orbitales de noyau gelées ne sont pas nécessaires par DETCI mais elles peuvent être nécessaires
pour MCSCF ou dégradés.

ICORE = entier
Spécifie comment gérer la mise en mémoire tampon des vecteurs CI. Une valeur de 0 rend le programme
effectuer des E/S un sous-bloc RAS à la fois ; 1 utilise des vecteurs CI entiers à la fois ; et 2
utilise un bloc irrep à la fois. Les valeurs de 0 ou 2 provoquent une certaine inefficacité dans le
E/S (nécessitant plusieurs lectures du vecteur C lors de la construction de H dans l'itération
sous-espace si DIAG_METHOD = SEM), mais nécessite moins de mémoire centrale.

J'ARRÊTE = booléen
If VRAI alors DETCI s'arrêtera après la formation des informations de chaîne et avant
les intégrales sont lues. Peut éventuellement changer en un entier afin que l'utilisateur puisse
choisir parmi plusieurs points d'arrêt.

MAXNVECT = entier
Donne le nombre maximum de vecteurs de sous-espace de Davidson qui peuvent être conservés sur le disque pendant
le coefficient CI et les vecteurs sigma. (Il y a un vecteur H(diag) et le nombre
de D vecteurs est égal au nombre de racines). Lorsque le nombre de vecteurs sur le disque
atteint la valeur de MAXNVECT, le sous-espace de Davidson sera réduit à
COLLAPSE_SIZE vecteurs pour chaque racine. Ceci est très utile pour économiser de l'espace disque.
La valeur par défaut est MAXITRE * NUM_ROOTS + NUM_INIT_VECS.

MIXTE = booléen
Ceci détermine si les excitations « mixtes » RAS II/RAS III sont autorisées dans le CI
espacer. Ceci est utile pour placer des contraintes supplémentaires sur un CI RAS.

MIXTE4 = booléen
Ceci est similaire au mot clé MIXED, mais fait référence aux excitations dans RAS IV.

NUITS = entier
Nombre de fichiers de travail à utiliser pour stocker les vecteurs C (et aussi pour le sigma
vecteurs).

OEI_ERASE = booléen
Ceci détermine si le programme efface le fichier d'intégrales à un électron après
a été lu. La valeur par défaut finira par être vraie, mais pendant le développement, le
la valeur par défaut est faux.

OEI_FILE = entier
Ce mot-clé permet à l'utilisateur de spécifier le fichier intégral à un électron transformé.
La valeur par défaut est 71.

PRINT_CIBLKS = booléen
Indique si le programme doit imprimer un résumé de tous les blocs du
Vecteur CI (qui peut être converti sous forme matricielle, voir les références.)

R4S = booléen
Restreint les chaînes RAS IV à l'ensemble minimal, économisant ainsi de la mémoire. Si vous êtes
préoccupé par cette option, vous devriez écrire à David pour obtenir des conseils, sauf si vous êtes un
Expert DETCI.

REF_SYM = entier
Cette option permet à l'utilisateur de rechercher des vecteurs CI d'un irrép différent que le
référence. Cela n'a probablement de sens que pour Full CI, et ce ne serait probablement pas
travailler avec des suppositions vectorielles unitaires. La numérotation commence à partir de zéro pour le total-
irrép symétrique.

REPL_OTF = booléen
Indique à DETCI s'il faut ou non effectuer des remplacements de chaîne à la volée. Peut sauver un
quantité de mémoire gigantesque (en particulier pour les CI tronqués) mais est quelque peu floconneuse et
n'a pas été testé depuis un moment. Si je me souviens bien, cela ne fonctionne que pour certaines classes de
calculs RAS. Contactez David pour obtenir de l'aide. Finalement, le vol
les éléments de remplacement doivent être refaits de manière beaucoup plus intelligente afin que cela ne prenne pas
des éons de temps CPU. Les travaux dans ce sens ont déjà commencé et peuvent être achevés
par la suite.

RESTART = booléen
Cette option permet à l'utilisateur de reprendre une itération DETCI qui s'est terminée
prématurément. Il suppose que les vecteurs CI et sigma sont sur le disque ; le nombre de
vecteurs spécifiés par RESTART_VECS est réduit à un vecteur par racine.

RESTART_VECS = entier
If RESTART = VRAI ceci spécifie le nombre de vecteurs CI (et sigma) à lire à partir
disque. Il s'agit généralement du nombre d'itérations terminées avec succès à partir d'un
l'exécution précédente multiplie le nombre de racines pour cette exécution.

TEI_ERASE = booléen
Ceci détermine si le programme efface le fichier d'intégrales à deux électrons après
a été lu. La valeur par défaut finira par être vraie, mais pendant le développement, le
la valeur par défaut est faux.

TEI_FILE = entier
Ce mot-clé permet à l'utilisateur de spécifier le fichier intégral à deux électrons transformé.
La valeur par défaut est 72.

MPN = booléen
Lorsque cette option est VRAIE, DETCI calculera la série MPn jusqu'au kième ordre où k
est déterminé par maxnvect. Pour les systèmes à coque ouverte (REF=ROHF, WFN = ZAPTN), DETCI
calculera la série ZAPTn. GUESS_VECTOR doit être défini sur TRUE. HD_OTF doit être
réglé sur VRAI. HD_AVE doit être défini sur orb_ener.

SAVE_MPN2 = entier
Lorsque MPN est VRAI et WIGNER est VRAI, cette option devient valide. Si réglé sur 1
alors l'énergie MP(2n-1) est économisée. S'il est réglé sur 2, l'énergie MP(2n-2) est économisée. Si seulement
l'autre valeur MPn énergie est économisée. La valeur par défaut est 0.

9 fév, 1996 détective(1)

Utilisez detci en ligne en utilisant les services onworks.net