<div dir="ltr"><div>The calculation of forces in DFT+U may take a significant amount of time when there are many atoms with U correction. Recent versions of QE are much faster than earlier ones in this respect. Which version of QE are you using, how many atoms (with and without Hubbard U) and on how many processors are you running?<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 21, 2015 at 4:21 AM, Mostafa Youssef <span dir="ltr"><<a href="mailto:myoussef@mit.edu" target="_blank">myoussef@mit.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Dear all,<br>
<br>
<br>
I performed ionic relaxation for 3x3x3 unit cells of a semiconducting oxide under electric field (lelfield =.true.) using ultrasoft PP and 1x1x3 MP K-point grid where 3 is in the direction of the electric field. These are GGA+U calculations. I noticed in the
 time report that force calculations time is becoming even more than the electrons time as shown below:<br>
<br>
     init_run     :      8.32s CPU      8.91s WALL (       1 calls)<br>
     electrons    :  42723.25s CPU  44995.21s WALL (       5 calls)<br>
     update_pot   :      2.62s CPU      3.19s WALL (       4 calls)<br>
     forces       :  66735.95s CPU  66943.64s WALL (       5 calls)<br>
<br>
This is first time in which I find PWscf use more time in force calculations than in electrons. For the same system and without electric field , forces time would be about 0.2 of electrons time.  Is this a normal behavior and is there a recommended possibility
 to reduce force calculation time?  <br>
<br>
I attempted using Norm Conserving PP's but unfortunately I could not get the code pass the 1 ionic step because of the S matrix not being positive definite.<br>
<br>
I appreciate your thoughts on this.<br>
<br>
Regards,<br>
Mostafa Youssef <br>
MIT<br>
<br>
</div>
</div>

<br>_______________________________________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div>