Create Presentation
Download Presentation

Download Presentation
## 5. Phonons Thermal Properties

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -

**5. Phonons Thermal Properties**• Phonon Heat Capacity • Anharmonic Crystal Interactions • Thermal Conductivity**Phonon Heat Capacity**• Planck Distribution • Normal Mode Enumeration • Density of States in One Dimension • Density of States in Three Dimensions • Debye Model for Density of States • Debye T3 Law • Einstein Model of the Density of States • General Result for D(ω)**where**= Thermal expansivity = isothermal compressibility = 1/ B B = Bulk modulus α = linear (1-D) thermal expansivity**Lattice heat capacity:**p = polarization Planck distribution:**Planck Distribution**System at constant T Canonical ensemble : Boltzmann factor For a set of identical harmonic oscillators Nn = number of oscillators in the nth excited state when system is in thermal equilibrium Probability of an oscillator in the nth excited state: Occupation number:**Density of States in One Dimension**Fixed boundary problem of N+1 particles. N = 10 → with Number of allowed K for non-stationary solutions is N–1 = Number of mobile atoms Polarization p : 1 long, 2 trans**Periodic boundary problem of N particles**N = 8 → with → Number of allowed K for non-stationary solutions is N = Number of mobile atoms**fixed B.C.**→ Periodic B.C. →**Density of States in Three Dimensions**Periodic B.C. ; N 3 cells in cube of side L → density of states in K-space is Number of modes per polarization lying between ωand ω + d ωis → density of states in ω-space is For isotropic materials, →**Debye Model for Density of States**Debye model: v velocity of sound (for a given type of polarization) For a crystal of N primitive cells: → Debye frequency**Thermal (vibrational) energy **Debye integrals: See Ex on Zeta functions, Arfken where θ Debye temperature for each acoustic branch**Debye T3 Law**For low T, xD → : → for each acoustic branch To account for all 3 acoustic branches, we set and so that Good for T < θ /50**Qualitative Explanation of the T3 Law**Of the 3N modes, only a fraction (KT /KD ) 3 = ( T / θ)3 is excited. →**Einstein Model of the Density of States**N oscillators of freq ω0. Diamond Classical statistical mechanics: Dulong-Petit valueCV = 3NkB**General Result for D(ω)**Si Debye solid vg ~ 0 Van Hove Singularities**Anharmonic Crystal Interactions**• Harmonic (Linear) Waves: • Normal modes do not decay. • Normal modes do not interact. • No thermal expansion. • Adiabatic & isothermal elastic constants are equal. • Elastic constants are independent of P and T. • C → constant for T > θ . Deviation from harmonic behavior → Anharmonic effects**Thermal Expansion**1-D anharmonic potential: c = 1 g = .2 f = .05 Boltzmann distribution High T: → Thermal Expansion**Thermal Conductivity**For phonons, JQ = JU . Heat current density: κ = Thermal conductivity coefficient • Key features of kinetic theory (see L.E.Reichl, “A Modern Course of Statistical Mechanics”, §13.4 ): • Quantities not conserved in particle collisions are quickly thermalized to (global) equilibrium values. (e.g., velocity directions & magnitudes ) • Conserved quantities can remain out of global equilibrium (e.g., stay in local equilibrium. They get transported spatially in the presence of a “gradient”. • MFP is determined by collisions that do not conserve the total momenta of particles. Net amount of A(z) transported across the x-y plane at z0in the +z direction per unit area per unit time: Δz = distance above/below plane at which particle suffered last collision. n = particle density, l = mean free path, a = some constant. bA = 1/3 is determined from self diffusion**For heat conduction, we set**where c = heat capacity per particle = sound velocity C = nc = heat capacity**Thermal Resistivity of Phonon Gas**Harmonic phonons: mfp l determined by collisions with boundaries & imperfections. Anharmonic phonons: only U-processes contribute.**Gas:**Elastic collisions. No T required. κ = Gas: No net mass flow. Inelastic collisions with walls sets up T & n gradients. Finite κ. Crystal: N-processes only. κ = . Crystal: U-processes. Finiteκ.**Umklapp Processes**2-D square lattice Normal processes: Umklapp processes: Energy is conserved: Condition: T > θ: all modes excited → no distinction between N- & U- processes → l 1/T. T < θ: probability of U- processes & hence l1 exp(–θ /2T).**Imperfections**Low T → umklapp processes negligible. Geometric effect dominates. Size effect: l > D = smallest dimension of specimen. Dielectric crystals can have thermal conductivities comparable to those of metals. Sapphire (Al2O3): κ ~ 200W cm–1 K–1 at 30K. Cu: max κ ~ 100W cm–1 K–1. Metalic Ga: κ ~ 845W cm–1 K–1 at 1.8K. ( Electronic contributions dominate in metals. ) Highly purified c-NaF**Isotope effect on Ge.**Enriched: 96% Ge74. Normal: 20% Ge70 , 8% Ge73 , 37% Ge74 , 8% Ge76.