当前位置: 首页 > 第一原理 >请教vasp+phonopy声子谱虚频问题

请教vasp+phonopy声子谱虚频问题

作者 幻恋小雨
来源: 小木虫 250 5 举报帖子
+关注

计算的是YAlO3出现了下面的情况
虚频从112开始太多了这里就贴119,120
band.conf
ATOM_NAME= Y Al O
DIM = 2 2 2
FORCE_CONSTANTS=READ
BAND = 0 0 0 0 0 0.5 -0.5 0 0.5 -0.5 0 0 -0.5 0.5 0 0 0.5 0 0 0.5 0.5 -0.5 0.5 0.5
BAND_POINTS = 101
BAND_LABELS = \Gamma Z T Y S X U R
对具体的虚频对应的坐标      
以及振动方向来调整原子坐标从而消除虚频不清楚
麻烦大家帮忙看看

119 f/i=    7.818909 THz    49.127654 2PiTHz  260.810722 cm-1    32.336420 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.002034   -0.000001    0.000000  
      3.640000  0.000000  0.000000    -0.002034    0.000001    0.000000  
      0.000000  3.640000  0.000000    -0.002034    0.000001    0.000000  
      3.640000  3.640000  0.000000     0.002034   -0.000001    0.000000  
      0.000000  0.000000  3.650000    -0.002034    0.000001    0.000000  
      3.640000  0.000000  3.650000     0.002034   -0.000001    0.000000  
      0.000000  3.640000  3.650000     0.002034   -0.000001    0.000000  
      3.640000  3.640000  3.650000    -0.002034    0.000001    0.000000  
      1.820000  1.820000  1.825000     0.000000    0.000000    0.000000  
      5.460000  1.820000  1.825000     0.000000    0.000000    0.000000  
      1.820000  5.460000  1.825000     0.000000    0.000000    0.000000  
      5.460000  5.460000  1.825000     0.000000    0.000000    0.000000  
      1.820000  1.820000  5.475000     0.000000    0.000000    0.000000  
      5.460000  1.820000  5.475000     0.000000    0.000000    0.000000  
      1.820000  5.460000  5.475000     0.000000    0.000000    0.000000  
      5.460000  5.460000  5.475000     0.000000    0.000000    0.000000  
      1.820000  1.820000  0.000000    -0.000153    0.250583    0.000000  
      5.460000  1.820000  0.000000     0.000153   -0.250583    0.000000  
      1.820000  5.460000  0.000000     0.000153   -0.250583    0.000000  
      5.460000  5.460000  0.000000    -0.000153    0.250583    0.000000  
      1.820000  1.820000  3.650000     0.000153   -0.250583    0.000000  
      5.460000  1.820000  3.650000    -0.000153    0.250583    0.000000  
      1.820000  5.460000  3.650000    -0.000153    0.250583    0.000000  
      5.460000  5.460000  3.650000     0.000153   -0.250583    0.000000  
      1.820000  0.000000  1.825000     0.000000    0.000000   -0.249407  
      5.460000  0.000000  1.825000     0.000000    0.000000    0.249407  
      1.820000  3.640000  1.825000     0.000000    0.000000    0.249407  
      5.460000  3.640000  1.825000     0.000000    0.000000   -0.249407  
      1.820000  0.000000  5.475000     0.000000    0.000000    0.249407  
      5.460000  0.000000  5.475000     0.000000    0.000000   -0.249407  
      1.820000  3.640000  5.475000     0.000000    0.000000   -0.249407  
      5.460000  3.640000  5.475000     0.000000    0.000000    0.249407  
      0.000000  1.820000  1.825000     0.000000    0.000000    0.000152  
      3.640000  1.820000  1.825000     0.000000    0.000000   -0.000152  
      0.000000  5.460000  1.825000     0.000000    0.000000   -0.000152  
      3.640000  5.460000  1.825000     0.000000    0.000000    0.000152  
      0.000000  1.820000  5.475000     0.000000    0.000000   -0.000152  
      3.640000  1.820000  5.475000     0.000000    0.000000    0.000152  
      0.000000  5.460000  5.475000     0.000000    0.000000    0.000152  
      3.640000  5.460000  5.475000     0.000000    0.000000   -0.000152  

120 f/i=    8.303368 THz    52.171601 2PiTHz  276.970541 cm-1    34.339983 meV
             X         Y         Z           dx          dy          dz
      0.000000  0.000000  0.000000     0.000000    0.000000    0.000000  
      3.640000  0.000000  0.000000     0.000000    0.000000    0.000000  
      0.000000  3.640000  0.000000     0.000000    0.000000    0.000000  
      3.640000  3.640000  0.000000     0.000000    0.000000    0.000000  
      0.000000  0.000000  3.650000     0.000000    0.000000    0.000000  
      3.640000  0.000000  3.650000     0.000000    0.000000    0.000000  
      0.000000  3.640000  3.650000     0.000000    0.000000    0.000000  
      3.640000  3.640000  3.650000     0.000000    0.000000    0.000000  
      1.820000  1.820000  1.825000     0.000000    0.000000    0.000000  
      5.460000  1.820000  1.825000     0.000000    0.000000    0.000000  
      1.820000  5.460000  1.825000     0.000000    0.000000    0.000000  
      5.460000  5.460000  1.825000     0.000000    0.000000    0.000000  
      1.820000  1.820000  5.475000     0.000000    0.000000    0.000000  
      5.460000  1.820000  5.475000     0.000000    0.000000    0.000000  
      1.820000  5.460000  5.475000     0.000000    0.000000    0.000000  
      5.460000  5.460000  5.475000     0.000000    0.000000    0.000000  
      1.820000  1.820000  0.000000     0.000000    0.000000    0.000000  
      5.460000  1.820000  0.000000     0.000000    0.000000    0.000000  
      1.820000  5.460000  0.000000     0.000000    0.000000    0.000000  
      5.460000  5.460000  0.000000     0.000000    0.000000    0.000000  
      1.820000  1.820000  3.650000     0.000000    0.000000    0.000000  
      5.460000  1.820000  3.650000     0.000000    0.000000    0.000000  
      1.820000  5.460000  3.650000     0.000000    0.000000    0.000000  
      5.460000  5.460000  3.650000     0.000000    0.000000    0.000000  
      1.820000  0.000000  1.825000     0.250000    0.000000    0.000000  
      5.460000  0.000000  1.825000    -0.250000    0.000000    0.000000  
      1.820000  3.640000  1.825000    -0.250000    0.000000    0.000000  
      5.460000  3.640000  1.825000     0.250000    0.000000    0.000000  
      1.820000  0.000000  5.475000    -0.250000    0.000000    0.000000  
      5.460000  0.000000  5.475000     0.250000    0.000000    0.000000  
      1.820000  3.640000  5.475000     0.250000    0.000000    0.000000  
      5.460000  3.640000  5.475000    -0.250000    0.000000    0.000000  
      0.000000  1.820000  1.825000     0.000000   -0.250000    0.000000  
      3.640000  1.820000  1.825000     0.000000    0.250000    0.000000  
      0.000000  5.460000  1.825000     0.000000    0.250000    0.000000  
      3.640000  5.460000  1.825000     0.000000   -0.250000    0.000000  
      0.000000  1.820000  5.475000     0.000000    0.250000    0.000000  
      3.640000  1.820000  5.475000     0.000000   -0.250000    0.000000  
      0.000000  5.460000  5.475000     0.000000   -0.250000    0.000000  
      3.640000  5.460000  5.475000     0.000000    0.250000    0.000000

请教vasp+phonopy声子谱虚频问题
1.png 返回小木虫查看更多

今日热帖
  • 精华评论
  • 幻恋小雨

    自己先顶一下

  • zyl雨田

    贴出来的是G点虚频吧。
    估计第一步的结构优化就没优化好

  • 幻恋小雨

    引用回帖:
    3楼: Originally posted by zyl雨田 at 2015-06-03 10:00:41
    贴出来的是G点虚频吧。
    估计第一步的结构优化就没优化好

    是g点虚频
    不过怎么找虚频对应的原子坐标呢

  • Samfain

    厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害厉害

  • prince1000

    什么叫g点虚频,从哪看出来是g点呢,不是很多的特殊k点的地方都有虚频么

猜你喜欢