Tutorial - Quantum Chemistry - Intro to Gaussian I

This is part one of a two part series of using Gaussian at SCS. Topics that are covered include: Generating an input file using Avogadro, Running a single-point Energy calculation, Calculation molecular orbitals, and Visualizing orbitals. Update: to submit Gaussian jobs on the cluster lop please use new script submit-g16.

 QUANTUM CHEMISTRY WITH GUASSIAN: A VERY BRIEF INTRODUCTION 

TARAS V. POGORELOV AND MIKE HALLOCK
SCHOOL OF CHEMICAL SCIENCES, UIUC
Created 2011. Updated/reviewed October 28, 2024.
This tutorial is designed to help getting started with Gaussian. Single point energy calculation is performed on a single water molecule. Molecular orbitals are calculated. Basic knowledge of Unix/Linux is assumed. Estimated time to complete this tutorial is 50 min.

To download and save this tutorial to your computer, print this page <Command>+<P> or <Ctrl>+<P>) and save the document as a PDF.  

SOFTWARE

You will need:
  • text editor
    • MacOS:     TextEdit, Vim
    • Windows:  Notepad, Vim
  • ssh client
    • MacOS:     built in Terminal, scp or Fugu for file transfers
    • Windows:   putty, WinSCP for file transfer
  • molecule builder

  1. PREPARE GAUSSIAN INPUT FILE FOR A WATER MOLECULE USING AVOGADRO

    1. Install Avogadro on your computer, if needed. 
    2. Open Avogadro and Select Element: Oxygen(8)
    3. Click once anywhere is the main window and watch hydrogens being added making a water molecule. 
    4. Prepare Gaussian input file: Open Extensions, Gaussian. 
    5. Note: the coordinates are already listed in the Gaussian Input plugin. 
    6. Enter title: water SP
    7. Select: Calculation: Single Point Energy, Theory level: RHF, Basis: 6-31G(d), Charge: 0, Multiplicity: 1. 
    8. Click: Generate button. 
    9. Save as: water.com

  1. Gaussian SINGLE POINT CALCULATION: INTERACTIVELY ON SCS' LINUX CLUSTER

    1. Copy Gaussian input file to the cluster you are using using scp command. 
    2. > ssh [netid]@lop.scs.illinois.edu 
    3. Make a directory waterSP in your home directory. 
    4. Copy Gaussian input file to cluster using scp command. (scp water.com [netid]@lop.scs.illinois.edu:~/waterSP)
    5. Note: please see the Unix/Linux Primer if you need to refresh your knowledge of the OS. 
    6. Load Gaussian module : 
    7. > module load gaussian/g16
    8. Start Gaussian job with prepared files: 
    9. > g16 < water.com
      • Note: < is the "pipe" flag to direct the input file to Gaussian
    10. Output will be printed on the screen. 
    11. To direct output into a file add > water.log in the end. 
      • Note: proceed to Appendix for an example of submitting a Gaussian calculation in queue. Please make sure to understand each line of the script. 

  1. Gaussian OUTPUT: BASIC FEATURES

    1. Open water.log in a text editor and locate the following: Input data, Standard orientation, Number of basis functions used, Single point energy, Symmetry, Mulliken charges.
      • Note: the energy is reported in A.U. please estimate how many significant digits are needed to guarantee sufficient accuracy. Define sufficient.

Exercise 3.1: Find data in the output without opening. Close the output file and use grep to find the energy value, without opening the file. This can become important for larger systems, as the size of the output file will grow very fast. 

  1. USE Gaussian TO CALCULATE MOLECULAR ORBITALS

    1. Copy water.com into waterMD.com and add in the end of the first line, after SP the following Pop=Reg formcheck. The former requests the data on molecular orbitals to be included into the output and the former asks to print what is called checkpoint file. It is used to restart calculations and to graph molecular orbitals. 
    2. Run the new calculation. Locate in the output file the following: the data mentioned above, molecular orbital coefficients, symmetries and energies of the MO's and HOMO/LUMO orbitals. 
      • Note: the checkpoint Test.FChk was produced. This is an ASCII file which has MO information. 
Exercise 4.1:  Display molecular orbitals using Avogadro. 
  1. Copy Test.FChk to your computer where Avogadro is installed. 
  2. Start Avogadro.
  3. File > OpenTest.FChk. Tow show surfaces go to Extensions > Create Surfaces to plot electrostatic potential select menu flags as shown on figure 1. 
Figure 1: "Create Surfaces" menu
The result will look similar to figure 2.
Figure 2: "Create Surface" output sample image
 
Plot HOMO for water molecule. Select menu flags as shown on figure 3.
Image 3 - "Create Surfaces" menu for water molecule
 
HOMO is shown on figure 4.
Figure 4: HOMO of a water molecule
 
Use Avogardo to graph LUMO of a water molecule.

  1. SUMMARY

This tutorial covered material to have one started using Gaussian for single point energy calculation, analyze the output, graph electrostatic potentials and molecular orbitals.

  1. CONTACT

If you found errors/typos or have suggestions or comments on material in this tutorial please contact us at the SCS Computer Center. We are looking forward to hearing from you.

  1. APPENDIX: RUNNING Gaussian IN QUEUE ON LOP (SCS LINUX CLUSTER)

Please make sure to understand each line of the shell script below. Edit it to match your file names.
  1. > module load gaussian    (Load Gaussian module)
  2. > submit-g16 water.com    (Submit the queue script to the queue)
    • Note: The following is probably not needed for the fast test job.
  3. > qstat     (Check the status of your job)
  4. > qdel [jobID]     (Delete the job if needed, where [jobID] is your job ID in the queue.

Keywords:
quantum chemistry, gaussian, intro
Doc ID:
103204
Owned by:
Taras P. in School of Chemical Sciences
UIUC
Created:
2020-06-18
Updated:
2024-10-28
Sites:
University of Illinois School of Chemical Sciences