| 查看: 1764 | 回复: 5 | |||||
[求助]
HF方程 已有1人参与
|
| 推导闭壳层分子体系的HF方程 |
回复此楼» 猜你喜欢
西南某211招27届博士
已经有33人回复
-
大连工业杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有99人回复
-
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有16人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
深圳大学高等研究院赵伟课题组招收化学调剂生
已经有3人回复
-
青岛大学固态电池研究团队招收材料、化学、物理等方向博士研究生
已经有29人回复
-
070300总分298求调剂
已经有1人回复
-
化学工程085602 305分求调剂
已经有17人回复
-
一志愿厦门大学化学学硕307求调剂
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3514.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★ ★ ★ ★
gmy1990: 金币+5 2014-05-29 05:45:56
gmy1990: 金币+5 2014-05-29 05:45:56
|
最简单的闭壳层模型就是两个电子，方便起见以这个为例子， 其中Hamilton量H=T+V(Ne)+V(ee)，动能算符T=-1/2▽^2，核对电子吸引V(Ne)=ΣZ/r1，电子-电子排斥V(ee)=1/r12 根据Pauli原理，体系波函数ψ写成Slater行列式，两个电子1和2有两个轨道i和j，分子体系的H不包含自旋变量并不考虑spin-orbital的耦合，可将自旋部分分离不考虑，ψ=1/√2[i(1)j(2)-i(2)j(1)]代入，方便起见先消去1/√2这个系数。 方便起见，把-1/2▽^2+ΣZ/r1写成hi算符，<[i(1)j(2)-i(2)j(1)]|h1+h2+V(ee)|[i(1)j(2)-i(2)j(1)]>-E<[i(1)j(2)-i(2)j(1)]|[i(1)j(2)-i(2)j(1)]>=0 h1对r1微分，根据i(1)j(1)正交及i(2)j(2)正交，<[i(1)j(2)-i(2)j(1)]|h1|[i(1)j(2)-i(2)j(1)]>=<i(1)|h1i(1)>+<j(1)|h(1)|j(1)> h2对r2微分，根据i(1)j(1)正交及i(2)j(2)正交，<[i(1)j(2)-i(2)j(1)]|h2|[i(1)j(2)-i(2)j(1)]>=<i(2)|h2i(2)>+<j(2)|h(2)|j(2)> <ψ|Vee|ψ>=<[i(1)j(2)-i(2)j(1)]|1/r12|[i(1)j(2)-i(2)j(1)]>=[<i(1)j(2)|1/r12|i(1)j(2)>+<i(2)j(1)|1/r12|i(2)j(1)>]-[<i(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|i(1)j(2)>] 令上式第一个[]=J=∫ρi(1)ρj(2)/r12dr1dr2+∫ρi(2)ρj(1)/r12dr1dr2，第二个[]=K=<i(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|i(1)j(2)> 对上式做i、j归一化约束条件的变分，构造泛函Ω=<i(1)|h1i(1)>+<j(1)|h(1)|j(1)>+<i(2)|h2i(2)>+<j(2)|h(2)|j(2)>+J-K-Ei<i(1)|i(1)>-Ej<j(1)|j(1)>-Ei<i(2)|i(2)>-Ej<j(2)|j(2)，其中Ei和Ej为Lagrange乘子。 方便起见，仅挑出对i(1)变分相关的项（增量i(1)无影响的项，变分中会被约去） Ω'=<i(1)|h1i(1)>+<i(1)j(2)|1/r12|i(1)j(2)>-[<i(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|i(1)j(2)>]-Ei<i(1)|i(1)>=0 对i(1)做作增量δi(1)，则Ω‘的增量为 δΩ'=<δi(1)|h1i(1)>+<i(1)|h1δi(1)]>+<δi(1)j(2)|1/r12|i(1)j(2)>+<i(1)j(2)|r12|δi(1)j(2)>+[<δi(1)j(2)|1/r12|i(2)j(1)>+<i(2)j(1)|1/r12|δi(1)j(2)>]-Ei<δi(1)|i(1)>-Ei<i(1)|δi(1)> 当δΩ'=0时，H为极值，挑出复共轭项<δi(1)| <δi(1)|h1i(1)>+<δi(1)j(2)|1/r12|i(1)j(2)>+<δi(1)j(2)|1/r12|i(2)j(1)>-Ei<δi(1)|i(1)>=0 由于变分过程中δi(1)的任意性，h1i(1)+<j(2)|1/r12|j(2)>i(1)+<j(2)|1/r12|i(2)>j(1)=Eii(1) 其中，第一项包含动能及核对电子吸引hi，第二项为电子-电子平均排斥Ji，第三项为交换势Ki，总体构成Fock算符Fi，右侧为i能级能量Ei。 打字打的比较费劲，其实直接参考一些量子力学和量子化学的书就可以了，上面公式是自己推导的，可能个别之处有不严谨甚至是出错的地方，由于本人数学及量子化学基础不太好，希望大家能够批评指正。 |
2楼2014-05-29 01:09:06
3楼2014-05-29 03:21:01
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3514.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2014-05-29 08:17:20
5楼2014-05-29 09:09:36
6楼2014-05-29 09:13:05
