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

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

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

谢谢版主的耐心回答，非常感谢

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

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

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

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

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

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

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