<div dir="ltr"><div><div><div>Dear Afshin,<br></div>it is not possible to convert total 
energy to stress; stress can be computed from the same ingredients 
(wavefunctions, self-consistent potential) as the total energy, but not 
just from the value of energy itself.<br><br>pw.x can compute the stress
 of any give structure, you have to set "tstress=.true." in the 
&control namelist, but of course you have to do it before starting 
the calculation.<br><br></div>kind regards.<br><br></div>p.s. repeating your question every 8 hours won't increase the chance of getting an answer, especially on Sunday<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 25 September 2016 at 17:49, Afshin Arjangmehr <span dir="ltr"><<a href="mailto:ms.arjangmehr@gmail.com" target="_blank">ms.arjangmehr@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"><span class=""><font style="font-size:12.8px" face="verdana, sans-serif">Dear QE developers and users,</font><div style="font-size:12.8px"><font face="verdana, sans-serif"><br></font></div><div style="font-size:12.8px"><font face="verdana, sans-serif">I have been able to calculate the total energy of crystalline structures under various strain tensors using pwscf. Now I intend to convert the total energy to stress tensors and then study the mechanical properties (either elastic or plastic) of different materials under strain (I have seen several papers that directly use E vs. Strain to calculate 2nd, 3rd, 4th, ... elastic constants and find yield & ultimate stress, but I believe that as Caro et al. <span style="background-color:rgb(245,245,245);color:rgb(51,51,51);white-space:pre-wrap"><a href="http://dx.doi.org/10.1088/0953-8984/25/2/025803" target="_blank">http://dx.doi.org/10.1088/<wbr>0953-8984/25/2/025803</a></span> mentione<wbr>d, using Stress-Strain curve produces much accurate results).</font></div><div style="font-size:12.8px"><font face="verdana, sans-serif">So, can anyone kindly give me some recommendation on how to convert Etot to Stress tensor? or Is there any tool to do this conversion? </font><span style="font-family:verdana,sans-serif">, since simply using (-dEtot/dVol) is not the solution</span></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Regards</div><div><br></div></span><span class="">-- <br><div><div dir="ltr"><div style="font-family:arial;font-size:small"><font size="1">With Best Regards</font></div><div style="font-family:arial;font-size:small"><font size="1">Afshin Arjhangemehr</font></div><div style="font-family:arial;font-size:small"><font size="1">PhD student in Radiation Application</font></div><div style="font-family:arial;font-size:small"><font size="1">Shahid Beheshti University G.C, Tehran, IRAN</font></div><div style="font-family:arial;font-size:small"><font size="1">(+98) 912 439 20 64</font></div></div></div>
</span></div>
<br>______________________________<wbr>_________________<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/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div>Dr. Lorenzo Paulatto </div><div>IdR @ IMPMC -- CNRS & Université Paris 6</div><div>phone: +33 (0)1 44275 084 / skype: paulatz</div><div>www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a></div><div>mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</div></div>
</div>