关于扩胞计算声子谱

jzjy1994

请教一下各位老师，在用castep进行声子谱计算，需要扩胞，但扩胞后晶体对称性会变成P1，倘若加入对称性，晶体又会变成原本单胞的形式。也就是说，倘若我想用超胞计算那只能用P1对称性？？？是因为我超胞扩展方式不对（比方说扩展需要平移？）还是加对称性方式不对？
另外请教一下老师如果用vasp扩胞计算声子谱加对称性会出现这类情况吗=0=。
谢谢各位老师了

卡开发发

castep有两种不同的方法做，线性响应方法是通过单胞直接计算多个q点的，但只能针对非金属；有限差分也有专门设置扩胞的选项，不需要在supercell的选项扩展超胞。vasp针对对称性有ibrion=6和8两个分别对应于考虑对称性的有限差分方法和线性响应方法，因为只能计算Gamma点的，通常需要使用超胞来计算然后做反折叠。

jzjy1994

卡开发发 发表于 2017-10-23 15:06
castep有两种不同的方法做，线性响应方法是通过单胞直接计算多个q点的，但只能针对非金属；有限差分也有专 ...

老师我用的是castep  finite displacement用原胞做出来虚频严重。=。=所以以为和vasp一样要扩胞。也就是我虚频严重的话（0以下有十几条格波）是无法通过扩胞解决的吗。=0=老师，你说的能在castep中扩胞的选项是在哪里？我把声子谱和设置发楼下-0-拜托老师看一下了，谢谢老师C:\Users\Administrator\Desktop\convergence test

jzjy1994

C:\Users\Administrator\Desktop\convergence test

卡开发发

jzjy1994 发表于 2017-10-23 16:30
老师我用的是castep  finite displacement用原胞做出来虚频严重。=。=所以以为和vasp一样要扩胞。也就是 ...

supercell defined by cutoff radius就是，你可以调大这个值。

jzjy1994

=0=谢谢老师，说起来老师已经帮我解答了5个问题了非常感谢老师的帮助，每次发帖总是能收到老师的回复。非常感谢

gog

卡开发发 发表于 2017-10-23 15:06
castep有两种不同的方法做，线性响应方法是通过单胞直接计算多个q点的，但只能针对非金属；有限差分也有专 ...

老师。您好！请教您几个问题
我这用castep算的声子。
线性响应方法，能量最小化报错，没有获得声子谱。
有限位置替换算出来，声子态密度是没有虚频。声子态密度可以改变显示单位，和声子谱单位一致。第一个问题是：声子态密度是声子谱或者声子谱等效谱图？
第二个问题是：声子态密度在零位处态密度值相比其他地方，要高出好几倍，这个怎么解释呢？
第三个问题是：网上有资料说，线性响应适合计算非磁性材料，有限替换法是适合磁性材料，这么讲，线性响应法没获得到声子谱、有限替换法获得了声子谱、声子态密度，我设计的材料该显示有磁性，这样理解正确不呢？
能带和态密度，是用非磁性材料去计算获得的。
如果是磁性材料，能带和态密度，就应该用磁性材料的特征去设置参数重新计算？

卡开发发

gog 发表于 2022-12-16 21:09
老师。您好！请教您几个问题
我这用castep算的声子。
线性响应方法，能量最小化报错，没有获得声子谱。 ...

1、

线性响应方法，能量最小化报错，没有获得声子谱。

极小化报错原因可能要找找，不过castep线性响应只能在模守恒赝势上做确实很麻烦。
2、

声子态密度在零位处态密度值相比其他地方，要高出好几倍，这个怎么解释呢？

这个有可能，但我也没办法确定，比方对于一些简单格子可能确实没有光学支，我对此了解有限，你可以看看其他人的讨论。
3、

网上有资料说，线性响应适合计算非磁性材料，有限替换法是适合磁性材料，这么讲，线性响应法没获得到声子谱、有限替换法获得了声子谱、声子态密度，我设计的材料该显示有磁性，这样理解正确不呢？

极小化报错没办法说明这件事。可能是没办法用于金属性的体系，我印象中对线性响应方法没办法使用展宽。但原理上金属性到底能否使用以及磁性能否使用，这点我没办法回答，得看castep是否实现，同类程序中qe肯定是能这么做，并且不强制模守恒赝势。

卡开发发

gog 发表于 2022-12-16 21:09
老师。您好！请教您几个问题
我这用castep算的声子。
线性响应方法，能量最小化报错，没有获得声子谱。 ...

4、

如果是磁性材料，能带和态密度，就应该用磁性材料的特征去设置参数重新计算？

磁性材料的能带和态密度应该开自旋极化去算，否则得到的结果就是非磁的。

gog

卡开发发 发表于 2022-12-17 10:49
4、
磁性材料的能带和态密度应该开自旋极化去算，否则得到的结果就是非磁的。

非常感谢老师详细解答。

Zim

卡开发发 发表于 2022-12-17 10:49
4、
磁性材料的能带和态密度应该开自旋极化去算，否则得到的结果就是非磁的。

老师您好，很抱歉这么晚打扰您。我想问一下算石墨烯的能带和态密度可以直接用MS软件自带的石墨然后切一层下来（1*1的晶格，里面有两个碳原子）直接算吗？还是说需要建立3*3的结构去算啊？因为我刚刚试了一下用1*1的算出来是正常的，但是用2*2的算出来能带就全是直线了

卡开发发

Zim 发表于 2023-4-19 00:26
老师您好，很抱歉这么晚打扰您。我想问一下算石墨烯的能带和态密度可以直接用MS软件自带的石墨然后切一层 ...

1、MS自动的路径不一定合理。
2、涉及到能带折叠问题。

Zim

卡开发发 发表于 2023-4-19 06:26
1、MS自动的路径不一定合理。
2、涉及到能带折叠问题。

老师，那我如果是算层数比较少石墨烯的能带和态密度的话可以不扩胞直接就用MS上面给的石墨切一层下来建立层数（形成1*1*层数的结构）直接做计算吗？

卡开发发

Zim 发表于 2023-4-19 11:41
老师，那我如果是算层数比较少石墨烯的能带和态密度的话可以不扩胞直接就用MS上面给的石墨切一层下来建立 ...

能带按道理该在相应的原胞上进行计算，态密度具备广延性就无所谓了。

Zim

卡开发发 发表于 2023-4-19 12:05
能带按道理该在相应的原胞上进行计算，态密度具备广延性就无所谓了。

谢谢卡老师