24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者zz11ss11zz将赠送您 50 个金币

zz11ss11zz

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 4666.4
散金: 50
红花: 1
帖子: 982
在线: 49.6小时
虫号: 7166386
注册: 2017-09-11
专业: 凝聚态物性 II ：电子结构

[求助] EPW phonon linewidth inconsistent with QE results

大家好，最近在用epw研究CDW问题，虽然epw算出的band与phonon和qe结果完美吻合，但是算出的声子线宽有很大差别如下图，求问如何解决这个矛盾：
epw.in：

CODE:

--

&inputepw

  prefix      = 'TaSe2'

  amass(1)    = 180.947

  amass(2)    = 78.96

  outdir      = './tmp/'

  iverbosity  = 3

  elph        = .true.

  epbwrite    = .true.

  epbread     = .false.

  epwwrite    = .true.

  epwread     = .false.

  nbndsub     =  11

  nbndskip    =  13

  wannierize  = .true.

  num_iter    = 300

  iprint      = 2

  dis_win_max = 6.5

  dis_win_min = -5

  dis_froz_max = 4

  dis_froz_min = -4

  proj(1)     = 'Ta:l=2'   

  proj(2)     = 'Se:l=1'

  wdata(1) = 'bands_plot = .true.'

  wdata(2) = 'begin kpoint_path'

  wdata(3) = 'G 0.00 0.00 0.00 M 0.00 0.50 0.00'

  wdata(4) = 'M 0.00 0.50 0.00 K 0.333 0.333 0.00'

  wdata(5) = 'K 0.333 0.333 0.00 G 0.00 0.00 0.00'

  wdata(6) = 'end kpoint_path'

  wdata(7) = 'bands_num_points = 30'

  wdata(8) = 'dis_num_iter = 10000'

  wdata(9) = 'guiding_centres = .true.'

  elecselfen  = .false.

  phonselfen  = .false.

  a2f         = .false.

  fsthick     = 3 ! eV 

  eptemp      = 30 ! K (same as PRB 76, 165108)

  degaussw    = 0.005 ! eV

  dvscf_dir   = '../../phonon/save'

!  filukk      = './diam.ukk'

!  filqf       = 'meshes/path.dat'

  nkf1        = 160

  nkf2        = 160

  nkf3        = 1

  nqf1        = 80

  nqf2        = 80

  nqf3        = 1

  nk1         = 16

  nk2         = 16

  nk3         = 1

  nq1         = 8

  nq2         = 8

  nq3         = 1

 /

      10 cartesian

   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00

   0.000000000000000E+00   0.144337567297387E+00   0.000000000000000E+00

   0.000000000000000E+00   0.288675134594775E+00   0.000000000000000E+00

   0.000000000000000E+00   0.433012701892162E+00   0.000000000000000E+00

   0.000000000000000E+00  -0.577350269189550E+00   0.000000000000000E+00

   0.124999999999992E+00   0.216506350946081E+00   0.000000000000000E+00

   0.124999999999992E+00   0.360843918243469E+00   0.000000000000000E+00

   0.124999999999992E+00   0.505181485540856E+00   0.000000000000000E+00

   0.249999999999985E+00   0.433012701892162E+00   0.000000000000000E+00

   0.249999999999985E+00   0.577350269189550E+00   0.000000000000000E+00

epw-linewidth.in

CODE:

--

&inputepw

  prefix      = 'TaSe2'

  amass(1)    = 180.947

  amass(2)    = 78.96

  outdir      = './tmp/'

  iverbosity  = 3

  elph        = .true.

  kmaps       = .true.

  epbwrite    = .false.

  epbread     = .false.

  epwwrite    = .false.

  epwread     = .true.

  nbndsub     =  11

  nbndskip    =  13

  wannierize  = .false.

  num_iter    = 300

  iprint      = 2

  dis_win_max = 6.5

  dis_win_min = -5

  dis_froz_max = 4

  dis_froz_min = -4

  proj(1)     = 'Ta:l=2'   

  proj(2)     = 'Se:l=1'

  wdata(1) = 'bands_plot = .true.'

  wdata(2) = 'begin kpoint_path'

  wdata(3) = 'G 0.00 0.00 0.00 M 0.00 0.50 0.00'

  wdata(4) = 'M 0.00 0.50 0.00 K 0.333 0.333 0.00'

  wdata(5) = 'K 0.333 0.333 0.00 G 0.00 0.00 0.00'

  wdata(6) = 'end kpoint_path'

  wdata(7) = 'bands_num_points = 30'

  wdata(8) = 'dis_num_iter = 10000'

  wdata(9) = 'guiding_centres = .true.'

  elecselfen  = .false.

  phonselfen  = .true.

  a2f         = .false.

  delta_approx= .true.

  fsthick     = 3 ! eV 

  eptemp      = 1 ! K (same as PRB 76, 165108)

  degaussw    = 0.1 ! eV

  dvscf_dir   = '../../phonon/save'

!  filukk      = './diam.ukk'

  filqf       = 'meshes/path.dat'

  nkf1        = 160

  nkf2        = 160

  nkf3        = 1

!  nqf1        = 80

!  nqf2        = 80

!  nqf3        = 1

  nk1         = 16

  nk2         = 16

  nk3         = 1

  nq1         = 8

  nq2         = 8

  nq3         = 1

 /

      10 cartesian

   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00

   0.000000000000000E+00   0.144337567297387E+00   0.000000000000000E+00

   0.000000000000000E+00   0.288675134594775E+00   0.000000000000000E+00

   0.000000000000000E+00   0.433012701892162E+00   0.000000000000000E+00

   0.000000000000000E+00  -0.577350269189550E+00   0.000000000000000E+00

   0.124999999999992E+00   0.216506350946081E+00   0.000000000000000E+00

   0.124999999999992E+00   0.360843918243469E+00   0.000000000000000E+00

   0.124999999999992E+00   0.505181485540856E+00   0.000000000000000E+00

   0.249999999999985E+00   0.433012701892162E+00   0.000000000000000E+00

   0.249999999999985E+00   0.577350269189550E+00   0.000000000000000E+00