<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 [like example 30</div>  <div>                                                     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 mixedterms, two different electric fields are required. This need apply electric fields of intensity  - 2h/1.414 , </div>  <div>-h/1.414, 0, h/1.414 and 2h/1.414 along the three different couples of cartesian directions, and apply equations.<BR>But how can I apply simultaneously finite fields along two different Cartesian directions 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 equations. <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>