Sandro, I am replying to this old e-mail now because the funds have become available to do some work in this area. I looked in Paolo Focher's thesis and there are some very helpful equations, so thank you. However, the eqn.'s are limited to LDA.

<br><br>Does anyone know where I could find the relevant equation(s) for the GGA functional contribution to the stress tensor? I imagine it should not be that different from LDA, but I want to be sure.<br><br><br><div><span class="gmail_quote">

On 7/13/06, <b class="gmail_sendername">Scandolo Sandro</b> <<a href="mailto:scandolo@ictp.it">scandolo@ictp.it</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;">

A number of useful formulae can be found in Paolo Focher's thesis (1994)<br>(dowloadable from <a href="http://www.sissa.it/cm/phd.php">http://www.sissa.it/cm/phd.php</a> ), as well as in a<br>subsequent publication:<br>

Paolo Focher and Guido L. Chiarotti Ab-initio Molecular Dynamics simulations of structural phase transitions in ``Progress in Computational Physics of Matter'', eds. L. Reatto and F. Manghi, p. 1-42, World Scientific, Singapore, (1995).

<br><br>(I can fax you the relevant pages of the latter when I'm back to Trieste,<br>late August).<br><br>Best regards,<br>Sandro<br><br><br>> Hi,<br>><br>> I and another post-doc are working on implementing a non-local

> functional into PWSCF. It is of the general form: > E_xc[n(r)] = \int dr dr' n(r) phi(r,r') n(r') > > where n(r) and phi(r,r') are the charge density and a kernel, > respectively.

<br>><br>> Forces require calculating V_xc which is analytically complicated, but<br>> has already been done. The next quantity that we want to compute is<br>> the stress tensor.<br>><br>> Does any happen to have any notes on calculating the stress tensor in

<br>> a PW basis set?<br>><br>> I looked at the original paper on calculating the stress:<br>><br>> O. H. Nielsen and R. M. Martin, Phys. Rev. Lett. 50, 697 (1983).<br>><br>> Eqn. 3 contains the V_xc (mu_xc in their notation) and it also

<br>> contains epsilon_xc which is the exchange-correlation energy per<br>> particle which cannot be easily written down in a closed analytic form<br>> for a non-local functional (I think). Otherwise, I don't see any

<br>> higher derivatives of E_xc in the the expression for the stress.<br>><br>> Does anyone have any useful references?<br>><br>> Thanks,<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>> _______________________________________________<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>><br><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>