G03 Sample Output
G03 Sample Output
The first part of the Gaussian output file states in considerable detail the contents of the license agreement.
This should be taken seriously. Gaussian 03 is no public domain software!!
Entering Gaussian System, Link 0=/scr1/g03/g03
Initial command:
/scr1/g03/l1.exe /scr1/zipse/Gau-26301.inp -scrdir=/scr1/zipse/
Entering Link 1 = /scr1/g03/l1.exe PID= 26302.
Gaussian, Inc.
Carnegie Office Park, Building 6, Pittsburgh, PA 15106 USA
Gaussian, Inc.
Carnegie Office Park, Building 6, Pittsburgh, PA 15106 USA
---------------------------------------------------------------
Warning -- This program may not be used in any manner that
competes with the business of Gaussian, Inc. or will provide
assistance to any competitor of Gaussian, Inc. The licensee
of this program is prohibited from giving any competitor of
Gaussian, Inc. access to this program. By using this program,
the user acknowledges that Gaussian, Inc. is engaged in the
business of creating and licensing software in the field of
computational chemistry and represents and warrants to the
licensee that it is not a competitor of Gaussian, Inc. and that
it will not use this program in any manner prohibited above.
---------------------------------------------------------------
Actual program output specific to a certain calculation starts with a statement of the program version,
Gaussian revision (here B.03), system software (here LINUX), and the current date. Subsequently the
keywords used in the input file are repeated together with other general settings such as the amount of main
memory needed for the calculations (here 48MB), and the location of a binary checkpoint file for storage of
important results (here /scratch/test1.chk). The "scf=tight" keyword used here specifies tight criteria for the
energy calculation and the quantum mechanical method used is "HF/6-31G(d)".
*********************************************
Gaussian 03: x86-Linux-G03RevB.03 4-May-2003
18-Oct-2004
*********************************************
%mem=6000000
%chk=/scratch/test1.chk
------------------------
#P HF/6-31G(d) scf=tight
------------------------
The keywords are transformed by Gaussian into a sequence of subroutine calls termed "links". The links are
given together with the corresponding options set for each link in a proprietary format. Provided that the "#P"
option is used in the input file, Gaussian prints out elapsed CPU times after leaving a link.
1/38=1/1;
2/17=6,18=5,40=1/2;
3/5=1,6=6,7=1,11=9,16=1,25=1,30=1/1,2,3;
4//1;
5/5=2,32=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
99/5=1,9=1/99;
Leave Link 1 at Mon Oct 18 10:33:45 2004, MaxMem= 6000000 cpu: 0.4
In link101 the program reads in or retrieves from the checkpoint file the structure of the system together with
other parameters and prints the structure (in a slightly modified format) together with overall charge and spin
multiplicity and the comments supplied in the input file. It is good practice to include the name of the input
file in the comments of the job. The system chosen here is formaldehyde in its electronic ground state.
(Enter /scr1/g03/l101.exe)
---------------------------------
test1 HF/6-31G(d) sp formaldehyde
---------------------------------
Symbolic Z-matrix:
Charge = 0 Multiplicity = 1
C1
O2 1 r2
H3 1 r3 2 a3
H4 1 r4 2 a4 3 d4 0
Variables:
r2 1.2
r3 1.
r4 1.
a3 120.
a4 120.
d4 180.
Atom 1 2 3 4
IAtWgt= 12 16 1 1
AtmWgt= 12.0000000 15.9949146 1.0078250 1.0078250
IAtSpn= 0 0 1 1
AtZEff= 0.0000000 0.0000000 0.0000000 0.0000000
AtQMom= 0.0000000 0.0000000 0.0000000 0.0000000
AtGFac= 0.0000000 0.0000000 2.7928460 2.7928460
Leave Link 101 at Mon Oct 18 10:33:46 2004, MaxMem= 6000000 cpu: 0.3
Link202 determines, among others, the symmetry of the system, decides on the symmetry properties that will
be used in the actual quantum mechanical calculations and rotates the molecule such that the center of mass
is located in the origin of the cartesian coordinate system, the principal axis (so it exists) points along the z-
axis, and the principal plane of symmetry (so it exists) is located in the yz-plane. The resulting orientation is
printed as "Standard orientation", which serves as the reference description for all information regarding the
wavefunction and first and second derivatives of the energy with respect to structural parameters.
(Enter /scr1/g03/l202.exe)
Input orientation:
---------------------------------------------------------------------
Center Atomic Atomic Coordinates (Angstroms)
Number Number Type X Y Z
---------------------------------------------------------------------
1 6 0 0.000000 0.000000 0.000000
2 8 0 0.000000 0.000000 1.200000
3 1 0 0.866025 0.000000 -0.500000
4 1 0 -0.866025 0.000000 -0.500000
---------------------------------------------------------------------
Distance matrix (angstroms):
1 2 3 4
1 C 0.000000
2 O 1.200000 0.000000
3 H 1.000000 1.907878 0.000000
4 H 1.000000 1.907878 1.732051 0.000000
Stoichiometry CH2O
Framework group C2V[C2(CO),SGV(H2)]
Deg. of freedom 3
Full point group C2V NOp 4
Largest Abelian subgroup C2V NOp 4
Largest concise Abelian subgroup C2 NOp 2
Standard orientation:
---------------------------------------------------------------------
Center Atomic Atomic Coordinates (Angstroms)
Number Number Type X Y Z
---------------------------------------------------------------------
1 6 0 0.000000 0.000000 -0.537500
2 8 0 0.000000 0.000000 0.662500
3 1 0 0.000000 0.866025 -1.037500
4 1 0 0.000000 -0.866025 -1.037500
---------------------------------------------------------------------
Rotational constants (GHZ): 334.3034006 40.3785408 36.0270459
Leave Link 202 at Mon Oct 18 10:33:47 2004, MaxMem= 6000000 cpu: 0.3
link301 loads all components necessary for the actual quantum mechanical part of the calculation.
(Enter /scr1/g03/l301.exe)
Standard basis: 6-31G(d) (6D, 7F)
There are 18 symmetry adapted basis functions of A1 symmetry.
There are 2 symmetry adapted basis functions of A2 symmetry.
There are 6 symmetry adapted basis functions of B1 symmetry.
There are 8 symmetry adapted basis functions of B2 symmetry.
Integral buffers will be 262144 words long.
Raffenetti 1 integral format.
Two-electron integral symmetry is turned on.
34 basis functions, 64 primitive gaussians, 34 cartesian basis functions
8 alpha electrons 8 beta electrons
nuclear repulsion energy 32.2605629825 Hartrees.
IExCor= 0 DFT=F Ex=HF Corr=None ExCW=0 ScaHFX= 1.000000
ScaDFX= 1.000000 1.000000 1.000000 1.000000
IRadAn= 0 IRanWt= -1 IRanGd= 0 ICorTp=0
NAtoms= 4 NActive= 4 NUniq= 3 SFac= 2.05D+00 NAtFMM= 60 Big=F
Leave Link 301 at Mon Oct 18 10:33:49 2004, MaxMem= 6000000 cpu: 0.1
link303 calculates a number of integrals necessary for the subsequent SCF (energy) calculation.
(Enter /scr1/g03/l302.exe)
NPDir=0 NMtPBC= 1 NCelOv= 1 NCel= 1 NClECP= 1 NCelD= 1
NCelK= 1 NCelE2= 1 NClLst= 1 CellRange= 0.0.
One-electron integrals computed using PRISM.
One-electron integral symmetry used in STVInt
NBasis= 34 RedAO= T NBF= 18 2 6 8
NBsUse= 34 1.00D-06 NBFU= 18 2 6 8
Leave Link 302 at Mon Oct 18 10:33:50 2004, MaxMem= 6000000 cpu: 0.8
(Enter /scr1/g03/l303.exe)
DipDrv: MaxL=1.
Leave Link 303 at Mon Oct 18 10:33:51 2004, MaxMem= 6000000 cpu: 0.0
Before the actual energy calculation is performed, a guess for the wavefunction is obtained using either the
Hueckel, the INDO, or the Harris functional method. Alternatively, a guess can also be read from the
checkpoint or the input file.
(Enter /scr1/g03/l401.exe)
Harris functional with IExCor= 205 diagonalized for initial guess.
ExpMin= 1.61D-01 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=1 IRadAn= 1 AccDes= 1.00D-06
HarFok: IExCor= 205 AccDes= 1.00D-06 IRadAn= 1 IDoV=1
ScaDFX= 1.000000 1.000000 1.000000 1.000000
Harris En= -113.945080996335
Initial guess orbital symmetries:
Occupied (A1) (A1) (A1) (A1) (B2) (A1) (B1) (B2)
Virtual (B1) (A1) (B2) (A1) (B1) (A1) (B2) (A1) (A1) (B2)
(B1) (A1) (B2) (A1) (A2) (B1) (A1) (A2) (B2) (A1)
(A1) (B1) (B2) (A1) (A1) (A1)
The electronic state of the initial guess is 1-A1.
Leave Link 401 at Mon Oct 18 10:33:52 2004, MaxMem= 6000000 cpu: 0.3
Calculation of the HF/6-31G(d) energy of the system is done in link 502. Some parameters such as the
currently selected convergence criteria are listed first. The final SCF energy given as E(RHF) =
-113.852967419 is the energy of the system with respect to its nuclei and electrons at infinite separation. The
energy is given in atomic units (Hartree).
(Enter /scr1/g03/l502.exe)
Closed shell SCF:
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on energy=1.00D-06.
No special actions if energy rises.
Using DIIS extrapolation, IDIIS= 1040.
Integral symmetry usage will be decided dynamically.
Keep R1 integrals in memory in canonical form, NReq= 649414.
IEnd= 21265 IEndB= 21265 NGot= 6000000 MDV= 5819242
LenX= 5819242
Symmetry not used in FoFDir.
MinBra= 0 MaxBra= 2 Meth= 1.
IRaf= 0 NMat= 1 IRICut= 1 DoRegI=T DoRafI=F ISym2E= 0 JSym2E=0.
Selected information on the optimized wavefunction is printed along with a Mulliken population analysis in
link601.
(Enter /scr1/g03/l601.exe)
Copying SCF densities to generalized density rwf, ISCF=0 IROHF=0.
**********************************************************************
**********************************************************************
Orbital symmetries:
Occupied (A1) (A1) (A1) (A1) (B2) (A1) (B1) (B2)
Virtual (B1) (A1) (B2) (A1) (B1) (A1) (B2) (A1) (A1) (B2)
(A1) (B1) (B2) (A1) (A2) (B1) (A1) (A2) (B2) (A1)
(A1) (B1) (B2) (A1) (A1) (A1)
The electronic state is 1-A1.
Alpha occ. eigenvalues -- -20.57603 -11.31546 -1.40250 -0.89790 -0.72365
Alpha occ. eigenvalues -- -0.64488 -0.52755 -0.44339
Alpha virt. eigenvalues -- 0.14325 0.27254 0.34861 0.37647 0.73913
Alpha virt. eigenvalues -- 0.80034 0.82453 0.97649 1.10532 1.18714
Alpha virt. eigenvalues -- 1.22933 1.27415 1.33708 1.62596 1.76972
Alpha virt. eigenvalues -- 1.80651 2.00802 2.24294 2.28103 2.51332
Alpha virt. eigenvalues -- 2.70989 2.89384 3.04381 3.30687 4.11684
Alpha virt. eigenvalues -- 4.52410
Condensed to atoms (all electrons):
1 2 3 4
1 C 4.535451 0.545272 0.388868 0.388868
2 O 0.545272 8.013683 -0.057728 -0.057728
3 H 0.388868 -0.057728 0.597808 -0.079928
4 H 0.388868 -0.057728 -0.079928 0.597808
Mulliken atomic charges:
1
1 C 0.141540
2 O -0.443500
3 H 0.150980
4 H 0.150980
Sum of Mulliken charges= 0.00000
Atomic charges with hydrogens summed into heavy atoms:
1
1 C 0.443500
2 O -0.443500
3 H 0.000000
4 H 0.000000
Sum of Mulliken charges= 0.00000
Electronic spatial extent (au): = 57.2540
Charge= 0.0000 electrons
Dipole moment (field-independent basis, Debye):
X= 0.0000 Y= 0.0000 Z= -2.8199 Tot= 2.8199
Quadrupole moment (field-independent basis, Debye-Ang):
XX= -11.2700 YY= -11.2435 ZZ= -11.7586
XY= 0.0000 XZ= 0.0000 YZ= 0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
XX= 0.1540 YY= 0.1805 ZZ= -0.3346
XY= 0.0000 XZ= 0.0000 YZ= 0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
XXX= 0.0000 YYY= 0.0000 ZZZ= -0.2636 XYY= 0.0000
XXY= 0.0000 XXZ= 0.3457 XZZ= 0.0000 YZZ= 0.0000
YYZ= -0.4711 XYZ= 0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX= -9.0469 YYYY= -15.1879 ZZZZ= -40.7625 XXXY= 0.0000
XXXZ= 0.0000 YYYX= 0.0000 YYYZ= 0.0000 ZZZX= 0.0000
ZZZY= 0.0000 XXYY= -4.2304 XXZZ= -8.2015 YYZZ= -9.1659
XXYZ= 0.0000 YYXZ= 0.0000 ZZXY= 0.0000
N-N= 3.226056298246D+01 E-N=-3.330532990315D+02 KE= 1.138889777761D+02
Symmetry A1 KE= 1.035048839320D+02
Symmetry A2 KE= 1.383951977151D-34
Symmetry B1 KE= 3.440292093306D+00
Symmetry B2 KE= 6.943801750813D+00
No NMR shielding tensors so no spin-rotation constants.
Leave Link 601 at Mon Oct 18 10:34:00 2004, MaxMem= 6000000 cpu: 1.1
At the very end of each Gaussian calculation, an archive entry in a very compact format is printed to
summarize the results. This archive entry is frequently used as supplemental material in publications of
theoretical results.
(Enter /scr1/g03/l9999.exe)
1\1\GINC-TERMINUS\SP\RHF\6-31G(d)\C1H2O1\ZIPSE\18-Oct-2004\0\\#P HF/6-
31G(D) SCF=TIGHT\\test1 HF/6-31G(d) sp formaldehyde\\0,1\C\O,1,1.2\H,1
,1.,2,120.\H,1,1.,2,120.,3,180.,0\\Version=x86-Linux-G03RevB.03\State=
1-A1\HF=-113.8529674\RMSD=2.649e-09\Dipole=0.,0.,-1.1094497\PG=C02V [C
2(C1O1),SGV(H2)]\\@
From a database of citations, Gaussian prints one entry together with some timing information.
MOST BUREAUCRACIES FUNCTION LIKE A SEPTIC TANK.
THE LARGE CHUNKS RISE TO THE TOP.
Job cpu time: 0 days 0 hours 0 minutes 9.3 seconds.
File lengths (MBytes): RWF= 11 Int= 0 D2E= 0 Chk= 7 Scr= 1
Normal termination of Gaussian 03 at Mon Oct 18 10:34:16 2004.