Introduction to IQmol

v2.13 (2019) Andrew Gilbert

Contents
1 Introduction

4

1.1

Installation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4

1.2

Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

52 Building Molecules

7

2.1

Adding Atoms and Fragments . . . . . . . . . . . . . . . . . . . . . . . . .

7

2.2

Optimizing Structures Using Molecular Mechanics . . . . . . . . . . . . . .

8

2.3

Symmetrizing Structures . . . . . . . . . . . . . . . . . . . . . . . . . . . .

9

2.4

Specifying Geometric Parameters . . . . . . . . . . . . . . . . . . . . . . .

9

2.5

Manipulating and Selecting Molecules . . . . . . . . . . . . . . . . . . . . . 10

2.6

Adding Additional Conformers . . . . . . . . . . . . . . . . . . . . . . . . . 11

3 Running Q-Chem Calculations

12

3.1

The QUI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

3.2

The Job Monitor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

3.3

Configuring Servers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

4 Analyzing Results

15

4.1

Molecular Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

4.2

Plotting Molecular Orbitals From a Checkpoint File . . . . . . . . . . . . . 16

4.3

Other Orbitals and Densities . . . . . . . . . . . . . . . . . . . . . . . . . . 17

4.4

Exporting Cube File Data . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

4.5

Visualizing Cube File Data . . . . . . . . . . . . . . . . . . . . . . . . . . . 19

4.6

Orbital Animations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20

4.7

Potential Energy Surfaces . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

2

4.8

Vibrational Frequencies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.8.1

4.9

Isotopic Substitution . . . . . . . . . . . . . . . . . . . . . . . . . . 23

NMR Spectra . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

5 Appearance

27

5.1

Viewer Camera . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

5.2

Clipping Planes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

5.3

Shaders . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28

5.4

Exporting POV-Ray Files . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

6 Sample Images

30

3

1

Introduction

IQmol is an open-source molecular editor and visualization package that runs underWindows, Mac OS X and Linux. It can read a variety of chemical file formats including
xyz, cml, pdb, mol, fchk, cube data and Q-Chem input/output. It also includes a freeform molecular builder that allows arbitrary molecular structures to be created. These
structures can be optimized using molecular mechanics force fields and symmetrized toensure the structure has the correct point group symmetry. A library of molecules andfunctional groups also exists, and these can be used to facilitate building more complicatedmolecules.
IQmol is capable of displaying a variety of molecular properties including atomic charges,
dipole moments and normal modes. Several surface types can be displayed including
molecular orbitals, (spin) densities and van der Waals surfaces. These surfaces can be colored according to an arbitrary scalar field, such as the electrostatic potential. Animations
are also available for vibrational frequencies, and reaction and optimization pathways.
IQmol can operate as a stand-alone package, but has also been written to work seamlessly with the Q-Chem (http://www.q-chem.com) computational chemistry package. A
comprehensive input file generator, the Q-Chem User Interface (QUI), provides accessto most of the available options in Q-Chem, and these options are presented in an intuitive, hierarchical fashion. The generated input files can be submitted to either local
or remote servers that have the Q-Chem software installed. In particular, a publicly accessible server is available that allows small (limited to approximately 10 mins) quantum
chemistry calculations to be run without having to purchase and install Q-Chem.

1.1

Installation

The latest version of IQmol can be downloaded from the website:
http://iqmol.org/downloads.html
For Mac (OS X) a disk image file is provided. After downloading, Simply double-click
the disk image to mount it and copy the application to the Applications directory, or any
other desired location.
For Windows an installer is provided that will guide you through the installation process
and will also create a shortcut to the application on your desktop.
For Linux .deb and .rpm packages are now available that can be installed using either
dpkg or yum:#> sudo dpkg -i iqmol 2.9.0.deb#> sudo apt-get install -f
4

The second command is required to resolve the dependency on the Qt libraries, which
you may not have installed.
#> sudo yum install iqmol-2.9.9-2.x86 64.rpm
Note that in either case you will need root permission.

1.2

Overview

Figure 1.1: The main IQmol window.

The main IQmol window is shown in Fig. 1.1 and comprises the following main parts:• The Viewer is the main part of the window and is where you can view and interact
with your molecule.
• The Toolbar at the top provides access to common commands and allows selection between different viewer modes including manipulation, selection

andbuilding.
• The Model View panel provides a hierarchical view of the data that are available
for the molecule.
• The History panel in the lower left shows a list of the most recent actions than canbe undone, either by clicking on them or by using the Edit I Undo menu option.
5

The Model View (MV) provides control over what objects are displayed in the Viewer,
and also allows access to configuration options for these objects. Visibility is controlled
by the associated check-box. Unchecking a check-box causes the item and all its children
to become hidden. If an item does not have a check-box, for example a bond, then its
visibility can only be controlled by items higher in the hierarchy.
The appearance of many objects can be configured by double-clicking the item in the MV.
For example, double-clicking the molecule name brings up the Configure Molecule dialog
which allows you to change the appearance of the molecular structure:

Note that several molecules can be viewed concurrently in the same Viewer and theappearance of each can be configured separately.

Figure 1.2: Rendering styles for molecular structures: balls and sticks, wire frame, spacefilling, tubes and plastic.

6

2

Building Molecules

2.1

Adding Atoms and Fragments

By default, IQmol opens in build mode. This is indicated by a red border around the
button in the Toolbar. The default build atom is indicated by the
button and canbe changed by clicking on this button. A pop-up periodic table will appear from which
the desired element type can be selected.
Clicking in the empty Viewer window will create an atom of the current build element.
Additional atoms can be added by clicking on an existing atom and dragging the mouse.
This creates a new atom bonded to the first. To create a disconnected atom, hold down
the alt modifier key when clicking in the Viewer window. (Note that some Linux window
managers use the alt modifier for other purposes. This behavior can be changed using
System Settings I Keyboard I Short-cuts menu).
Bond orders can be increased by clicking and dragging between two existing atoms. If nobond exists between the two atoms, one is created. Otherwise the bond order is increased.
To decrease the bond order, the bond must first be deleted and a new bond created.

Figure 2.1: Increasing the bond order.

button, ensuring the Functional Group
Functional groups can be added by clicking theradio button is clicked and selecting the desired group from the menu. Groups are addedin the same way as atoms, i.e. clicking and dragging from an existing atom. The empty
valence is indicated by the yellow bond and shows where the group will be connected.
Entire molecules can also be added to your system by clicking the
button, ensuring
the Molecules radio
buttonisclicked,andselectingthedesiredmolecule
from the menu.
Be sure to click Select or the build selection will not be updated. Unlike the other build
modes, clicking anywhere in the viewer window will add the selected molecule (no mousemodifier is required). This allows several molecules of the same type to be added quickly,(e.g. for solvation), but changes the usual mouse behavior. If you accidentally add too
many molecules, use the Edit I Undo menu option. When adding molecules, if you click
and hold it is possible to alter the local orientation of the newly added molecule before
being added to the global frame.
7

Figure 2.2: The build fragment pop-up that appears when the add fragment button isclicked.

Once the backbone of the molecule has been drawn, the Add Hydrogens button
be clicked to automatically add hydrogen atoms to any unfilled valencies.

2.2

can

Optimizing Structures Using Molecular Mechanics

The builder in IQmol is free-form, so your initial structure may look a bit wonky. To
improve the geometry, click the
button. This optimizes the geometry using a molecularmechanics (MM) force field. The default force field is the Universal Force Field (UFF)[2]
which has the advantage of being defined for almost all of the periodic table. However,the UFF does not perform well for systems that contain hydrogen bonds and, in thesecases, it is recommended that the force field be changed using the Build I Select Force
Field menu option.

Figure 2.3: Changing the molecular mechanics force field.

8

2.3

Symmetrizing Structures

If your molecule has symmetry, the MM optimization is unlikely to find a structure withthe desired symmetry. If this is the case, a nearly symmetric structure can be symmetrizedusing the Build I Symmetrize Molecule menu option. Finding very high symmetry may
require relaxing the tolerance using the Build I Set Symmetry Tolerance menu option.
Relaxing the tolerance allows the program to move the nuclear coordinates more in order
to find a symmetric structure.
IQmol uses a modified version of the Symmol[1] program written by Tullio Pilati andAlessandra Forni to symmetrize molecular structures.

2.4

Specifying Geometric Parameters

Specific values for geometric parameters can be set by first selecting the atoms involved
and using the Build I Set Geometric Constraint menu option. A dialog will appear thatallows the parameter to be either set, constrained or scanned. Constrained parameters
apply to any subsequent MM optimization and are also passed through to the Q-Chem
input file, if a optimization job is requested. Scan options are also passed through to
Q-Chem for scan jobs.

Figure 2.4: Dialog for setting geometric constraints
The type of constraint depends on the number of atoms selected:
1. Fixed atom position
2. Inter-atomic distance (or select a single bond)
3. Bond angle
4. Torsion (dihedral) angle
Active constraints are visible in the Viewer and can be deactivated by clicking the adjacent
check-box in the MV.
92.5

Manipulating and Selecting Molecules

IQmol has been designed to work best with a three-button mouse or track pad. If you
are using a track pad on a Mac, it is recommended that you enable the secondary clickby going to System Preferences I Trackpad.
Manipulate mode is activated by clicking thethe following mouse functionality:

button in the Toolbar, and implements• Left click and drag: Rotate the view of the molecule.
• Middle click and drag: Zoom in and out.
• Right click and drag: Translate the view of the molecule.
It is also possible to manipulate part of the molecule independently from the rest. To dothis, make a selection and press and hold the ctrl modifier (command key on Mac). The
mouse movements will affect only the selected atoms as follows:• Left click and drag: Rotate the selected atoms about their center.
• Right click and drag: Translate the selected atoms.
If only a single bond is selected, then the mouse movements have the following effects• Left click and drag: Rotate around the axis of the bond.
• Right click and drag: Change the length of the bond.
Select mode is activated by clicking thefollowing mouse functionality:

button in the Toolbar, and implements the• Left click: Adds atom or bond to selection.
• Click and drag: Creates a selection rectangle, all atoms and bonds within theselection rectangle are added to the selection.
• Right click: Removes atom or bond from selection.
It is also possible to select all, select none and invert the selection via the options under
the Edit menu.

10

2.6

Adding Additional Conformers

For setting up jobs that use the Freezing String Method (FSM), two conformers are
required corresponding to the initial and final geometries of the string. Having built thefirst geometry, a second geometry can be added to the molecule by right-clicking themolecule name in the MV. This brings up a context menu with the option ‘Duplicate
Geometry’. Selecting this menu option creates a Geometries item in the MV with two(identical) geometries. The second of these can be modified by selecting it in the MV andmanipulating selected atoms in the viewer using the ctrl (command on Mac) modifier key.
If ‘Freezing String’ is chosen for the Calculate option in the QUI (see the next section)then both geometries will be included in the $molecule section of the input deck, separated
by ****.

11

3

Running Q-Chem Calculations

3.1

The QUI

IQmol has a built-in input file generator for Q-Chem calculations, the QUI, that can beaccessed via the Calculation I Q-Chem Setup menu.

Figure 3.1: The Q-Chem User Interface (QUI) dialog.

The left hand side of the QUI dialog contains controls for setting up the calculation. Theseare presented in a hierarchical fashion so that most commonly used options are presented
in the top part of the panel, with other relevant options appearing in the lower section
depending on what type of job is selected. More advanced options can be accessed via
the Advanced tab. The generated input file is echoed in the panel on the right hand sideof the QUI.
Clicking the Submit button will start the job running on the selected server. Once the job
has completed you will be prompted to copy the results (for remote servers) before they
are automatically loaded into IQmol. The icon appears next to the molecule name in
the MV when the molecule has been updated with the results of a calculation.

12

Figure 3.2: The Job Monitor dialog. Right-clicking on a selected job brings up a context
menu with additional options.

3.2

The Job Monitor

Submitted jobs can be monitored via the Calculation I Job Monitor menu option. This
brings up the Job Monitor dialog that provides information on the progress of jobs. Rightclicking a row in the Job Monitor brings up a context menu which allows running jobs tobe killed or queried, and finished jobs to be opened.
If a job has the Error status, hovering the mouse over the status will show the error
message. Alternatively, double clicking the job will open the output file, if available, andmay provide additional information about why the job failed. Double-clicking a completed
job will copy the results from the server (if they have not already been copied) and reload
them into the main IQmol window.

3.3

Configuring Servers

By default IQmol is configured to submit jobs to the Q-Chem server in Pleasanton,
California. This is a publicly available server than can be used by anyone wishing to run
test calculations before purchasing the Q-Chem software. Jobs submitted to this server
can access the full suite of electronic structure methods available in Q-Chem, but are
time-limited to 10 minutes.
Additional servers can be configured to access other computers that have the Q-Chem
software installed. These computers can be either local servers (i.e. the same machineas the one running IQmol) or remote servers that can be connected to via SSH. To add
13

Figure 3.3: The server configuration dialog.

an additional server, go to the Calculation I Edit Servers menu option and click the
button.
The information required to configure a server depends on its type, but default options
exist for each case. The following indicates the minimum that should be considered:• Local: Set the Queue System to ‘Basic’ unless you know you are running queuing
software on the machine. Be sure the QC and QCSCRATCH variables are set to theircorrect values in the Run File Template, which is accessed by clicking the Configure
button.
• SSH: The host name and account on the target machine will be required in order toconnect via SSH. Set the Authentication combo-box to ‘Password Prompt’ unless
you have set up an alternative authentication protocol. Depending on how thetarget host has been set up, the Run File Template may require some editing to get
it to work, but this will differ from machine to machine.
• HTTP: The default options should be suitable. Note that on Windows HTTP
should be used (and is the default in v2.8) whereas for Mac and Linux HTTPS canbe used and will be the default in future versions.

14

4Analyzing Results

IQmol reads in a wide variety of file types and allows the user to visualize many of theresults they contain. Individual files can be opened via the File I Open menu option, or
alternatively by dragging and dropping the file onto the Viewer window. Directories canalso be opened via the File I Open Dir menu option, or by dragging and dropping thedirectory onto the Viewer.
Opening directories allows more than one file associated with a molecule to be loaded
together, such as a Q-Chem output file and a formatted checkpoint file. The directory
name determines the base name for the molecule and IQmol loads all files contained
within the directory that have a matching base name.

Figure 4.1: Example: Opening the Aniline directory will load the Aniline.ESP.cube,
Aniline.FChk, Aniline.inp and Aniline.out files into the one molecule in the MV. The
notes.txt file will be ignored as it has a different base name.

4.1

Molecular Surfaces

Several molecular surfaces can be generated without first performing a quantum chemical
calculation. These include: van der Waals, pro-molecule and superposition of ionic densities (SID). The SID surface is similar to the pro-molecule surface except that it scales
the densities based on the atomic charges of the atoms and may provide a more accurate
representation of the density of charged systems.
To plot one of these pseudo-density surfaces, double-click the ‘Surfaces’ item in the MV
associated with the molecule.

15

4.2

Plotting Molecular Orbitals From a Checkpoint File

Plotting molecular orbitals (MOs) requires a formatted checkpoint file (.fchk file extension) which is generated by default when running Q-Chem calculations from IQmol(GUI rem variable should be set to 2). After opening the fchk file, a ‘Canonical Orbitals’
item will appear in the MV under the Surfaces item. Double-clicking this brings up the
Add Surface dialog:

Figure 4.2: The Add Surface dialog allows MOs and (spin-)densities to be plotted.

The Add Surface dialog allows MOs, total densities, and spin densities to be computed.
Several surfaces can be queued and computed at the same time, which is more efficientas the shell data only need to be computed once. The quality, colors and opacity ofthe surfaces can be be set in the dialog, and changes to these settings are saved to thepreferences as default values for any subsequent surfaces generated.
Note that each tick on the quality scale corresponds to roughly four times the numberof grid points as the previous one, and will therefore take roughly four times as long tocompute. Also note that densities require evaluating shell-pair values and are significantly
more expensive to compute that MOs, which only require shell values. For most purposes,
the pseudo densities (such as the pro-molecule density) available by double-clicking theSurfaces item provide an almost identical density surface, but at a much cheaper cost.
Once computed, the individual surfaces appear as sub-items in the MV, and these can befurther configured by double-clicking the item in the MV.
The Add Surface dialog also contains an interactive energy level diagram on the righthand side panel. The vertical scale can be zoomed in and out using the scroll wheel
16

Figure 4.3: An example molecular orbital for aniline.

(or equivalent) on the mouse. The scale can also be translated by a left click-and-drag.
Individual orbitals can be selected with a left-click and this will cause the energy of theselected orbital to appear below the diagram, as shown in Fig. 4.2.

4.3

Other Orbitals and Densities

Orbitals and densities other than those based on the canonical orbitals can also be visualized, if the appropriate calculation has been performed. Localized orbitals, naturaltransition orbitals, natural bond orbitals and attachment/detachment densities are all
possible.

Figure 4.4: Localised σ-bonding orbitals in ethene.

Dyson orbitals can be also displayed, either from the checkpoint file (requires Q-Chem17

≥ v5.1) or from the output file using an older version of Q-Chem. In the latter case, the
PRINT GENERAL BASIS option should be set to true in order for IQmol to generate
the basis function data.
$rem
METHODBASIS
EOM_EA_BETA
CC_DO_DYSONGUIPRINT_GENERAL_BASIS$end= EOM_CCSD= cc-pVTZ= [1,1,1,0]= TRUE= 2= TRUE! Ony required for QChem < v5.1

If available, Dyson Orbitals appear under the Surfaces item in the Model View and behave
similarly to the Canonical Orbitals.

4.4

Exporting Cube File Data

Each surface requires the generation of data on a 3D grid and these data are storedinternally so that subsequent calculations of the same surface (with different isovalues,for example) are much faster. To see what grid data is being stored, right-click on the
MO Surfaces item in the MV to bring up the context menu. Selecting the Show Grid Infomenu brings up the Grid Information dialog:

Figure 4.5: The Grid Information dialog shows what data is being stored and offers theoption of exporting the data in a cube file format.

From this dialog it is possible to export cube files containing the grid data. Right-clicking
on the desired grid brings up a context menu with the Export Cube File option. The
Switch Phase option swaps the sign of the data, and may be useful when comparing
different MOs where the (arbitrary) phases differ. Cube files can be saved for future18

plotting to avoid having to recompute the data, or for reading into another plotting
package. This is described in the following section.

4.5

Visualizing Cube File Data

Cube files contain volumetric data such as electron densities, molecular orbitals or electrostatic potentials (ESP). Because the data has been pre-computed, generating surfaces
using them is very quick. After opening a cube file, a Cube Data item will appear in theMV, double-clicking this item brings up the Add Surface dialog which allows an isovalue
surface to be requested (Fig. 4.6). The Signed check-box causes two isosurfaces to begenerated corresponding to ± the specified isovalue, and this should be checked for datasuch as MOs and spin densities.

Figure 4.6: The Add Surface dialog for cube data files

Cube file data can also be used to color a surface. This will require either two cube files(one containing the surface data and the second containing the property used to color thesurface) or a cube file and a checkpoint file. In either case, two files need to be loaded
into the one molecule using the File I Open Dir menu option, described above.
After creating a surface, double-clicking the surface item will cause the Configure Surface
dialog to appear:
The Property combo-box will have the cube data as one of the options, and selecting thiswill cause the surface to be colored according to the data in the cube file. The gradientcolors can be altered by clicking inside the gradient box.

19

Figure 4.7: Configure Surface dialog showing the Cube Data as an option for the property.

Figure 4.8: Electron density plot of aniline colored according to the electrostatic potential.

4.6

Orbital Animations

Orbital animations can be created e.g., for watching the evolution of the frontier orbitals
during a reaction, or watching the relaxation of an orbital during a self-consistent field(SCF) calculation. Key frames should be saved as cube file data in a single directory.
The base name of the cube files needs to match the directory name, e.g.
relaxation/relaxation.Frame.1.cube
relaxation/relaxation.Frame.2.cube
Open the directory, as described in Sec. 4, and right-click on the first cube file item whichappears in the MV, this will allow you to access the Surface Animator dialog:

20

The dialog is shown in Fig. 4.9 and can be used to adjust the key frame order and the settings for the animation, including the numberof interpolation
frames (higher values givesmoother results, but take longer). Pressing Calculate generates each frame and individual frames can be accessed in the MV. Use the playback options to run the animation,
recording it with the

button, if desired.

Figure 4.9: The surface animation dialog.

21

4.7

Potential Energy Surfaces

Several calculation types provide information about the potential energy surface (PES).
These include geometry optimizations, PES scans, reaction pathways and transition state
searches. On opening a Q-Chem output file containing one of these job types a Geometries
item appears in the MV and can be double-clicked to bring up the Geometries dialog.

Figure 4.10: The Geometries dialog showing a selected geometry.

The pathway can be animated in the viewer using the play button, and individual frames
can be selected from either the table or by selecting a dot on the energy plot. The Viewer
window will be automatically updated with the corresponding structure. As with the MOenergy diagrams, the energy plot can also be zoomed and scrolled.
Pathways can also be read in from an XYZ file. The format for this is simply a concatenation of regular XYZ files:
5-291.77177
Si 0.71979
H
0.73262H
1.10451
H -1.32980
H -1.22713
5-291.76961
Si 0.67218
H
0.73053
H
1.10815
H -1.29502
H -1.19409-0.08082-1.27801
0.936730.238940.18317-0.76577
0.229900.34825-0.282400.47002-0.07483-1.29428
0.952700.23567
0.17584-0.72935
0.22634
0.34714-0.33016
0.49156

22

....

The first line gives the number of atoms, the second (comment) line must have the corresponding energy and then the following lines have the xyz coordinates of each atom inthe system. This pattern repeats for as many frames as you have available.

4.8

Vibrational Frequencies

Vibrational frequencies can be read in from a Q-Chem output file and appear as aFrequencies item in the MV. Normal mode vectors can be visualized by selecting theassociated frequency item in the MV.

Figure 4.11: Arrows indicating a normal mode in acetic acid.

Double-clicking a frequency in the MV will cause the molecule to vibrate according to
the selected mode. Double-clicking the Frequencies item in the MV will bring up theVibrational Frequencies dialog:
This dialog contains an impulse spectrum showing the positions and relative intensities ofthe frequencies. Individual modes can be selected on the spectrum by clicking the hollow
circles, and this will also update the Viewer window with the selected mode. The impulses
can be broadened using either Gaussians or Lorentzians to give a more realistic looking
spectrum, and the image can be exported by righ-clicking on the spectrum to bring up acontext menu.
The horizontal scale of the spectrum can be zoomed (using the scroll wheel on the mouse)
and translated (left-click and drag) to give greater detail.

4.8.1

Isotopic Substitution

IQmol makes it easy to set up isotopic substitution loops for a Q-Chem vibrationalfrequency calculation. First select the atoms you wish to substitute in the Viewer and23

Figure 4.12: The Vibrational Frequencies dialog.

then select the Build I Set Isotopes menu option. A dialog will appear that allows you
to change the isotope for the selected atoms, along with the ability to change the pressure
and temperature used for calculating the thermochemical data.. Once the masses have

Figure 4.13: Isotope substitution dialog allowing non-default masses to be specified forthe selected atoms. Note that only one non-default mass is allowed for each element type.
been selected, click OK and an Isotopes item appears in the MV containing the substitution information. This can be edited or checked by double-clicking on the appropriate
Substitution item in the MV. Different substitutions can be set up by selecting different
atoms in the Viewer and/or different masses and these will be looped over in the Q-Chem
frequency calculation. Each substitution loop only requires a re-weighting of the Hessian
before diagonalization and so multiple loops can be obtained very quickly once the Hessian
has been calculated.
24

4.9

NMR Spectra

NMR shieldings and chemical shifts can be visualized using IQmol. A Q-Chem calculation (JOB TYPE = NMR) must first be run and the output file loaded into IQmol.
Shielding constants can be displayed as atom labels in the viewer by selecting the Display
I Atom Labels I NMR option.

Figure 4.14: Propanol showing the NMR shieldings

An NMR item will appear on in the MV and double-clicking this brings up the NMR
Spectrum dialog.

Figure 4.15: The NMR Spectrum dialog

Various spectra can be displayed (1 H, 13 C etc.) and if a reference is available (e.g. TMS) for
the given level of theory, chemical shifts can also be displayed. Selecting rows in the table
will also highlight the corresponding atoms in the viewer for easy identification. If impulses

25are plotted, then these can also be selected to identify which atoms are responsible forthe selected signal.
As with the vibrational frequencies spectra, NMR spectra can be exported by rightclicking on the image to bring up the context menu.

26

5Appearance

5.1

Viewer Camera

When viewing a molecule, the coordinates of the atoms and surfaces are fixed in the worldreference frame. When manipulating the molecule using the mouse (see § 2.5) it is actually
the camera that changes position and orientation. If more precise control of the camera isrequired (e.g. for reproducing a particular image), then the camera can be configured via
the Display I Camera menu option. By default the camera uses a perspective projection
to give better depth perception, but an orthographic projection can also be selected, if
required.

Figure 5.1: The Configure Camera dialog allows precise control over position and orientation of the Viewer camera.

5.2

Clipping Planes

A Global clipping plane is available and can be viewed by enabling the Clipping Plane item
under Global in the MV. Clipping planes can be useful for showing surface information
without obscuring the underlying molecular structure. While selected, the clipping plane
can be arbitrarily translated and rotated using the manipulate selection mode describedin § 2.5. Alternatively, a precise orientation of the clipping plane can be set by doubleclicking the Clipping Plane item in the MV to open the configuration dialog.
Individual surfaces can be clipped, by selecting the Clip check-box in in the surface configurator dialog (see Fig. 4.7) for each surface. Once a surface has been clipped, the clipping
plane remains active even if hidden by unchecking the Clipping Plane item check-box in
the MV.

27

Figure 5.2: Image showing a clipped density surface allowing the underlying molecular
structure to be more easily seen.

5.3

Shaders

IQmol uses GLSL shaders to enhance the appearance of the Viewer. Shaders can be madeactive by going to the Display I Appearance menu option, and selecting the Shader panel.

Figure 5.3: Configuring shader options

Changes to the various
 shader options
 will take immediate effect and can be saved as thedefault by clicking Save As Default . Some shaders have several lights associated withthem, and these can be turned on or off to give different lighting effects.
28

5.4

Exporting POV-Ray Files

IQmol can generate scene files (.pov) for use with an external POV-Ray package. POVRay uses ray tracing algorithms to generate images of very high quality and resolution. POV-Ray is open source, but several official and unofficial packages exist withpre-compiled binaries, in particular testing was done using the Mac version of MegaPOV.
Binaries can be downloaded from the following sites:• Mac: http://megapov.inetart.net/povrayunofficial_mac• Windows: http://www.povray.org/download
To export a scene file use the Edit I Generate POV-Ray Input menu option. Scene files
are plain text with instructions for the POV-Ray program, and can be edited by the user
for fine tuning before processing. IQmol can apply several effects (e.g. surface textures)to the molecule, and these can be configured by opening the Appearance dialog (Display
I Appearance menu option).

Figure 5.4: Appearance dialog showing Options for POV-Ray scene files.

Note that IQmol attempts to generate a scene file that matches as closely as possible
what is shown in the Viewer window. However, there may be slight variations in lighting,color and camera angle, and some trial and error may be required to obtain a pleasing
image. In particular, the Gamma value may need experimenting with, and changes to theambient and diffuse lighting (via the shader options) may be required depending on if alight or dark background is selected for the scene.

29

6

Sample Images

Here are some examples of the kinds of images that can be generated using IQmol’s
shader and rendering options. Most of these are simply screen shots, however, tiles 3, 4
and 5 in the middle row were obtained by running POV-Ray on the IQmol generated
scene file.

30

References[1] Pilati, T., and Forni, A. Symmol: a program to find the maximum symmetrygroup in an atom cluster given a prefixed tolerance. J. Appl. Cryst. 31 (1998), 503.
[2] Rappe, A. K., Casewit, C. J., Colwell, K. S., Goddard III, W. A., andSkiff, W. M. Uff, a full periodic table force field for molecular mechanics andmolecular dynamics simulations. J. Am. Chem. Soc. 114 (1992), 10024.

31