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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于密度泛函计算的疑惑（初学）
汕头大学海洋科学接受调剂
53/1返回列表
查看: 4198  |  回复: 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的回帖

KalaShayminS

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
25楼: Originally posted by tsglss at 2014-10-14 23:09:12
首先波函数怎么求到，而且它应该还是个数值，更糟的是算符本身就不是个数，就算模拟出来了，也不是很准确，因为T=0.感觉无重要意义
...

求解KS方程自然就得到KS波函数了，求解过程是个复杂的数学问题，DFT程序就是做这个的。上面已经说过了，波函数用基组来表示。算符也可以用基组来表示。因为KS波函数的形式是Slater行列式，所以其产生湮灭算符是显然的。
从这个和下面那个回复来看，里面提到的对波函数、算符、基态、混沌、温度、费米子、相互作用的概念的理解全都是不正确的。
赞一下 回复此楼
30楼2014-10-15 21:59:13
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +5 关忆北. 2026-04-14 5/250 2026-04-15 09:30 by 逆水乘风
[考研] 求调剂学校 +12 不会吃肉 2026-04-13 14/700 2026-04-15 09:15 by 逆水乘风
[考研] 药学求调剂 +11 喽哈加油 2026-04-14 13/650 2026-04-14 21:14 by qingfeng258
[考研] 290求调剂 +21 luoziheng 2026-04-10 23/1150 2026-04-14 15:49 by zs92450
[考研] 一志愿沪9，326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 机械工程313分找工科调剂 +4 双一流本科机械 2026-04-08 4/200 2026-04-14 07:32 by Abskk
[考研] 302求调剂 +10 易！? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 346分，工科0854求调剂，专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 0860004 求调剂 309分 +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 332，085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
信息提示
关闭
请填处理意见
关闭