detcas - En ligne dans le Cloud

PROGRAMME:

Nom

detcas - Interaction de configuration du déterminant / Actif complet
Programme de champs cohérents dans l'espace

DESCRIPTION

Le programme détcas effectue un champ auto-cohérent d'espace actif complet (CASSCF)
optimisation des orbitales moléculaires via une procédure en deux étapes en conjonction avec le
programme d'interaction de configuration déterminante DETCI. Le programme est assez simple et
utilise actuellement une approche Newton-Raphson pour mettre à jour les orbitales, en utilisant un simple
Hessian orbital approximatif. La convergence est accélérée en utilisant l'inversion directe de Pulay de
la procédure de sous-espace itérative (DIIS). Le code a été écrit pour permettre une utilisation plus générale
des fonctions d'onde qui ne comportent pas nécessairement un traitement CI complet de l'espace actif. Dans
En particulier, toute fonction d'onde CI d'espace actif restreint (RAS) prise en charge par DETCI peut être
utilisé, permettant des fonctions d'onde MCSCF de type RASSCF.

Références

Hesse orbitale approximative :

1. G. Chaban, MW Schmidt, et MS Gordon, Theor. Chim. Actes 97, 88-95 (1997).

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).

Espace actif restreint SCF :

1. P.-A. Malmqvist, A. Rendell et BO Roos, J. Phys. Chem. 94, 5477 (1990).

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.

2. CD Sherrill et HF Schaefer, La méthode d'interaction de configuration : avancées
dans Approches hautement corrélées, Adv. Chimie quantique. 34, 143-269 (1999).

DES DOSSIERS EST REQUIS

input.dat - Fichier d'entrée
file78 - MO intégrales à un électron (opérateur fzc)
file72 - MO intégrales à deux électrons
file73 - Matrice de densité à une particule MO
file74 - Matrice de densité à deux particules MO
file75 - MO Lagrangien

DES DOSSIERS ACTUALISÉ

output.dat - Fichier de sortie
file14.dat - Enregistrement des énergies et gradients orbitaux

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.

Une entrée supplémentaire pour ce programme est lue à partir du fichier input.dat. Le plus couramment utilisé
les mots clés sont :

WFNPlus = un magnifique
Les types de fonction d'onde pris en charge sont CASSCF et RASSCF.

CONVERGENCE = entier
Convergence souhaitée sur le gradient orbital. La convergence est atteinte lorsque le RMS
de l'erreur dans le gradient orbital est inférieure à 10**(-n). La valeur par défaut est 4 pour
calculs énergétiques et 7 pour les gradients.

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.

ENERGIE_CONVERGENCE = entier
Convergence souhaitée sur l'énergie totale MCSCF. La valeur par défaut est 7.

CONGELÉ_DOCC = tableau_entiers
Le nombre d'orbitales les plus énergétiques doublement occupées dans chaque irréductible
représentation qui sera littéralement figée (non mise à jour dans le MCSCF). Les
L'ordre du coton des 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
qui sera littéralement gelé (non mis à jour dans le MCSCF). La valeur par défaut est le zéro
vecteur.

RESTRICTED_DOCC = tableau_entiers
Le nombre d'orbitales les plus énergétiques doublement occupées dans chaque irréductible
représentation qui sera optimisée mais gardée doublement occupée dans le MCSCF.
Ces orbitales viennent après le CONGELÉ_DOCC orbitales. La valeur par défaut est le zéro
vecteur.

RESTRICTED_UOCC = tableau_entiers
Le nombre d'orbitales inoccupées les plus énergétiques dans chaque représentation irréductible
qui sera optimisé mais maintenu inoccupé dans le MCSCF. Ces orbitales viennent
avant la CONGELÉ_UOCC orbitales. La valeur par défaut est le vecteur zéro.

NCASITER = entier
Nombre maximum d'itérations pour optimiser les orbitales. Cette option devrait be
spécifié in le DEFAULT of contribution, car il doit être visible pour le
programme de contrôle PSI. La valeur par défaut est 1.

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).

Les mots-clés les moins utilisés sont :

DIIS_FREQ = entier
La procédure d'extrapolation DIIS sera tentée chaque n itérations. Le défaut
est 1.

DIIS_MAX_VECS = entier
Le nombre maximal de vecteurs de sous-espace pour la procédure DIIS. Après ce numéro
des vecteurs sont atteints, les vecteurs plus anciens seront supprimés du sous-espace comme
nécessaire. La valeur par défaut est 8.

DIIS_MIN_VECS = entier
Le nombre minimum de vecteurs de sous-espace avant une interpolation DIIS peut être
effectué. La valeur par défaut est 2. Cela n'a pas de sens d'avoir des valeurs inférieures à 2.

DIIS_START = entier
Le numéro d'itération lors de l'inversion directe du sous-espace itératif de Pulay (DIIS)
la procédure d'accélération de la convergence doit être activée. La valeur par défaut est 3.
Avant cette itération, les vecteurs ne sont pas ajoutés au sous-espace DIIS. La première
L'étape DIIS ne sera pas prise avant DIIS_MIN_VECS les vecteurs sont dans le sous-espace DIIS.

LAG_FILE = entier
Fichier (numéro d'unité) de lecture de la matrice lagrangienne. La valeur par défaut est
actuellement 75.

OPDM_FILE = entier
Fichier (numéro d'unité) pour la lecture de la matrice de densité à une particule. La valeur par défaut
a actuellement 73 ans.

ÉCHELLE_GRAD = booléen
Indique s'il faut mettre à l'échelle le gradient orbital par l'orbitale approximative (diagonale)
Toile de jute. La valeur par défaut est VRAI.

SCALE_STEP = réal
Facteur d'échelle pour le pas de rotation orbitale. La valeur par défaut est 1.0.

TPDM_FILE = entier
Fichier (numéro d'unité) pour la lecture de la matrice de densité à deux particules. La valeur par défaut
a actuellement 74 ans.

LEVEL_SHIFT = booléen
Indique s'il faut autoriser le décalage de niveau de la matrice de jute. Cela peut être utilisé
pour s'assurer que la toile de jute est définie positive pour les itérations de début. Si
la toile de jute n'est pas définie positive près de la convergence, alors cela peut être un signe de
instabilités numériques dans le MCSCF. La valeur par défaut est VRAI.

SHIFT = réal
Indique de combien les éléments diagonaux de la toile de jute doivent être décalés pour le niveau
déplacement. La valeur par défaut est 0.01.

DETERM_MIN = réal
Minimum autorisé valorisé pour le déterminant de la matrice hessienne en cas de changement de niveau
est sur. Ceci n'est pas utilisé si LEVEL_SHIFT=FALSE. La valeur par défaut est 0.00001.

ÉTAPE_MAX = réal
Il s'agit de la rotation orbitale unique maximale autorisée. La valeur par défaut est 0.30.

FORCE_ÉTAPE = booléen
Cela permet à l'utilisateur d'outrepasser le pas calculé et de forcer un pas dans un
direction particulière. Cela peut être utile si vous essayez de forcer le calcul
d'un point de selle. Le défaut est faux.

FORCE_PAIR = entier
C'est l'indice de la paire indépendante qui doit être tournée si
FORCE_STEP=TRUE. Ceci est ignoré sinon. Il n'y a pas de défaut.

FORCE_VALUE = réal
C'est la valeur de rotation orbitale pour faire tourner les orbitales spécifiées par
FORCE_PAIR. Ceci n'est utilisé que si FORCE_STEP est défini sur TRUE. La valeur par défaut est
0.0.

CHECK_HESSIEN = booléen
Calculer les valeurs propres de la toile de jute orbitale. C'est bon pour vérifier
solutions de point de selle. Doit être utilisé sans changement de niveau. La valeur par défaut est
FAUX.

EIGEN_VECTEURS = booléen
Voulez-vous également calculer les vecteurs propres de la toile de jute orbitale. C'est
ignoré si CHECK_HESSIAN=FALSE. La valeur par défaut est FAUX.

8 mai, 1998 détcas(1)

Utilisez detcas en ligne en utilisant les services onworks.net