Gaussian Input Files


The following is a typical Gaussian input file, using a single point energy calculation on formaldehyde (HF/6-31G(d) level) as an example:
%mem=6000000
%chk=/scratch/test1.chk
#P HF/6-31G(d) scf=tight

test1 HF/6-31G(d) sp formaldehyde

0 1
C1
O2  1  r2
H3  1  r3  2  a3
H4  1  r4  2  a4  3  d4

r2=1.20
r3=1.0
r4=1.0
a3=120.
a4=120.
d4=180.
      

The first line of the input file:


%mem=6000000 
 
specifies how much memory is available for this calculation (in words). The default value (if nothing is specified) amounts to 6000000 words (corresponding to 48MB on a 32bit operating system). This information can also be specified with

%mem=48MB 
 
The second line of the input file:

%chk=/scratch/test1.chk

defines a binary "checkpoint" file to which important results are written. If no file name is supplied in the input file, the checkpoint files are generated during program execution, but deleted on successful termination of the program. The line directly following these two lines is the first one absolutely necessary for the program to execute at all and contains keywords defining the major steps of the calculation to be performed.

#P HF/6-31G(d) scf=tight

Here a HF/6-31G(d) single point calculation is done. The keywords must be preceded by the "#" sign. Inclusion of "P" provides additional output. The keyword line can actually be longer than one physical line and is concluded by one additional blank line. Following the blank line starts the comments section (in free format) providing space for any type of comment. It is good practice to include the name of the input file also here together with a short description of the job. The comments section is concluded by one additional blank line.

test1 HF/6-31G(d) sp formaldehyde
 
Before specifying the structure of the system, the overall charge of the system and the overall spin multiplicity must be given. The spin multiplicity "1" describes a singlet state, "2" a doublet state, "3" a triplet state . . . .

0 1

The structure of the system is described in the lines directly following the charge/spin multiplicity information. Beware of additional blank lines in this case! In the current example the structure of formaldehyde is described in an internal coordinate system on 11 lines. The geometry specification section is concluded by at least one blank line (more doesn`t hurt).
C1
O2  1  r2
H3  1  r3  2  a3
H4  1  r4  2  a4  3  d4

r2=1.20
r3=1.0
r4=1.0
a3=120.
a4=120.
d4=180.
      


A series of separate calculations can be combined into one input file by separating the various input sections with the Link1 comand. The following example demonstrates how the energy of formaldehyde can be calculated with two different theoretical methods (here HF/6-31G(d) and AM1):

%mem=6000000
%chk=/scratch/test1.chk
#P HF/6-31G(d) scf=tight

test1 HF/6-31G(d) sp formaldehyde

0 1
C1
O2  1  r2
H3  1  r3  2  a3
H4  1  r4  2  a4  3  d4

r2=1.20
r3=1.0
r4=1.0
a3=120.
a4=120.
d4=180.


--Link1--
%mem=6000000
%chk=/scratch/test1.chk
#P AM1 scf=tight geom=check

test1 AM1 sp formaldehyde

0 1


      


Please observe that the geometry of the system is only given explicitly for the first of these calculations. The second job step (AM1 energy calculation) reads the geometry from the checkpoint file /scratch/test1.chk using the geom=check command. This will only work correctly if there are TWO BLANK LINES separating the two parts of this compound job.


last changes: 16.10.2004, HZ
questions & comments to: zipse@cup.uni-muenchen.de