计算比旋光度的基本流程:
- 首先优化分子结构
- 使用优化好的分子作为初始结构,用更大的基组进行比旋光度计算
- 计算结果分析
计算比旋光度的基本操作:
以甲基环氧乙烷作为例子,通过计算两种手性分子构型的比旋光度,并与实验值对比,以确定分子的绝对构型。实验数据我们选用气相条件下,在355nm和633nm下测得的比旋光度。
1、对映体的结构优化:
S构型甲基环氧乙烷结构优化的输入文件如下
添加图片注释,不超过 140 字(可选)
R构型甲基环氧乙烷结构优化的输入文件:
添加图片注释,不超过 140 字(可选)
2、对映体的比旋光度计算:
使用第一步优化得到的结构作为初始结构,计算S和R型甲基环氧乙烷的比旋光度。计算S和R型甲基环氧乙烷在气相条件下,355nm和633nm入射光下比旋光度的输入文件分别如下。
S构型甲基环氧乙烷比旋光度计算的输入文件:
添加图片注释,不超过 140 字(可选)
R构型甲基环氧乙烷比旋光度计算的输入文件:
添加图片注释,不超过 140 字(可选)
如果上一步优化步骤已经保存了check文件,比如opt.chk,则可以简化为下面的输入文件:
添加图片注释,不超过 140 字(可选)
3、Gaussview查看结果:
利用GaussView查看分子的比旋光度,我们以S构型甲基环氧乙烷比旋光度计算结果为例,步骤如下:
1、GaussView 6 | File | Open
打开文件:S-enantiomer-OR.log
2、GaussView 6 | Results | ORD
在Results下拉菜单中点击ORD选项,将显示不同入射波长下的比旋光度。用鼠标点击图中的点,在下方将显示指定波长下比旋光度。
添加图片注释,不超过 140 字(可选)
添加图片注释,不超过 140 字(可选)
4. 从Gaussian输出文件查看
除了用GaussView查看外,还可以直接从输出文件里查看,找[Alpha]行就可以,比如下面这个样子:
添加图片注释,不超过 140 字(可选)
其中,[Alpha] ( 5893.0 A) = -366.99 deg就是[α]589.30=-366.99°的意思。
5、结果分析:
利用同样的方法我们可以查看每个构型在不同入射频率下的的比旋光度,下表列出了所有计算和实验得到的比旋光度,可以看出在这个例子中,S构型计算得到的比旋光度与实验测得的更为接近,因此我们认为实验测量的分子为S构型,这也与事实一致。
添加图片注释,不超过 140 字(可选)
公众号“320科技工作室