<div dir="ltr"><div><div>Dear PW community, <br><br></div>I have the following error <br><br>     Error in routine set_irr_sym_new (511):<br>     wrong representation<br><br></div><div>That I cannot get rid off while trying to do a phonon calculation.<br></div><div><br>My PH input<br></div><div>--<br>&inputph<br>  recover  = .true.,<br>  prefix   = 'x',<br>  fildyn   = 'x.dyn',<br>  ldisp    = .true.,<br>  epsil    = .true.,<br>  trans    = .true.,<br>  fildvscf = 'dvscf',<br>  niter_ph = 800,<br>  nq1=8,<br>  nq2=8,<br>  nq3=2,<br>  tr2_ph   =  1.0d-14,<br>  max_seconds = 86000,<br> /<br><br></div><div>Prior to that I made a vc-relax (cell and ions) calculation with<br>    tprnfor         = .true.,<br>    tstress         = .true.,<br>    forc_conv_thr   = 1d-5, ! Ry/Bohr<br></div><div>which converged (with the correct sym). I then used these lattice paramters for the ph calculation (through a first scf run).<br></div><div><br></div><div>The q-point that crashes is the following one:<br>     Calculation of q =    0.0000000  -0.5773503  -0.1312539<br><br></div><div>The q-points before that one are fine. <br><br></div><div>Things I've checked so far:<br><br></div><div>1) change the q-grid to 8x8x3. The problematic q-point is then <br></div><div>q =    0.0000000   0.1443376   0.0875026<br></div><div>In that case the problem is not the same (no error) but the calculation never converged<br>      iter # 361 total cpu time : 81102.7 secs   <a href="http://av.it">av.it</a>.:  16.5<br>      thresh= 1.419E-03 alpha_mix =  0.700 |ddv_scf|^2 =  3.809E-04<br><br>      iter # 362 total cpu time : 81178.6 secs   <a href="http://av.it">av.it</a>.:  13.8<br>      thresh= 1.952E-03 alpha_mix =  0.700 |ddv_scf|^2 =  9.521E-05<br>....<br><br><br></div><div>2) I've increased ecutwfc from 90 to 120 in the scf run<br></div><div>same error<br></div><div><br></div><div>3) I've increase kgrid from 10x10x3 to 16x16x4<br></div><div>same error<br><br></div><div>4) I went into /PHonon/PH/set_irr_sym.f90<br></div><div>And changed <br></div><div>              if (jpert.ne.ipert .and. abs(wrk)> 1.d-6 ) &<br>                     call errore('set_irr_sym_new','wrong representation',100*irr+10*jpert+ipert)<br>              if (jpert.eq.ipert .and. abs(wrk-1.d0)> 1.d-6 ) &<br>                     call errore('set_irr_sym_new','wrong representation',100*irr+10*jpert+ipert)<br></div><div>into<br>              if (jpert.ne.ipert .and. abs(wrk)> 1.d-3 ) &<br>                     call errore('set_irr_sym_new','wrong representation',100*irr+10*jpert+ipert)<br>              if (jpert.eq.ipert .and. abs(wrk-1.d0)> 1.d-3 ) &<br>                     call errore('set_irr_sym_new','wrong representation',100*irr+10*jpert+ipert)<br></div><div>Same error<br><br><br></div><div>What else can I test?<br><br></div><div>PS: I'm using the latest trunk of QE. <br><br></div><div>Cheers, <br><br></div><div>Samuel Ponce<br></div><div><br><br><br><br><br><br><br><br><br><br></div></div>