<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> </div>
<div>P.S</div>
<div>The pseudopotentials i'm using are all PBE 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 with 180 Ry cutoff for wavefunction.</div>
<div>I've opened the structure in Xcrysden and it looks OK, no weird things.</div>
<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></div>