<div>Dear Paolo U,</div>  <div> </div>  <div>Thanks for your reply. I'm very interested for your finite electric field method.<BR>The dielectric polarizability tensor can be calculated by taking dinite differences of  atomic forces.<BR>I have several questions to ask you.</div>  <div> </div>  <div>1. We obtained the diagonal terms, by considering values for the electric field of<BR>-2h,-h, 0, h and 2h in the cp.x code [tefield=.true. <BR>                                                        tcg =
 .true.,<BR>                                                      epol = 1 or 2 or 3<BR>                                                     efield = -2h,-h, 0, h or 2h],<BR>and by using the five-point formular . Is it right?</div>  <div> </div>  <div>2. For the calculation of the off-diagonal terms, two different electric fields <BR>are required. This need apply electric fields of intensity -2h/1.414,-h/1.414, <BR>0, h/1.414 and 2h/1.414 along the three different
 couples of cartesian directions, and apply eauation.<BR>But how can I apply simultaneously finite fields along two different Cartesian directions<BR>in cp.x code?[ESSPRESSO-3.2]<BR> <BR>3. In this way, we could obtain all the dielectric polarizability tensor. Then the<BR>Raman susceptibilities associated to the normal mode n can be obtianed according to the equation. <BR>How to implement it ? By a small code?</div>  <div> </div>  <div>Best!</div>  <div> </div>  <div>I'm counting on the new version.</div>  <div> </div>  <div>Niu Li</div><p>



      <hr size=1><a href="http://cn.mail.yahoo.com" target=blank>ÇÀ×¢ÑÅ»¢Ãâ·ÑÓÊÏä3.5GÈÝÁ¿£¬20M¸½¼þ£¡</a>