24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4052  |  回复: 25

tinkie0411

银虫 (小有名气)

送红花一朵
引用回帖:
9楼: Originally posted by chuanghua304 at 2014-04-01 08:48:16
其实理论理论上两种晶胞做优化都可以,手册的例子里好像就是用的立方结构,但是要在band.conf里设置,primitive axial。用原胞也可以,那就不用设置primitive axial了。具体体系不同扩展的倍数也不同,你128个原子应 ...

是的 128原子是用立方基扩的 16x2x2x2=128 如果换成3x3x3的话16x3x3x3=432个原子 程序太大了 跑不起来 我也试过用primitive cell 的基矢扩成3x3x3 的4x3x3x3=108个原子 结果差别不大
我用的是PAW-GGA91 Ecut最大值是282 我取的500ev
好的 以后的计算我会把EDIFFG = -1.0e-06 和NPAR=4删除 我现在再用primitive cell扩成4x4x4的超胞 256个原子 已经跑的很慢很慢了 而且我
也不确定算这么长时间最后能否有个可以用的结果 很困惑
10楼dx0620说DFPT只能算Gamma点 你怎么看 还有你告诉我的IBRION=6的差分法 是不是优化完之后 再算超胞的时候加这个 这样的算法超胞计算只需要一次 我的理解对吗?非常感谢!
11楼2014-04-01 16:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)

引用回帖:
10楼: Originally posted by dx0620 at 2014-04-01 09:19:29
我怎么感觉出现这样的虚频不是计算方法的问题啊,你用的什么赝势?如果出现这种问题,建议先不要用考虑半芯态赝势。你如果惯用胞的超胞只有128个原子的话,就用惯用胞吧,对原胞的设置如果不熟悉的话很容易出错。据 ...

感谢你的回复!我用的是PAW-GGA91,你是说是赝势的问题吗?那我试试PAW-PBE。惯用胞2x2x2的超胞是128个原子,但是我不确定我这种离奇的虚频是否是因为超胞不够大而引起的。我也试过用primitive cell的基矢扩胞扩成4x3x3x3=108个原子的,也有这种虚频的情况。我现在用的是直接位移法,也打算试试DFPT和6楼chuanghua说的差分法。我现在最大的问题就是不知道产生这种离奇虚频的原因,所以就只能不停地尝试各种方法但也不确定能否解决。
12楼2014-04-01 16:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)

引用回帖:
10楼: Originally posted by dx0620 at 2014-04-01 09:19:29
我怎么感觉出现这样的虚频不是计算方法的问题啊,你用的什么赝势?如果出现这种问题,建议先不要用考虑半芯态赝势。你如果惯用胞的超胞只有128个原子的话,就用惯用胞吧,对原胞的设置如果不熟悉的话很容易出错。据 ...

我看phonopy manual 1.7.4手册上62页用DFPT方法画的NaCl声子谱与第4页直接位移法是一致的 是不是DFPT方法算的也不只是Gamma点了呢 期待您的回复 谢谢!
13楼2014-04-01 18:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
12楼: Originally posted by tinkie0411 at 2014-04-01 16:32:17
感谢你的回复!我用的是PAW-GGA91,你是说是赝势的问题吗?那我试试PAW-PBE。惯用胞2x2x2的超胞是128个原子,但是我不确定我这种离奇的虚频是否是因为超胞不够大而引起的。我也试过用primitive cell的基矢扩胞扩成 ...

你就用惯用胞2x2x2的超胞就行,IBRION=-1。我还是觉得出现你这种虚频的话跟赝势有很大关系,我之前算过的就有这种情况,你可以看看,更奇怪。赝势你都先别用考虑半芯态的
Vasp+phonopy 声子谱出现离奇虚频!

» 本帖已获得的红花(最新10朵)

14楼2014-04-01 19:48:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
13楼: Originally posted by tinkie0411 at 2014-04-01 18:00:34
我看phonopy manual 1.7.4手册上62页用DFPT方法画的NaCl声子谱与第4页直接位移法是一致的 是不是DFPT方法算的也不只是Gamma点了呢 期待您的回复 谢谢!...

NaCl对称性那么高,算的声子谱跟前面的一样很正常的,想出问题都难。1.7.4手册上第80页说的很清楚,只能计算Gamma点上的。
15楼2014-04-01 19:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)

送红花一朵
引用回帖:
14楼: Originally posted by dx0620 at 2014-04-01 19:48:35
你就用惯用胞2x2x2的超胞就行,IBRION=-1。我还是觉得出现你这种虚频的话跟赝势有很大关系,我之前算过的就有这种情况,你可以看看,更奇怪。赝势你都先别用考虑半芯态的
...

非常感谢你的回复!你说的半芯态赝势是什么意思,是不是指赝势文件中带—sv —pv 和—d的赝势呢?
我找到的解释
有些轨道不属于价轨道,但是能量又不是太低(通常位于HOMO以下5-10eV的范围)。如果在post-HF计算中做冻结处理,会对分子结构和特性通常有不太大也不太小的影响(对有些量,影响是非常大的,比如四极矩)。这些轨道就是semi-core轨道,其中的电子称为semi-core电子。例如3d金属的3s3p电子,碘的4d电子,金的4f5s5p电子。
如果追求较高的精度,需要对这部分电子进行相关能处理,而不能冻结。当然对于HF计算,这么做是没有任何意义的。 是这样理解吗?
我现在做的是Ni2MnGa,画的是fcc奥氏体的,空间群225,按理说对称性也是相当好啊,不知道为什么会出现这种情况。POSCAR我在1楼贴出来了,POTCAR我用的是PAW-GGA91,Ni10Mn7Ga13 就是取的最简单的那种(Ga考虑了d),但是我也做过Ga3的,也出现这种情况,这应该没用你说的考虑半芯态的赝势吧? 顺带一句,我也试过USPP-GGA91,结果一样。我现在打算试一下PAW-GGAPBE,如果还是有这种问题,我不知道还能换什么赝势了。另外我最近做了四方结构马氏体相的声子谱,就没有这种离奇虚频出现了。虽说奥氏体相低温不是稳定相,但是我现在做的是重复很多已经发表的文章,人家的都没有我这种情况,我到底错哪了呢???
16楼2014-04-01 21:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)

引用回帖:
15楼: Originally posted by dx0620 at 2014-04-01 19:55:03
NaCl对称性那么高,算的声子谱跟前面的一样很正常的,想出问题都难。1.7.4手册上第80页说的很清楚,只能计算Gamma点上的。...

好的,看到了,多谢!
17楼2014-04-01 21:27:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
16楼: Originally posted by tinkie0411 at 2014-04-01 21:25:51
非常感谢你的回复!你说的半芯态赝势是什么意思,是不是指赝势文件中带—sv —pv 和—d的赝势呢?
我找到的解释
有些轨道不属于价轨道,但是能量又不是太低(通常位于HOMO以下5-10eV的范围)。如果在post-HF计算 ...

差不多这么理解,简单来说就是不单单只考虑最外侧电子。那你先试试PBE的赝势看看吧
18楼2014-04-01 21:35:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)

引用回帖:
14楼: Originally posted by dx0620 at 2014-04-01 19:48:35
你就用惯用胞2x2x2的超胞就行,IBRION=-1。我还是觉得出现你这种虚频的话跟赝势有很大关系,我之前算过的就有这种情况,你可以看看,更奇怪。赝势你都先别用考虑半芯态的
...

这是文献中的声子谱,我画的是Gamma到X那段。奥氏体是不稳定,但是虚频我不是我画出来的那种。他们用的软件是PHONON software,我用的phonopy,虽说前者是收费的可能更好用,但是phonopy也发了很多好文章。我现在的情况总不会是因为软件的问题吧?很困惑!
Vasp+phonopy 声子谱出现离奇虚频!-1
QQ图片20140401153253.jpg

19楼2014-04-01 21:35:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tinkie0411

银虫 (小有名气)

引用回帖:
18楼: Originally posted by dx0620 at 2014-04-01 21:35:13
差不多这么理解,简单来说就是不单单只考虑最外侧电子。那你先试试PBE的赝势看看吧...

回复的好快!实在抱歉跟你抱怨一大堆,这件事情真是让我困惑了好久!我会试试PBE的,有结果了我还会贴上来,衷心感谢!!!
20楼2014-04-01 21:40:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tinkie0411 的主题更新
信息提示
请填处理意见