<div dir="ltr"><div><div><div><div>It is quite natural if you optimize the cell parameter along the z-direction also in your vc-relax. The system will try to minimize its energy by interacting with the periodic images along z, i.e., tending to form bulk. In order to do so, the value of the cell parameter will reduce along z. <br><br></div>If you are interested in finding the lattice parameters in the plane of your sheet, you should do a vc-relax calculation keeping the z fixed.<br><br></div>Since you have a 4x4 unit cell, you can also try to do manually by varying in-plane lattice parameter, plot the total energy vs a and find the minima of the curve.<br><br></div>With regards,<br></div>Prasenjit<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 29 August 2017 at 10:37,  <span dir="ltr"><<a href="mailto:l.eslami@srbiau.ac.ir" target="_blank">l.eslami@srbiau.ac.ir</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">On 2017-08-28 11:40, Lorenzo Paulatto wrote:<br>
> Dear Leila,<br>
><br>
> if the total energy keeps decreasing, that it is completely normal: it<br>
> means that you started from an energy saddle point. It may however<br>
> indicate that your initial crystal structure was wrong.<br>
><br>
> kind regards<br>
><br>
><br>
> On 28/08/17 07:22, <a href="mailto:l.eslami@srbiau.ac.ir">l.eslami@srbiau.ac.ir</a> wrote:<br>
>> Hi all,<br>
>><br>
>> We are trying to perform vc-relax calculation for a single layer ZnO<br>
>> sheet in the presence of an oxygen vacancy in a 4*4 suppercell. We<br>
>> have<br>
>> used espresso-5.4.0. with different psudopotentials such as GGA,<br>
>> NORMCONS and non relativistic Pseudopotentials: "O.pbe-mt_fhi.upf and<br>
>> Zn.pbe-mt_fhi.upf,  as well as GBRV pseudopotentials:<br>
>> “o_pbe_v1.2.uspp.F.UPF” and “zn_pbe_v1.uspp.F.UPF”. The calculations<br>
>> do<br>
>> not return any error. However, convergences of force and pressure do<br>
>> not<br>
>> occur even after long time running of the calculations. In fact, the<br>
>> code starts to decrease force and pressure but after a couple of<br>
>> relaxation steps they increase. I would be grateful if anyone could<br>
>> kindly give me any piece of advice to solve the problem.<br>
>> many thanks,<br>
>><br>
>> Leila Eslami<br>
>> Department of physics,<br>
>> Islamic Azad University,<br>
>> Tehran, Iran<br>
>> ______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
<br>
<br>
<br>
<br>
</div></div>Dear Lorenzo,<br>
thank you for your kind attention. I checked the total energy of the<br>
system and it keeps decreasing. However, there is another ambiguous<br>
point in the results. In fact, as our system is a sheet and we set 15<br>
Angstrom distance between subsequent layers in z-direction, I expect<br>
that the z-components of atomic positions do not change during relation<br>
process. But after a couple of relation steps this happens. To me it<br>
doesn't seem normal!!!!<br>
<br>
regards,<br>
<div class="HOEnZb"><div class="h5">Leila Eslami<br>
  Department of physics,<br>
  Islamic Azad University,<br>
  Tehran, Iran<br>
<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a></div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">PRASENJIT GHOSH,<br>IISER Pune,<br>Dr. Homi Bhabha Road, Pashan<br>Pune, Maharashtra 411008, India<br><br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790</div></div>
</div>