<br clear="all"><div><div style="text-align: left;"><span style="font-family:'Comic Sans MS';font-size:12px;line-height:14px">Dear </span>Hongsheng</div><div style="text-align: left;"><span class="Apple-style-span" style="font-family: 'Comic Sans MS'; font-size: 12px; line-height: 14px; "><br>
</span></div><div style="text-align: left;"><span class="Apple-style-span" style="font-family: 'Comic Sans MS'; font-size: 12px; line-height: 14px; ">In this case I think that the more important issue to be noticed is considering spin polarization for most of atoms. Many mistakes in this issue arise from not considering spin polarization.</span></div>
<div style="text-align: left;"><span class="Apple-style-span" style="font-family: 'Comic Sans MS'; font-size: 12px; line-height: 14px; ">I guess </span><font class="Apple-style-span" face="'Comic Sans MS'"><span class="Apple-style-span" style="font-size: 12px; line-height: 14px;">Considering "Symmetry" may symmetrize the charge density and causes some changes in electronic structure of atom (population of orbitals) or may do reversely with nosym='true'. </span></font></div>
<div style="text-align: left;"><font class="Apple-style-span" face="'Comic Sans MS'"><span class="Apple-style-span" style="font-size: 12px; line-height: 14px;">Also, it is impossible for me </span></font><span class="Apple-style-span" style="font-family: 'Comic Sans MS'; font-size: 12px; line-height: 14px; ">to describe this issue based on DFT and the exchange-correlation parts in details.</span></div>
<div style="text-align: left;"><span class="Apple-style-span" style="font-family: 'Comic Sans MS'; font-size: 12px; line-height: 14px; "><br></span></div><div style="text-align: left;"><font class="Apple-style-span" face="'Comic Sans MS'"><span class="Apple-style-span" style="font-size: 12px; line-height: 14px;">In quantum chemistry software based on Gaussian type orbitals basis set (not plane wave basis set) we compare energy of an atom with considering different spin multiplicities then find the lowest energy and call it "ground state". I can do it for some atoms (whatever you select) by Firefly software (former PC-GAMESS) and then you can compare it with the corresponding calculation by QE. </span></font></div>
<div style="text-align: left;"><font class="Apple-style-span" face="'Comic Sans MS'"><span class="Apple-style-span" style="font-size: 12px; line-height: 14px;">These days I am very busy and going on a journey so let me do it in a due course; contact me :) </span></font></div>
<div style="text-align: left;"><span class="Apple-style-span" style="font-family: 'Comic Sans MS'; font-size: 12px; line-height: 14px; "><br></span></div><div style="text-align:center"><span style="font-family:'Comic Sans MS';font-size:12px;line-height:14px">Best Wishes</span></div>
<div style="text-align:center"><span style="font-family:'Comic Sans MS';font-size:12px;line-height:14px"><br></span></div><span style="color:black;font-family:'Comic Sans MS';font-size:9pt;line-height:14px"><div style="text-align:center">
Masoud</div><div style="text-align:center"><br></div><div style="text-align:center"><br></div><div style="text-align:center"><br></div></span><span style="color:black;font-family:'Comic Sans MS';font-size:9pt;line-height:14px">--------------</span><br>
<span style="color:black;font-family:'Comic Sans MS';font-size:9pt;line-height:14px">Masoud Nahali, Sharif University of Technology</span></div><div><span style="color:black;font-family:'Comic Sans MS';font-size:9pt;line-height:14px"><a href="mailto:masoud.nahali@gmail.com" target="_blank">masoud.nahali@gmail.com</a>; <a href="http://alum.sharif.edu/%7Em_nahali" target="_blank">a</a><a href="http://lum.sharif.edu/~m_nahali" target="_blank">lum.sharif.edu/~m_nahali</a></span></div>
<div><span style="color:black;font-family:'Comic Sans MS';font-size:9pt;line-height:14px"><p><span><span style="font-size:24.0pt;line-height:115%;font-family:Webdings;color:green">P</span></span><span><span style="font-size:10pt;line-height:115%;color:black"> </span></span><span style="color:black;font-family:'Comic Sans MS';font-size:9pt;line-height:14px"><span><span style="font-size:10.0pt;line-height:115%;color:black"></span></span><span><span style="font-size:7.5pt;line-height:115%;color:green">Save a Tree . . . </span></span></span><span><span style="font-size:7.5pt;line-height:115%;color:green">Please don't print this e-mail
unless you really need to.</span></span></p></span></div>
<div></div><div><br></div><div><br></div>Hongsheng Zhao Wrote :<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"> Thanks a lot for all of the helps on this topic. I've checked and read<br>
most the topics on this ML which are related to "calculating the energy<br>
of isolated atom". It seems this is a tricky thing based on most of<br>
the replies on this similar topics on this ML.<br>
<br>
Based on my understanding, till now, there are still some issues<br>
puzzling me on this topic:<br>
<br>
1- On the manual of pwscf (p37), we can find the following descriptions:<br>
<br>
--------<br>
...<br>
Specify nosym=.true. to avoid generation of additional k-points in low<br>
symmetry<br>
...<br>
-----------<br>
<br>
I want to know, by using nosym=.true., what's the different on the<br>
k-points set obtained with the ones without using this option? For<br>
example, when using a bcc box with a atom at its center, what's the<br>
corresponding k-points set when using nosym=.true.?<br>
<br>
Furthermore, when we say a isolated atom, we mean a zero dimension<br>
system. So the corresponding BZ sampling by using gamma only should be<br>
OK. If so, there is no k-points issues at all - there only one k point<br>
should be considered. Hence I'm confused on the description given by<br>
the manual *avoid generation of additional k-points in low symmetry*.<br>
Any hints?<br>
<br>
2- From this ML, I found the following notes:<br>
<br>
---------<br>
while doing calculations for any zero dimensions you should always set<br>
nosym =.true.<br>
-----------<br>
<br>
Does this mean: when I do the calculations for isolated<br>
atom/molecule/nanocluster and so on, I should always set nosym =.true.?<br>
<br>
3- Prof. Paolo Giannozzi given the following descriptions on some<br>
similar topic in this ML:<br>
<br>
---------<br>
calculations on isolated atoms with plane waves are tricky.<br>
Compare your results with what you get with the atomic code,<br>
in which you can fix occupancies at will<br>
------------<br>
<br>
My issue: how to fix the occupancies based on the results? Any examples?<br>
<br>
What do you mean by saying *at will*? If really is a thing *at will*,<br>
why should we do it?<br>
<br>
Regards<br>
<br>
><br>
> On Thu, Jul 28, 2011 at 10:12 AM, Gabriele Sclauzero<br>
> <<a href="mailto:gabriele.sclauzero@epfl.ch">gabriele.sclauzero@epfl.ch</a> <mailto:<a href="mailto:gabriele.sclauzero@epfl.ch">gabriele.sclauzero@epfl.ch</a>>> wrote:<br>
><br>
><br>
> Il giorno 28/lug/2011, alle ore 01.55, Hongsheng Zhao ha scritto:<br>
><br>
>> Hi all,<br>
>><br>
>> I've learned that when using the a large-cell ( i.e., a big box ) to<br>
>> calculate the energy of isolated atom, it should be done without<br>
>> symmetry. But I'm not sure what's the meaning of the "without<br>
>> symmetry" here and how to ensure this when we using a a big box<br>
>> model<br>
>> for this purpose. Any hints on this issue?<br>
><br>
> This is a FAQ on this forum: have you spent some time browsing the<br>
> archives?<br>
><br>
> OK, I'll give you a hint: "nosym"... but use it with care!<br>
><br>
> HTH<br>
><br>
> GS<br>
><br>
>><br>
>> Regards<br>
>> --<br>
>> Hongsheng Zhao <<a href="mailto:zhaohscas@yahoo.com.cn">zhaohscas@yahoo.com.cn</a><br>
>> <mailto:<a href="mailto:zhaohscas@yahoo.com.cn">zhaohscas@yahoo.com.cn</a>>><br>
>> School of Physics and Electrical Information Science,<br>
>> Ningxia University, Yinchuan 750021, China<br>
>> _______________________________________________<br>
>> Pw_forum mailing list<br>
>> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a> <mailto:<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>
><br>
> ? Gabriele Sclauzero, EPFL SB ITP CSEA<br>
> / PH H2 462, Station 3, CH-1015 Lausanne/<br>
><br><br>
</blockquote></div><br>