Cerius2 Modeling Environment



3       Tutorials

This chapter contains several tutorials covering the basics of working within the Cerius2 modeling environment. The tutorials in this chapter are:

Basic GUI Functions

Using the 3D Sketcher

Managing Models and More

Managing and Modifying Graphs

Basic Molecular Simulation: Minimization & Dynamics


Basic GUI Functions

This section contains an introductory Cerius2 Visualizer tutorial. It shows you the basics of the GUI including loading, viewing, and display of simple models; manipulating control panels; and exiting the program. You learn about functions on the menu bar and the tool bar.

This chapter is divided into three parts:

Before you begin

To complete this tutorial, you need a licensed copy of Cerius2 and a directory into which you can save the files you create during the tutorial.

In this tutorial, you familiarize yourself with the look and feel of the Cerius2 GUI and Visualizer. You use the Cerius2 menu bar (File, Edit, Build, etc.) and the tool bar immediately below it. You learn how to load a model and then practice viewing, displaying, and selecting it. The tutorial concludes with a brief exercise about control panels.

Summary of tutorial

The tutorial is divided into four short exercises:

A. Loading and displaying a molecular model

B. Viewing and labeling a model

C. Selecting atoms

D. Manipulating control panels

Step-by-step tutorial

This section takes you through a series of short progressive exercises designed to familiarize you with some basics of the Cerius2 Visualizer.

A. Loading and displaying a molecular model

In the following steps, you load an msi format structure file of trimethylbenzene and display it in both Stick and Ball styles.

1.   Start Cerius2 and load the model.

If Cerius2 isn't already open, type

>	cerius2 
in a UNIX shell to run the program.

The program takes a few seconds to initialize; a welcome screen appears briefly and then the Cerius2 GUI appears.

Troubleshooting

If typing cerius2 does not start the program, consult the Cerius2 Installation and Administration Guide for help on setting up the correct environment for running Cerius2.  

Select File/Load Model... to open the Load Model control panel.

Remember to use the left mouse button for choosing menu items.

Load the model 135benz.msi from the directory Cerius2-Models/organics: double-click the Cerius2-Models/ directory name in the browser box to open that directory; double-click the organics/ directory to open the organics directory; then double-click the 135benz.msi filename to load the file.

Note

You may need to use the scroll bar at the side of the file list to locate the organics/ directory.  

Your model window should look similar to this:

Troubleshooting

If you need to return to the Cerius2-Models/ directory, choose the Run directory selector at the top right of the browser box.     If you can't find the file in the organics/ directory, make sure that the File Format popup reads MSI and that the text box at the foot of the panel reads *.msi.  

2.   Change the display style of the model.

Use the Display Style popup menu on the left end of the tool bar to change the display style to Ball.

The atoms are displayed as spheres. By default they are filled to 0.7 of their van der Waals radii.

Troubleshooting

If some of the atoms stayed in stick display? If not all the atoms appear in Ball display?     You may have inadvertently selected some of the atoms in the model window. If any atoms were in a selected state, then only those selected atoms will appear as balls. To correct this, click in some empty space in the model window: this deselects any previously selected atoms. Then reselect Ball from the Display Style popup.  

Select Stick from the Display Style popup menu.

The model returns to stick display where all bonds drawn as lines.

B. Viewing and labeling a model

The following steps show you how to rotate and translate the trimethylbenzene model. You also learn how to label atoms in the model by element name.

1.   Manipulate the model's presentation using the mouse.

Rotate the model by placing the cursor over the model window and then dragging while holding down the right mouse button. Familiarize yourself with the action of the mouse.

If the cursor is near the center of the screen, it rotates the model about an axis in the plane of the screen and perpendicular to the direction in which you drag. If the cursor is towards the edge of the screen, it rotates the model around an axis perpendicular to the screen.

Translate the model by dragging over the model window with the middle mouse button depressed. Scale the model by dragging over the model window with the middle mouse button and the <Shift> key depressed, or with both the middle and right mouse buttons depressed. Try rotating the model in some other display styles. Then return the display style to Stick.

Note that rotation is quickest in the stick display.

[Aside: The Ellipsoid style is for models for which temperature factors have been defined. The Polyhedra style is used to display coordinated cations, for example, in zeolites.]

2.   Label the atoms in the model with their element type.

Change the Atom Labeling popup on the tool bar from No Label to Elements.

All the atoms in the model window should now be labeled as hydrogens or carbons, either C or H.

If the display style isn't stick, the labels may be difficult or impossible to see.

Troubleshooting

If some of the atoms are not labeled, you may have inadvertently selected some of the atoms in the model window. If any atoms are selected, then only those selected will be labeled. To correct this, click the Deselect All icon in the tool bar, and then reselect Elements from the Display Style popup.  

C. Selecting atoms

In the following steps, you learn several techniques for selecting and deselecting atoms.

1.   Select an atom.

Use the left mouse button to pick one of the atoms in the model.

In stick display, when the atom is selected it is highlighted. Cerius2 beeps to let you know that you've selected the atom. The atom is also identified in the text window.

Select several more atoms by holding down the <Shift> key as you select them.

You should now have three or four atoms highlighted in the model window.

2.   Change the display style.

Use the Display Style popup to change to Ball & Stick.

Note that the new display style only applies to the selected atoms. In the ball and the cylinder display styles, selected atoms are indicated with thin stripes.

3.   Select and deselect all atoms.

Deselect all atoms by clicking on empty space in the model window.

Find the Select All and Deselect All icons on the tool bar. They look like the icons in the margin.

Next, you use these buttons to select and deselect atoms in the model window.

Choose the Select All button then changethe Display Style popup to Ball.

All the atoms in the window are highlighted and then appear as balls wrapped in thin blue lines.

Choose the Deselect All button.

All the atoms are deselected -- the thin blue lines disappear.

D. Manipulating control panels

In the following steps, you learn about control panels: how to open, close, scale, and stack them. You use on-line help to get information about the controls on the panels. And finally you learn how to re-initialize and/or exit from Cerius2.

1.   Close and open control panels.

Close the Load Model control panel by clicking on the close box in the right corner.

All control panels can be closed this way.

Open the Save Model control panel by choosing File/Save Model... .

Notice that those pulldown menu items that open control panels all have names that end in ellipses (...).

2.   Open more control panels.

Choose Build/Element Defaults to open the Edit Elements control panel.

You would use the Edit Elements control panel if you wanted to alter any element properties such as the display color used in Cerius2.

Choose View/Display Attributes... to open the Display Attributes control panel.

The Display Attributes control panel has the same functionality as the Display Style and Atom Labeling tool bar popups, which you have already used.

You can select the Preferences button to change the defaults for the different displays. For example, you can change the size of the spheres in the Ball display from 0.7 to 1.0 in van der Waals radii of the atoms.

3.   Access help for a command in a control panel.

Click the Preferences... button on the Display Attributes control panel to open the Style Preferences control panel.

Again, control panel buttons with names ending in ellipses (...) open other control panels.

Click with the right mouse button on a command in one of the control panels.

This brings up a help text panel that describes the function.

Click in the help text panel to close it.

Try clicking with the right mouse button in other items in the GUI: check boxes, buttons, popups, sliders, menu items, and so on.

4.   Access the Stereo control panel from a pullright menu.

Select View/Graphics> and then move the mouse to the right to select the Stereo... item on the pullright menu.

The Stereo control panel opens. You would use the Stereo view control panel if you wanted to view the model as a stereo pair or with the Crystaleyes stereo hardware.

Right arrows at the end of pulldown items indicate that there are pullright menus for these items.

5.   Move and resize control panels.

Click and drag some control panel title bars to different screen locations.

Resize some control panels by clicking and dragging their edges or corners. Click the Cascade Panels button on the tool bar.

This stacks all the open control panels with their title bars visible.

Bring the last panel to the front by clicking on its title bar, and close it by selecting its close box (in the top right corner).

6.   Close the control panels.

Select the Clear Panels button on the tool bar.

The remaining control panels close.

7.   If you are going to continue immediately with another tutorial:

Re-initialize the program by selecting File/New Session and selecting the Confirm button on the Re-initialize dialog box.

Re-initialization resets all the Cerius2 control panels to their original defaults.

Alternatively, if you want to stop now:

Choose File/Exit to quit from Cerius2, and select the Exit button on the Exit Option... dialog box.

Review

In this tutorial you learned how to:

You may now want to proceed to the next tutorial, "Using the 3D Sketcher", where you learn how to build your own models.

Alternatively, you may want to learn more about the Visualizer by reading the user manual sections listed below.

Related material

The Visualizer is discussed in detail in the Cerius2 Modeling Environment. The sections listed here will provide you with more detailed information about atom display and selection than are included in this tutorial book.

1.   Sections in the Introducing Cerius2 and Viewing and Displaying Models chapters of the Cerius2 Modeling Environment.

2.   The section, Selecting atoms and groups of atoms in the Building Models chapter of the Cerius2 Modeling Environment

3.   The Managing Models and Sessions chapter of the Cerius2 Modeling Environment.


Using the 3D Sketcher

This chapter contains a tutorial on molecular sketching. It introduces the Cerius2 3D Sketcher. You learn how to draw and edit models in 3D and how to use the Clean function to do simple optimization of model geometry. You also learn how to save the models, which you build, into structure files.

This chapter is divided into three parts:

Before you begin

To complete this tutorial, you will need a licensed copy of Cerius2 and a directory into which you can save the files that you create during the tutorial.

Before you begin the tutorial on page 50, you should be familiar with the contents of the section Basic GUI Functions. If you have not read this section, you should ensure that you know how to do the review items listed in the Review on page 61 before beginning this tutorial.

In Basic GUI Functions, you loaded a model from file into the model window. The purpose of this tutorial chapter is to show you how you can build your own models using the 3D molecular Sketcher. The Sketcher is the principal building tool of the Cerius2 Visualizer.

Note

There are other ways to build models in Cerius2: using the Analog, Crystal, Surface, Interface, Polymer, and Amorphous builders. You may have licensed some or all of these modules with your configuration of Cerius2. The tutorials in the Cerius2 Tutorials: Materials Science and Life Science demonstrate the use of some of these sketchers.  

Below is an illustration of the Sketcher control panel. The Sketcher provides tools for mouse-driven sketching, editing, and cleaning of models. The operation of the Sketcher is modal, that is, the function of the left mouse button in the model window is determined by which icon is currently selected in the Sketcher control panel. The default setting is to have the Selection mode icon highlighted. This allows normal Cerius2 operation, with the mouse button selecting atoms (just as in the previous tutorial). It is a good habit to always highlight the Selection mode icon before closing the Sketcher control panel.

Summary of tutorial

The tutorial is divided into five short exercises:

A. Building cyclohexane and saving the structure to file

B. Converting cyclohexane into the boat conformation

C. Editing cyclohexane to make benzene

D. Editing benzene to make phenol

E. Building xylene from template fragments

Step-by-step tutorial

This section takes you through a series of short progressive exercises designed to familiarize you with some basics of the Sketcher.

A. Building cyclohexane and saving the structure to file

In the following steps, you open the Sketcher panel and use its Sketch, H Adjust, and Clean tools to build a model of cyclohexane. After you build cyclohexane, you save it into a structure file.

1.   Start Cerius2 and open the Sketcher.

If Cerius2 isn't already open, type

>	cerius2 
in a UNIX shell to run the program.

The program takes a few seconds to initialize; a welcome screen appears briefly at first, and then the Cerius2 GUI appears.

Open the Sketcher control panel by selecting the Build/3D-Sketcher... item.

The Sketcher control panel appears.

2.   Sketch the cyclohexane model.

Choose the Sketch icon: x

Note the Sketch With element is set to C for carbon. Notice also that in the model window the mouse cursor appears as an inverted V; this indicates that the mouse is now in sketching mode.

Click once in the model window.

A single carbon atom appears, represented as a cross, and a dotted "guideline" follows the mouse cursor. Keep the cursor within the model window -- if the cursor moves out of the model window, you lose the guideline.

Troubleshooting

If you lose the guideline, click once on the first carbon atom (with the mouse still acting in Sketch mode), and the guide line will reappear. (If you click directly on the existing carbon, you won't generate a second atom.)  

Click somewhere else to put down a second atom and continue to build a ring of six carbon atoms --terminate the ring by clicking on the original atom.

You should now have a hexagonal form somewhat similar to this in the model window:

You can correct any mistakes with the Undo button.

3.   Set the sketching mode back to Selection.

Highlight the Selection mode icon once more so that all future clicks in the Model window do not result in new atoms!

4.   Adjust the hydrogens and minimize the structure.

Click the H Adjust button.

The H Adjust command adds (or deletes) hydrogens to all atoms in order to satisfy their valences. In this case, H Adjust adds two hydrogen atoms to each of the carbon atoms.

Hold down the Clean button. Release the button when the atoms stop moving.

Clean attempts to adjust bond lengths and angles to reasonable values. Note that recleaning after manipulating or deforming the structure may find a different structure.

You have built a cyclohexane ring. The cyclohexane ring that you see on the screen is probably in the "chair" conformation.

5.   Experiment with the display and view of the model.

Try altering display styles and rotating and rescaling the model to view it. (If you need hints on how to do this, refer to the Basic GUI Functions tutorial lesson.)

Reset the view -- use the icon beside the label popup on the tool bar.

6.   Save the cyclohexane model into a structure file.

Select File/Save Model... to open the Save Model control panel. Type cyclohexane.msi into the filename text box, which is next to the Save button, then select the Save button to save the structure to a file. Use the directory selector popup and the file browser on the panel to find the correct directory before you save the file.

B. Converting cyclohexane into the boat conformation

In the following steps, you learn how to move atoms to convert the chair conformer of cyclohexane into the boat conformer.

1.   Adjust your view of the model.

Use the mouse to rotate the model to view the "chair" from the side.

Select the carbon on one end of the ring as you look at it, and select the two hydrogens attached to that carbon.

(Remember you need to hold down the <Shift> key to select multiple atoms.)

Hold down the <Control> key and drag with the middle mouse button to move the selected atoms. Drag these atoms well to the other side of the plane defined by the four carbons in the seat of the chair.

Clean the model.

You should now have either the boat form or the twist-boat form of cyclohexane. The boat form is shown above.

Note

If cleaning returns the model back to the chair form, drag the CH2 group even farther above the seat of the chair than it was in step 3 and even closer to the back of the chair. Then repeat the Clean step.  

Rotate and view the model.

C. Editing cyclohexane to make benzene

In the following steps, you learn how to edit bonds and to edit and delete atoms, and so build benzene from cyclohexane.

1.   Delete the hydrogens from cyclohexane.

Find the two popups on the tool bar that currently read Atom. Change the second of these to El for element.

Both of these popups control the selection rules for atom picking. Setting the second popup to El means that all atoms of the same element type are selected when one atom of that type is picked.

Select one hydrogen atom.

All hydrogens appear selected in the model window.

Double-click the Delete Atom icon in the Sketcher control panel.

All the hydrogens in the model are deleted.

This delete procedure has illustrated the alternative way of applying the Sketcher icons -- double clicking on an icon carries out its action on all selected atoms.

2.   Make the double bonds.

Change the second selection rule popup back to its default Atom setting. Also, highlight the Selection mode icon once more.

Change the Edit Bond popup from Single to Double and click any carbon atom. Pick an adjoining carbon atom.

A dotted guideline appears and then becomes the double bond between the carbon atoms.

Repeat to make alternate bonds in the ring double.

3.   Adjust the hydrogens and minimize the model.

Use the H Adjust and Clean buttons, then rotate and view this planar benzene.

D. Editing benzene to make phenol

In the following steps, you add an OH group to benzene to make phenol. Also, you learn how to use the Sketcher tools to alter torsion angles, bond angles, and bond lengths.

1.   Add an oxygen.

Use the popup adjacent to the Edit Element text box to select O for oxygen.

The popup is a fast way to choose commonly used element symbols.

Select the periodic table icon after this popup.

A periodic table appears.

Select any element in the table.

The element is written to the text box.

You can also change the element by typing a new element symbol into the Periodic Table text box.

Close the periodic table. If necessary, restore the Edit Element symbol to O.

Pick any one of the hydrogens to change it to an oxygen.

In Edit Element mode, the mouse cursor is a circle.

Use the H Adjust and Clean buttons. Highlight the Selection mode icon once more. Rotate and view the model; try some other display styles, then return to the Stick display. Choose the Twist icon.

In the model window, the mouse cursor for twist mode is shaped like a bow tie.

You must select four atoms to define a twist torsion.

Pick the hydrogen attached to the oxygen, the oxygen itself, the carbon attached to the oxygen, and one of the adjoining ring carbons.

In this case, don't hold the <Shift> key down; the atoms don't become highlighted when they are selected in Twist mode.

Drag the mouse over the model window with the left mouse button depressed to increase/decrease the torsion.

The change in the torsion angle is interactively updated. Note that the atom picked last is the one that moves.

2.   Try the other geometry manipulation tools.

Use the Stretch and Angle tools, which work similarly to Twist.

3.   Save the new structure.

Save the phenol structure into a file called phenol.msi. If necessary, see Step 6 on page 68 for a reminder of how to save.

Clear the model window.

E. Building xylene from template fragments

In the following steps, you use the template tool to create a model of p-xylene. The template tools let you sketch with predefined molecular fragments (called templates) -- for common structures, this is often more convenient than sketching atom-by-atom.

1.   Create a model from the phenyl template.

Click the Template icon on the Sketcher -- this looks like a "T" and is located about halfway down the box.

The current template is phenyl.

Click in the model window.

A phenyl ring appears.

2.   Add a methyl fragment from a template to make toluene.

Click the Templates... button in the Sketcher. Double-click the organic/ directory and find the template file methyl. Double-click methyl.

A file browser opens on a directory called templates/ and the template name in the Sketcher changes from phenol to methyl, indicating that the methyl template has been successfully loaded.

Pick one of the phenyl hydrogens.

The hydrogen atom is replaced by a methyl group to make toluene.

3.   Add another methyl group.

Select another of the phenyl hydrogens -- choose the one farthest away from the methyl group to make p-xylene.

4.   Save the model.

Save the p-xylene structure into a file called xylene.msi. If necessary, see Step 6 on page 68 for a reminder of how to save.

5.   If you are going to continue immediately with another tutorial:

Re-initialize the program by selecting File/New Session and selecting the Confirm button on the Re-initialize dialog box.

Re-initialization resets all the Cerius2 control panels to their original defaults.

Alternatively, if you want to stop now:

Choose File/Exit to quit from Cerius2 and select the Exit button on the Exit Option... dialog box.

Review

In this tutorial you learned how to ...

o Sketch atoms

o Move and delete atoms

o Edit atom and bond types

o Edit bond angles and lengths

o Adjust hydrogens in a model

o Use the template tool

o Clean a model

o Save a model to file

You may now want to learn more about the Sketcher by:

Or you may proceed immediately to the next tutorial, Managing Models and More.


Managing Models and More

This chapter contains a progressive tutorial designed to show you some more functions of the Visualizer including managing multiple models, copying between model spaces, measuring models, and printing. You also get additional practice at loading, saving, and sketching models.

This chapter is divided into three parts:

Before you begin

To complete this tutorial, you need a licensed copy of Cerius2 and a directory into which you can save the files that you create during the tutorial.

Before you start this tutorial lesson (page 50), you should have read (or be familiar with the contents of) Basic GUI Functions, and have done the tutorial exercises in Using the 3D Sketcher.

In this tutorial, you load in and use the models that you created in Using the 3D Sketcher.

Summary of tutorial

The tutorial is divided into five short exercises:

A. Using the model manager

B. Making picric acid from phenol (transferring data between model spaces)

C. Building a model of indigo

D. Measuring the indigo model

E. Generating hard copy

Step-by-step tutorial

This section includes a series of short progressive exercises describing various aspects of the Cerius2 Visualizer: model management, measuring models, and creating Postscript images for printing. You will use the controls in the model manager to complete the steps of this tutorial. The components of the model manager are named in the illustration below.

A. Using the model manager

The model manager is used to view collections of models, sort models by name and by identity number, and discard models from the Visualizer.

1.   Load the models.

If Cerius2 isn't already running, type cerius2 in a UNIX shell to start the program.

Choose Load Model... from the File pulldown menu. Use the Load Model control panel to load models from the previous tutorials: phenol.msi, cyclohexane.msi, and xylene.msi.

The file names, without the ".msi" extensions, appear in the model manager list in the same order that the files are loaded in.

Load butane.msi and napthalene.msi from the -Cerius2-Models/organics directory.

There are now five models in the system. Each is named and numbered on the model manager list. The diamond after the model number is the current model indicator (i.e., the displayed model on which calculations will be performed). The current model is marked with red.

Click the diamond for the different models in the list. After trying the other models, select phenol to be displayed.

The model displayed in the window changes with each selection. Phenol should be the last one displayed.

It is possible to display multiple models on the screen at any one time in a variety of display styles. The Model Display Selector boxes to the right of the model names (under the heading Visible) show which models are displayed. Displayed models are indicated by a filled blue box, the visibility indicator.

The current model is always marked visible.

Check the Model Display Selector box for xylene.

The model window divides into border mode, and the xylene model appears in the top right corner of the window.

Select two more models.

The four models are displayed in border mode, with the current model (phenol) as a large model surrounded by smaller images of the other models.

2.   Try rotating and translating the models.

Use the mouse to rotate and translate the models in various ways.

Notice how all are rotated together.

To the left of the Cerius2 logo are three icons that control display style.

Click on the middle Display Mode icon to get a grid display.

Click the top icon.

It overlays all currently visible models.

Return to border display mode by clicking the bottom icon. Click the filled blue box on the bar above the Model Manager list.

All models should now be displayed.

Click the adjacent empty box icon.

Only the current model is displayed.

Click the Name button immediately above the list.

This orders the models according to the first letter of the model names. (This is case sensitive: names beginning with uppercase letters are listed before those beginning with lowercase letters.)

Try clicking on the ID, A, and Visible buttons to sort models in other ways.

Make napthalene the current model. Click the "-" icon above the Model Manager list and choose the OK button when prompted by the dialog box.

This deletes the model.

Make cyclohexane the current model. Click the "[ ]" icon and choose the OK button when prompted by the dialog box.

This clears the model -- leaving it in the list but removing all data from it.

B. Making picric acid from phenol (transferring data between model spaces)

In the following steps, you learn how to create a new model space and to copy atoms from one model space and paste them into another. You copy the phenol model into the new model space and then use the Sketcher to modify the phenol model to form picric acid (2,4,6-trinitrophenol).

1.   Copy the phenol model.

Set the current model to be phenol. Select all the atoms in the phenol model.

(You can do this with the Select All item from the ccpulldown or use theicon.) All the atoms in the current model will be highlighted.

Choose the Copy option from the Edit pulldown.

This copies the phenol model into the buffer.

Choose the "+" icon on the Model Manager.

A new model space, called Model6, is created.

Choose Paste from the Edit menu.

The phenol model is copied into the new model space.

2.   Edit the model to create 2,4,6-trinitrophenol.

Open the Sketcher, and click the Fragments button. Load the nitro template -- the templates/organic/nitro file. Use the nitro template to create 2,4,6-trinitrophenol.

This molecule, 2,4,6-trinitrophenol, is also called picric acid.

Return the highlight to the selection mode icon once more. Clean the model and give it the name ferrocene in the Model Manager list.
.

C. Building a model of indigo

In the following steps, you use the Sketcher to build a model of indigo. In doing this, you review sketching while also learning some important new things:

1.   Sketch the indigo model.

Create a new model space by choosing the "+" button on the Model Manager.

Rename this new model space "indigo".

You're going to sketch the indigo model in this space.

Because indigo is centro-symmetric, save time by sketching only half of the molecule and then doubling it.

2.   Sketch this fragment.

Hint: you may want to begin by using the template tool with the phenyl fragment.

Choose the Select All item from the Edit pulldown menu or the icon.

All atoms are highlighted.

Choose the Duplicate item from the Edit pulldown menu.

A second copy of the molecule is pasted into the model window.

Pick in the Model window and hold down the left mouse button while you drag the mouse to create a box enclosing the whole of one molecule. Hold down the <Ctrl> key, and use the rotate and translate mouse controls.

Only the selected molecule moves.

Move one fragment relative to the other until the two carbon atoms at the ends of the five-membered rings are about a bond length apart. Link the two fragments by a double bond using the Edit Bond option in the Sketcher. Clean to complete the indigo molecule.

Return the Sketcher to selection mode, and then close the Sketcher control panel.

D. Measuring the indigo model

In the following steps, the Measurements control panel is used to measure the geometry of the indigo model from C. Building a model of indigo.

1.   Measure the distance between two atoms.

Open the Measurements control panel from the Geometry pulldown menu. Click on the Distance icon.

As in the Sketcher, these icons determine the function of the mouse in the model window.

Pick any two atoms.

The distance between them will be marked on the model. You may need to rotate the model or zoom in on it to see the measurement clearly.

2.   Try some of the other measurement tools.

Pick three atoms to measure angles, four to measure torsions or inversions, and six to measure the angle between two planes.

Choose the List Measurements button.

Pick three atoms to measure angles, four to measure torsions or inversions, and six to measure the angle between two planes.

All measurements you have made will be written out to the text window. (You may want to label the model by Numbers to match up the atom serial numbers in the text window with the atoms in the model window.)

3.   Optional step:

Load the model Cerius2-Models/molecular-crystals/pigments/indigo_1.msi.

This is a crystal structure of indigo.

Perform the same measurements on this model as you made on the indigo that you sketched. Compare the results.

Because the model that you made was only crudely minimized by the Clean tool and not affected by crystal packing constraints, you'll notice significant differences between the two models.

Finish by choosing the Delete Measurements button.

This removes all the measure lines and values from the models.

E. Generating hard copy

In the following steps, you output the contents of the model window to a PostScript file, and then print a view of the trinitrophenol model in ball display.

Note

This exercise assumes that you are able to print grayscale PostScript files from your workstation.  

Make trinitrophenol the current model and the only model visible. Set the display style to Ball.

If necessary, rotate, scale, or translate the model in order to get a pleasing view of it in the model window.

Open the Print control panel by selecting Print... from the File pulldown menu. Choose the Page Setup... button to open the Page Setup control panel. On the Page Setup control panel set the background Color popup to White.

The model window background turns to white. (A white background will save toner when the image is printed.)

If you want your name or any other text to appear on the image, type into the User Name and Title Text boxes.

Make sure that no windows or control panels are overlapping the Model Window. Choose File to be the PostScript Destination and enter trinitrophenol.ps in the text box that currently reads CERIUS.ps. Choose the Print button.

A message is written to the text window:


Starting to output image ...
..................Image Captured
The PostScript image will be written to a file called -rinitrophenol.ps.

If you are running Cerius2 on an SGI, choose the Preview button. (If you are running on another platform skip to step 4, below.)

This calls the xpsview program to give a preview of how the PostScript image will look on the page. The preview gives you a chance to alter the model view or the page setup controls before printing.

Close the xpsview window by choosing Exit from the xpsview File menu in the top left corner.

4.   Select one of two methods to print the file.

There are two ways that you can print the file:

The first is to change the PostScript Destination back to Printer, and choose the Print button. If the laser printer that you want is not the system default, you will also need to edit the print command lp -c.

Or

The second is to print the trinitrophenol.ps file from a UNIX shell in the way you would normally print any PostScript file. In the simplest configuration, the print command would be lp trinitrophenol.ps.

On the Page Setup control panel, set the background Color popup back to Black.

The model window background returns to black. (You may have found that the white background is hard on your eyes.)

Note

It is also possible to use the SGI snapshot utility or the xv program to capture rgb images from Cerius2.  

5.   Exit or continue to another tutorial.

Either:

Re-initialize the program by selecting New Session from the File pulldown menu and choosing the Confirm button on the Re-initialize dialog box.

Re-initialization resets all the Cerius2 control panels to their original defaults.

Alternatively, if you want to stop now:

Choose Exit from the File pulldown menu to quit from Cerius2, and choose the Exit button on the Exit Option... dialog box.

Review

In this tutorial you learned how to:

o Display multiple models using the Model Manager

o Remove and clear models from the Model Manager list

o Create a new model space

o Copy atoms between model spaces

o Duplicate selected atoms

o Move one model fragment relative to another

o Measure model geometry

o Generate a PostScript image from the model window

Now, you may want to proceed to the next tutorial, Managing and Modifying Graphs.

Related material

The model manager is discussed in detail in the Cerius2 Modeling Environment. If you haven't already done so, you may want to read the introductory tutorial in Chapter 2 of that book. Other sections of that book that relate to the activities in this lesson include:

Using the main control panel in Chapter 2

Controls and control panels in Chapter 2

Basic building in Chapter 4

Printing models and graphs in Chapter 7


Managing and Modifying Graphs

This chapter contains a tutorial exercise about graphs in Cerius2. It introduces the basics of viewing, modifying, saving, and loading graphs that are created by Cerius2.

This chapter is divided into four parts:

Before you begin

To complete this tutorial, you need a licensed copy of Cerius2 and a directory into which you can save the files you create during the tutorial.

In this tutorial, you learn about the Cerius2 Graphs card. Nearly all the Cerius2 materials science modules produce some output to the graph window. It's important that you understand how the Graphs card works so that you can get the most from your graphical results.

You can use the functions on the Graphs card:

You have opportunity to do all these things in this tutorial exercise.

Cerius2 Graphs terminology

The Cerius2 graph system is based on a hierarchy of three structures:

These are defined below and shown symbolically in the figure.

Data-set

At the bottom of the hierarchy is the data-set. The data-set must be one of three types:

Graph

At the middle level of the hierarchy is the graph. A graph is made up of one or more data-sets displayed on axes.

Gallery

At the top of the hierarchy is the gallery. The gallery is the name given to the group of graphs on display in the graph window at a given time.

Summary of tutorial

The tutorial is divided into six exercises:

A. Loading graphs and using the Gallery Manager

B. Using the Graph Manager and setting plotting attributes

C. Scaling graphs

D. Annotating graphs

E. Saving graph files

F. Creating new graphs and using expert functions

Step-by-step tutorial

This section takes you through a series of short progressive exercises designed to familiarize you with the Cerius2 Graphs card.

A. Loading graphs and using the Gallery Manager

Graphs are saved in files with .grf extensions; this format is described in the "File Formats" appendix of the Cerius2 Modeling Environment.

Files generated using Insight II can also be loaded; these have the extension .tbl and are described at the MSI website, www.msi.com.

The sample graphs that you load for this tutorial were created by the Blends module and the Sorption module. The contents in themselves are irrelevant for this tutorial, and you needn't have licensed the Blends or the Sorption modules to load or manipulate these graphs.

1.   Load the graph files.

Find the Graphs card (on the Tables & Graphs card deck).

If the card is not already at the front of the deck, click on the card title to bring it to the front.

Open the Load Graphs control panel by choosing the
File -> Load Graphs item on the Graphs card. Find the directory Cerius2-Resources /EXAMPLES/data.

Be sure that the filename text box at the bottom of the browser box reads *.grf.

Load the file EM_MM_DG.grf: pick its name on the browser box and then click the Load button.

A confirmation message like the one shown below appears in the text window. The graph appears in the graph window.


Loading Grapher file: ./Cerius2-Resources/EXAMPLES/data/EM_MM_DG.grf
Loading Done.

Place the graph window in the top left corner of the screen over the Model window. Enlarge the graph window (click and drag a corner of the window) until the window is almost as large as the Model window.

You do not need the Model window for this tutorial.

Now load two more graph files from the EXAMPLES/data directory: EM_MM_Emix.grf and Loading.grf. Close the Load Graphs control panel.

2.   Arrange the graphs in the Gallery.

Open the Gallery Manager control panel by choosing Gallery -> Manager from the Graphs card. In the Graph Chooser list box, you will see the names of the three graphs that you've loaded into the system.

They are: Sorption Loading Plot v1, Emix Fit v1, and Thermodynamics v1. Only the Sorption Loading Plot is on display.

Click on the Display check boxes next to Emix Fit v1 and Thermodynamics v1.

Three graphs should now be on display in the graph window.

Notice that the Gallery Width and the Gallery Height popups on the Gallery Manager control panel are both set to Auto. This automatically arranges the graphs to best effect in the graph window.

Change the Gallery Width to 3 and the Gallery Height to 1, and note the result. Try some other width/height combinations before returning both settings to Auto.

Note

If you specify a gallery size that is too small for the number of graphs, graphs will disappear from the Gallery window. To redisplay a lost graph, set the gallery size large enough for all the graphs and then select again the missing graphs on the Graph Chooser list box of the Gallery Manager control panel.  

Set the Label Graphs check box to on. Note the appearance of the graph labels in the graph window. Uncheck the Emix Fit v1 and the Thermodynamics v1 graphs on the Graph Chooser list box to leave only the Sorption Loading Plot on display in the graph window. Be sure that the Edit box for the Sorption Loading Plot graph is toggled on.

B. Using the Graph Manager and setting plotting attributes

Now that you've looked at the Gallery Manager and used it to edit the layout of graphs in the graph window, you learn how to edit the graphs themselves using the Graph Manager.

1.   Open the Graph Manager

Open the Graph Manager control panel by clicking the Graph Manager... button at the bottom of the Gallery Manager control panel. (Alternatively open the control panel by choosing Graph -> Manager from the Graphs card.)

The Graph Manager control panel has a Data-Set Chooser list box that parallels the Graph Chooser browser box in the Gallery Manager. In the Data-Set Chooser are listed the data-sets of all three of the loaded graphs, but only the data-sets currently displayed on the Sorption Loading Plot are checked.

Uncheck the Display button for the Sorption Loading Data 01 data-set in the Data-set Chooser list box to leave the Sorption Loading Data 02 and the Sorption Loading Data 03 data-sets on display in the graph window. Set the check boxes in the Edit column so that the Sorption Loading Data 02 data-set can be edited but the Sorption Loading Data 03 data-set cannot be.

The top portion of your Data-Set Chooser should look like this:

Click the Plotting Attributes... button at the base of the Graph Manager control panel to open the Plotting Attributes control panel. Change the Style popup (currently reading Point) to Line.

Observe that the data plot style changes from points to a line.

Use the Color popup (currently reading YEL) to change the line color to light green (LGR).

Edit the Label text box.

Observe that the text for the green line is modified.

Edit the Y Offset text box: for example, change it to 30.

Observe the Y shift in the data plotted.

2.   Editing more than one dataset at a time.

On the Graph Manager control panel, select all the Sorption data-sets for both Display and Edit.

Return to the Plotting Attributes control panel.

Notice how the Style and Color popups show the -> symbol and the Y Offset and Label text boxes are blank. This is because the attributes are not the same for all of the data-sets. The 01 and 03 data-sets are in point style, and the 02 is in line. Furthermore, the colors of all three are different, and the key labels and Y offsets also differ.

Set the Y Offset to 0 (zero).

Experiment with the Style popup: set it to Histogram, to Delta, to Point, and finally to Line.

You could also use the Color popup to set all the lines to the same color or the Label text box to make the labels the same for all of the datasets -- but these would not be very useful.

Close the Plotting Attributes control panel.

C. Scaling graphs

In the following steps, you learn four ways to alter the scale of a graph:

1.   Use the text boxes to scale a graph.

Open the Graph Scaling control panel by choosing the Graph -> Scaling item on the Graphs card or by choosing the Graph Scaling... button at the foot of the Graph Manager.

Enter a value of -1e5 in the Min X text box and a value of
-20 in the Min Y text box.

Notice the effect on the graph axes.

2.   Try the zoom and scroll buttons and observe their effect.

increases the zoom level
decreases the zoom level
resets the graph to its original scale

You can use the % text box associated with the zoom buttons to specify the magnitude of the zoom buttons' effect.

Alter its value from 5 to 20% and notice the effect on the zoom rate.

Use the zoom buttons to adjust the graph scale so that the origin of the graph is not visible. Now use the scroll buttons until you can see the origin back on the graph.

The scroll buttons work very similarly to the zoom ones. The arrows indicate the scroll directions, and the central button resets the scale. The scroll % text box is used to set the percentage of the axis scrolled in one button click.

Next, you try the mouse/keyboard shortcuts for scrolling and zooming.

3.   Manipulate the graphs with the mouse.

Try mouse scrolling: put the cursor within the axes of the graph, and move the mouse while holding down the middle mouse button.

By moving the mouse diagonally, you can scroll in both the x and y directions at once.

Try mouse zooming: put the cursor within the axes of the graph, and move the mouse while holding down both the <Shift> key and the middle mouse button.

Click the Zoom Reset button to return the graph to its original scale.

Click the Axis Scaling popup to view the choice of scale settings. Then click again using the right mouse button to display the help text that describes the five sorts of scaling.

When you have finished reading the help text, click it to make it disappear.

Try the Auto-XY setting, and then the Nice-XY setting.

Notice the change in the Max Y value when you do this.

Exit from the Graph Scaling control panel.

D. Annotating graphs

In the following steps, you learn how to use the Annotation Option popups to display (or not display) graph axes and various text items. You also edit the text for the graph title and axes.

(Remember, you've already learned how to edit the text for the graph key, or legend, in the Edit the Label text box. step on page 102.)

1.   Experiment with the Graph Annotations.

Open the Graph Annotation control panel by choosing Graph -> Annotation from the Graphs card or by clicking the Graph Annotation... button at the foot of the Graph Manager.

Change the popup options from Yes to No and back again so you can identify the various components of the graph: Axes, Label, Title, and Key.

The Cell option only applies to 3D graphs.

Note

You will never want to select --> because it has no effect. However, Cerius2 may set a popup to --> to indicate that this option is not unique for the graphs in the edit list.  

Edit the X Axis Label, the Y Axis Label, and Graph Title text boxes.

Notice that the text on the graph is updated with the new labels and title.

Close the Graph Annotation control panel.

E. Saving graph files

A graph (.grf) file can contain a single graph or multiple graphs. Each of the graph files that you loaded in section A (page 96) contained only one graph. In the following steps, you save a gallery of two graphs into one .grf file. Then, you reset the Graphs card and load back in the two-graph file that you saved.

1.   Save multiple graphs into one file.

On the Gallery Manager control panel select the Emix Fit and the Thermodynamics graphs to be displayed (and do not display the Sorption Loading Plot graph).

Open the Save Graphs control panel by choosing File -> Save Graphs from the Graphs card. Type two_graphs.grf into the filename text box and click the Save button.

A confirmation message like the one shown below is written to the text window.


Saving Gallery File: ./two_graphs.grf
Saved GRAPHS file.
2.   Delete graphs from the Gallery.

Choose Reset from the Graphs card, and choose Re-initialize Graphs from the dialog box to confirm the reset.

Reset deletes all graphs from the Gallery Manager and returns all Graphs controls to their default settings.

The Gallery and graph window should now be empty.

3.   Load the saved file of multiple graphs.

Open the Load Graphs control panel by choosing File -> Load Graphs from the Graphs card.

You may have to move up several directories to return to your working directory where the two_graphs.grf file is located.

Load back the two_graphs.grf file.

F. Creating new graphs and using expert functions

In this final section of the Graphs tutorial, you create a new User Graph and copy datasets from other graphs onto it. You also learn how to use the Expert Functions: Clean, Data Placement, and Keep Same.

Clear

In part E (page 108), you deleted all the graphs from the Gallery Manager by resetting the Graphs card. In the following steps, you learn how to selectively remove graphs from the gallery.

Data Placement & Keep Same

There is another version of the Emix Fit graph saved in the file Cerius2-Resources/EXAMPLES/data/EM_PM_Emix.grf. Here, you load this graph to illustrate the usefulness of the Keep Same and Data Placement options of the Graphs Expert Functions control panel. In this exercise, you load the second version from file, but, in practice, it will usually be a Cerius2 module that creates a new version of a graph and pipes it into Graphs. However, the principle is the same regardless of whether the graph is loaded or piped. The Data Placement control lets you decide how to treat new versions of the same graph. The Keep Same controls let you apply the graph formatting that you've specified for one graph to subsequent versions of the graph.

1.   Create a User graph.

Open the Gallery Manager control panel by choosing Gallery -> Manager from the Graphs card. Create a new User Graph by checking the third (empty) Display box in the Graph Chooser list box.

This causes the list box to update with the name User Graph v1, and an empty set of axes appears in the graph window.

Next, you copy a dataset from the Thermodynamics graph to the new user graph.

Be sure that the Edit box of the User Graph is toggled on and that the Display boxes for the Emix Fit and the Thermodynamics graph are toggled off. Open the Graph Manager control panel by clicking the Graph Manager... button on the Gallery Manager control panel. Toggle on the Display boxes for the Enthalpy and the Free Energy data-sets.

These two datasets are plotted on the User Graph.

Optional: You may want to review what you learned in the previous exercises (pages 100-106) by customizing this graph: altering plot style, changing colors, scaling axes, and editing text.

2.   Delete extra graphs from the Gallery.

Open the Graphs Expert Functions control panel by choosing Gallery -> Expert Functions from the Graphs card. In the Graph Chooser list box of the Gallery Manager control panel, choose to display only the Emix Fit graph.

Click the CLEAR button on the Gallery Manager/Graphs Expert Functions control panel.

Clicking the CLEAR button deletes from the Gallery Manager all graphs that are not currently on display in the graph window. The Graph Chooser list box should now look like the figure below (the Thermodynamics and the User Graphs graphs are gone).

3.   Add a second graph.

Set the Data Placement popup of the Graphs Expert Functions control panel to Additional Graph. Open the Load Graphs control panel (File -> Load from the Graphs card), and load the file Cerius2-Resources/EXAMPLES/data/EM_PM_Emix.grf.

There are now two Emix Fit graphs listed in the Graph Chooser list box, and two Emix Fit graphs on display in the graph window.

4.   Delete the second graph.

On the Gallery Manager control panel, set the second version of the Emix Fit graph to off. Then choose the CLEAR button of the Graphs Expert Functions control panel.

This deletes the second version of the Emix Fit graph.

5.   Combine the data-sets from the two graphs.

Set the Data Placement option to Combine With Existing. Re-load the Cerius2-Resources/ EXAMPLES/data/EM_PM_Emix.grf graph file.

This time, when you load the second version of the Emix Fit graph, the data-sets from the two versions are combined and displayed on the same axes.

Now set the Data Placement popup to Replace Existing and reload the Cerius2-Resources/ EXAMPLES/data/EM_PM_Emix.grf.

Now, the older version of the graph is overwritten by the new one. Replace Existing is the default for Data Placement; when a computational module creates versions of the same graph, the Replace Existing option prevents these many obsolete versions from accumulating in the gallery.

6.   Make two graphs with identical axes.

Open the Graph Scaling control panel (Graph -> Scaling from Graphs card), and set the scale of the X axis to 300 - 600 and the scale of the Y axis to -1.3 - 2.0.

Open the Plotting Attributes control panel (Data -> Plotting Attributes from Graphs card), and set the Color popup to your favorite color.

Set the following option values in the Graphs Expert Functions control panel: Data Placement Additional Graph
Keep Same: Graph Axes On
Keep Same: Plot Attributes On

Reload the file Cerius2-Resources/EXAMPLES/ data/EM_MM_Emix.grf.

Because the Keep Same options are toggled on, the scale and plot attributes for the new graph are not the ones stored in the file, but instead, they are the same as the previous version of the graph.

7.   Exit or continue to new tutorial.

If you are going to continue immediately with another tutorial:

Re-initialize the program by selecting New Session from the File pulldown menu and choosing the Confirm button on the Re-initialize dialog box.

Re-initialization resets all the Cerius2 control panels to their original defaults.

Alternatively, if you want to stop now:

Choose Exit from the File pulldown menu to quit from Cerius2 and then choose the Exit button on the Exit Option... dialog box.

Checklist review

In this tutorial you learned how to...

o Load and save graph files

o Display graphs using the Gallery Manager

o Display data-sets using the Graph Manager

o Alter the data plot style and color

o Specify a Y-offset value for data-sets

o Edit text for graph titles, axes, and keys

o Use the Expert Functions to selectively remove graphs and to control version overwrites

What next?

In this tutorial exercise, you have encountered only the xy linear graph. Many Cerius2 modules also create 2D cartesian and polar graphs. These graphs can be manipulated (loaded, saved, modified) in much the same way as the xy graphs. When a polar or cartesian contour graph is selected for editing, the appearance of the Graph Scaling and Plotting Attributes control panels are automatically altered to show the special controls you need to modify 2D plots. You may now want to run one of the Cerius2 example files that creates a 2D graph and practice scaling it and modifying its plot attributes. For more information about 2D plots, see the Graphs chapter of the Cerius2 Modeling Environment.

These are some of logfiles that generate 2D plots:

ANL_exafs_SiAs.log ANL_hrtem_particle.log
ANL_gsas_inorganic.log QM_mopac_SN2.log
POL_diff_crys.log POL_diff_am.log

Alternatively, you may want to learn more about the graphs in Cerius2 by reading the user manual selection listed under the Related material heading below.

You may now want to proceed to the tutorials in Cerius2 Materials Sciences Tutorials or Cerius2 Life Sciences Tutorials, or you may prefer to explore other Cerius2 modules on your own with the aid of the Cerius2 user's reference books and the on-line help text.

Related material

1.   The Graphs chapter of Cerius2 Modeling Environment.

2.   The File Formats appendix of Cerius2 Modeling Environment.


Basic Molecular Simulation: Minimization & Dynamics

This chapter shows how to use molecular mechanics techniques to simulate the behavior of a system on an atomistic scale. The example used is a simulation of beta-D-glucose in water.

The tutorial has four parts:

A. Building an initial model of the glucose molecule using the 3D Sketcher

B. Refining the model using molecular mechanics minimization to optimize bond lengths, angles etc.

C. Adding water (for the purposes of this tutorial only a small amount of water is included)

D. Simulating the behavior of the glucose/water system over time using molecular dynamics.

The purpose of this lesson is to introduce the most basic tools for atomistic level simulation. More advanced simulation and analysis options are described in Cerius2 Tutorials--Materials Science and --Life Science.

A. Build a model of the glucose molecule

1.   Construct a glucose model

To start Cerius2, enter cerius2 in a UNIX window. After Cerius2 starts, open the 3D Sketcher from the Build pulldown and construct a model of glucose.

See Using the 3D Sketcher, for guidance on using the 3D Sketcher.

Click and hold the CLEAN button in the Sketcher control panel to clean the structure.

B. Minimize the glucose energy

During a molecular mechanics simulation the movement of the atoms (the changing of bond lengths, angles, etc., and the interactions between non-bonded atoms and molecules) are controlled by a mathematical function usually referred to as the forcefield. The accuracy of a molecular mechanics simulation depends critically on the choice of forcefield. The PCFF forcefield has been shown to give accurate models of a wide range of organic molecules and polymers and is used in this tutorial. Information about other forcefields and theoretical background about atomistic-level simulation using molecular mechanics can be found in the Forcefield-based Simulations book.

1.   Load the forcefield.

Open the OFF SETUP card deck and locate the OPEN FORCEFIELD card. Load the pcff_300_1.01 forcefield.

You may need to scroll down through a number of other forcefields to locate pcff. Note that references are written to the textport for your information.

A molecular mechanics forcefield contains different sets of parameters (different bond lengths, etc.) for atoms in different bonding environments. For example, the oxygen in a carbonyl group will be modeled using different parameters than those used for oxygen in an ether group or in water; carbon in benzene will be given different parameters to carbon in ethanol or in ethene. Before running the simulation it is necessary to identify these different atom types.

2.   Identify the atom types.

Open the Atom Typing menu by holding down the Typing option on the OPEN FORCEFIELD card and selecting Atoms. Select the Run button next to Calculate using typing rules to ask Cerius2 to identify the atoms types for all the atoms in the current model (the glucose molecule).

Change the Atom Labeling popup on the tool bar from NO LABEL to FFTYPE, and then to CHARGES.

Note that the oxygen atoms and carbon atoms in different bonding environments have been identified as different atom types. Once the atom types and bonding environment were identified, partial charges were automatically allocated to each atom.

3.   Minimize the energy of the model.

Select the MINIMIZER card from the OFF METHODS card deck (it is usually the first card in this deck).

During molecular mechanics minimization the positions of the atoms in the model are altered to find lower energy conformations of the molecule.

Open the Energy Minimization menu by selecting Run from the Minimizer card.

The Smart Minimizer works through a series of minimization algorithms that have been found to give the most efficient minimization for molecular models. Pick the Smart Minimizer with the right-hand mouse button to read about the methods used. With the options set as shown, up to 500 iterations (changes in the molecular model which reduce its energy) will be carried out. The minimization will stop either at 500 interactions, or earlier if a minimum energy structure is located.

Select the Preferences... button to see the convergence criteria used.

Select the Run button to minimize the energy.

C. Add water to the system

1.   Create water molecules with the Sketcher

If the 3D-Sketcher menu is not still active, open by selecting it from the Build pulldown menu. Change the Sketch with element to oxygen, and double-click at five or six random locations around the glucose molecule.

This sketches five or six isolated oxygen atoms --shown by default as small red stars.

Select the H ADJUST button on the Sketcher menu to automatically add hydrogens to satisfy the oxygen valences.

This produces five or six water molecules.

D. Simulating the behavior of the glucose/water system

1.   Identify the atom types

Return to the OFF SETUP card deck and open the Forcefield Atom Typing menu (select Typing/Atoms on the OPEN FORCEFIELD card). Calculate the atom types using typing rules.

This identifies the atom types in the water molecules and assigns appropriate partial charges.

You can select CHARGES from the Atom labelling popup menu to view these charges. Change back to NO LABEL before continuing.

Before carrying out a molecular dynamics simulation it is always wise to run a few steps of molecular mechanics minimization to relax any "hot spots".

Select the MINIMIZER card from the OFF METHODS card deck. Open the Energy Minimization card by selecting Run from this card. Reduce the Maximum Number of Interactions to 200 and Minimize the Energy.

The most significant interactions between the glucose molecule and the surrounding water will be hydrogen bonds. The PCFF forcefield does not contain an explicit hydrogen bond term, but correctly models hydrogen bonding as a result of electrostatic and van der Waals interactions between non-bonded atoms. You can use the Edit H-bonds command under the Build pulldown menu to view the predicted hydrogen bonding.

2.   Calculate hydrogen bonds

Select the Build/Edit H-bond from the top menu bar. Click the Preferences... button to see the criteria that will be used to judge the presence of hydrogen bonds. Click the CALCULATE button to display the pairs of atoms that fall within the H-bond criteria.

Turn on the Enable Automated Recalculation so that the changing hydrogen bonding pattern will be displayed during the molecular dynamics simulation.

3.   Simulate the system dynamics

Select the Dynamics Simulation card from the OFF METHODS card deck and open the Dynamics Simulation menu by selecting Run.

This model is fairly small so the number of dynamics steps can be increased to 2000 without making the simulation too time consuming. Each step simulates 0.001 ps, so the entire run will simulate the behavior of the glucose/water system for 2 ps.

A certain amount of analysis is carried out during the simulation and results will be displayed in the graphics and graph windows, but further analysis, after the simulation completes, requires data to be saved.

Select the Trajectory... button at the bottom of the Dynamics Simulation menu, or select Trajectory -> from the Dynamics Simulation card to open the Dynamics Trajectory Output menu.

Select to Create a Dynamics Trajectory file with an update every 100 steps (Trajectory Update Frequency). This means that during the 2000-step simulation, 20 trajectory frames will be saved for later analysis. Type in glucose_water as the Filename Prefix.

Return to the Dynamics Simulation menu and select the Run Dynamics button.

As the simulation continues, the structure on the graphics screen is updated every 25 steps.Notice the changing H-bond pattern. The graphs show that the temperature of the simulation fluctuates around the requested 300K.


What to do next

This simple simulation might be useful to provide insight into hydrogen bonding patterns and dynamic behavior. To obtain a more quantitative measure of hydrogen bonding between water and beta-B-glucose the calculation could be repeated with a larger number of water molecules using periodic boundary conditions. The resulting trajectory file could be analyzed to predict the pair correlation function between the hydroxyl groups on the glucose and the water.

The area under the h-bond peaks is proportional to the number of hydrogen bonds in the system and can be compared to results for other systems (for example to other carbohydrates).

It is known that monosaccharides in water undergo structural isomerization. Processes that involve bond breaking and formation will not be observed using molecular mechanics, but could be studied using quantum mechanics models (for example, using FastStructure).




Last updated April 13, 1999 at 06:36PM Pacific Daylight Time.
Copyright © 1999, Molecular Simulations Inc. All rights reserved.