当前位置: 首页 > 第一原理 >VASP-DFPT & phonopy calculation 计算声子谱求助

VASP-DFPT & phonopy calculation 计算声子谱求助

作者 yangna188
来源: 小木虫 500 10 举报帖子
+关注

1. 第一次计算声子谱,想重复vasp手册上面的NaCl例子,但是做出来的图线有点多,
以下是我的band,conf里面的参数
ATOM_NAME = Na Cl                           
FORCE_CONSTANTS = READ
DIM = 2 2 2                                 
BAND = 0.00 0.00 0.00 0.50 0.00 0.00 0.50 0.50 0.00 0.00 0.00 0.00 0.50 0.50 0.50 波矢路径
BAND_POINTS = 200
INCAR里面的参数如下:
PREC = Accurate
  ENCUT = 500
IBRION = 8
  EDIFF = 1.0e-08
  IALGO = 38
ISMEAR = 0; SIGMA = 0.1
  LREAL = .FALSE.
ADDGRID = .TRUE.
  LWAVE = .FALSE.
LCHARG = .FALSE.
2. 文献上面说如果晶体中有N个原胞,每个原胞中有n个原子,那么一共有3nN种不同频率模式。这3nN个独立振动分为3n支,其中3支是声学波,其余3(n-1)支是光学波,请问波支是什么,是波的条数?还是什么?还有声学波与光学波怎么区分?

VASP-DFPT & phonopy calculation 计算声子谱求助
TZp.png

  • 精华评论
  • zyl雨田

    算例中的NaCl 一个原胞里面有三个原子吧,所以有6个振动模式。但是你这个怎么会这么多振动模式?看输入参数没看出有什么问题?NaCl 分不分原胞跟惯用胞?
    声子谱最后画出来的是原胞的色散关系 所以振动模式应该是3n 支( n表示原胞中的原子数)至于怎么区分声学波 光学波 建议你看看黄昆的固体物理  晶格振动那一章

  • yangna188

    引用回帖:
    2楼: Originally posted by zyl雨田 at 2018-01-08 16:50:37
    算例中的NaCl 一个原胞里面有三个原子吧,所以有6个振动模式。但是你这个怎么会这么多振动模式?看输入参数没看出有什么问题?NaCl 分不分原胞跟惯用胞?
    声子谱最后画出来的是原胞的色散关系 所以振动模式应该是 ...

    我的原胞和教程里面的一模一样,就是如下:
    Na Cl                                   
       1.00000000000000     
         5.6903014761756712    0.0000000000000000    0.0000000000000000
         0.0000000000000000    5.6903014761756712    0.0000000000000000
         0.0000000000000000    0.0000000000000000    5.6903014761756712
       Na   Cl
         4     4
    Direct
      0.0000000000000000  0.0000000000000000  0.0000000000000000
      0.0000000000000000  0.5000000000000000  0.5000000000000000
      0.5000000000000000  0.0000000000000000  0.5000000000000000
      0.5000000000000000  0.5000000000000000  0.0000000000000000
      0.5000000000000000  0.5000000000000000  0.5000000000000000
      0.5000000000000000  0.0000000000000000  0.0000000000000000
      0.0000000000000000  0.5000000000000000  0.0000000000000000
      0.0000000000000000  0.0000000000000000  0.5000000000000000
    还有,您说的NaCl胞有三个原子,那应该画出来是九条线,而不是六条呀

  • 小虫虫LG

    引用回帖:
    3楼: Originally posted by yangna188 at 2018-01-08 16:57:08
    我的原胞和教程里面的一模一样,就是如下:
    Na Cl                                   
       1.00000000000000     
         5.6903014761756712    0.0000000000000000    0.0000000000000000
         0.000000000000 ...

    这不得了,你用的是惯用胞。8个原子。

  • yangna188

    引用回帖:
    5楼: Originally posted by zyl雨田 at 2018-01-09 10:32:50
    啊我说错了,不好意思啊,是两个原子。楼下有回答你了,你的POScar不是有8个原子吗?算出来有24条振动模式,你数数你画出来的振动模式是不是有这么多。我这边没有phonopy的手册,所以你需要确定一下你用的POScar跟算 ...

    第一次计算声子谱,所以问题比较小白:
    1. PRIMITIVE这个参数是在band.conf中设置吗,加入PRIMITIVE= 0.0 0.5 0.5  0.5 0.0  0.5 0.5 0.5 0.0,后执行% phonopy --dim="2 2 2" -c POSCAR-unitcell band.conf,老是报错。
    2. 还有POSCAR中初胞原子的坐标格式也是上面的POSCAR一致吗?只是相应的原子数减小,坐标个数减少?
    3. 最主要一点,第一次求助,如何给你金币

  • zyl雨田

    你应该根据phonopy手册上面的算例来输入命令。PRIMITIVE是将惯用胞转换成原胞的矩阵,这个矩阵你是否有输对?  另外我知道为什么处理band.conf的时候还需要用到--dim这个扩胞命令?  其次 我不知道你第二个问题是什么意思  最后  不用给我金币了   好好看看手册,根据手册上面说的一步步来应该没有问题才是啊。

应助之星
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓