<div dir="ltr"><div>Hi</div><div><br></div><div>
<div>just guessing: the space-group machinery takes the provided Wyckoff positions, applies crystal symmetries, throws out translationally equivalent positions. It works also if positions are provided in the form (x,y,z) (wrt crystal axis), but you have to ensure that the Wyckoff site symmetry is correct. I suspect that one or more of the Ge positions are not exact Wyckoff positions but very close to them. <br></div></div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Sep 11, 2020 at 2:25 PM José Carlos Conesa Cegarra <<a href="mailto:jcconesa@icp.csic.es">jcconesa@icp.csic.es</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all,<br>
<br>
I have found a problem with the number of atoms resulting in a <br>
computation with pw.x. The essentials of the geometry specification are:<br>
<br>
&SYSTEM<br>
space_group=155, rhombohedral=.TRUE.<br>
A=8.7018, B=8.7018, C=8.7018<br>
cosAB=-0.0032742, cosAC=-0.0032742, cosBC=-0.0032742<br>
nat=13, ntyp=4<br>
...<br>
<br>
/<br>
<br>
ATOMIC_POSITIONS crystal_sg<br>
V 0.00000 0.00000 0.00000<br>
Ge -0.75399 -0.99800 -0.74920<br>
N -0.73084 -0.51728 -0.00123<br>
N -0.01635 -0.01149 -0.23613<br>
N -0.52047 -0.23372 -0.50080<br>
N -0.22951 -0.73368 -0.74864<br>
N -0.26226 -0.48504 -0.98625<br>
Sn -0.12491 -0.62544 -0.12283<br>
Ge 0.00000 -0.49511 -0.50488<br>
Ge -0.50000 -0.75534 -0.24466<br>
Sn -0.50000 -0.24993 -0.75007<br>
N -0.76700 -0.76700 -0.76700<br>
Sn -0.62767 -0.62767 -0.62767<br>
<br>
With this input I expect to have 56 atoms, as corresponds to a V-doped <br>
(Ge,Sn) spinel nitride having almost the same number of Ge and Sn atoms <br>
(in fact, there is one Sn atom less) as I have verified independently. <br>
However, I obtain at the beginning of the pw.x calculation the following <br>
list of atoms:<br>
<br>
site n. atom positions (alat units)<br>
1 V tau( 1) = ( 0.0000000 0.0000000 0.0000000 )<br>
2 Ge tau( 2) = ( -0.0033926 -0.2015179 0.2870436 )<br>
3 Ge tau( 3) = ( 0.1762159 0.0978209 0.2870436 )<br>
4 Ge tau( 4) = ( -0.1728234 0.1036970 0.2870436 )<br>
5 Ge tau( 5) = ( -0.0033926 0.2015179 1.4393268 )<br>
6 Ge tau( 6) = ( 0.1762159 -0.0978209 1.4393268 )<br>
7 Ge tau( 7) = ( -0.1728234 -0.1036970 1.4393268 )<br>
8 N tau( 8) = ( -0.5167561 -0.1236930 1.0074235 )<br>
9 N tau( 9) = ( 0.3654993 -0.3856774 1.0074235 )<br>
10 N tau( 10) = ( 0.1512567 0.5093704 1.0074235 )<br>
11 N tau( 11) = ( -0.5167561 0.1236930 0.7189470 )<br>
12 N tau( 12) = ( 0.3654993 0.3856774 0.7189470 )<br>
13 N tau( 13) = ( 0.1512567 -0.5093704 0.7189470 )<br>
14 N tau( 14) = ( 0.1556621 0.0938462 1.5744671 )<br>
15 N tau( 15) = ( -0.1591043 0.0878842 1.5744671 )<br>
16 N tau( 16) = ( 0.0034422 -0.1817305 1.5744671 )<br>
17 N tau( 17) = ( 0.1556621 -0.0938462 0.1519033 )<br>
18 N tau( 18) = ( -0.1591043 -0.0878842 0.1519033 )<br>
19 N tau( 19) = ( 0.0034422 0.1817305 0.1519033 )<br>
20 N tau( 20) = ( -0.0139315 0.2264700 1.0041779 )<br>
21 N tau( 21) = ( -0.1891630 -0.1253001 1.0041779 )<br>
22 N tau( 22) = ( 0.2030945 -0.1011699 1.0041779 )<br>
23 N tau( 23) = ( -0.0139315 -0.2264700 0.7221925 )<br>
24 N tau( 24) = ( -0.1891630 0.1253001 0.7221925 )<br>
25 N tau( 25) = ( 0.2030945 0.1011699 0.7221925 )<br>
26 N tau( 26) = ( 0.3676808 -0.2000458 0.7412862 )<br>
27 N tau( 27) = ( -0.0105956 0.4184438 0.7412862 )<br>
28 N tau( 28) = ( -0.3570852 -0.2183980 0.7412862 )<br>
29 N tau( 29) = ( 0.3676808 0.2000458 0.9850842 )<br>
30 N tau( 30) = ( -0.0105956 -0.4184438 0.9850842 )<br>
31 N tau( 31) = ( -0.3570852 0.2183980 0.9850842 )<br>
32 N tau( 32) = ( 0.5127756 0.1138545 0.7287873 )<br>
33 N tau( 33) = ( -0.3549887 0.3871495 0.7287873 )<br>
34 N tau( 34) = ( -0.1577869 -0.5010040 0.7287873 )<br>
35 N tau( 35) = ( 0.5127756 -0.1138545 0.9975831 )<br>
36 N tau( 36) = ( -0.3549887 -0.3871495 0.9975831 )<br>
37 N tau( 37) = ( -0.1577869 0.5010040 0.9975831 )<br>
38 Sn tau( 38) = ( -0.0014732 -0.4102001 1.2238930 )<br>
39 Sn tau( 39) = ( 0.3559803 0.2038242 1.2238930 )<br>
40 Sn tau( 40) = ( -0.3545071 0.2063759 1.2238930 )<br>
41 Sn tau( 41) = ( -0.0014732 0.4102001 0.5024774 )<br>
42 Sn tau( 42) = ( 0.3559803 -0.2038242 0.5024774 )<br>
43 Sn tau( 43) = ( -0.3545071 -0.2063759 0.5024774 )<br>
44 Ge tau( 44) = ( -0.3506754 0.2104528 0.5754626 )<br>
45 Ge tau( 45) = ( -0.0069197 -0.4089202 0.5754626 )<br>
46 Ge tau( 46) = ( 0.3575951 0.1984674 0.5754626 )<br>
47 Ge tau( 47) = ( 0.3575880 0.1984633 0.5754510 )<br>
48 Ge tau( 48) = ( -0.0069197 -0.4089120 0.5754510 )<br>
49 Ge tau( 49) = ( -0.3506683 0.2104487 0.5754510 )<br>
50 Ge tau( 50) = ( -0.1808480 0.3132379 0.8631852 )<br>
51 Ge tau( 51) = ( 0.3616960 -0.0000000 0.8631852 )<br>
52 Ge tau( 52) = ( -0.1808480 -0.3132379 0.8631852 )<br>
53 Sn tau( 53) = ( 0.1771154 -0.3067729 0.8631852 )<br>
54 Sn tau( 54) = ( -0.3542309 -0.0000000 0.8631852 )<br>
55 Sn tau( 55) = ( 0.1771154 0.3067729 0.8631852 )<br>
56 N tau( 56) = ( 0.0000000 -0.0000000 0.4022443 )<br>
57 N tau( 57) = ( -0.0000000 -0.0000000 1.3241261 )<br>
58 Sn tau( 58) = ( -0.0000000 -0.0000000 0.6427795 )<br>
59 Sn tau( 59) = ( 0.0000000 -0.0000000 1.0835909 )<br>
<br>
i.e. there are 3 Ge atoms more than expected (the count for the other <br>
atoms is OK). And one can see that atoms 44/49, 45/48 and 46/47 are at <br>
nearly the same positions. But of course with these positions added <br>
automatically by pw.x the calculation cannot progress. How can this have <br>
happened? There is something wrong here.<br>
<br>
Best regards,<br>
<br>
-- <br>
José C. Conesa<br>
Research Professor<br>
Instituto de Catálisis y Petroleoquímica, CSIC<br>
Marie Curie, 2; Campus de Cantoblanco<br>
28028 Madrid (Spain)<br>
Telef. +34 915854766<br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><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>