[Pw_forum] pw bands calculation crash

Davide Tiana d.tiana at bath.ac.uk
Fri Apr 27 16:04:46 CEST 2012


Dear all,
After optimised a molecule I tried to calculate the band structure but  
pw crashes and I can't figure out why.

The optimisation step should be fine(here it converge immediately  
since I rerun my previous calculation):

     BFGS Geometry Optimization

      bfgs converged in   1 scf cycles and   0 bfgs steps
      (criteria: energy < 0.10E-03, force < 0.10E-02)

      Final energy   =    -646.2249695321 Ry
Begin final coordinates

ATOMIC_POSITIONS (crystal)
C        0.953315090   0.000413805   0.235484404
..
..
H        1.009939685   0.169442779   0.056058275
H        1.044496546   0.158711398  -0.056214459
End final coordinates



      Writing output data file mol_a.save

..
..
..
    This run was terminated on:  13: 4:41  27Apr2012

=------------------------------------------------------------------------------=
    JOB DONE.
=------------------------------------------------------------------------------=
"

but when I run pw.x to calculate bands it crashes

"
      Atomic positions and unit cell read from directory:
      ./mol_a.save/


      Subspace diagonalization in iterative solution of the eigenvalue problem:
      a serial algorithm will be used

Application 1981542 exit codes: 139
Application 1981542 exit signals: Killed
Application 1981542 resources: utime ~2s, stime ~2s
"
My input should not contain any errors. I used it for other systems  
similar to this one and it works. just in case:
"
&control
     calculation='bands',
     title='mol_a',
     wf_collect=.true.,
     outdir='./',
     prefix='mol_a'
     etot_conv_thr=1.0d-4,
     forc_conv_thr=1.0d-3,
  /
  &system
    ibrav=14,
    celldm(1)=28.4693677,
    celldm(2)=1.058648879,
    celldm(3)=0.600431923,
    celldm(4)=0,
    celldm(5)=0,
    celldm(6)=0,
    nat=23,
    ntyp=6,
    nbnd=70
    tot_charge=0.0,
    ecutwfc=45.0,
    ecutrho=450.0,
    nosym=.false.,
    nosym_evc=.false.,
    london=.true.,
    london_s6=0.75,
    london_rcut=200,
  /
  &electrons
    electron_maxstep=80,
    conv_thr=1.0d-6,
    mixing_mode='plain',
    mixing_beta=0.7,
    mixing_ndim=8,
    mixing_fixed_ns=0,
    diagonalization='david',
  /
  ATOMIC_SPECIES
    H  1.007   H.pw91-van_ak.UPF
  ...
    Zn 65.409  Zn.pw91-n-van.UPF
  ATOMIC_POSITIONS {crystal}
    C   0.952835      0.500000      0.826917
....
Zn   0.0000        0.500000      0.191590
K_POINTS {crystal}
11
0 0  0 1
0 0 .1 2
0 0 .2 3
0 0 .3 4
0 0 .4 5
0 0 .5 6
0 0 .6 7
0 0 .7 8
0 0 .8 9
0 0 .9 10
0 0  1 11
"


Any suggestion?
thanks a lot
Davide




More information about the users mailing list