Dear Prof. Baroni<br> <br>Thanks again<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div style="word-wrap:break-word"><div><div class="im">
<blockquote type="cite"><div>Dear Prof. Baroni</div>
<div> </div>
<div>Thanks for your reply.</div>
<div>I wonder shouldn't the plot show the same frequency for a range of displacements?<span dir="ltr"> </span></div></blockquote><div><br></div></div>what range? why should it?</div><div><div class="im"><br></div></div>
</div></blockquote><div><br style="background-color: rgb(255, 255, 255);"><font style="background-color: rgb(255, 255, 255);">-- Because if we are in the
neighberhood of a minimum "where a quadratic could be fitted", any
displacement in that range should result the same frequency. By
displacement, I mean the same pattern but different in amplitude.</font><br> </div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap: break-word;">
<div><div class="im"><blockquote type="cite"><div><span dir="ltr">And the same freq. for the range to be the harmonic frequency? or the zero displacement limit, is the harmonic frequency?</span><br></div></blockquote><div>
<br></div></div>yes, by definition</div><div><div class="im"><br><blockquote type="cite"><div>
</div>
<div>And when I test the displacements with the same amplitude, but
different in sign (regarding the inward and outward atomic movements in
the basis) I do not get the same results. </div></blockquote><div><br></div></div>why should you?</div><div><div class="im"><br><blockquote type="cite"><div>Shouldn't the plot have
symmetry vs. the vertical axis?</div></blockquote><div><br></div></div>why?</div></div></blockquote><div><br><font style="background-color: rgb(255, 255, 255);">-- Consider for example a linear
chain of atoms, with 2 atoms per basis. The optical mode is the result
of the outward movement of the atoms (The one in +a/2 moves towards +x
and the one in -a/2 moves towards -x). I expect that if I displace the
atoms in an inward movement with the same amplitude as the outward, I
should see the same frequency (The one in +a/2 moves towards -x and the
one in -a/2 moves towards +x) because the 2 patterns present the same
mode, right? That is what I've done in the case of Diamond. Is my
expectation wrong?</font><br style="background-color: rgb(255, 255, 255);"> </div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap: break-word;">
<div><br></div><div>please, try you best to answer my questions, and I will make my best to let you understand</div><div><br></div><font color="#888888"><div>SB</div></font><div><div></div><div class="h5"><div><br></div>
</div>
</div></div></blockquote><div><br>Mohaddeseh <br><br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap: break-word;">
<div><div class="h5"><div><blockquote type="cite">
<div> </div>
<div>Thank You again<br>Mohaddeseh<br></div><br><br><div class="gmail_quote">On Mon, Aug 1, 2011 at 6:33 PM, Stefano Baroni <span dir="ltr"><<a href="mailto:baroni@sissa.it" target="_blank">baroni@sissa.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">Hi Mohaddesh. I do not understand why you say that the calculated frequency (which, I assume you define as 2 * Delta E / u^2 ???) does not converge as u->0. According to your plot, it does, and to a reasonable value. The finite slope of your plot around u=0 seems to me nothing but a manifestation of cubic anharmonicity ... What's wrong in this? SB<div>
<br><div><div><div></div><div><div>On Aug 1, 2011, at 3:01 PM, mohaddeseh abbasnejad wrote:</div><br></div></div><blockquote type="cite"><div><div></div><div><br>Dear users,<br clear="all">
<div style="text-align:justify;direction:ltr" dir="ltr"><br>I have
plotted the harmonic frequency (obtained within the frozen phonon
approach) of Diamond as a function of atomic displacements (small
displacements in the [1 1 1] direction, u(1)=±(u,u,u) and
u(2)=±(-u,-u,-u) are the displacements of two atoms in the basis). I
expect to see a range of displacements that tend to the similar
frequency, but there is no convergence in the plot. Also, the
displacements with the same amplitude but different in sign, do not tend
to the same frequency, which is not what I expect, too. The plot is attached. <br>I was wondering
if anyone could help me to find the source of error. I should expect the
mentioned symmetry, right? Why there is no convergence?</div>
<div style="text-align:justify;direction:ltr" dir="ltr"> </div>
<div style="text-align:justify;direction:ltr" dir="ltr">Thanks in advance.</div>
<div style="text-align:justify;direction:ltr" dir="ltr"> </div>
<div style="text-align:justify;direction:ltr" dir="ltr">Yours,</div>
<div style="text-align:justify;direction:ltr" dir="ltr">Mohaddeseh</div>---------------------------------------------------------<br><br>Mohaddeseh Abbasnejad, <br>Room No. 323, Department of Physics, <br>University of Tehran, North Karegar Ave.,<br>
Tehran, P.O. Box: 14395-547- IRAN <br>E-Mail: <a href="mailto:m.abbasnejad@gmail.com" target="_blank">m.abbasnejad@gmail.com</a><br>
Website: <a href="http://physics.ut.ac.ir/" target="_blank">http://physics.ut.ac.ir</a><br><br>---------------------------------------------------------<br><br>
</div></div><span><fig..pdf></span>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><div>
<span style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px"><font size="3"><span style="font-size:12px">---</span></font></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px">
<font size="3"><span style="font-size:12px">Stefano Baroni - SISSA</span></font><font size="3"><span style="font-size:12px"> </span></font><font size="3"><span style="font-size:12px">&</span></font><font size="3"><span style="font-size:12px"> </span></font><font size="3"><span style="font-size:12px">DEMOCRITOS National Simulation Center - Trieste</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px"><font size="3"><span style="font-size:12px"><a href="http://stefano.baroni.me/" target="_blank">http://stefano.baroni.me</a> <a href="tel:%5B%2B39%5D%20040%203787%20406" value="+390403787406" target="_blank">[+39] 040 3787 406</a> (tel) -528 (fax) / stefanobaroni (skype)</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;min-height:14px;font-size:12px"><br style="font-size:12px"></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;min-height:14px;font-size:12px">
<span style="font-size:14px"><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><font size="3"><span style="font-size:12px">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><font size="3"><span style="font-size:12px"><br></span></font></div></span></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px">
<font size="3"><span style="font-size:12px">Please, if possible, don't</span></font><font size="3"><span style="font-size:12px"> </span></font><font size="3"><span style="font-size:12px">send me MS Word or PowerPoint attachments</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px"><font size="3"><span style="font-size:12px">Why? See:</span></font><font size="3"><span style="font-size:12px"> </span></font><font color="#0023E9" size="3"><span style="font-size:12px"><a href="http://www.gnu.org/philosophy/no-word-attachments.html" target="_blank">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div>
</div><div><font color="#0023E9" size="3"><br></font></div></span><br>
</div>
<br></div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>---------------------------------------------------------<br><br>Mohaddeseh Abbasnejad, <br>Room No. 323, Department of Physics, <br>University of Tehran, North Karegar Ave.,<br>
Tehran, P.O. Box: 14395-547- IRAN <br>Tel. No.: <a href="tel:%2B98%2021%206111%208634" value="+982161118634" target="_blank">+98 21 6111 8634</a> & Fax No.: <a href="tel:%2B98%2021%208800%204781" value="+982188004781" target="_blank">+98 21 8800 4781</a><br>
Cellphone: <a href="tel:%2B98%20917%20731%207514" value="+989177317514" target="_blank">+98 917 731 7514</a><br>E-Mail: <a href="mailto:m.abbasnejad@gmail.com" target="_blank">m.abbasnejad@gmail.com</a><br>
Website: <a href="http://physics.ut.ac.ir/" target="_blank">http://physics.ut.ac.ir</a><br><br>---------------------------------------------------------<br><br>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><div>
<span style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px"><font size="3"><span style="font-size:12px">---</span></font></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px">
<font size="3"><span style="font-size:12px">Stefano Baroni - SISSA</span></font><font size="3"><span style="font-size:12px"> </span></font><font size="3"><span style="font-size:12px">&</span></font><font size="3"><span style="font-size:12px"> </span></font><font size="3"><span style="font-size:12px">DEMOCRITOS National Simulation Center - Trieste</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px"><font size="3"><span style="font-size:12px"><a href="http://stefano.baroni.me/" target="_blank">http://stefano.baroni.me</a> <a href="tel:%5B%2B39%5D%20040%203787%20406" value="+390403787406" target="_blank">[+39] 040 3787 406</a> (tel) -528 (fax) / stefanobaroni (skype)</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;min-height:14px;font-size:12px"><br style="font-size:12px"></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;min-height:14px;font-size:12px">
<span style="font-size:14px"><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><font size="3"><span style="font-size:12px">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px"><font size="3"><span style="font-size:12px"><br></span></font></div></span></div><div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px">
<font size="3"><span style="font-size:12px">Please, if possible, don't</span></font><font size="3"><span style="font-size:12px"> </span></font><font size="3"><span style="font-size:12px">send me MS Word or PowerPoint attachments</span></font></div>
<div style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;font-size:14px"><font size="3"><span style="font-size:12px">Why? See:</span></font><font size="3"><span style="font-size:12px"> </span></font><font color="#0023E9" size="3"><span style="font-size:12px"><a href="http://www.gnu.org/philosophy/no-word-attachments.html" target="_blank">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div>
</div><div><font color="#0023E9" size="3"><br></font></div></span><br>
</div>
<br></div></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://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br><br>