Debye modeli

Termodinamik ve katı hal fiziğinde **Debye modeli;&nbsp;**Peter Debye&nbsp;tarafından 1912 yılında geliştirilen, katılarda&nbsp;özgül ısıya (ısı kapasitesi) olan fonon&nbsp;katkısını tahmin etmek için kullanılan metottur.&nbsp;Atomik kristal yapının salınımlarını, bir kutu içerisindeki fononlar gibi düşünerek ele alır.&nbsp;Bu; katıya ayrı ayrı kuantum harmonik osilatörlerden oluşmuş olarak davranan Einstein modelinin tam tersidir. Debye modeli;&nbsp;&nbsp;– **Debye T<sup>3</sup>&nbsp;yasası - **ısı kapasitesini düşük sıcaklıklarda doğru bir şekilde tahmin eder., düşük sıcaklıklarda olan.&nbsp;Tıpkı Einstein modeli gibi,&nbsp;yüksek sıcaklıklarda&nbsp;Dulong–Petit Yasasını da doğru bir şekilde kapsar.&nbsp;Ama, ara sıcaklıklarda basitleştirmek için yapılan varsayımlar nedeniyle doğruluğu kusurludur.

Debye modelinin en doğru kullanımları için&nbsp;M. Sunada ve T. Sunada kaynak olarak kullanılabilir.

## Türetme
Debye modeli; elektromanyetik radyasyona kapalı bir kutudaki fotonlardan oluşmuş gaz olarak davranan&nbsp;Planck'ın siyah cisim ışımasının katı-hal muadilir.&nbsp;Debye modeli, atomik titreşimlere kutudaki (kutu = katı cisim)&nbsp;fononlar olarak yaklaşır. Hesaplamaların çoğu aynıdır.

Bir kenarı&nbsp;uzunluğunda bir küp düşünün. Kutudaki parçacık başlığında, kutunun içindeki (sadece bir eksende hizalanmış) sonik bozuklukların rezonans modlarının dalga boyları:
: 
olarak verilir.&nbsp; bir tam sayıdır.&nbsp;Bir fononun enerjisi ise:
: 
 Planck sabiti&nbsp;ve&nbsp; bu fononun frekansıdır. Eğer frekansın dalgaboyu ile ters orantılı olduğunu varsayarsak, elimizde:
: 
Burada&nbsp; katı içindeki ses hızıdır.
Üç boyutta ele alırsak:
: 
&nbsp;fononun üç boyutta skaler momentum büyüklüğüdür.

Frekansın dalga boyuyla (ses hızını sabit kabul edersek) ters orantılı olduğu varsayımı yüksek-enerjili fononlar (detaylı bilgi için fonon başlığına bakınız) için olmasa da düşük-enerjili fononlar için iyi bir yaklaşımdır. Bu; Debye modelinin sınırlandırmalarından birisidir ve ara sıcaklıklarda sonuçların hatalı olmasına karşılık gelir. Oysaki düşük ve yüksek sıcaklıklarda hatasızlardır.

Şimdi kutudaki toplam enerjiyi hesaplarsak,
: 
Burada&nbsp; kutudaki&nbsp;. Diğer bir deyişle, toplam enerji, enerji ve o enerjiye sahip (tek boyutta) fonon sayısının çarpımının toplamıdır.&nbsp;3 boyutta:
: 
Şimdi, burası&nbsp;**Debye modeli** ve Planck'ın siyah cisim ışıması farklılaştığı yerdir. Kutudaki&nbsp;elektromanyetik radyasyonun aksine, fononun enerji seviyesi sayısı sınırlıdır. Çünkü, bir fonon sınırsız frekansa sahip olamaz. Frekansı, ilerlediği ortamla—katının atomik örgüsü (atomic lattice of the solid) sınırlanmıştır. Aşağıdaki; enine(transverse) bir fononun tasvirine bakınız.
::::: 
Alttaki figürden de anlaşılacağı üzere, bir fononun sahip olabileceği minimum dalga boyunun atom aralığının iki katı olduğunu varsaymak makul görünüyor. Bir katının içerisinde&nbsp; atom var. Bizim katımızın şekli bir küp, bu da demektir ki kenar başına&nbsp;&nbsp;atom var. Öyleyse atom aralığı&nbsp;ve minimum dalga boyu:
: 
Böylelikle; maksimum mod sayısı&nbsp;&nbsp;(foton&nbsp;için sonsuz):
: 
Bu, üçlü enerji toplamının üst sınırıdır.
: 
Yavaşça değişen, sürekli (well-behaved) fonksiyonlar için, toplam integral olarak yazabilir. (Thomas-Fermi yaklaşımı):
: 
Şimdiye kadar,&nbsp;&nbsp;Fononlar Bose–Einstein istatistiklerine uyarlar. Dağılımları ünlü&nbsp;Bose–Einstein formülü ile verilir:
: 
Bir fononun üç olası kutuplaşma durumu (bir boyuna&nbsp;ve yaklaşık olarak enerjisini etkilemeyen iki enine) onun
yukarıdaki enerji) olduğu için yukarıdaki formül 3 ile çarpılır,
: 
(Aslında&nbsp;*efektif ses hızı kullanırsak* &nbsp;Debye sıcaklığı&nbsp; (aşağıya bakınız) ile daha doğrusu 

Bunu da enerji integralinde yerine yazarsak,
: 
Bu integrallerin fotonlar için hesaplanmasının kolay olmasının sebebi, ışığın frekansının, en azından yarı klasik olarak, bağlanmamış olması gerçeğidir. Yukarıdaki resimde görüldüğü gibi, bu fononlar için geçerli değildir. Debye, bu üçlü integrali yaklaşık olarak hesaplayabilmek için küresel koordinatları kullandı.
: 
ve küpü, kürenin sekizde biri olarak varsaydı.
: 
&nbsp;kürenin yarıçapı, küpteki ve kürenin sekizde birindeki parçacık sayının korunumu ile bulunabilir. Küpün hacmi,&nbsp;&nbsp;tane birim-hücre(unit-cell) hacmi toplamına eşit.
: 
Yarıçap için çözüldüğünde:
: 
Doğru integral için bir küreye göre integralin yerine yazılması, modele başka bir yanlışlık kaynağı getirir.

Öyleyse enerji integrali
: 
İntegral için değişken değiştirirsek&nbsp;,
: 
Görünümü basitleştirmek için,&nbsp;**Debye Sıcaklığı** 
: 
Birçok referans Debye sıcaklığını bazı sabitler ve materyale bağlı değişkenler için kısaltma olarak tanımlar.&nbsp;Ancak, aşağıda gösterildiği gibi, &nbsp;kabaca minimum dalga boyu modunun fonon enerjisine eşittir. Böylelikle Debye sıcaklığını, en yüksek frekans modunun (ve dolayısıyla her modun) uyarıldığı sıcaklık olarak yorumlayabiliriz.

Devam edersek, elimizde özgül bir iç enerji var:
: 
&nbsp;: (üçüncü) Debye fonksiyonu

&nbsp;'ye göre türevini aldığımızda, boyutsuz ısı kapasitesini elde ederiz:
: 
Bu formüller, Debye modelini tüm sıcaklıklarda sağlar. Aşağıda verilen daha temel formüller, düşük ve yüksek sıcaklıkların sınırında asimptotik davranışı verir. Daha önce belirtildiği gibi, bu davranış, ara davranışa aksi olarak daha kesindir.&nbsp;Düşük ve yüksek enerjilerdeki kesinliğin temel nedeni, sırasıyla: (i) kesin&nbsp;*dağılım ilişkisi&nbsp;(dispersion relation)*&nbsp; düşük frekanslar ve (ii) kesin&nbsp;*seviye yoğunlukları(density of states)*&nbsp;frekans aralığı başına titreşim sayısı ile ilgili olarak.

## Debye'nın derivasyonu
Aslında, Debye denklemini biraz daha farklı ve daha basit türetti. Sürekli ortamın katı mekaniğini kullanarak, belirli bir değerden daha düşük frekanstaki titreşim seviyesi(state) sayısının aşağıdaki ifadeye asimptotik olduğunu buldu.
: 
 hacim ve esneklik katsayıları&nbsp;ve yoğunluk ile hesaplanan bir çarpan (faktör).&nbsp;Bunu, T sıcaklığında bir harmonik osilatörün beklenen enerjisiyle(expected energy) birleştirirsek ve (Einstein'ın modelinde zaten kullandığı gibi) titreşim frekansları sonsuza giderse enerji:
: 
Bu form düşük sıcaklıklarda doğru olan&nbsp;&nbsp;davranışını verir. düşük sıcaklıklarda doğru olan davranış. Ama Debye N tane atom için&nbsp;Yaptığı varsayıma göre; bir atomik katı, titreşim seviyelerinin frekanslarının spektrumu yukarıdaki kurala uymaya devam etmek zorundaydı. Maksimum frekans&nbsp;:
: 
Debye, bu varsayımın gerçekten doğru olmadığını biliyordu (yüksek frekanslar birbirine varsayılandan daha yakındı), ancak yüksek sıcaklıktaki doğru davranışı garanti eder (Dulong-Petit yasası). Enerji şu şekildedir:
: 

:: 

:: &nbsp;burada&nbsp;.'ye eşit.

:: 

:: 
Debye fonksiyonu&nbsp;olarak isimlendirilen fonksiyondur.

## Başka bir derivasyon
Önce titreşim frekansı dağılımını türetiriz. Aşağıdaki türetme Appendix VI kısmından alınmıştır.&nbsp;Yan uzunlukları . Elastik dalga&nbsp;dalga denklemine uyacak ve&nbsp;düzlem dalgalardan oluşacaktır; dalga vektörü &nbsp;ele alalım.



Dalga denkleminin çözümü:
:: 
ve sınırlayıcı şartları &nbsp;,

 pozitif tam sayılar. ()'yi ()'in içerisinde yerine yazarsak ve&nbsp;dağılım ilişkisi 
:: 
Yukarıdaki denklem, sabit frekans&nbsp; pozitif). Frekansı&nbsp; (çok büyük bir paralel prizma) için elipsin hacmine yaklaşık verilebilir.&nbsp;Dolayısıyla, modların sayısı&nbsp;,&nbsp;&nbsp;frekans aralığında:

öylece tanımlarsak:&nbsp;.

Aşağıdaki derivasyon kaynağından alınmıştır.&nbsp;Titreşim frekansı&nbsp;N tane atom olduğu için, 3N tane kuantum harmonik osilatör (her x-, y-, z-, yönü için 3 tane) vardır ve&nbsp;olayısıyla&nbsp;:

.

k:&nbsp;Boltzmann sabiti ve h:&nbsp;Planck sabiti&nbsp;ve ()'ü ()'te yerine yazdığımızda, elimizde:

bu tanım daha standarttır. . Maxwell-Boltzmann istatistiklerini kullanarak,&nbsp;
:: .


.

 (çünkü &nbsp;tane mod,&nbsp;. Elimizde:
:: 
Yukarıdaki ifadelerden 1/A için bir tanım elde edebiliriz. Bunu () içerisinde yazarsak:
:: 

: 
Frekansa göre integral aldığımızda:
:: 

## Düşük sıcaklık sınırı
Eğer&nbsp;
: 
Bu belirli integral tam olarak hesaplanabilir.
: 
Düşük sıcaklık sınırında, yukarıda bahsedilen Debye modelinin sınırlamaları geçerli değildir ve (fononik) ısı kapasitesi, sıcaklık, elastik katsayılar ve atom başına düşen hacim arasında doğru bir ilişki verir. (bu son miktarlar, Debye sıcaklığı içerisinde bulunur)

## Yüksek sıcaklık sınırı
Debye katısının sıcaklığı&nbsp;. ğer ise
: 

: 
Bu Dulong-Petit yasasıdır ve ısı kapasitesinin daha da yükselmesine neden olan anharmoniği dikkate almıyor olmasına rağmen oldukça doğrudur. Bir iletken veya yarıiletken katı maddenin toplam ısı kapasitesi elektronlardan da ihmal edilemez bir katkı içerebilir.

## Debye vs Einstein

Peki Debye ve Einstein modelleri deneyle ne kadar örtüşüyor? Şaşırtıcı derecede yakın, ancak Debye düşük sıcaklıklarda doğru olmasına rağmen Einstein doğru değildir.

Modeller ne kadar farklı? Bu soruyu cevaplamak için, doğal olarak, aynı eksenlerde ikisini de çizmek gerekli.&nbsp;Tek sorun, çizilemeyecek olması. Hem Einstein modeli hem de Debye modeli, ısı kapasitesi için fonksiyonel bir form sağlar. Bunlar modellerdir ve ölçeksiz bir model olmaz. Ölçek, modelin gerçek dünyadaki eşdeğeri ile ilgilidir. Aşağıda verilen Einstein modelinin ölçeğini görebilirsiniz:
: 
Debye modelinin ölçeği,&nbsp;
: 
Bu da onları aynı eksenler dizisine koymanın mantıklı olmadığını gösterir. Bunlar aynı şeyin iki modeli fakat farklı ölçeklerde. Eğer biri Einstein'ın sıcaklığını şöyle tanımlarsa:
: T_E \ \stackrel{\mathrm{def}}{=}\  {\epsilon\over k}\,,
Öyleyse:
: T_E \ne T_D\,,
İkisini kıyaslayabilmek için, bir oran aramalıyız.
: \frac{T_E}{ T_D} = ?
Einstein katısı&nbsp;**tek**-frekanslı&nbsp;kuantum harmonik oscilatörden oluşur, \epsilon = \hbar\omega = h\nu. Bu frekans, eğer gerçekten var olsaydı, katıdaki sesin hızıyla bağlantılı olurdu. Eğer sesin ilerlemesini birbirine çarpan bir dizi atom olarak düşünürsek, atomik örgü tarafından sürdürülebilir minimum dalga boyuna&nbsp;\lambda_{min}
: \nu = {c_s\over\lambda} = {c_s\sqrt[3]{N}\over 2L} = {c_s\over 2}\sqrt[3]{N\over V}
Böylelikle&nbsp;**Einstein sıcaklığı:**
: T_E = {\epsilon\over k} = {h\nu\over k} = {h c_s\over 2k}\sqrt[3]{N\over V}\,,
Ve aradığımız oran:
: {T_E\over T_D} = \sqrt[3]{\pi\over6}\ = 0.805995977...
Şimdi her iki model de aynı grafik üzerinde çizilebilir. Bu oran, bir 3-boyutlu küre hacminin, onu içeren küpün hacmine oranının küp kökü olduğuna dikkat edelim; bu sadece yukarıdaki enerji integralini hesaplarken Debye tarafından kullanılan düzeltme faktörüdür.

Alternatif olarak, 2 sıcaklığın oranı, tüm osilatörlerin salınacağı Einstein'ın tek frekansının ve Debye'nin maksimum frekansının oranı olarak görülebilir. Einstein'ın tek frekansının Debye modeli için mevcut olan frekansların ortalaması olduğu görülebilir.

## Debye sıcaklık tablosu
Debye modeli tamamen doğru olmasa da, diğer katkıların (oldukça mobil iletim elektronları gibi) ihmal edilebilir olduğu yalıtkan, kristal katıların düşük sıcaklık ısı kapasitesi için iyi bir yaklaşım verir. Metaller için ısıya elektron katkısı** T** ile orantılıdır, bu da, alçak sıcaklıklarda atomik örgü titreşimleri iT^3 tarafından domine edilir.&nbsp;Bu durumda, Debye modeli, yalnızca belirli ısının örgü katkısı için yaklaşık olarak söylenebilir. Aşağıdaki tabloda çeşitli saf elementler için Debye sıcaklıkları listelenmiştir::

Debye modelinin deneysel verilere uyması, genellikle Debye sıcaklığının sıcaklığa bağlı olmasına izin vererek fenomenolojik olarak geliştirilir; örneğin, sıcaklık Mutlak sıfırdan yaklaşık 100 K'ye yükseldiğinde, su buzunun değeri yaklaşık 222 K'den 300 K'ye yükselir.

## Diğer yarı-parçacıklar için kapsama
Diğer bozonik yarı-parçacıklar için, mesela fononlar (kuantize olmuş ses dalgaları) yerine ferromıknatıslardaki magnonlar için benzer sonuçlar kolayca elde edilebilir. Bu durumda, düşük frekanslarda farklı dağılma ilişkilerine (dispersion relation) sahiptir. Örneğin, magnonlar için&nbsp;E(\nu )\propto k^2&nbsp;, fonon olsaydı&nbsp;E(\nu )\propto k (dalga vektörü&nbsp;k=2\pi /\lambda ). Ayrıca seviyelerin yoğunluğu farklılık gösterir.&nbsp;(Örneğin&nbsp;\int g(\nu ){\rm d}\nu \equiv N\,). Sonuç olarak ferromıknatıslarda, ısı kapasitesine magnon katkısı olur.&nbsp;\Delta C_{\,{\rm V|\,magnon}}\,\propto T^{3/2}, bu, yeterince düşük sıcaklıklarda fonon katkısını domine eder, \,\Delta C_{\,{\rm V|\,phonon}}\propto T^3. Buna karşın metallerde, düşük sıcaklıklardaki asıl ısı kapasitesine yapılan katkı&nbsp;\propto T, elektronlardan gelmektedir. Fermiyoniklerdir,&nbsp;,Arnold Sommerfeld'e kadar uzanan farklı metotlarla hesaplanırlar.

## Sıvılar için kapsama
Sıvılar sadece boyuna fononları destekleyebileceği, enine fononları desteklemeyeceği için, uzun bir süre boyunca, fonon teorisinin sıvıların ısı kapasitesini açıklamayacağını düşünülüyordu, katılardaki ısı kapasitesinin 2/3'ünden sorumlu olduklarını düşünürsek haksız bir çıkarım değildi.&nbsp;Ancak, nötronlar ve X-Işınları ile yapılan&nbsp;Brillouin saçılması deneyleri&nbsp;Yakov Frenkel'in sezgilerini destekliyordu.&nbsp;Bununla; her ne kadar Frenkel frekansı denilen, eşik frekansının üzerindeki frekanslarla sınırlı olsa da, enine fononların sıvılarda da mevcut olduğu anlaşıldı. Enerjinin çoğu bu yüksek-frekans modlarında bulunduğu için,&nbsp;basit sıvıların deneysel ısı kapasitelerine iyi bir yakınlık elde etmek için Debye modelinin basit bir modifikasyonu yeterlidir.

## Ayrıca bakınız
* Bose gas
* Debye frequency
* Gas in a box
* Kinetic theory of solids
* Grüneisen parameter

## Başvurular

## Konuyla ilgili yayınlar
* *CRC Kitabının Kimya ve Fizik*, 56 Edition (1975-1976)
* Schroeder, Daniel V. *Giriş için Termal Fizik*. Addison-Wesley, San Francisco (2000). Bölüm 7.5.

## Dış bağlantılar
* Özgül ısı, bir hata arama ve kullanma kuvars termal ve ısı iletkenliği deneysel olarak belirlenmesi.