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 <pw_forum@pwscf.org><br>
<br>
> <br>
> Dear Pieremanuele,<br>
> <br>
> if you don't have the crystal coordinates of the atoms you can <br>
> have the <br>
> code print them.<br>
> <br>
> <br>
> On May 27 2009, Pieremanuele Canepa wrote:<br>
> <br>
> >Dear all,<br>
> >I was  able to compile  the program resp_mat.f90 in <br>
> order to  carry out the<br>
> >U in a self consistent way. Now I am trying to write up <br>
> correctly the input<br>
> >for resp_mat.x. I have realized reading trough the guide (found <br>
> at link<br>
> >http://vlab.msi.umn.edu/events/lecture.shtml)  that I do <br>
> need to have an<br>
> >external file   with the fractional coordinate of the <br>
> atom belonging to the<br>
> >supercell.<br>
> >Actually, if I look at those external file (call for example <br>
> pos....) that<br>
> >come togheter with the guide found online, the second part is <br>
> concerning to<br>
> >the fractional coordinate of the atoms in the supercell, while <br>
> what is the<br>
> > first part? Are these the vector of the supercell in Cartesian <br>
> > coordinates ?<br>
> <br>
> yes.<br>
> <br>
> >Why, when  for instance  the   Cartesian <br>
> coordinate of the  supercell<br>
> >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>
> >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>
> <br>
> it doesn't matter the length of the vectors. only the ratio <br>
> between length <br>
> does.<br>
> This is used to compute the distance between atoms ut the scale <br>
> doesn't <br>
> matter at all.<br>
> <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 = <br>
> 'neutral' mean  in the<br>
> >input file for resp_mat.x?<br>
> ><br>
> <br>
> neutral means that you are enforcing the total response of the <br>
> system to be <br>
> neutral. so there will be a compensating "background" term in <br>
> the response <br>
> matrix that results in an additional column and row. see PRB 71 35105.<br>
> <br>
> hope this helps. <br>
> <br>
> regards,<br>
> <br>
> Matteo<br>
> <br>
> <br>
> I hope  you can get my point!<br>
> >Best Regards, Piero<br>
> ><br>
> _______________________________________________<br>
> Pw_forum mailing list<br>
> Pw_forum@pwscf.org<br>
> http://www.democritos.it/mailman/listinfo/pw_forum<br>
> <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>