Hi all,<br><br>   I posted this issue one week ago. Now I have relaxed the structure and got total force less than 5.0D-5. I used a large cubic supercell with 30bohr alat and ecutwfc up to 50Ry. However, things only went worse. The highest frequency from 
ph.x dropped from 3187 cm^(-1) to 1849, while gaussian03 shows it should be around 3300cm^(-1).<br><br>I also checked H2 molecule. In this case, I got three negative values and three positive ones( values are 898,898,3679cm^-1 ). But there should be only ONE vibrational frequency which corresponds to the stretching mode along bonding direction. And the highest freq. is also too low (should be around 43xx).
<br><br>I really couldn't figure out what's wrong. Can someone help me? Thanks in advance!<br><br>The input files for H2 are: <br>---------------------------scf_input-----------------(run two times just like example07 )-----
<br> &control<br>    calculation='scf'<br>    restart_mode='from_scratch',<br>    prefix='h2',<br>    pseudo_dir = '$PSEUDO_DIR/',<br>    outdir='$TMP_DIR/'<br> /<br> &system<br>    ibrav=  8, celldm(1) =18, celldm(2)=
1.0, celldm(3)=1.0,<br>    nat= 2, ntyp= 1,<br>    ecutwfc = 50.0,   <br>    occupations='smearing', smearing='gausssian', degauss=0.01,<br>    la2F = .false., nspin=1,<br> /<br> &electrons<br>    conv_thr =  1.0d-8<br>
    mixing_beta = 0.7<br> /<br>ATOMIC_SPECIES<br>  H     1.008    H.blyp-vbc.UPF<br>ATOMIC_POSITIONS (bohr)<br>H        0.000000000   0.000000000  -0.707244771<br>H        0.000000000   0.000000000   0.707244771<br>K_POINTS {automatic}
<br>1 1 1 0 0 0<br>------------------------------ph_input---------------------------<br> &inputph<br>  tr2_ph=1.0d-8,<br>  prefix='h2',<br>  fildvscf='h2dv',<br>  amass(1)=1.008,<br>  outdir='$TMP_DIR/',<br>  fildyn='
h2.dyn',<br>  elph=.false.,<br>  trans=.true.,<br>  ldisp=.true.,<br>  nq1=1, nq2=1, nq3=1<br>--------------------------------------------------------------------<br><br>
Regards,<br>
Jin Zhang<br>
<br>