stop at "Using .ukk from disk"

General discussion around the EPW software

Moderator: stiwari

Post Reply
YueGu
Posts: 1
Joined: Tue Jun 04, 2019 2:13 am
Affiliation:

stop at "Using .ukk from disk"

Post by YueGu »

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
Top
Sylvester10
Posts: 15
Joined: Sat Jul 18, 2020 9:35 pm
Affiliation: Kenyatta University

Re: stop at "Using .ukk from disk"

Post by Sylvester10 »

I got the same error, Any solution so far
Top
hlee
Posts: 415
Joined: Thu Aug 03, 2017 12:24 pm
Affiliation: The University of Texas at Austin

Re: stop at "Using .ukk from disk"

Post by hlee »

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
Top
hpaudya1
Posts: 190
Joined: Tue Mar 21, 2017 7:11 pm
Affiliation:

Re: stop at "Using .ukk from disk"

Post by hpaudya1 »

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
Top
Post Reply

Return to “General discussion”