Solving Phonon Excitation in Diatomic Lattices

  • Thread starter Thread starter aabb009
  • Start date Start date
  • Tags Tags
    Excitation Phonon
aabb009
Messages
1
Reaction score
0

Homework Statement


Ok, I need to show that in an acoustic mode of vibration in a diatomic lattice, for small k, \omega \propto k, and find the constant of proportionality.

Homework Equations


A_1\left(\omega^2M-\frac{2T}{a}\right)+A_2\left(\frac{2T}{a}cos(ka)\right)=0
, and
A_1\left(\frac{2T}{a}cos(ka)\right)+A_2\left(\omega^2m-\frac{2T}{a}\right)=0
hence:
\omega^2 = \frac{T}{a}\left[\frac{1}{M} + \frac{1}{m}\right] - \frac{T}{a}\left[\left(\frac{1}{M}+\frac{1}{m}\right)^2-\frac{4sin^2(ka)}{Mm}\right]^{1/2}


The Attempt at a Solution


I work through it, but repeatedly find that \omega^2 \propto k, and I can't see anyway of getting a k^2 factor on the right.
\omega^2 = \frac{T}{a}\left[\left(\frac{1}{M} + \frac{1}{m}\right) - \sqrt{\left(\frac{1}{M}+\frac{1}{m}\right)^2-\frac{4sin^2(ka)}{Mm}}\right]
\omega^2 = \frac{T}{a}\left[\left(\frac{1}{M} + \frac{1}{m}\right) - \sqrt{\left(\frac{1}{M}+\frac{1}{m}\right)^2\left[1-\frac{Mm4sin^2(ka)}{(M+m)^2}\right]}\right]
\omega^2 = \frac{T}{a}\left(\frac{M+m}{Mm}\right)\left[1 - \sqrt{\left[1-\frac{Mm4sin^2(ka)}{(M+m)^2}\right]}\right]
with small angle approximation we get:
\omega^2 = \frac{T}{a}\left(\frac{M+m}{Mm}\right)\left[1 - \sqrt{1+\frac{4Mmk^2a^2}{(M+m)^2}}\right]
\omega^2 = \frac{T}{a}\left(\frac{M+m}{Mm}\right)\left(1-1+\frac{2\sqrt{Mm}ka}{m+M}\right)
hence
\omega^2 = \frac{2T}{\sqrt{Mm}}k

Where am I going wrong? I don't see any way to prove this.
 
Physics news on Phys.org
aabb009 said:

Homework Statement


Ok, I need to show that in an acoustic mode of vibration in a diatomic lattice, for small k, \omega \propto k, and find the constant of proportionality.

Homework Equations


A_1\left(\omega^2M-\frac{2T}{a}\right)+A_2\left(\frac{2T}{a}cos(ka)\right)=0
, and
A_1\left(\frac{2T}{a}cos(ka)\right)+A_2\left(\omega^2m-\frac{2T}{a}\right)=0
hence:
\omega^2 = \frac{T}{a}\left[\frac{1}{M} + \frac{1}{m}\right] - \frac{T}{a}\left[\left(\frac{1}{M}+\frac{1}{m}\right)^2-\frac{4sin^2(ka)}{Mm}\right]^{1/2}

The Attempt at a Solution


I work through it, but repeatedly find that \omega^2 \propto k, and I can't see anyway of getting a k^2 factor on the right.
\omega^2 = \frac{T}{a}\left[\left(\frac{1}{M} + \frac{1}{m}\right) - \sqrt{\left(\frac{1}{M}+\frac{1}{m}\right)^2-\frac{4sin^2(ka)}{Mm}}\right]
\omega^2 = \frac{T}{a}\left[\left(\frac{1}{M} + \frac{1}{m}\right) - \sqrt{\left(\frac{1}{M}+\frac{1}{m}\right)^2\left[1-\frac{Mm4sin^2(ka)}{(M+m)^2}\right]}\right]
\omega^2 = \frac{T}{a}\left(\frac{M+m}{Mm}\right)\left[1 - \sqrt{\left[1-\frac{Mm4sin^2(ka)}{(M+m)^2}\right]}\right]
the next step is wrong, you changed a minus into a plus between the terms in the sqrt
with small angle approximation we get:
\omega^2 = \frac{T}{a}\left(\frac{M+m}{Mm}\right)\left[1 - \sqrt{1+\frac{4Mmk^2a^2}{(M+m)^2}}\right]
the next step is wrong. you didn't use the right expansion of sqrt(1-x). use
<br /> \sqrt(1-x)\approx 1-x/2<br />
\omega^2 = \frac{T}{a}\left(\frac{M+m}{Mm}\right)\left(1-1+\frac{2\sqrt{Mm}ka}{m+M}\right)
hence
\omega^2 = \frac{2T}{\sqrt{Mm}}k

Where am I going wrong? I don't see any way to prove this.
 
Last edited:
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top