卡开发发 发表于 2024-12-27 00:45 感谢卡卡老师百忙之中抽空回答! |
going6667 发表于 2024-12-26 23:08 1、金属仍然可能需要,甚至可以说对于小晶胞的情况绝大多数情况都得扩。只是有可能不见得要那么大的超胞,其他情况有可能需要的超胞更大,按到也需要通过足够的测试,不过因为声子的计算量太大大家没人愿意这样做。除此之外,即便对DFPT来说VASP也没办法处理声子BZ采样q点为Gamma以上的情况,得通过超晶胞的声子谱进行折叠。 2、有限位移没有这方面限制吧,但是确实需要超晶胞上面已经提到,因为原子位移导致的密度变化也可以看成是一种响应。DFPT原则上可以不扩胞,这方面大概是VASP自身算法的原因,具体如何实现我没有深究,如果有人能讨论得清楚可以进行补充。 3、严格说得测试,科学的角度说主观经验不太可取,但碍于算不动,没人愿意这么做。 |
参与人数Participants 1 | eV +3 | 收起 理由Reason |
---|---|---|
| + 3 | 谢谢 |
卡开发发 发表于 2024-10-24 17:46 卡卡老师,您好。我是一名刚刚接触声子色散计算的学生,我从很多途径中了解到声子色散计算需要构建超胞,从vaspwiki了解到声子色散需要构建超胞的原因是:半导体和绝缘体的电子屏蔽作用比较弱,需要考虑长程相互作用对力常数矩阵的影响。目前我有几个比较简单的问题一直想不出答案,能麻烦您指导一下吗: (1)声子色散计算需要构建超胞的原因只是因为vaspwiki上说的这个原因吗?那么金属体系是不是不需要进行超胞? (2)vasp+phonopy计算动力学矩阵有有限位移法和DFPT两种方法,我了解到有限位移法要求声子波矢与原胞边界supersize正交,而DFPT原则上不需要扩胞,但因为vasp原因,计算DFPT还是需要扩胞,请问我了解的地方有错误的地方吗? (3)上述两种方法扩胞大小有经验条件吗?我从一些DFT数据库网站上了解到声子谱计算扩胞后原则上不应小于12埃,到20埃比较保险。是不是还得去测试超胞大小直到声子谱变化收敛? |
卡开发发 发表于 2024-10-24 17:46 好的,感谢! |
松鼠 发表于 2024-10-24 15:37 IBRION=-1本身没办法算有限差分,只是通过外部程序例如Pymatgen或者ASE产生一批扩超晶格原子位移的结构,然后通过IBRION=-1+NSW=0进行单点计算获得力,然后再通过后处理得到Hessian矩阵。IBRION=5其实相当于上面的过程全部在VASP内部进行,IBRION=6用DFPT替代有限差分。 |
卡开发发 发表于 2021-11-20 21:05 请问一下有限差分法用IBRION=-1和用IBRION=5或6具体的区别是啥吗? |
本帖最后由 松鼠 于 2024-10-24 15:30 编辑 请问一下有限差分法用IBRION=-1和用IBRION=5或6具体的区别是啥啊? |
谢谢老师 |
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“用DFPT的方法计算声子谱”改了,以后务必注意 |
ShuaiCui 发表于 2021-11-20 21:13 声子谱计算与优化过程的计算方法和精度一致即可。 |
卡开发发 发表于 2021-11-20 21:05 那我还需要重新结构优化嘛? |
VASP应该没做DFPT部分对IVDW的支持,你要不换有限差分做(IBRION=5或6)。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-17 00:52 , Processed in 1.677348 second(s), 26 queries , Gzip On.