[Pw_forum] error of calculation for electeon-phonon using the saved fildvscf and fildyn

oyxf328 oyxf328 at 126.com
Thu Jun 2 08:54:25 CEST 2016


Dear  developer:
      I am calculating the electron_phonon coupling and the lambda of superconduct case(FeSe) using the version 5.1. 
     First step, I have done the phonon calculation to obtain the fildvscf, fildrho and fildyn。 These files saved in the directory. The second step, I need to calculate the electron_phonon using the saved fildvscf and fildyn. But I always was prevented by the  following errors:


     forrtl: severe (59): list-directed I/O syntax error, unit 26, file /data2/ouyang/work/fese-qe/NM+CAFM/3/fese.dyn3
Image              PC                Routine            Line        Source
ph.x               0000000000B6985D  Unknown               Unknown  Unknown
ph.x               0000000000B68365  Unknown               Unknown  Unknown
ph.x               0000000000B053D9  Unknown               Unknown  Unknown
ph.x               0000000000A8E45A  Unknown               Unknown  Unknown
ph.x               0000000000A8DC50  Unknown               Unknown  Unknown
ph.x               0000000000ABFE3B  Unknown               Unknown  Unknown
ph.x               0000000000ABD9C8  Unknown               Unknown  Unknown
ph.x               0000000000459FE0  readmat_                  192  elphon.f90
ph.x               000000000045659A  elphon_                   110  elphon.f90
ph.x               0000000000447166  do_phonon_                 95  do_phonon.f90
ph.x               000000000044429C  MAIN__                     75  phonon.f90
ph.x               00000000004441DC  Unknown               Unknown  Unknown
libc.so.6          000000328BE1D994  Unknown               Unknown  Unknown
ph.x               00000000004440E9  Unknown               Unknown  Unknown
I will list the input files as following:
the first input file for scf:
&control
    calculation        = 'scf'
    restart_mode       = 'from_scratch'
    prefix             = 'fese'
    pseudo_dir         = '/data2/ouyang/software/espresso-5.1/espresso-5.1/pseudo/'
    outdir             = './tempdir/'
/
&system
    ibrav              = 0
    nat                = 4
    ntyp               = 4
    ecutwfc            = 60
    ecutrho            = 600
    occupations        = 'smearing'
    smearing           = 'methfessel-paxton'
    degauss            =  0.02
    nspin              = 2
    starting_magnetization(1) =1
    starting_magnetization(2) =-1
    starting_magnetization(3) =0
    starting_magnetization(4) =0
/
&electrons
    electron_maxstep   = 50
    conv_thr           = 1.0D-14
    mixing_beta        = 0.7
    mixing_mode        = 'plain'
/
ATOMIC_SPECIES
   Fe1    55.845        Fe.pbe-nd-rrkjus.UPF
   Fe2    55.845        Fe.pbe-nd-rrkjus.UPF
   Se1    78.960        Se.pbe-n-rrkjus_psl.0.2.UPF
   Se2    78.960        Se.pbe-n-rrkjus_psl.0.2.UPF
CELL_PARAMETERS {angstrom}
3.754198    0.000000   0.000000
0.000000    3.754198   0.000000
0.000000    0.000000   6.545954
ATOMIC_POSITIONS {crystal}
Fe1   0.750000      0.250000      0.000000
Fe2   0.250000      0.750000      0.000000
Se1   0.750000      0.750000      0.774837
Se2   0.250000      0.250000      0.225163
K_POINTS (automatic)
12 12 6 0 0 0




The second file for phonon calculation:
Electron-phonon coefficients for fese
 &inputph
  tr2_ph=1.0e-14
  alpha_mix=0.3
  prefix='fese'
  fildvscf='fesedv'
  fildrho='drho'
  amass(1)=55.845
  amass(2)=55.845
  amass(3)=78.96
  amass(4)=78.96
  outdir='tempdir'
  fildyn='fese.dyn'
  start_q   = 3
  last_q    = 3
  trans=.true.
  ldisp=.true.
  nq1=2, nq2=2, nq3=2




The third file for calculating the electron_phonon using the saved fildyn and fildvscf :
Electron-phonon coefficients for fese
 &inputph
  tr2_ph=1.0e-14
  alpha_mix=0.3
  prefix='fese'
  fildvscf='fesedv'
  fildrho='drho'
  amass(1)=55.845
  amass(2)=55.845
  amass(3)=78.96
  amass(4)=78.96
  electron_phonon='interpolated'
  el_ph_sigma=0.005
  el_ph_nsigma=20
  outdir='tempdir'
  fildyn='fese.dyn'
  start_q   = 3
  last_q    = 3
  trans=.false.
  ldisp=.true.
  nq1=2, nq2=2, nq3=2




I want to know how to solve the problem. I will be very appreciative for your answer..





 
 
 
 
 
 






 





 





 





 





 





 





 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20160602/65a2f21d/attachment.html>


More information about the users mailing list