Startseite |
If you trudge somewhere, you walk there with slow, heavy steps.
Input (nur die für TRUDGE relevanten Zeilen sind angegeben):
$CONTRL
SCFTYP=RHF MULT=1
RUNTYP=TRUDGE
COORD=HINT CITYP=GUGA $END
$TRUDGE
OPTMIZ=GEOMETRY NPAR=1 IEX(1)=11 $END
H2-Molekuel, CI
DNH 2
H
1.0 LC
0.500 0.0
0.0 +
O
$END
NPAR=1: number of parameters to be optimized. Das ist klar, nur der Bindungsabstand wird variiert.
IEX= defines the parameters to be
optimized. Pointers to the HINT internal coordinates which will be optimized.
The pointers to the internal coordinates are defined as: (the number of atom in
the input list)*10 + (the number of internal coordinate for that atom). For
each atom, the HINT internal coordinates are numbered as 1, 2, and 3 for BOND,
ALPHA, and BETA, respectively.
IEX(1)=11: Die erste 1 bedeutet Atom 1 (es gibt nur
dieses in $DATA). Die zweite 1 steht für Bindungslänge.
P= Defines the initial values of the parameters to be optimized. If omitted, the $DATA values are used. Da P in $TRUDGE nicht auftaucht, ist der Startwert für die Bindungslänge 0.5 (wegen der Symmetrie bedeutet das Bindungsabstand 1.0 Å).
Im Outputfile wird das wie folgt wiedergegeben:
NON-GRADIENT ENERGY MINIMIZATION
-------------------------------------------------------------
$TRURST KSTART=-1 JSTART=
0 TOLF=
0.001000 TOLR= 0.050000
FNOISE=
0.000500 FNOT=
0.000000000
$END
-----------------------------------------------------------
INTERNAL COORDINATES (ANGSTROMS/DEGREE)
I BOND
ALPHA BETA
SIGN ICONX IATCON INATOM
1 0.5000000 0.00000
0.00000 1.0
1501502503 2
-----------------------------------------------------------
----- TRUDGE RESTART DATA AT NSTEP
0 --------
$TRUDGE OPTIMIZE=GEOMETRY NPAR=
1
IEX(1)=
11,
P(1)= 0.500000,
TRUDGE ENERGY VALUE AT NSTEP= 0
IS -0.7467206299
$TRUDGE OPTMIZ=GEOMETRY NPAR=2 IEX(1)=21,22 P(1)=0.96 $END
$DATA
Water CI-SD Geometry Optimization
CNV 2
O 8.0 LC 0.0 0.0 0.0 - O K
H 1.0 PCC 0.94 53.0 0.0 + O K I
NPAR= 2: Bindungslänge und Bindungswinkel sollen optimiert werden.
IEX(1)=21,22: 21: Bindungslänge des H-Atoms. 22: Bindungswinkel, den das H-Atom bildet. Ob die 1 in Klammern das Bezugsatom darstellt ist noch unklar.
P(1)=0.96: Startwert für Bindungslänge.
N2O4-Molekuel
CNH 2
N 7.0 LC 0.90
0.0 0.0 + O
O 8.0 PCC 1.20 112.0
0.0 + 1
O 8.0 PCC 1.20 -112.0
0.0 + 1
$TRUDGE OPTMIZ=GEOMETRY NPAR=1 IEX(1)=11 $END
Damit wird nur der N-N-Abstand variiert.
$TRUDGE OPTMIZ=GEOMETRY NPAR=5 IEX(1)=11,21,22,31,32 $END
Damit werden alle Parameter optimiert.
Seitenanfang | Startseite |