[Pw_forum] About GWW

[shudong wang]

Dear developers:
I use GWW to calculate my 2D unit cell with 2C and 2F atoms,first I have
performed  scf_head.in with pw.x code ,it works well.Then I have a ph
calculation at Gamma with head.x but it appears:

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         0
     from solve_e : error #         1
     called in the wrong case
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
My input file is:

&control
    title='gf'
    calculation='scf',
    restart_mode='from_scratch',
    prefix='gf'
    pseudo_dir = '/home/sdwang/espresso-4.1.2/pseudo',
    outdir='/home/sdwang/gw/tmp'
/
&system
    ibrav = 4, celldm(1)=4.839,celldm(3)=7.023,nat=4, ntyp=2,
    ecutwfc = 60.0,nosym=.false.,
occupations='smearing',degauss=0.01,smearing='gaussian'
/
&electrons
    diagonalization='david',
    conv_thr =  1.0d-10,
    mixing_beta = 0.5,
    startingwfc='random',
/
ATOMIC_SPECIES
 C   12.0  C.pz-vbc.UPF
 F   19.0  F.pz-rrkj.UPF
ATOMIC_POSITIONS (crystal)
C        0.331724716   0.667136313   0.096514182
C        0.665075100   0.333783657   0.069110636
F        0.331574264   0.667153962   0.171629429
F        0.665085920   0.333706069  -0.005994247
K_POINTS (automatic)
9 9 1  0 0 0


phonons of gf at Gamma
 &inputph
  tr2_ph=1.d-4,
  prefix='gf',
  outdir='/home/sdwang/gw/tmp/',
  epsil=.false.,
  amass(1)=12.0,
  amass(2)=19.0
  fildyn='gf.dynG',
  recover=.false.
  trans=.false.
  l_head=.true.
  omega_gauss=20.
  n_gauss=79
  grid_type=2
  niter_ph=1
 /
0.0 0.0 0.0
What is my problem?
Thanks!!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20100705/750612d9/attachment.html>