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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 六角晶系结构参数的优化问题
351/1返回列表
查看: 1701  |  回复: 34
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

奋力拼搏

至尊木虫 (文坛精英)

[交流] 六角晶系结构参数的优化问题

由于初学VASP，不知道如何优化六角晶系结构参数，看侯老师讲义上的内容，还是不明白。真诚地希望各位高手们帮忙指教，谢谢！！
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2012-08-12 22:57:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

坦然随风

木虫 (著名写手)
★ ★
奋力拼搏(金币+1): 谢谢参与
奋力拼搏: 金币+1 2012-08-15 19:44:18
这个不懂，只能来领金币了
赞一下 回复此楼
9楼2012-08-13 07:10:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lizaxuzhu

新虫 (著名写手)
★
奋力拼搏(金币+1): 谢谢参与
xueht987: 恭喜什么呢？这里不是道贺的地方亲！！！ 2012-08-13 14:17:55
恭喜~~~~~
回复此楼
18楼2012-08-13 09:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)
★ ★ ★ ★
奋力拼搏(金币+1): 谢谢参与
WDD880227: 金币+1, 多谢提示 2012-08-13 19:07:31
奋力拼搏: 金币+1 2012-08-15 19:39:17
奋力拼搏: 金币+1 2012-08-15 19:44:26
楼主这问题和“vasp怎么用啊？”性质差不多啊。。。。。你不说具体点，谁都没办法回答，只好灌水了
INCAR，事实上只要写上ISIF=3，其他几乎所有参数采用默认值都不会有很大问题
KPOINTS注意用Gammacenter
POSCAR自己写吧…………
暂时想到这么多
赞一下 回复此楼
22楼2012-08-13 15:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cclw

新虫 (小有名气)
★ ★ ★
奋力拼搏(金币+1): 谢谢参与
WDD880227: 金币+1, 多谢指教 2012-08-15 08:59:46
奋力拼搏: 金币+1, 非常感谢您的回复，谢谢您！！ 2012-08-15 19:43:09
六角就两个参数a和c，可以分别优化。建议用a作为scale,先固定c/a，优化a,直至能量最低；然后优化c/a. 优化的过程可以设isif=2查看stress,然后确定增大或减小参数。
赞一下 回复此楼
23楼2012-08-14 09:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

魔鬼中的天使

金虫 (小有名气)
★ ★ ★
奋力拼搏(金币+1): 谢谢参与
WDD880227: 金币+1, 多谢指教 2012-08-15 08:59:59
奋力拼搏: 金币+1, 谢谢您的回复！ 2012-08-15 19:43:53
用isif=3就是简单优化
可以先把精度设置低一些（prec、ecut、k-mesh这些参数）
然后逐渐提高精度直至收敛
赞一下 回复此楼
25楼2012-08-14 20:04:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)
引用回帖:
23楼: Originally posted by cclw at 2012-08-14 09:56:41
六角就两个参数a和c，可以分别优化。建议用a作为scale,先固定c/a，优化a,直至能量最低；然后优化c/a. 优化的过程可以设isif=2查看stress,然后确定增大或减小参数。

麻烦您给我细讲一下怎样做才能固定c/a优化a,能量最低之后又怎样做才能优化c/a呢？
麻烦您给我讲一下吧，谢谢
赞一下 回复此楼
27楼2012-08-15 19:42:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaomeiyan

木虫 (正式写手)
★
奋力拼搏(金币+1): 谢谢参与
祝福楼主
回复此楼
28楼2012-08-15 20:18:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaomeiyan

木虫 (正式写手)
★
奋力拼搏(金币+1): 谢谢参与
祝福楼主
回复此楼
31楼2012-08-15 20:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

成广

铜虫 (初入文坛)
★ ★ ★ ★ ★
奋力拼搏(金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
goingabroad2楼
2012-08-12 23:35   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 非水区，拒绝灌水！ 2012-08-13 00:07:11
xueht987: 金币+1, 操作失误，见谅 2012-08-13 14:24:05
zhaowei35633楼
2012-08-13 00:01   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 00:07:31
zhaowei35634楼
2012-08-13 00:01   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 00:07:42
zhaowei35635楼
2012-08-13 00:01   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 00:07:51
zhaowei35636楼
2012-08-13 00:02   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-4, 拒绝灌水！ 2012-08-13 00:08:08
xueht987: 金币+4, 请为本版的发展做出您应有的贡献，大家一起来营造积极、和谐的学术交流氛围！处罚过于严厉了些，故再返4金。 2012-08-13 14:30:37
zhaowei35637楼
2012-08-13 00:02   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-5, 拒绝灌水！ 2012-08-13 00:08:21
坦然随风8楼
2012-08-13 07:10   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-1, 拒绝灌水！ 2012-08-13 14:15:22
坦然随风10楼
2012-08-13 07:10   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 14:15:33
坦然随风11楼
2012-08-13 07:11   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 14:15:47
坦然随风12楼
2012-08-13 07:11   回复  
奋力拼搏(金币+1): 谢谢参与
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
xiaoyi-198816楼
2012-08-13 09:22   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-2, 拒绝灌水！ 2012-08-13 14:16:42
xiaoyi-198817楼
2012-08-13 09:22   回复  
奋力拼搏(金币+1): 谢谢参与
xueht987: 金币-3, 拒绝灌水！ 2012-08-13 14:17:01
xiaoyi-198819楼
2012-08-13 09:22   回复  
奋力拼搏(金币+1): 谢谢参与
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
lihong198924楼
2012-08-14 10:32   回复  
奋力拼搏(金币+1): 谢谢参与
奋力拼搏: 金币+1 2012-08-15 19:43:33
郑克宁26楼
2012-08-14 20:56   回复  
奋力拼搏(金币+1): 谢谢参与
奋力拼搏: 金币+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): 谢谢参与
相关版块跳转 我要订阅楼主 奋力拼搏 的主题更新
351/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭