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}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)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](https://chem.libretexts.org/@api/deki/files/469747/Figure2A.png?revision=1&size=bestfit&width=715&height=393)
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](https://chem.libretexts.org/@api/deki/files/469748/Figure2B.png?revision=1&size=bestfit&width=828&height=342)
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](https://chem.libretexts.org/@api/deki/files/469749/Figure3.png?revision=1)
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](https://chem.libretexts.org/@api/deki/files/469750/Figure4.png?revision=1&size=bestfit&width=451&height=202)
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](https://chem.libretexts.org/@api/deki/files/469754/Figure6A.png?revision=1&size=bestfit&width=705&height=199)
![Figure6B.png](https://chem.libretexts.org/@api/deki/files/469769/Figure6B.png?revision=1&size=bestfit&width=418&height=503)
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](https://chem.libretexts.org/@api/deki/files/469770/Figure7.png?revision=1&size=bestfit&width=1133&height=309)
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](https://chem.libretexts.org/@api/deki/files/469771/Figure8.png?revision=1&size=bestfit&width=857&height=580)
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
- Neese, F. The ORCA Program System. WIREs Computational Molecular Science 2012, 2 (1), 73–78. https://doi.org/10.1002/wcms.81.
- 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.
- 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.
- 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.
- Avogadro: An Open-Source Molecular Builder and Visualization Tool. http://avogadro.cc/.