<DIV>dear all,</DIV>  <DIV>   I calcultated the Raman tensor  using the pwscf version 3.0. The scf part run correctly. But the the running of the ph.x will stop unexpectedly. The fowllowing part is the ph.out.</DIV>  <DIV> </DIV>  <DIV> Program PHONON    v.3.0    starts ...<BR>     Today is 26Mar2006 at  0:46:20 </DIV>  <DIV>     Ultrasoft (Vanderbilt) Pseudopotentials</DIV>  <DIV>     Reading file SnO.save ...      only dimensions<BR>     read complete</DIV>  <DIV>     Reading file SnO.save ...      all except wavefuctions<BR>     read complete</DIV>  <DIV>     nbndx  =    10  nbnd   =    10  natomwfc =    64  npwx  
 =    9600<BR>     nelec  =  20.00  nkb   =    52  ngl    =    4494<BR>     WRITING PATTERNS TO FILE drho.pat</DIV>  <DIV>     phonons of SnO at Gamma                                                    </DIV>  <DIV>                  ........................</DIV>  <DIV>                  .........................</DIV> 
 <DIV>                  ........................</DIV>  <DIV>     Dielectric constant from finite-differences</DIV>  <DIV>          (       9.668658880       0.000000000       0.000000000 )<BR>          (       0.000000000       9.668658880       0.000000000 )<BR>          (       0.000000000       0.000000000       6.794783712 )</DIV>  <DIV>              Electro-optic tensor is defined as
 <BR>            the derivative of the dielectric tensor <BR>              with respect to one electric field <BR>                 units are Rydberg a.u. </DIV>  <DIV>            to obtain the static chi^2 multiply by 1/2  </DIV>  <DIV>            to convert to pm/Volt multiply per 2.7502  </DIV>  <DIV><BR>          Electro-optic tensor in cartesian axis: </DIV>  <DIV>          (       0.000000000       0.000000000       0.000000000 )<BR>         
 (       0.000000000       0.000000000       0.000000000 )<BR>          (       0.000000000       0.000000000       0.000000000 )</DIV>  <DIV>          (       0.000000000       0.000000000       0.000000000 )<BR>          (       0.000000000       0.000000000       0.000000000 )<BR>          (       0.000000000       0.000000000       0.000000000 )</DIV> 
 <DIV>          (       0.000000000       0.000000000       0.000000000 )<BR>          (       0.000000000       0.000000000       0.000000000 )<BR>          (       0.000000000       0.000000000       0.000000000 )</DIV>  <DIV><BR>     Computing Second order response </DIV>  <DIV><BR>      iter #   1   av.it.:  13.9<BR>      thresh= 0.100E-01 alpha_mix =  0.700 |ddv_scf|^2 =  0.476E-02 "</DIV>  <DIV> </DIV>  <DIV>It stopped here without any error warning.</DIV>  <DIV>So I
 the change parameters such as tr2_ph, kpoints and ecut to the lower limitition. It run correctly. But  "Electro-optic tensor in cartesian axis" is the same as above. So please help me with this.</DIV>  <DIV> <BR>ps: the  phonon frequencies at Gama point  can be calculated freely.</DIV><p>
                <hr size=1><a href="http://cn.mail.yahoo.com" target=blank> 
ÑÅ»¢1GÃâ·ÑÓÊÏä°Ù·Ö°Ù·ÀÀ¬»øÐÅ</a><br> 
<a href="http://cn.zs.yahoo.com" target=blank>ÑÅ»¢ÖúÊÖ-ËÑË÷¡¢É±¶¾¡¢·ÀɧÈÅ</a>