原子軌道線性組合

本條目中，向量與標量分別用粗體與斜體顯示。例如，位置向量通常用 $\mathbf{r}\,\!$ 表示；而其大小則用 $r\,\!$ 來表示。

原子軌域線性組合（Linear combination of atomic orbitals，或者簡寫為LCAO），是量子化學中用於求解分子軌域的一種方法，這種方法是通過對原子軌域進行線性疊加來構造分子軌域。因為它屬於分子軌域方法的一種，所以又稱原子軌域線性組合的分子軌域方法，或者叫LCAO-MO。它於1929年由約翰·蘭納-瓊斯爵士引入用於描述元素周期表第一行上原子構成的雙原子分子的成鍵，並且經由Ugo Fano進行了擴展。

在量子力學裡，原子的電子組態由波函數來描述。從數學上來看，這些波函數構成了函數基組。在化學反應過程中，軌道波函數會發生改變，根據原子所參與形成的化學鍵的類型，電子雲的形狀會相應改變。

LCAO的數學形式為：

$$\ \Psi_i = \sum_{j}^n c_{ji} \varphi_j$$

其中$\Psi_i$為第$i$條分子軌道，它被表示為$n$個原子基函數（原子軌道）$\varphi_j$的線性疊加。係數$c_{ji}$表示了第$j$條原子軌道對該分子軌道$i$的貢獻大小。

作為基函數的原子軌道$\varphi_j$通常是在（核）中心場作用下的單電子波函數。所使用的基函數通常是類氫原子，因為類氫原子波函數已知有解析的表達式。當然，基函數也可以選擇如高斯函數的其他形式。

通過變分法求系統總能量的最低值，人們可以獲得線性展開式前每項的係數$c_{ji}$。這種定量方法稱為Hartee-Fock方法。但隨著計算化學的發展，人們一般不用LCAO做波函數的實際優化，只用其作定性估測，以衡量或預測其他計算方法的結果。

## 基本計算過程

假設分子系統的哈密頓量為$\hat H$，其定態薛丁格方程為$\hat H\Psi = E\Psi$。

其中$\Psi$為分子軌道（分子波函數），$E$分子體系的能量。

LCAO的基本思想就是用原子軌道$\varphi$的線性組合來表示分子軌道$\Psi$：

$$|\Psi\rangle = \sum\limits_k {{c_k}\left| {{\varphi _k}} \right\rangle }$$

將其代入到定態薛丁格方程中。

$$\sum\limits_k {{c_k}\hat H\left| {{\varphi _k}} \right\rangle } = E\sum\limits_k {{c_k}\left| {{\varphi _k}} \right\rangle }$$

$$\sum\limits_k {{c_k}\underbrace {\left\langle {\varphi _i} \right|\hat H\left| \varphi _k \right\rangle }_{H_{ik}}} = E\sum\limits_k {{c_k}\underbrace {\left\langle {\varphi _i} | {\varphi _k} \right\rangle }_{S_{ik}}}$$

$$<pre> \sum\limits_k {{c_k}\left( {{H_{ik}} - E{S_{ik}}} \right) = 0} </pre>$$

所得到的線性方程組系統為久期方程。注意，在LCAO中，$\left\langle {\varphi _i} | {\varphi _k} \right\rangle \ne \delta_{i,k}$，這是因為這裡的$i,k$代表的不再是同一原子的波函數，而是處於不同位置的原子的波函數，它們一般不滿足正交歸一性。$S_{ik}$與原子間的位置相關，原子間相距近，則波函數間交疊大；若原子相距很遠，$S_{ik}$則趨於零，因此$S_{ik}$被稱作重疊積分（overlap integral）。

## 實例：雙原子分子

記雙原子分子中兩個原子的波函數分別為$\varphi_A$與$\varphi_B$，根據LCAO，分子波函數可以寫作線性組合：

$$\Psi = {c_A}{\varphi _A} + {c_B}{\varphi _B}$$

代入到定態薛丁格方程$\hat H\Psi = E\Psi$中。

$$\hat H\left( {{c_A}{\varphi _A} + {c_B}{\varphi _B}} \right) = E\left( {{c_A}{\varphi _A} + {c_B}{\varphi _B}} \right)$$

分別用兩個原子波函數與上式做內積。

$$\int {d\tau \;\varphi _A^*\hat H\left( {{c_A}{\varphi _A} + {c_B}{\varphi _B}} \right)} = E\int {d\tau \;\left( {{c_A}\varphi _A^*{\varphi _A} + {c_B}\varphi _A^*{\varphi _B}} \right)}$$

$$\int {d\tau \;\varphi _B^*\hat H\left( {{c_A}{\varphi _A} + {c_B}{\varphi _B}} \right)} = E\int {d\tau \;\left( {{c_A}\varphi _B^*{\varphi _A} + {c_B}\varphi _B^*{\varphi _B}} \right)}$$

展開。

$${c_A}\underbrace {\int {d\tau \;\varphi _A^*\hat H{\varphi _A}} }_{{H_{AA}}} + {c_B}\underbrace {\int {d\tau \;\varphi _A^*\hat H{\varphi _B}} }_{{H_{AB}}} = E{c_A}\underbrace {\int {d\tau \;\varphi _A^*{\varphi _A}} }_1 + E{c_B}\underbrace {\int {d\tau \;\varphi _A^*{\varphi _B}} }_{{S_{AB}}}$$

$${c_A}\underbrace {\int {d\tau \;\varphi _B^*\hat H{\varphi _A}} }_{{H_{BA}}} + {c_B}\underbrace {\int {d\tau \;\varphi _B^*\hat H{\varphi _B}} }_{{H_{BB}}} = E{c_A}\underbrace {\int {d\tau \;\varphi _B^*{\varphi _A}} }_{{S_{BA}}} + E{c_B}\underbrace {\int {d\tau \;\varphi _B^*{\varphi _B}} }_1$$

因此得到。

$${c_A}\left( {{H_{AA}} - E} \right) + {c_B}\left( {{H_{AB}} - E{S_{AB}}} \right) = 0$$

$${c_A}\left( {{H_{BA}} - E{S_{BA}}} \right) + {c_B}\left( {{H_{BB}} - E} \right) = 0$$

相應的久期方程矩陣形式為

$$\begin{bmatrix} <pre> {{H_{AA}} - E} & {{H_{AB}} - E{S_{AB}}} \\ {{H_{BA}} - E{S_{BA}}} & {{H_{BB}} - E} \\ </pre> \end{bmatrix} \begin{bmatrix} <pre> {c_A} \\ {c_B} \\ </pre> \end{bmatrix} = 0$$

線性組合的係數由此可求得。

雙原子分子體系的能量$E$可由兩個方程之比求得。

$$\frac{{{H_{AA}} - E}}{{{H_{BA}} - E{S_{BA}}}} = \frac{{{H_{AB}} - E{S_{AB}}}}{{{H_{BB}} - E}}$$

### 最簡單的分子： H$_2^+$

H$_2^+$是由兩個質子與一個電子組成的同核雙原子分子，是最簡單的分子形式。設想H$_2^+$的分子軌道可以由兩個氫原子的基態波函數1s線性疊加而成。此時滿足${H_{AA}} = {H_{BB}} = \alpha ,{H_{AB}} = {H_{BA}} = \beta ,{S_{AB}} = {S_{BA}} = S$，其中α為庫侖積分，β為交換積分，S為重疊積分。於是，代入用於求能量的比值式：

$$\frac{{\alpha - E}}{{\beta - ES}} = \frac{{\beta - ES}}{{\alpha - E}}$$

可得到兩個可能的能量值；回代入久期方程，可得到係數$c_A$與$c_B$的關係。

${E_ + } = \frac{{\alpha + \beta }}{{1 + S}}$，此時有$c_A=c_B$

${E_ - } = \frac{{\alpha - \beta }}{{1 - S}}$，此時有$c_A=-c_B$

因此，令$c_A=c$，可得到兩個分子軌道

$${\Psi _ + } = c\left( {{\varphi _A} + {\varphi _B}} \right)$$

$${\Psi _ - } = c\left( {{\varphi _A} - {\varphi _B}} \right)$$

c可由歸一化條件最終確定。

已知氫原子基態波函數（1s）在空間中表示為$e^{-\frac{\mathbf{r}}{a_0}}$，考慮二維情況$\mathbf{r}=(x,y)$，設一個處於$x=0$處的氫原子基態波函數為$\varphi_A(\mathbf{r})=e^{-\frac{\sqrt{x^2+y^2}}{a_0}}$，另一個處於$x=x_0$處的氫原子基態波函數為$\varphi_B(\mathbf{r})=e^{-\frac{\sqrt{(x-x_0)^2+y^2}}{a_0}}$，對波函數按上面得到的分子軌道表達式進行線性疊加可得。

$${\Psi _ + }(x,y) = c\left( e^{-\frac{\sqrt{x^2+y^2}}{a_0}} + e^{-\frac{\sqrt{(x-x_0)^2+y^2}}{a_0}} \right)$$

$${\Psi _ - }(x,y) = c\left( e^{-\frac{\sqrt{x^2+y^2}}{a_0}} - e^{-\frac{\sqrt{(x-x_0)^2+y^2}}{a_0}} \right)$$