Dear prasenajit<br>thank you for your valuable inputs <br><br><div class="gmail_quote">On Fri, Sep 11, 2009 at 2:43 PM, Prasenjit Ghosh <span dir="ltr"><<a href="mailto:prasenjit.jnc@gmail.com">prasenjit.jnc@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">For Ti atom, the electronic configuration is [<b>Ar</b>].<b>3d</b><sup>2</sup>.<b>4s</b><sup>2</sup>. So the magnetization should be 2.0 Bohr magneton due to the two electrons in the d-orbital, which is not what you are getting from your calculation.<br>
<br>Also your smearing contribution is -0.00972779 Ry (about 0.1 eV), which I would say very high. If you would have find the correct magnetic ground state, then this should have been 0.0 Ry.<br><br>As far as the bulk system is concerned, I'm not sure. It depends how Ti bulk behaves. You can search in the literature to find out what values of degauss other people have used.<br>
<br>Mixing beta doesn't effect the energy of the system. It just helps you in convergence.<div><div></div><div class="h5"><br><br>Prasenjit.<br><br><div class="gmail_quote">2009/9/11 udayagiri sai babu <span dir="ltr"><<a href="mailto:udayagiri3@gmail.com" target="_blank">udayagiri3@gmail.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;">Dear prasenajit this is the part of my output for the posted input file. (smearing contribution is there! though it is small ) is the degauss value too high for isolated atom or it is high for bulk calculation also? because i use around 0.01 Ry value for all my calculations. Does the mixing beta value effect the energy of the system?<br>
<br> one-electron contribution = -155.91791626 Ry<br> hartree contribution = 74.49216944 Ry<br> xc contribution = -14.37611953 Ry<br> ewald contribution = -20.42854186 Ry<br> smearing contrib. (-TS) = -0.00972779 Ry<br>
<br> total magnetization = 3.24 Bohr mag/cell<br> absolute magnetization = 3.24 Bohr mag/cell<br><br> convergence has been achieved in 34 iterations<div><div></div><div><br><br><div class="gmail_quote">
On Fri, Sep 11, 2009 at 2:06 PM, Prasenjit Ghosh <span dir="ltr"><<a href="mailto:prasenjit.jnc@gmail.com" target="_blank">prasenjit.jnc@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear Udayagiri,<br><br>First of all it's not clear whether you have any question other than whether your input file is correct.<br>
<br>Now some comments about your input file:<br><br>(a) while doing calculations for any zero dimensions you should always set nosym =.true.<br>
(b) Your value of degauss is too high.....you should use some value around 0.001 Ry. According to my experience, doing cal. for a single atom is tricky, specially when you are using a smearing. Even that doesn't confirm you will get the correct ground state for the atom. So what you should do is at the end of the calculation check the values of the total magnetization, absolute magnetization and the smearing contribution. The last term should be exactly zero.<br>
(c) You don't need to set tstress and tprnfor true. <br>(c) mixing_beta= 0.7 too large.....I would put something smaller....0.3<br>(d) I hope you have checked the convergence of ecut rho & ecutwfc.<br><br>Prasenjit. <br>
<br><div class="gmail_quote">2009/9/11 udayagiri sai babu <span dir="ltr"><<a href="mailto:udayagiri3@gmail.com" target="_blank">udayagiri3@gmail.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;">
<div><div></div><div>
Hi everybody<br>I want to calculate the cohesive energy of Titanium for that i need to calculate the energy of an isolated atom previously in the forum somebody has suggested me to do spin polarisation calculation for the isolated atom. I have written an input file accordingly some one can see and tell me if my input script is true. I have read the forum archives on this topic but they are not exactly answering my question. Can somebody please help me. <br>
<b><br>Input file for isolated atom </b> <br><br>&control<br>calculation = 'scf',<br>restart_mode='from_scratch',<br>outdir='/home/sai/temp',<br>prefix='titanium'<br>pseudo_dir='/home/sai/softwares/espresso-4.0.1/pseudo/',<br>
tstress = .TRUE.,<br>tprnfor = .TRUE.,<br>/<br>&system<br>ibrav=1<br>celldm(1)=20,<br>nat=1,<br>ntyp=1<br>ecutwfc=25,<br>ecutrho=250,<br>occupations='smearing',<br>smearing='mv',<br>degauss=0.010000<br>
nspin=2<br>starting_magnetization(1)=0.2<br>/<br>&electrons<br>conv_thr=1.D-8,<br>diagonalization='david',<br>mixing_mode='plain',<br>mixing_beta= 0.7,<br>/<br>ATOMIC_SPECIES<br>Ti 47.867 Ti.pbe-sp-van_ak.UPF<br>
ATOMIC_POSITIONS (crystal)<br>Ti 0.00000000 0.00000000 0.00000000<br>K_POINTS GAMMA<br clear="all"><font color="#888888"><br>-- <br>U.Saibabu<br>PhD student,<br>Deformation mechanisms modeling group,<br>Materials engineering department,<br>
IISc Bangalore,<br>India.<br><br>
</font><br></div></div>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">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><br clear="all"><br>-- <br>PRASENJIT GHOSH,<br>POST-DOC,<br>ROOM NO: 265, MAIN BUILDING,<br>CM SECTION, ICTP,<br>STRADA COSTERIA 11,<br>TRIESTE, 34104,<br>ITALY<br>PHONE: +39 040 2240 369 (O)<br>
+39 3807528672 (M)<br>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">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><br clear="all"><br>-- <br>U.Saibabu<br>PhD student,<br>Deformation mechanisms modeling group,<br>Materials engineering department,<br>IISc Bangalore,<br>India.<br><br>
</div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">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><br clear="all"><br>-- <br>PRASENJIT GHOSH,<br>POST-DOC,<br>ROOM NO: 265, MAIN BUILDING,<br>CM SECTION, ICTP,<br>STRADA COSTERIA 11,<br>TRIESTE, 34104,<br>ITALY<br>PHONE: +39 040 2240 369 (O)<br>
+39 3807528672 (M)<br>
</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://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>U.Saibabu<br>PhD student,<br>Deformation mechanisms modeling group,<br>Materials engineering department,<br>IISc Bangalore,<br>India.<br><br>