<div dir="ltr"><div>Answering myself: later, for a specific k-point. Quick example showing the problem attached.<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 8, 2016 at 7:12 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"><div dir="ltr"><div>When does this happen? at the beginning of the phonon calculation or later?<br><br></div>Paolo<br></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Tue, Nov 8, 2016 at 1:22 PM, Antoine JAY <span dir="ltr"><<a href="mailto:antoine.jay@polytechnique.edu" target="_blank">antoine.jay@polytechnique.edu</a><wbr>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello community & developers,<br>
<br>
Performing phonons on gamma-boron at 30 GPa (using GGA), I have this message:<br>
  Â  Â from set_irr_sym_new : error #  Â  Â  2111<br>
  Â  Â wrong representation<br>
<br>
I send down the atomic positions as I know from previous mail (Lorenzo<br>
Paulatto 2012-07-23) that you need some tests to find a reliable way to<br>
generate the initial random pseudo-dynamical matrix (random_matrix.f90).<br>
Lattice and atomic positions are already converged.<br>
<br>
Increasing the threshold 10^-6 to 10^-4 line 156 and 158  of set_irr_sym.f90<br>
does not solve this issue.<br>
<br>
Do you have any other suggestion?<br>
<br>
Thank you very much<br>
<br>
Antoine Jay<br>
<br>
&control<br>
  Â  calculation = 'scf',<br>
  Â  prefix  Â  Â  = '$name.$Pression',<br>
  Â  tstress  Â  Â = .true.,<br>
  Â  tprnfor  Â  Â = .true.,<br>
  Â  pseudo_dir  = '$PSEUDO_DIR/',<br>
  Â  outdir  Â  Â  = '$TMP_DIR/',<br>
  Â  etot_conv_thr=1.0d-6,<br>
  Â  forc_conv_thr=1.0d-4,<br>
  Â  nstep=1400,<br>
  Â  wf_collect=.true.<br>
/<br>
&system<br>
  Â  ibrav =8<br>
  Â  celldm(1) = 9.213227<br>
  Â  celldm(2) = 1.113897<br>
  Â  celldm(3) = 1.365520<br>
  Â  nat=  28,<br>
  Â  ntyp= 5,<br>
  Â  ecutwfc =80,<br>
/<br>
 &electrons<br>
  Â  diagonalization='david'<br>
  Â  mixing_beta = 0.7<br>
  Â  conv_thr =  1.0d-13<br>
  Â  electron_maxstep=800<br>
/<br>
ATOMIC_SPECIES<br>
B1  Â 10.810  B.gga.pw91.UPF<br>
B2  Â 10.810  B.gga.pw91.UPF<br>
B3  Â 10.810  B.gga.pw91.UPF<br>
B4  Â 10.810  B.gga.pw91.UPF<br>
B5  Â 10.810  B.gga.pw91.UPF<br>
ATOMIC_POSITIONS {crystal}<br>
B1  Â  Â  Â 0.1704979  Â 0.475506132  Â 0.000000000<br>
B2  Â  Â  Â 0.1597472  Â 0.719824335  Â 0.627390639<br>
B3  Â  Â  Â 0.3462299  Â 0.908257987  Â 0.787690908<br>
B4  Â  Â  Â 0.3505384  Â 0.726605058  Â 0.000000000<br>
B5  Â  Â  Â 0.1635420  Â 0.991258129  Â 0.000000000<br>
B1  Â  Â  Â 0.8295020  Â 0.524493868  Â 0.000000000<br>
B1  Â  Â  Â 0.3295030  Â 0.975507731  Â 0.500000000<br>
B1  Â  Â  Â 0.6704969  Â 0.024492269  Â 0.500000000<br>
B2  Â  Â  Â 0.8402527  Â 0.280175665  Â 0.627390639<br>
B2  Â  Â  Â 0.3402485  Â 0.219820769  Â 0.872609128<br>
B2  Â  Â  Â 0.6597514  Â 0.780179231  Â 0.872609128<br>
B2  Â  Â  Â 0.8402527  Â 0.280175665  Â 0.372609361<br>
B2  Â  Â  Â 0.1597472  Â 0.719824335  Â 0.372609361<br>
B2  Â  Â  Â 0.6597514  Â 0.780179231  Â 0.127390872<br>
B2  Â  Â  Â 0.3402485  Â 0.219820769  Â 0.127390872<br>
B3  Â  Â  Â 0.6537700  Â 0.091742013  Â 0.787690908<br>
B3  Â  Â  Â 0.1537717  Â 0.408261072  Â 0.712309649<br>
B3  Â  Â  Â 0.8462282  Â 0.591738928  Â 0.712309649<br>
B3  Â  Â  Â 0.6537700  Â 0.091742013  Â 0.212309092<br>
B3  Â  Â  Â 0.3462299  Â 0.908257987  Â 0.212309092<br>
B3  Â  Â  Â 0.8462282  Â 0.591738928  Â 0.287690351<br>
B3  Â  Â  Â 0.1537717  Â 0.408261072  Â 0.287690351<br>
B4  Â  Â  Â 0.6494615  Â 0.273394942  Â 0.000000000<br>
B4  Â  Â  Â 0.1494677  Â 0.226608664  Â 0.500000000<br>
B4  Â  Â  Â 0.8505322  Â 0.773391336  Â 0.500000000<br>
B5  Â  Â  Â 0.8364579  Â 0.008741871  Â 0.000000000<br>
B5  Â  Â  Â 0.3364577  Â 0.491264003  Â 0.500000000<br>
B5  Â  Â  Â 0.6635422  Â 0.508735997  Â 0.500000000<br>
K_POINTS {automatic}<br>
6 6 6 1 1 1<br>
EOF<br>
<br>
and the ph run:<br>
<br>
phonon of $name<br>
 &inputph<br>
  recover=.true.<br>
  tr2_ph=1.0d-14,<br>
  prefix='$name.$Pression',<br>
  amass(1)=10.810,<br>
  amass(2)=10.810,<br>
  amass(3)=10.810,<br>
  amass(4)=10.810,<br>
  amass(5)=10.810,<br>
  outdir='$TMP_DIR/'<br>
  fildyn='$name.$Pression.dyn',<br>
  fildrho='$name.$Pression.drho'<wbr>,<br>
  ldisp=.true.,<br>
  nq1=4, nq2=4, nq3=4<br>
 /<br>
EOF<br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listi<wbr>nfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br></div></div><span class="HOEnZb"><font color="#888888">-- <br><div class="m_350171698652569959gmail_signature" data-smartmail="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 <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a><br><br></div></div></div></div></div>
</font></span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="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>
</div>