<br><blockquote><ol><li><div>Using a k-points grid that includes the K point can be a source of troubles, I would recommend using 16x16x1 or even 32x32x1 for better results.</div></li></ol></blockquote><div>p.s. and treat it as a metal, adding some smearing, to avoid the absurd dielectric constant</div><br><div><signature id="local-80497166-1ed1">--<br><span style="font-size:0.9em;color:gray;">Lorenzo Paulatto - Paris<span></span></span></signature></div><div class="gmail_quote_attribution">On Mar 30 2021, at 5:00 pm, Lorenzo Paulatto <paulatz@gmail.com> wrote:</div><blockquote><div>You have specified the q-points in crystal fractionary units, hence you have to specify <font style="font-family:monospace"><span style="color:rgb(24, 24, 178)">q_in_cryst_coord=.true.</span></font></div><div>This can be safely used in conjunction with q_in_band_form, which is a completely unrelated option.</div><br><div>Also:</div><ol><li><div>You way of specifying the cell parameters is discouraged and may cause problems with symmetry detection, but it should work.</div></li><li><br><blockquote><ol><li><div>Using a k-points grid that includes the K point can be a source of troubles, I would recommend using 16x16x1 or even 32x32x1 for better results.</div></li></ol></blockquote></li></ol><br><div>hth</div><br><div><div>--</div><div><span style="color:gray"><font style="font-size:0.9em">Lorenzo Paulatto - Paris</font></span></div></div><div class="gmail_quote_attribution">On Mar 30 2021, at 4:19 pm, Evgenii <eisysoev4@stud.eltech.ru> wrote:</div><blockquote><div><div><div>Dear Quantum Espresso experts, </div><div> </div><div>I am having trouble calculating the phonon bands of graphene.</div><div> </div><div>First, I took the cif-file of the already relaxed graphite cell from the Materials Project site https://next-gen.materialsproject.org/ and made graphene cell from it by hand.</div><div>Then I relaxed the cell obtained with vc-relax and 12x12x1 k-point set.</div><div>I used the same parameters to perform a scf calculation:</div><div> </div><div>&CONTROL</div><div>calculation='scf'</div><div>title='graphene'</div><div>prefix='graphene'</div><div>verbosity='high'</div><div>restart_mode='from_scratch'</div><div>nstep=1000</div><div>iprint=1</div><div>tstress=.true.</div><div>tprnfor=.true.</div><div>disk_io='default'</div><div>pseudo_dir = '/home/mxm2/Desktop/Graphene',</div><div>outdir='/home/mxm2/Desktop/Graphene/temp'</div><div>/</div><div>&SYSTEM</div><div>ibrav = 0,</div><div>celldm(1)=1.889726,</div><div>nat = 2,</div><div>ntyp = 1,</div><div>ecutwfc = 60.0 ,</div><div>ecutrho = 480.0 ,</div><div>/</div><div>&ELECTRONS</div><div>electron_maxstep = 1000,</div><div>conv_thr = 1.0d-12 ,</div><div>mixing_mode = 'plain' ,</div><div>mixing_beta = 0.7d0 ,</div><div>/</div><div>ATOMIC_SPECIES</div><div>C 12.011 C.pz-rrkjus.UPF</div><div>ATOMIC_POSITIONS (crystal)</div><div>C        0.000000000   0.000000000   0.000861580</div><div>C        0.333333000   0.666667000   0.000838420</div><div>CELL_PARAMETERS (alat)</div><div>   2.438515837  -0.000000000   0.000000000</div><div>  -1.219257919   2.111816662   0.000000000</div><div>  -0.000000000  -0.000000000  19.995916314</div><div> </div><div>K_POINTS automatic</div><div>12 12 1 0 0 0</div><div> </div><div>Then, I run ph.in calculation:</div><div> </div><div>&inputph</div><div>recover=.true.</div><div>tr2_ph=1.0d-14,  </div><div>prefix='graphene',</div><div>ldisp=.true. </div><div>nq1=12, </div><div>nq2=12, </div><div>nq3=1,  </div><div>amass(1)=12.011,   </div><div>outdir='/home/mxm2/Desktop/Graphene/temp', </div><div>fildyn='graphene.dyn',    </div><div>/</div><div> </div><div>My q2r:</div><div> </div><div>&input</div><div>fildyn='graphene.dyn',</div><div>zasr='simple',</div><div>flfrc='graphene.12x12x1.fc'</div><div>/</div><div> </div><div>Matdyn input:</div><div> </div><div>&input</div><div>asr='simple',</div><div>amass(1)=12.011,</div><div>flfrc='graphene.12x12x1.fc', </div><div>flfrq='graphene.disp.freq' </div><div>q_in_band_form=.true.,</div><div>/</div><div>4 ! number of q-points</div><div>0.000 0.000 0.000 10 ! gamma point</div><div>0.500 0.000 0.000 10 ! M</div><div>0.333 0.333 0.000 10 ! K</div><div>0.000 0.000 0.000  1 ! gamma point</div><div> </div><div>Plotband:</div><div> </div><div>graphene.disp.freq</div><div>0 1600</div><div>graphene.disp.freq.dat</div><div>graphene.disp.freq.ps</div><div>0</div><div>100 0</div><div> </div><div>Why the phonon dispersion does not match the pictures from the papers and for calculations with e.g. 18x18x1 k-point and nq-point sets the output file of graphene.dyn1 look like:</div><div> </div><div>Dielectric Tensor:</div><div> </div><div>************************        166.332153320312          0.000000000000</div><div>        166.332275390625************************          0.000000000000</div><div>          0.000000000000          0.000000000000          1.116060253916</div><div> </div><div>     Effective Charges E-U: Z_{alpha}{s,beta}</div><div> </div><div>     atom #    1</div><div>  -11117471.575740195811         -0.002632206306          0.000000000000</div><div>         -0.002632205375-11117471.579664016142            0.000000000000</div><div>          0.000000000000          0.000000000000         -0.003734609171</div><div>          </div><div><div>Also, what is the best way to set the lattice parameters, via Ibrav or CELL_PARAMETERS?</div></div><div>Any help is appreciated.</div><div> </div><div>Best regards,</div><div>Evgeniy Sysoev.</div></div></div><div>_______________________________________________</div><div>Quantum ESPRESSO is supported by MaX (www.max-centre.eu)</div><div>users mailing list users@lists.quantum-espresso.org</div><div>https://lists.quantum-espresso.org/mailman/listinfo/users</div></blockquote></blockquote><img class="mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border:0; width:0; height:0;" src="https://link.getmailspring.com/open/9A3DB2B1-03BC-41F6-8DA9-32A82466AF83@getmailspring.com?me=ca1a67e4&recipient=dXNlcnNAbGlzdHMucXVhbnR1bS1lc3ByZXNzby5vcmc%3D">