Javad,<br><br>It is unlikely that you can converge the stress more than a couple of kilobar. (It is calculated in a complete different way than in classical MD.) As hinted at in the previous e-mail, the pressure will converge slower than the energy within the PW pseudopotential scheme. In other words, if the calculation needs only 40 Ry to converge the energy well, don't be surprised if you have to go 60 - 80 Ry to converge you stress tensor.
<br><br>You can also calculate the pressure (roughly) as<br>P = -dE/dV<br><br>where E and V are the total energy and volume. You would need to work out the unit conversion to kbar.<br><br>It is good to look back in the literature to papers that discuss this matter in more detail:
<br>Phys. Rev. B 32, 3780 (1985).<br>Phys. Rev. B 52, 8160 (1995).<br><br><br><div><span class="gmail_quote">On 12/28/06, <b class="gmail_sendername">Katalin Gaal-Nagy</b> <<a href="mailto:katalin.gaal-nagy@physik.uni-regensburg.de">
katalin.gaal-nagy@physik.uni-regensburg.de</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Dear Javad,<br><br>one reason for the discrepancy could be that you are not converged enough
<br>with respect to the number of plane waves, since the stress is very<br>sensitive to that quantity. If you have calculated also the stress during<br>your convergence tests with respect to the ecut, you can estimate the
<br>error in the stress.<br>Of course one expects that the stress vanishes at the minimum of the total<br>energy.<br><br>Another reason coul be that you are using the experimental lattice<br>constant, which could yield a energy which is not at the minimum of the
<br>total energy due to the over/underestimation of the lattice constant using<br>LDA/GGA. Thus, I suggest to relax also the lattice constant together with<br>the c/a in order to obtain the values at the energy minimum of your
<br>calculation.<br><br>However, 10 kbar is a rather small value which could be somehow<br>acceptable.<br><br>But maybe someone else could give you further hints.<br><br>Best wishes,<br>Katalin<br><br>On Thu, 28 Dec 2006, Javad hashemifar wrote:
<br><br>> Dear PWscf users<br>><br>> I am trying to relax the c lattice parameter of a tetragonal supercell<br>> that its a lattice parameter is fixed to some experimental value.<br>> Therefore I am calculating this supercell at different values of c
<br>> (celldm(3)) and looking at the total energy and stress.<br>> My feeling is that these two parameter should have consistent behavior<br>> and both could be used to find the optimized c. I mean I expect that<br>
> when total energy is minimized, then the stress should vanish but it<br>> is not the case for my calculation and I observe that when the total<br>> energy is minimum, there is considerable value of stress (~ -10 kbar)
<br>> in the system.<br>> I would appreciate if some body explain what is the problem and<br>> if it is really true that stress and total energy have not the same<br>> behavior, which one should be used to find out relaxed lattice
<br>> parameters. I expect that total energy is the main criteria for<br>> relaxation of the lattice parameters ??<br>><br>> Best regards<br>> Javad Hashemifar<br>><br>_______________________________________________
<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote>
</div><br><br clear="all"><br>-- <br>Nichols A. Romero, Ph.D.<br>1613 Denise Dr. Apt. D<br>Forest Hill, MD 21050<br>443-567-8328 (C)<br>410-306-0709 (O)<br>