<div dir="auto"><div>Sorry, I forgot to mention that the results or dispersion were obtained using qe.6.3-bkp.</div><div dir="auto"><br></div><div dir="auto">Best<br><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Wed, Mar 13, 2019, 15:04 Prof. Raphael Longuinhos Monteiro Lobato <<a href="mailto:raphael.lobato@ufla.br" target="_blank" rel="noreferrer">raphael.lobato@ufla.br</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
Back on the celldm post (sorry for the six months delay to reply ... I <br>
will bring back some parts to remember):<br>
<a href="http://lists.quantum-espresso.org/pipermail/users/2018-October/041418.html" rel="noreferrer noreferrer noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/2018-October/041418.html</a><br>
<br>
and<br>
<br>
<a href="http://lists.quantum-espresso.org/pipermail/users/2018-October/041419.html" rel="noreferrer noreferrer noreferrer" target="_blank">http://lists.quantum-espresso.org/pipermail/users/2018-October/041419.html</a><br>
<br>
Previously, I wrote:<br>
+==<br>
ph.x from pre-release q-e-qe-6.3Max reads from outdir unexpected values <br>
for celldm and writes these in the header of .dyn . It leads to <br>
"rescaled" distances in .axsf and polarizability in the dynmat.x out <br>
(omega, IR, Raman tensor and depolarization factor are as expected). <br>
Hand-correcting the .dyn header (correct values for celldm) solves it.<br>
<br>
Tests in v6.2 and 6.2.2 leads to expected results (they read from outdir <br>
the expected values for celldm). ph.x prints wrong celldm while reading <br>
from outdir ... it is not related to io_dyn_mat write_dyn_mat_header ... <br>
also, the .xml files are the same.<br>
+==<br>
<br>
Prof. Giannozzi mentioned that<br>
<a href="https://gitlab.com/QEF/q-e/commit/9a579539b65bd1502d1f28747e483a1d1059f5f8" rel="noreferrer noreferrer noreferrer" target="_blank">https://gitlab.com/QEF/q-e/commit/9a579539b65bd1502d1f28747e483a1d1059f5f8</a><br>
<br>
could have solved the issue. Maybe I'm missing something, but it appears <br>
it did not.<br>
<br>
The issue on the celldm can lead to other unexpected differences on <br>
post-processing the dyn.<br>
<br>
In my case, the high-frequency (low frequency is okay) phonon dispersion <br>
near Gamma calculated with the original (wrong) celldm values and the <br>
ones using hand-correcting values are somewhat different. As example, <br>
the high-frequency LO mode obtained with the wrong cell_dm seems <br>
quadratic (or some d^2f/dk^2 >0 dependency) like using no LR <br>
electrostatics, while it becomes linear-like (I would expect some thing <br>
like that) when using the correct cell_dm. It is show in the images in <br>
appended. Also, there are some "discontinuities" in the branches I not <br>
used to see on results. Also, there are some extra inflection points in <br>
the region I show in the graphs that I would not expect (these appears <br>
using no LR electrostatics).<br>
<br>
Note: The system is expected to be polar semiconductor. In the graphs nq <br>
07x07x01 and 14x14x01 uses nk is 14x14x1 (spacing of ~0.03*2pi/A). <br>
Later, I used both nk and nq of 16x16x01. Orange diamonds are direct 2D <br>
DFPT calculations (nk 14x14x01 or 16x16x01 leads to identical results in <br>
these direct 2D DFPT calculations). The highest frequency modes from <br>
direct 2D DFPT and from interpolation of dynmat are some what different <br>
(maybe the interpolation is not converged). Data with cell_dm corrected <br>
values are also indicated.<br>
<br>
I will repeat all these using those versions I reported correctly <br>
reading cell_dm values.<br>
<br>
Hope all it is not a "false alarm" (some unnoticed error of my part). <br>
Sorry bring problems, not solutions.<br>
<br>
Best,<br>
<br>
-- <br>
Associate Professor Raphael Longuinhos Monteiro Lobato<br>
CV: <a href="http://lattes.cnpq.br/1475483971471395" rel="noreferrer noreferrer noreferrer" target="_blank">http://lattes.cnpq.br/1475483971471395</a>.<br>
Research Group site: <a href="http://jenainassoares2.wixsite.com/nanomat" rel="noreferrer noreferrer noreferrer" target="_blank">http://jenainassoares2.wixsite.com/nanomat</a> (Portuguese only - English version on the way)<br>
Departamento de Fisica, Universidade Federal de Lavras<br>
Av. Dr. Silvio Menicucci 1001, bairro Kenedy<br>
Post office box 3037<br>
37200-000 - Lavras, MG, Brazil<br>
<br>
</blockquote></div></div></div>