<div dir="ltr">Dear Forzad,<div><br></div><div>if you search the QE webpage (in the "resources" section, I believe) you will find the notes and lectures from several past QE schools, workshops, tutorials and there should be some about the calculation of U. One is in the Santa Barbara school from 2009. A second one should be among the notes from a school in Pune last year.</div><div><br></div><div>hope this will help you.</div><div><br></div><div>best regards,</div><div><br></div><div>Matteo</div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 8, 2015 at 12:11 PM, Farzad Molani <span dir="ltr"><<a href="mailto:farzad_c81@yahoo.com" target="_blank">farzad_c81@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:HelveticaNeue-Light,Helvetica Neue Light,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"><div><span>Dear QE users,</span></div><div dir="ltr"><span>I want to calculate Hubbard parameter. Has any body  a tutorial "step by step" for calculating the hubbard parameter? Is this input file correct? How can I use the output file?</span></div><div dir="ltr"><br><span></span></div><div dir="ltr"><span><br></span></div><div dir="ltr"><span>&CONTROL<br>  calculation  ="scf",<br>  pseudo_dir   = "./",<br>  restart_mode = 'restart',<br>  outdir       = "./",<br>  prefix       = "geraphene-Fe-ldaU",<br>  verbosity    = 'high',<br>  tprnfor       = .true.,<br>  tstress       = .true.,<br>/<br>&SYSTEM<br>  ibrav       = 4,<br>  celldm(1)   = 9.282,<br>  celldm(3)   = 2.035,<br>  nat         = 9,<br>  ntyp        = 2,<br>  ecutwfc     = 44.2,<br>  ecutrho     =300,<br>  nspin       = 2,<br>  starting_magnetization(1)=0.5,<br>  starting_magnetization(2)=0.5,<br>  occupations='smearing',<br>  smearing   ='m-v',<br>  degauss    =0.003,<br>  lda_plus_u = .true.<br>  Hubbard_alpha(2)=-0.04,<br>/<br>&ELECTRONS<br>  conv_thr  = 1.0d-10,<br>  mixing_beta = 0.5d0,<br>  diagonalization  = 'cg',<br>/<br>&IONS<br>/<br>&CELL<br>/<br>ATOMIC_SPECIES<br>C   12.010   C.pz-rrkjus.UPF<br>Fe  55.8452  Fe.pz-nd-rrkjus.UPF<br>ATOMIC_POSITIONS {angstrom}<br>C       -0.034035640   1.396956433   4.962804121<br>C        1.208811211   0.679456444   4.961853267<br>C        2.451565749   1.396929605   4.962300355<br>C        3.664662075   0.696596615   5.017209108<br>C       -1.247079599   3.532441567   5.017139042<br>C       -0.034018302   2.832082543   4.961775758<br>C        1.208707850   3.549535163   4.962481191<br>C        2.451470952   2.831954512   4.963586323<br>Fe       1.207584375   2.113742807   6.427714456<br>K_POINTS automatic<br> 5 5 1 0 0 0</span></div><div><div> </div><div><br></div><div>Best regards<span class="HOEnZb"><font color="#888888"><br></font></span></div></div><span class="HOEnZb"><font color="#888888"><div>Farzad Molani, <br>Ph.D in Computational Chemistry<br>Department of Chemistry, <br>Islamic Azad University of Sanadaj<br>Tel.: 009891 4442 3308</div></font></span></div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div></div>