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<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">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 class="h5">
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">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>