<div dir="ltr"><div><div>Dear Lorenzo Dona,<br><br></div>Perhaps you could try to use a visualization software to verify the structure and export it in .xyz format. You could use VESTA to make the structure and visualize it. VESTA can read Wyckoff positions. Furthermore, it can export the positions of the atoms in .xyz format. You can use the coordinates of the .xyz file as input. This way, you avoid using space_group and crystal_sg.<br><br></div>Dae Kwang Jun<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 6, 2016 at 7:24 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">You need:<br>
- the space group number<br>
- the cell parameters (for the rhombohedral lattice, not the hexagonal one)<br>
- the Wyckoff positions. If the latter are given in hexagonal axis, as<br>
they are here:<br>
  <a href="http://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-wp-list" rel="noreferrer" target="_blank">http://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-wp-list</a><br>
  use "rhombohedral=.false."<br>
This is all I know about the specific case of rhombohedral groups<br>
<span class="HOEnZb"><font color="#888888"><br>
Paolo<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
On Sat, Jun 4, 2016 at 8:01 PM, Lorenzo Donà <<a href="mailto:lorechimica91@hotmail.it">lorechimica91@hotmail.it</a>> wrote:<br>
> Dear Paolo I have tried a lot of inputs with space_group=148 keyword but the<br>
> structure found is incorrect…<br>
> Or better is correct the number of atom 54 in the cell but the position is<br>
> incorrect.<br>
> Please can you help me to understand the correct input for my system, can<br>
> you tell me how to create the correct structure input using the keyword<br>
> space_group.<br>
> I must tell you that i found problem only for this system for cubic or<br>
> ORTHORHOMBIC system I found NO problem.<br>
> Thanks a lot if you can help me with my problem.<br>
> Dearly<br>
> lorenzo<br>
><br>
> Il giorno 31 mag 2016, alle ore 21:58, Paolo Giannozzi<br>
> <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> ha scritto:<br>
><br>
> For the rhombohedral (or trigonal) trigonal lattice, cosab=-0.5 means<br>
> that the angle between any two of the three primitive vectors is 120,<br>
> that is, they lie on a plane. This is why the code complains. By the<br>
> way, here a converter between rhombohedral and hexagonal cells:<br>
> <a href="http://qpeng.org/tools/r2h.html" rel="noreferrer" target="_blank">http://qpeng.org/tools/r2h.html</a>, and here a  picture (fig.5):<br>
> <a href="http://xrayweb.chem.ou.edu/notes/symmetry.html" rel="noreferrer" target="_blank">http://xrayweb.chem.ou.edu/notes/symmetry.html</a><br>
><br>
><br>
> On Tue, May 31, 2016 at 9:57 AM, Lori 91 <<a href="mailto:lorechimica91@hotmail.it">lorechimica91@hotmail.it</a>> wrote:<br>
><br>
> Thanks to replay Davide.<br>
> cosab=-0.5 give me an error.<br>
> Thanks a lot dearly<br>
> Lorenzo<br>
><br>
> Inviato da iPhone<br>
><br>
> Il giorno 31 mag 2016, alle ore 09:16, Tiana Davide <<a href="mailto:davide.tiana@epfl.ch">davide.tiana@epfl.ch</a>><br>
> ha scritto:<br>
><br>
> Ciao Lorenzo<br>
><br>
> I am not sure but one problem could be that you're using cosab 0.5 (60<br>
> degrees) instead of cosab=-0.5 (120 degrees)<br>
> another one (again just a guess) is your cut-off seems to be really low even<br>
> for an ultrasoft pseudopotential.<br>
><br>
> anyway, you should try with this electron setting<br>
><br>
>  mixing_mode='local-TF',<br>
>  mixing_beta=0.6,<br>
><br>
> I've not problem in converging it (using different pseudo though)<br>
><br>
> Cheers<br>
> Davide<br>
><br>
> P.S. the last test you can do is to use a bit of smearing, for instance:<br>
><br>
>  occupations='smearing',<br>
>  smearing='marzari-vanderbilt'<br>
>  degauss=0.05<br>
> ________________________________________<br>
> ------------------------------<br>
><br>
> Message: 2<br>
> Date: Sun, 29 May 2016 21:11:49 +0200<br>
> From: Lorenzo Don? <<a href="mailto:lorechimica91@hotmail.it">lorechimica91@hotmail.it</a>><br>
> Subject: [Pw_forum] CPO-27-Zn scf not converge<br>
> To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br>
> Message-ID: <BLU436-SMTP367637BD78842FCA03B145CB440@phx.gbl><br>
> Content-Type: text/plain; charset="us-ascii"<br>
><br>
> Dear all can you help me or give me some tips to make scf convergence on<br>
> this calculation:<br>
><br>
> &control<br>
>   pseudo_dir='/home/utente/Scaricati/espresso-5.4.0/pseudo/'<br>
>   outdir = './',<br>
>   wf_collect = .true.,<br>
>   verbosity= high,<br>
> /<br>
> &system<br>
>   ibrav=  5,<br>
>   a= 25.226, cosab= 0.5,  space_group = 148,<br>
>   nat=  9, ntyp= 4,<br>
>   ecutwfc =15.0,<br>
>   occupations='fixed', smearing='gauss',<br>
>    input_dft= PBE,<br>
> /<br>
> &electrons<br>
>   mixing_mode='plain'<br>
>   diagonalization='david'<br>
>   mixing_beta = 0.01<br>
>   conv_thr= 1.0e-4<br>
>   electron_maxstep=500,<br>
> /<br>
> ATOMIC_SPECIES<br>
> Zn 65.409 Zn.pbe-van.UPF<br>
> C 12.010  C.pbe-van_ak.UPF<br>
> O 16.00   O.pbe-van_ak.UPF<br>
> H 1.00    H.pbe-van_ak.UPF<br>
> ATOMIC_POSITIONS crystal_sg<br>
> C   7.127430447398E-02 -2.455696438145E-01  7.991093030472E-02<br>
> C   1.197468901477E-01 -2.066445045066E-01  2.152449225957E-01<br>
> C   2.117353032357E-01 -1.117780631842E-01  2.553368090727E-01<br>
> C   1.646071563067E-01 -1.528403540257E-01  1.440875844425E-01<br>
> H   1.629458066146E-01 -1.419696182443E-01 -4.760745519579E-03<br>
> Zn  -1.943259345214E-02  2.834354147594E-01  3.060064447988E-01<br>
> O   7.229455751443E-02 -2.290060932141E-01 -9.085450716092E-02<br>
> O   3.018234203718E-02 -2.947952345384E-01  1.422041978763E-01<br>
> O   2.519030126218E-01 -6.097111538485E-02  1.721974260943E-01<br>
> K_POINTS gamma<br>
><br>
> Thanks a lot to help me<br>
><br>
> dearly<br>
><br>
> lorenzo<br>
><br>
><br>
><br>
> _______________________________________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
><br>
><br>
><br>
> --<br>
> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
> Phone +39-0432-558216, fax +39-0432-558222<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
><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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
<br>
<br>
--<br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></div></div></blockquote></div><br></div>