ONIOM tips

go back to Main Page, Computational Resources, Chemistry & More, Computational Codes, GAUSSIAN, G09

ONIOM tips:

• You can use iop(4/33=2) to print out all MM contributions.
• To get out the model system from an oniom calculation include iop(1/33=2). In the output file, search for: "ONIOM: generating point 2 high level on model system". All atoms without the label -1 belong to the model system (this list includes also the link hydrogens).
• There are two options if you want to check other force fields which are not implemented in Gaussian09;

Option1: You can do one-dimensional potential energy surface scans over a reaction coordinate, where the constrained structure optimisations at the potential energy surface can be performed with ONIOM(DFT:UFF), ONIOM(DFT:DREIDING) or ONIOM (DFT:AMBER), where the iop(1/33=2) option helps to obtain QM and MM parts of the molecule. Then you can do MM optimisations on the QM and MM parts of the molecule separately with MACROMODEL or TINKER, where they have various force filed parameters.

Option2: ´External´ keyword.

There are two examples in the G09 test suite which make use of the ´External´ keyword.

(1) Test726 uses a shell script to run g09 as an external process with UFF. The script ´extuffex´ is stored in the g09/bsd directory.

(2) Test769 uses a perl script, mygau, which is stored in the g09/bsd subdirectory. Note that you will likely need to edit mygau to change the path to perl, rather than /usr/local/bin, because most machines store it in /usr/bin.

In both cases you could use ´External´ keyword in the input file.

Eg:- ONIOM(B3LYP/6-31G(d):AM1:External="Script_name") Opt

How to use external MM program:

This is not fully implemented in G09. The Gau_External script in the G09 home is an example where we use Gaussian as an external program, specifically 'Gau_External' script helps to run MOPAC as external procedure for ONIOM. You would need to make a copy of this script and modify it to create an input file for other MM programes (Eg. TINKER), and run the MM program with this input, and then use the output from the MM program to generate result files suitable for processing by Gaussian.

• In Gaussview: When selecting MM regions in QM/MM methods you can change the status (H/L) of several atoms at same time, by:i) select the upper atom in the list ii) go to the last one and press control+shift 3) change the last atom status to low level.
• GaussView can help assigning atom types and it will also add point charges according to the AMBER convention for atoms in standard protein residues. For other molecules that are not standard residues, there is no information about appropriate charges (GaussView will try to assign reasonable atom types), so one would need to derive charges for the atoms in such molecule, or look in the literature to see if other researchers have derived charges and parameters for such molecule.
• Only the "amber96" version is internally hardwired in the G09 code and that is the one selected when using the AMBER keyword. The "amber.prm" file in the "g09" directory is simply a text file that shows this same parameter file that is internally hardwired in the code, but the program does not read that file. It is included for illustrative purposes so one see examples of how to specify MM parameters manually. The files that you mentioned (amber98.prm and oplsaa.prm) unfortunately are not ready to be used.

You can make modifications to the internally hardwired parameter file by adding the modifications to the input file, for instance: