mirror of https://github.com/LCPQ/Cost_package.git
254 lines
8.0 KiB
Plaintext
Executable File
254 lines
8.0 KiB
Plaintext
Executable File
CALCULS CHOLESKY CALCULS CHOLESKY CALCULS CHOLESKY CALCULS CHOLESKY
|
|
16 dec 2010
|
|
|
|
|
|
I Calculs CASDI
|
|
|
|
I.1. Les integrales bi-electroniques sont obtenues dans motra.
|
|
=========================================================
|
|
La seule difference dans les fichiers de donnees est d'ajouter CHOL ou
|
|
HIGH ou LOW dans les donnees de seward.
|
|
|
|
- Le calcul motra est tres long, et cette methode n'est pas la meilleure.
|
|
|
|
|
|
I.2. Motra transforme les vecteurs de Cholesky, qui sont transformes en
|
|
integrales par un programme ulterieur (chol_X).
|
|
=========================================================
|
|
chol_X ne les calcule
|
|
pas toutes. Il commence par calculer les integrales d'echange, puis,
|
|
des seuils etant donnes, ne calcule que les integrales entre orbitales
|
|
en "interaction forte"
|
|
|
|
seward donnee supplementaire "CHOL" ou "HIGH" ou "LOW"
|
|
scf
|
|
(rasscf)
|
|
|
|
motra donnees supplementaires (2 lignes): CTONLY
|
|
(dans cet ordre !) KPQ
|
|
(Cholesky Transformation ONLY ; kpq : dans les
|
|
vecteurs cholesky V(k,pq) k sont les lignes
|
|
et pq les colonnes
|
|
PAS DE GEL EN DONNEE !
|
|
|
|
molcost donnees supplementaires : ychol=T
|
|
frozen, delete, fermi (apres gel) sont
|
|
donnes ICI
|
|
si on veut reordonner les orbitales en
|
|
Occ Act Virt (pour Nevpt2) donner ICI
|
|
act, ninact
|
|
|
|
chol_X donnees : &chol sli1=1.d-7,sli2=1.d-7 /
|
|
dans l'integrale (ij|kl):
|
|
sli1 : seuil pour considerer la distance i-j
|
|
sli2 : seuil pour considerer la distance ij-kl
|
|
|
|
casdet
|
|
casdi donnee obligatoire : itest(63)=1 dans &dav
|
|
(les integrales sont lues sur ijcl1,ijcl2,ijcl3)
|
|
|
|
|
|
II. Calculs en Orbitales Localisees
|
|
|
|
II.1. Calculs CASDI
|
|
|
|
seward donnee supplementaire "CHOL" ou "HIGH" ou "LOW"
|
|
scf
|
|
(rasscf)
|
|
|
|
dolo donner &smu ychol=T dans DOLOIN
|
|
(proj_scf)
|
|
schmudort
|
|
cp LOCORB INPORB
|
|
|
|
motra DOLO_MOTRAIN
|
|
DOLO_MOTRAIN a ete cree dans dolo
|
|
par rapport a motra, il a les 2 lignes
|
|
supplementaires :
|
|
CTONLY
|
|
KPQ
|
|
(Cholesky Transformation ONLY ; kpq : dans les
|
|
vecteurs cholesky V(k,pq) k sont les lignes
|
|
et pq les colonnes
|
|
pas de gel ni de del en donnee
|
|
|
|
molcost < DOLO_MOLCOSTIN
|
|
DOLO_MOLCOSTIN a ete cree dans dolo
|
|
avec la donnee supplementaire : ychol=T
|
|
frozen, delete, fermi (apres gel) sont
|
|
donnes ici
|
|
|
|
chol_X donnees : &chol sli1=1.d-7,sli2=1.d-7 /
|
|
dans l'integrale (ij|kl):
|
|
sli1 : seuil pour considerer la distance i-j
|
|
sli2 : seuil pour considerer la distance ij-kl
|
|
casdet
|
|
casdi donnee obligatoire : itest(63)=1
|
|
(les integrales sont lues sur ijcl1,ijcl2,ijcl3)
|
|
|
|
|
|
II.2. Calculs CASDILOC
|
|
|
|
Motra transforme les vecteurs de Cholesky. Puis molcost les reecrit
|
|
dans un fichier lu par ijkloc. ijkloc transforme les vecteurs de Cholesky
|
|
de la base atomique a la base moleculaire, puis calcule les integrales en
|
|
bas locale, conformement a ds seuils en donnee.
|
|
|
|
seward donnee supplementaire "CHOL" ou "HIGH" ou "LOW"
|
|
scf
|
|
(rasscf)
|
|
|
|
si localisation dolo :
|
|
dolo donner &smu ychol=T dans DOLOIN
|
|
(proj_scf)
|
|
schmudort
|
|
cp LOCORB INPORB
|
|
si localisation lewis :
|
|
lewis
|
|
cp LOCORB INPORB
|
|
donner ychol=T dans &lewis
|
|
|
|
|
|
|
|
motra DOLO_MOTRAIN (DOLO_MOTRAIN a ete genere dans dolo ou lewis)
|
|
DOLO_MOTRAIN a ete cree dans dolo (ou lewis)
|
|
par rapport a motra, il a les 2 lignes
|
|
supplementaires (a cause de la donnee ychol=T
|
|
dans dolo ou lewis) :
|
|
CTONLY
|
|
KPQ
|
|
(Cholesky Transformation ONLY ; kpq : dans les
|
|
vecteurs cholesky V(k,pq) k sont les lignes
|
|
et pq les colonnes
|
|
pas de gel ni de del en donnee
|
|
|
|
molcost < DOLO_MOLCOSTIN (DOLO_MOLCOSTIN a ete genere dans dolo ou lewis)
|
|
DOLO_MOLCOSTIN a ete cree dans dolo (ou lewis)
|
|
avec la donnee supplementaire : ychol=T
|
|
frozen, delete, fermi (apres gel) sont
|
|
donnes ici
|
|
|
|
|
|
ijkloc donnee supplementaire : ychol=T
|
|
sli1 : seuil, dans (ijkl) pour les
|
|
distributions i-j et k-l
|
|
sli2 : seuil d'interaction de ij avec kl
|
|
|
|
casdiloc
|
|
|
|
|
|
exemple de fichier casdiloc
|
|
|
|
#!/bin/bash
|
|
Project='acetone'
|
|
CurrDir=$PWD
|
|
WorkDir=$CurrDir/TMP
|
|
export Project WorkDir CurrDir
|
|
export PATH=$PATH:/home/daniel/cost2/bin
|
|
export MOLCAS=/home/daniel/7.7.dev
|
|
echo 'Start of job:' $Project
|
|
echo 'Current directory:' $CurrDir
|
|
echo 'Scratch directory:' $WorkDir
|
|
cd $WorkDir
|
|
#------------------------------------------------------------------------------#
|
|
# Start executing molcas job #
|
|
#------------------------------------------------------------------------------#
|
|
|
|
|
|
######################seward scf############
|
|
rm -fr $WorkDir
|
|
mkdir $WorkDir
|
|
cd $WorkDir
|
|
molcas $CurrDir/sewin
|
|
molcas $CurrDir/scfin
|
|
cp ../TMP0/acetone.ScfOrb .
|
|
############################################
|
|
|
|
##############interface#####################
|
|
molcost << EOF
|
|
&cost prefix='acetone.', ymo=f /
|
|
EOF
|
|
###########################################
|
|
|
|
######## LOCALISATION ###################################
|
|
dolo
|
|
proj_scf << EOF
|
|
&pscf prefix='acetone.' /
|
|
EOF
|
|
mv NONORLOC_scf NONORLOC
|
|
schmudort << EOF
|
|
&ort /
|
|
EOF
|
|
###########################################
|
|
|
|
############# motra ########################
|
|
cp LOCORB INPORB
|
|
molcas DOLO_MOTRAIN
|
|
###########################################
|
|
|
|
############# molcost ########################
|
|
molcost < DOLO_MOLCOSTIN
|
|
###########################################
|
|
|
|
############# ijkloc ########################
|
|
ijkloc << EOF
|
|
&ijkloc prefix='acetone.',ychol=t,nprint=1
|
|
symm='XY X',sli=0.d0 /
|
|
EOF
|
|
###########################################
|
|
|
|
############# casdiloc ##################
|
|
casdiloc << EOF
|
|
&locci prefix='acetone.',gener='CAS+SD'
|
|
/
|
|
EOF
|
|
###########################################
|
|
1,1 Haut
|
|
DOLOIN:
|
|
|
|
&smufil prefix='acetone.' /
|
|
&smu
|
|
nprint=0,ychol=t
|
|
orb='acetone.ScfOrb' /
|
|
/
|
|
&oao /
|
|
C* 1s(1) pr=1
|
|
O* 1s(1) pr=1
|
|
C* 1s(2) pr=2
|
|
O* 1s(2) pr=2
|
|
H* 1s(1) pr=2
|
|
C* 2p(1) pr=2
|
|
O* 2p(1) pr=2
|
|
fin
|
|
&orb symm='XY X' /
|
|
&smufil prefix='acetone.' /
|
|
&smu
|
|
nprint=0,ychol=t
|
|
orb='acetone.ScfOrb' /
|
|
/
|
|
&oao /
|
|
C* 1s(1) pr=1
|
|
O* 1s(1) pr=1
|
|
C* 1s(2) pr=2
|
|
O* 1s(2) pr=2
|
|
H* 1s(1) pr=2
|
|
C* 2p(1) pr=2
|
|
O* 2p(1) pr=2
|
|
fin
|
|
&orb symm='XY X' /
|
|
|
|
|
|
OLD_DATA:
|
|
|
|
BOND_ BASREF='O:1S(2) 2P(1),C:1S(2) 2P(1),H:1S(1)'
|
|
BOND_ CHAIN='C1=O1(LP2)'
|
|
BOND_ CHAIN='C1-C2'
|
|
BOND_ CHAIN='C2-H1'
|
|
BOND_ CHAIN='C2-H2'
|
|
CORE_ C* 1S(1) label='G'
|
|
CORE_ O* 1S(1) label='G'
|
|
|
|
|
|
|
|
|