Dear all, <br>I was  able to compile  the program resp_mat.f90 in order to  carry out the U in a self consistent way. Now I am trying to write up correctly the input for resp_mat.x. I have realized reading trough the guide (found at link <a href="http://vlab.msi.umn.edu/events/lecture.shtml" target="_blank">http://vlab.msi.umn.edu/events/lecture.shtml</a>)  that I do need to have an external file   with the fractional coordinate of the atom belonging to the supercell. <br>
Actually, if I look at those external file (call for example pos....) that come togheter with the guide found online, the second part is concerning to the fractional coordinate of the atoms in the supercell, while what is the first part? Are these the vector of the supercell in Cartesian coordinates ? Why, when  for instance  the   Cartesian coordinate of the  supercell vectors printed in the early part of the PWscf output are:<br>
  crystal axes: (cart. coord. in units of a_0)<br>               a(1) = (  1.000000  0.000000  0.000000 )  <br>               a(2) = (  0.000000  1.000000  0.000000 )  <br>               a(3) = (  0.000000  0.000000  1.000000 )  <br clear="all">
in the file pos  these become<br>20.d0 0.d0 0.d0<br>0.d0 20.d0 0.d0<br>0.d0 0.d0 20.d0<br>?<br>The cell parameter  in this case is 10.8400 (FeNi in the guide)<br><br>Another question is : what does the line   back = 'neutral' mean  in the input file for resp_mat.x? <br>
I hope  you can get my point!<br>Best Regards, Piero<br>-- <br>Pieremanuele Canepa<br>Room 230<br>School of Physical Sciences, Ingram Building,<br>University of Kent, Canterbury, Kent,<br>CT2 7NH<br>United Kingdom<br>-----------------------------------------------------------<br>
<br>