<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p></p>
<pre>Thank you for your reply Dr. Lorenzo.<br><br>As you say, my system is isolated at every direction, <br>just a single finite-sized nanoribbon constructed along the x direction (with a definite y-width).<br><br>Then, should i use just a single k-point (gamma-point) for a 4x4x1 supercell, in a band calculation?<br><br>And again, in a Thermal conductivity calculation, which requires 3 and 4 order forces<br>should I use this same input?<br><br>Thank you again for your time and patience.<br><br>Pablo.<br><br>On Wednesday, 26 April 2017 14:26:54 CEST Pablo Álvarez Zapatero wrote:
><i> So the question is, is this input fine? Do will solve the energy problem
</i>><i> (DFT ) accurately?
</i>><i>
</i>
Dear Pablo,
if your system is isolated (i.e. not periodic) along a certain direction, than
you should use only one k-point in that direction. If the results are absurd
the problem is not the lack of k-points.
If you are building a supercell, then you can safely divide the numbr of k-
points in the unit cell by the number of unit cells in the supercell. This is
perfectly accurate.
However, I do not understand really what you are trying to do. Phonopy
computes 2-body force constants by finite differences, it does not "solve"
anything. It then uses these force constants to interpolate the phonon
dispersion. In order to have decent phonons in Graphene, you need to compute
2-body force constants up to at the very least four unit cells, i.e. using a
4x4x1 supercell.
Incidentally this can be done with the ph.x code at a fraction of the
computational cost.
cheers
--
Dr. Lorenzo Paulatto </pre>
<br>
<p></p>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr style="display:inline-block; width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>De:</b> Pablo Álvarez Zapatero <pabloalvazapa@hotmail.com><br>
<b>Enviado:</b> miércoles, 26 de abril de 2017 14:26<br>
<b>Para:</b> pw_forum@pwscf.org<br>
<b>Asunto:</b> About graphene nanoribbon input.</font>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hi everyone,</p>
<p><br>
</p>
<p>I am trying to calculate thermal conductivity of a graphene nanoribbon throught other program called "phono3py", which uses QE to perform ab initio calculations. Then, in order to understand and validating my results, i am currently trying to calculate the
phonon bands of a different sized graphene nanoribbon, using the so called phonopy.</p>
<p>Firstly, i succesfully solved the infinite (phonon bands) grapehene system with phonopy+QE. Now, I have looked for propper inputs in QE forums for a nanorribon, but i do not know if the input i am using is well posed. Using a single gamma k point (as i read
in this forum), results are ridicolous. Now, using these k-points it lasts for 24 hours for a sigle .in calculation (there are 11 more in a band calculation).</p>
<p><br>
</p>
<p></p>
<div>&control<br>
calculation = 'scf'<br>
prefix='grapheneph',<br>
tprnfor = .true.<br>
tstress = .true.<br>
pseudo_dir = './',<br>
outdir='./saleph51/'<br>
/<br>
&system<br>
ibrav = 0,<br>
nat = 84,<br>
ntyp = 1,<br>
occupations = 'smearing'<br>
smearing = 'm-p'<br>
degauss = 0.01<br>
ecutwfc = 30.0,<br>
ecutrho = 150.0,<br>
/<br>
&electrons<br>
mixing_beta= 0.7<br>
conv_thr = 1.0d-8<br>
diagonalization = 'david'<br>
mixing_mode = 'plain'<br>
/<br>
<b>K_POINTS automatic<br>
12 12 1 0 0 0</b><br>
! ibrav = 0, nat = 84, ntyp = 1<br>
CELL_PARAMETERS bohr<br>
29.8200000000000003 0.0000000000000000 0.0000000000000000<br>
0.0000000000000000 16.1490000000000009 0.0000000000000000<br>
0.0000000000000000 0.0000000000000000 10.0000000000000000<br>
ATOMIC_SPECIES<br>
C 12.01070 C.pbe-rrkjus.UPF</div>
<span>ATOMIC_POSITIONS crystal</span>
<p></p>
<p>[...]<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p>I read that is convenient to use a multiple of 3 k-point mesh to achieve better results in any direction (i read for instance 12 or even 60!). Since my system is flat, i set that k-point mesh in x-y.
<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p>So the question is, is this input fine? Do will solve the energy problem (DFT ) accurately?</p>
<p><br>
</p>
<p>Thank you for your time and patience!</p>
<p>Pablo.<br>
</p>
</div>
</div>
</div>
</div>
</body>
</html>