<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="">
Dear Evgeniy,
<div class=""><br class="">
</div>
<div class="">I use the following for my graphene dispersion. Please note that the actual coordinates depend on the PSP.</div>
<div class=""><br class="">
</div>
<div class="">4<br class="">
0.0 0.0 0.0 150 ! G<br class="">
0.0000000 -0.577588 0.0 150 ! M<br class="">
-0.3334706 -0.577588 0.0 150 ! K<br class="">
0.0 0.0 0.0 1</div>
<div class=""><br class="">
</div>
<div class="">So the matdyn.x code will automatically set 150 points between the first 2 high-symmetry points and so on.</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Vahid<br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Mar 30, 2021, at 11:58 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="">Dear Vahid,</div>
<div class=""> </div>
<div class="">
<div class="">If we use the default values of the <span style="background-color: rgb(255, 255, 255); float: none; font-family: arial, sans-serif; font-size: 15px; font-style: normal; font-weight: 400; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;" class="">
q_in_band_form</span> variables, then how do we find intermediate points between the points of high symmetry, as for example on this site
<a href="https://www.afs.enea.it/buonocor/graphane-with-phonon.html" rel="noopener noreferrer" target="_blank" class="">
https://www.afs.enea.it/buonocor/graphane-with-phonon.html</a> ? Are intermediate points always needed, we noticed that in some examples they are missing.</div>
<div class=""> </div>
<div class="">
<div class="">Thanks for the help.</div>
<div class=""> </div>
<div class="">Evgeniy Sysoev</div>
</div>
<div class=""> </div>
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">30.03.2021, 17:28, "Vahid Askarpour" <<a href="mailto:vh261281@dal.ca" class="">vh261281@dal.ca</a>>:</div>
<blockquote class="">
<div class="" style="word-wrap:break-word">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 class=""><br class="">
<blockquote class="">
<div class="">On Mar 30, 2021, at 11:19 AM, Evgenii <<a class="" href="mailto:eisysoev4@stud.eltech.ru">eisysoev4@stud.eltech.ru</a>> wrote:</div>
<br class="34bd394265c9c48ce7af3175baf9655dApple-interchange-newline">
<div class="">
<div class="">
<table class="" align="center" border="0" cellpadding="0" cellspacing="0" width="100%" style="padding:10px 0 10px 0">
<tbody class="">
<tr class="">
<td class="" style="font-size:0px;line-height:0px">
<table class="" align="center" border="0" cellpadding="0" cellspacing="0" style="background:#707372;border-radius:5px;overflow:hidden;width:100%">
<tbody class="">
<tr class="">
<td class="" style="border-top-color:#fbe122;border-top-style:solid;border-top-width:8px;padding:4px 8px 4px 8px;text-align:left;vertical-align:top">
<table class="" border="0" cellpadding="0" cellspacing="0">
<tbody class="">
<tr class="">
<td class="" align="left">
<div class="" style="color:#ffffff;font-family:'arial' , sans-serif;font-size:12px;line-height:16px;text-align:left">
<span class="" style="font-size:12px;font-weight:bold">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>
<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 class="" href="https://next-gen.materialsproject.org/">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.<span class="1f1ea193f6735cf0wmi-callto">000000000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000861580</span></div>
<div class="">C 0.<span class="1f1ea193f6735cf0wmi-callto">333333000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">666667000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000838420</span></div>
<div class="">CELL_PARAMETERS (alat)</div>
<div class=""> 2.<span class="1f1ea193f6735cf0wmi-callto">438515837</span> -0.<span class="1f1ea193f6735cf0wmi-callto">000000000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000</span></div>
<div class=""> -1.<span class="1f1ea193f6735cf0wmi-callto">219257919</span> 2.<span class="1f1ea193f6735cf0wmi-callto">111816662</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000</span></div>
<div class=""> -0.<span class="1f1ea193f6735cf0wmi-callto">000000000</span> -0.<span class="1f1ea193f6735cf0wmi-callto">000000000</span> 19.<span class="1f1ea193f6735cf0wmi-callto">995916314</span></div>
<div class=""> </div>
<div class="">K_POINTS automatic</div>
<div class=""><span class="1f1ea193f6735cf0wmi-callto">12 12 1 0 0 0</span></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.<span class="1f1ea193f6735cf0wmi-callto">332153320312</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span></div>
<div class=""> 166.<span class="1f1ea193f6735cf0wmi-callto">332275390625</span>************************ 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span></div>
<div class=""> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span> 1.<span class="1f1ea193f6735cf0wmi-callto">116060253916</span></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=""> -<span class="1f1ea193f6735cf0wmi-callto">11117471</span>.<span class="1f1ea193f6735cf0wmi-callto">575740195811</span> -0.<span class="1f1ea193f6735cf0wmi-callto">002632206306</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span></div>
<div class=""> -0.<span class="1f1ea193f6735cf0wmi-callto">002632205375-11117471</span>.<span class="1f1ea193f6735cf0wmi-callto">579664016142</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span></div>
<div class=""> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span> 0.<span class="1f1ea193f6735cf0wmi-callto">000000000000</span> -0.<span class="1f1ea193f6735cf0wmi-callto">003734609171</span></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="5d6269b6be4d49e0cfc96bf333b4085cid:0B30C32A-4939-45E7-9E9F-5CE3D4B64253" class=""><graphene.disp.freq.jpg></span>_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a class="" href="http://www.max-centre.eu/">www.max-centre.eu</a>)<br class="">
users mailing list <a class="" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br class="">
<a class="" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
,
<p class="">_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="https://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><br class="">
</p>
</blockquote>
</div>
</div>
_______________________________________________<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>