Indeed! Thanks very much!<br>BTW, is there a "tr1_ph"? This name convention is really interesting:)<br><br>Best,<br>Jin Zhang<br><br><div><span class="gmail_quote">On 11/27/06, <b class="gmail_sendername">Xunlei Ding

</b> <<a href="mailto:ding@sissa.it">ding@sissa.it</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear Zhang,<br><br>

The tr2_ph is too large in you case.<br>Just change the tr2_ph=1.0d-14, you will get:<br>    omega( 1) =       2.198936 [THz] =      73.349094 [cm-1]<br>     omega( 2) =       3.192379 [THz] =     106.487008 [cm-1]<br>     omega( 3) =       

3.192379 [THz] =     106.487008 [cm-1]<br>     omega( 4) =       4.915826 [THz] =     163.975388 [cm-1]<br>     omega( 5) =       4.915826 [THz] =     163.975388 [cm-1]<br>     omega( 6) =     131.771213 [THz] =    4395.443762

 [cm-1]<br><br>Yours,<br>Ding<br><br>Jin Zhang wrote:<br><br>> Hi all,<br>><br>>    I posted this issue one week ago. Now I have relaxed the structure<br>> 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<br>> and three positive ones( values are 898,898,3679cm^-1 ). But there<br>> 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). > > 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<br>> 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>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum">

http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br>Dep. of Physics, Peking University<br>Beijing, P.R.China<br>Phone: 86-10-6275-3190