26.8 Example: Inserting rigid bodies

For this, we will approximately duplicate the structure model used in the MOF structure from Silvana Urcia-Romero et al. (2022), Crystal Growth & Design, 2382–2391.1. The CIFs from this paper can be downloaded for free from the CCDC 2. Any of the CIFs from that paper will do. After the coordinates are read, and the “Draw Atoms” or “Draw Options” phase tab is selected, the structure will appear as shown in 26.22.

PIC

Figure 26.22: The asymmetric unit for the MOF structure being used to demonstrate use of rigid bodies. The pink atom is Co, which is bonded to a bipyridine (half shown), a pyrazinedicarboxylic acid unit and a water molecule. Also present is a \(\rm CO_2\) molecule. Note that the view point has been shifted away from the center of the unit cell.

26.8.1 Example 1: Inserting the pyrazinedicarboxylic acid ion rigid body

To insert the pyrazinedicarboxylic acid ion rigid body that was previously created, select the “RB Models” phase tab. If for some reason the structure is not plotted, select the “Draw Atoms” or “Draw Options” phase tab and then return to the ‘RB Models” tab. The ‘RB Models” tab posts the menu labeled “Edit Body” and from that menu select the “Locate & Insert Rigid Body” menu command. This will open a window where one can select from the three rigid bodies that were previously defined. The body we want was called “pzdc”; select this and press OK. The data window will appear, as seen in Fig. 26.19 and the structure visualization will appear as shown in 26.23, where now the rigid body is also displayed.

PIC

Figure 26.23: The asymmetric unit for the MOF structure with the addition of the pyrazinedicarboxylic acid rigid body in the process of being inserted. The stick figure with green lines for bonds shows the rigid body location and orientation.

Our next task will be to orient the pyrazinedicarboxylic acid rigid body on the approximate positions of the atoms that are in that group. This will be very hard to see with these atoms shown as space-filling. Press the “Ball & Sticks” or “Sticks” button to change the display according to your preference. As noted before (in §26.5) we can identify atoms to be paired by highlighting them. The most easily identified atoms to me are the ring nitrogens and the carboxyilic carbons. In Table 26.1 I have provided a list of atoms to be paired.

rigid body structure
N2 N13
N11 N14
C4 C6
C8 C7
Table 26.1: Easily identifiable paired atoms for the pyrazinedicarboxylic acid rigid body and corresponding atoms in the structure.

Assign these four atoms, as seen in Fig. 26.24. Then press the “Set both” button and note that the rigid body is now largely superimposed, as seen in Fig. 26.25 and the assignments table (Fig. 26.26) now shows the C and N paired atoms are all within a few tenths of an Ångstrom. The O atoms can be brought into better agreement by adjusting the torsion angles with the sliders.

PIC

Figure 26.24: The assignments table for the pyrazinedicarboxylic acid rigid body, which shows how atoms are paired. This is after four atom assignments have been made, but before the rigid body position has been updated.

PIC

Figure 26.25: The pyrazinedicarboxylic acid rigid body superimposed on the structure, in the process of being inserted.

PIC

Figure 26.26: The assignments table for the pyrazinedicarboxylic acid rigid body after four atom assignments have been made and the “Set both” button has been used to superimpose the bodies. Note the close distances between the C and N atom pairs.

Finally, complete the rigid body insertion process by pressing the “Add” button. The two sets of superimposed atoms are replaced by a single set, where the bonds are now colored in orange, to indicate that this is a rigid body, as seen in 26.27. The positions of the rigid body atoms are updated in the Atoms tab and the coordinates and \(\rm U_{iso}\) values are shaded to indicate that these values are generated from the rigid body, as seen in Fig. 26.28.

PIC

Figure 26.27: The pyrazinedicarboxylic acid ion rigid body after insertion. The orange bonds indicate that this is a rigid body.

PIC

Figure 26.28: The contents of the Atoms tab after the pyrazinedicarboxylic acid ion rigid body has been inserted. Note the shading of the coordinates and \(\rm U_{iso}\) values to indicate that these values are generated from the rigid body.

26.8.2 Example 2A: Inserting the bipyridine molecule rigid body

Before beginning this example, you should save the GSAS-II project so that we can return to this point, where only the first rigid body has been inserted. This rigid body is a bit different from the previous because it has internal symmetry. It is duplicated by a center of symmetry. This \(\overline {1}\) located at (0,\(\frac {1}{2}\),0), so this is where the origin of this rigid body needs to be placed.

Again, select the “RB Models” phase tab, if needed. To insert the bipyridine molecule rigid body from the “half-bipy” rigid body that was previously created; again use the “Edit Body”/“Locate & Insert Rigid Body” menu command. This will open a window where one can select from the three rigid bodies that were previously defined. The body we want was called “half-bipy”; select this and press OK. The data window will appear similar to what was seen before (Fig. 26.19) and the half-molecule bipyridine rigid body is also displayed.

To set the origin location, change the “Origin” y value to 0.5 at the top of the data window and select the “Lock” selection so this does not get changed accidentally. At this point, given that the origin is fixed, it is not that hard to use the Alt+Left mouse and Alt+Middle mouse drag motions to align the rigid body to the atoms in the structure. This will not be easy to explain in writing how to do this, but let me try. Note that the following paragraph is optional as the following one duplicates the alignment process.

The initial plot appears with the rigid body approximately opposite where we want it, as seen in Fig. 26.29. Holding the Alt key down, drag the middle mouse button (Alt+middle mouse) vertically until the two rings are approximately aligned, as seen in Fig. 26.30. While this looks to be fairly reasonable alignment, rotating the view by approximately \(90^\circ \), by dragging with the left mouse button (no Alt key) produces a view that shows the rings are not in the same plane, as seen in Fig. 26.31. Rotation of the rigid body using Alt+left mouse, will bring the rings to be closely overlapped.

PIC

Figure 26.29: Initial location of the half-bipyridine molecule rigid body during the insertion process.

PIC

Figure 26.30: Location of the half-bipyridine molecule rigid body during the insertion process after rotation using the Alt+middle mouse button.

PIC

Figure 26.31: Location of the half-bipyridine molecule rigid body during the insertion process after rotation using the Alt+middle mouse button and rotation of the view.

The other way to align the rigid bodies is to assign atoms, as was done before. Note that the rigid body N1 atom must be paired to the N15 atom in the structure. The two rigid body atoms bonded to N1 are C2 and C3, which can be seen from the atom labels in the plot. In the structure, it can be determined that C8 and C12 are bonded to N15 using the “Crystal Highlight” option. As before, assign these atoms. Assign rigid body atoms C2, N1, and C3, respectively, to structure atoms C8, N15, and C12, respectively, and press the “Set Orientation” button. (Note that it does not matter if C8 is assigned to C2 or to C3, so the two C atoms assignments could be reversed.)

Finally, press “Add” to complete the insertion of this rigid body. The half of the bipyridine molecule will now appear with orange bonds, as was seen before for the pyrazinedicarboxylic acid.

26.8.3 Example 2B: Inserting the bipyridine molecule rigid body with a duplicate atom

The last example will be to repeat the previous rigid body insertion, but to use a version of the rigid body that includes atoms that are duplicated by symmetry. The rigid body insertion in previous step needs to be reversed before this step can be performed. Either reload the GSAS-II project (.gpx) file that was saved in the last step or delete the rigid body insertion by selecting the half-bipy:0 body on the “RB Models” tab and then press the “Delete” button.

Follow approximately the same process as before to insert the bipyridine molecule rigid body. Again use the “Edit Body”/“Locate & Insert Rigid Body” menu command, but this time use the “half-bipy+1” rigid body that was previously created, . This will open a window where one can select from the three rigid bodies that were previously defined. Again, set the origin location by changing the “Origin” y value to 0.5 at the top of the data window and select the “Lock” selection. Assign rigid body atoms C2, N1, and C3, respectively, to structure atoms C8, N15, and C12, respectively, and press the “Set Orientation” button.

One important extra step is needed relative to the previous example. While atoms C1-C5, N1, and H1-H4 match atoms in the structure, there is no matching atom for C6. The default action would be to pair this atom with one of the \(\rm CO_2\) C atoms, but here we need is for GSAS-II to add this atom to the structure when the rigid body is inserted. From the pull-down menu in the “Assign as atom” for C6, select the “Create new” option. Then press the “Update Assignments” button and the assignments table should appear, as in seen in Fig. 26.32.

PIC

Figure 26.32: The assignments table for the bipyridine molecule+1 rigid body. Note that atom C6 has been designated as a new atom.

After the “Add” button has been pressed, the rigid body information is updated in the Atoms tab, which appears as in seen in Fig. 26.33. Note that the C6 atom in the rigid body has been added to the atoms list as atom #29 (the second atom labeled as RbC29.) Note that since this atom is a symmetry duplicate, the occupancy has been set as zero during the insertion process.

PIC

Figure 26.33: The contents of the Atoms tab after the bipyridine molecule+1 rigid body has been inserted. The new atom from the rigid body has been inserted as the last atom.