Links to
www.Gaussian.com/
G98 Technical information G98 Release notes
-------------------------------------------------------------------------------- How to split large RWF files (>2GB)
%rwf=f1,1500MB,f2,1900MB,f3,1900MB,f4,1900MB,f5,1900MB,f6,1900MB,f7,1900MB %NoSave
%Chk=myoo2.hus.chk %Mem=800MB
#P UB3LYP/Gen Freq SCF(Vshift=0,Conver=6,MaxCycle=300) Geom=AllCheck Guess=Read MaxDisk=13000MB ! Gen basis set
@/teo/ulf/Myo/Vac/O2/B3lyp/myobas/N
-------------------------------------------------------------------------------- How to run excited states with TD
Also open-shell ground state. Use keyword TD (G98 manual p. 155). Chapter 9 (7, p. 173) in Exploring chemistry with electronic structure methods gives a detailed description of the CIS method which has similar input and output. Use a normal SCF calculation input file and add in the final step: TD(Nstates=12,Root=1,50-50) Density Pop=Full
where Nstates is the number of interesting states and Root is the root of interest (for population analysis and geometry optimization), 50-50 indicates that you are interested in both singlet and triplet states. -------------------------------------------------------------------------------- Typical output of a Frequency calculation:
SCF Done: E(RB+HF-LYP) = -76.4070229620 A.U. after 10 cycles Convg = 0.9037D-08 -V/T = 2.0052 S**2 = 0.0000
KE= 7.601163519654D+01 PE=-1.989603930100D+02 EE= 3.745590692757D+01
Leave Link 502 at Wed Mar 28 04:59:08 2001, MaxMem= 6291456 cpu: 7.3 (Enter /usr/local/lib/Gaussian/G98new/g98/l801.exe) Range of M.O.s used for correlation: 1 18
NBasis= 18 NAE= 5 NBE= 5 NFC= 0 NFV= 0 NROrb= 18 NOA= 5 NOB= 5 NVA= 13 NVB= 13
Leave Link 801 at Wed Mar 28 04:59:08 2001, MaxMem= 6291456 cpu: 0.0 (Enter /usr/local/lib/Gaussian/G98new/g98/l1101.exe) Using compressed storage. Will process 3 atoms per pass.
Leave Link 1101 at Wed Mar 28 04:59:10 2001, MaxMem= 6291456 cpu: 1.7 (Enter /usr/local/lib/Gaussian/G98new/g98/l1102.exe)
41
Use density number 0.
Leave Link 1102 at Wed Mar 28 04:59:11 2001, MaxMem= 6291456 cpu: 0.4 (Enter /usr/local/lib/Gaussian/G98new/g98/l1110.exe) forming Gx(P) for the SCF density.
Integral derivatives from FoFDir, PRISM(SPDF).
Do as many integral derivatives as possible in FoFDir. G2DrvN: MDV= 6291456.
G2DrvN: will do 3 atoms at a time, making 1 passes doing MaxLOS=2. Petite list used in FoFDir. MinBra= 0 MaxBra= 2 Meth= 1.
IRaf= 0 NMat= 1 IRICut= 1 DoRegI=T DoRafI=F ISym2E= 1 JSym2E=1. FoFDir used for L=0 through L=2.
Leave Link 1110 at Wed Mar 28 04:59:19 2001, MaxMem= 6291456 cpu: 6.1 (Enter /usr/local/lib/Gaussian/G98new/g98/l1002.exe) Minotr: Closed-shell wavefunction. Direct CPHF calculation.
Solving linear equations simultaneously. Using symmetry in CPHF.
Requested convergence is 1.0D-08 RMS, and 1.0D-07 maximum. Secondary convergence is 1.0D-12 RMS, and 1.0D-12 maximum. Differentiating once with respect to electric field. with respect to dipole field.
Differentiating once with respect to nuclear coordinates. NewPWx=T KeepS1=F KeepF1=F KeepIn=F MapXYZ=F. MDV= 6291456 Using IRadAn= 2.
Store integrals in memory, NReq= 876501. 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. There are 12 degrees of freedom in the 1st order CPHF. 10 vectors were produced by pass 0.
AX will form 10 AO Fock derivatives at one time. 9 vectors were produced by pass 1. 9 vectors were produced by pass 2. 9 vectors were produced by pass 3. 5 vectors were produced by pass 4. 3 vectors were produced by pass 5.
Inv2: IOpt= 1 Iter= 1 AM= 1.09D-15 Conv= 1.00D-12.
Inverted reduced A of dimension 45 with in-core refinement.
Leave Link 1002 at Wed Mar 28 04:59:27 2001, MaxMem= 6291456 cpu: 5.5 (Enter /usr/local/lib/Gaussian/G98new/g98/l601.exe)
Copying SCF densities to generalized density rwf, ISCF=0 IROHF=0.
42
********************************************************************** Population analysis using the SCF density.
********************************************************************** ...
Leave Link 601 at Wed Mar 28 04:59:28 2001, MaxMem= 6291456 cpu: 0.5 (Enter /usr/local/lib/Gaussian/G98new/g98/l701.exe) Compute integral second derivatives.
... and contract with generalized density number 0.
Leave Link 701 at Wed Mar 28 04:59:31 2001, MaxMem= 6291456 cpu: 2.5 (Enter /usr/local/lib/Gaussian/G98new/g98/l702.exe)
L702 exits ... SP integral derivatives will be done elsewhere.
Leave Link 702 at Wed Mar 28 04:59:32 2001, MaxMem= 6291456 cpu: 0.0 (Enter /usr/local/lib/Gaussian/G98new/g98/l703.exe) Compute integral second derivatives.
Integral derivatives from FoFDir, PRISM(SPDF). Petite list used in FoFDir. MinBra= 0 MaxBra= 2 Meth= 1.
IRaf= 0 NMat= 1 IRICut= 1 DoRegI=T DoRafI=F ISym2E= 1 JSym2E=1. Leave Link 703 at Wed Mar 28 04:59:52 2001, MaxMem= 6291456 cpu: 13.9 (Enter /usr/local/lib/Gaussian/G98new/g98/l716.exe)
Dipole = 1.97314156D-05-8.22036304D-01 0.00000000D+00 Polarizability= 7.33514026D+00 1.70957942D-03 5.34051953D+00 2.28341401D-04 1.91679577D-05 2.73823266D+00 Full mass-weighted force constant matrix:
Low frequencies --- -54.4316 -0.0023 -0.0019 -0.0016 46.0086 70.3001 Low frequencies --- 1712.8326 3718.8261 3841.6059 Harmonic frequencies (cm**-1), IR intensities (KM/Mole),
Raman scattering activities (A**4/AMU), Raman depolarization ratios, reduced masses (AMU), force constants (mDyne/A) and normal coordinates: 1 2 3 A' A' A' Frequencies -- 1712.8326 3718.8261 3841.6059
Gaussian软件应用第一到第十章
http://210.34.15.126/cgi-bin/jinghua.cgi?action=list&forum=3&show=75 Gaussian98使用经验谈1－4
Gaussian 98 – some hints for doing calculations
Gaussian 98 is an excellent program that is able to do complex calculations and produce output that can be readily understood. However, like all programs, it can have its little idiosyncrasies that irritate and confuse the user. The best sources of information on how to use the program are the two manuals that are shipped with the program, if you cannot find what you need within those the next best source is the Gaussian website. Unfortunately, there are many
43
simple procedures which are either missing from the manuals or where the instructions are either incorrect or incomprehensible (or both). Therefore you may find the following hints and tips useful. We stress that the information here is intended to supplement the manuals and web-based material rather than replace them.
Contents
Some Things That You Really Should Do Before You Start Getting Rid Of The Error Message When Opening Gaussview Adding A New Basis Set ·For G98W ·For G98 Linux Freezing A Parameter Setting Up A Batch File
Dealing with Negative Frequencies Inexpensive Integrals Geometry Optimizations
How To Produce A Potential Energy Surface How To Get A Good Guess for FC
To Add Extra Basis Functions To Only Certain Specified Atoms Using freqchk
Visualising Molecular Orbitals and Potential Energy Surfaces. Other Useful Things
Some Things That You Really Should Do Before You Start!
Before you start to use Gaussian 98 there are just a few things you should do to try to minimize problems later on. ·Using Windows Explorer, left click on the G98W directory. If G98W is installed in drive C then you should see seven folders, these are; 1.data 2.examples 3.exercise 4.help 5.quick 6.scratch 7.tests
·You should save all your work files into the 'Scratch' directory, doing this means that the program knows where they are. At the end of a successful job you can then move your files into a folder of your choice. ·Don't put spaces in file names or folders as your computer will get lost!
·If a job crashes or if you kill it manually then some files are left behind in the Scratch directory - these all start with the letters 'gxx.'. For example, there will be a file that is named gxx.rwf. These left-over files can take up a lot of space and you can safely delete them. When a job finishes successfully then these files are automatically deleted by G98W.
·When you run a job you will find that you have created three different types of file, these are;
1.the *.gjf file - this is the file that you created by drawing the molecule and setting up the calculation. 2.the *.log file (or *.out file) - this contains all the results and shows the progress of the calculation. You can see an updated log file during the course of a calculation by opening Gaussview, file, open, navigate to the Scratch
44
directory (where you have saved your files!!) and double click on *.log file with the appropriate name. This opens a Gaussview window with your molecule in it, so go to 'results', 'view file' to see how the job is progressing. 3.the *.chk file - this is the most extensive file; it contains all of the detailed electronic information about the molecule and the information that the computer needs to perform different tasks. It is important to save this file, you can use it to restart jobs, visualize MO's as well as using it as the input for a second calculation. This file cannot be read by opening it in WordPad or other such programs, you open it using Gaussview. The key message here is; ALWAYS save the checkpoint file! The default in G98W is not to save this file, so you have to actively tell the program to save it. In the calculation setup box there is a small window called 'Link 0 Commands'. Make sure that you put the file name beside '%chk=', otherwise the .chk file will not be saved.
Getting Rid Of The Error Message When Opening Gaussview
GaussViewW (GVW) was released before Windows 2000 and there is a minor change in the Windows Registry entry for WordPad so that GVW cannot locate it. There are two options, live with the error and modify a GVW script or, two, modify the registry so that GVW can locate WordPad. Both work but if you are uncomfortable modifying the registry the first is safer.
·Option 1 Edit the script g98w\\editor.bat and change the line start/wait notepad %2 to start/wait notepad %1 For large files notepad can activate WordPad.
·Option 2 Use the RegEdit utility in Windows to locate the entry for word pad,
HKEY_LOCAL_MACHINE\\SOFTWARE\\Microsoft\\Windows\\CurrentVersion\\AppPaths\\WORDPAD.EXE and change the type from REG_EXPAND_SZ to REG_SZ. Then the string \and %ProgramFiles% replaced by the actual path to the WindowsNT directory. Exit RegEdit and GVW should then locate WordPad correctly. As should other applications.
Adding A New Basis Set
When adding a new basis set you should only include definitions for atoms that are actually in the molecule. For G98W
1.Draw molecule as normal - using gaussview 2.Select 'Gaussian' on the 'calculate' menu
3.Select your usual parameters in the calculation setup dialog box, but select a 'blank' for the basis set. 4.In the 'additional keywords' box type in the model chemistry together with the keyword gen, e.g. rb3lyp/gen 5.Click on edit, you'll be asked to save the file first.
6.Add the text for the new basis set in the 'molecule specification' box at the end of the existing text - leave a blank line between the sets of text. Check that the route section is correct, there may be two instances of the model chemistry. Save 7.Click okay 8.Submit For G98 Linux
1.Draw molecule as normal - using gaussview 2.Select 'Gaussian' on the 'calculate' menu
3.Select the usual parameters in the calculation setup dialog box, but select a 'blank' for the basis set. 4.In the 'additional keywords' box type in the model chemistry together with the keyword gen, e.g. rb3lyp/gen 5.Click 'edit' then go to the bottom of the existing text and place the cursor there. Then select 'include file' from the 'file' menu and navigate to the new basis set, select it - it will be added to the input file - remove the '!' symbol and do not leave a blank line between text. Make sure that in the route section there is only one instance of
45