Στερεό του Αϊνστάιν

Το **στερεό του Αϊνστάιν** είναι ένα θεωρητικό μοντέλο στερεών που προτάθηκε για πρώτη φορά από τον Άλμπερτ Αϊνστάιν το 1907, με σκοπό τη μελέτη της εξάρτησης της θερμοχωρητικότητας των στερεών από τη θερμοκρασία. 

Η σημασία του (απλουστευμένου) αυτού μοντέλου είναι μεγάλη, καθώς είναι σε θέση να εξηγήσει, σε αδρές γραμμές, τα πειραματικά δεδομένα που ήταν γνωστά εκείνη την εποχή όσον αφορά την εξάρτηση της θερμοχωρητικότητας των στερεών από τη θερμοκρασία, κάτι που δεν είναι δυνατόν να γίνει με τους νόμους της Κλασικής Φυσικής.

Το μοντέλο του στερεού του Αϊνστάιν είναι κατ' ουσίαν η εφαρμογή των νόμων της στατιστικής μηχανικής για κανονικές κατανομές. Το μοντέλο αυτό βασίζεται στις εξής υποθέσεις:

* Το στερεό αποτελείται από Ν *ταυτόσημους*, *μη αλληλεπιδρόντες* κβαντικούς αρμονικούς ταλαντωτές (ΚΑΤ) που βρίσκονται σε ισορροπία με περιβάλλον θερμοκρασίας Τ
* Όλοι οι ταλαντωτές έχουν την ίδια κυκλική συχνότητα ω

Αν και αμφότερες υποθέσεις σε ρεαλιστικά συστήματα είναι εσφαλμένες, μπορεί κανείς να αντλήσει σημαντικά συμπεράσματα που είναι σε θέση να εξηγήσουν χονδρικά την εμπειρικά γνωστή εξάρτηση της θερμοχωρητικότητας των στερεών από τη θερμοκρασία.

## Πειραματικά δεδομένα

Το βασικό εμπειρικό δεδομένο της εποχής για πολλά στερεά ήταν ότι η θερμοχωρητικότητά τους είναι περίπου ίση με **3R** ανά γραμμομόριο ουσίας, όπου R η παγκόσμια σταθερά των αερίων (νόμος των Dulong-Petit). Αυτό βρισκόταν σε συμφωνία με την Κλασική Φυσική, σύμφωνα με την οποία η ειδική θερμότητα πρέπει να είναι ανεξάρτητη της θερμοκρασίας. Ανεξάρτητα πειραματικά δεδομένα έδειξαν όμως ότι η θερμοχωρητικότητα μεταβάλλεται σύμφωνα με τη θερμοκρασία, πλησιάζοντας το μηδέν καθώς η θερμοκρασία τείνει στο απόλυτο μηδέν.

## Η θεωρητική εξήγηση

Έστω σύστημα Ν, μη αλληλεπιδρώντων, ΚΑΤ συχνότητας ω το καθένα, σε θερμοδυναμική ισορροπία με περιβάλλον θερμοκρασίας Τ. Η ενέργεια κάθε ΚΑΤ ισούται με

:  \epsilon_n=\left(n+\frac{1}{2}\right)\hbar\omega \ , 

όπου n=0,1,2,... ο κβαντικός αριθμός που χαρακτηρίζει την κατάσταση διέγερσης στην οποία βρίσκεται ο ταλαντωτής. Αν ε1,ε2,...,εN η ενέργεια του υπ' αριθμόν 1,2,...,Ν σωματιδίου με αντίστοιχο κβαντικό αριθμό n1,n2,...,nN, τότε η ενέργεια μιας *τυχαίας* μικροκατάστασης r του συστήματος θα ισούται με το *άθροισμα* των επιμέρους ενεργειών του κάθε ταλαντωτή. Συνεπώς,

:  E_r=\left(n_1+\frac{1}{2}\right)\hbar\omega+\left(n_2+\frac{1}{2}\right)\hbar\omega+...+\left(n_N+\frac{1}{2}\right)\hbar\omega 

Για παράδειγμα, μία από τις μικροκαταστάσεις του συστήματος θα μπορούσε να αντιστοιχεί σε ταλαντωτές που βρίσκονται όλοι στη θεμελιώδη τους κατάσταση (n1,n2,...,nN=0). Στη περίπτωση αυτή, η ενέργεια θα ισούται προφανώς με Νħω/2.

Εν γένει, μπορούμε να πούμε ότι η τυχαία μικροκατάσταση r του συστήματος θα καθορίζεται *πλήρως* αν γνωρίζουμε ακριβώς σε ποια κατάσταση βρίσκεται ο κάθε ταλαντωτής. Δηλαδή, η τυχαία μικροκατάσταση r αντιστοιχεί σε μία συγκεκριμένη αλληλουχία κβαντικών αριθμών n1,n2,...,nN. Μαθηματικά,

:  r=\{n_1,n_2,...,n_N\} \ \ \ 

Για να υπολογίσουμε τις διάφορες θερμοδυναμικές ποσότητες του συστήματος, θα χρειαστεί να υπολογίσουμε τη συνάρτηση επιμερισμού

:  Z=\sum_r e^{-\beta E_r}=\sum_r \exp \left\lbrace -\beta\left[\left(n_1+\frac{1}{2}\right)\hbar\omega+\left(n_2+\frac{1}{2}\right)\hbar\omega+...+\left(n_N+\frac{1}{2}\right)\hbar\omega\right] \right\rbrace 

Ο υπολογισμός της παραπάνω ποσότητας εκ πρώτης όψεως είναι, αν όχι πολύ δύσκολη, πρακτικά ακατόρθωτη καθώς υπάρχουν άπειρες μικροκαταστάσεις για το σύστημά μας (αφού ο κβαντικός αριθμός n κάθε ταλαντωτή μπορεί να πάρει τιμές από το μηδέν έως το άπειρο). Παρ' όλα αυτά, μπορούμε να κάνουμε χρήση της υπόθεσης που κάναμε ότι οι ταλαντωτές δεν αλληλεπιδρούν μεταξύ τους και ότι ταλαντώνονται με την ίδια συχνότητα για να υπολογίσουμε τη συνάρτηση επιμερισμού σε κλειστή μορφή.

Ας θεωρήσουμε αρχικά την απλουστευμένη περίπτωση των 2 ταλαντωτών (Ν=2). Οι δυνατές μικροκαταστάσεις ενός τέτοιου συστήματος είναι οι εξής:

:  r={(n_1,n_2)}=\{ (0,0),(0,1),...,(0,\infty),(1,0),(1,1),...,(1,\infty),...,(\infty,0),(\infty,1),...,(\infty,\infty) \} 

Για να αθροίσουμε πάνω σε όλες τις καταστάσεις λοιπόν,

:  \sum_r(n_1,n_2)=\sum_{n_2=0}^{\infty}(0,n_2)+\sum_{n_2=0}^{\infty}(1,n_2)+...+\sum_{n_2=0}^{\infty}(\infty,n_2)=\sum_{n_1=0}^{\infty}\left(\sum_{n_2=0}^{\infty}(n_1,n_2)\right) 

Συνεπώς,

:  \sum_r \ \xrightarrow \ \sum_{n_1=0}^{\infty}\sum_{n_2=0}^{\infty} 

Αντίστοιχα για Ν ταλαντωτές,

:  \sum_r \ \xrightarrow \ \sum_{n_1=0}^{\infty}\sum_{n_2=0}^{\infty}\cdot\cdot\cdot\sum_{n_N=0}^{\infty} 

Άρα λοιπόν, η συνάρτηση επιμερισμού θα ισούται με

:  Z=\sum_{n_1=0}^{\infty}\sum_{n_2=0}^{\infty}\cdot\cdot\cdot\sum_{n_N=0}^{\infty} e^{-\beta \left[ \left(n_1+\frac{1}{2}\right)\hbar\omega+\left(n_2+\frac{1}{2}\right)\hbar\omega+...+\left(n_N+\frac{1}{2}\right)\hbar\omega \right]}=\left(\sum_{n_1=0}^{\infty}e^{-\beta\left(n_1+\frac{1}{2}\right)\hbar\omega}\right)\left(\sum_{n_2=0}^{\infty}e^{-\beta\left(n_2+\frac{1}{2}\right)\hbar\omega}\right)\cdot\cdot\cdot\left(\sum_{n_N=0}^{\infty}e^{-\beta\left(n_N+\frac{1}{2}\right)\hbar\omega}\right) 

Όμως, οι δείκτες άθροισης n1,n2,...,nN είναι βωβοί. Αυτό σημαίνει ότι όλα τα αθροίσματα της παραπάνω σχέσης είναι πανομοιότυπα. Θέτοντας λοιπόν n1=n2=...=nN=n, βρίσκουμε ότι

:  Z=\left(\sum_{n=0}^{\infty}e^{-\beta\left(n+\frac{1}{2}\right)\hbar\omega}\right)^N 

Ο υπολογισμός της συνάρτησης επιμερισμού του συστήματος ανάγεται λοιπόν στον υπολογισμό ενός μόνο αθροίσματος. Έχουμε:

:  \sum_{n=0}^{\infty}e^{-\beta\left(n+\frac{1}{2}\right)\hbar\omega}=e^{-\beta\hbar\omega/2}\sum_{n=0}^{\infty}\left(e^{-\beta\hbar\omega}\right)^n 

Για 0 \sum_{n=0}^{\infty}x^n=\frac{1}{1-x} 

συνεπώς

:  \sum_{n=0}^{\infty}e^{-\beta\left(n+\frac{1}{2}\right)\hbar\omega}=\frac{e^{-\beta\hbar\omega/2}}{1-e^{-\beta\hbar\omega}} 

Η συνάρτηση επιμερισμού όμως ισούται με την παραπάνω ποσότητα υψωμένη στη δύναμη N:

:  Z=\frac{e^{-\frac{N}{2}\beta\hbar\omega}}{\left(1-e^{-\beta\hbar\omega}\right)^N}\Leftrightarrow \ln{Z}=-\frac{N}{2}\beta\hbar\omega-N\ln\left(1-e^{-\beta\hbar\omega}\right) 

Η εσωτερική ενέργεια, U, του συστήματος ισούται με

:  U=-\frac{\partial \ln{Z}}{\partial\beta}=N\hbar\omega\left(\frac{1}{2}+\frac{1}{e^{\beta\hbar\omega}-1}\right) 

Είναι εύκολο να δειχθεί ότι η θερμοχωρητικότητα του συστήματος υπό σταθερό όγκο θα είναι

:  C_V=\frac{\partial U}{\partial T}=Nk_B(\beta\hbar\omega)^2\frac{e^{\beta\hbar\omega}}{(e^{\beta\hbar\omega}-1)^2} 

Η παραπάνω έκφραση απλοποιείται σημαντικά αν θεωρήσουμε ότι kBT>>ħω (βħω C_V\approx Nk_B(1+\beta\hbar\omega)\approx Nk_B 

όπου αγνοήσαμε όρους της τάξης του βħω. Μέχρι στιγμής, θεωρήσαμε ότι οι ταλαντωτές του στερεού είναι μονοδιάστατοι. Τα ίδια επιχειρήματα ακριβώς ισχύουν και στην τριδιάστατη περίπτωση, με την προφανή διαφορά ότι θα πρέπει να πολλαπλασιάσουμε το αποτέλεσμά μας με τον αριθμό των βαθμών ελευθερίας του συστήματος, δηλαδή με το 3. Η θεωρητική πρόβλεψη του μοντέλου του στερεού του Αϊνστάιν είναι λοιπόν ότι

:  C_V(k_BT\gg\hbar\omega)\approx 3Nk_B=3nR 

Που είναι ακριβώς το γνωστό εμπειρικό αποτέλεσμα στο οποίο κατέληξαν οι Dulong-Petit για τη θερμοχωρητικότητα στερεών.

## Ρεαλιστικότερα μοντέλα

Αν και η ακριβής μοντελοποίηση της εξάρτησης της θερμοχωρητικότητας των διαφόρων στερεών από τη θερμοκρασία είναι ένα περίπλοκο ζήτημα που απαιτεί εκτενή χρήση υπολογιστικών τεχνικών (οι λεπτομέρειες των οποίων εξαρτώνται από τη χημική σύσταση του κάθε στερεού), υπάρχουν ορισμένα ακριβέστερα αναλυτικά θεωρητικά μοντέλα από το εκείνο του Αϊνστάιν. Ένα τέτοιο μοντέλο είναι εκείνο του Debye, το οποίο είναι σε θέση να εξηγήσει τη συμπεριφορά της θερμοχωρητικότητας σε χαμηλές θερμοκρασίες (σε αντίθεση με εκείνο του Αϊνστάιν).

## Αναφορές


## Βιβλιογραφία

* 
* Mandl, F. *«Statistical Physics»*. Chichester: John Wiley & Sons Ltd.
* Οικονόμου, Ε.Ν. (2002). *«Στατιστική Φυσική & Θερμοδυναμική»*. Πανεπιστημιακές Εκδόσεις Κρήτης, Ηράκλειο.

## Πρόσθετα
* Σημειώσεις Στατιστικής Μηχανικής του καθηγητή Eric Poisson του Πανεπιστημίου Guelph στον Καναδά σε μορφή PDF https://web.archive.org/web/20111225232318/http://www.physics.uoguelph.ca/poisson/research/spii.pdf

Κατηγορία:Στατιστική μηχανική
Κατηγορία:Φυσική συμπυκνωμένης ύλης