| 查看: 1726 | 回复: 5 | |||||
[求助]
HF方程 已有1人参与
|
| 推导闭壳层分子体系的HF方程 |
» 猜你喜欢
量子点电致发光器件(QLED):蓝光量子点(QDs)和纳米氧化锌电子传输材料
已经有0人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有27人回复
物理化学论文润色/翻译怎么收费?
已经有163人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生(还有1个名额)
已经有0人回复
需要一个圆偏振光诱导化学反应,请问使用圆偏振光作为光源,应该需要哪些设备搭建
已经有0人回复
华北电力大学(北京)(第一性原理计算)博士招生——学博,专博各1人
已经有45人回复
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习 催化剂设计方向
已经有0人回复
法国里昂大学招生CSC 资助博士生
已经有0人回复
法国里昂大学接收CSC 资助访问学者
已经有0人回复
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生(计算优先)
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-

专家经验: +224 - 应助: 47 (小学生)
- 金币: 3492.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1253.7小时
- 虫号: 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 (小学生)
- 金币: 3492.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1253.7小时
- 虫号: 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












回复此楼