Geometrieoptimierung unter Beibehaltung einiger Geometrieparameter. Z. B. Berechnung der Potentialkurve (Reaktionsprofil) der Dissoziation von N2O4 in 2 NO2.
!
! N2O4-Molekuel
!
$CONTRL SCFTYP=RHF MULT=1 RUNTYP=OPTIMIZE NZVAR=12 COORD=ZMT $END
$SYSTEM TIMLIM=1000 MEMORY=5000000 $END
$ZMAT IZMAT(1)=1,2,1 1,3,2 1,4,1 1,5,1 1,6,2 2,3,2,1 2,4,1,2
2,5,1,2 2,6,2,1 3,4,1,2,3 3,5,1,2,3 3,6,2,1,4 $END
$STATPT IFREEZ(1)=1 $END
$BASIS GBASIS=STO NGAUSS=6 $END
$GUESS GUESS=HUCKEL $END
$DATA
N2O4-Molekuel
DNH 2
N1
N2 1 1.800
O3 2 1.200 1 112.0
O4 1 1.200 2 112.0 3 180.0
O5 1 1.200 2 112.0 3 0.0
O6 2 1.200 1 112.0 4 0.0
$END
Die N-N-Bindung wird eingefroren: IFREEZ(1)=1: erster Geometrieparameter in IZMAT(1). Alle übrigen Parameter werden optimiert.
Output
--------------------
INTERNAL COORDINATES
--------------------
- - ATOMS - - COORDINATE COORDINATE
NO. TYPE I J K L M N (BOHR,RAD)
(ANG,DEG)
----------------------------------------------------------------
1 STRETCH 2 1
3.4015068 1.8000000
2 STRETCH 3 2
4.8455146 2.5641361
3 STRETCH 4 1
4.8455146 2.5641361
4 STRETCH 5 1
4.8455146 2.5641361
5 STRETCH 6 2
4.8455146 2.5641361
6 BEND 3 2 1
0.4643552 26.6055941
7 BEND 4 1 2
0.4643552 26.6055941
8 BEND 5 1 2
0.4643552 26.6055941
9 BEND 6 2 1
0.4643552 26.6055941
10 TORSION 4 1 2 3 3.1415927
180.0000000
11 TORSION 5 1 2 3 0.0000000
0.0000000
12 TORSION 6 2 1 4 0.0000000
0.0000000
Achtung! Die Nummerierung hat sich
im Output gegenüber dem Input geändert!
z
|
O3 | O6
|
N1
|
----------------|---------------> x
|
N2
|
O5 | O4