<div dir="auto">Dear Paolo,<div dir="auto">Thanks a lot for the clarification. </div><div dir="auto">Best regards, </div><div dir="auto">Mauro. </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il ven 29 mag 2020, 15:14 Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr">On Fri, May 29, 2020 at 1:56 PM Mauro Sgroi <<a href="mailto:maurofrancesco.sgroi@gmail.com" target="_blank" rel="noreferrer">maurofrancesco.sgroi@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="auto"><br><div>My space group is 15. I've the coordinates of the not equivalent atoms from a cif experimental file.</div><div>Those Wyckoff positions are referred to the non-primitive conventional cell with 48 atoms.</div><div>[...] I imagine that pw.x automatically understands that those are referred to the conventional cell and transforms them in the equivalent positions in the primitive (24-atoms) cell.</div><div>Is this correct?</div></div></div></blockquote><div><br></div><div>I think that the code generates all atomic positions  by applying all symmetry operations to the provided Wyckoff positions. Then, it throws away atoms that, modulo lattice translations, overlap. The code always assumes the smallest (primitive) unit cell.<br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="auto"><div>Or should I transform manually the coordinates referring them to the primitive monoclinic axes?</div></div></div></blockquote><div><br></div><div>I don't think you need to do anything (unless you don't get the correct number of atoms, of course): the code expects Wyckoff positions as they are defined in big crystallography books, produces atomic positions in the primitive unit cells as they are internally used.</div><div><br></div><div>Paolo<br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="auto"><div>Thanks a lot and best regards,</div><div>Mauro Sgroi.</div><div dir="auto"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il dom 17 mag 2020, 19:03 Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" target="_blank" rel="noreferrer">p.giannozzi@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Sat, May 16, 2020 at 10:23 PM Mauro Sgroi <<a href="mailto:maurofrancesco.sgroi@gmail.com" rel="noreferrer noreferrer" target="_blank">maurofrancesco.sgroi@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Checking the output of ph.x I found that the celldm parameters are changed with respect to the initial scf calculation (and also the number of atoms in the cell pass from 24 to 30).</div></div></div></div></div></div></blockquote><div><br></div><div class="gmail_quote">With your input and the latest development version, the self-consistent calculation also produces 30 atoms from your Wyckoff positions. I get exactly the same lattice parameters in the phonon and in the scf code. I also get this interesting message:</div></div><div class="gmail_quote">     BEWARE: axis for ibrav=-13 changed, see documentation!</div><div class="gmail_quote">related to this change I did some time ago upon suggestion by don't remember who:<br></div><div class="gmail_quote"><a href="https://gitlab.com/QEF/q-e/-/commit/962a723a9d1e79244c8a2d7468937ab9f29982f4" rel="noreferrer noreferrer" target="_blank">https://gitlab.com/QEF/q-e/-/commit/962a723a9d1e79244c8a2d7468937ab9f29982f4</a></div><div class="gmail_quote">Are you by any chance running different version of the scf and phonon code?<br></div><div class="gmail_quote"><br></div><div class="gmail_quote">Paolo</div><div class="gmail_quote"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Could you please help me to solve this problem?</div><div>Below I'm attaching the input of pw.x, ph.x and the initial output of ph.x.</div><div>Thanks a lot in advance and best regards,</div><div>Mauro Sgroi.</div><div><br></div><div>SCF</div><div><div>&CONTROL</div><div>                       title = 'Li2TiS3_monoclinic'</div><div>                  pseudo_dir = '/workhpc/FCA/FCA_CRF_STRUT/sgroi/DATABASE/ESPRESSO'</div><div>                      prefix = 'LTS_mono'</div><div>                     outdir  = '/workhpc/FCA/FCA_CRF_STRUT/sgroi/tmp/LTS_mono_low_kpt_cell'</div><div>                 calculation = 'scf'</div><div>                   verbosity = 'high'</div><div> /</div><div> &SYSTEM</div><div>                  space_group = 15</div><div>                     uniqueb = .true.</div><div>                   celldm(1) = 11.552652</div><div>                   celldm(2) = 1.73103</div><div>                   celldm(3) = 1.952171</div><div>                   celldm(5) = -0.172475</div><div>                         nat = 8</div><div>                        ntyp = 8</div><div>                     ecutwfc = 52</div><div>                     ecutrho = 575</div><div>                        nbnd = 120</div><div>                 occupations = 'smearing'</div><div>                     degauss = 0.005</div><div>                    smearing = 'cold'</div><div> /</div><div> &ELECTRONS</div><div>            electron_maxstep = 200</div><div>                    conv_thr = 1.0D-12</div><div> /</div><div><br></div><div><br></div><div>ATOMIC_SPECIES</div><div> Li1 6.941  Li.pbe-sl-kjpaw_psl.1.0.0.UPF</div><div> Li2 6.941  Li.pbe-sl-kjpaw_psl.1.0.0.UPF</div><div> Li3 6.941  Li.pbe-sl-kjpaw_psl.1.0.0.UPF</div><div> Ti1 47.867  Ti.pbe-spn-kjpaw_psl.1.0.0.UPF</div><div> Ti2 47.867  Ti.pbe-spn-kjpaw_psl.1.0.0.UPF</div><div> S1  32.06   S.pbe-nl-kjpaw_psl.1.0.0.UPF</div><div> S2  32.06   S.pbe-nl-kjpaw_psl.1.0.0.UPF</div><div> S3  32.06   S.pbe-nl-kjpaw_psl.1.0.0.UPF</div><div><br></div><div>ATOMIC_POSITIONS {crystal_sg}</div><div>Li1 0.3368954  0.8267566  0.9998022</div><div>Li2 0.4164759  0.4164759  0.2500000</div><div>Li3 0.5000000  0.0000000  0.5000000</div><div>Ti1 0.0836168  0.0836168  0.2500000</div><div>Ti2 0.7514196  0.7514196  0.2500000</div><div>S1  0.1929125  0.9744265  0.6332331</div><div>S2  0.3771191  0.1029193  0.1331970</div><div>S3  0.4557540  0.7282681  0.3676504</div><div><br></div><div>K_POINTS automatic</div><div>6  3  3   0 0 0</div></div><div><br></div><div>PHonon</div><div><div>Normal modes for LTS</div><div> &inputph</div><div>  prefix='LTS_mono'</div><div>  outdir  = '/workhpc/FCA/FCA_CRF_STRUT/sgroi/tmp/LTS_mono_low_kpt_cell'</div><div>  tr2_ph=1.0d-14</div><div>  amass(1)=6.941</div><div>  amass(2)=6.941</div><div>  amass(3)=6.941</div><div>  amass(4)=47.867</div><div>  amass(5)=47.867</div><div>  amass(6)=32.06</div><div>  amass(7)=32.06</div><div>  amass(8)=32.06</div><div>  epsil=.false.</div><div>  !lraman=.true.</div><div>  trans=.true.</div><div>  asr=.true.</div><div>  fildyn='dmat.lts'</div><div> /</div><div> 0.0 0.0 0.0</div></div><div><br></div><div>Part of Phonon output</div><div><div>bravais-lattice index     =           13</div><div>     lattice parameter (alat)  =      11.5527  a.u.</div><div>     unit-cell volume          =    2566.1373 (a.u.)^3</div><div>     number of atoms/cell      =           30</div><div>     number of atomic types    =            8</div><div>     kinetic-energy cut-off    =      52.0000  Ry</div><div>     charge density cut-off    =     575.0000  Ry</div><div>     convergence threshold     =      1.0E-14</div><div>     beta                      =       0.7000</div><div>     number of iterations used =            4</div><div>     Exchange-correlation      = PBE ( 1  4  3  4 0 0)</div><div><br></div><div><br></div><div>     celldm(1)=   11.55265  celldm(2)=    0.99956  celldm(3)=    0.00000</div><div>     celldm(4)=   -0.50022  celldm(5)=    0.00000  celldm(6)=    0.00000</div><div><br></div><div>     crystal axes: (cart. coord. in units of alat)</div><div>               a(1) = (  0.5000  0.8655  0.0000 )</div><div>               a(2) = ( -0.5000  0.8655  0.0000 )</div><div>               a(3) = ( -0.3367  0.0000  1.9229 )</div><div><br></div><div>     reciprocal axes: (cart. coord. in units 2 pi/alat)</div><div>               b(1) = (  1.0000  0.5777  0.1751 )</div><div>               b(2) = ( -1.0000  0.5777 -0.1751 )</div><div>               b(3) = (  0.0000 -0.0000  0.5200 )</div><div><br></div><div><br></div><div>     Atoms inside the unit cell:</div><div><br></div><div>     Cartesian axes</div><div><br></div><div>     site n.  atom      mass           positions (alat units)</div><div>        1     Li1  6.9410   tau(    1) = (   -0.49974    0.56563    1.92254  )</div><div>        2     Li1  6.9410   tau(    2) = (   -0.00531    0.56563    0.96184  )</div><div>        3     Li1  6.9410   tau(    3) = (    0.16304    1.16540    0.00038  )</div><div>        4     Li1  6.9410   tau(    4) = (   -0.33139    1.16540    0.96108  )</div><div>        5     Li2  6.9410   tau(    5) = (    0.33230    0.72093    0.48073  )</div><div>        6     Li2  6.9410   tau(    6) = (   -0.50065    0.72093    0.48073  )</div><div>        7     Li2  6.9410   tau(    7) = (   -0.16900    0.14458    1.44219  )</div><div>        8     Li2  6.9410   tau(    8) = (   -0.33605    0.14458    1.44219  )</div><div>        9     Li3  6.9410   tau(    9) = (   -0.16835    0.86552    0.96146  )</div><div>       10     Li3  6.9410   tau(   10) = (    0.00000    0.86552    0.00000  )</div><div>       11     Ti1 47.8670   tau(   11) = (   -0.00056    0.14474    0.48073  )</div><div>       12     Ti1 47.8670   tau(   12) = (   -0.16779    0.14474    0.48073  )</div><div>       13     Ti1 47.8670   tau(   13) = (    0.16386    0.72077    1.44219  )</div><div>       14     Ti1 47.8670   tau(   14) = (   -0.66891    0.72077    1.44219  )</div><div>       15     Ti2 47.8670   tau(   15) = (    0.16724    0.43521    0.48073  )</div><div>       16     Ti2 47.8670   tau(   16) = (   -0.33559    0.43521    0.48073  )</div><div>       17     Ti2 47.8670   tau(   17) = (   -0.00395    0.43030    1.44219  )</div><div>       18     Ti2 47.8670   tau(   18) = (   -0.50111    0.43030    1.44219  )</div><div>       19     S1  32.0600   tau(   19) = (   -0.52030    0.82125    1.21765  )</div><div>       20     S1  32.0600   tau(   20) = (    0.01525    0.82125    1.66672  )</div><div>       21     S1  32.0600   tau(   21) = (    0.18360    0.90978    0.70526  )</div><div>       22     S1  32.0600   tau(   22) = (   -0.35195    0.90978    0.25620  )</div><div>       23     S2  32.0600   tau(   23) = (   -0.16773    1.04367    0.25613  )</div><div>       24     S2  32.0600   tau(   24) = (   -0.00062    1.04367    0.70533  )</div><div>       25     S2  32.0600   tau(   25) = (   -0.16897    0.68736    1.66679  )</div><div>       26     S2  32.0600   tau(   26) = (   -0.33608    0.68736    1.21758  )</div><div>       27     S3  32.0600   tau(   27) = (   -0.16803    0.39514    0.70696  )</div><div>       28     S3  32.0600   tau(   28) = (   -0.00032    0.39514    0.25450  )</div><div>       29     S3  32.0600   tau(   29) = (   -0.16867    1.33589    1.21595  )</div><div>       30     S3  32.0600   tau(   30) = (   -0.33638    1.33589    1.66842  )</div></div><div><br></div><div><br></div></div></div></div></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>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></div></div></div></div></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>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></div></div></div></div></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>