Qumica computacional

PRÁCTICA 1
Objetivo: Determinar las energía de la función de onda con un calculo tipo RHF Y MSCF , a diferentes radios para encontrar el mínimo. Comparar la diferencia entre uno y otro para ver como el calculo mejoró.
Ejemplo del input
$CONTRL MAXIT=200 UNITS=ANGS SCFTYP=RHF $END
$CONTRL ICHARG=0 MULT=1 RUNTYP=ENERGY $END
$GUESS GUESS=HCORE $END
$SCF $END
$DATA
H2 singuleteEquipo MARVISA CALCULO RHF distacia = 0.5 A.
C1
H1 1.0 0.0 0.0 0.0
S 3
1 18.7311370 0.0334946000
2 2.8253937 0.2347269500
3 0.6401217 0.8137573300
S 1
1 0.1612778 1.0000000000
S 1
1 0.0360000 1.0000000000
P 1
1 1.1000000 1.0000000000
H21.0 0.0 0.0 0.5
S 3
1 18.7311370 0.0334946000
2 2.8253937 0.2347269500
3 0.6401217 0.8137573300
S 1
1 0.1612778 1.0000000000
S 1
1 0.0360000 1.0000000000
P 1
1 1.1000000 1.0000000000
$END

r | E (Hartrees)RHF | E (Hartrees)MSCF |
6 | -0.75946573| -0.99760279 |
5 | -0.7705299 | -0.99760441 |
4 | -0.79035735 | -0.99764455 |
3 | -0.83133089 | -0.99859909 |
2 | -0.92280993 | |
1 | -1.09990147 | -1.13255598 |
0.75 | -1.13119984 | -1.15744525 |
0.5 | -1.06149218 | -1.08286367 |
0.25 | -0.37521748 | -0.39178527 |
















PRACTICA 2 Molécula H2O
Objetivo: Calcular la energía del H20con RHF, RHF optimizado, usar simetría CN2V, y hacerlo con coordenadas internas.
Ej. Input con simetría
$CONTRL MAXIT=200 UNITS=ANGS SCFTYP=RHF $END
$CONTRL ICHARG=0 MULT=1 RUNTYP=OPTIMIZE $END
$GUESS GUESS=HCORE $END
$SCF $END
$DATA
H2O MARVISA CALCULO RHF OPTIMIZADO GEOMETRIA C2V
CNV 2
O 8.0 0.0 0.0 0.0
S 3
1 322.03700000.0592394
2 48.4308000 0.3515000
3 10.4206000 0.7076580
L 2
1 7.4029400 -0.4044530 0.2445860
2 1.5762000 1.2215600 0.8539550
L 1
1 0.3736840 1.0000000 1.0000000
L 1
1 0.0845000 1.00000001.0000000
H1 1.0 0.0 0.76 0.585
S 3
1 18.7311370 0.0334946000
2 2.8253937 0.2347269500
3 0.6401217 0.8137573300
S 1
1 0.1612778 1.0000000000
S 1
1 0.0360000 1.0000000000
P 1
1 1.1000000 1.0000000000
$END


Ej input concoordenadas internas
$CONTRL MAXIT=200 UNITS=ANGS SCFTYP=RHF COORD=ZMT $END
$CONTRL ICHARG=0 MULT=1 RUNTYP=ENERGY $END
$BASIS GBASIS=STO NGAUSS=2 $END
$GUESS GUESS=HCORE $END
$SCF $END
$DATA
H2O MARVISA CALCULO RHF ENERGIA COORDENADAS ZMT
CNV 2

O
H 1 rOH
H 1 rOH 2 aHOH

rOH=0.959
aHOH=104.8

$END
También se obtuvieron los orbitales molecularesdel H2O

PRACTICA SF4
Objetivo: usar coordenadas internas con bases STOG 3-G y 4-31G
PRACTICA BENCENO
Objetivo: usar las comodidades de simetría con RHF y una base STO
PRACTICA CF4, CH4, CF3H
Objetivo analizar las poblaciones de Mulliken usando optimización.
CF4
TOTAL MULLIKEN AND LOWDIN ATOMIC POPULATIONS |
ATOM MULL.POP. CHARGE | |
1C1 4.802899 1.197101 | | |
2 F2 9.299276 -0.299276 | | |
3 F3 9.299282 -0.299282 | | |
4 F4 9.299277 -0.299277 | | |
5 F5 9.299266 -0.299266 CF3H | | |
TOTAL MULLIKEN AND LOWDIN ATOMIC POPULATIONS |
ATOM MULL.POP. CHARGE...