<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 style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span id="" title="" class="elementToProof">Dear Lorenzo,</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span id="" title="" class=""><span id="" title="" class=""><br>
</span></span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span id="" title="" class=""><span id="" title="" class="">Thank you</span></span> very much much for your detailed explanation and your valuable time.</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
I will continue to perform some tests and if the problem of the of a negative line width persists, I will let you know.</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
I will be also cheking the store_lw flag, as I am using the SMA method for the time being.
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Regards <br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Lorenzo Paulatto <lorenzo.paulatto@cnrs.fr><br>
<b>Sent:</b> Friday, June 16, 2023 8:06 AM<br>
<b>To:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br>
<b>Subject:</b> Re: [QE-users] A question on the volume_factor flag in d3_tk.x</font>
<div> </div>
</div>
<div style="padding-bottom:1px">
<p><br>
</p>
<div class="x_moz-cite-prefix">On 6/15/23 21:12, Elio Physics wrote:<br>
</div>
<blockquote type="cite"><style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Dear professor Paulatto,</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Thanks for the information. <br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
1)Would that divergence produce negative line widths? I have been trying to obtain some decent results for graphene but when the smearing is above 1.0 , the code crashes stating that the linewidth is negative! i tried to increase the grid dimension but in vain.<br>
</div>
</blockquote>
<p><br>
</p>
<p>To get a negative linewidth it takes a seriously broken calculation, as the last therm can mathematically be negative</p>
<p><img alt="" data-outlook-trace="F:1|T:1" src="cid:part1.KFZ5q240.nA6oCBjt@cnrs.fr"></p>
<p><br>
</p>
<p>But that would require that the population of the higher phonon (q'') is higher than the lower phonon (q') as w'' = w+w' > w'.</p>
<p><br>
</p>
<p>Also a broken compilation of the code, or a bug, but to say anything I would need a specific case to test and reproduce.<br>
</p>
<p><br>
</p>
<blockquote type="cite">
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
2) Is there a flag in the input that ensures regularization or one that imposes a cutoff?</div>
</blockquote>
<p><br>
</p>
<p>No, it has been done ad-hoc for the specific paper about 2D systems, although I did not do those calculations, I'm not sure of the details. It is very easily done in post-processing, if store_lw=.true. is used to save the linewidths to file, using the matlab
script recompute_sma.m, or any home-made script. A bit more complicated in the exact minimization method, setting a reasonable Casimir scattering length may be the best bet. Exact algorithm is very expensive to converge anyway, because one cannot use shifted
grid (and symmetry is not exploited at the moment)<br>
</p>
<p><br>
</p>
<p><br>
</p>
<blockquote type="cite">
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div class="x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Thank you <br>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> users
<a class="x_moz-txt-link-rfc2396E" href="mailto:users-bounces@lists.quantum-espresso.org">
<users-bounces@lists.quantum-espresso.org></a> on behalf of Lorenzo Paulatto <a class="x_moz-txt-link-rfc2396E" href="mailto:lorenzo.paulatto@cnrs.fr">
<lorenzo.paulatto@cnrs.fr></a><br>
<b>Sent:</b> Thursday, June 15, 2023 5:02 PM<br>
<b>To:</b> Quantum ESPRESSO users Forum <a class="x_moz-txt-link-rfc2396E" href="mailto:users@lists.quantum-espresso.org">
<users@lists.quantum-espresso.org></a><br>
<b>Subject:</b> Re: [QE-users] A question on the volume_factor flag in d3_tk.x</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div dir="auto">Hello Elio,<br>
The problem is that tk is given for bulk per unit of volume. When you do a 2D material in plane waves, you need to put a lot of vacuum. The volume factor is used to make the value independent on the amount of vacuum. It can is the ratio of the bulk uni cell
to the 2d unit cell, or the ratio between the c axis in bulk and in the film, it is the same thing.<br>
<br>
Please not that 2d materials are very tricky: the thermal conductivity from the quadratic band is integrable but diverging. Using some kind of regularization or a cutoff may be required.<br>
<br>
Cheets<br>
<br>
</div>
<br>
<br>
<div class="x_x_gmail_quote">
<div dir="auto">On June 15, 2023 8:48:16 p.m. GMT+02:00, Elio Physics <a class="x_moz-txt-link-rfc2396E" href="mailto:Elio-Physics@live.com">
<Elio-Physics@live.com></a> wrote:</div>
<blockquote class="x_x_gmail_quote" style="margin:0pt 0pt 0pt
0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div class="x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<span class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Dear all (and specifically Professor Paulatto),</span>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br class="x_x_ContentPasted0">
</div>
<div class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Sorry to bother you again with another question regarding the d3q code. To calculate the lattice thermal conductivity for 2D materials, we need to define a volume_factor =H/V. H is defined as the interlayer distance for the corresponding bulk material. If I
understood correctly, let us say we have a 1H MoS2, we know that the bulk can be modeled as two layers (a total of six atoms per unit cell). H would then be the distance between the upper S layer of the bottom MoS2 and the lower S layer of the MoS2 on top.
Correct?</div>
<div class="x_x_x_elementToProof"><span><br class="x_x_ContentPasted0">
</span></div>
<div class="x_x_x_elementToProof"><span><br class="x_x_ContentPasted0">
</span></div>
<div class="x_x_x_elementToProof"><span class="x_x_ContentPasted0">Now what if we have a 1T 2D structure. Up to my knowledge, the bulk is modeled as one layer ( 3 atoms per unit cell). In this case, what would H be ? The thickness of the structure?</span></div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<div><br class="x_x_ContentPasted0">
</div>
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br class="x_x_ContentPasted0">
</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br class="x_x_ContentPasted0">
</div>
<div class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Regards</div>
<div class="x_x_x_elementToProof" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br class="x_x_ContentPasted0">
</div>
<div class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Elie Moujaes</div>
<div class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Federal University of Rondonia</div>
<div class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Brazil</div>
<span class="x_x_x_elementToProof x_x_ContentPasted0" style="font-family:Calibri,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">Porto Velho</span><br>
</div>
</blockquote>
</div>
</div>
<br>
<fieldset class="x_moz-mime-attachment-header"></fieldset>
<pre class="x_moz-quote-pre">_______________________________________________
The Quantum ESPRESSO community stands by the Ukrainian
people and expresses its concerns about the devastating
effects that the Russian military offensive has on their
country and on the free and peaceful scientific, cultural,
and economic cooperation amongst peoples
_______________________________________________
Quantum ESPRESSO is supported by MaX (<a class="x_moz-txt-link-abbreviated" href="http://www.max-centre.eu">www.max-centre.eu</a>)
users mailing list <a class="x_moz-txt-link-abbreviated" href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>
<a class="x_moz-txt-link-freetext" href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></pre>
</blockquote>
<div class="x_moz-signature">-- <br>
<small>Dr. Lorenzo Paulatto<br>
IdR @ IMPMC - CNRS UMR 7590 & Sorbonne Université<br>
phone: +33 (0)1 442 79822 / skype: paulatz<br>
<a href="http://www.impmc.upmc.fr/~paulatto/" class="x_moz-txt-link-freetext">http://www.impmc.upmc.fr/~paulatto/</a> -
<a href="https://anharmonic.github.io/" class="x_moz-txt-link-freetext">https://anharmonic.github.io/</a><br>
23-24/423 B115, 4 place Jussieu 75252 Paris CX 05<small></small></small></div>
</div>
</body>
</html>