Dear  Ivanov,<br><br>if you have a look to the file resp_mat.f90 you can check the reason to have atomic positions read from file.<br>there are needed to calculate the matrix elements of chi and chi0.<br><br>atomic positions are given in fractional coordinates (direct lattice).<br>
what about primitive vectors?<br><br>if you are considering a FeO unit cell with 8 ions you can perturb only one Fe site (you don`t need to look at occupancies on O sites) if  all the other Fe are equivalent to this one.<br>
one have to perturb once all the ions which are not equivalent.<br><br>cheers,<br><br>Gianluca<br><br><div class="gmail_quote">2010/2/16 Максим Иванов <span dir="ltr"><<a href="mailto:maaaxim@yahoo.com">maaaxim@yahoo.com</a>></span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit;" valign="top">
Dear PW-users,<br>I have question about calculating the response matrice in lda+u calculation. At the tutorial with examples at page<a href="http://media.quantum-espresso.org/santa_barbara_2009_07/index.php" target="_blank"> http://media.quantum-espresso.org/santa_barbara_2009_07/index.php</a> is used input file to calculate the matrice. In this input file you need file containing the atomic positions and the vectors of the unit cell plus number of unit cells in each direction to construct the supercell. As I understand response matrice is resiprocal matrix of matrix containing dn/dalfa derivatives. So, why does author include atomic positions in his program? <br>
In the case when there are more than two Fe atoms in FeO unit cell (let's say 8) is it correct to make perturbation only for one Fe atom? Or should I make perturbation for 8 Fe atoms and calculate 8x8
 matrix?<br>Thanks in advance,<br>Maxim Ivanov<br><div class="hm"><table border="0" cellpadding="0" cellspacing="0" width="100%"><tbody><tr><td><br></td><td><br></td></tr></tbody></table><br><br></div></td></tr></tbody></table>
<div class="hm"><br>





      <hr size="1"><br>Вы уже с Yahoo!? 
Испытайте обновленную и улучшенную. <a href="http://ru.mail.yahoo.com" target="_blank">Yahoo! Почту</a>!</div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br>