Startseite    

TRUDGE, non-gradient energy minimization

If you trudge somewhere, you walk there with slow, heavy steps.

H2

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

H2O

 $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

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