[Pw_forum] phonon eigenvectors
Eric Abel
etabel at hotmail.com
Thu Aug 31 15:22:21 CEST 2006
>From: Stefano Baroni <baroni at sissa.it>
>Reply-To: pw_forum at pwscf.org
>To: pw_forum at pwscf.org
>Subject: Re: [Pw_forum] phonon eigenvectors
>Date: Thu, 31 Aug 2006 10:05:04 +0200
>
>Eric: do not be satisfied by meaningful answer to your meaningful and
>legitimate questions. You have to *understand* them (i.e. do not stop
>thinking and striving to understand, until answers given by others sound
>to you as YOUR answers. Keep asking yourself questions until you are sure
>that you can answer all of them).
>
>Signed: Grandad
>
>Coming now to your present doubt. You have to understand here that what is
>arbitrary here is the overall _phase_ of the eigenvectors. This obviously
>implies (in some sense) that the imaginary part is arbitrary, but *NOT*
>that it can (always) be neglected, as this specic case shows. What has to
>be nonzero is NOT the real part of the eigenvector, but its modulus.
Stefano, I apologize if it sounded as though I was trying to trivialize your
(previous) answer. I read most of the postings to this users forum, and
understand there are a great many who prefer to simply ask questions rather
than seek out the answers themselves. As for me, I am too proud to ask for
help, so when I do, it's because I feel I have no where else to turn. That
said, I probably didn't ask the question in the proper way. My bad. My gut
instinct was that the modulus, rather than the real part, was the important
value...I just wanted confirmation. The questions which I asked myself
before posting to the forum were the following:
Question: What is the meaning of an eigenvector?
Answer: An eigenvector tells how the atoms are displaced in the vibration.
Question: What is the imaginary part?
Answer: The imaginary part is a phase factor. If one atom has a significant
imaginary part with respect to the rest, then then it's displacement in the
vibration will be phase shifted with respect to the others. In the case of
a completely imaginary eigenvector, the displacement will be completely out
of phase.
Question: If the imaginary part is a phase factor, then what does it mean
if all of the components are imaginary?
Answer: First answer: nothing. If all of the atoms are "out of phase",
then they are "in phase" with respect to eachother, therefore, having a
completely imaginary or completely real eigenvector should be equivalent.
Question: Great. If this is the case, then why does this eigenvector come
up imaginary when all of the other vectors come up real?
Answer: Hmmm. That's a good one. There must be some reason that the code
chooses an imaginary eigenvector for this mode...time to get help on this
one. We are fresh out of answers.
And with this I come first to the pwscf discussion archive. I didn't find
any discussion of imaginary eigenvectors with real eigenvalues. I looked to
the literature, with the same result. There comes a time when one asks
himself questions to which he doesn't have the answers, then it comes time
to discuss with his peers. With that I come to you with my question:
Why did the code "choose" to make this displacement imaginary. Is this
simply an artifact of the matrix diagonalization, or is there some physical
implications to this?
Thank you for your time and patience.
Eric
>
>Is it clearer now? Could you answer similar questions by yourself in the
>future?
>
>Signed: Stefano B.
>
>
>On Aug 31, 2006, at 6:07 AM, Eric Abel wrote:
>
>>Hello PWSCF users,
>>
>>I have a question which is sort of a follow-up to a similar question
>>which I asked previously. Anyway, I was wondering about the imaginary
>>part of the phonon eigenvectors, which at the time, Stefano had informed
>>me are more or less arbitrary. At the time that made sense. But now I
>>am trying to make sense of the following eigenvectors:
>>
>>q = 0.0000 0.0563 0.0000
>>**********************************************************************
>>****
>> omega( 1) = -0.308036 [THz] = -10.275037 [cm-1]
>>( 0.000000 0.000000 0.000000 -0.142674 -0.338669 0.000000
>>)
>>( 0.000000 0.000000 0.031918 -0.201463 -0.388896 -0.061613
>>)
>>( 0.000000 0.000000 0.019612 -0.123787 -0.367551 -0.058231
>>)
>>( 0.000000 0.000000 0.000000 -0.226888 -0.354765 0.000000
>>)
>>( 0.000000 0.000000 0.036445 -0.230039 -0.328252 -0.052005
>>)
>>( 0.000000 0.000000 0.000000 -0.096682 -0.403396 0.000000
>>)
>> omega( 2) = 0.617574 [THz] = 20.600196 [cm-1]
>>( -0.403925 -0.012145 0.000000 0.000000 0.000000 0.000000
>>)
>>( -0.397050 -0.075201 0.000000 0.000000 0.000000 0.000000
>>)
>>( -0.397199 -0.075229 0.000000 0.000000 0.000000 0.000000
>>)
>>( -0.404073 -0.012149 0.000000 0.000000 0.000000 0.000000
>>)
>>( -0.408949 -0.077455 0.000000 0.000000 0.000000 0.000000
>>)
>>( -0.416112 -0.012511 0.000000 0.000000 0.000000 0.000000
>>)
>> omega( 3) = 0.788199 [THz] = 26.291668 [cm-1]
>>( 0.000000 0.000000 0.000000 0.389819 -0.233273 0.000000
>>)
>>( 0.000000 0.000000 -0.055122 0.347926 -0.110454 -0.017499
>>)
>>( 0.000000 0.000000 -0.061881 0.390589 -0.148121 -0.023467
>>)
>>( 0.000000 0.000000 0.000000 0.342277 -0.187111 0.000000
>>)
>>( 0.000000 0.000000 -0.050913 0.321357 -0.249835 -0.039581
>>)
>>( 0.000000 0.000000 0.000000 0.384275 -0.092029 0.000000
>>)
>>
>>These are the first 3 (accoustic) modes of the spectrum. If the
>>imaginary part doesn't matter, then it would appear that the accoustic
>>displacement along the b-direction is practically zero. However, the
>>imaginary part is the same magnitude as the displacements in the other
>>orthogonal directions. Is there any physical meaning to this, or is it
>>just an artifact of the calculation?
>>
>>Eric,
>>Ph.D. Student, M.I.T.
>>
>>
>>_______________________________________________
>>Pw_forum mailing list
>>Pw_forum at pwscf.org
>>http://www.democritos.it/mailman/listinfo/pw_forum
>
>---
>Stefano Baroni - SISSA & DEMOCRITOS National Simulation Center - Trieste
>[+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)
>
>Please, if possible, don't send me MS Word or PowerPoint attachments
>Why? See: http://www.gnu.org/philosophy/no-word-attachments.html
>
>
>
More information about the Pw_forum
mailing list