Skip to main content
Chemistry LibreTexts

12.3: Computational Instructions

  • Page ID
    470406
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    To calculate the energetics of the reaction we will need to calculate the starting materials and products for both the opening of oxetane and propylene oxide, as shown in Figure 2. While we have learned to calculate the transition states of substitution reactions in a previous computational exercise, the transition states for these reactions will be provided to help save time.

    Figure2A.png
    Figure 2A. Opening of propylene oxide with the hydrosulfide anion. B. Opening of oxetane with the hydrosulfide anion.

    Start by creating a folder on the desktop of your computer and label it as Ether Substitution. Within this folder, please create the following subfolders: SM_Epox, PR_Epox, SM_Oxet, and PR_Oxet. Next you should download the supporting files for this exercise. These files will include a generic input file (denoted by a .inp file extension), and molecular coordinates for the products and starting materials that you will run your calculations upon. Place a copy of the generic input file into each of the nested subfolders and place the appropriate coordinates file into each subfolder. A description of this file structure is shown in Figure 3.

    Figure2B.png
    Figure \(\PageIndex{1}\): Copy and Paste Caption here. (Copyright; author via source)

    The process for calculating the energy of the products or reactants of the oxetane or epoxide opening is the same for each calculation. To illustrate this process, we will walk through the calculation of the starting materials for the epoxide reaction. Begin by opening the starting coordinates file sm_epox_coord.xyz in Avogadro4,5 and make sure that it has the starting propylene oxide and hydrosulfide anion as shown in Figure 4.

    Figure3.png
    Figure 4. Expoxide opening reaction starting material viewed in Avogadro.

    You should then open the generic input file in notepad and save it as sm_epoxide.inp. As shown in Figure 5, please change the name of the file from coord.xyz to sm_epox_coord.xyz so that the computer knows to look for the epoxide starting material coordinates in your folder. Be sure to save your changes.

    Figure4.png
    Figure 5. Generic input script for determining the energy values of starting materials or products.

    We can now run our calculation using Orca via the command line as we did in previous exercises. Briefly, open the command prompt to your PC by right clicking on the start button and searching for command prompt. First, we need to tell the computer to look on the C drive and we do this by typing C: and hitting enter. Next, we need to tell the computer where the input script and the coordinates file are to run the calculation. We do this by typing cd (space) and pasting the file path. When you hit enter, the computer will paste a new line indicating that the current directory has changed, as shown in Figure 6A. To find the file path of your input script, right click on the input script (sm_epoxide.inp) and select properties. The file path will appear under location, and you can highlight and copy this file path (Figure 6B).

    Figure6A.png
    Figure6B.png
    Figure 6A. (Left) Changing of the file path in the command prompt to match the location of our input script. 6B. (Right) Locating the file path on the properties window of the input script (Red Arrow).

    Next, we will run the calculation by typing orca sm_epoxide.inp > sm_epoxide.out and pressing enter. After pressing enter it may not appear like much is happening, but the computer is working on your computation and depositing the results in the SM_Epox folder that you created. Depending upon the speed of your computer, the calculation will take about 1-3 minutes, and upon completion the command prompt will print another line indicating that it is ready for the next command (Figure 7).

    Figure7.png
    Figure 7. Running the calculation using the command line. The line indicated by the red arrow shows the computer that we want to use Orca to calculate the commands in sm_epoxide.inp and that the results of this calculation should be placed in the output file sm_epoxide.out. The line indicated by the blue arrow is the computer indicating that the calculation is complete, and the command prompt is ready for the next command.

    After the Orca job has completed you can access the energy values by opening the output file (sm_expoxide.out) in notepad. At the very end of the file (scroll to the bottom) will be the thermodynamic values that Orca calculated for the starting materials of the epoxide opening reaction. As shown in Figure 8, the value of Gibbs free energy in Hartree (Eh) can be found under the heading of Final Gibbs free energy.

    Figure8.png
    Figure 8. The Gibbs free energy of the starting material of the epoxide opening reaction (indicated by a red arrow).

    Please repeat this process for the products of the epoxide opening, the starting materials of the oxetane opening, and the products of the oxetane opening. After completing these computations use these results to complete the questions at the end of this assignment.

    References

    1. Neese, F. The ORCA Program System. WIREs Computational Molecular Science 2012, 2 (1), 73–78. https://doi.org/10.1002/wcms.81.
    2. Neese, F. Software Update: The ORCA Program System, Version 4.0. WIREs Computational Molecular Science 2018, 8 (1), e1327. https://doi.org/10.1002/wcms.1327.
    3. Neese, F.; Wennmohs, F.; Becker, U.; Riplinger, C. The ORCA Quantum Chemistry Program Package. J. Chem. Phys. 2020, 152 (22), 224108. https://doi.org/10.1063/5.0004608.
    4. Hanwell, M. D.; Curtis, D. E.; Lonie, D. C.; Vandermeersch, T.; Zurek, E.; Hutchison, G. R. Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J Cheminform 2012, 4 (1), 17. https://doi.org/10.1186/1758-2946-4-17.
    5. Avogadro: An Open-Source Molecular Builder and Visualization Tool. http://avogadro.cc/.

    This page titled 12.3: Computational Instructions is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Nicholas Boaz and Orion Pearce.

    • Was this article helpful?