EPW Forum

Posted: **Mon Jun 17, 2019 1:25 am**

Hello Everyone,
I am a beginner of epw. A few days ago I tried the exercise1 in Wed.4. (http://epw.org.uk/School2018/School2018 ... .Verdi.pdf), and it ran successfully.
However, when I changed the calculated material from Pb to 2D-As, the calculation stopped at the last step. In the epw4.out file, the last content is "Using As.ukk from disk".
My input files:

scf.in

Code: Select all

 &control
    calculation='scf'
    restart_mode='from_scratch',
    prefix='As',
    pseudo_dir = './',
    outdir='./'
 /
 &system
    ibrav=  0, 
    nat= 2,
    ntyp= 1,
    ecutwfc = 30.0
    occupations='smearing',
    smearing='marzari-vanderbilt',
    degauss=0.05
 /
 &electrons
    conv_thr =  1.0d-10
    mixing_beta = 0.7
 /
ATOMIC_SPECIES
 As 74.9216 As-fr.upf

CELL_PARAMETERS (angstrom)
     3.6091988604408050   -0.0000000000356632   -0.0000000000000000
    -1.8045993107470519    3.1256576930176023   -0.0000000000000000
     0.0000000000000000    0.0000000000000000   11.9844661236618624

ATOMIC_POSITIONS (crystal)
As         0.6666669849999991  0.3333330150000009  0.5565192311861956
As         0.3333329559999996  0.6666669250000012  0.4399620088138056

K_POINTS { automatic }
9  9  1 0 0 0

ph.in

Code: Select all

--
&inputph
  prefix   = 'As',
! epsil    = .false.,
  fildyn   = 'As.dyn',
  ldisp    = .true.
  fildvscf = 'dvscf'
  search_sym = .FALSE.
  nq1=3,
  nq2=3,
  nq3=1,
  tr2_ph   =  1.0d-12,
 /

nscf.in

Code: Select all

 &control
    calculation='nscf'
    restart_mode='from_scratch',
    prefix='As',
    pseudo_dir = './',
    outdir='./'
 /
 &system
    ibrav=  0, 
    nat= 2,
    ntyp= 1,
    ecutwfc = 30.0
    occupations='smearing',
    smearing='marzari-vanderbilt',
    degauss=0.05
    nbnd=24
    nosym=.true.
 /
 &electrons
    conv_thr =  1.0d-10
    mixing_beta = 0.7
 /
ATOMIC_SPECIES
 As 74.9216 As-fr.upf
 
CELL_PARAMETERS (angstrom)
     3.6091988604408050   -0.0000000000356632   -0.0000000000000000
    -1.8045993107470519    3.1256576930176023   -0.0000000000000000
     0.0000000000000000    0.0000000000000000   11.9844661236618624

ATOMIC_POSITIONS (crystal)
As         0.6666669849999991  0.3333330150000009  0.5565192311861956
As         0.3333329559999996  0.6666669250000012  0.4399620088138056

K_POINTS crystal
36
0.0 0.0 0.0 0.0277777777778
0.0 0.166666666667 0.0 0.0277777777778
0.0 0.333333333333 0.0 0.0277777777778
0.0 0.5 0.0 0.0277777777778
0.0 0.666666666667 0.0 0.0277777777778
0.0 0.833333333333 0.0 0.0277777777778
0.166666666667 0.0 0.0 0.0277777777778
0.166666666667 0.166666666667 0.0 0.0277777777778
0.166666666667 0.333333333333 0.0 0.0277777777778
0.166666666667 0.5 0.0 0.0277777777778
0.166666666667 0.666666666667 0.0 0.0277777777778
0.166666666667 0.833333333333 0.0 0.0277777777778
0.333333333333 0.0 0.0 0.0277777777778
0.333333333333 0.166666666667 0.0 0.0277777777778
0.333333333333 0.333333333333 0.0 0.0277777777778
0.333333333333 0.5 0.0 0.0277777777778
0.333333333333 0.666666666667 0.0 0.0277777777778
0.333333333333 0.833333333333 0.0 0.0277777777778
0.5 0.0 0.0 0.0277777777778
0.5 0.166666666667 0.0 0.0277777777778
0.5 0.333333333333 0.0 0.0277777777778
0.5 0.5 0.0 0.0277777777778
0.5 0.666666666667 0.0 0.0277777777778
0.5 0.833333333333 0.0 0.0277777777778
0.666666666667 0.0 0.0 0.0277777777778
0.666666666667 0.166666666667 0.0 0.0277777777778
0.666666666667 0.333333333333 0.0 0.0277777777778
0.666666666667 0.5 0.0 0.0277777777778
0.666666666667 0.666666666667 0.0 0.0277777777778
0.666666666667 0.833333333333 0.0 0.0277777777778
0.833333333333 0.0 0.0 0.0277777777778
0.833333333333 0.166666666667 0.0 0.0277777777778
0.833333333333 0.333333333333 0.0 0.0277777777778
0.833333333333 0.5 0.0 0.0277777777778
0.833333333333 0.666666666667 0.0 0.0277777777778
0.833333333333 0.833333333333 0.0 0.0277777777778

epw.in

Code: Select all

--
&inputepw
  prefix      = 'As',
  amass(1)    = 74.9216
  outdir      = './'
  dvscf_dir   = './save'

  elph        = .true.
  kmaps       = .false.
  epbwrite    = .true.
  epbread     = .false.
  epwwrite    = .true.
  epwread     = .false.
  
 efermi_read = .true.           
 fermi_energy = -0.9348 

  wannierize  = .true.
  nbndsub     =  6
  nbndskip    =  2
  num_iter    = 300
  dis_froz_min= -7
  dis_froz_max= 2
  proj(1)     = 'As:p'
  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.333333    0.666667  0'
  wdata(5) = 'K 0.333333    0.666667  0 G 0.00 0.00 0.00'
  wdata(6) = 'end kpoint_path'
  wdata(7) = 'bands_num_points = 100'

  elecselfen  = .false. 
  phonselfen  = .false.
  a2f         = .false.

  fsthick     = 1 ! eV
  eptemp      = 0.075 ! K
  degaussw    = 0.05 ! eV

  nkf1         = 6
  nkf2         = 6
  nkf3         = 1
  nqf1         = 3
  nqf2         = 3
  nqf3         = 1

  nk1         = 6
  nk2         = 6
  nk3         = 1
  nq1         = 3
  nq2         = 3
  nq3         = 1
 /
      3 cartesian
   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00
   0.380327422288020E-11   0.384900205005807E+00   0.000000000000000E+00
   0.333333333339035E+00   0.577350294767568E+00   0.000000000000000E+00

epw2.in

Code: Select all

--
&inputepw
  prefix      = 'As',
  amass(1)    = 74.9216
  outdir      = './'
  dvscf_dir   = './save'

  elph        = .true.
  kmaps       = .true.
  epbwrite    = .false.
  epbread     = .false.
  epwwrite    = .false.
  epwread     = .true.
  
 efermi_read = .true.          
 fermi_energy = -0.9348    
  
  wannierize  = .false.
  nbndsub     =  6
  nbndskip    =  2
  num_iter    = 300
  dis_froz_min= -7
  dis_froz_max= 2
  proj(1)     = 'As:p'
  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.333333    0.666667  0'
  wdata(5) = 'K 0.333333    0.666667  0 G 0.00 0.00 0.00'
  wdata(6) = 'end kpoint_path'
  wdata(7) = 'bands_num_points = 100'

  elecselfen  = .false. 
  phonselfen  = .false.
  a2f         = .false.
  band_plot   = .true.

  fsthick     = 1 ! eV
  eptemp      = 0.075 ! K
  degaussw    = 0.1 ! eV

  filqf       = 'As_band.kpt'
  filkf       = 'As_band.kpt'
  !nkf1         = 6
  !nkf2         = 6
  !nkf3         = 1
  !nqf1         = 3
  !nqf2         = 3
  !nqf3         = 1

  nk1         = 6
  nk2         = 6
  nk3         = 1
  nq1         = 3
  nq2         = 3
  nq3         = 1
 /
      3 cartesian
   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00
   0.380327422288020E-11   0.384900205005807E+00   0.000000000000000E+00
   0.333333333339035E+00   0.577350294767568E+00   0.000000000000000E+00

epw3.in

Code: Select all

--
&inputepw
  prefix      = 'As',
  amass(1)    = 74.9216
  outdir      = './'
  dvscf_dir   = './save'

  elph        = .true.
  kmaps       = .true.
  epbwrite    = .false.
  epbread     = .false.
  epwwrite    = .false.
  epwread     = .true.
  
 efermi_read = .true.       
 fermi_energy = -0.9348    
  
  wannierize  = .false.
  nbndsub     =  6
  nbndskip    =  2
  num_iter    = 300
  dis_froz_min= -7
  dis_froz_max= 2
  proj(1)     = 'As:p'
  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.333333    0.666667  0'
  wdata(5) = 'K 0.333333    0.666667  0 G 0.00 0.00 0.00'
  wdata(6) = 'end kpoint_path'
  wdata(7) = 'bands_num_points = 100'

  elecselfen  = .false. 
  phonselfen  = .true.
  a2f         = .false.
  band_plot   = .false.
  delta_approx= .true.

  fsthick     = 1 ! eV
  eptemp      = 0.075 ! K
  degaussw    = 0.1 ! eV

  filqf       = 'As_band.kpt'
  !filkf       = 'As_band.kpt'
  nkf1         = 20
  nkf2         = 20
  nkf3         = 1
  !nqf1         = 3
  !nqf2         = 3
  !nqf3         = 1

  nk1         = 6
  nk2         = 6
  nk3         = 1
  nq1         = 3
  nq2         = 3
  nq3         = 1
 /
      3 cartesian
   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00
   0.380327422288020E-11   0.384900205005807E+00   0.000000000000000E+00
   0.333333333339035E+00   0.577350294767568E+00   0.000000000000000E+00

epw4.in

Code: Select all

--
&inputepw
  prefix      = 'As',
  amass(1)    = 74.9216
  outdir      = './'
  dvscf_dir   = './save'

  elph        = .true.
  kmaps       = .true.
  epbwrite    = .false.
  epbread     = .false.
  epwwrite    = .false.
  epwread     = .true.
  
 efermi_read = .true.         
 fermi_energy = -0.9348     
  wannierize  = .false.
  nbndsub     =  6
  nbndskip    =  2
  num_iter    = 300
  dis_froz_min= -7
  dis_froz_max= 2
  proj(1)     = 'As:p'
  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.333333    0.666667  0'
  wdata(5) = 'K 0.333333    0.666667  0 G 0.00 0.00 0.00'
  wdata(6) = 'end kpoint_path'
  wdata(7) = 'bands_num_points = 100'

  elecselfen  = .false. 
  phonselfen  = .true.
  a2f         = .true.
  band_plot   = .false.
  delta_approx= .true.

  fsthick     = 1 !eV
  eptemp      = 0.075 ! K
  degaussw    = 0.1 ! eV

  !filqf       = 'As_band.kpt'
  !filkf       = 'As_band.kpt'
  nkf1         = 20
  nkf2         = 20
  nkf3         = 1
  nqf1         = 12
  nqf2         = 12
  nqf3         = 1

  nk1         = 6
  nk2         = 6
  nk3         = 1
  nq1         = 3
  nq2         = 3
  nq3         = 1
 /
      3 cartesian
   0.000000000000000E+00   0.000000000000000E+00   0.000000000000000E+00
   0.380327422288020E-11   0.384900205005807E+00   0.000000000000000E+00
   0.333333333339035E+00   0.577350294767568E+00   0.000000000000000E+00

and the last of epw4.out looks like

Code: Select all

     Atoms inside the unit cell: 

   Cartesian axes

     site n.  atom      mass           positions (a_0 units)

 
     No symmetry!

     G cutoff =    0.0000  (      0 G-vectors)     FFT grid: (  0,  0,  0)
     number of k points=    0
                       cart. coord. in units 2pi/a_0
     EPW          :     0.00s CPU         0.01s WALL

     EPW          :     0.00s CPU         0.01s WALL

     No wavefunction gauge setting applied

     -------------------------------------------------------------------
     Using As.ukk from disk
     -------------------------------------------------------------------

Can anybody tell me how to change the input files to make the calculations go smoothly.

Looking forward for the replies.

Best wishes
Yue Gu

Posted: **Sun Jul 26, 2020 4:58 pm**

I got the same error, Any solution so far

Posted: **Sun Jul 26, 2020 5:58 pm**

Dear Sylvester10:

Could you be more specific to the error you encountered?
Did you encounter the error when doing calculations for 2D-As as in YueGu's post?

I can't give you the clear answer if you don't provide the information on the calculation in which you encountered the error.

I need the followings:
(1) the version of EPW you are using
(2) all inputs and outputs with the information on pseudopotentials.

Sincerely,

H. Lee

Posted: **Mon Jul 27, 2020 11:30 am**

Hi Sylvester10 and Yue Gu,

Are you sure, you have 'As.ukk' (rotation matrix) saved from the 1st calculation (epw1.in)? I wonder if you have removed it by any chance.

Best,
Hari Paudyal

EPW Forum

stop at "Using .ukk from disk"

stop at "Using .ukk from disk"

Re: stop at "Using .ukk from disk"

Re: stop at "Using .ukk from disk"

Re: stop at "Using .ukk from disk"