<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
The coordinates of the high-symmetry points in matdyn.in should be in cartesian not crystal because you are specifying q_in_band_form=.true.
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Vahid</div>
<div class=""><br class="">
</div>
<div class="">Vahid Askarpour<br class="">
Department of physics and atmospheric science<br class="">
Dalhousie University<br class="">
Halifax, NS<br class="">
CANADA<br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Mar 30, 2021, at 11:19 AM, Evgenii <<a href="mailto:eisysoev4@stud.eltech.ru" class="">eisysoev4@stud.eltech.ru</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class=""><!-- START CAUTION Box Code -->
<table align="center" border="0" cellpadding="0" cellspacing="0" style="padding:10px 0 10px 0" width="100%" class="">
<tbody class="">
<tr class="">
<td style="line-height:0px;font-size:0px;mso-line-height-rule:exactly;" class="">
<table align="center" border="0" cellpadding="0" cellspacing="0" style="background:#707372;background-color:#707372;width:100%;border-radius:5px;overflow:hidden;" class="">
<tbody class="">
<tr class="">
<td style="border-top:solid 8px #fbe122;padding:4px 8px;text-align:left;vertical-align:top;" class="">
<table role="presentation" border="0" cellpadding="0" cellspacing="0" class="">
<tbody class="">
<tr class="">
<td align="left" class="">
<div style="font-family:Arial,sans-serif;font-size:12px;line-height:16px;text-align:left;color:#ffffff;" class="">
<span style="font-weight:bold;font-size:12px;" class="">CAUTION:</span> The Sender of this email is not from within Dalhousie.</div>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
<!-- END CAUTION Box Code -->
<div class="">
<div class="">
<div class="">
<div class="">Dear Quantum Espresso experts, </div>
<div class=""> </div>
<div class="">I am having trouble calculating the phonon bands of graphene.</div>
<div class=""> </div>
<div class="">First, I took the cif-file of the already relaxed graphite cell from the Materials Project site
<a href="https://next-gen.materialsproject.org/" class="">https://next-gen.materialsproject.org/</a> and made graphene cell from it by hand.</div>
<div class="">Then I relaxed the cell obtained with vc-relax and 12x12x1 k-point set.</div>
<div class="">I used the same parameters to perform a scf calculation:</div>
<div class=""> </div>
<div class="">&CONTROL</div>
<div class="">calculation='scf'</div>
<div class="">title='graphene'</div>
<div class="">prefix='graphene'</div>
<div class="">verbosity='high'</div>
<div class="">restart_mode='from_scratch'</div>
<div class="">nstep=1000</div>
<div class="">iprint=1</div>
<div class="">tstress=.true.</div>
<div class="">tprnfor=.true.</div>
<div class="">disk_io='default'</div>
<div class="">pseudo_dir = '/home/mxm2/Desktop/Graphene',</div>
<div class="">outdir='/home/mxm2/Desktop/Graphene/temp'</div>
<div class="">/</div>
<div class="">&SYSTEM</div>
<div class="">ibrav = 0,</div>
<div class="">celldm(1)=1.889726,</div>
<div class="">nat = 2,</div>
<div class="">ntyp = 1,</div>
<div class="">ecutwfc = 60.0 ,</div>
<div class="">ecutrho = 480.0 ,</div>
<div class="">/</div>
<div class="">&ELECTRONS</div>
<div class="">electron_maxstep = 1000,</div>
<div class="">conv_thr = 1.0d-12 ,</div>
<div class="">mixing_mode = 'plain' ,</div>
<div class="">mixing_beta = 0.7d0 ,</div>
<div class="">/</div>
<div class="">ATOMIC_SPECIES</div>
<div class="">C 12.011 C.pz-rrkjus.UPF</div>
<div class="">ATOMIC_POSITIONS (crystal)</div>
<div class="">C 0.000000000 0.000000000 0.000861580</div>
<div class="">C 0.333333000 0.666667000 0.000838420</div>
<div class="">CELL_PARAMETERS (alat)</div>
<div class=""> 2.438515837 -0.000000000 0.000000000</div>
<div class=""> -1.219257919 2.111816662 0.000000000</div>
<div class=""> -0.000000000 -0.000000000 19.995916314</div>
<div class=""> </div>
<div class="">K_POINTS automatic</div>
<div class="">12 12 1 0 0 0</div>
<div class=""> </div>
<div class="">Then, I run ph.in calculation:</div>
<div class=""> </div>
<div class="">&inputph</div>
<div class="">recover=.true.</div>
<div class="">tr2_ph=1.0d-14, </div>
<div class="">prefix='graphene',</div>
<div class="">ldisp=.true. </div>
<div class="">nq1=12, </div>
<div class="">nq2=12, </div>
<div class="">nq3=1, </div>
<div class="">amass(1)=12.011, </div>
<div class="">outdir='/home/mxm2/Desktop/Graphene/temp', </div>
<div class="">fildyn='graphene.dyn', </div>
<div class="">/</div>
<div class=""> </div>
<div class="">My q2r:</div>
<div class=""> </div>
<div class="">&input</div>
<div class="">fildyn='graphene.dyn',</div>
<div class="">zasr='simple',</div>
<div class="">flfrc='graphene.12x12x1.fc'</div>
<div class="">/</div>
<div class=""> </div>
<div class="">Matdyn input:</div>
<div class=""> </div>
<div class="">&input</div>
<div class="">asr='simple',</div>
<div class="">amass(1)=12.011,</div>
<div class="">flfrc='graphene.12x12x1.fc', </div>
<div class="">flfrq='graphene.disp.freq' </div>
<div class="">q_in_band_form=.true.,</div>
<div class="">/</div>
<div class="">4 ! number of q-points</div>
<div class="">0.000 0.000 0.000 10 ! gamma point</div>
<div class="">0.500 0.000 0.000 10 ! M</div>
<div class="">0.333 0.333 0.000 10 ! K</div>
<div class="">0.000 0.000 0.000 1 ! gamma point</div>
<div class=""> </div>
<div class="">Plotband:</div>
<div class=""> </div>
<div class="">graphene.disp.freq</div>
<div class="">0 1600</div>
<div class="">graphene.disp.freq.dat</div>
<div class="">graphene.disp.freq.ps</div>
<div class="">0</div>
<div class="">100 0</div>
<div class=""> </div>
<div class="">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 class=""> </div>
<div class="">Dielectric Tensor:</div>
<div class=""> </div>
<div class="">************************ 166.332153320312 0.000000000000</div>
<div class=""> 166.332275390625************************ 0.000000000000</div>
<div class=""> 0.000000000000 0.000000000000 1.116060253916</div>
<div class=""> </div>
<div class=""> Effective Charges E-U: Z_{alpha}{s,beta}</div>
<div class=""> </div>
<div class=""> atom # 1</div>
<div class=""> -11117471.575740195811 -0.002632206306 0.000000000000</div>
<div class=""> -0.002632205375-11117471.579664016142 0.000000000000</div>
<div class=""> 0.000000000000 0.000000000000 -0.003734609171</div>
<div class=""> </div>
<div class="">
<div class="">Also, what is the best way to set the lattice parameters, via Ibrav or CELL_PARAMETERS?</div>
</div>
<div class="">Any help is appreciated.</div>
<div class=""> </div>
<div class="">Best regards,</div>
<div class="">Evgeniy Sysoev.</div>
</div>
</div>
</div>
</div>
<span id="cid:0B30C32A-4939-45E7-9E9F-5CE3D4B64253"><graphene.disp.freq.jpg></span>_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" class="">www.max-centre.eu</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>