<div dir="ltr">Thanks Giuseppe for the advises.<div><br></div><div>I had a closer look at the {alat} and {crystal} option for the atomic position input. For the {crystal} option, the basis vectors are the primitive lattice vectors (a_i), i.e. x a_1 + y a_2 +z a_3. For the {alat} option, the basis are the unit vectors in the Cartesian coordinates (i^{\hat}), i.e. x(a*x^{\hat})+y(b*y^{\hat})+z(c*z^{\hat}), where a, b, c are the lattice constants. Is my understanding correct?</div><div><br></div><div>For the current case, in which I obtain the atomic positions from the Wyckoff positions in the rhombohedral axes, I should use the {crystal} option, right?</div><div><br></div><div>Regards,</div><div>Fermin</div><div><div><br></div></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-06-01 15:55 GMT+08:00 Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Fermin<br>
You are using ultrasoft pseudopotential (I hope that you are aware of<br>
the difference between norm-conserving and ultrasoft...). The ecutrho<br>
variable does not appear in your input, so you are using the default<br>
ecutwfc*4 value which is suitable for NC but not for US. Reasonably<br>
safe values for your pseudopotentials should be ecutwfc=40.0 and<br>
ecutrho =320.0~480.0, but always check! In difficult cases (which<br>
yours should not be) you can lower mixing_beta down to 0.05 0r 0.02,<br>
or also carefully try to use the scf_must_converge=.false. setup.<br>
Remember also that error on stress is one order of magnitude higher<br>
than error on forces, and that it is therefore safe to start a<br>
vc-relax calculation with at least conv_thr=1.0D-7. *BUT* your<br>
calculation does not converge because *you did not check* the<br>
structure... If you try to open your input file with xcrysden you will<br>
touch it with your hands! A hint: do you know the difference between<br>
{alat} and {crystal} coordinates?<br>
HTH<br>
Giuseppe<br>
<br>
Giuseppe Mattioli<br>
ISM-CNR<br>
Italy<br>
<br>
lung Fermin <<a href="mailto:ferminlung@gmail.com">ferminlung@gmail.com</a>> ha scritto:<br>
<div><div class="h5"><br>
> Dear QE community,<br>
><br>
> I want to calculate the phonon diseprsion relation of PbTiO3 in the<br>
> strucutre of R-3c spacegroup. First I did the vc-relax using the<br>
> experimental lattice parameters. However, the scf calculation fails to<br>
> converge within 100 iteration and the job exits with error. I have tried<br>
> changing the mixing_beta to 0.3 and mixing_ndim to 10, but still fails. I<br>
> am a newbie to QE. May someone please give me some advices on this issue?<br>
><br>
> INPUT file:<br>
>  &control<br>
>     calculation='vc-relax',<br>
>     restart_mode='from_scratch',<br>
>     prefix='PbTiO3',<br>
>     outdir = './tmp/',<br>
>     pseudo_dir = '$PSD_DIR',<br>
>     max_seconds = 4.1D+6<br>
>  /<br>
>  &system<br>
>     nat = 10<br>
>     ntyp = 3<br>
>     ibrav = 5<br>
>     celldm(1) = 10.037,<br>
>     celldm(4) = 0.53,<br>
>     ecutwfc =60<br>
>  /<br>
>  &electrons<br>
>    mixing_beta=0.3D0,<br>
>    mixing_ndim=10<br>
>  /<br>
>  &ions<br>
>  /<br>
> &cell<br>
>  /<br>
> ATOMIC_SPECIES<br>
> Pb  207.2  Pb.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
> Ti  47.867  Ti.pbe-spn-rrkjus_psl.1.0.0.UPF<br>
> O  15.999  O.pbe-n-rrkjus_psl.1.0.0.UPF<br>
> ATOMIC_POSITIONS {alat}<br>
> Pb  0.25 0.25 0.25<br>
> Pb  0.75 0.75 0.75<br>
> Ti  0.00 0.00 0.00<br>
> Ti  0.50 0.50 0.50<br>
> O  0.62351014  0.87648986<br>
> 0.25000000<br>
> O  0.37648986  0.12351014  0.75000000<br>
> O  0.87648986  0.25000000  0.62351014<br>
> O  0.12351014  0.75000000  0.37648986<br>
> O  0.25000000  0.62351014  0.87648986<br>
> O  0.75000000  0.37648986  0.12351014<br>
> K_POINTS AUTOMATIC<br>
> 15 15 15 0 0 0<br>
><br>
> Thanks in advance.<br>
> Fermin<br>
<br>
<br>
</div></div>--<br>
********************************************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et ègaux en droits. Les distinctions sociales<br>
ne peuvent être fondèes que sur l'utilitè commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l'homme. Ces droits sont la libertè,<br>
la propriètè, la sùretè et la rèsistance à l'oppression.<br>
********************************************************<br>
<br>
    Giuseppe Mattioli<br>
    CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
    v. Salaria Km 29,300 - C.P. 10<br>
    I 00015 - Monterotondo Stazione (RM)<br>
    Tel <a href="tel:%2B%2039%2006%2090672836" value="+390690672836">+ 39 06 90672836</a> - Fax <a href="tel:%2B39%2006%2090672316" value="+390690672316">+39 06 90672316</a><br>
    E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>><br>
<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></blockquote></div><br></div>