计算化学公社

标题: 关于声子谱的一些问题 [打印本页]

作者
Author: jzjy1994 时间: 2017-10-17 15:58
标题: 关于声子谱的一些问题
本帖最后由 jzjy1994 于 2017-10-17 17:10 编辑

首先第一个问题是声子谱中格波的数量和原胞的原子数有关。那么用超胞算或者某一个面是会产生声子折叠？？？还有超胞计算的时间是要比原胞长很多吗
第二个问题是用castep结构要收敛到何种程度才能用于声子计算
第三个问题就是声子计算的时间了，比方说28个原子，18个k点，我用18个核来计算大概要多久
最后一个问题是比方说我要计算A/B（各2层）超胞的声子，那我想看B声子谱在结合A之后的变化，是要计算一个4层的B的超胞吗？
在这里向各位老师请教
非常感谢

作者
Author: jzjy1994 时间: 2017-10-17 17:14
另外想再请教一下各位老师一个问题，如用石墨烯和和金属Ni的111面进行结合，那么理论上来讲这个结构的布里渊区是要改变的，那么在对结合前和结合后石墨烯的能带和声子怎么进行比较呢，因为k路径都变化了。

作者
Author: 卡开发发 时间: 2017-10-17 19:16
1、会产生色散曲线的折叠。超胞比原胞计算确实慢很多，对于有限差分，首先等同于计算6N的单点能，然后每个单点能和体系大小大致也是O(N^3)的，但是有限差分的方法必须使用超胞，线性响应方法我不太清楚具体的标度。不过你后面提到的石墨烯和Ni复合的问题应该是金属性倾向的体系，因为Fermi面附近的能级近简并，可能一般的线性响应方法会有点问题，我相信castep本身也会给出同样的warning。
2、原则上这个没有限制，但总归不能太宽松。至于程序有时候确实对一些优化好的结构也同样warning，我不太清楚程序怎么判断的。
3、有限差分标度在问题1的回答我已经给出来了，但每个单点算多快还得看你服务器的配置，你可以估计一下计算时间的。不过对于Ni这种有磁性的体系，有限差分过程有些scf可能会很长甚至难收敛。
4、你需要把两个layer复合到一块再计算，而且还涉及到匹配性的问题，对于其中至少一个layer可能不是单胞。一样的六方格子，路径应该差别不大，复合之后涉及到折叠程度不同的问题，应该没办法参照各自单胞的情形。

作者
Author: jzjy1994 时间: 2017-10-17 20:07

卡开发发发表于 2017-10-17 19:16
1、会产生色散曲线的折叠。超胞比原胞计算确实慢很多，对于有限差分，首先等同于计算6N的单点能，然后每个 ...

谢谢版主的耐心回答，因为在文献中有看到这种比较方式，但是不清楚怎么操作，这是文献中的图C:\Users\Administrator\Desktop\convergence test

作者
Author: 卡开发发 时间: 2017-10-17 20:43

jzjy1994 发表于 2017-10-17 20:07
谢谢版主的耐心回答，因为在文献中有看到这种比较方式，但是不清楚怎么操作，这是文献中的图

我看了下原因是Ni111和石墨烯的晶格匹配比较接近，所以可以用单胞的layer直接摞在一块就可以。这样直接算石墨烯单胞在Ni111单胞晶格尺寸下和反过来即可。如果需要晶格常数公倍数才能重叠的两个layer，用单胞就没办法解决了。

作者
Author: jzjy1994 时间: 2017-10-17 21:10
版主的意思是用Ni111单胞的晶格尺寸来作为石墨烯的尺寸，然后在此基础上再计算单是石墨烯的声子谱。再拿这个声子谱和结合后的声子谱相比较吗？文献中貌似用了2层的石墨烯，中间有三层的Ni在这种状况下原胞内的石墨烯原子数也和会比单胞多出一倍，这样声子模式不会变多吗？还是说是因为他选择的k路径只考虑了一个平面的路径？（另外结合之后布里渊区应该变化了吧，那路径=。=应该也变化了吧）

作者
Author: 卡开发发 时间: 2017-10-17 21:46

jzjy1994 发表于 2017-10-17 21:10
版主的意思是用Ni111单胞的晶格尺寸来作为石墨烯的尺寸，然后在此基础上再计算单是石墨烯的声子谱。再拿这 ...

是的，在Ni111的晶格下面分别算石墨烯、Ni111和结合的，再是石墨烯晶格下面分别算石墨烯、Ni111和结合的。声子模式因为原子数目变多肯定会变多。这个体系本来就是2D的系统，上下应该是敞开的，这样的路径才会是G-M-K-G，因为两者都是六方格子，所以q点的路径不会不同。

作者
Author: jzjy1994 时间: 2017-10-17 21:54
谢谢版主的耐心回答，非常感谢

作者
Author: hf_hua 时间: 2021-8-8 16:16
您好lz，想要请教一下声子折叠问题。就是如何查看折叠，就是比如扩胞之后再计算声子，他的谱线会增加，格波由于多次折射或者反射什么的，就会和未扩胞的声子谱有一些不同，那么怎么能够一一对应起来呐？就是如何找到想对应的模式呐

作者
Author: hf_hua 时间: 2021-9-1 17:11

jzjy1994 发表于 2017-10-17 20:07
谢谢版主的耐心回答，因为在文献中有看到这种比较方式，但是不清楚怎么操作，这是文献中的图

你好，请问这图中，gamma点处的简并点是如何理解呐？由于某种对称性会导致简并的情况？

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