GAUSSIAN: Difference between revisions
| Line 114: | Line 114: | ||
http://aliga.iciq.es/wiki/images/files/g98/index.htm | http://aliga.iciq.es/wiki/images/files/g98/index.htm | ||
== Gaussian09 vs Gaussian03 Scaling == | == Gaussian09 vs Gaussian03 Scaling and New Functional Testing == | ||
Tests calculations performed at the CESCA supercomputer. | Tests calculations performed at the CESCA supercomputer. | ||
'''Calcfc, opt to TS and frequency.''' | |||
Job has 687 basis functions with B3LYP method. | Job has 687 basis functions with B3LYP method. | ||
Step G03 G09 | Step G03 G09 | ||
| Line 130: | Line 131: | ||
SCF takes longer, but it is using the new GEDIIS algorithm (vs GDIIS in G03 which now doesn't exist). | SCF takes longer, but it is using the new GEDIIS algorithm (vs GDIIS in G03 which now doesn't exist). | ||
They say GEDIIS should give better performance, especially for calculations not so close to convergence as this example. | They say GEDIIS should give better performance, especially for calculations not so close to convergence as this example. | ||
''' | |||
Testing new functionals with G09''' | |||
Job has 687 basis functions, starting structure from B3LYP method in G03. | |||
Step Functional | |||
M06 B97D | |||
calcfc to 2nd l103 4h 49m 34s 1h 59m 42s | |||
2nd opt step 13m 36s 7m 42s | |||
3rd opt step 13m 37s 7m 15s | |||
Total opt time 19h 05m 32s 3h 36m 35s | |||
No. steps 89 32 | |||
avg opt step time 772s 406s | |||
freq 5h 27m 59s 2h 11m 48s | |||
Total calc time 29h 50m 18s 7h 48m 5s | |||
The M06 functional is 33% slower than B3LYP (based on the frequency times), but still faster than B3LYP in G03. | |||
Grimme's B97D functional is very fast, and also takes less steps to optimize in this case too. | |||
Revision as of 16:27, 8 October 2009
go back to Main Page, Computational Resources, Chemistry & More, Computational Codes
Gaussian 09 is the latest in the Gaussian series of electronic structure programs. Gaussian is used by chemists, chemical engineers, biochemists, physicists and others for research in established and emerging areas of chemical interest.
Starting from the basic laws of quantum mechanics, Gaussian predicts the energies, molecular structures, and vibrational frequencies of molecular systems, along with numerous molecular properties derived from these basic computation types. It can be used to study molecules and reactions under a wide range of conditions, including both stable species and compounds which are difficult or impossible to observe experimentally such as short-lived intermediates and transition structures. This article introduces several of its new and enhanced features.
GAUSSIAN was first written by John Pople and released in 1970. According to the manual, it can deal with different computational approaches: molecular mechanics (AMBER, UFF, DREIDING), semi-empirical calculations (AM1,PM3,CNDO,INDO,MINDO/3,MNDO), SCF methods (Restricted, Unrestricted, and Restricted Open-shell Hartree-Fock), Møller-Plesset perturbation theory, DFT methods (Hybrid functionals, exchange functionals and correlation functionals), ONIOM (QM/MM method), Complete Active Space (CAS) and Multi-Configurational Self-Consistent Field calculations, Coupled Cluster calculations, QCI methods and Quantum chemistry composite methods.
Nearly everything you need to know can be found on the Gaussian's homepage, 1 and the on-line manual, [2]]
Installing Gaussian in your local computer
It is sometimes useful to be able to run short calculations on your local machine. Although this should not be done systematically to avoid burning out the machine.
One option is to compile the code (potentially painful).
Another option is to copy a pre-compiled version following the instructions below. It will not run optimally for your machine, but ...
Installing Gaussian on your local computer
Gaussian for beginners
If you are not familiar with Gaussian and Gaussview, it might be a good idea to run some "simple" calculations before starting your "own" project. Following the ling below you will find some exercices prepared for the Master on Computational Chemistry, which are a good starting point.
Graphical Interfaces
link to the GaussView help : [[1]] and Gauss View in the wiki Gauss View
Input Examples
In spite of the highly informative Gaussian homepage, it can sometimes be difficult for newcomers to get a quick overview of the program. In order to make the introduction to the program easier, input examples for several kinds of calculations will be presented in the following. In this way users new to the program can just copy, paste and modify the input examples and they will have a calculation ready to submit.
Optimisation from transition state geometries
Gaussian questions & answers
Links
http://www.gaussian.com/tech_top_level.htm
http://en.wikipedia.org/wiki/GAUSSIAN
Introduction and input/output explanation of PCM
http://www.cup.uni-muenchen.de/oc/zipse/compchem/solv/pcm.html
ABOUT SENDING CALCULATIONS :
- Checkpoints: They are big files which contain information such as the wave function or the geometry of finished calculations. They can not be edited with vi or nedit... (as they are written in binary) but the information they contain can be read by Gaussian. Although they might be very useful, to store them is not recommended if we are not sure about their utility.
In kimik, when the calculation finishes they are saved in the same directory as the output. For practical reasons while the calculation is running the checkpoint file should be saved in the Scratch. The submision script qs adds the lines
%chk=/scratch/filename.chk %nprocs=XX (XX corresponds the type of node (single or quad core) you requested in qs, where "qs 1 (s/n/l) g03 filename.in" -> %nprocs=1, and "qs 1 (q) g03 filename.in" -> %nprocs=4) %nproclinda=XX (XX corresponds to the number of nodes you requested in qs, where "qs 1 l g03 filename.in" will not add the line, but "qs 4 l g03 filename.in" -> %nproclinda=4)
automatically at the beginning of your input file creating a new input called filename.in.in. You should not include these lines in your input file when using qs, if not you will end up with two %chk= lines and two %nproc=.
- Memory: The qs script does not add a line to specify memory usage, which must be done manually with the "%Mem=" command. The single core nodes have 1GB of RAM, and the quads have 4GB (but only 2GB is accessible due to the 32bit version of Gaussian 03 currently compiled).
As a rule you should specify no more than 85% of a node's memory capacity (single=850MB, quad=1950MB) to allow for normal OS operations, like disk writing and network access, otherwise you will start using swap space and slow down the calculation drastically (especially in frequency calculations)
- --link1--; you can link calculations by writing the inputs one after the other and adding --link1-- between them (one blank line before link1 and no blank line after). With this and the keywords geom=(checkpoint), guess=(read)... you could save time when doing solvent calculations or other. The outputs will be written in a single file.
ABOUT NBO: Natural Bond Orbital Analysis
It is recommended to use the keyword pop=nboread in the command line and at the end of the file $NBO 3CBONDS RESONANCE $END
http://www.chem.wisc.edu/~nbo5/ch3nh2.html
ABOUT PCM (CPCM, IPCM,..) :
- If you have one or more hydrogens "bonded" to two atoms (close to two atoms), the default parameters of PCM will not be able to compute the molecular cavity for them. There are two solutions for this problem:
-1- Use a different model to build up the cavity. Change the default UAO which does not explicitly consider hydrogens, by UHF, PAULING, BONDI.. which consider hydrogens explicitly.
-2- Use the keyword SPHEREONH = atom number of the hydrogen "bonded" to two atoms.
- Often Gaussian has problems computing the cavity in specific structures, giving an error message like the one below...
AdVTs1: ISph= 3760 is engulfed by JSph= 3762 but Ae( 3760) is not yet zero!
There are two methods to solve this:
-1- Increase the RMIN value used in the calculation. Increasing the RMIN value smooths the surface of the spheres that make the cavity and helps prevent overlapping. The default value in Gaussian 03 is RMIN=0.20, but you can use any value up to 0.40. You should try to keep the new RMIN within 0.1 vs the default, as a larger difference will give errors in the solvation free energy to the 4th decimal place (from 0.05 kcal/mol upwards).
-2- Try to modify slightly your coordinates. A very small change of the coordinates may be enough to eliminate the problem. A good way to do that is to round the coordinates (i.e. from 2.347895 to 2.348) or re-optimize the molecule (the coordinates you will get will be slightly different). Even though that sounds quite stupid some times it is enough for the program to work.
Scripts
Old version of Gaussian, G98
http://aliga.iciq.es/wiki/images/files/g98/index.htm
Gaussian09 vs Gaussian03 Scaling and New Functional Testing
Tests calculations performed at the CESCA supercomputer.
Calcfc, opt to TS and frequency. Job has 687 basis functions with B3LYP method.
Step G03 G09 calcfc to 2nd l103 7h 54m 08s 3h 42m 10s 2nd opt step 8m 00s 9m 16s 3rd opt step 6m 06s 7m 06s freq 8h 30m 25s 4h 6m 15s Total calc time 16h 38m 39s 8h 4m 47s
So frequency is about 2x faster in the new Gaussian. SCF takes longer, but it is using the new GEDIIS algorithm (vs GDIIS in G03 which now doesn't exist). They say GEDIIS should give better performance, especially for calculations not so close to convergence as this example. Testing new functionals with G09 Job has 687 basis functions, starting structure from B3LYP method in G03.
Step Functional
M06 B97D
calcfc to 2nd l103 4h 49m 34s 1h 59m 42s 2nd opt step 13m 36s 7m 42s 3rd opt step 13m 37s 7m 15s Total opt time 19h 05m 32s 3h 36m 35s No. steps 89 32 avg opt step time 772s 406s freq 5h 27m 59s 2h 11m 48s Total calc time 29h 50m 18s 7h 48m 5s
The M06 functional is 33% slower than B3LYP (based on the frequency times), but still faster than B3LYP in G03. Grimme's B97D functional is very fast, and also takes less steps to optimize in this case too.