<div dir="ltr">Dear prof. Hutter<div><br></div><div>thank you for your reply. </div><div>I emailed our issue to the libxc mailing list.</div><div>I have just received what is stated below:</div><div><br></div><div><i><span style="font-family: arial, sans-serif; font-size: 12.8px;">Hi,</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;">I also had such problem with TMx. In addition to the usual</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;">tau = max(tau,tauw)</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;">(tau and tauw are the KS and Weizsäcker kinetic energies)</span></i></div><div><i><span style="font-family: arial, sans-serif; font-size: 12.8px;">that should be always used in general, I also needed</span></i><br></div><div><i><span style="font-family: arial, sans-serif; font-size: 12.8px;">to add something like</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;">if (abs(tau) .lt. 1.0d-18) tau = 1.0d-18</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;">for TMx because of divisons tauw/tau. It should be the same</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;">problem for you.</span><br style="font-family: arial, sans-serif; font-size: 12.8px;"><span style="font-family: arial, sans-serif; font-size: 12.8px;"><br></span></i></div><div><i><span style="font-family: arial, sans-serif; font-size: 12.8px;">Best regards,</span></i><br style="font-family: arial, sans-serif; font-size: 12.8px;"><br>Do I need to change the CP2K code?</div><div>If yes, could you please point me where?</div><div><br></div><div>Bests</div><div>Luca </div><div><br><br><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><br> 
<br></blockquote></div></div>