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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于密度泛函计算的疑惑（初学）
51/1返回列表
查看: 4618  |  回复: 32
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

物理毁一生

新虫 (初入文坛)

[交流] 关于密度泛函计算的疑惑（初学） 已有5人参与

一、原来求解薛定谔方程的波函数（本证函数）和能量（本征值）有3N个变量，KS证明通过密度（仅3个变量）也能求出总能量（不过是基态的），所以现在要对含有电子密度的能量泛函理进行求解。
二、给出K-S方程，能量可有密度表示，E[ro] = Ts[ro] + Exc[ro] + EH[ro] + Eext[ro]，其中交换关联项可采用LDA，GGA等多种近似，而其他几项都有清晰的表达式。
三、利用泛函分析，利用SCF求解使能量最低的密度函数。（其中数学方法暂时不讨论，因为我看不懂）。
四、利用得到的基态能量和电子密度计算体系的光学、电学等性质。
我的理解是这样的，那么，问题来了，挖掘机技术哪家强......
Q1：在计算过程中，我只能求得基态能量和电子密度，那么体系的波函数我能够求得吗？如果不知道波函数，怎么可能了解其它的物理性质呢？只有波函数才能反映体系的状态啊!
Q2:K-S方程中动能项Ts[ro]的表达式需要各电子的波函数，而不仅仅是电子密度，求解过程中能够得到吗？
Q3：LDA的物理图像是什么，能用LDA近似需要满足什么样的条件？
Q4：在计算过程中（使用软件时），自洽求解时方法是不是唯一的？是不是我们只用选择交换关联参数和边界限制？
（任意一个问题都可，求大神指导）
回复此楼

» 收录本帖的淘帖专辑推荐

光催化

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
尼玛，学物理，毁一生
1楼 2014-10-11 10:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

物理毁一生

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 卡开发发 at 2014-10-11 11:11:55
Q1：KS的框架虽然是“密度”泛函，但是求解过程还是借助了“（赝）波函数”来进行计算。但是这个波函数在Exc非精确的情况下是非物理的，但目前没有一个Exc是真正意义上的精确。
Q2：如上面所说，自洽过程借助的是波 ...

请问一下如果求解一个体系，我不知道它里面真实的势，在用KS方程求解时，势如何选择？貌似在DFT中，只进行了LDA（或其他的），但是对势场也没有一个明确的说明。
赞一下 回复此楼
尼玛，学物理，毁一生
5楼2014-10-12 11:19:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 33 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
Q1：KS的框架虽然是“密度”泛函，但是求解过程还是借助了“（赝）波函数”来进行计算。但是这个波函数在Exc非精确的情况下是非物理的，但目前没有一个Exc是真正意义上的精确。
Q2：如上面所说，自洽过程借助的是波函数，计算密度则通过<ψ|ψ>进行，都是基于波函数进行的。
Q3：LDA的导出基于均匀电子气，一般比较适合电子密度变化较为缓慢的体系。
Q4：这个没明白你指的方法是？自洽求解可以采用不同的展开基，比如LCAO或者是平面波等，求解的方式从数学的角度稍有差异，广义来说没什么区别。
赞一下(2人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
2楼2014-10-11 11:11:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
LZ对密度泛函和K-S方程的关系理解错了。
密度泛函是H-K的定理，K-S方程是“在密度泛函理论框架下求解问题”的方法。

如楼上已经说了，KS方法借用了一个非相互作用系统的波函数，用来表达动能。这个波函数和真实体系的波函数没啥必然联系。但是，这个波函数仍然可以用来计算各种物理性质。

Q4：自洽是个数学问题，有多种可行的数学方法，一般不需要操心。程序一般会给选项来选方法。
赞一下 回复此楼
3楼2014-10-11 23:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

物理毁一生

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by KalaShayminS at 2014-10-11 23:08:14
LZ对密度泛函和K-S方程的关系理解错了。
密度泛函是H-K的定理，K-S方程是“在密度泛函理论框架下求解问题”的方法。

如楼上已经说了，KS方法借用了一个非相互作用系统的波函数，用来表达动能。这个波函数和真实 ...

请问在自洽过程中，最后得到的是非相互作用的波函数，还是真实系统的波函数？
赞一下 回复此楼
尼玛，学物理，毁一生
4楼2014-10-12 11:13:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 33 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[教师之家] 跳槽后在研项目怎么办？ +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
信息提示
关闭
请填处理意见
关闭