Cost_package/Manuals/noscf_emploi

98 lines
2.8 KiB
Plaintext
Executable File

noscf : obtention d'orbitales optimisees par iterations sur la
matrice densite.
genere les donnees et appelle les programmes :
- motra
- molcost ou molcsd
- casdet
- casdi
- locnats
puis itere jusqu'a convergence
Utilisation :
cp "Orbitales" INPORB (avec schmud : cp LOCORB INPORB)
export CurrDir WorkDir Project
export PATH=$PATH:"lieu_de_cost"/bin
noscf donnees (ET NON i"noscf < donnees !!!)
orbitales de sortie : LOCORB2
Donnees :
&nofil PREFIX=
/
&no
data=' ' 'AUTO' lit quelques donnees d'un fichier de schmudorb
les donnees de schmudorb pilotent ainsi une grande partie de
la suite du calcul (defaut ' ')
FERMI=0 par sym, nb d'orb occupees
(lu de schmudorb si data='AUTO')
MOLCAS=5 (ou 4)
frozen=0 donne par sym
(lu de schmudorb si data='AUTO')
delete=0 donne par sym
(lu de schmudorb si data='AUTO')
maxit=50 nb max de cycles motra...locnats
r=' > ' (les fichiers de sortie sont dans la
directory de travail (sinon r=' # ')
nmat=1 nb de fichiers de mat densite
netat= nb d'etats par fichier (nmat donnees)
metat= numeros des etats demandes (Somme(netat)) donnees
coef=1.d0 coefs des etats (Somme(netat)) donnees
stop=' ' ('MOLCOST' pour arret apres molcsd, donne la numerotation
des orbitales)
seuil=1.d-4 seuil de convergence du processus (sur les elements des
blocs exterieurs de la matrice densit)
/
&casdi (nmat fois)
syspin='0',
nessai=0,
iprec=6
noac=0, (lu de schmudorb si data='AUTO')
numac= (lu de schmudorb si data='AUTO')
NELAC=0, (Attention, non lu ! A donner imperativement)
gener='CAS+S',
is0=1
ms2=0
/
les variables en majuscules cont celles obligatoires
les valeurs donnees sont celles par defaut
Donnees supplementaires:
noscf genere, a partir de &nofil, &no et &casdi les donnees pour
motra, molcost, casdet, casdi et locnats.
Il peut arriver qu'on veuille donner a un de ces programmes un donnee
supplementaire. Par exemple, si on veut donner la donnee twvec a &dav
pour le 2eme calcul casdi (alors nmat >=2).
On ajoute a la fin du fichier :
DAV2
twvec=0.01
FIN
explication : DAV pour &dav, 2 pour 2eme casdi
la (les) donnee(s)
FIN (de l'ajout a &dav)
la ligne est ajoutee avant la fin de namelist "/"
Rem : si on veut ajouter un vecteur d'essai en donnee
DAV1
nessai=1
/
&ess file='ESSOUT0'
FIN
le fichier de donnees devient apres rajout des 3 lignes :
&dav
...
...
nessai=1 \
/ > 3 lignes supplementaires
&ess file='ESSOUT0' /
/ (deja present, donc pas a rajouter)