查看“︁科恩-沈吕九方程”︁的源代码

{{Refimprove|time=2012-06-25T07:14:31+00:00 }}
'''科恩－沈吕九方程'''（{{lang-en|'''Kohn–Sham equation'''}}，简称'''科恩－沈方程'''）在[[密度泛函理论]]里面指的是与真实体系相关的虚拟体系所满足的[[薛定谔方程]]。该虚拟体系中的粒子（通常是[[电子]]）在无相互作用的有效势场中运动，粒子密度在空间各点均与真实系统相同。<ref>
{{cite journal
 | last1 = Kohn | first1 = Walter
 | last2 = Sham | first2 = Lu Jeu
 | year = 1965
 | title = Self-Consistent Equations Including Exchange and Correlation Effects
 | journal = [[Physical Review]]
 | volume = 140 | pages = A1133–A1138 | issue = 4A
 | doi = 10.1103/PhysRev.140.A1133
|bibcode = 1965PhRv..140.1133K }}</ref><ref>
{{cite book
 | last1 = Parr | first1 = Robert G.
 | last2 = Yang | first2 = Weitao
 | year = 1994
 | title = Density-Functional Theory of Atoms and Molecules
 | publisher = [[Oxford University Press]]
 | isbn = 978-0-19-509276-9
}}</ref>科恩－沈吕九方程中的有效势通常用 <math>v_{\rm s}(\mathbf r)</math> 或 <math> v_{\rm eff}(\mathbf r)</math>) 来表示，称为科恩－沈势。虚拟系统中的粒子是彼此无相互作用的[[费米子]]，因此科恩－沈方程的精确解为单个斯莱特行列式，行列式中的轨道则称为科恩－沈轨道，每一个科恩－沈轨道都可以表示为原子轨道的线性组合，也可以按照基函数展开。科恩－沈方程的形式如下：

<center><math>\left(-\frac{\hbar^2}{2m}\nabla^2+v_{\rm eff}(\mathbf r)\right)\phi_{i}(\mathbf r)=\varepsilon_{i}\phi_{i}(\mathbf r)</math></center>

式中 <math>\varepsilon_i</math> 为科恩－沈轨道 <math>\phi_i</math> 的轨道能。含有 <math>N</math> 个粒子的科恩－沈系统的电子密度则由下式给出：

<center><math>\rho(\mathbf r)=\sum_i^N |\phi_{i}(\mathbf r)|^2.</math></center>

科恩－沈方程于 1965 年由[[加利福尼亚大学圣迭戈分校]]的[[沃尔特·科恩]]与[[沈吕九]]提出并以他们的名字命名。

==科恩－沈势==
[[密度泛函理论]]中，体系的能量是电子密度的[[泛函]]：

<center><math> E[\rho]  = T_s[\rho] + \int d\mathbf r\ v_{\rm ext}(\mathbf r)\rho(\mathbf r) + V_{H}[\rho] + E_{\rm xc}[\rho]</math></center>

式中 <math>T_s</math> 是科恩－沈[[动能]]项，可以用科恩－沈轨道表出如下：

<center><math>T_s[\rho]=\sum_{i=1}^N\int d\mathbf r\ \phi_i^*(\mathbf r)\left(-\frac{\hbar^2}{2m}\nabla^2\right)\phi_i(\mathbf r),</math></center>

<math> v_{\rm ext}</math> 是作用在真实系统上的外势（至少包括原子核与电子之间的相互作用势），<math> V_H </math> 是哈特里（库仑）能：

<center><math> V_{H}={e^2\over2}\int d\mathbf r\int d\mathbf{r}'\  {\rho(\mathbf r)\rho(\mathbf r')\over|\mathbf r-\mathbf r'|}.</math></center>

<math>E_{\rm xc}</math> 是交换相关能量。对虚拟体系总能量表达式右端除动能项之外的部分取电子密度的泛函微商<ref>{{Cite web |url=http://muchomas.lassp.cornell.edu/P480/Notes/dft/node11.html |title=存档副本 |access-date=2019-01-10 |archive-date=2020-02-18 |archive-url=https://web.archive.org/web/20200218032643/http://muchomas.lassp.cornell.edu/P480/Notes/dft/node11.html }}</ref>，就得到科恩－沈势的表达式：

<center><math>v_{\rm eff}(\mathbf r) = v_{\rm ext}(\mathbf{r}) + e^2\int {\rho(\mathbf{r}')\over|\mathbf r-\mathbf r'|}d\mathbf{r}' + {\delta E_{\rm xc}[\rho]\over\delta\rho(\mathbf r)}.</math></center>

上式中最后一项

<center><math>v_{\rm xc}(\mathbf r)\equiv{\delta E_{\rm xc}[\rho]\over\delta\rho(\mathbf r)}</math></center>

是交换相关势项。在整个密度泛函理论中只有这一项（及与之相关的能量）是未知的。

科恩－沈轨道能 <math>\varepsilon_i</math> 并没有明确的物理含义。它与体系总能量的关系由下式给出（参见[[库普曼斯定理]]）：

<center><math>E = \sum_{i}^N \varepsilon_i - V_{H}[\rho] + E_{\rm xc}[\rho] - \int {\delta E_{\rm xc}[\rho]\over\delta\rho(\mathbf r)} \rho(\mathbf{r}) d\mathbf{r}</math></center>

在限制性开壳层计算中，因为科恩－沈轨道的选取不唯一，上式仅对某些轨道能的选取成立。

== 参考文献 ==
{{reflist}}

[[Category:密度泛函理论]]