| 查看: 5014 | 回复: 18 | ||
future_wl木虫 (著名写手)
|
[求助]
VASP计算介电张量的结果文件
|
|
计算完介电张量后（LPEAD=.TRUE LEPSILON=.TRUE), 在结果文件中搜索DIELECTRIC 发现两部分矩阵： HEAD OF MICROSCOPIC STATIC DIELECTRIC TENSOR (independent particle, excluding Hartree and local field effects) ------------------------------------------------------ 7.531256 0.000000 0.000000 0.000000 7.531256 0.000000 0.000000 0.000000 6.922124 ------------------------------------------------------ MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects in DFT) ------------------------------------------------------ 7.314593 0.000000 0.000000 0.000000 7.314593 0.000000 0.000000 0.000000 6.607187 ------------------------------------------------------ 请问这两部分分别表示什么？我如果要取介电张量的话应该用哪部分？ 谢谢大家！ |
回复此楼» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于G09计算的极化率的比较
已经有3人回复
- VASP光学性质计算步骤
已经有10人回复
- 如何用VASP计算体系的介电常数
已经有24人回复
- 请教压电原理的计算（有点类似电容）。
已经有7人回复
- 電磁場能動量?埩康挠嬎銌栴}，?V義相對論下的
已经有26人回复
- 【求助】用vasp算压电张量（系数），弹性张量需要注意的问题？
已经有9人回复
- 【求助】请问有人用WIEN2K算过光电导张量吗？
已经有4人回复
- 【求助】算VASP弹性常数时的能量-应变出错
已经有18人回复
- 【求助】关于Dielectric Tensor的问题
已经有5人回复
- 【求助】Pwscf现在可以算哪些力学性质？
已经有3人回复
- 【求助】求教VASP计算介电常数步骤
已经有16人回复
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-11-05 10:03:32
future_wl
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 3295.4
- 散金: 1881
- 红花: 5
- 帖子: 1404
- 在线: 676.6小时
- 虫号: 950891
- 注册: 2010-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2013-11-05 10:13:45
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
identation: 金币+2, 谢谢交流 2013-11-06 13:45:47
identation: 金币+2, 谢谢交流 2013-11-06 13:45:47
|
体系微观介电常数epsilon(r,r')是一个二变量的函数，对其作傅立叶变换后，傅立叶系数可以写成一个矩阵的形式epsilon( G, G' )，G和G'是倒格矢。 宏观介电常数由微观介电常数在原胞内积分取均值获得，可以证明最终结果等于其傅立叶变换系数的(0,0)项，也就是epsilon(0,0)。 但是一般程序在计算时算的不是epsilon(G,G')，而是epsilon^-1(G,G')。这样必须先取整个epsilon^-1(G,G')矩阵的逆，再取其(0,0)元，这样做时就包含了局域场效应。 如果为了计算简便，仅对epsilon^-1(G,G')的对角元取逆，就忽略了局域场效应。 一般来说，局域场效应对inhomogeneous系统影响较大。 |
» 本帖已获得的红花（最新10朵）
future_wl4楼2013-11-05 14:21:04
future_wl
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 3295.4
- 散金: 1881
- 红花: 5
- 帖子: 1404
- 在线: 676.6小时
- 虫号: 950891
- 注册: 2010-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2013-11-06 09:05:42
future_wl
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 3295.4
- 散金: 1881
- 红花: 5
- 帖子: 1404
- 在线: 676.6小时
- 虫号: 950891
- 注册: 2010-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
6楼2013-11-06 09:11:12
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2013-11-06 09:35:56
future_wl
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 3295.4
- 散金: 1881
- 红花: 5
- 帖子: 1404
- 在线: 676.6小时
- 虫号: 950891
- 注册: 2010-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
8楼2013-11-06 20:31:06
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
fzx2008: 金币+3, 谢谢指导 2013-11-07 12:34:30
fzx2008: 金币+3, 谢谢指导 2013-11-07 12:34:30
|
RPA是一种计算节点函数的理论，局域场效应是计算介电函数过程中所涉及的一种偏数学的近似手段。 这里VASP手册里写得我也看得不是太明白，按我理解应该是这样的： 算介电函数时先算的是极化函数X，然后由X计算eps^-1。由前面DFT计算得到的波函数和本征值可以根据RPA算出X0，更精确的应该是由X0算出的X。考不考虑RPA的区别主要在X上。 你的理解应该是对的。 推荐去看下DP和exciting这两个TDDFT程序的主页，上面有讲解TDDFT理论的PPT，里面对极化率节点函数讲得很清楚。 |
9楼2013-11-07 09:45:41
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
抱歉，上面打错了。应该是DP和yambo。 10楼2013-11-07 09:46:59
