No subject

http://www.democritos.it/pipermail/pw_forum/2008-September/010123.html
http://www.democritos.it/pipermail/pw_forum/2010-February.txt

I understand that the third derivatives were never implemented and
that only the LDA part was calculated, for the Perdew Zunger
parameterization. Is this third derivative the third derivative that
is mentioned in the routine PH/raman_mat.f90 (i.e, the third
derivative of the energy with respect to the atomic displacements and
with respect to two electric fields, as described in Lazzeri and
Mauri, PRL 90,036401)?

I am interested in getting a feel for what would be needed to
implement the third derivatives with GGA into the code, as well as
what sort of errors one might expect when using the third-order term
as currently coded with LDA pseudopotentials (comments in the forums
have hinted that the errors are likely small).  The shift in the Raman
spectra for cd-Si is pretty minor when using a GGA pseudo compared to
that of an LDA calculation, but I wanted to be able to quantify the
error better for other systems where comparing the GGA raman spectra
to one calculated within LDA was not possible (e.g., when the system
doesn't exhibit a gap in LDA).

Thanks for the help!

Best,
Brad
UC Berkeley

On Tue, Nov 9, 2010 at 5:00 PM, Brad Malone <brad.malone at gmail.com> wrote:

> Hi, I am seeking more information on the error
>
> > from phq_setup : error # 1
> > third order derivatives not implemented with GGA
>
> From the forums:
>
> http://www.democritos.it/pipermail/pw_forum/2008-September/010123.html
> http://www.democritos.it/pipermail/pw_forum/2010-February.txt
>
> I understand that the third derivatives were never implemented and
> that only the LDA part was calculated, for the Perdew Zunger
> parameterization. Is this third derivative the third derivative that
> is mentioned in the routine PH/raman_mat.f90 (i.e, the third
> derivative of the energy with respect to the atomic displacements and
> with respect to two electric fields, as described in Lazzeri and
> Mauri, PRL 90,036401)?
>
> I am interested in getting a feel for what would be needed to
> implement the third derivatives with GGA into the code, as well as
> what sort of errors one might expect when using the third-order term
> as currently coded with LDA pseudopotentials (comments in the forums
> have hinted that the errors are likely small).  The shift in the Raman
> spectra for cd-Si is pretty minor when using a GGA pseudo compared to
> that of an LDA calculation, but I wanted to be able to quantify the
> error better for other systems where comparing the GGA raman spectra
> to one calculated within LDA was not possible (e.g., when the system
> doesn't exhibit a gap in LDA).
>
> Thanks for the help!
>
> Best,
> Brad
> UC Berkeley
>

--000e0cd1e7da1a89630494a936f2
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable

It appears that my last email was severely truncated on the mailing list (e=
ven though it looks fine in my outbox). Below is what is should have been:<=
div>-----------------------</div><div><meta http-equiv=3D"content-type" con=
tent=3D"text/html; charset=3Dutf-8"><span class=3D"Apple-style-span" style=
=3D"font-family: arial, sans-serif; font-size: 13px; border-collapse: colla=
pse; color: rgb(51, 51, 51); ">Hi, I am seeking more information on the err=
or<br>
<br>&gt; from phq_setup : error # 1<br>&gt; third order derivatives not imp=
lemented with GGA<br><br>From the forums:<br><br><a href=3D"http://www.demo=
critos.it/pipermail/pw_forum/2008-September/010123.html" target=3D"_blank" =
style=3D"color: rgb(61, 84, 89); ">http://www.democritos.it/pipermail/pw_fo=
rum/2008-September/010123.html</a><br>
<a href=3D"http://www.democritos.it/pipermail/pw_forum/2010-February.txt" t=
arget=3D"_blank" style=3D"color: rgb(61, 84, 89); ">http://www.democritos.i=
t/pipermail/pw_forum/2010-February.txt</a><br><br>I understand that the thi=
rd derivatives were never implemented and<br>
that only the LDA part was calculated, for the Perdew Zunger<br>parameteriz=
ation. Is this third derivative the third derivative that<br>is mentioned i=
n the routine PH/raman_mat.f90 (i.e, the third<br>derivative of the energy =
with respect to the atomic displacements and<br>
with respect to two electric fields, as described in Lazzeri and<br>Mauri, =
PRL 90,036401)?<br><br>I am interested in getting a feel for what would be =
needed to<br>implement the third derivatives with GGA into the code, as wel=
l as<br>
what sort of errors one might expect when using the third-order term<br>as =
currently coded with LDA pseudopotentials (comments in the forums<br>have h=
inted that the errors are likely small). =A0The shift in the Raman<br>spect=
ra for cd-Si is pretty minor when using a GGA pseudo compared to<br>
that of an LDA calculation, but I wanted to be able to quantify the<br>erro=
r better for other systems where comparing the GGA raman spectra<br>to one =
calculated within LDA was not possible (e.g., when the system<br>doesn&#39;=
t exhibit a gap in LDA).<br>
<br>Thanks for the help!<br><br>Best,<br>Brad<br>UC Berkeley</span><br><br>=
<div class=3D"gmail_quote">On Tue, Nov 9, 2010 at 5:00 PM, Brad Malone <spa=
n dir=3D"ltr">&lt;<a href=3D"mailto:brad.malone at gmail.com">brad.malone at gmai=
l.com</a>&gt;</span> wrote:<br>
<blockquote class=3D"gmail_quote" style=3D"margin:0 0 0 .8ex;border-left:1p=
x #ccc solid;padding-left:1ex;">Hi, I am seeking more information on the er=
ror<br>
<br>
&gt; from phq_setup : error # 1<br>
&gt; third order derivatives not implemented with GGA<br>
<br>