【答案】应助回帖

选择科学就要耐得住寂寞！

11楼2012-03-18 14:34:28

chenyuchen

木虫 (正式写手)

应助: 32 (小学生)
贵宾: 0.005
金币: 4590.1
散金: 102
红花: 3
帖子: 832
在线: 208.4小时
虫号: 406410
注册: 2007-06-18
性别: GG
专业: 凝聚态物性 II ：电子结构

★
fzx2008: 金币+1, 谢谢回帖交流！ 2012-03-18 14:38:37

引用回帖:

10楼: Originally posted by fanarsenal at 2012-03-17 22:22:42:
这是涉及激发态的计算，好像只要是这样的都和wien2k有关，求科普。wien2k和有什么不同啊……

wien2k是全电子而非赝势方法，所以激发态计算更准一些

选择科学就要耐得住寂寞！

12楼2012-03-18 14:36:49

锐利的碎片

木虫 (正式写手)

star watcher

1ST强帖: 1
应助: 136 (高中生)
金币: 3637.1
散金: 252
红花: 22
帖子: 988
在线: 1224.9小时
虫号: 961933
注册: 2010-03-05
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 谢谢指教 2012-03-18 19:27:18

我记得wien2k不能算激发态。楼主不用纠结动量矩阵，那论文里也没有分析只是提了下，实际这个没法分析。因为矩阵非常大。

13楼2012-03-18 18:59:40

fanarsenal

木虫 (正式写手)

应助: 47 (小学生)
金币: 2021.7
红花: 23
帖子: 515
在线: 659.4小时
虫号: 1069215
注册: 2010-08-04
专业: 金属结构材料

★ ★
franch: 金币+2, 鼓励交流。。。。。呵呵 2012-03-19 09:40:09

引用回帖:

12楼: Originally posted by chenyuchen at 2012-03-18 14:36:49:
wien2k是全电子而非赝势方法，所以激发态计算更准一些

谢谢哈……了解了……
全电子对激发态应该也只是比赝势好些而已吧，激发态还是应该用格林函数吧？

14楼2012-03-18 22:21:25

bananawuzhu

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 384.4
散金: 9
帖子: 27
在线: 64.8小时
虫号: 892507
注册: 2009-11-03
专业: 凝聚态物性 II ：电子结构

引用回帖:

13楼: Originally posted by 锐利的碎片 at 2012-03-18 18:59:40:
我记得wien2k不能算激发态。楼主不用纠结动量矩阵，那论文里也没有分析只是提了下，实际这个没法分析。因为矩阵非常大。

谢谢您的回复！
之前是觉得如果可以定量算出相应数值能更好的进行分析。

15楼2012-03-19 09:19:19

chwipoc

木虫 (小有名气)

应助: 4 (幼儿园)
金币: 2214.1
散金: 518
帖子: 148
在线: 172.2小时
虫号: 1669095
注册: 2012-03-06
性别: GG
专业: 计算机体系结构

果断求指导，目前打算做GaAs纳米线有关光学性质的这一块，目前有些迷茫

没有

16楼2012-11-06 14:58:56

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

\int_{k}dk<\psi_k^v|p_i|\psi_k^c>^2*\delta(E_{\psi_k^c}-E_{\psi_k^v}-\hbar\omega)

<\psi_k^v|p_i|\psi_k^c> 这一项对应的就是跃迁矩阵，应该决定了电子跃迁时的动量守恒
(E_{\psi_k^c}-E_{\psi_k^v}-\hbar\omega) 这一项是能量守恒。

5.2版本的确没有动量矩阵元的输出，在4.6版本中有输出文件为MME momentum matrix element

但是5.2版本有介电函数的输出，有了介电函数，线性光学性质就可以代入计算了。

但如果要计算更高阶的光性，则需要动量矩阵元的单独输出。

楼主刚学习，可以对照4.6
尝试修改vasp的main.f程序，以及相应的光学计算程序，比如optic.f，编译后计算，估计能单独得到MME。
本人编程很差，尝试过几次，未成功。

17楼2012-11-08 09:47:25

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

14楼: Originally posted by fanarsenal at 2012-03-18 22:21:25
谢谢哈……了解了……
全电子对激发态应该也只是比赝势好些而已吧，激发态还是应该用格林函数吧？...

激发态的计算是需要格林函数，
格林函数G描述的只是准粒子的运动特性，
准粒子的物理作用用动态屏蔽库伦势W表速
iGW是自能，替代DFT里的Exc，则可以计算激发态。

18楼2012-11-08 09:49:56

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

10楼: Originally posted by fanarsenal at 2012-03-17 22:22:42
这是涉及激发态的计算，好像只要是这样的都和wien2k有关，求科普。wien2k和有什么不同啊……

光学计算不一定必须考虑激发态，但是考虑激发态计算结果更可靠。
激发态的考虑，则需要计算双准粒子的BSE方程，目前vasp还不能完成。

vasp目前激发态的计算也只能用GW近似算下能带隙（个人见解）

19楼2012-11-08 09:55:13

707462558

银虫 (著名写手)

应助: 2 (幼儿园)
金币: 358
散金: 41
红花: 37
帖子: 1096
在线: 431.8小时
虫号: 1422788
注册: 2011-09-29
性别: MM
专业: 凝聚态物性 II ：电子结构

引用回帖:

17楼: Originally posted by ApolloYang at 2012-11-08 09:47:25
\int_{k}dk<\psi_k^v|p_i|\psi_k^c>^2*\delta(E_{\psi_k^c}-E_{\psi_k^v}-\hbar\omega)

<\psi_k^v|p_i|\psi_k^c> 这一项对应的就是跃迁矩阵，应该决定了电子跃迁时的动量守恒
(E_{\psi_k^c}-E_{\p ...

请问一个介电函数矩阵怎么样才能带进去这个公式里呀

ND5T@`T9`HVF{`K%]FPTYVG.jpg

无论生活和你开了多大的玩笑我们都得继续生活

20楼2013-03-10 14:53:48