Um eine Bindungslänge, einen Winkel oder ein Dieder zu fixieren, ist es wichtig, die Direktive Opt=ModRedundant zu setzen.

%nprocshared=20
%mem=40GB
%chk=/global/archive/phchm/ir-goettingen/s0opt2/ir_opt.chk
#P b3lyp/def2SVP opt=ModRedundant scf=(tight,MaxCycle=500) pop=full iop(6/7=3) gfinput gfoldprint gfprint EmpiricalDispersion=GD3 SCRF=(PCM,Solvent=Acetonitrile)

Ir complex with B3LYP

1 1
C -2.38008 0.32587 0.65976
S -4.86842 0.34670 -0.74159
S -5.07214 1.14093 0.59786
C -3.51219 -0.45064 -1.51185
...
Ir 0.59352 0.00001 -0.00005

D 1 2 3 4 0.00 B
D 1 2 3 4 F

Das D 1 2 3 4 0.00 B gibt ein Dieder an (weitere Möglichkeiten: X für kartesische Koordinate, B für Bindungslänge, A für Winkel oder L für lineare Biegung), die vier beteiligten Atome (angegeben durch ihre Nummer in der xyz-Liste), den Wert des Dieders (hier 0.00) und B für "build", d.h. diese redundante Definition erzeugen.

Die folgende Zeile, D 1 2 3 4 F, friert das Dieder auf den zuvor angegebenen Wert ein. Hier wird das Dieder auf 0 Grad fixiert.

Hier ist auch ein SLURM-Submission-Script:

#!/bin/bash
#SBATCH -J opt
#SBATCH -N 1
#SBATCH --ntasks-per-node=20
#SBATCH -p compute
#SBATCH --mem=42G
#SBATCH --constraint beegfs

module use /usr/license/modulefiles
module load gaussian
export g16root=/usr/license/gaussian
source $g16root/g16/bsd/g16.profile
export TMPDIR=/public/`hostname -s`/scratch/tmp/$LOGNAME.$$
mkdir $TMPDIR

g16 < ir_opt.inp > ir_opt.log

rm -rf $TMPDIR