Thanks a lot Matteo! <br>
<br>
----- Original Message -----<br>
From: matteo@umn.edu<br>
Date: Thursday, May 28, 2009 10:24 am<br>
Subject: Re: [Pw_forum] input for resp_mat.x in to optimize the U for LDA+U calculations<br>
To: PWSCF Forum &lt;pw_forum@pwscf.org&gt;<br>
<br>
&gt; <br>
&gt; Dear Pieremanuele,<br>
&gt; <br>
&gt; if you don't have the crystal coordinates of the atoms you can <br>
&gt; have the <br>
&gt; code print them.<br>
&gt; <br>
&gt; <br>
&gt; On May 27 2009, Pieremanuele Canepa wrote:<br>
&gt; <br>
&gt; &gt;Dear all,<br>
&gt; &gt;I was&nbsp; able to compile&nbsp; the program resp_mat.f90 in <br>
&gt; order to&nbsp; carry out the<br>
&gt; &gt;U in a self consistent way. Now I am trying to write up <br>
&gt; correctly the input<br>
&gt; &gt;for resp_mat.x. I have realized reading trough the guide (found <br>
&gt; at link<br>
&gt; &gt;http://vlab.msi.umn.edu/events/lecture.shtml)&nbsp; that I do <br>
&gt; need to have an<br>
&gt; &gt;external file&nbsp;&nbsp; with the fractional coordinate of the <br>
&gt; atom belonging to the<br>
&gt; &gt;supercell.<br>
&gt; &gt;Actually, if I look at those external file (call for example <br>
&gt; pos....) that<br>
&gt; &gt;come togheter with the guide found online, the second part is <br>
&gt; concerning to<br>
&gt; &gt;the fractional coordinate of the atoms in the supercell, while <br>
&gt; what is the<br>
&gt; &gt; first part? Are these the vector of the supercell in Cartesian <br>
&gt; &gt; coordinates ?<br>
&gt; <br>
&gt; yes.<br>
&gt; <br>
&gt; &gt;Why, when&nbsp; for instance&nbsp; the&nbsp;&nbsp; Cartesian <br>
&gt; coordinate of the&nbsp; supercell<br>
&gt; &gt;vectors printed in the early part of the PWscf output are:<br>
&gt; &gt;&nbsp; crystal axes: (cart. coord. in units of a_0)<br>
&gt; &gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(1) = (&nbsp; 1.000000&nbsp; 0.000000&nbsp; 0.000000 )<br>
&gt; &gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(2) = (&nbsp; 0.000000&nbsp; 1.000000&nbsp; 0.000000 )<br>
&gt; &gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(3) = (&nbsp; 0.000000&nbsp; 0.000000&nbsp; 1.000000 )<br>
&gt; &gt;in the file pos&nbsp; these become<br>
&gt; &gt;20.d0 0.d0 0.d0<br>
&gt; &gt;0.d0 20.d0 0.d0<br>
&gt; &gt;0.d0 0.d0 20.d0<br>
&gt; &gt;?<br>
&gt; <br>
&gt; it doesn't matter the length of the vectors. only the ratio <br>
&gt; between length <br>
&gt; does.<br>
&gt; This is used to compute the distance between atoms ut the scale <br>
&gt; doesn't <br>
&gt; matter at all.<br>
&gt; <br>
&gt; <br>
&gt; &gt;The cell parameter&nbsp; in this case is 10.8400 (FeNi in the guide)<br>
&gt; &gt;<br>
&gt; &gt;Another question is : what does the line&nbsp;&nbsp; back = <br>
&gt; 'neutral' mean&nbsp; in the<br>
&gt; &gt;input file for resp_mat.x?<br>
&gt; &gt;<br>
&gt; <br>
&gt; neutral means that you are enforcing the total response of the <br>
&gt; system to be <br>
&gt; neutral. so there will be a compensating "background" term in <br>
&gt; the response <br>
&gt; matrix that results in an additional column and row. see PRB 71 35105.<br>
&gt; <br>
&gt; hope this helps. <br>
&gt; <br>
&gt; regards,<br>
&gt; <br>
&gt; Matteo<br>
&gt; <br>
&gt; <br>
&gt; I hope&nbsp; you can get my point!<br>
&gt; &gt;Best Regards, Piero<br>
&gt; &gt;<br>
&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; Pw_forum@pwscf.org<br>
&gt; http://www.democritos.it/mailman/listinfo/pw_forum<br>
&gt; <br>
<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>