<div dir="ltr">Meta GGA's have been in Wien2k for I think at least a decade.<div><br></div><div>What exactly are you trying to do? If you want to see sections of the kinetic energy density then, by hand, you can use lapw5 for real space and (probably) lapw3 for reciprocal space. You will need to look in the UG and ensure that you are using the correct normalizations.</div><div><br></div><div>If, somehow, you want critic2 to do something with it then you need to create case.inm_vresp (just copy from case.inm) and run a single iteration with a converged density to create the relevant case.vrespsum/up/dn files.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 11, 2018 at 6:54 AM, Gavin Abo <span dir="ltr"><<a href="mailto:gsabo@crimson.ua.edu" target="_blank">gsabo@crimson.ua.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div text="#000000" bgcolor="#FFFFFF">
<p>At <a class="m_-6057817230025761251moz-txt-link-freetext" href="https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_aoterodelaroza_critic2_blob_master_doc_user-2Dguide.txt&d=DwMDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=0sYGUIbhGzz1AlZTZCpwptGMlbSTSo2riKIQW1vDd0U&s=BXt7dvSmqQYoV6czdOq45IWbnZHvqvmK2DjZroyGNMA&e=" target="_blank">
https://github.com/<wbr>aoterodelaroza/critic2/blob/<wbr>master/doc/user-guide.txt</a>, there is:</p>
<p>* meta-GGA: xc(rho,grad,lapl,tau,idx)<br>
</p>
<p>where rho is the electron density expression, grad is its gradient, lapl is its Laplacian and tau is the kinetic energy density<br>
</p>
In the WIEN2k 17.1 usersguide, there is on page 62 in the section "5.1.1 Main execution script (x lapw)":<br>
<br>
-vresp-> runs lapw2/mixer and creates case.vrespval/sum (TAU/meta-GGA)<br>
<br>
The WIEN2k updates page [ <a class="m_-6057817230025761251moz-txt-link-freetext" href="https://urldefense.proofpoint.com/v2/url?u=http-3A__susi.theochem.tuwien.ac.at_reg-5Fuser_updates_&d=DwMDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=0sYGUIbhGzz1AlZTZCpwptGMlbSTSo2riKIQW1vDd0U&s=-1-IXIlwjXwaobfAeW_nfp4Md9SH12jgcGSE7dNJdTU&e=" target="_blank">
http://susi.theochem.tuwien.<wbr>ac.at/reg_user/updates/</a> ] has:<br>
<br>
<b>VERSION_16.1: 1.12.2016 </b><br>
<br>
SRC_lapw0: new GGA and meta-GGA functionals <br>
<br>
So maybe the kinetic energy density can be created using the -vresp switch.  However, it looks like it may not be possible with WIEN2k 13.1 as it looks like the meta-GAA was not implemented until WIEN2k 16.1.<br>
<br>
Might as well upgrade to WIEN2k 17.1 if that is the case, or perhaps to avoid many WIEN2k 17.1 patches [
<a class="m_-6057817230025761251moz-txt-link-freetext" href="https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_gsabo_WIEN2k-2DPatches_tree_master_17.1&d=DwMDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=0sYGUIbhGzz1AlZTZCpwptGMlbSTSo2riKIQW1vDd0U&s=vTej_ew7uxsLzT8wNQnv55GnufMneEUlj65Ec0qdQNo&e=" target="_blank">
https://github.com/gsabo/<wbr>WIEN2k-Patches/tree/master/17.<wbr>1</a> ], wait until WIEN2k 18 is released [
<a class="m_-6057817230025761251moz-txt-link-freetext" href="https://urldefense.proofpoint.com/v2/url?u=https-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_msg17501.html&d=DwMDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=0sYGUIbhGzz1AlZTZCpwptGMlbSTSo2riKIQW1vDd0U&s=95T3kuUJtoQrdjeXxumZZlbGHzLTakf0U-IthfFuW4I&e=" target="_blank">
https://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/msg17501.html</a> ].<span class=""><br>
<br>
<div class="m_-6057817230025761251moz-cite-prefix">On 6/11/2018 4:56 AM, Bon, Marta wrote:<br>
</div>
<blockquote type="cite">
<pre>Tau is the kinetic energy density. It should be possible to print it out with Wien2K (but I don't understand how).
I think with lapw0, but again I am stuck with the preparation of the case.in0 file.


Marta</pre>
</blockquote>
<br>
</span></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>
</div>