Building and Manipulating Crystal Structures

Tutorial Created with Software Release: 2024-3
Topics: Catalysis & Reactivity, Energy Capture & Storage, Metals, Alloys & Ceramics, Thin Film Processing
Methodology: Periodic Quantum Mechanics
Products Used: MS Maestro

Tutorial files

0.4 MB

This tutorial is written for use with a 3-button mouse with a scroll wheel.
Words found in the Glossary of Terms are shown like this: Workspacethe 3D display area in the center of the main window, where molecular structures are displayed

 

Tip: You can hover over a glossary term to display its definition. You can click on an image to expand it in the page.
Abstract:

 

In this tutorial, we will learn how to import, manipulate and build crystal structures in MS Maestro.

 

Tutorial Content
  1. Introduction to Crystals and Periodic Boundary Conditions 

  1. Creating Projects and Importing Structures

  1. Manipulating an Inorganic Crystal Structure

  1. Manipulating an Organic Crystal Structure

  1. Building a Crystal Structure with Crystal Builder

  1. Introducing Vacancies and Defects into a Crystal Structure

  1. Handling Atomic Disorder

  1. Conclusion and References

  1. Glossary of Terms        

1. Introduction to Crystals and Periodic Boundary Conditions

In order to predict bulk and surface properties of infinite or semi-infinite crystalline materials, it is first necessary to build and manipulate a crystal structure using the MS Maestro interface. The procedures described in this tutorial are therefore fundamental to a variety of applications, including but not limited to:

  • Calculating band structure, mechanical, dielectric, magnetic and thermodynamic properties of inorganic materials
  • Creating compositional phase diagrams
  • Examining material defects
  • Probing surface degradation and reactivity
  • Constructing multi-interfacial systems containing a crystalline component
  • Modeling bulk and surface solid alloys

A crystal is a periodic array of atoms. As such, it is defined by the crystal unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors - three translation vectors a,b,c and the atomic site coordinates (basis of atoms). Translation vectors (also known as lattice parametersthe lengths of the unit cell vectors (a, b, c) and angles between them (ɑ, β, ɣ) which define the unit cell), can be also specified by six numbers - the length of the translation vectors (a, b, c) and the angles (ɑ, β, ɣ) between them (Figure 1). The primitive unit cellthe smallest unit cell is the smallest repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along a, b, and c.

Figure 1. The unit cell with periodic boundary conditions

The smallest asymmetric subset of atomic sites in the unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors constitutes the asymmetric unitthe smallest asymmetric subset of atomic sites in the unit cell. In conjunction with a space group it uniquely defines the unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors . The space groupthe point group operations (reflections and rotations) and translations that define the unit cell comprises the point group operations (reflections and rotations) and translations. Materials Science Suite supports all three-dimensional space groups as defined in the International Tables for Crystallography. For a full description of space groups, see the International Tables for Crystallography.

For more detail on introductory solid state concepts, see the References section at the end of this tutorial.

In this tutorial, we will learn how to import crystal structures into MS Maestro for viewing and manipulation. We will use the crystal tools and our understanding of periodic boundary conditions to prepare model systems. In addition, we will explore building crystal structures directly, as well as procedures for introducing defects into a system.

2. Creating Projects and Importing Structures

At the start of the session, change the file path to your chosen Working Directorythe location where files are saved in MS Maestro to make file navigation easier. Each session in MS Maestro begins with a default Scratch Projecta temporary project in which work is not saved, closing a scratch project removes all current work and begins a new scratch project, which is not saved. A MS Maestro project stores all your data and has a .prj extension. A project may contain numerous entries corresponding to imported structures, as well as the output of modeling-related tasks. Once a project is saved, the project is automatically saved each time a change is made.

Structures can be built in MS Maestro or can be imported using File > Import Structures (or drag-and-dropped), and are added to the Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and Project Tabledisplays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data. The Entry Lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion is located to the left of the Workspacethe 3D display area in the center of the main window, where molecular structures are displayed. The Project Tabledisplays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data can be accessed by Ctrl+T (Cmd+T) or Window > Project Table if you would like to see an expanded view of your project data.

  1. Double-click the Materials Science icon

Figure 2-1. Change Working Directory option.

  1. Go to File > Change Working Directory
  2. Find your directory, and click Choose
  3. Pre-generated input and results files are included for running jobs or examining output. Download the zip file here: schrodinger.com/sites/default/files/s3/release/current/Tutorials/zip/builders_crystals.zip
  4. After downloading the zip file, unzip the contents in your Working Directorythe location where files are saved for ease of access throughout the tutorial

Figure 2-2. Save Project panel.

  1. Go to File > Save Project As
  2. Change the File name to crystals_tutorial, click Save
    • The project is now named crystals_tutorial.prj

Figure 2-3. Importing .cif files.

While building crystal structures from scratch is do-able (demonstrated in Section 5), it is typically much simpler to begin by importing a structure file. The Materials Science suite supports .cif and .pdb crystal file formats. Some common sources of structure files are:

In this tutorial, the provided files are .cif and should now be located in your Working Directorythe location where files are saved following Step 7 above

  1. Go to File > Import Structures
    • Import panel opens
  2. Navigate to the tutorial files and select (Cmd + Click or Ctrl + Click) both alpha-glycine.cif and SiO2_alpha_quartz.cif
  3. Click Open
    • The two crystal structures are imported

Figure 2-4. Renamed and reordered imported entries.

Depending on the origin of a structure file, it may have a nondescript name. Such is the case here. Let’s rename and reorder our files

  1. Double-click on each name in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion to rename the entries as follows:
    • 301K → alpha_glycine
    • a-quartz → alpha_quartz
  2. Click and drag the newly named files in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion to a new order: alpha_quartz then alpha_glycine corresponding to the order that we will utilize these files in the tutorial
    • The names and order should correspond to the Figure.
  3. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue alpha_quartz
    • The alpha_quartz entry will be selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and will appear in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Note: Please refer to the Glossary of Terms for the difference between includedthe entry is represented in the Workspace, the circle in the In column is blue and selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries

Additional Material for Importing Structures:

Though not covered explicitly in this tutorial, a relevant panel to be familiar with is:

The panel allows you to query the database of the Materials Project for structures that match a set of criteria. The structures are downloaded and can be incorporated into the project.

3. Manipulating an Inorganic Crystal Structure

The crystal tools available in MS Maestro are routinely used to visualize and then prepare an inorganic crystal structure for subsequent calculations. Here, we will first consider alpha-quartz, the most common crystalline form (polymorph) of SiO2. In this section we will familiarize ourselves with the Periodic Structure Tools, prepare a supercell, and then create a cell representing the system under a stress or strain using the Redefine Lattice panel.

Figure 3-1. Applying style settings.

  1. With the alpha_quartz entry selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includedthe entry is represented in the Workspace, the circle in the In column is blue, in the Style () menu, click ball-and-stick representation and Color Atoms
    • The atoms now appear as spheres rather than points
    • The coloring is updated

Note: The ball-and-stick representation is recommended for viewing solid state structures. As such, this step will always be done first when importing and analyzing crystal structures

Figure 3-2. Menus directly available from the workspace that are useful for crystal structures.

When initially visualizing crystal structures, we will typically use two menus directly available from the workspacethe 3D display area in the center of the main window, where molecular structures are displayed: the Periodic Structure Tool Window () and the Workspace Configuration Panel ()

Note: Whenever a periodic structure is includedthe entry is represented in the Workspace, the circle in the In column is blue in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed, next to the button, you will see the associated space group as well as whether you are visualizing the asymmetric unit (ASU). Here the SiO2 structure in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed is currently the asymmetric unit and the space group is P3121 (space group 152)

Figure 3-3. Building a unit cell.

In the next few steps, we will build a unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors from the asymmetric unit. Note that in practice, any of these steps (2-17) can be done concurrently, but we will do them separately to demonstrate the individual role of each selection

  1. Click on the button to open the Periodic Structure Tool Window
  2. Go to Build Cell
  3. Click Apply
    • All of the atoms comprising one unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors are shown
    • Bonding between expected atoms will not necessarily be shown

Figure 3-4. Adding bonds to a unit cell.

  1. Click on the button to open the Periodic Structure Tool Window
  2. Go to Build Cell
  3. Click Recalculate Connectivity
  4. Click Apply
    • In addition to all of the atoms comprising one unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors being shown, bonds are included

 

Note: If unexpected bonds are shown (or conversely, expected bonds are not shown), you can manually set the tolerance for bonds from the Workspace Configuration Panel () by adjusting the bond options here:

Figure 3-5. Translating the cell to fit into the unit cell box in the workspace.

  1. Click on the button to open the Periodic Structure Tool Window
  2. Go to Build Cell
    • Note that Recalculate Connectivity is still selected. Let’s maintain that selection
  3. Click Translate to First Unit Cell
  4. Click Apply
    • Now, all of the atoms comprising one unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors are shown (again with bonds) and they are translated to the unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors represented by the blue box in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Figure 3-6. Visualizing extents.

 

We will now explore the extents tool. The typical use case for the extents tool is to visualize how the unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors translates in 3D space to generate the infinite crystalline material

  1. Click on the button to open the Periodic Structure Tool Window
  2. Go to Build Cell
    • Note that Translate to First Unit Cell and Recalculate Connectivity are still selected. Let’s maintain those selections.
  3. Click Extents
  4. Insert 1 for the +A, +B and +C directions
  5. Click Apply
    • An extended view of the crystal is displayed in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed
    • Cell 2 x 2 x 2 is shown in green in the bottom right

Note: We have not changed anything about the crystal structure or its periodic boundary conditions in these steps. Stated another way, these steps are purely visual so far. However, extents should be removed before running a calculation

Note: If you select Create new entry from the Build Cell menu, rather than altering the current entry in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed, a new entry will be added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion with the associated build parameters

Note: After doing these steps, you can always deconstruct the cell and return the asymmetric unit by clicking Revert to ASU from the Periodic Structure Tool window

Figure 3-7. Toggling the unit cell box on or off.

A related useful feature is the ability to show or hide the unit cell box

  1.   From the Workspace Configuration Panel (), click Unit Cell
    • The unit cell disappears from the workspacethe 3D display area in the center of the main window, where molecular structures are displayed
  2. Click Unit Cell again
    • The unit cell reappears
    • Visualizing the unit cell is useful moving forward, so let’s keep this workspacethe 3D display area in the center of the main window, where molecular structures are displayed toggle selected

 

Note: Further display options for the unit cell box are available in the Preferences menu in Workspace > Periodicity

Figure 3-8. Toggling polyhedra on or off.

The polyhedral representation of the model can also be visualized.

  1.   From the Workspace Configuration Panel (), click Polyhedrons.
    • Polyhedra appear

For our purposes, we do not need to keep the polyhedra displayed.

  1. Click Polyhedrons again
    • The polyhedra disappear

Figure 3-9. Preparing a P1 cell after using the extents tool.

Now suppose we want to transform the original unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors to a new unit cellthe repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors .  There are three typical cases that might warrant doing so:

  • To create a larger cell (possibly of a different shape).
  • To create a cell representing the crystal system under stress or strain.
  • To create or manipulate a cell with a vacuum gap (for example for slabs or molecule in the box models)

Each case requires careful thought before redefining the lattice. The first case, creating a larger cell, is the most common transformation (for example, to prepare an appropriate sized unit cell to study adsorption of a small molecule) and is straightforward. For integer extent expansions, this can be done directly in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed.

  1. Click on the button to open the Periodic Structure Tool Window
  2. Select Make P1 Cell
    • The unit cell is now what was the 2 x 2 x 2 supercell displayed earlier
    • The bottom right indicates that this is now the ASU in the P1 space group
    • A new entry has been added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion titled alpha_quartz P1 and is includedthe entry is represented in the Workspace, the circle in the In column is blue by default

Note: Because we are now fundamentally altering the periodic boundary conditions associated with our original crystal, a new entry is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion (alpha_quartz P1) so that we also retain our original structure (alpha_quartz). In this case, the structures are fundamentally the same in the context of the infinite crystal, but they are defined by different sized unit cells

Figure 3-10. The Redefine Lattice panel.

Optional: To conclude this section, we will learn how to use the Redefine Lattice panel to transform the cell to represent a system under stress or strain. Feel free to skip to Section 4 if this topic is not relevant

  1. Go to Tasks > Materials > Tools > Redefine Lattice
    • The Redefine Lattice panel opens
    • So long as the alpha_quartz P1 entry is includedthe entry is represented in the Workspace, the circle in the In column is blue, the cell parameters are loaded into the panel

Suppose we wanted to prepare a crystal that is uniaxially strained along the a-axis

  1. Keep the Set new cell parameters radio button checked, adjust a to 10.02, and press enter (representing a hypothetical ~2% strain from the original a = 9.828 Å)
    • The Transformation Matrix automatically updates
  2. For Transformation mode, select Strain (fractional coordinates are fixed)
  3. Click Run
  4.  
    • A new entry is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion. It is both selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includedthe entry is represented in the Workspace, the circle in the In column is blue, and is named alpha_quartz P1

Figure 3-11. The P1 cell after the cell parameters have been redefined.

  1. Close the Redefine Lattice panel
  2. Rename this entry alpha_quartz P1_strained

Note: We have now fundamentally altered the periodic boundary conditions and cartesian coordinates associated with our original crystal. Doing so requires caution. An optimization is always suggested before performing any property calculations on a crystal structure. This is especially true after altering periodic boundary conditions

Here are a few additional general notes for the advanced user regarding the Redefine Lattice panel:

  • The transformation matrix can be adjusted manually. Note that the transformation matrix is the more general tool analogous to the extents tool.
  • Three transformation modes are available:
    • Regular: expects an integer valued transformation matrix with a positive integer determinant (D); creates a supercell with D times the number of original atoms and volume of the loaded cell parameters.
    • Strain: The number of atoms in the cell stays the same, and the fractional coordinates of the atoms stay the same. The translation vectors and coordinates are scaled proportionally to the transformation.
    • Frame: The number of atoms in the cell stays the same, and the cartesian coordinates of the atoms stay the same. The atoms stay fixed and the cell itself is transformed to generate a gap between the atoms and the cell vectors.

4. Manipulating an Organic Crystal Structure

All of the capabilities for handling inorganic periodic systems also apply to organic crystal structures. In this section we’ll briefly look at an organic crystal structure.

Figure 4-1. Select and include alpha glycine.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue alpha_glycine from the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
    • The alpha glycine crystal structure appears in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Figure 4-2. Updating style.

  1. In the Style menu, click ball-and-stick representation and Color Atoms
    • The atoms now appear as spheres rather than points
    • The coloring is updated
  2. From the Color Atoms dropdown, select Element
    • The atoms are colored by element

Figure 4-3. Building a cell.

As was the case for inorganic structures, the Periodic Structure Tool Window () allows visualization of the periodic system:

  1. Click on the button to open the Periodic Structure Tool Window
  2. Select Build Cell
  3. Select all four of the options
    • Within Translate to First Unit Cell, select Intact Molecules
    • For extents, add 1 for +A and +B to create a 2 x 2 x 1 supercell.
  4. Click Apply
    • A new supercell appears in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed

 

Note: Choosing Intact Molecules in the Translate to First Unit Cell section translates the molecules so that all of their centroids are in the unit cell, allowing visualization of complete molecules rather than fragments

Feel free to explore visualizing and manipulating the organic system analogously to our steps in Section 3.

5. Building a Crystal Structure with Crystal Builder

In some instances, you may wish to build your own crystal structure from scratch rather than importing an experimental structure. Let’s learn how to do so for a simple NaCl crystal structure. Building a crystal structure with the Crystal Builder panel requires that you define a space group and the fractional coordinates and occupancies of the atoms in the cell.

Figure 5-1. Reset the builder panel.

  1. Go to Tasks > Materials > Structure Builders > Crystal Structure
  2. If the panel is populated with data from one of the previous steps, simply click Reset in the bottom left corner of the panel
    • All inputs return to their default settings

Figure 5-2. Defining NaCl in the builder.

The space groupthe point group operations (reflections and rotations) and translations that define the unit cell, periodic boundary conditions and fractional coordinatesa coordinate system for the position of atoms in the basis of the unit cell vectors associated with NaCl are provided in the following steps. Populate the panel with the following information:

  1. Space Group Number: 225
    • F 4/m -3 2/m appears
  2. Lattice Parameters:
    • a = b = c = 5.64 Å
    • ɑ = β = ɣ = 90º
    • Some parameters are fixed, as dictated by the space group
  3. Within Asymmetric unit, click Add Atom twice
    • Two C atoms are added by default
  4. Click on the first C to update the first atom to Na and repeat to update the second atom to Cl
  5. The Na ion should have the fractional coordinates (0 0 0) and the Cl ion should have the fractional coordinatesa coordinate system for the position of atoms in the basis of the unit cell vectors (0.5 0 0). Update the f1, f2 and f3 entries accordingly as needed
    • Fractional coordinatesa coordinate system for the position of atoms in the basis of the unit cell vectors is defined in the Glossary of Terms
  6. Maintain the remaining defaults. The panel should match that which is shown in the Figure
  7. Click Create Crystal Structure
  8. Title it NaCl
    • A new entry is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and is includedthe entry is represented in the Workspace, the circle in the In column is blue in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed. Although only two atoms were specified (the asymmetric unit), eight atoms are generated (the full unit cell) due to the symmetry operators applied to the asymmetric unit
  9. Click OK
  10. Close the panel

Figure 5-3. Updating the style.

  1. In the Style () menu, click ball-and-stick representation
    • The atoms now appear as spheres rather than points
  2. If anything is not to your visual liking with respect to bonds and connectivity, refer to Section 3 and feel free to experiment with viewing and manipulating the NaCl structure

Figure 5-4. A snippet of some of the data available in the Project Table.

We may be interested in the volume or density of this cell. We also may also wish to see the lattice parametersthe lengths of the unit cell vectors (a, b, c) and angles between them (ɑ, β, ɣ) which define the unit cell for any crystals in the project all in one place

  1. Open the Project Table ()

Notice that for all of the structures, the unit cell density, unit cell volume, space group and lattice parameters are shown

Note: To add or remove columns from the Project Table, the Property Tree () is a useful tool

  1. Close the Project Table ()

6. Introducing Vacancies and Defects into a Crystal Structure

We might be interested in how properties of a solid state structure change with the introduction of vacancies or defects into the crystal structure. The first step would typically be optimization of the cell parameters and atomic coordinates of the perfect crystal at the quantum mechanical level (see the Electronic Structure Calculations of Bulk Crystals Using Quantum ESPRESSO tutorial). In this example, however, we will simply proceed with the unoptimized NaCl structure from Section 5 and introduce vacancies or defects.

Figure 6-1. Building a larger cell.

  1. With the NaCl entry selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includedthe entry is represented in the Workspace, the circle in the In column is blue, click on the button to open the Periodic Structure Tool Window
  2. Click Build Cell 
  3. Select Recalculate Connectivity
  4. Select Extents and input 1 in the +A, +B +C, -A, -B, and -C directions
  5. Check Create new entry and then click Apply.
    • A new entry is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion named NaCl by default
    • The entry is selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includedthe entry is represented in the Workspace, the circle in the In column is blue
    • The structure in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed now includes 216 atoms.

Figure 6-2. Renaming the entry.

  1. Rename the entry by double-clicking on the new NaCl in the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion and typing NaCl_vac

Figure 6-3. Duplicate the entry.

  1. Right click on NaCl_vac and click Duplicate > In Place
    • Another entry named NaCl_vac is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion.
  2. Rename this entry (as in Step 6): NaCl_dop
    • Your entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion should now includethe entry is represented in the Workspace, the circle in the In column is blue three entries related to the NaCl crystal:
      • The original unit cell with 8 atoms (NaCl)
      • The 3 x 3 x 3 cell with 216 atoms (NaCl_vac)
      • The same cell with 216 atoms (NaCl_dop)

Figure 6-4. Setting the selection scope to atoms.

First let’s introduce vacancies manually.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue NaCl_vac

Make sure your current selection scope is Atoms (it may still be Molecules depending on which previous steps you have completed):

  1. Set the selection scope to Atoms using the menu in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed or by typing a on your keyboard
    • The cursor updates to notify you that atom selection is active

Figure 6-5. Selecting and deleting a Na and Cl ion.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries one Na ion and one Cl ion in the lattice (using Shift + Click)
    • A selection box appears around the selected ions
  2. From the Build menu (in the Toolbar), click delete selected atoms
    • The selected ions are deleted from the structure. You have now created a Schottky defect

Note: In practice, from here you must create a new P1 cell before proceeding to calculate properties of the crystal structure with these vacancies.

Note: Vacancies can be introduced more systematically in the Crystal Builder panel (without deleting the original coordinates) by changing the occupancy of an atom to 0%

Figure 6-6. New P1 cell with a centered unit cell.

Let’s look at one final useful case, to introduce a substitutional dopant

  1. Now, select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue NaCl_dop

If you wanted to replace one Na ion with a K ion, you could use a similar process to the previous steps with the Build menu. In this case, let’s make this alteration slightly differently to demonstrate using the Crystal Builder panel:

  1. Click on the button to open the Periodic Structure Tool Window
  2. Select Make P1 Cell
    • A new entry is added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion named NaCl_dop P1.
  3. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue NaCl_dop P1.
  4. If the unit cell box is not properly centered, center the cell in the Workspace Configuration Panel ()

Figure 6-7. Substituting one Na ion for a K ion.

  1. Go to Tasks > Materials > Structure Builders > Crystal Structure
    • The Crystal Structure panel opens
  2. Click Load from Workspace
    • The asymmetric unit is populated with the atoms comprising the P1 cell from the workspacethe 3D display area in the center of the main window, where molecular structures are displayed
  3. Change one of the Na atoms to K by clicking on the Na element button in the table and using the Choose Element panel
  4. Click OK
  5. Click Create Crystal Structure
  6. Name the structure NaCl_dop_K
  7. Click OK
  8. Close the Crystal Structure panel

 

Figure 6-8. Selecting the K ion.

  1. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue Na_Cl_dop_K
  2. In the Style menu, click ball-and-stick representation
    • The atoms now appear as spheres rather than points

The K default color happens to be the same as Na. Let’s use a convenient selection method to locate the K ion and change its color to purple

  1. From the main menu, go to Select > Define
    • The atom selection panel opens
  2. Click Element, K and then Add
    • K appears in the ASL box
  3. Click OK
    • The potassium ion will now be selected in the crystal structure

Note: Selecting atoms via this method is quite useful when seeking particular selections in large structures

Figure 6-9. Coloring the K ion.

  1. In the Style menu, use the paint tool to color the selected atom purple
    • The K ion turns purple

Figure 6-10. The Assign Space Group panel.

After making a new P1 cell and manipulating the structure (in this case substituting one Na ion for a K ion), it may be practical to assign a space group to the new structure - the assignment of the space group is mainly useful for reducing the cell to a primitive cell.

  1. With NaCl_dop_K selected(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includedthe entry is represented in the Workspace, the circle in the In column is blue, go to Tasks > Materials > Tools > Assign Space Group
  2. Click Run
    • New entries are added to the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion which can be stylized again as needed
  3. Close the Assign Space Group panel
  4. Open the Project Table ()
    • You should notice that the new structure has Space Group ID 221

7. Handling Atomic Disorder

Finally, we’ll now learn how to handle an imported crystal structure with disorder. We will walk through an example for a zinc complex. Optionally, there are two additional .cif files included in the tutorial that can be imported for extra practice with handling typical disorder: an organic crystal (lovastatin.cif) and a MOF crystal (MOF-801.cif).

In nature, crystals are often characterized by atomic disorder. There is a form of disorder of specific interest, where a crystal unit cell is still well-defined but:

  1. Some lattice sites may be occupied by different atomic species with certain probability;
  2. Two or more symmetry related lattice sites exist, only one of which may be occupied by an atom. In the crystal structure (cif) files such a structure would be typically overspecified having both possible sites listed, each with fractional occupancy. For example, consider an H atom that may occupy one of the two atomic sites: just above or just below the mirror plane.      

In such cases, it is customary to characterize the crystal with its full space symmetry, but assigning some lattice sites fractional occupancies. We note that in atomic modeling, we must use structures with unique positions of atoms with complete occupancy. So, in order to consider a disorder, we need to generate an ensemble of all possible disordered structures (or at least their representative ensemble). Luckily, in many cases various realizations of disordered structures are equivalent, so only one or very few disordered structures are necessary to describe the whole ensemble. The following steps allow us to generate such representative structures. For subsequent computation, fully occupied atom sites are necessary.

Figure 7-1. Import, rename, select and include the zinc complex.

  1. Go to File > Import Structures
    • Import panel opens
  2. Select Section_07 > 1510909.cif
  3. Click Open
    • The crystal structure is imported
  4. Rename the entry:
    • Cm48usa → Zn_Disorder
  5. Select(1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries and includethe entry is represented in the Workspace, the circle in the In column is blue Zn_Disorder from the entry lista simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion
    • The Zn_Disorder crystal structure appears in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed

Figure 7-2. Select, include and stylize Zn_Disorder.

  1. In the Style menu, click ball-and-stick representation and Color Atoms
    • The atoms now appear as spheres rather than points
    • The coloring is updated

Notice that the oxygen atoms on the perchlorate ligand and anion were modeled with disorder. Also notice that the hydrogen atoms are not necessarily bonded by default

 

Figure 7-3. Building the unit cell for the Zn system.

  1. Click on the button to open the Periodic Structure Tool Window
  2. Select Build Cell > Recalculate Connectivity and Recalculate Bond Orders
  3. Click Apply
    • An error appears “Atoms with disorder group found”
    • Click OK
  4. Click Open Crystal Builder

Figure 7-4. Resolving the disorder.

If you scroll through the atoms in the asymmetric unit, you will notice that some of the oxygen atoms have partial occupancies (< 1.0)

  1. Click Create Crystal Structure
    • A drop-down appears for handling the disorder
  2. Select which group of atoms you would like to include with full occupancy, and in doing so, which group will be erased. A radio button appears alongside either group with which to make the selection. In this case, let’s choose the top radio button because the occupancies are a bit higher
  3. Click OK
  4. Again click Create Crystal Structure
    • A warning appears that the atoms in the excluded group will now have zero occupancy
  5. Click OK
  6. Rename the entry Zn_Disorder_occ
  7. Click OK
  8. Close the panel
    • The crystal structure includedthe entry is represented in the Workspace, the circle in the In column is blue in the workspacethe 3D display area in the center of the main window, where molecular structures are displayed (Zn_Disorder_occ) now only contains one set of oxygen atoms at the disordered perchlorate sites and can now be stylized, visualized and used for subsequent calculations

8. Conclusion and References

In this tutorial, we learned how to import, build and manipulate crystal structures in MS Maestro. These steps are elemental for proceeding with any calculations on bulk crystalline materials.

For further learning:

For introductory content, focused on navigating the Schrödinger Materials Science interface, an Introduction to Materials Science Maestro tutorial is available. Please visit the materials science training website for access to 70+ tutorials. For scientific inquiries or technical troubleshooting, submit a ticket to our Technical Support Scientists at help@schrodinger.com.

For self-paced, asynchronous, online courses in Materials Science modeling, including access to Schrödinger software, please visit the Schrödinger Online Learning portal on our website.

For some related practice, proceed to explore other relevant tutorials:

For further reading:

9. Glossary of Terms

Unit Cell - the repeating pattern in the material having the full symmetry of the crystal structure, which can be infinitely repeated by translation along unit cell vectors a, b, and c

Primitive Unit Cell - the smallest unit cell

Asymmetric Unit - the smallest asymmetric subset of atomic sites in the unit cell

Fractional Coordinates - a coordinate system for the position of atoms in the basis of the unit cell vectors

Lattice Parameters - the lengths of the unit cell vectors (a, b, c) and angles between them (ɑ, β, ɣ) which define the unit cell

Space Group - the point group operations (reflections and rotations) and translations that define the unit cell

Entry List - a simplified view of the Project Table that allows you to perform basic operations such as selection and inclusion

Included - the entry is represented in the Workspace, the circle in the In column is blue

Project Table - displays the contents of a project and is also an interface for performing operations on selected entries, viewing properties, and organizing structures and data

Recent actions - This is a list of your recent actions, which you can use to reopen a panel, displayed below the Browse row. (Right-click to delete.)

Scratch Project - a temporary project in which work is not saved, closing a scratch project removes all current work and begins a new scratch project

Selected - (1) the atoms are chosen in the Workspace. These atoms are referred to as "the selection" or "the atom selection". Workspace operations are performed on the selected atoms. (2) The entry is chosen in the Entry List (and Project Table) and the row for the entry is highlighted. Project operations are performed on all selected entries

Working Directory - the location where files are saved

Workspace - the 3D display area in the center of the main window, where molecular structures are displayed