<br><br><div class="gmail_quote">On Wed, Oct 31, 2012 at 6:14 PM, david furman <span dir="ltr"><<a href="mailto:sirok4@gmail.com" target="_blank">sirok4@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>Hi all,</div>
<div> </div>
<div>I sucsessfully calculated equation of state for CaF2 by varying celldm(1) in a loop.</div>
<div>then fitting the EV data to murnaghan EOS to get a parabola like plot.</div>
<div> </div>
<div>Now I try to calculate EV data for CaTiO3 (Pnma space group) by the same method.</div>
<div>But, no matter how large the celldm(1) is (keeping the b/a and c/a fixed), the energy still goes down instead of going up, with no parabola like behavior.</div>
<div>I assume this is wrong, because when distances between atoms are large (expansion of unit cell), this state should be higher</div>
<div>in energy than the equilibrium volume state, just like when compressing the unitcell the energy goes up because of repulsion. </div>
<div> </div>
<div>Can anybody tell me why is this happening ? Could it be a bug ?</div></div></blockquote><div>Before reaching any excellent conclusion can you provide your input file for system. <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div> </div>
<div>P.S</div>
<div>The pseudopotentials i'm using are all PBE</div></div></blockquote><div>whether it is NCPP or USPP? <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div> and verified with other systems so they should be OK.</div>
<div>There are 20 atoms in the primitive unit cell. I'm using a 8x8x8 k-point grid</div></div></blockquote><div>Have you check your k-point convergence with the total energy? <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div> with 180 Ry cutoff for wavefunction.</div></div></blockquote><div>It is very large for USPP and check it  for both pseudopotential convergence .<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div>I've opened the structure in Xcrysden and it looks OK, no weird things.</div><span class="HOEnZb"><font color="#888888">
<div> </div>
<div><br clear="all"><br>-- <br></div>
<div dir="ltr">
<p style="TEXT-ALIGN:left;DIRECTION:ltr"><b><span style="FONT-FAMILY:'Candara','sans-serif';COLOR:#00b0f0">David Furman</span></b><span dir="ltr"></span></p>
<div><font size="3" face="'Times New Roman'"><i><br></i></font>
<div><span style="FONT-FAMILY:'Times New Roman';FONT-SIZE:medium"><br></span></div></div></div><br></font></span></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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br>Thanks and Regards<br>Bramha Prasad Pandey<br>Indian School of Mines(ISM)<br>
Dhanbad, INDIA.<br><br>