www.毛片在线观看_怡红院亚洲红怡院天堂麻豆_欧美jizz18欧美_国产精品免费久久久久影院

Document

當前位置：文庫百科 ? 文章詳情

Gaussian量化模擬入門教程（三）：實例操作之自旋多重度的判斷

來源：本站時間：2020-11-04 09:59:03 瀏覽：17645次

注：在未特別說明的情況下，本文中涉及到的所有量化模擬實例均是在本地電腦（Windows系統(tǒng)）建模并生成輸入文件，在服務器（Linux系統(tǒng)）中進行計算，最后在本地電腦進行數(shù)據(jù)處理和分析。

在修改Gaussian輸入文件時，自旋多重度往往是困擾入門者的難題。正如前文所述，其計算公式為：2S+1= Nα-Nβ+1，其中S為自旋量子數(shù)，Nα為alpha（即自旋向上）電子數(shù)，Nβ為beta（即自旋向下）電子數(shù)。換言之，只要能夠知道模型分子的電子占據(jù)情況，即可算出該體系的自旋多重度。

簡單體系

1、對于比較簡單的體系，可根據(jù)原子或分子軌道理論來判斷。下面通過一些例子來進行解釋。

① 氫原子H
氫原子H核外電子排布為1s1，顯然只有1個alpha電子，無beta電子。因此其自旋多重度Nα-Nβ+1=2。
② 氧原子O
  氧原子O核外電子排布為[O] 2s22p4，其價電子共4個alpha電子，2個beta電子，因此其自旋多重度Nα-Nβ+1=3。
③ 氯原子Cl
  氯原子Cl核外電子排布為[Ne]3s23p5，其價電子共4個alpha電子，3個beta電子，因此其自旋多重度Nα-Nβ+1=2。
④ 鐵原子Fe
  鐵原子Fe核外電子排布為[Ar]3d64s2，其價電子共6個alpha電子，2個beta電子，因此其自旋多重度Nα-Nβ+1=5。
⑤ 二價鐵離子Fe2+
  二價鐵離子Fe2+核外電子排布為[Ar]3d6，其價電子共5個alpha電子，1個beta電子，因此其自旋多重度Nα-Nβ+1=5。
⑥ 三價鐵離子Fe3+
  三價鐵離子Fe3+核外電子排布為[Ar]3d5，其價電子共5個alpha電子，無beta電子，因此其自旋多重度Nα-Nβ+1=6。

⑦ 三價鐵離子Fe3+的其他自旋多重度

三價鐵離子Fe3+核外電子排布為[Ar]3d5。由于3d軌道共有5個伸展方向，根據(jù)洪特規(guī)則，當5個電子以自旋相同的方式分別占據(jù)不同軌道時，其總能量最低，自旋多重度為6，對應于Fe3+的基態(tài)。但倘若不考慮能量最低的穩(wěn)態(tài)結構，F(xiàn)e3+還會存在另外兩種電子排布形式，如下圖所示，其自旋多重度分別為4和2。根據(jù)洪特規(guī)則，這兩種電子結構的能量都應該比基態(tài)高，可以被視作Fe3+的某種激發(fā)態(tài)。據(jù)此，我們可以得知Fe3+所有可能的自旋多重度僅有2、4和6三種可能，此時自旋多重度=2所對應的電子結構被稱作低自旋態(tài)，而自旋多重度=6所對應的電子結構被稱作高自旋態(tài)。

事實上，量子化學模擬所計算的多是基態(tài)結構。因此，自旋多重度的設置是量化計算中非常重要的一步。若自旋多重度設定不合理，最終計算得到的結果必然不會是基態(tài)結構的相關信息，實例⑧則進一步展示了這一點。

⑧ 氧氣分子O2

對于分子而言，需要考慮分子軌道及電子占據(jù)情況。以O2分子為例，當兩個氧原子的原子軌道組合形成分子軌道時，其軌道及電子占據(jù)情況如下（內層軌道及電子省略），顯然可以計算出Nα-Nβ+1=3，此即為三線態(tài)氧，也就是通常所說的氧氣分子。

但分子氧中π2p*軌道上的電子還可以以另外兩種形式排布，此時的自旋多重度均為Nα-Nβ+1=1，即單線態(tài)氧。根據(jù)洪特規(guī)則，顯然三線態(tài)氧的電子排布所帶來的總能量要低于兩種單線態(tài)氧。因此，在進行基態(tài)的氧氣分子計算時，其自旋多重度應設定為3。

大對數(shù)體系

2、對于難以通過軌道理論進行判斷的體系，可以采取計算不同自旋多重度下的單點能，找出能量最低時對應的自旋多重度，即為基態(tài)電子結構。若只是為了找到體系基態(tài)的自旋多重度，對于大多數(shù)體系均可適用。

⑨ 尋找基態(tài)Fe3+的自旋多重度

倘若我們要對一個Fe3+進行模擬，但并不知道自旋多重度應取何值，且不擅長或無法通過軌道理論進行判斷，那么可以通過計算Fe3+在不同自旋多重度下的單點能的方式來尋找其基態(tài)自旋多重度。具體來說，由于Fe3+總電子數(shù)為奇數(shù)（共計23個電子），因此其自旋多重度必然是偶數(shù)值。既然如此，我們分別設定自旋多重度為2、4、6、8、10等等進行多次單點能計算。對比最終得到的能量可以發(fā)現(xiàn)，當自旋多重度取6時Fe3+能量最低，表明此時為其基態(tài)。

倘若自旋多重度設定值不合理，計算很可能會失敗。例如若Fe3+自旋多重度設定為7，在提交任務開始計算后很快就會報錯并中斷。從相應輸出文件末尾的錯誤原因提示可以看到，對于電子總數(shù)為23的Fe3+而言，自旋多重度為7是不可能出現(xiàn)的情況。此時應回過頭來重新修改輸入文件，解決問題后重新進行計算。

⑩ 尋找基態(tài)O2的自旋多重度

在實例⑧中我們展示了分子氧的三種電子結構。若不確定哪種才是基態(tài)，也可以分別計算三線態(tài)和單線態(tài)氧的能量進行對比。從下圖所示的結果可以看到，三線態(tài)氧的能量更低，因此判斷自旋多重度為3時的O2為基態(tài)分子氧。

Gaussian量化模擬入門教程（三）：實例操作之自旋多重度的判斷+參考圖6

Gaussian量化模擬入門教程（三）：實例操作之自旋多重度的判斷+參考圖7

在這個模擬中有一點必須說明，根據(jù)實例⑧中的分子氧電子結構示意圖可以判斷存在兩種單線態(tài)氧，其自旋多重度均為1。若想通過Gaussian輸入文件區(qū)分這種自旋多重度重合的電子結構，通常需要做激發(fā)態(tài)計算。常用的方法主要為CIS和TDDFT方法，可求解出在同一自旋多重度下的多個電子結構。由于本例中并未計算激發(fā)態(tài)，此處所得的是所有單線態(tài)氧中能量最低的情況，即1Δg型單線態(tài)氧。1Σg*型單線態(tài)氧屬于較為特殊的非限制性開殼層單重態(tài)，其能量要略高于1Δg型單線態(tài)氧，在本例中并未作為結果被輸出。

即便如此，上述事實也并不影響我們尋找基態(tài)分子氧的自旋多重度。因為基態(tài)的能量必然低于激發(fā)態(tài)，因此不管對于何種體系，只要找到能量最低的態(tài)即可斷定為基態(tài)。

上一篇：XRD精修教程（二）—CMPR軟件介紹

下一篇：Gaussian量化模擬入門教程（四）：實例操作之多巴胺分子的結構優(yōu)化及基本性質分析