This image presents the dispersion relation of phonon polaritons according to a simple theory^[1], giving the two following dispersion relations for the coupled optical phonon/electromagnetic field system :

1. longitudinal optical phonons (LO phonons)
   ${\displaystyle \omega ^{2}=\omega _{0}^{2}+{\frac {\Omega ^{2}}{1+\chi _{\infty }}}={\text{const}};}$
2. polaritons (transverse em field + TO phonons)
   ${\displaystyle \mathbf {k} ^{2}={\frac {\omega ^{2}}{c^{2}}}\left[1+\chi _{\infty }+{\frac {\Omega ^{2}}{\omega _{0}^{2}-\omega ^{2}}}\right].}$

The dashed red curves correspond to the uncoupled modes (photons and phonons: Ω = 0), and the black curves represent solutions of the above equations with parameters chosen to fit the dispersion relation of polaritons in GaP^[2], giving

Uncoupled optical phonon frequency

ω₀ = 45.4 meV

Coupling parameter

Ω = 61 meV

Electric susceptibility

χ_∞ = 7.5

Note that k is very small compared to π/a.

1. ↑ Max Born, Kun Huang, Dynamical Theory of Crystal Lattices, Oxford University Press, 1954 (ISBN 9780198503699), « Chapter 8 : Infrared Dispersion and the Retardation Effect on Lattice Vibrations »
2. ↑ C. H. Henry and J. J. Hopfield, Raman Scattering by Polaritons, Phys. Rev. Lett. 15, 964 - 966 (1965).

(further minor editing and conversion to svg done with inkscape)

描述Phonon polaritons.svg	English: Dispersion relation of phonon polaritons with parameters corresponding to GaP. Français : Relation de dispersion des polaritons avec des paramètres correspondant au GaP.
日期	2008年11月15日
來源	自己的作品
作者	Mathieu Perrin
SVG開發 InfoField	本SVG檔案的原始碼通過W3C驗證。   本vector image使用MATLAB創作。
原始碼 InfoField	MATLAB code e=1.602e-19; hbar = 1.05457168e-34; c=2.9978e8; ang=1e-10; omega0=45.4e-3*e/hbar; %conversion meV->rd/s Om=61e-3*e/hbar; chi=7.5; k=linspace(0,0.0004/ang)'; y1=omega0*ones(size(k)); %uncoupled phonon dispersion y2=c*k/sqrt(1+chi); %uncoupled photon dispersion y3=sqrt(y1.^2+Om^2/(1+chi));%longitudinal phonons %solve for omega in %k^2*c^2=omega^2*(1+chi+Om^2/(omega0^2-omega^2)) a1 = (1+chi)*ones(size(k)); a2 = -((1+chi)*omega0^2+k.^2*c^2+Om^2); a3 = k.^2*c^2*omega0^2; D = sqrt(a2.^2-4*a1.*a3); y4=sqrt((-a2+D)./(2*a1));%upper polaritons y5=sqrt((-a2-D)./(2*a1));%lower polaritons plot(k*ang, [y1 y2]*1000*hbar/e,'r--','Linewidth',1.8); hold on plot(k*ang, [y3 y4 y5]*1000*hbar/e,'k','Linewidth',1.8); hold off axis([0 0.0004 35 60]); set(gca,'FontSize',18,'TickLength',[.02 .02],'LineWidth',1.5); xlabel('k (A^{-1})'); %A to be replaced ylabel('h\omega (meV)'); %h to be replaced

我，本作品的著作權持有者，決定用以下授權條款發佈本作品：

已授權您依據自由軟體基金會發行的無固定段落、封面文字和封底文字GNU自由文件授權條款1.2版或任意後續版本，對本檔進行複製、傳播和／或修改。該協議的副本列在GNU自由文件授權條款中。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

此檔案採用創用CC 姓名標示-相同方式分享 3.0 未在地化版本、2.5 通用版、2.0 通用版以及1.0 通用版授權條款。

您可以自由：

• 分享 – 複製、發佈和傳播本作品
• 重新修改 – 創作演繹作品

惟需遵照下列條件：

• 姓名標示 – 您必須指名出正確的製作者，和提供授權條款的連結，以及表示是否有對內容上做出變更。您可以用任何合理的方式來行動，但不得以任何方式表明授權條款是對您許可或是由您所使用。
• 相同方式分享 – 如果您利用本素材進行再混合、轉換或創作，您必須基於如同原先的相同或兼容的條款，來分布您的貢獻成品。

https://creativecommons.org/licenses/by-sa/3.0CC BY-SA 3.0 Creative Commons Attribution-Share Alike 3.0 truetrue

您可以選擇您需要的授權條款。