wiki.chpc.ac.za

User Tools

Site Tools

You are here: start » howto » gaussian
Trace:
howto:gaussian

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision 		Previous revision
howto:gaussian [2013/07/16 13:18]
wikiadmin [Example moab job script for g09.A01] 		howto:gaussian [2025/10/17 09:51] (current)
alopis [Easy Gaussian09 and Gaussian16 PBS creation and submission]
Line 1: Line 1:
 ====== Gaussian at CHPC ====== ====== Gaussian at CHPC ======
  
-  * **NOTE : You should always run your jobs from ''scratch5''**+This license can only be used by CHPC registered users, which are associated with either **South African** or **SKA partner** higher education institutions
  
 +Gaussian can be run over 1 or 2 nodes (see example below) 
  
-We have two versions of Gaussian 09 installed at CHPC. +====Easy Gaussian09 and Gaussian16 PBS creation and submission====
-Here is an example on how to access them:+
  
-  embele@login01:~/scratch5 $ module avail                   ### list available modules +In order to see all available Easy Scripts, first include/load this module - even better if your add this line to your .bashrc:
-  embele@login01:~/scratch5 $ module add gaussian/g09.A01    ### load g09 version A01 (older) or +
-  embele@login01:~/scratch5 $ module add gaussian/g09.D01    ### load g09 version D01 (new) +
-  +
-  +
-====Example moab job script for g09.A01====+
  
 +  module load chpc/easy_scripts
  
-  #/bin/csh +Execute **qg09_E01** or **qg16_A03** or **qg16_B01** or **qg16_C01**
-  #MSUB -l nodes=1:ppn=12 +
-  #MSUB -l feature=dell +
-  #MSUB -l walltime=2:00:00 +
-  #MSUB -m be +
-  #MSUB -V +
-  #MSUB -o /export/home/username/scratch5/log.out +
-  #MSUB -e /export/home/username/scratch5/log.err +
-  #MSUB -d /export/home/username/scratch5 +
-    +
-  source /opt/gridware/applications/gaussian/old/g09/g09setup  +
-  source /etc/profile.d/modules.sh +
-  module add gaussian/g09.A01 +
-  g09 < input.com > output.log +
-  +
  
-====Example moab job script for g09.D01==== +This script handles PBS setup and submission. 
- +Also takes care of jobs running over more than one node with LINDA and does cleanup of temporary GAU files when job is completed/killed. Examples of these scripts are provided below.
  
-  #/bin/csh +              EXAMPLE1 
-  #MSUB -l nodes=1:ppn=12 +  Enter project name/shortname 
-  #MSUB -l feature=dell +  CHEM0100 
-  #MSUB -l walltime=2:00:00 +  Enter input file (with .com or .gjf extension) 
-  #MSUB -m be +  test.gjf 
-  #MSUB -V +  Enter number of nodes on which to run job 
-  #MSUB -o /export/home/username/scratch5/log.out +  1 
-  #MSUB -e /export/home/username/scratch5/log.err +  Enter total walltime (hour:minute) 
-  #MSUB -d /export/home/username/scratch5 +  2:00 
-    +  Enter email address 
-      +  testing@gmail.com 
-  source /opt/gridware/applications/gaussian/g09/g09setup +  Generated pbs file for test 
-  source /etc/profile.d/modules.sh +  Do you wish to submit job to cluster (y/n) 
-  module add gaussian/g09.D01 +  y
-  g09 < input.com > output.log+
  
  
-Both examples above use this example of a gaussian input file:+    EXAMPLE2 [PLEASE TAKE NOTE OF EMPTY SPACES] 
 +  Enter project name/shortname 
 +  CHEM0100 
 +  Enter input file (with .com or .gjf extension) 
 +  test.gjf 
 +  Enter number of nodes on which to run job
      
-  %nprocshared=12 +  Enter total walltime (hour:minute) 
-  %nprocl=1 +   
-  #P HF/6-31G*  IOP(6/33=2,6/41=10,6/42=17) SCF=Tight Pop=MK        +  Enter email address 
-  Title Card Required +  testing@gmail.com 
-  0 1 +  Generated pbs file for test 
-   C                 -3.19550100    0.11344600   -0.18511100 +  Do you wish to submit job to cluster (y/n) 
-   O                 -3.05859100    0.83554400   -1.15941100 +  y 
-   N                 -2.14200500   -0.25446800    0.60848900 +   
-   H                 -2.32841200   -0.78506500    1.44398900 +   
-   C                 -0.75900100    0.12341400    0.36048900 +If you wish to use the utility **formchk** (to create readable format files for visualisers) please do the following. 
-   C                  0.14379200   -0.42259800    1.48678900 +   
-   C                  0.59047400   -1.84630300    1.03238900 +  module load chpc/gaussian16/C01 
-   C                  1.59689900    0.10938400    1.28248900 +  source $g16root/g16/bsd/g16.profile 
-   C                  0.08697100   -2.05959700   -0.42431100 +   
-   C                  2.05278100   -1.31402200    0.82788900 +Then the command **formcheck** should be available for use. Other Gaussian commands should also be available, such as **freqchk** and **chkchk**. 
-   C                 -0.14011000   -0.59319400   -0.87441100 + 
-   C                  1.35866600   -2.48501300   -1.17731100 + 
-   C                  1.32989700   -0.05621300   -1.07501100 + 
-   C                  1.66721100    0.99508300    0.02128900 +You can use a **checkpoint file (.chk)** which will allow you to continue with a 2nd (or subsequent) calculation from close to where the initial/previous calculation ended. This is especially helpful for jobs which exceed walltime limits or have stopped for another reason. There is info online including from Gaussian website, these links may be helpful too: 
-   C                  2.20428100   -1.28752400   -0.71551100 + 
-   C                  0.60162500    2.08159700    0.02358900 +https://gaussian.com/faq2/ 
-   N                 -0.63668200    1.58111300    0.28218900 + 
-   O                  0.86964000    3.24909300   -0.21751100 +https://gaussian.com/restart/ 
-   H                 -1.43757500    2.10622300   -0.04521100 + 
-   O                  2.93691800    1.56466600   -0.15981100 +https://gaussian.com/tip4/ 
-   H                  2.76333000    2.51016900   -0.30861100 + 
-   H                 -0.28050600   -0.28619200    2.48278900 +https://www.researchgate.net/post/How-can-I-res... 
-   H                  0.43696300   -2.68340100    1.71298900 + 
-   H                  2.07960500    0.56667700    2.14378900 +Very occassionally  a checkpoint file may be corrupted, so keeping the most recent backup copy of a good checkpoint (.chk) file can be helpful. Otherwise in these circumstances the calculation may need to be started from scratch - or the most recent molecular geometry (the coordinated of all the atoms) could be used from the Gaussian output (.log) file. 
-   H                 -0.78403700   -2.70728500   -0.52611100 +
-   H                  2.86017500   -1.79913300    1.37408900 +
-   H                 -0.75990800   -0.47138600   -1.76101100 +
-   H                  1.22276500   -2.52601100   -2.26181100 +
-   H                  1.75375300   -3.44771800   -0.83891100 +
-   H                  1.50280200    0.34398500   -2.07391100 +
-   H                  3.22968200   -1.22473700   -1.07621100 +
-   C                 -4.55170800   -0.42883700    0.22928900 +
-   H                 -5.23089700    0.41157200    0.38338900 +
-   H                 -4.52881600   -1.04383700    1.13088900 +
-   H                 -4.94951500   -1.02253200   -0.59631100+
  
 +Please **recheck** your molecular **geometry very carefully** to ensure there are no issues - this will result in calculations which take much longer than necessary.
/app/dokuwiki/data/attic/howto/gaussian.1373973490.txt.gz · Last modified: 2021/12/09 16:42 (external edit)

Page Tools

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Noncommercial-Share Alike 3.0 Unported
CC Attribution-Noncommercial-Share Alike 3.0 Unported Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki