<div dir="ltr">Dear all,<br><br>I have been trying to perform an electron-phonon calculation for a hexagonal structure but the calculation crashes at second point on a 2x2x2 q-mesh in the elphon subroutine with the error &quot;cannot remap grid on k-point list&quot;&nbsp; do n=1,nkh<br>
<br>&nbsp;&nbsp;&nbsp;&nbsp; do nk=1,nkBZ<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (eqBZ(nk) == n) go to 20<br>&nbsp;&nbsp;&nbsp;&nbsp; end do<br>&nbsp;&nbsp;&nbsp;&nbsp; !&nbsp; this failure of the algorithm may indicate that the displaced grid<br>&nbsp;&nbsp;&nbsp;&nbsp; !&nbsp; (with k1,k2,k3.ne.0) does not have the full symmetry of the lattice<br>
&nbsp;&nbsp;&nbsp;&nbsp; call errore(&#39;lint&#39;,&#39;cannot remap grid on k-point list&#39;,n)<br>20&nbsp;&nbsp; continue<br>&nbsp; end do<br><br>The problem is that the code finds more irreducible k-points than there should be. In this example k(27) is the same as k(25) and k(29) is the same as k(9). So, the correct number of k points should be 26 not 30.<br>
&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp; number of k points=&nbsp;&nbsp; 30&nbsp; gaussian broad. (Ry)=&nbsp; 0.0050&nbsp;&nbsp;&nbsp;&nbsp; ngauss =&nbsp;&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cart. coord. in units 2pi/a_0<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 1) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000), wk =&nbsp;&nbsp; 0.0312500<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 2) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 3) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.0625000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 4) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 5) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.0312500<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 6) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.0000000&nbsp; -0.6097561), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 7) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp;&nbsp; 0.0000000), wk =&nbsp;&nbsp; 0.1875000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 8) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp;&nbsp; 9) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.0937500<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 10) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 11) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.1875000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 12) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp; -0.6097561), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 13) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.5773503&nbsp;&nbsp; 0.0000000), wk =&nbsp;&nbsp; 0.0937500<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 14) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.5773503&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 15) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.5773503&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.1875000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 16) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.5773503&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 17) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.5773503&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.0937500<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 18) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.5773503&nbsp; -0.6097561), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 19) = (&nbsp;&nbsp; 0.2500000&nbsp;&nbsp; 0.4330127&nbsp;&nbsp; 0.0000000), wk =&nbsp;&nbsp; 0.1875000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 20) = (&nbsp;&nbsp; 0.2500000&nbsp;&nbsp; 0.4330127&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 21) = (&nbsp;&nbsp; 0.2500000&nbsp;&nbsp; 0.4330127&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.3750000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 22) = (&nbsp;&nbsp; 0.2500000&nbsp;&nbsp; 0.4330127&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 23) = (&nbsp;&nbsp; 0.2500000&nbsp;&nbsp; 0.4330127&nbsp; -0.3048780), wk =&nbsp;&nbsp; 0.1875000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 24) = (&nbsp;&nbsp; 0.2500000&nbsp;&nbsp; 0.4330127&nbsp; -0.6097561), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 25) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.2886751&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.0937500<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 26) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.2886751&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 27) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.2886751&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.0937500<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 28) = (&nbsp;&nbsp; 0.0000000&nbsp; -0.2886751&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 29) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp;&nbsp; 0.1524390), wk =&nbsp;&nbsp; 0.0937500<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k(&nbsp;&nbsp; 30) = (&nbsp;&nbsp; 0.0000000&nbsp;&nbsp; 0.2886751&nbsp; -0.1524390), wk =&nbsp;&nbsp; 0.0000000<br><br>These are the input files for this simple test that I can be easily reproduces.<br>
<br><a href="http://sic.scf.fit.in">sic.scf.fit.in</a><br>---------------------------<br>&nbsp;&nbsp; &amp;control<br>&nbsp;&nbsp;&nbsp; calculation = &#39;scf&#39;<br>&nbsp;&nbsp;&nbsp; tstress = .true.<br>&nbsp;&nbsp;&nbsp; tprnfor = .true.<br>&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;./&#39;,<br>
&nbsp;&nbsp;&nbsp; prefix=&#39;sic2&#39;,<br>&nbsp;&nbsp;&nbsp; outdir=&#39;./&#39;<br>&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0d-5,<br>&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav=&nbsp; 4, celldm(1) = 5.81, celldm(3) = 1.64,<br>&nbsp;&nbsp;&nbsp; nat=&nbsp; 4, ntyp=3,<br>&nbsp;&nbsp;&nbsp; ecutwfc =25.0,ecutrho=200.0,<br>
&nbsp;&nbsp;&nbsp; occupations=&#39;smearing&#39;, smearing=&#39;methfessel-paxton&#39;, degauss=0.005,<br>&nbsp;&nbsp;&nbsp; la2F=.true.<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; diagonalization=&#39;david&#39;<br>&nbsp;&nbsp;&nbsp; mixing_mode = &#39;plain&#39;<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<br>
&nbsp;&nbsp;&nbsp; conv_thr&nbsp;&nbsp; = 1.0d-8<br>/<br>ATOMIC_SPECIES<br>C&nbsp;&nbsp; 12.011 C.pbe-rrkjus.UPF<br>Si&nbsp; 28.086 Si.pbe-n-van.UPF<br>B&nbsp;&nbsp; 10.811 B.pbe-n-van.UPF<br>ATOMIC_POSITIONS crystal<br>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.500000000<br>
Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.375000000<br>Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.875000000<br>B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000<br>K_POINTS automatic<br>&nbsp;&nbsp; 8 8 8 0 0 0<br><br><a href="http://sic.scf.in">sic.scf.in</a><br>
----------------<br>&amp;control<br>&nbsp;&nbsp;&nbsp; calculation = &#39;scf&#39;<br>&nbsp;&nbsp;&nbsp; tstress = .true.<br>&nbsp;&nbsp;&nbsp; tprnfor = .true.<br>&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;./&#39;,<br>&nbsp;&nbsp;&nbsp; prefix=&#39;sic2&#39;,<br>&nbsp;&nbsp;&nbsp; outdir=&#39;./&#39;<br>&nbsp;&nbsp;&nbsp; etot_conv_thr = 1.0d-5,<br>
&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav=&nbsp; 4, celldm(1) = 5.81, celldm(3) = 1.64,<br>&nbsp;&nbsp;&nbsp; nat=&nbsp; 4, ntyp=3,<br>&nbsp;&nbsp;&nbsp; ecutwfc =25.0,ecutrho=200.0,<br>&nbsp;&nbsp;&nbsp; occupations=&#39;smearing&#39;, smearing=&#39;methfessel-paxton&#39;, degauss=0.005,<br>
&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; diagonalization=&#39;david&#39;<br>&nbsp;&nbsp;&nbsp; mixing_mode = &#39;plain&#39;<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<br>&nbsp;&nbsp;&nbsp; conv_thr&nbsp;&nbsp; = 1.0d-8<br>/<br>ATOMIC_SPECIES<br>C&nbsp;&nbsp; 12.011 C.pbe-rrkjus.UPF<br>Si&nbsp; 28.086 Si.pbe-n-van.UPF<br>
B&nbsp;&nbsp; 10.811 B.pbe-n-van.UPF<br>ATOMIC_POSITIONS crystal<br>C&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.500000000<br>Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.375000000<br>Si&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.666666667&nbsp;&nbsp; 0.333333333&nbsp;&nbsp; 0.875000000<br>B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000<br>
K_POINTS automatic<br>&nbsp;&nbsp; 4 4 4 0 0 0<br><br><a href="http://sic.elph.in">sic.elph.in</a><br>----------------<br>Electron-phonon coefficients<br>&nbsp;&amp;inputph<br>&nbsp; tr2_ph=1.0d-10,<br>&nbsp; prefix=&#39;sic2&#39;,<br>&nbsp; fildvscf=&#39;dv&#39;,<br>
&nbsp; amass(1)=12.011,<br>&nbsp; amass(2)=28.086,<br>&nbsp; amass(3)=10.811,<br>&nbsp; outdir=&#39;./&#39;,<br>&nbsp; fildyn=&#39;sic2.dyn&#39;,<br>&nbsp; elph=.true.,<br>&nbsp; trans=.true.,<br>&nbsp; ldisp=.true.,<br>&nbsp; nq1=2, nq2=2,nq3=2<br>&nbsp;/<br><br>Thank you,<br>
&nbsp; Roxana Margine<br><br>====================================================<br>Elena Roxana Margine, Ph.D.<br>Laboratoire de Physique de la Matière Condensée et Nanostructures<br>Université Claude Bernard Lyon 1, 43 bd du 11 novembre 1918<br>
69622 Villeurbanne - France&nbsp;&nbsp;&nbsp; E-mail: <a href="mailto:rmargine@lpmcn.univ-lyon1.fr">rmargine@lpmcn.univ-lyon1.fr</a><br>Phone: +33 (0) 4 724 482 37&nbsp; <br>====================================================<br><br><br><br>
<br><br></div>