| 查看: 1904 | 回复: 34 | |||
[交流]
六角晶系结构参数的优化问题
|
|||
| 由于初学VASP，不知道如何优化六角晶系结构参数，看侯老师讲义上的内容，还是不明白。真诚地希望各位高手们帮忙指教，谢谢！！ |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于根据电子衍射图样求六角晶系轴长比c/a的问题
已经有3人回复
- 结构优化参数设置
已经有6人回复
- 结构优化中问题
已经有6人回复
- VASP优化结构遇到的问题。
已经有9人回复
- 结构优化中遇到的问题
已经有9人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】六角晶系中[11-20]与(11-20)为什么不垂直？
已经有11人回复
- 【求助】一个结构优化的问题
已经有11人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【求助成功】计算能带结构时出现错误，请教
已经有19人回复
- 【求助】利用侯博的方法优化结构出现的问题（求验证）
已经有11人回复
- 【请教】请教一个六方晶系的问题
已经有12人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
» 抢金币啦！回帖就可以得到:
-
3-丁烯-1-醇（CAS：627-27-0）
+1/87
-
3-羟基丙腈（CAS：109-78-4）
+1/87
-
1-甲基-1,2-双(三甲基硅氧基)乙烷（CAS：17887-27-3）
+1/86
-
十二烷基溴化吡啶（CAS：104-73-4）
+1/84
-
香港城市大学郑星课题组诚聘博士后 (生态学、环境科学、植物学、大气科学背景)
+1/78
-
实验室富余
+1/73
-
无机非金属材料科学相关 实习研发人员招聘
+1/36
-
请问谁有2025最新影响因子Excel表？金币求助！
+1/33
-
哈工大（深圳）何自开团队诚招2027级化学工程与技术专业博士生
+1/28
-
浙江大学医学院附属第一医院招聘科研助理
+1/28
-
哈工大（深圳）何自开团队诚招化学工程与技术专业博士生
+1/27
-
中山大学计算数学方向招2027届考核制博士生
+1/18
-
高耐腐蚀性的镁合金微弧氧化工艺
+1/15
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/12
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/9
-
中山大学李国强课题组招收博士生（柔性电子方向）
+1/8
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
夏尽秋来，蓄力沉淀，静待硕果成文
+1/4
-
高弯曲强度玻璃制品的特殊加工工艺
+1/3
-
请大家抽空欣赏一下我们网站的配图
+1/3
9楼2012-08-13 07:10:47
22楼2012-08-13 15:23:31
23楼2012-08-14 09:56:41
25楼2012-08-14 20:04:31
27楼2012-08-15 19:42:43
★ ★ ★ ★ ★
奋力拼搏(金币+1): 谢谢参与
WDD880227: 金币+1, 多谢指教 2012-08-17 08:43:59
奋力拼搏: 金币+3, 非常感谢您的回复！！谢谢！1 2012-08-19 17:39:42
奋力拼搏(金币+1): 谢谢参与
WDD880227: 金币+1, 多谢指教 2012-08-17 08:43:59
奋力拼搏: 金币+3, 非常感谢您的回复！！谢谢！1 2012-08-19 17:39:42
|
看你对精度有啥要求，如果要求不是太高的话，直接设置ISIF=3，弛豫就可，最后的结构就是优化后的结果 如果要求精度较高，就得做V-E曲线，vasp manual上的： Accurate bulk relaxations with internal parameters (one) The general message is: whenever possible avoid volume relaxation with the default energy cutoff. Either increase the basis set by setting ENCUT manually in the INCAR file, or use method two suggested below, which avoids doing volume relaxations at all. If volume relaxations are the only possible and feasible option please use the following step by step procedure (which minimizes errors to a minimum): Relax from starting structure (ISMEAR should be 0 or 1, see section 6.38). Start a second relaxation from previous CONTCAR file (re-relaxation). As a final step perform one more energy calculation with the tetrahedron method switched on (i.e. ISMEAR=-5), to get very accurate and unambiguous energies (no relaxation for the final run). The final calculation should be done with PREC=High, to get very accurate energies. A few things should be remarked here: Never take the energy obtained at the end of a relaxation run, if you allow for cell shape relaxations (the final basis set might not be isotropic). Instead perform one additional static run at the end. The relaxation will give a structure which is correct to first order, the final error in the energy of step 3 is of second order (with respect to the structural errors). If you take the energy directly from the relaxation run, errors are usually significantly larger. Another important point is that the most accurate results for the relaxation will be obtained if the starting cell parameters are very close to the final cell parameters. If different runs yield different results, then the run which started from the configuration which was closest to the relaxed structure, is the most reliable one. We strongly recommend to do any volume (and to lesser extend cell shape) relaxation with an increased basis set. ENCUT=1.3 default cutoff is reasonable accurate in most cases. PREC=High does also increase the energy cutoff by a factor 1.25. At an increased cutoff the Pulay stress correction are usually not required. Accurate bulk relaxations with internal parameters (two) It is possible to avoid volume relaxation in many cases: The method we have used quite often in the past, is to relax the structure (cell shape and internal parameters) for a set of fixed volumes (ISIF=4). The final equilibrium volume and the groundstate energy can be obtained by a fit to an equation of state. The reason why this method is better than volume relaxation is that the Pulay stress is almost isotropic, and thus adds only a constant value to the diagonal elements of the stress tensor. Therefore, the relaxation for a fixed volume will give an almost correct structure. The outline for such a calculation is almost the same as in the previous section. But in this case, one has to do the calculations for a set of fixed volumes. At first sight this seems to be much more expensive than method number one (outlined in the previous section). But in many cases the additional costs are only small, because the internal parameters do not change very much from volume to volume. Select one volume and relax from starting structure keeping the volume fixed (ISIF=4 see sec. 6.24; ISMEAR=0 or 1, see section 6.38). Start a second relaxation from previous CONTCAR file (if the initial cell shape was reasonable this step can be skipped, if the cell shape is kept fixed, you never have run VASP twice). As a final step perform one more energy calculation with the tetrahedron method switched on (ISMEAR=-5), to get very accurate unambiguous energies (no relaxation for the final run). The method has also other advantages, for instance the bulk modulus is readily available. We have found in the past that this method can be used safely with the default cutoff. (see also section 9.2). |
35楼2012-08-16 23:16:05
18楼2012-08-13 09:22:44
28楼2012-08-15 20:18:30
31楼2012-08-15 20:19:11
简单回复
2012-08-12 23:35
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 非水区，拒绝灌水！ 2012-08-13 00:07:11
xueht987: 金币+1, 操作失误，见谅 2012-08-13 14:24:05
xueht987: 金币-1, 非水区，拒绝灌水！ 2012-08-13 00:07:11
xueht987: 金币+1, 操作失误，见谅 2012-08-13 14:24:05
2012-08-13 00:01
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 00:07:31
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 00:07:31
2012-08-13 00:01
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 00:07:42
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 00:07:42
2012-08-13 00:01
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 00:07:51
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 00:07:51
2012-08-13 00:02
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 00:08:08
xueht987: 金币+4, 请为本版的发展做出您应有的贡献，大家一起来营造积极、和谐的学术交流氛围！处罚过于严厉了些，故再返4金。 2012-08-13 14:30:37
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 00:08:08
xueht987: 金币+4, 请为本版的发展做出您应有的贡献，大家一起来营造积极、和谐的学术交流氛围！处罚过于严厉了些，故再返4金。 2012-08-13 14:30:37
2012-08-13 00:02
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-5, 拒绝灌水！ 2012-08-13 00:08:21
xueht987: 金币-5, 拒绝灌水！ 2012-08-13 00:08:21
坦然随风8楼
2012-08-13 07:10
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 14:15:22
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 14:15:22
坦然随风10楼
2012-08-13 07:10
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 14:15:33
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 14:15:33
坦然随风11楼
2012-08-13 07:11
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 14:15:47
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 14:15:47
坦然随风12楼
2012-08-13 07:11
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 14:15:58
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 14:15:58
liu2004m13楼
2012-08-13 08:53
回复
奋力拼搏(金币+1): 谢谢参与
xiaoyi-198814楼
2012-08-13 09:21
回复
奋力拼搏(金币+1): 谢谢参与
xiaoyi-198815楼
2012-08-13 09:21
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 14:16:35
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 14:16:35
xiaoyi-198816楼
2012-08-13 09:22
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 14:16:42
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 14:16:42
xiaoyi-198817楼
2012-08-13 09:22
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 14:17:01
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 14:17:01
xiaoyi-198819楼
2012-08-13 09:22
回复
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 14:17:12
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 14:17:12
迷路的游侠20楼
2012-08-13 10:24
回复
奋力拼搏(金币+1): 谢谢参与
zhengkning21楼
2012-08-13 15:03
回复
奋力拼搏(金币+1): 谢谢参与
奋力拼搏: 金币+1 2012-08-15 19:39:12
奋力拼搏: 金币+1 2012-08-15 19:39:12
lihong198924楼
2012-08-14 10:32
回复
奋力拼搏(金币+1): 谢谢参与
奋力拼搏: 金币+1 2012-08-15 19:43:33
奋力拼搏: 金币+1 2012-08-15 19:43:33
郑克宁26楼
2012-08-14 20:56
回复
奋力拼搏(金币+1): 谢谢参与
奋力拼搏: 金币+1 2012-08-15 19:43:59
奋力拼搏: 金币+1 2012-08-15 19:43:59
gaomeiyan29楼
2012-08-15 20:18
回复
奋力拼搏(金币+1): 谢谢参与
gaomeiyan30楼
2012-08-15 20:19
回复
奋力拼搏(金币+1): 谢谢参与
gaomeiyan32楼
2012-08-15 20:19
回复
奋力拼搏(金币+1): 谢谢参与
strive12312333楼
2012-08-16 15:03
回复
奋力拼搏(金币+1): 谢谢参与
hoborchem34楼
2012-08-16 17:07
回复
奋力拼搏(金币+1): 谢谢参与
