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">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;">
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>_______________________________________________<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>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>