<div dir="ltr">Thank you both for your responses. <div><br></div><div>Paolo to clarify, only non-symmetric systems can be relaxed? </div><div><br></div><div>Regards </div><div><br></div><div>Zarah</div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Fri, Feb 28, 2014 at 3:07 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
If forces are zero by symmetry, there is nothing to relax<br>
<br>
P.<br>
<div class="HOEnZb"><div class="h5"><br>
On Fri, 2014-02-28 at 13:06 +0800, Zarah Baiyee wrote:<br>
> Hi All,<br>
><br>
><br>
> I am a new user to Quantum Espresso, therefore this may be a bit<br>
> trivial but I any help would be very welcome.<br>
><br>
><br>
> I am attempting a relax calculation using espresso version 5.0.2, for<br>
> the geometry optimisation of a perovskite structure. For the input<br>
> below, the calculation runs, however does not compute any forces and<br>
> convergence without geometry optimisation to the original structure,<br>
> with 0 bfgs steps. I am able to perform an expected relaxation if<br>
> alter the initial atomic positions to asymmetric: (i.e.)  an oxygen<br>
> molecule from:  O (0.500  0.5  0.0) to O(0.499 0.5 0.0).<br>
><br>
><br>
> Please may someone shed some light on to why this is the case, and if<br>
> there is anything wrong with my input.<br>
><br>
><br>
> Thanks<br>
><br>
><br>
> Zarah<br>
> <a href="mailto:Zbaiyee@gmail.com">Zbaiyee@gmail.com</a><br>
><br>
><br>
> &CONTROL<br>
>   calculation  = "relax",<br>
>   prefix       = "CO",<br>
>   pseudo_dir   = "/Users/zarahbaiyee/espresso-5.0.2/pseudo",<br>
>   outdir       = "/Users/zarahbaiyee/tmp",<br>
>   forc_conv_thr = 1.0D-5<br>
> /<br>
> &SYSTEM<br>
>   ibrav     = 0,<br>
>   nat       = 5,<br>
>   ntyp      = 3,<br>
>   ecutwfc   = 60,<br>
>   ecutrho   = 420,<br>
> /<br>
> &ELECTRONS<br>
>   conv_thr    = 1.D-7,<br>
>   mixing_beta = 0.7D0,<br>
>   electron_maxstep = 400<br>
> /<br>
> &IONS<br>
> /<br>
> CELL_PARAMETERS bohr<br>
>  7.60  0.0  0.0<br>
>  0.0  7.60  0.0<br>
>  0.0  0.0  7.60<br>
> ATOMIC_SPECIES<br>
>  Ba  137.32 Ba.pbe-mt_fhi.UPF<br>
>  Fe  55.84 Fe.pbe-mt_fhi.UPF<br>
>  O   15.99 O.pbe-kjpaw.UPF<br>
> ATOMIC_POSITIONS {crystal}<br>
> Ba  0.000  0.0  0.0  0 0 0<br>
> Fe  0.500  0.5  0.5<br>
> O   0.000  0.5  0.5<br>
> O   0.500  0.0  0.5<br>
> O   0.500  0.5  0.0<br>
><br>
> K_POINTS {automatic}<br>
> 9 9 9 0 0 0<br>
><br>
><br>
</div></div><div class="im HOEnZb">> _______________________________________________<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>
<br>
</div><span class="HOEnZb"><font color="#888888">--<br>
Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</a><br>
</font></span><div class="HOEnZb"><div class="h5"><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>
</div></div></blockquote></div><br></div>