<HTML><HEAD></HEAD>
<BODY dir=ltr>
<DIV dir=ltr>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">
<DIV>Hi Pascal</DIV>
<DIV> </DIV>
<DIV>Thanks for replying. I agree with you that building a PP consistent with
the metaGGA functional would solve most of the instability problems. But I do
not have *any* experience building PPs (which, by the way, seems to be not
trivial to me, even less trivial for metaGGA), so I followed <A
title=http://dx.doi.org/10.1016/j.cpc.2013.02.020 style="TEXT-ALIGN: justify"
href="http://dx.doi.org/10.1016/j.cpc.2013.02.020">http://dx.doi.org/10.1016/j.cpc.2013.02.020</A>
and tried the trick. I was just wondering if any general strategy for
convergence exists. </DIV>
<DIV> </DIV>
<DIV>Regarding point 2, I asked one of the authors of the aforementioned paper
and he declared not the have seen such a problem before, so I guess that the
kinetic terms are not responsible. I have overcome the problem by using a dense
kpoints path already in scf calculations and calculating the bands without the
standard intermediate nscf run, but it would be nice if the developers could
have a look and check. </DIV>
<DIV> </DIV>
<DIV>Again, thank you very much for your interest.</DIV>
<DIV> </DIV>
<DIV>Juanjo</DIV>
<DIV> </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">Juan J.
Meléndez <BR>Associate Professor<BR>Department of Physics · University of
Extremadura<BR>Avda. de Elvas, s/n 06006 Badajoz (Spain)<BR>Phone: +34 924 28 96
55<BR>Fax: +34 924 28 96 51<BR>Email: melendez@unex.es<BR>Web:
http://materiales.unex.es/miembros/personal/jj-melendez/Index.html</DIV>
<DIV
style="FONT-SIZE: small; FONT-FAMILY: 'Calibri'; FONT-WEIGHT: normal; COLOR: #000000; FONT-STYLE: normal; TEXT-DECORATION: none; DISPLAY: inline">
<DIV style="FONT: 10pt tahoma">
<DIV><FONT size=3 face=Calibri></FONT> </DIV>
<DIV style="BACKGROUND: #f5f5f5">
<DIV style="font-color: black"><B>From:</B> <A title=pascal.boulet@univ-amu.fr
href="mailto:pascal.boulet@univ-amu.fr">Pascal BOULET</A> </DIV>
<DIV><B>Sent:</B> Thursday, October 30, 2014 7:15 PM</DIV>
<DIV><B>To:</B> <A title=pw_forum@pwscf.org
href="mailto:pw_forum@pwscf.org">PWSCF Forum</A> </DIV>
<DIV><B>Subject:</B> Re: [Pw_forum] TB09 metaGGA</DIV></DIV></DIV>
<DIV> </DIV></DIV>
<DIV
style="FONT-SIZE: small; FONT-FAMILY: 'Calibri'; FONT-WEIGHT: normal; COLOR: #000000; FONT-STYLE: normal; TEXT-DECORATION: none; DISPLAY: inline">
<DIV>
<DIV style="FONT-SIZE: 13px; FONT-FAMILY: georgia; COLOR: #000000">hello,
<DIV> </DIV>
<DIV>I have never tried meta-GGA functionals with QE, so I am guessing. </DIV>
<DIV> </DIV>
<DIV>For point 1, one possible reason is the inadequacy between the pseudo
potential and the functional. I guess you have build a special pseudopotential
for this functional.</DIV>
<DIV> </DIV>
<DIV>Point 2: Could it be because the meta-GGA kinetic terms are not saved in
the checkpoint file and hence not readable when doing the
nscf?<BR><BR><BR>Pascal</DIV>
<DIV> </DIV>
<DIV>
<DIV> </DIV>
<DIV dir=ltr>"Juanjo Meléndez" <melendez@unex.es> wrote:
<BLOCKQUOTE cite="<BA0A7E01AFE046D08178A05800725F9D@Otello>" type="cite">
<DIV dir=ltr>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">
<DIV>Hi all</DIV>
<DIV> </DIV>
<DIV align=justify>I am starting to work with metaGGA functionals in QE (v.
5.1 + libXc), following Èric Germaneau and colleagues’ work (<A
title=http://dx.doi.org/10.1016/j.cpc.2013.02.020
href="http://dx.doi.org/10.1016/j.cpc.2013.02.020">http://dx.doi.org/10.1016/j.cpc.2013.02.020</A>).
I know that metaGGA is little tested, but I would be pleased if somebody could
answer a couple of questions:</DIV>
<DIV align=justify> </DIV>
<DIV align=justify>1) Convergence is quite tricky. Up to know, I got some
results in simple systems using a very small mixing_beta (within 0.01–0.05),
not too many bands (just one or two conduction bands, never converged for more
bands) and restarting from a previously fully converged calculation. I got
results only for carbon, silicon and germanium, never got any convergence for
any binary simple compound (like LiF). Does anyone have additional hints for
metaGGA calculations to converge?</DIV>
<DIV align=justify> </DIV>
<DIV align=justify>2) In addition, I get something wrong with nscf
calculations. After convergence of a scf run, I get a list of bands at each
k-point, as usual –nothing strange here. Then I made a path of k-points to get
the band structure and run a nscf calculation. This finishes fine as well, but
the bands are not only different from those from the scf calculations, but
also unrealistically high. I am attaching the input and output files for both
the scf and nscf runs, as well as a couple of eigenval.xml files taken after
nscf. Does anybody know how to manage this? Or could this be a bug in the
code?</DIV>
<DIV> </DIV>
<DIV>Thanks a lot in advance</DIV>
<DIV> </DIV>
<DIV>Take care</DIV>
<DIV> </DIV>
<DIV>Juanjo</DIV>
<DIV> </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">Juan J.
Meléndez <BR>Associate Professor<BR>Department of Physics · University of
Extremadura<BR>Avda. de Elvas, s/n 06006 Badajoz (Spain)<BR>Phone: +34 924 28
96 55<BR>Fax: <A href="tel:+34924289651" target=_blank><A
href="tel:+34924289651" target=_blank>+34 924 28 96 51</A></A><BR>Email:
melendez@unex.es<BR>Web: <A
href="http://materiales.unex.es/miembros/personal/jj-melendez/Index.html">http://materiales.unex.es/miembros/personal/jj-melendez/Index.html</A></DIV>
<DIV
style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000"> </DIV>
<DIV
style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">****************************</DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">“All
those who look stupid actually are, as also are half of those who do not look
like” (F. de Quevedo)</DIV>
<DIV
style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000"> </DIV>
<DIV
style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">****************************</DIV></DIV></DIV><BR><BR>
<HR
style="BORDER-TOP: medium none; HEIGHT: 1px; BORDER-RIGHT: medium none; BORDER-BOTTOM: medium none; COLOR: #909090; BORDER-LEFT: medium none; WIDTH: 99%; BACKGROUND-COLOR: #b0b0b0">
<TABLE
style="BORDER-TOP: medium none; BORDER-RIGHT: medium none; BORDER-COLLAPSE: collapse; BORDER-BOTTOM: medium none; COLOR: #000000; BORDER-LEFT: medium none">
<TBODY>
<TR>
<TD
style="BORDER-TOP: medium none; BORDER-RIGHT: medium none; BORDER-BOTTOM: medium none; PADDING-BOTTOM: 0px; PADDING-TOP: 0px; PADDING-LEFT: 8px; BORDER-LEFT: medium none; PADDING-RIGHT: 15px"><A
href="http://www.avast.com/"><IMG border=0
src="http://static.avast.com/emails/avast-mail-stamp.png"> </A></TD>
<TD>
<P
style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri','Verdana','Arial','Helvetica'; COLOR: #3d4d5a">Este
mensaje no contiene virus ni malware porque la protección de <A
href="http://www.avast.com/">avast! Antivirus</A> está activa.
</P></TD></TR></TBODY></TABLE><BR><PRE style="FONT-SIZE: 13px">_______________________________________________<BR>Pw_forum mailing list<BR><A href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</A><BR><A href="http://pwscf.org/mailman/listinfo/pw_forum" target=_blank>http://pwscf.org/mailman/listinfo/pw_forum</A></PRE></BLOCKQUOTE></DIV><BR><BR>
<DIV style="FONT-SIZE: 10px; FONT-FAMILY: arial; COLOR: #000000"><FONT
size=2><FONT size=1></FONT><BR>>-----------------<BR>Pascal
Boulet<BR>Aix-Marseille University <BR>MADIREL Laboratory<BR>Avenue
Normandie-Niemen <BR>13397 Marseille Cedex 20 <BR>Email:
pascal.boulet@univ-amu.fr <BR>Tel. <A href="tel:+33413551810" target=_blank>+33
413 55 18 10</A> <BR>Fax <A href="tel:+33413551850" target=_blank>+33 413 55 18
50</A></FONT><BR><BR></DIV></DIV></DIV></DIV>
<P>
<HR>
_______________________________________________<BR>Pw_forum mailing
list<BR>Pw_forum@pwscf.org<BR>http://pwscf.org/mailman/listinfo/pw_forum</DIV></DIV></DIV>
<br /><br />
<hr style='border:none; color:#909090; background-color:#B0B0B0; height: 1px; width: 99%;' />
<table style='border-collapse:collapse;border:none;'>
<tr>
<td style='border:none;padding:0px 15px 0px 8px'>
<a href="http://www.avast.com/">
<img border=0 src="http://static.avast.com/emails/avast-mail-stamp.png" />
</a>
</td>
<td>
<p style='color:#3d4d5a; font-family:"Calibri","Verdana","Arial","Helvetica"; font-size:12pt;'>
Este mensaje no contiene virus ni malware porque la protección de <a href="http://www.avast.com/">avast! Antivirus</a> está activa.
</p>
</td>
</tr>
</table>
<br />
</BODY></HTML>