It is currently Thu Dec 12, 2019 12:14 am

Tell a friend!


Post new topic Reply to topic  [ 1 post ] 
stop at "Using .ukk from disk" 
Author Message

Joined: Tue Jun 04, 2019 2:13 am
Posts: 1
University: tsinghua
Post stop at "Using .ukk from disk"
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:
&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:
--
&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:
&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:
--
&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:
--
&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:
--
&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:
--
&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:
     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


Mon Jun 17, 2019 1:25 am
Profile E-mail
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum

Search for:
Jump to:  
cron
Powered by phpBB © phpBB Group.
Designed by Vjacheslav Trushkin

(All content on this board is governed by and is the sole responsibility of the board administrator.)


Gratis forum Free forum hosting| gratis phpbb3 forum | phpbb3 styles