版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (6934)
>虫友互识 (1836)
>休闲灌水 (815)
>基金申请 (662)
>文献求助 (355)
>硕博家园 (128)
>招聘信息布告栏 (126)
>博后之家 (124)
>考博 (114)
>有机交流 (67)
>论文投稿 (64)
>复合材料 (62)
>论文道贺祈福 (61)
>导师招生 (56)
>育儿交流 (54)
>催化 (49)

返回列表

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 iamikaruk 的 40 个金币

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

[交流] 请教密度泛函理论中屏蔽库伦势的问题

最近我在做高能电子与固体弹性散射方面的一些研究，基本思路就是通过DFT求解得到固体的屏蔽库伦势 $V\left(\vec{r}\right)$ ：
$V\left(\vec{r}\right)=V_{\text{Coulomb}}\left(\vec{r}\right)+\int \frac{e^2 n\left(\vec{r}'\right)}{\left|\vec{r}-\vec{r}'\right|} d^3\vec{r}'$
原子核的库伦势可以解析的写出来，通过FP-LAPW方法可以得到总的电子密度 $n\left(\vec{r}\right)$ 。
那么我的问题是
1. 这样得到的屏蔽库伦势应该是绝对的吧？即可否认为这里的屏蔽库伦势的绝对值是相对于无穷远处真空为0得到的？
2. 如果1是正确的，那么计算得到电子的屏蔽库伦势>0该怎么理解？固体中的电子会自动发射到真空中？

回复此楼

» 猜你喜欢

申请2024或2025年博士研究生已经有12人回复
在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠已经有13人回复
双非本科毕业论文，气人已经有12人回复
化学B02口青基代表作都是什么水平的？向大佬求助已经有10人回复
每到中夜，情难自抑已经有57人回复
化学口B0109(高分子合成)，拿青年基金一般需要怎样的文章水平？已经有22人回复
每天学术时间不能保证，能保证的只有：已经有10人回复
研0二导师分到新来的博士后靠谱吗已经有7人回复
研究生在毕业答辩时挂了，遗憾已经有13人回复
24年博士招生已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

密度泛函理论PPT 已经有16人回复
密度泛函理论新进展及应用已经有397人回复
鸟瞰密度泛函理论中文版(A_Bird's-Eye_View_of_Density-Functional_Theory) 已经有952人回复
一篇关于密度泛函理论的Review 已经有217人回复
有没有基于spin-polarized 理论的密度泛函方法啊？已经有5人回复
留学去不去（密度泛函理论的问题）已经有8人回复
【求助】求论文《應用第一原理密度泛函理論研究鍍在過渡金屬上的單層石墨烯的能帶結構已经有4人回复
【求助】含时密度泛函计算金团簇的激发态输出出错已经有8人回复
【讨论】密度泛函理论与薛定谔方程的关系？已经有7人回复
【求助】请教会用“密度泛函理论”的师兄师姐已经有11人回复
【求助】第一性原理与密度泛函理论的关系是什么已经有15人回复
求密度泛函理论的离散变分法在化学和材料物理学中的应用 pdf 已经有7人回复
【求助】审稿人关于密度泛函的问题，请帮忙回答一下！已经有4人回复
【讨论】密度泛函理论中“密度”究竟指什么已经有16人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2012-12-12 13:29:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+3, 谢谢指教！ 2012-12-12 19:20:39
iamikaruk: 金币+5, 多谢指教 2012-12-12 19:33:33

我的理解是这样的：
由于势函数加上或者减去一个常数不影响相对能级的大小，所以这里的屏蔽库伦势的数值不是一个绝对值
由于在凝聚态中，电子的能级是小于真空能级的，所以如果想按照真空能级对得到的屏蔽库伦势进行标度的话，可以沿着远离原子核或者固体表面的方向计算静电势，离的足够远的地方势能趋于定值，这个值就是真空能级的能量，然后以这个能量作为参考态就可以

赞一下

回复此楼

3楼2012-12-12 19:19:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

引用回帖:

3楼: Originally posted by uuv2010 at 2012-12-12 19:19:23
我的理解是这样的：
由于势函数加上或者减去一个常数不影响相对能级的大小，所以这里的屏蔽库伦势的数值不是一个绝对值
由于在凝聚态中，电子的能级是小于真空能级的，所以如果想按照真空能级对得到的屏蔽库伦势进 ...

你的意思是固体本身相对于真空有个功函数？然后在DFT的框架下这个得要通过slab方式的计算得到功函数？
我不理解的一点是如果用的势函数是 $V_{\text{Coulomb}}\left(\vec{r}\right)=-Z\frac{e^2}{r}$ 的形式，以无穷远处真空为0，这样处理结果得到的也不是绝对的静电势？
不过实际处理中我倒是可以用 $V_{\text{Coulomb}}\left(\vec{r}\right)=-Z\frac{e^2}{r}+C$ 去拟合很靠近原子核处的DFT结果，然后把这个常数C去掉。

赞一下

回复此楼

4楼2012-12-12 19:46:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
帖子: 390
在线: 236.7小时
虫号: 582657

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

屏蔽库伦势 GW方法理的W不就是吗？不可以直接得出？

赞一下

回复此楼

5楼2012-12-12 23:34:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

引用回帖:

5楼: Originally posted by ApolloYang at 2012-12-12 23:34:23
屏蔽库伦势 GW方法理的W不就是吗？不可以直接得出？

我关心这个得到的是绝对值还是相对值，是否要把得到的库伦势做平移
GW方法的W是绝对的？我没有用过GW方法，请教一下。

赞一下

回复此楼

6楼2012-12-13 08:44:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

4楼: Originally posted by iamikaruk at 2012-12-12 19:46:42
你的意思是固体本身相对于真空有个功函数？然后在DFT的框架下这个得要通过slab方式的计算得到功函数？
我不理解的一点是如果用的势函数是V_{\text{Coulomb}}\left(\vec{r}\right)=-Z\frac{e^2}{r}的形式，以 ...

是的，由于这个常数的存在，所以加上一个常数拟合的方法可能比较合适

赞一下

回复此楼

7楼2012-12-13 10:22:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
帖子: 390
在线: 236.7小时
虫号: 582657

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

5楼: Originally posted by ApolloYang at 2012-12-12 23:34:23
屏蔽库伦势 GW方法理的W不就是吗？不可以直接得出？

给个邮件
我发几个GW的介绍论文给你看看

GW方法里屏蔽库伦势我没理解透彻。

VASP里的“自能＝iGW” G是格林函数，
要得出W，估计还要修改main.F

赞一下

回复此楼

8楼2012-12-13 10:22:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

引用回帖:

7楼: Originally posted by uuv2010 at 2012-12-13 10:22:52
是的，由于这个常数的存在，所以加上一个常数拟合的方法可能比较合适...

试了一下，实际上不可操作
屏蔽库伦势在很接近于原子核的地方才逼近于库伦势，这个范围<0.05埃，实际上很难达到那么高的采样精度，就算达到数据点也不够用于拟合。
昨晚上推导了一个小时的公式，完全从电荷密度出发得到屏蔽库伦势。早上高兴的过来写了个代码验证，发现一个单胞内的电荷密度在0.01埃的sampling下总电荷密度还不收敛

赞一下

回复此楼

9楼2012-12-13 11:16:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

引用回帖:

8楼: Originally posted by ApolloYang at 2012-12-13 10:22:56
给个邮件
我发几个GW的介绍论文给你看看

GW方法里屏蔽库伦势我没理解透彻。

VASP里的“自能＝iGW” G是格林函数，
要得出W，估计还要修改main.F...

我的邮件地址是thuxielin@gmail.com
vasp不是全电子方法，没有考虑core electrons，它得到的势真实么？

赞一下

回复此楼

10楼2012-12-13 11:16:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
帖子: 390
在线: 236.7小时
虫号: 582657

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
iamikaruk: 金币+5, 多谢指教 2012-12-13 12:01:07

引用回帖:

9楼: Originally posted by iamikaruk at 2012-12-13 11:16:01
试了一下，实际上不可操作
屏蔽库伦势在很接近于原子核的地方才逼近于库伦势，这个范围<0.05埃，实际上很难达到那么高的采样精度，就算达到数据点也不够用于拟合。
昨晚上推导了一个小时的公式，完全从电荷密 ...

你这里的屏蔽库伦势和GW里的屏蔽库伦势是不一样的概念。

GW的屏蔽库伦势是相对裸露库伦势而言的。
GW讨论的是激发态。假设电子被激发到VB，在CB上留下一个空穴，这个CB上的空穴会吸引其它的电子形成一个准粒子。不恰当地说，就向空穴外围分布了许多的电子。如此而言，空穴和CB里电子的作用就由裸露的库伦作用弱化为屏蔽库仑作用。

赞一下

回复此楼

12楼2012-12-13 11:54:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

引用回帖:

12楼: Originally posted by ApolloYang at 2012-12-13 11:54:18
你这里的屏蔽库伦势和GW里的屏蔽库伦势是不一样的概念。

GW的屏蔽库伦势是相对裸露库伦势而言的。
GW讨论的是激发态。假设电子被激发到VB，在CB上留下一个空穴，这个CB上的空穴会吸引其它的电子形成一个准粒 ...

这是进一步的东西了，要考虑高能电子激发固体内部电子。
我现在仅仅研究高能电子与固体的弹性散射作用，所以普通的DFT就可以胜任了。

赞一下

回复此楼

13楼2012-12-13 12:00:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

1ST强帖: 1
应助: 11 (小学生)
金币: 4950.3
帖子: 1362
在线: 543.8小时
虫号: 1315985

引用回帖:

debug了一下，发现fortran在读入电子密度的xsf文件时出错了，重新改代码

赞一下

回复此楼

14楼2012-12-13 12:02:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

简单回复

zfl891052楼

2012-12-12 15:15 回复

likai113911楼

2012-12-13 11:34 回复

uuv201015楼

2012-12-13 12:10 回复

引用回帖:

9楼: Originally posted by iamikaruk at 2012-12-13 11:16:01 试了一下，实际上不可操作屏蔽库伦势在很接近于原子核的地方才逼近于库伦势，这个范围<0.05埃，实际上很难达到那么高的采样精度，就算达到数据点也不够用于拟合。昨晚上推导了一个小时的公式，完全从电荷密 ...

相关版块跳转我要订阅楼主 iamikaruk 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[硕博家园] 每天学术时间不能保证，能保证的只有： +10	hahamyid 2024-05-27	10/500	2024-06-01 21:11 by 小小芝麻影
[教师之家] 研究生在毕业答辩时挂了，遗憾 +10	spacr2022 2024-05-27	13/650	2024-06-01 17:45 by keyaner23
[硕博家园] 论大家对6070后普通教授导师的看法 +6	SNaiL1995 2024-05-28	10/500	2024-06-01 16:55 by 篮球易飞
[基金申请] 九部门发文：不得将专利授权数量作为人才评价、项目评审、职称评定、高校评价等的条件 +16	sjtu2012 2024-05-28	19/950	2024-06-01 15:24 by 天外飞去来
[基金申请] B口人才项目 +9	WOWO159357 2024-05-29	19/950	2024-06-01 14:24 by linxuhuizj
[硕博家园] 实验室太吵闹，无法安静学习，怎么办？ +6	utahh 2024-05-31	11/550	2024-06-01 13:27 by utahh
[基金申请] 博后特别资助状态变化 +24	随梦而飞2017 2024-05-30	35/1750	2024-06-01 10:10 by 青岛阳仔
[基金申请] 数理的人才答辩通知发了吗 +7	yzy3327 2024-05-30	7/350	2024-06-01 09:15 by 撒哈拉西部
[考博] 广东以理材料系碳点与功能材料课题组 — 2博士名额 / 科研助理 +4	小城夜很美 2024-05-27	11/550	2024-05-31 21:26 by 小城夜很美
[硕博家园] 求问endnote上Wiley引用格式在哪下载 +3	Elon. 2024-05-30	3/150	2024-05-31 10:13 by crazy peng
[文学芳草园] 对对子啊 +5	天若孤独 2024-05-29	7/350	2024-05-31 09:00 by wjykycg
[基金申请] 入职高校3年发表10+SCI，尽人事听天命 +29	kaoyan250 2024-05-27	40/2000	2024-05-31 08:44 by Xiaolin81
[电化学] 有锂金属软包电池代工公司代做软包嘛 +5	yyywj995 2024-05-26	6/300	2024-05-30 08:33 by 那年的旧城
[论文投稿] 审稿专家比较坚定的让补充实验，但实在没法补充实验，修回还有希望吗？ (EPI+1) 3+3	qweasd12345 2024-05-29	6/300	2024-05-30 08:11 by qweasd12345
[硕博家园] 又想换工作 +15	brightmj 2024-05-27	20/1000	2024-05-29 23:25 by zyqchem
[基金申请] 信息学部函评结束了吗？ +6	ducan21 2024-05-28	7/350	2024-05-29 12:10 by WORLD0256
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊，不投稿收费SC，对分区没有要求 5+3	xintangren 2024-05-28	4/200	2024-05-29 10:46 by xintangren
[基金申请] E05青基有几个评审 +4	KYXY123 2024-05-28	4/200	2024-05-28 19:25 by popt2t
[有机交流] 机理求助 120+4	15147165026 2024-05-26	10/500	2024-05-28 14:42 by 江东闲人
[硕博家园] 周日 +6	1加油哦棒 2024-05-26	9/450	2024-05-27 10:30 by hahamyid

24小时热门版块排行榜

[交流] 请教密度泛函理论中屏蔽库伦势的问题

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:查看全部散金贴

» 抢金币啦！回帖就可以得到:
查看全部散金贴