In der Festkörperphysik beschreibt das Einstein-Modell (nach Albert Einstein) eine Methode, um den Beitrag der Gitterschwingungen (Phononen) zur Wärmekapazität eines kristallinen Festkörpers zu berechnen. Da sich das Einstein-Modell ausschließlich auf optische Phononen anwenden lässt, ist es nicht so erfolgreich wie das Debye-Modell, das auch akustische Phononen beschreibt.
Die Gitterschwingungen des Kristalls werden gequantelt, d. h. der Festkörper kann Schwingungsenergie nur in diskreten Quanten $ \hbar \cdot \omega _{E} $ aufnehmen. Diese Quanten nennt man auch Phononen. Man beschreibt den Festkörper dann als aus N quantenharmonischen Oszillatoren bestehend, die jeweils in drei Richtungen unabhängig schwingen können. Die Besetzungswahrscheinlichkeit $ \langle n\rangle $ einer solchen Schwingungsmode (eines Phonons) hängt von der Temperatur T ab und folgt (da Phononen Bosonen sind) der Bose-Einstein-Verteilung:
mit
Damit ergibt sich die innere Energie U im Festkörper zu (Es wurde die Quantisierungsbedingung des harmonischen Oszillators verwendet):
mit
Der Beitrag $ {\frac {\hbar \cdot \omega _{E}}{2}} $ gibt die Nullpunktenergie an.
Der Beitrag der Phononen zur Wärmekapazität ist dann:
mit
Mit der Einstein-Temperatur $ \Theta _{E}={\frac {\hbar \cdot \omega _{E}}{k_{B}}} $ ergibt sich eine einfachere Schreibweise:
Wie das Debye-Modell liefert das Einstein-Modell das korrekte Hochtemperaturlimit nach dem Dulong-Petit-Gesetz:
Im Limes kleiner Temperaturen ergibt sich:
Dieser Verlauf von CV(T) für kleine Temperaturen weicht allerdings erheblich von Messungen ab. Dies hängt mit der Annahme zusammen, alle harmonischen Oszillatoren im Festkörper würden mit einer einheitlichen Frequenz schwingen. Die Verhältnisse im realen Festkörper sind jedoch deutlich komplizierter.