<div dir="ltr"><div dir="ltr"><div><div>Dear Paolo: <br><br>Thank you very much for helping. So there is no way I can do the lattice
 constance optimization right now with metaGGA functionals in QE, unless
 somebody can re-design the pp files. Is my understanding correct?<br><br></div>Best regards,<br><br></div>Hong Tang<br><br><div><div><pre><br><br>______________________________<wbr>____________<br><br>There isn't much to explain: nobody has ever written down and
implemented  the formula for stress in the case of meta-GGA. In QE
meta-GGA is implemented (sort of) only for norm-conserving
pseudopotential, so there is no need (bot that I know) to store the
kinetic energy density in the PP file<br><br><br><br>Paolo

On Wed, Jun 28, 2017 at 4:34 AM, Hong Tang <<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">tanghong992 at gmail.com</a>> wrote:
><i> Dear Paolo:
</i>><i>
</i>><i>
</i>><i>
</i>><i> Thank you so much for helping.
</i>><i>
</i>><i> I tried grep -n “stress” Si.rx.out
</i>><i>
</i>><i> What I got is below:
</i>><i>
</i>><i> [<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">hongtang at login4</a> results]$ grep -n "stress" Si.rx.out
</i>><i>
</i>><i> 1797:     entering subroutine stress ...
</i>><i>
</i>><i> 1799:     Message from routine stress:
</i>><i>
</i>><i> 1800:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 3128:     entering subroutine stress ...
</i>><i>
</i>><i> 3130:     Message from routine stress:
</i>><i>
</i>><i> 3131:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 4442:     entering subroutine stress ...
</i>><i>
</i>><i> 4444:     Message from routine stress:
</i>><i>
</i>><i> 4445:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 5764:     entering subroutine stress ...
</i>><i>
</i>><i> 5766:     Message from routine stress:
</i>><i>
</i>><i> 5767:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 7555:     entering subroutine stress ...
</i>><i>
</i>><i> 7557:     Message from routine stress:
</i>><i>
</i>><i> 7558:     Meta-GGA and stress not implemented
</i>><i>
</i>><i>
</i>><i>
</i>><i> It shows MetaGGA and stress not implemented.  Maybe this is the problem
</i>><i> there.
</i>><i>
</i>><i> I want to do structure optimization with metagga functional (TPSS here).  In
</i>><i> VASP, I know the pseudo potentials (pps) of the elements should have
</i>><i> “kinetic energy density”.  In quantum espresso, I don’t know what should be
</i>><i> in the pps. The pps I chose for Si is for tpss, as you can see from the
</i>><i> input file.
</i>><i>
</i>><i>
</i>><i>
</i>><i> Please explain to me a bit more clear how to fix the problem. I appreciate
</i>><i> your helps very much.
</i>><i>
</i>><i>
</i>><i>
</i>><i> Hong Tang
</i>><i>
</i>><i> Temple University
</i>><i>
</i>><i> Philadelphia, PA 19122 US</i></pre></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 28, 2017 at 10:37 PM, Hong Tang <span dir="ltr"><<a href="mailto:tanghong992@gmail.com" target="_blank">tanghong992@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Dear Paolo: <br><br>Thank you very much for helping. So there is no way I can do the lattice
 constance optimization right now with metaGGA functionals in QE, unless
 somebody can re-design the pp files. Is my understanding correct?<br><br></div>Best regards,<br><br></div>Hong Tang<br><br><div><div><pre><br><br>______________________________<wbr>____________<br><br>There isn't much to explain: nobody has ever written down and
implemented  the formula for stress in the case of meta-GGA. In QE
meta-GGA is implemented (sort of) only for norm-conserving
pseudopotential, so there is no need (bot that I know) to store the
kinetic energy density in the PP file<br><br><br><br>Paolo

On Wed, Jun 28, 2017 at 4:34 AM, Hong Tang <<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">tanghong992 at gmail.com</a>> wrote:
><i> Dear Paolo:
</i>><i>
</i>><i>
</i>><i>
</i>><i> Thank you so much for helping.
</i>><i>
</i>><i> I tried grep -n “stress” Si.rx.out
</i>><i>
</i>><i> What I got is below:
</i>><i>
</i>><i> [<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">hongtang at login4</a> results]$ grep -n "stress" Si.rx.out
</i>><i>
</i>><i> 1797:     entering subroutine stress ...
</i>><i>
</i>><i> 1799:     Message from routine stress:
</i>><i>
</i>><i> 1800:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 3128:     entering subroutine stress ...
</i>><i>
</i>><i> 3130:     Message from routine stress:
</i>><i>
</i>><i> 3131:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 4442:     entering subroutine stress ...
</i>><i>
</i>><i> 4444:     Message from routine stress:
</i>><i>
</i>><i> 4445:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 5764:     entering subroutine stress ...
</i>><i>
</i>><i> 5766:     Message from routine stress:
</i>><i>
</i>><i> 5767:     Meta-GGA and stress not implemented
</i>><i>
</i>><i> 7555:     entering subroutine stress ...
</i>><i>
</i>><i> 7557:     Message from routine stress:
</i>><i>
</i>><i> 7558:     Meta-GGA and stress not implemented
</i>><i>
</i>><i>
</i>><i>
</i>><i> It shows MetaGGA and stress not implemented.  Maybe this is the problem
</i>><i> there.
</i>><i>
</i>><i> I want to do structure optimization with metagga functional (TPSS here).  In
</i>><i> VASP, I know the pseudo potentials (pps) of the elements should have
</i>><i> “kinetic energy density”.  In quantum espresso, I don’t know what should be
</i>><i> in the pps. The pps I chose for Si is for tpss, as you can see from the
</i>><i> input file.
</i>><i>
</i>><i>
</i>><i>
</i>><i> Please explain to me a bit more clear how to fix the problem. I appreciate
</i>><i> your helps very much.
</i>><i>
</i>><i>
</i>><i>
</i>><i> Hong Tang
</i>><i>
</i>><i> Temple University
</i>><i>
</i>><i> Philadelphia, PA 19122 US</i></pre></div></div></div>
</blockquote></div><br></div>