Follow the first link above for the latest version of MOPlot |
  MOPlot
program for Plotting Molecular Orbitals (and other stuff)
Contents 1. Introduction and a bit of history2. Program overview 2.1. Features of MOPlot 2.2. Rules for drawing MO's 2.3. Limitations and caveats 3. MOPlot program interface 3.1. Overview 3.2. Camera tab 3.3. View tab 3.4. MO tab 3.5. Vibr tab 3.6. Geom tab 3.7. Extra tab 3.8. Misc tab 4. Gathering postscript pictures for output 5. Help options 6. Referencing
Introduction and a bit of history Although many programs are available that can produce very pretty and quantitatively accurate electron density plots from the outputs of quantum chemical calculations, the computation of each such plot is rather time-consuming, even on current hardware, so viewing several molecular orbitals, e.g. to choose a set to include in an active space, may require a considerable degree of patience. In addition, the resulting images often do not immediately reveal the "nature" of a molecular orbital very clearly, and images must be turned and twisted around before one gains insight, e.g. into the nodal structure (note that people who use such plots in publications usually have to display several views to make them intelligible to the readers). In our work we felt the need for a program that would allow us to flip rapidly through a set of MOs, such as to have an immediate view of their nodal structure and an approximate idea of their spatial distribution over a molecule. The graphic tools implemented in GaussView, Spartan, CaChe, Hyperchem etc. visualize MO's as three-dimensional electron density plots which take some time to calculate and do not always readily yield insight into the nodal structure of MO's. Also, we were tired of trying to adapt the voluminous pixel images produced by the above electron density display programs for use in publications, so we wanted a utility that produces vector output. HOMO of admantylidene displayed by MOPLOT fills this gap by producing 2-D projected images of geometries, MO's and normal vibrations produced immediately from the output of various standard quantum chemical programs. An solution that fulfils all of the above requirements had actually been proposed and implemented on a PDP 8 computer over 30 years ago by Haselbach and Schmelzer with the purpose of displaying MOs from semiempirical calculations in the ZDO approximation. In the 1980s we adapted a copy of this Fortran program to PCs running DOS, and later to the Windows operating system. Since the 1990s we used it also to display MOs from ab initio calculations, although we had to cut some corners to do that (see below). However, the advent of Windows NT/2000 (which lacks a DOS kernel), tolled the death bell to MOPlot, so it began to fall into disuse, and students went back to make guesses about the shapes of orbitals instead of viewing them. Then we had the good fortune to gain the collaboration of Dr. Rouslan Olkhov who, while waiting for parts for a laser system to arrive, agreed to look into the MOPlot problem. After thinking about this for a moment, Rouslan proposed to rewrite MOPlot from scratch, using the LabView programming environment. At first we were a bit incredulous, but it was not long before Rouslan produced the first MOPlot-style pictures of MOs that looked much prettier than anything we had had before. In addition it turned out that the user interface was much easier to program than with any version of Fortran, so we decided to scrap all the code we had and build on what Rouslan was producing. After about six months, he had the first full working version of MOPlot, and ever since he has been ironing out bugs, adding features (see below), and refining the user interface until he arrived at the present version which is not only much better than anything we had before, but has the advantage that it runs on any platform and operating system for which LabView exists, i.e. Windows, Linux, MacOS. As our University holds a site licence to the LabView professional development system this makes it possible for us to offer compiled versions of MOPlot for all these platforms which run independently of the LabView programming environment. Overview of MOPlot program
MOPlot can display In addition, MOplot can MOPlot has a help mode where all the features in the different menus are explained in pop-up windows when one moves with the mouse over the feature on which one would like to get information. Rules for drawing MO's An LCAO-representation of MO's is given in the ZDO approximation (no overlap densities). AO contributions are displayed as single spheres (pure s-AO's), or pairs of contiguous spheres (pure p-AO's) whose (combined) area is proportional to the contribution of an AO to the electron density at the atom on which this AO is centered In s/p hybrid MO's, the area of one sphere is enlarged relative to the other in proportion to the size and sign of the s-contribution. If c(s)>c(p), a single sphere is drawn whose center is displaced in the direction of the (px, py, pz) vector in proportion to the p-contribution. (in MOPlot, the s/p ratio threshold where the display switches from pairs of spheres to single spheres can be changed. The default value is 1) An illustration of this on a progressively pyramidalized methyl radical is shown below: Limitations and caveats It is important to remember that MOPlot was originally developed for viewing the results of semiempirical calculations to which the ZDO approximation is inherent. Its MO display mode reflects these theoretical models and MOPlot should therefore not be regarded as a program that yields accurate renditions of ab-initio wavefunctions. If this is what one wants, one should resort to one of the many fancy commercial or shareware electron density plotting packages (e.g. Spartan, GaussView, Molekel ...) which do an admirable job at that. The main purpose of MOPlot is to gain quick insight into the nodal properties of MOs and their approximate distribution over a molecule. These features can be more easily recognized in a ZDO-type rendition, such as it is provided by MOPlot. MOPlot knows only about AO-coefficients, it ignores that fact that the same AOs have different exponents in different atoms. As long as one stays in the CHON group of atoms and within the realm of valence MOs, the errors introduced by this are usually quite tolerable and certainly not misleading, but if one does e.g. hydrogen fluoride, the results have very little relation to reality with regard to the electron density distribution. Also, to keep things simple, the contributions of inner and outer split valence basis AO's are simply summed up and polarization functions are ignored. After this, MOs are renormalized to one, ignoring overlap, to make them comparable amongst each other. MOPLOT is not a commercial piece of software, and we can assume no responsibility whatsoever with regard to its performance. We have tested it to the extent that we can, and it looks as though it runs quite reliable. Nevertheless users are very likely to encounter bugs here and there and we would be very indebted if users would take the time to communicate these bugs to us. MOPlot program interface When you first open the program you are presented with a stark gray square on the left hand side of the program window (the picture area), and with a collection of menus on the right hand side. Before you can do anything, you must read a file by clicking on the blue "open file" button whereupon you are presented with a standard file selection window. MOPlot can only open Gaussian outputs, Molden inputs, and files in its proprietary standard format. If you try for example to open a Word document, MOPlot will complain. Gaussian outputs can assume gargantuan sizes. Since MOPlot must plough through the entire output until it can be sure that is has found everything it needs, opening a file can take a little while. Don't despair, eventually it will get there, and it will display in the area below the menu bar what it found. Once a file has been successfully read, the screen should look something like this: The buttons in the menu bar allows the user to [open (a new) file], to prepare postscript [output] of pictures observed in picture area (see instructions below), to get some information [about] the authors of the program, and, finally, to [exit] the program and free the precious memory of the user’s computer. The picture area is simply the graphics display and has no active interface elements., apart from the fact that you can twist and turn the molecule in all (hopefully intuitive) directions by dragging the mouse over it while you hold down the(left) mouse button. Pressing the shift key while you drag the mouse allows you to rotate the molcule in the plane of the display. Note the little coordinate system on the left bottom: it indicates the axes with regard to which the cartesians of the molecule were defined in the output. The tab menu contains all the control elements, which allow user to manipulate the data to be shown in the picture area. Six pages are accessible in the tab menu:   |
  [camera] is provided for more precise control over the rotation of the molecule around three perpendicular axes. A slide allows an adjustment of the [rotation speed]. At any moment the set of Eulerian angles, which correspond to the current orientation of the molecule in the picture area, can be memorized via [get angles from picture] into the a, b and f fields. The [apply angles to picture] button allows to instantly change the orientation according to the rotational angles entered in the a, b, and f fields (this feature can be used to re-establish a particular orientation). Finally there is a [scale] slide to increase or decrease the size of the molecule, a [light gradient] slide, which adjusts the depth hue, and [3D shading] button to invoke pseudo-3D rendering of the picture.   Note that the shaded spheres look much nicer in the postscript output than they look on the screen where the number of hues is kept small to maintain the quickness that is the prime quality of MOPlot.   |
  [view] allows to select what one wants to view (except for MOs and vibrations, those are controlled in their own menus). Atoms are represented as circles whose radii are proportional (including the factors from the scale slides) to:
Note however that, once created, dummies cannot be removed! But, of course, one can always reload the original file.   |
  [mo] controls the display of MOs:
  |
  [vibr] controls the display of normal modes (vibrations)
 
"spherical vectors" after: W. Hug, Visualizing Raman and Raman
optical activity generation in polyatomic molecules. Chemical Physics (2001),
264(1), pp. 53-69   Note: one has to leave the loop mode before doing anything else.   Note: If the frequencies of mixed vibrations are different, the program treats them as if the ratio of their frequencies were integer from 1 to 6 (six is a maximum since the cartoon is limited by 20 frames and we need at least 3 frames per vibrational period).   |
  [geom] is a utility to calculate geomety parameters (bondlengths and –angles, dihedrals) between 1-4 points in the molecule that is currently displayed. The unique feature of MOPlot is, that these points may but need not coincide with atoms, they can also be midpoints between atoms (in which case, instead of a single number, the input to the boxes on the right consists of two numbers separated by a "-"). This has often proven to be very useful to calculate, e.g. the angle between the two benzene rings in 9,10 dihydroanthracene. When entering numbers, one may toggle to the next field with Tab. However, after entering the last number, one must click with the mouse somewhere outside the boxes to activate the calculation (just hitting "Enter" at the end will not do!) [List feasible bonds] button creates a list of atoms that appear to be linked by chemical bonds along with bond lengths. This list can be saved to an ASCII file or copied into the clipboard.   |
  [extra] The first extra feature is the option to view the course of geometry optimisations, similar to Molden. Initially an autoscaled plot of all the steps is displayed, with the relative energy set to zero for the first point. With the magnifying glass one can zoom in on different parts of this plot with various options (or undo a zoom operation). With the hand one can move the plot around, and with the cursor tool one can grab the cursor and place it on any desired point (initially it is on the last point which is where one has to fetch it first), whereupon the geometry of this point is displayed (and can be saved to disk by save step geometry. One can also move the cursor incrementally to the left with the green diamond and to the right with the red diamond at the top right of the plot. At any time, for example if one is lost, one can autoscale the plot again. The second option offers to [remap atomic labels], which overrides numbering of atoms from input data file according to the custom table which can be edited via [edit label map]. The rest of the controls deal with plotting "gradient difference" and "derivative coupling" vectors at a conical intersection point if such a calculation was performed. Initially, the vectors as they appear in the output are displayed. Upon pressing the [orthonormalize] these are normalized to 1Å and orthogonalized. The vectors can be displayed a [single] one at a time or [both] together. Changing the [phase] angle results in a mixing of the vectors to produce two new vectors in the same plane (the orientation of a [single] vector at zero phase corresponds to the original orientation of the derivative coupling vector). [copy vectors to vibr] converts the current vectors into a pair of vibrational modes. This allows to use all commands from [vibr] submenu on these vectors, including animation. [save discoid coords] creates a file with a set of geometries that lie on a circle in the plane of the derivative coupling and the gradient difference vectors, centered at the conical intersection. The radiu of the cirle and angular increment at which geometries are evaluated should be specified. On pressing the "write file" menu, a file selection box pops up.   |
 
[misc] contains miscellaneous controls of the way things are displayed in MOPlot, such as the fonts used to label atoms or to write the information on vibrations or MOs at the top of the display area. Also, if you don't like the colors of MOs, charges, spins or whatever, the colors of just about everything can be changed in this menu. The drop-down menu [picture colour style] allows you to switch between color, black-and-white, and a black/ white/grey picture style, which are provided for the purposes of creating publication output for journals which do not accept color figures. The [perspective] slide adds varying degrees of perspective with regard to the axis perpendicular to the screen. The [external ps viewer] option lets you specify the location of a postscript viewer program (such as ghostview) on your computer (click on the litte folder to the left). The[write ini file upon exit] button allows you to save the present program parameters in the moplot.ini file. The saved parameters will be taken as default values next time you open the MOPlot program. One may also edit the moplot.ini file (which is in plain text format) manually. If the .ini file becomes corrupted simply delete it, and restart MOPlot which will return to its original default values and create new ini file with those parameters if you activate the [write ini file upon exit] option. Currently supported program parameters are:
  Note for Linux users: normally, the MOPlot ini file can only be updated by its original owner. If other users want to use their own ini files they should use private copies of MOPlot, only two files are required: moplot and lvanlys.so   |
  [output] If you click that button in the menu bar at the top, MOPlot kicks into a mode that lets you prepare postscript output. The label of the [output] button changes to [gathering] and a new menu appears below: The first thing you must decide is, how many pictures you want on a page (default: one, you can go up to twelve by clicking on the bar below "page layout") Then, every time you have a picture in the display area that you want to add to your page, click the blue bar on the right side and a little field with a label will pop up in the schematic page layout. At and time you can delete one of the cells by entering the cell number in the blue field and press [delete cell]. By default the first cell (on the top left) is deleted. You can also delete everything by [clear page]. If you add more pictures than fit on your layout, it will tell you how many pictures do not fit. As you delete cells, these "surplus" pictures will be added one by one to the layout. The picture area remains active as long as you do not "clear" it, or exit the program. Thus you can also read another molecule and add views of that to the current output page. When your page is full, or if you decide that you have everything you need, you should write your page to disk. Before you do that you must decide whether you want interpreted postscript [ps] (which is what you spool to your printer) or encapsulated postscript [eps] (which is often used to import into graphics programs). When you click on [save page] you are presented with a standard file dialog box that lets you decide how you want to name the file and where you want to save it. There is an option to [preview (the) page] with an external postscript viewer program, which must be specified on the [misc] menu. Clicking on the [gathering] button while you are in the output mode gets you out of this mode, i.e. the output menu and the blue bar on the side of the picture area disappear. However, the output page is maintained.   Example of the MOPlot postscript output with 3D shading turned on   |
  [help] clicking this button puts MOPlot into help mode which means that each time you move your mouse over a feature in a menu or an area where some information is displayed, some information on that feature pops up in a little window. If you are tired of this, click [help] again, and the help menus disappear. [manual] starts the defauld web browser and loads the online manual of MOPlot. Naturally the feature does not work if the host computer is not connected to the internet. Referencing Program reference to be used in publications is: Rouslan V. Olkhov, Thomas Bally, MOPlot v.1.86, http://www-chem.unifr.ch/tb/moplot/moplot.html Thank you for reading this far! |