Dear Paolo,<br>         Thanks for your answer. However, I still cannot quite understand it. In my script I defined the plane as (001), which is ,<br>e1(1) =1.0, e1(2)=0.0, e1(3) = 0.0;   <br>e2(1) =0.0, e2(2)=1.0, e2(3) = 0.0;<br>
I didn't define x(i). In my scf file, the unit cell is just like quarter, is that the reason for my quarter like charge density contour?<br><br>           Thanks for your patience. <br><br>Chengyu Yang<br>UCF PhD student.<br>
<br><br>