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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请大牛指教。关于截断动能、原子数目，K点密度
51/1返回列表
查看: 1321  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dyc_2008

捐助贵宾 (著名写手)

[交流] 请大牛指教。关于截断动能、原子数目，K点密度

基于周期性边界条件的VASP，要用超胞堆叠整个3维空间，这样的话，原子数目是无限的，而实际计算中，为什么原子数目多了就算不动？

长话短说，不提公式，只说原理：
1、类比于原子轨道。晶态材料中有‘晶体轨道’。也就是单电子函数。
2、原子中的电子，按泡利不相容原理，从低能级的轨道开始填充。  VS   晶体中的电子，也按泡利不相容原理从低能量的‘晶体轨道’开始填起。
3、怎样得到体系的能量呢？（我这里说的不是密度泛函版，这个太难描述了。只说一个最简单的原型）
a、求出体系的所有‘晶体轨道’，当然个轨道所对能量也顺便得出来了。
b、电子一个一个，从低能轨道排到高能轨道。各电子的能量，就是它所在轨道的能量。
c、各电子能量求和就OK了。

4、可见，关键在于求解各‘晶体轨道’。基于Bloch地理，计算机是这样处理的。
a、用超胞中各原子的原子轨道的线性叠加来表示‘晶体轨道’。
b、利用变分原理。把该问题转化为求解一个齐次线性方程组。也就是对角化一个矩阵，就可求出各系数。
c、可见。如果超胞中的原子越多，则参与计算的原子轨道及越多。需要对角化的矩阵就越大。所有不能处理太大的超胞。
5、另外，每个原子的原子轨道是无限多条。晶体轨道也是无限多条。如果不做处理，那个矩阵是无穷维的。‘无穷×无穷’的矩阵
a、引入截断动能。能量高于这个值的原子轨道则不参与计算。
这样一来‘无穷×无穷’简化为 ‘无穷×有穷’
b、还有一个方向是无限维。该K空间踩点来处理了。你不是无穷吗？我只踩有限个点来计算。
于是： ‘无穷×有穷’再次简化为 ‘有穷×有穷’

总结：胞越大，截断动能越大，K点越大。要对角化的矩阵就越大。

水平有限，必有疏漏。恭请指正。
回复此楼

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

VASP

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-07-06 11:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
★ ★
dyc_2008(金币+1): 谢谢参与
dyc_2008: 金币+1 2012-09-11 09:10:26
只是说其中一点，也就是3c所说的“c、各电子能量求和就OK了。”
恐怕是不可以的，即简单地将各个电子的能量求和并不能描述体系的真实的能量状态，至少交换能没有考虑进去。
赞一下 回复此楼
2楼2012-07-06 12:05:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)
引用回帖:
2楼: Originally posted by uuv2010 at 2012-07-06 12:05:27
只是说其中一点，也就是3c所说的“c、各电子能量求和就OK了。”
恐怕是不可以的，即简单地将各个电子的能量求和并不能描述体系的真实的能量状态，至少交换能没有考虑进去。

我已经说了这是最基本的模型
不是DFT
没有考虑交互作用
仅仅就是单电子方程

谢谢参与
赞一下 回复此楼
3楼2012-07-06 12:33:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
★ ★ ★
dyc_2008: 金币+3 2012-09-11 09:10:35
无论DFT还是HF，都是需要考虑交换能的，呵呵
如果是想计算氢原子的能级的话(这个可能是LZ所说的最简单模型)，没有必要考虑原子个数和多k点计算
另外，关于“为什么原子数目多了就算不动”的问题，是因为DFT的计算量是随着原子数目的立方的关系增加的
赞一下 回复此楼
4楼2012-07-06 15:50:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)
引用回帖:
4楼: Originally posted by uuv2010 at 2012-07-06 15:50:41
无论DFT还是HF，都是需要考虑交换能的，呵呵
如果是想计算氢原子的能级的话(这个可能是LZ所说的最简单模型)，没有必要考虑原子个数和多k点计算
另外，关于“为什么原子数目多了就算不动”的问题，是因为DFT的计 ...

感谢您老的热情关注。（老，不是真的老，而是一种尊称）
我就是想问问，VASP的骨架是不是这样的？


交互作用什么的，就是血肉了。不可忽略。
赞一下 回复此楼
5楼2012-07-06 16:27:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dyc_2008 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 283求调剂，材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
关闭
请填处理意见
关闭