| 查看: 1557 | 回复: 4 | ||||
[交流]
请大牛指教。关于截断动能、原子数目，K点密度
|
|
基于周期性边界条件的VASP，要用超胞堆叠整个3维空间，这样的话，原子数目是无限的，而实际计算中，为什么原子数目多了就算不动？ 长话短说，不提公式，只说原理： 1、类比于原子轨道。晶态材料中有‘晶体轨道’。也就是单电子函数。 2、原子中的电子，按泡利不相容原理，从低能级的轨道开始填充。 VS 晶体中的电子，也按泡利不相容原理从低能量的‘晶体轨道’开始填起。 3、怎样得到体系的能量呢？（我这里说的不是密度泛函版，这个太难描述了。只说一个最简单的原型） a、求出体系的所有‘晶体轨道’，当然个轨道所对能量也顺便得出来了。 b、电子一个一个，从低能轨道排到高能轨道。各电子的能量，就是它所在轨道的能量。 c、各电子能量求和就OK了。 4、可见，关键在于求解各‘晶体轨道’。基于Bloch地理，计算机是这样处理的。 a、用超胞中各原子的原子轨道的线性叠加来表示‘晶体轨道’。 b、利用变分原理。把该问题转化为求解一个齐次线性方程组。也就是对角化一个矩阵，就可求出各系数。 c、可见。如果超胞中的原子越多，则参与计算的原子轨道及越多。需要对角化的矩阵就越大。所有不能处理太大的超胞。 5、另外，每个原子的原子轨道是无限多条。晶体轨道也是无限多条。如果不做处理，那个矩阵是无穷维的。‘无穷×无穷’的矩阵 a、引入截断动能。能量高于这个值的原子轨道则不参与计算。 这样一来‘无穷×无穷’简化为 ‘无穷×有穷’ b、还有一个方向是无限维。该K空间踩点来处理了。你不是无穷吗？我只踩有限个点来计算。 于是： ‘无穷×有穷’再次简化为 ‘有穷×有穷’ 总结：胞越大，截断动能越大，K点越大。要对角化的矩阵就越大。 水平有限，必有疏漏。恭请指正。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 我的差分电荷密度正常吗？能不能给个提示。谢谢啦！
已经有19人回复
- TEM 原子序数大的原子对比度更大吗
已经有8人回复
- 求大牛指教：怎么根据DTA-TG曲线制定烧结制度~
已经有4人回复
- 两个原子的局域态密度在费米处重叠
已经有4人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 关于K点选择与原子间的力
已经有9人回复
- 【求助完毕】这篇JACS的communicaton是如何从电荷密度图看出共价键来的？
已经有12人回复
- 【求助】请教几个有关表面活性剂的问题(标题要长，金币要多，虚心求指教，额外金币奉�
已经有13人回复
- 【请教】向大侠请教，已知空间群，并且已经有粉末衍射数据，怎么确定原子位置？
已经有4人回复
- 【求助】谁做ZSM-5分子筛，请高手指教
已经有24人回复
- 【求助】电子密度图中背景电子密度的概念？
已经有12人回复
2楼2012-07-06 12:05:27
3楼2012-07-06 12:33:33
4楼2012-07-06 15:50:41
5楼2012-07-06 16:27:30