计算化学公社

标题: 关于扩胞计算声子谱 [打印本页]

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

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

作者
Author: jzjy1994 时间: 2017-10-23 16:30

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

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

作者
Author: jzjy1994 时间: 2017-10-23 16:30
本帖最后由 jzjy1994 于 2017-10-23 16:35 编辑

C:\Users\Administrator\Desktop\convergence test

作者
Author: 卡开发发 时间: 2017-10-23 17:00

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

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

作者
Author: jzjy1994 时间: 2017-10-23 17:16
=0=谢谢老师，说起来老师已经帮我解答了5个问题了非常感谢老师的帮助，每次发帖总是能收到老师的回复。非常感谢

作者
Author: gog 时间: 2022-12-16 21:09
本帖最后由 gog 于 2022-12-16 21:19 编辑

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

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

作者
Author: 卡开发发 时间: 2022-12-17 10:46

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

1、

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

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

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

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

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

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

作者
Author: 卡开发发 时间: 2022-12-17 10:49

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

4、

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

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

作者
Author: gog 时间: 2022-12-18 08:36

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

非常感谢老师详细解答。

作者
Author: Zim 时间: 2023-4-19 00:26

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

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

作者
Author: 卡开发发 时间: 2023-4-19 06:26

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

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

作者
Author: Zim 时间: 2023-4-19 11:41

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

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

作者
Author: 卡开发发 时间: 2023-4-19 12:05

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

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

作者
Author: Zim 时间: 2023-4-19 14:14

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

谢谢卡老师

作者
Author: Zim 时间: 2023-4-20 12:48

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

卡老师那石墨烯计算的时候要加自旋极化吗？还有就是我画出来的态密度图和别人的不一样是因为参数选择不一样吗？我需要怎么判断我的态密度图正不正确呢

作者
Author: 卡开发发 时间: 2023-4-20 14:31

Zim 发表于 2023-4-20 12:48
卡老师那石墨烯计算的时候要加自旋极化吗？还有就是我画出来的态密度图和别人的不一样是因为参数选择不一 ...

完整的纯石墨烯应该没有自旋极化这个问题吧。
你可能需要把你计算的细节和你参考资料的细节给出来。

作者
Author: Zim 时间: 2023-4-20 16:00

卡开发发发表于 2023-4-20 14:31
完整的纯石墨烯应该没有自旋极化这个问题吧。
你可能需要把你计算的细节和你参考资料的细节给出来。

我是准备算石墨烯的能带，然后我用的是石墨烯单胞。需要勾选自旋极化吗

作者
Author: 卡开发发 时间: 2023-4-20 19:03

Zim 发表于 2023-4-20 16:00
我是准备算石墨烯的能带，然后我用的是石墨烯单胞。需要勾选自旋极化吗

不需要。

作者
Author: Zim 时间: 2023-4-20 19:24

卡开发发发表于 2023-4-20 19:03
不需要。

好的，非常感谢卡老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)