| 查看: 5140 | 回复: 9 | |||
1243946679至尊木虫 (正式写手)
逍遥游
|
[求助]
VASP各种能量的选择与区别
|
|
请问：1.在计算结合能，测试ENCUT、KPOINTS、真空层厚度的收敛时，都是选择energy without entropy 进行比较吗？（看到vasp论坛上的帖子，迷茫了。链接如下http://cms.mpi.univie.ac.at/vasp ... en+subtracted#p5621） 2.vasp的MD计算时，OSZICAR中能量举例如下：T=296 E= -1424.9819 F= -1542.1294 E0= -1542.1346 EK= 7.6232 SP= 110 SK= 1.4E-4。手册上说，E is the total free energy (including the kinetic energy of the ions and the energy of the Nos ′ e thermostat)， F is the total free energy (at this point the energy of the reference atom has been subtracted)。我的理解是，F是赫姆霍茨自由能，F=U-T*S，E是所计算体系的内能加上热浴的内能，即，E=U（体系）+U（热浴）=F+T*S+SP+SK，但等式两侧并不相等，当等式右侧加上EK后，才基本相等（并不完全相等）。如何理解E和F这两个概念，重点是，如何得到所计算体系的势能，是 F-EK（或者，F+T*S-EK）吗？ |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
0856材料专硕353求调剂
已经有5人回复
-
求调剂
已经有6人回复
-
293求调剂
已经有4人回复
-
306求0703调剂一志愿华中师范
已经有10人回复
-
【双一流院校新能源、环境材料，材料加工与模拟招收大量调剂】
已经有7人回复
-
08工学调剂
已经有12人回复
-
材料与化工328分调剂
已经有4人回复
-
279分求调剂 一志愿211
已经有20人回复
-
一志愿华东理工大学081700，初试分数271
已经有6人回复
-
085404电子信息284分求调剂
已经有3人回复
wuli8
荣誉版主 (知名作家)
…………
- 1ST强帖: 2
- 应助: 35 (小学生)
- 贵宾: 12.924
- 金币: 20190.4
- 散金: 15888
- 红花: 88
- 沙发: 4
- 帖子: 7840
- 在线: 1114.6小时
- 虫号: 465889
- 注册: 2007-11-23
- 专业: 物理学I
- 管辖: 计算模拟
|
让你彻底明白所有vasp里能量的物理意义 （一）首先我们应明白，固体的结合能就是固体的内能E（结合）=U（内能）， 原因如下： 一般情况都把孤立原子的能量作为能量参考点。前段时间有个同学问VASP中得出的绝对能量是相对于什么的，其实就是相对孤立原子得。 （二）其次我们根据自由能与内能之间的关系F=U-TS 而且我们都知道VASP的所有计算都是在绝对0度下的情况，T=0代入上式，有F=U。所以结合就等于内能等于自由能。肯定有Free energy TOTEN=energy without entropy恒成立... 这时候肯定有人会说不对啊，可以看VASP手册，候博的参考书作证，肯定不对得。 现在我告诉你确实它们二者确实有区别，区别在下面的情况 （1）当我们用ISMEAR=-5时，费米能这儿没有展宽，它算出来的就是完全在绝对0度的能量。Free energy TOTEN=energy without entropy恒成立。 （2）有时为了在数学上处理的方便，为了更容易积分，我们也用ISMEAR！=-5（！=是不等于的意思）的方法，这个时候费米能这儿有一定的展宽。此时，我们容易想到，有展宽不就是相当有一定的熵值吗？所以这个时候虽然算的是绝对0度的情况，但是有一定的熵值（我们应明白，这个熵值不是由一定的温度带来的，而是数学处理的结果）。所以在SMEAR！=-5的方法我们会发现Free energy TOTEN和energy without entropy有一定的差别。此时energy without entropy是Free energy TOTEN在SIGMA趋于0的极限。 注意:（1）有人在算单个原子的能量时会发现单个原子的能量虽然很小但并不是0,但是按我上面的推导,固体中的结合能是相对孤立体系的能量而来的,所以单个原子得到的TOTEN肯定是0啊,原因在于我们的POTCAR不可能绝对合理,而且我们也知道计算单个原子的能量就是为了检测赝势,单原子得到的TOTEN越小说明赝势越好。但一般不会正好是0. （2）如果你注意的话，energy without entropy与Free energy TOTEN在SIGMA趋于0也不是完全相等，但是也会发现它们之间的差别在10E-3左右，原因在于计算机求积分、求极限不能像我们人一样达到任意的精度。 还有，指出你几个概念上的问题，供探讨 第一，关于结合能。结合能是定义为相距无穷远的原子结合形成一定结构的物质所放出的能量 第二，关于单点能。它是第一性原理计算直接得到的能量，或者说是赝能，是一个空间点阵平均每阵点上采用赝势计算所得到的能量，其中包含了结合能的贡献，但是更多的，也包含了靠近芯区附近的电子在采用赝势近似下的能量，这一部分能量既不是原子芯区附近电子能量的真实反应，也不会影响化学键性质，不会对结合能有所贡献。 第三，Free energy TOTEN是体系总能，要减去阵点上分布的原子的能量再除以平均原子数才是结合能（当然，这个和你的计算脚本的设计有关），而且这还没有考虑不加展宽时没有被计算到的能带的因素。 第四，是否考虑能级展宽，和结合能的定义没有关系。 第五，结合能计算时对单个原子能量的计算应该只计算Gamma点能量，且用消除简并。 关于你说的这几点做个讨论 （一）关于结合能。你说“结合能是定义为相距无穷远的原子结合形成一定结构的物质所放出的能量” 你和我说的没区别，我说的是结合能是相对于“孤立原子做参考点的”，也就是它与周围任何原子没有相互作用，和你所说的相距无穷远一回事，我这个好像没有任何错误。 （二）关于单点能。你说“它是第一性原理计算直接得到的能量，或者说是赝能，是一个空间点阵平均每阵点上采用赝势计算所得到的能量，其中包含了结合能的贡献，但是更多的，也包含了靠近芯区附近的电子在采用赝势近似下的能量，这一部分能量既不是原子芯区附近电子能量的真实反应，也不会影响化学键性质，不会对结合能有所贡献”。 我赞同你的大部分观点，也提出你说的几点错误，单点能准确的来说它包含了所有哈密顿的量，而且这儿的单点能不是你所说的“平均每个原子的能量”，而是你计算的整个原胞的能量。但是这个能量有一个参考点。你可以看候博得，也可以看我回的下一个贴子，至于“影响不影响成键之类的内容”固体力学上已经说的很清楚了。 （三）你说“Free energy TOTEN是体系总能，要减去阵点上分布的原子的能量再除以平均原子数才是结合能（当然，这个和你的计算脚本的设计有关），而且这还没有考虑不加展宽时没有被计算到的能带的因素” 我不赞同你后面说的几点。Free energy TOTEN从字面意思上我们也知道它的结果是自由能，你可以说它是总能，因为根据我上面的推导，它们至少在数值是相等得。不在于你把它说成什么，你就是把它说成总能，其实它还是等于结合能，等于自由能，等于内能。至于除不除原子总数在于你想得到的是平均每个原子的还是总体系的，这在于个人。考虑不考虑展宽，那要看ISMEAR等于几，做几个实例就会感觉到它考虑没考虑了。 （四）你说“是否考虑展宽和结合能的定义没有关系” 我也没说和它的定义有什么关系啊，但是由于数学处理带来的误差，它对结合能的结果有一定的影响啊。 （五）对单个原子的结合能的计算应该只计算Gamma点的能量，且用削除简并。 你说得第五点我不懂，我算单个原子的能量时一直是按三个表面的构造方法来算的，也没想过什么简并。还望有高手给我帮助，第五点怎么理解。 我给你写一遍候柱锋书上怎么说结合能得，然后和我说得对照一下：“VASP计算得到的总能己经减去了在以原子参考组态计算得到的总能（也就是构造赝势时，得到的总能，对应于POTCAR文件中的EATOM）”要得到准确的结合能，还需要减去前面单个原子计算中的第2种情况计算得到的修正值。 所以：：：：：我们可以认为vasp得到的总能就是结合能，但是由于赝势带来的影响，它存在一定的误差。为什么要减去修正值。我在这儿做个解释，由于我们在计算单原子时也用这个赝势了，计算总能时也用这个赝势了，两都有误差，但对于同一种元素来说，它们的误差相近，这样相减的话，正好能把误差消去（不可能完全消去，但是使我们的结果更加精确）。 这只是为了提高精度，就像候博说得“为了得到准确的结合能”，候博不是说从前面的计算来得到结合能，这只是为了考虑POTCAR带来的误差，而不是概念上差别。 关于你说的这几点做个讨论 （一）关于结合能。你说“结合能是定义为相距无穷远的原子结合形成一定结构的物质所放出的能量” 你和我说的没区别，我说的是结合能是相对于“孤立原子做参考点的”，也就是它与周围任何原子没有相互作用，和你所说的相距无穷远一回事，我这个好像没有任何错误。 =============== 这里你说的是没有错误，但是我觉得有必要先澄清一下。 （二）关于单点能。你说“它是第一性原理计算直接得到的能量，或者说是赝能，是一个空间点阵平均每阵点上采用赝势计算所得到的能量，其中包含了结合能的贡献，但是更多的，也包含了靠近芯区附近的电子在采用赝势近似下的能量，这一部分能量既不是原子芯区附近电子能量的真实反应，也不会影响化学键性质，不会对结合能有所贡献”。 我赞同你的大部分观点，也提出你说的几点错误，单点能准确的来说它包含了所有哈密顿的量，而且这儿的单点能不是你所说的“平均每个原子的能量”，而是你计算的整个原胞的能量。但是这个能量有一个参考点。你可以看候博得，也可以看我回的下一个贴子，至于“影响不影响成键之类的内容”固体力学上已经说的很清楚了。 ============== 从你的回复中，我可以知道你肯定没有学过晶体学或者空间群理论，你应该看看晶体学国际表中对于阵点的定义，阵点并不是每个原子，这里你的理解有问题，阵点是一个抽象点，一个晶体中包含所有对称性的可以仅通过平移来构造整个晶体的结构所占据的位置就是一个阵点，换句话说，一个阵点，就是一个满足平移对称性的原子集团，且该集团内部的位置满足该晶体结构的全部对称性，而且它不仅仅是“原胞” （三）你说“Free energy TOTEN是体系总能，要减去阵点上分布的原子的能量再除以平均原子数才是结合能（当然，这个和你的计算脚本的设计有关），而且这还没有考虑不加展宽时没有被计算到的能带的因素” 我不赞同你后面说的几点。Free energy TOTEN从字面意思上我们也知道它的结果是自由能，你可以说它是总能，因为根据我上面的推导，它们至少在数值是相等得。不在于你把它说成什么，你就是把它说成总能，其实它还是等于结合能，等于自由能，等于内能。至于除不除原子总数在于你想得到的是平均每个原子的还是总体系的，这在于个人。考虑不考虑展宽，那要看ISMEAR等于几，做几个实例就会感觉到它考虑没考虑了。 =============== 这个能量确切的说应该是叫做考虑电子振动熵的体系总自由能，当不考虑展宽的时候，它是等于总能的，如果你读过Vasp的代码，就知道TOTEN在vasp的计算中就是总能，这个和结合能不是一个概念，还包含有非成键部分的贡献，至于内能的定义，如果你阅读过塞兹的现代固体理论，或者Pauling的书，或者读过Morse当时提出morse势的那篇文献，就应该知道，固体物理中所使用的内能，指的是离子实的动能和原子的“结合能”之和——这里结合能之所以要打引号，是因为按定义，是要形成稳定结构或者亚稳态结构时才能称之为结合能，内能定义并没有考虑粒子芯区附近电子能量的影响，正如你所说，是“相对于“孤立原子做参考点的””，在0K下已经不考虑动能，因此就应是总能减去孤立原子的能量和才行，至于结合能，则是稳定状态下结构的这个能量。 （四）你说“是否考虑展宽和结合能的定义没有关系” 我也没说和它的定义有什么关系啊，但是由于数学处理带来的误差，它对结合能的结果有一定的影响啊。 （五）对单个原子的结合能的计算应该只计算Gamma点的能量，且用削除简并。 你说得第五点我不懂，我算单个原子的能量时一直是按三个表面的构造方法来算的，也没想过什么简并。还望有高手给我帮助，第五点怎么理解。 ============= 简单的操作是计算单个原子能量只考虑Gamma点，然后三边都设置在10A以上，且不相等 至于原因，应该去查量子力学的书，记得本科的时候，老师都会讲到的。 根据F=U-TS，E（结合）=E（bulk）-nE(单个离散)。而E（bulk）就是U，通常我们取E（离散）为参考点。也就是把E（离散）看为0，这样推出来，在T=0时，F=E（结合）。 它是包含你所说得那些，而且就像我在前面的贴子中说得，它就是总得H得到的能量，但是它有一个参考，而这个参考就是离散原子的能量 |
2楼2015-11-03 16:01:50
wuli8
荣誉版主 (知名作家)
…………
- 1ST强帖: 2
- 应助: 35 (小学生)
- 贵宾: 12.924
- 金币: 20190.4
- 散金: 15888
- 红花: 88
- 沙发: 4
- 帖子: 7840
- 在线: 1114.6小时
- 虫号: 465889
- 注册: 2007-11-23
- 专业: 物理学I
- 管辖: 计算模拟
|
让你彻底明白所有vasp里能量的物理意义 http://muchong.com/bbs/viewthread.php?tid=2628844&fpage=1&target=blank |
3楼2015-11-03 16:02:41
1243946679
至尊木虫 (正式写手)
逍遥游
- 应助: 0 (幼儿园)
- 金币: 10509.1
- 散金: 299
- 帖子: 611
- 在线: 454.1小时
- 虫号: 2071207
- 注册: 2012-10-18
- 专业: 金属非晶态、准晶和纳米晶
4楼2015-11-03 19:05:08
凌云志123
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1203.9
- 散金: 454
- 红花: 6
- 帖子: 452
- 在线: 98.7小时
- 虫号: 2073907
- 注册: 2012-10-20
- 性别: GG
- 专业: 金属结构材料
5楼2015-12-18 10:55:38
1243946679
至尊木虫 (正式写手)
逍遥游
- 应助: 0 (幼儿园)
- 金币: 10509.1
- 散金: 299
- 帖子: 611
- 在线: 454.1小时
- 虫号: 2071207
- 注册: 2012-10-18
- 专业: 金属非晶态、准晶和纳米晶
6楼2015-12-20 11:32:38
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
7楼2015-12-21 02:57:15
1243946679
至尊木虫 (正式写手)
逍遥游
- 应助: 0 (幼儿园)
- 金币: 10509.1
- 散金: 299
- 帖子: 611
- 在线: 454.1小时
- 虫号: 2071207
- 注册: 2012-10-18
- 专业: 金属非晶态、准晶和纳米晶
8楼2015-12-21 09:19:16
9楼2016-04-20 11:17:01
1243946679
至尊木虫 (正式写手)
逍遥游
- 应助: 0 (幼儿园)
- 金币: 10509.1
- 散金: 299
- 帖子: 611
- 在线: 454.1小时
- 虫号: 2071207
- 注册: 2012-10-18
- 专业: 金属非晶态、准晶和纳米晶
10楼2016-04-20 12:20:17
