Statisztikus termodinamikai alapok ismtlse A rszletes anyag elrhet

  • Slides: 73
Download presentation
Statisztikus termodinamikai alapok ismétlése A részletes anyag elérhető: http: //keszei. chem. elte. hu/Stat. Th/

Statisztikus termodinamikai alapok ismétlése A részletes anyag elérhető: http: //keszei. chem. elte. hu/Stat. Th/

Tankönyv – Ernő Keszei: Chemical Thermodynamics – An Introduction, Springer, 2011 (Magyar változata: http:

Tankönyv – Ernő Keszei: Chemical Thermodynamics – An Introduction, Springer, 2011 (Magyar változata: http: //keszei. chem. elte. hu/fizkem 1/Tankonyv. pdf) Statisztikus termodinamikáról a 10. fejezet szól Források – Szépfalusi Péter, Tél Tamás: Statisztikus fizika, Egyetemi jegyzet, ELTE Elméleti Fizikai Tanszék, 1976 – Herbert B. Callen: Thermodynamics and an Introduction to Thermostatistics, 2 nd Edition, John Wiley & Sons, 1985 – W. Greiner, L. Neise, H. Stöcker: Thermodynamics and Statistical Mechanics, 2 nd Edition, Springer, 1994 – P. W. Atkins: Physical Chemistry, 6 th Edition, Oxford University Press, 1998 – R. A. Alberty, R. J. Silbey: Physical Chemistry, John Wiley & Sons, 1992

Bevezetés 1 A TERMODINAMIKA alkalmas kémiailag érdekes anyagi rendszerek egyensúlyainak leírására. Összefüggései általánosak. Konkrét

Bevezetés 1 A TERMODINAMIKA alkalmas kémiailag érdekes anyagi rendszerek egyensúlyainak leírására. Összefüggései általánosak. Konkrét anyagi rendszerek esetén viszont mindig szükség van az egyedi fundamentális egyenletre, állapotegyenletekre, vagy táblázatokra. A használt alapmennyiségek (U, S, V, ni , T, P, stb. . . ) között csak egy van, amelyet újként kellett a termodinamikában megtanulni: az entrópia. (A többi új mennyiség ebből már származtatható. ) A termodinamikából azt is tudjuk, hogy ha ismerjük egy adott rendszer fundamentális egyenletét, akkor összes lehetséges állapotát ismerjük. A fundamentális egyenlet alakja lehet pl. S = S (U, V, n ) A STATISZTIKUS TERMODINAMIKA fő célja: a rendszert alkotó részecskék (molekulák) egyszerű fizikai tulajdonságaiból meghatározni a fundamentális egyenletet. Ezzel egyúttal közelebbről meghatározzuk az entrópiát is: Értelmezzük az entrópiát a molekulák sokaságának tulajdonságai alapján.

Bevezetés 2 A STATISZTIKUS TERMODINAMIKAI LEÍRÁST (mint ahogy a klasszikus termodinamikában is tettük) egy

Bevezetés 2 A STATISZTIKUS TERMODINAMIKAI LEÍRÁST (mint ahogy a klasszikus termodinamikában is tettük) egy EGYSZERŰ, IZOLÁLT rendszerrel kezdjük. EGYSZERŰ rendszer: makroszkopikusan homogén, izotróp, (elektromosan) töltetlen, (kémiailag) inert, (térfogatához képest) kis felületű, elektromos, mágneses és gravitációs tértől mentes. IZOLÁLT rendszer: U, V, és n állandó. Az EGYSZERŰ, IZOLÁLT rendszer legyen EGYKOMPONENSŰ is. Az n mólszám (ez persze a IUPAC kedvéért anyagmennyiség) helyett használjuk a molekulák (részecskék) N számát. Egy ilyen rendszer modelljét nevezzük mikrokanonikus sokaságnak.

MIKROKANONIKUS SOKASÁG 1 MIKROKANONIKUS SOKASÁG: N részecske E összenergiával V térfogatban Mit jelent pontosan

MIKROKANONIKUS SOKASÁG 1 MIKROKANONIKUS SOKASÁG: N részecske E összenergiával V térfogatban Mit jelent pontosan a sokaság ? ? Az összes olyan elrendezést, amely N részecskét tartalmaz V térfogatban, amelyeknek az összenergiája E a sokasághoz tartozónak tekintünk. A KVANTUMMECHANIKA szerint nagy részecskeszám (N 1020) esetén a rendszernek igen sok olyan kvantumállapota lehet, amely az előírt N, E, V értékekkel összeegyeztethető – ezek alkotják a sokaságot.

Mikrokanonikus sokaság 2 Tegyük fel, hogy a rendszer éppen egy adott kvantumállapotban van. Ez

Mikrokanonikus sokaság 2 Tegyük fel, hogy a rendszer éppen egy adott kvantumállapotban van. Ez az állapot azonban nem változatlan: a kvantummechanika szerint még egyetlen részecske különböző állapotai között is van véges átmeneti valószínűség. 1020 nagyságrendű részecskeszám esetén a rendszer kvantumállapotai nagyon sűrűn helyezkednek el, és az állapotok között igen nagy valószínűséggel, gyorsan történnek átmenetek. Ezeket az átmeneteket véletlen folyamatok váltják ki, ezért várható, hogy a részecskerendszer bármely állapotból bármely másikba átmehet. Ezt fejezi ki a statisztikus fizika egyetlen alapfeltevése (axiómája): Mikrokanonikus sokaságban az állandó E, V, N feltétellel összeegyeztethető minden egyes kvantumállapot ugyanakkora valószínűséggel fordul elő, ha a sokaság a rendszer egyensúlyát reprezentálja.

Egyenletes eloszlás A mikrokanonikus sokaság állapotai tehát egyensúlyban egyenletes eloszlásúak az eseményteret alkotó kvantumállapotok

Egyenletes eloszlás A mikrokanonikus sokaság állapotai tehát egyensúlyban egyenletes eloszlásúak az eseményteret alkotó kvantumállapotok fölött. Mit jelent ez? Ha az (E, V, N ) mikrokanonikus sokaságnak megfelelő rendszert egy adott kvantumállapotból indítjuk, egyensúly akkor lehet, ha a sokaság az összes lehetséges állapotot – tehát az elérhető állapotok maximumát – azonos valószínűséggel járja be. Vizsgáljuk meg az analógiát az állapotok maximumának bejárása és az izolált rendszerek entrópiája között. Osszunk ketté gondolatban egy adott (E, V, N ) mikrokanonikus sokaságnak megfelelő termodinamikai rendszert: E 1 , V 1 , N 1 E 2 , V 2 , N 2

Mikrokanonikus entrópia 1 E 1 , V 1 , N 1 E 2 ,

Mikrokanonikus entrópia 1 E 1 , V 1 , N 1 E 2 , V 2 , N 2 (Ez megfelel egy gondolatban két alrendszerre osztott izolált rendszernek. ) Erre felírhatjuk az E = E 1 + E 2 V = V 1 + V 2 N = N 1 + N 2 összefüggéseket. Tudjuk, hogy a két alrendszer fölött az entrópia is additív: S (E, V, N ) = S 1 (E 1 , V 1 , N 1 ) + S 2 (E 2 , V 2 , N 2 )

Mikrokanonikus entrópia 2 Az állapotok számáról viszont azt tudjuk, hogy a teljes rendszerben a

Mikrokanonikus entrópia 2 Az állapotok számáról viszont azt tudjuk, hogy a teljes rendszerben a két rendszer állapotainak száma összeszorzódik: (E, V, N ) = 1 (E 1 , V 1 , N 1 ) 2 (E 2 , V 2 , N 2 ) Írjuk fel a fenti szorzatfüggvény teljes differenciálját: d = 2 d 1 + 1 d 2 , alakban, majd ezt az egyenletet osszuk el az előzővel: Az egyesített rendszerben az S = S 1 + S 2 additivitás miatt az S teljes differenciálja felírható d. S = d. S 1 + d. S 2 alakban.

Mikrokanonikus entrópia 3 és d. S = d. S 1 + d. S 2

Mikrokanonikus entrópia 3 és d. S = d. S 1 + d. S 2 egyszerre csak akkor lehet igaz, ha fennáll az alábbi összefüggés: S ( E, V, N ) = k ln ( E, V, N ) Vegyük észre, hogy ez az entrópia kifejezése a mikrokanonikus rendszer összes lehetséges állapotának számával. Ez akkor lesz összhangban a termodinamikában használt Kelvin-féle hőmérsékleti skálával, ha a k együttható éppen a Boltzmann konstans: k = 1. 3807 · 10– 23 J / K (A konstans értéke könnyebben megjegyezhető k = R / NA alakban. ) k értékének ilyen megválasztása mellett a termodinamikából jól ismert deriváltra éppen a Kelvin hőmérséklet reciprokát kapjuk.

A MIKROKANONIKUS FORMALIZMUS ÖSSZEFOGLALÁSA A sokaság i állapotainak valószínűségi sűrűségfüggvénye: N részecske E összenergiával

A MIKROKANONIKUS FORMALIZMUS ÖSSZEFOGLALÁSA A sokaság i állapotainak valószínűségi sűrűségfüggvénye: N részecske E összenergiával V térfogatban Az mennyiség neve: állapotösszeg, vagy partíciós függvény. (Zustandssumme) (partition function) A statisztikus termodinamika mikrokanonikus formalizmusa az entrópia-bázisú fundamentális egyenlethez vezet, melynek alakja: S ( E, V, N ) = k ln ( E, V, N ) A mikrokanonikus „recept” tehát nagyon egyszerű: Számoljuk össze a lehetséges állapotokat, és abból egyszerűen kiszámíthatjuk az entrópiafüggvényt. Ezzel egyúttal meghatároztuk a termodinamikai rendszer minden lehetséges makroszkopikus állapotát is.

Einstein kristálymodell Nézzünk az állapotszámlálásra egy példát: Einstein kristálymodelljét (1907) Elhanyagolta: a magenergia változását,

Einstein kristálymodell Nézzünk az állapotszámlálásra egy példát: Einstein kristálymodelljét (1907) Elhanyagolta: a magenergia változását, az elektronok mozgását, és a rezgéseken kívül minden más gerjesztést. Számításba vette: a kristály rezgéseit. Az N atom a kristály rácspontjaihoz a kitéréssel arányos (harmonikus) erővel van kötve. Így a tér mindhárom irányában 0 sajátfrekvenciával rezeghet. A szomszédos rácspontok rezgései csatoltak, így az egész kristály 3 N normál rezgési módusban rezeghet. A kialakuló csatolt rezgések alsó frekvenciahatára 0 (ez a legnagyobb hullámhosszú rezgés, ami „belefér” a kristályba), felső frekvenciahatára a rácspontok távolságának nagyságrendjébe esik. Einstein ekvivalens modellje: 3 N harmonikus oszcillátor, 0 sajátfrekvenciával.

Einstein kristálymodell 2 Az egyes oszcillátorok energiája lehet. Célszerű az energiaskálát úgy választani, hogy

Einstein kristálymodell 2 Az egyes oszcillátorok energiája lehet. Célszerű az energiaskálát úgy választani, hogy n értéke 0, 1, 2, . . . legyen. a Planck állandó Az állapotszámlálásnál felhasználjuk, hogy összesen U / rezgési kvantumot kell szétosztani 3 N rezgési mód között. 0 Itt érkeztünk el a konkrét számolás elvégzéséhez: az a szám, ahányféleképpen szétoszthatunk U / energiakvantumot 3 N oszcillátor között. 0 Ez a (kombinatorikai) feladat azonos azzal, hogyan lehet U / 0 db megkülönböztethetetlen korongot elhelyezni 3 N számozott (megkülönböztethető) dobozban.

Einstein kristálymodell 3 Használjunk dobozok helyett pálcikákat (ezekből már csak 3 N – 1

Einstein kristálymodell 3 Használjunk dobozok helyett pálcikákat (ezekből már csak 3 N – 1 kell), és tegyük le sorba egymás mellé a pálcikákat és a korongokat!. . . 1 2 3 4 3 N – 1 3 N A lehetséges elrendezések száma U / 0 + 3 N – 1 elem ismétléses permutációja, ahol egyik elemből U / 0 azonos, a másikból pedig 3 N – 1 azonos: kifejezése tehát alakú.

Einstein kristálymodell 4 Ez éppen az lehetséges állapotszám, amiből az entrópia számítható. A számításhoz

Einstein kristálymodell 4 Ez éppen az lehetséges állapotszám, amiből az entrópia számítható. A számításhoz használjuk fel a faktoriálisra vonatkozó egyik Stirling-formulát: ln (N !) N ln N – N , ha N 1 Itt a feltétel teljesül, mivel 1020 nagyságrendű szám faktoriálisát kell számítani. Az S = k ln -hoz számítsuk ki előbb ln -t. A feladat tehát az előbbiek alapján kiszámítása:

Einstein kristálymodell 5 a=3 N b = U / 0 1 mol kristályra: Tudjuk:

Einstein kristálymodell 5 a=3 N b = U / 0 1 mol kristályra: Tudjuk: 0 N = NA , S = s , U = u k NA = R , legyen: (innentől intenzív s ) u 0 = 3 NA

Einstein kristálymodell 6 Ezzel máris kész az Einstein kristálymodell fundamentális egyenlete: (mivel s intenzív,

Einstein kristálymodell 6 Ezzel máris kész az Einstein kristálymodell fundamentális egyenlete: (mivel s intenzív, ez már csak u és v függvénye) Az ebből számolt fajhő T = 0 K-nél zérus, onnan exponenciálisan emelkedik, majd telítésbe megy át, ami kvalitatíve megfelel a kísérleti hőkapacitások lefutásának. A modell primitívségéről árulkodnak annak mechanikai tulajdonságai: A egyenlet alapján a nyomás alakban számítható. A fundamentális egyenlet viszont nem mutat V-függést, így a nyomásra mindig zérust kapnánk. (A 3. főtétel viszont teljesül. )

Kétállapotú molekulák Mechanikai szempontból tehát az Einstein kristálymodell nem használható – még kvalitatíve sem.

Kétállapotú molekulák Mechanikai szempontból tehát az Einstein kristálymodell nem használható – még kvalitatíve sem. Kétállapotú molekulák Tekintsünk egy másik egyszerű példát, az ún. kétállapotú molekulákat. Ebben az esetben minden molekulának van egy „alapállapota”, amit zérus energiájúnak tekintünk, valamint egy „gerjesztett állapota”, aminek energiája legyen . (Azaz, feltesszük, hogy a többi energiaszint annyival magasabb -nál, hogy azok nincsenek betöltve soha. ) Az U, N, V mikrokanonikus sokaságban ekkor U / darab molekula van gerjesztett állapotban, N – (U / ) pedig alapállapotban. Az analóg kombinatorikai probléma U / molekula kiválasztása az összes molekula közül. Ez éppen N elem közül U / elem ismétlés nélküli kombinációinak száma: ( alakú kifejezés)

Kétállapotú molekulák 2 Az entrópia kiszámításához ismét használjuk a Stirling képletet: Visszahelyettesítve: Láthatóan ez

Kétállapotú molekulák 2 Az entrópia kiszámításához ismét használjuk a Stirling képletet: Visszahelyettesítve: Láthatóan ez az entrópia is független a térfogattól. A modell termikus tulajdonságai azonban jók: ; Mivel U < N , ezért

Kétállapotú molekulák 3 Ha a fenti állapotegyenletet U -ra megoldjuk, akkor az energiafüggvényt kapjuk.

Kétállapotú molekulák 3 Ha a fenti állapotegyenletet U -ra megoldjuk, akkor az energiafüggvényt kapjuk. Ennek alapján azaz nagyon magas hőmérsékleten a molekulák fele gerjesztett állapotban, fele alapállapotban lenne – ha ilyenkor nem töltődnének be magasabb energiájú gerjesztett állapotok is. Az energia deriválásából kapható a c. V moláris hőkapacitás:

Kétállapotú molekulák 4 Ez a függvény lefutású Valóban, a kísérleti adatok azt bizonyítják, hogy

Kétállapotú molekulák 4 Ez a függvény lefutású Valóban, a kísérleti adatok azt bizonyítják, hogy ha a c. V ilyen maximum szerint változó függvény, akkor két alacsony energiájú állapot fölött nagyon messze van a következő energiaszint. Nem véletlenül foglalkoztunk ilyen egyszerű rendszerekkel: bonyolultabb esetben a kombinatorikus analógiák már nem vezetnek kiszámítható -hoz. Van ugyan egy másik lehetőség is az összeszámolás helyett, ami sokdimenziós terekben az azonos energiákhoz tartozó felületi pontok számának integrálással történő közelítő számítása, de ez sem igazán célravezető bonyolultabb rendszerek esetén. Megoldás: az entrópia helyett más, entrópia-jellegű függvények számítása más sokaságokon.

KANONIKUS SOKASÁG 1 Példaként nézzük a kanonikus sokaságot. KANONIKUS SOKASÁG: N részecske V térfogatban

KANONIKUS SOKASÁG 1 Példaként nézzük a kanonikus sokaságot. KANONIKUS SOKASÁG: N részecske V térfogatban T hőmérsékleten A rendszer, amit ez a sokaság modellez, régi ismerősünk: N részecske V térfogatban, amely belemerül egy T hőmérsékletű hőtartályba (fala pedig hővezető) (A termodinamikából már tudjuk, hogy ez az egyensúly az F szabadenergiával (Helmholtz potenciállal) jellemezhető könnyen. ) Ebben az esetben az energia a sokaságban nincs lerögzítve, így bármely E energia elérhető a sokaság számára. A különböző állapotok azonban ebben a sokaságban nem azonos valószínűségűek.

Kanonikus sokaság 2 Az állapotok feletti eloszlás kiszámításához a rendszert és a hőtartályt egyetlen

Kanonikus sokaság 2 Az állapotok feletti eloszlás kiszámításához a rendszert és a hőtartályt egyetlen izolált összetett rendszernek tekintjük. (Így jártunk el a termodinamikában is. ) Tekintsünk egy analógiát! Legyen egy piros és dobókockánk. (A fehér kockák jelentik a tartályt. ) Nézzük először azt, mekkora a valószínűsége a három kocka egy adott „konfigurációjának”. Nyilván minden konfiguráció valószínűsége azonos. (1 / 63 = 1 / 216) Ha csak a piros kockát dobjuk egyedül, ott is azonos minden szám valószínűsége. Így az 1, 2, 3, 4, 5, 6 számok mindegyike 1/6 valószínűséggel fordul elő. Számítsuk ki, mekkora a piros kockán az 1, 2, 3, 4, 5, 6 számok előfordulásának feltételes valószínűsége, ha a három kocka dobásakor csak azokat az eseteket vesszük figyelembe, ahol a dobott számok összege 12. (Az „összenergia” megfelelője az egyesített rendszerben. )

Kanonikus sokaság 3 A piros kockán az 1, 2, 3, 4, 5, 6 számok

Kanonikus sokaság 3 A piros kockán az 1, 2, 3, 4, 5, 6 számok előfordulásának valószínűsége a három kocka dobásakor, ha az összeg 12, a következőképpen számítható: piros 1 2 3 56, 65 4 5 6 44, 26, 62, 35, 53 55, 46, 64 36, 63, 45, 54 16, 61, 25, 52, 34, 43 33, 15, 51, 24, 42 összesen lehetőségek száma valószínűség 2 3 4 2 / 25 3 / 25 4 / 25 5 6 5 25 5 / 25 6 / 25 5 / 25 1 Miközben tehát a 3 kocka összességére igaz az egyenletes eloszlás, ugyanakkor, ha lekötjük a három szám összegét, abból egy kiemelt kocka különböző állapotainak feltételes valószínűségei nem azonosak.

Kanonikus sokaság 4 Az előbbi analógia alapján levonhatjuk az alábbi következtetést: A termosztált alrendszer

Kanonikus sokaság 4 Az előbbi analógia alapján levonhatjuk az alábbi következtetést: A termosztált alrendszer adott i állapotának valószínűsége kiszámítható az alrendszer Ei energiájával összeegyeztethető tartály-állapotok számának és a teljes rendszer (tartály + alrendszer) összes állapotai számának hányadosaként: Etot : az összenergia Ei : a termosztált rendszer energiája res (Etot – Ei ) : az elérhető állapotok száma a tartályban, ha a termosztált rendszer energiája Ei tot (Etot) : az egyesített rendszer összes elérhető állapotának száma

Kanonikus sokaság 5 A kifejezésbe írjuk be az állapotok számának helyébe annak entrópiával kifejezett

Kanonikus sokaság 5 A kifejezésbe írjuk be az állapotok számának helyébe annak entrópiával kifejezett értékét a megfelelő mikrokanonikus összefüggés alapján: Ha a termosztált alrendszer energiájának egyensúlyi értéke (az energia várható értéke) U, akkor az entrópia additivitása miatt felírható: S (Etot) = S (U ) + Sres (Etot – U ) Fejtsük sorba az Sres (Etot – Ei) entrópiát az U egyensúlyi energia körül: Sres (Etot – Ei) = Sres (Etot – U + U – Ei) = Sres (Etot – U ) +

Kanonikus sokaság 6 A termodinamikából tudjuk, hogy és azt is, hogy a tartállyal cserélt

Kanonikus sokaság 6 A termodinamikából tudjuk, hogy és azt is, hogy a tartállyal cserélt energia csak , , így a sorfejtésben elegendő az első tag: Az entrópia additivitása miatt felírható: Sres és Stot fenti két kifejezését helyettesítsük be a pi kitevőjébe:

Kanonikus sokaság 7 alapján már látszik az, amit a termodinamikai előzmények ismeretében sejthetünk: az

Kanonikus sokaság 7 alapján már látszik az, amit a termodinamikai előzmények ismeretében sejthetünk: az eredményben a tartályra vagy az egyesített rendszerre jellemző mennyiségek már nem, csak a termosztált rendszer (a kanonikus sokaság) jellemzői szerepelnek. A keresett sűrűségfüggvény a fentiek alapján az alábbi alakba írható: Az egyszerűség kedvéért vezessük be a jelölést, és a már ismert függvényt. Az F függvény statisztikus kifejezését természetesen nem ismerjük. Éppen azt szeretnénk megkapni.

Kanonikus sokaság 8 F kiszámításához felhasználjuk, hogy , azaz Ebből: Ezt logaritmálva és átrendezve:

Kanonikus sokaság 8 F kiszámításához felhasználjuk, hogy , azaz Ebből: Ezt logaritmálva és átrendezve: visszahelyettesítésével: ahol Q a kanonikus partíciós függvény: Vegyük észre: itt Q a T, V, N függvénye!

Kanonikus sokaság 9 Kanonikus sokaságon tehát az szabadenergia számítható egyszerűen, ahol Q a kanonikus

Kanonikus sokaság 9 Kanonikus sokaságon tehát az szabadenergia számítható egyszerűen, ahol Q a kanonikus partíciós függvény: Az F függvény átírható a mikrokanonikus entrópiához hasonló alakba: A kanonikus partíciós függvény logaritmusának k-szorosa tehát a Massieu-függvényt adja meg, ami éppen az F-nek megfelelő entrópia-reprezentációjú potenciálfüggvény.

Kanonikus eredmények Megfogalmazhatjuk tehát a kanonikus sokaságból kiszámítható fundamentális egyenletre vonatkozó szabályt: KANONIKUS SOKASÁG:

Kanonikus eredmények Megfogalmazhatjuk tehát a kanonikus sokaságból kiszámítható fundamentális egyenletre vonatkozó szabályt: KANONIKUS SOKASÁG: N részecske V térfogatban T hőmérsékleten A termodinamikai rendszer, amit ez a sokaság modellez: N részecske V térfogatban, amely belemerül egy T hőmérsékletű hőtartályba (a rendszer fala hővezető) A kanonikus sokaság i állapotainak és a hozzájuk tartozó Ei energiáknak függvényében kiszámíthatjuk a Q kanonikus partíciós függvényt. Ebből egyszerűen számítható az F (T, V, N ) függvény – a megfelelő termodinamikai rendszer fundamentális egyenlete.

A KANONIKUS FORMALIZMUS ÖSSZEFOGLALÁSA KANONIKUS SOKASÁG: N részecske V térfogatban T hőmérsékleten a sokaság

A KANONIKUS FORMALIZMUS ÖSSZEFOGLALÁSA KANONIKUS SOKASÁG: N részecske V térfogatban T hőmérsékleten a sokaság i állapotainak valószínűségi sűrűségfüggvénye a kanonikus állapotösszeg, vagy kanonikus partíciós függvény A statisztikus termodinamika kanonikus formalizmusa a szabadenergia-bázisú fundamentális egyenlethez vezet, melynek alakja: A kanonikus „recept” tehát a következő: számítsuk ki az állapotösszeget, és abból a szabadenergia-függvényt.

Kanonikus energia A kanonikus sűrűségfüggvényből kiszámítható az energia várható értéke: Tehát a Q kanonikus

Kanonikus energia A kanonikus sűrűségfüggvényből kiszámítható az energia várható értéke: Tehát a Q kanonikus partíciós függvényből kiszámítható a belső energia is alakban.

Kanonikus partíciós függvény 1 Eddig azt reméltük, hogy majd egyszerűbb lesz a Q kiszámítása,

Kanonikus partíciós függvény 1 Eddig azt reméltük, hogy majd egyszerűbb lesz a Q kiszámítása, mint volt az kiszámítása. Ehhez az alábbi megfontolások segítenek hozzá. A kanonikus sokaság Ei energiája felírható az egymástól független (molekuláris) módusok i energiájának összegeként. Ennek feltétele, hogy a módusok között ne legyen kölcsönhatás. Ilyen esetben (a zárójelbe tett betűk különböző molekulák molekuláris módusait jelentik) A partíciós függvényre ezért felírhatjuk:

Kanonikus partíciós függvény 2 Jelöljük qj -vel a j -edik molekula molekuláris partíciós függvényét:

Kanonikus partíciós függvény 2 Jelöljük qj -vel a j -edik molekula molekuláris partíciós függvényét: Az l index a molekulák állapotát jelenti, a j pedig a molekulák sorszámát. Ezzel a jelöléssel a szorzat, és így a (teljes) partíciós függvény egyszerűbb alakban felírható: A sokaság kanonikus partíciós függvénye a független molekuláris mozgásformákhoz tartozó molekuláris partíciós függvények szorzataként állítható elő.

IDEÁLIS GÁZOK ÁLLAPOTEGYENLETE Ennyi előzmény után már nekivághatunk (ideális) gázok állapotegyenlete statisztikus kiszámításának is.

IDEÁLIS GÁZOK ÁLLAPOTEGYENLETE Ennyi előzmény után már nekivághatunk (ideális) gázok állapotegyenlete statisztikus kiszámításának is. Gázmolekulákban a következő módusok fordulhatnak elő: – három transzlációs módus – rotációs módusok (forgás) – vibrációs módusok (rezgés) – elektronállapot-módusok Ez azt jelenti, hogy a partíciós függvény faktorizálható: alakba.

Ideális gázok 2 Mindegyik tényezőre felírható a szorzatra bontás. Hasonlóképpen, mivel a szorzások sorrendje

Ideális gázok 2 Mindegyik tényezőre felírható a szorzatra bontás. Hasonlóképpen, mivel a szorzások sorrendje felcserélhető, ezt úgy is írhatjuk, hogy először a molekuláris partíciós függvényeket írjuk fel alakban, majd ezeket szorozzuk össze az N molekula partíciós függvényéhez. Kezdjük az (ideális) gáz transzlációs partíciós függvényével. Számítsuk ki ehhez először egy gázmolekula egydimenziós transzlációs mozgásának energiáit. Az egydimenziós „részecske a dobozban” probléma megoldásaként kapjuk az energiát (L hosszúságú dobozra):

Ideális gázok 3 A legalacsonyabb energiaszintet az energiaskála zéruspontjának tekintve az energiára írható: Az

Ideális gázok 3 A legalacsonyabb energiaszintet az energiaskála zéruspontjának tekintve az energiára írható: Az összes energiaszintre számított transzlációs partíciós függvény: Mivel ezek az energiaszintek nagyon közel vannak egymáshoz a makroszkopikus méretű rendszerben, ezért az összegzést helyettesíthetjük egy integrálással:

Ideális gázok 4 A integrál kiszámításához legyen A változócseréhez használjuk a következő összefüggéseket: így

Ideális gázok 4 A integrál kiszámításához legyen A változócseréhez használjuk a következő összefüggéseket: így Az egy dimenzióban: energia beírása után az eredmény:

Ideális gázok 5 Ha a mozgás háromdimenziós, X , Y és Z dobozmérettel: mivel

Ideális gázok 5 Ha a mozgás háromdimenziós, X , Y és Z dobozmérettel: mivel X Y Z = V, a doboz térfogata. Vezessük be a termikus hullámhosszat. Ennek segítségével a háromdimenziós transzlációs partíciós függvény alakba írható.

Ideális gázok 6 Nézzük meg, mekkora a termikus hullámhossz és a molekuláris transzlációs partíciós

Ideális gázok 6 Nézzük meg, mekkora a termikus hullámhossz és a molekuláris transzlációs partíciós függvény 25 ºC-on , 100 cm 3 O 2 molekulát tartalmazó gázban. Szobahőmérsékletű gázokban ez a tipikus nagyságrendje a termikus hullámhossznak és a molekuláris transzlációs partíciós függvénynek.

Ideális trans gázok 7 q -ból már kiszámíthatjuk egyatomos gázok állapotegyenletét is. (Ezeknek nincs

Ideális trans gázok 7 q -ból már kiszámíthatjuk egyatomos gázok állapotegyenletét is. (Ezeknek nincs forgási és rezgési állapotuk, gerjesztett elektronállapotukat pedig tekintsük sokkal nagyobb energiájúnak az alapállapotnál. ) Ekkor N atomból álló gázra: Ez a formula azonban nem jó így, mert a valódi állapotok száma az N elemű rendszerben kisebb, mint ami a fenti képlettel számítható. A gázmolekulák ugyanis nem megkülönböztethetők. Ha tehát a „megszámozott” molekulákat összekeverjük, nem jutunk más konfigurációhoz. A fenti számot ezért el kell osztani N!-sal, mert ennyiféleképpen „átszámozhatók” a molekulák. (N elem ismétlés nélküli permutációinak száma. ) A „jó” partíciós függvény tehát:

Ideális gázok 8 Ebből most már ki tudjuk számítani az ideális gáz fundamentális egyenletét.

Ideális gázok 8 Ebből most már ki tudjuk számítani az ideális gáz fundamentális egyenletét. Tudjuk: Ebből

Ideális gázok 9 Az U = F + TS alapján kiszámítható a belső energia

Ideális gázok 9 Az U = F + TS alapján kiszámítható a belső energia is: Ebből kifejezhető T az U függvényében: Ezt beírva a hőmérséklet helyére, az S(U, V, N ) számítható: Sackur-Tetrode egyenlet (fundamentális)

Rotációs partíciós függvény Az U – ból könnyen kapható CV is: Itt is látszik

Rotációs partíciós függvény Az U – ból könnyen kapható CV is: Itt is látszik pl. az, hogy a statisztikus állapotegyenletek konstansának k = R / NA választása vezet a termodinamikából ismert képletekhez. Láttuk tehát, milyen hatékony eszköz a kanonikus sokaság. Nézzük most a partíciós függvény többi tényezőjét! A forgási partíciós függvényhez a forgási energia kifejezése kell. Vizsgáljunk meg egy heteronukleáris lineáris rotort (mint pl. a HCl). A lineáris rotor rotációs állandója: ( I: tehetetlenségi nyomaték) A lineáris rotor energiája: (J: forgási kvantumszám) a rotációs partíciós függvény:

Rotációs partíciós függvény 2 ÚJ ELEM: A (2 J + 1)-es faktor. (Eddig ilyen

Rotációs partíciós függvény 2 ÚJ ELEM: A (2 J + 1)-es faktor. (Eddig ilyen még nem volt. ) Tudjuk: a J forgási kvantumszámhoz tartozó energiaszint (2 J + 1)-szeresen degenerált. Ezért kell (2 J + 1)-el szorozni J szerinti (állapot szerinti, és nem energia szerinti) összegzés esetén. Az összeg kiszámítása ismét integrálásra vezethető vissza, ha h c B << k T (h c B a J = 1 -hez tartozó energia) Helyettesítés:

Rotációs partíciós függvény 3 A heteronukleáris lineáris rotornak (mint pl. a HCl) 2 szabadsági

Rotációs partíciós függvény 3 A heteronukleáris lineáris rotornak (mint pl. a HCl) 2 szabadsági foka van. Ennek megfelelően két független forgástengelye van, történetesen mindegyikhez ugyanaz a rotációs állandó, így ugyanazok a forgási energiák tartozhatnak. Az általános rotornak 3 szabadsági foka, ennek megfelelően 3 forgástengelye és 3 rotációs állandója van. Jelöljük ezeket A, B és C betűkkel. Egy ilyen általános (nemlineáris) molekula rotációs partíciós függvénye: ( : a forgási szimmetriacsoport degenerációfoka) Szobahőmérsékleten sok forgási állapot energiája elérhető. Ebben a hőmérséklettartományban a tipikus érték q R 1000.

Rezgési partíciós függvény A rezgési partíciós függvényhez írjuk fel a rezgési energiát. Ez harmonikus

Rezgési partíciós függvény A rezgési partíciós függvényhez írjuk fel a rezgési energiát. Ez harmonikus oszcillátorra: v : latin „vé” : görög „nű” Legyen az energiaskála zéruspontja Ekkor: Ez éppen az quotiensű mértani sor összege:

Rezgési / elektronikus partíciós függvény A molekula minden (harmonikus) normálrezgésére ilyen alakú a partíciós

Rezgési / elektronikus partíciós függvény A molekula minden (harmonikus) normálrezgésére ilyen alakú a partíciós függvény. Így a teljes rezgési partíciós függvény: Az egyes normálrezgésekhez tartozó q. V(i ) -kről lehet tudni, hogy nem túl magas hőmérsékleten (szobahőmérséklet körül) értékük 1 -3 körüli. Hátra van még az elektronikus partíciós függvény. Gyakori eset, hogy szobahőmérséklet környékén csak az alapállapot elérhető: Kivétel: degenerált alapállapotok (pl. az alkáli atomok spin-degenerációja) Ekkor: g. E a degenerációfok

Ideális gázok 10 Pl. az NO molekula esetében: 2 alapállapot + 2 „kis energiájú”

Ideális gázok 10 Pl. az NO molekula esetében: 2 alapállapot + 2 „kis energiájú” gerjesztett állapot: Ez a „kétszeresen degenerált” kétállapotú rendszer. Ezzel készen van az (ideális) gáz teljes partíciós függvényéhez szükséges minden tényező: Ennek segítségével már nemcsak egyatomos gázok állapotegyenleteit írhatjuk fel, hanem bonyolultabb molekulákból álló gázokét is.

A kanonikus sűrűségfüggvény tulajdonságai Mielőtt áttérünk a kanonikus formalizmus szerinti fundamentális egyenletek termodinamikai alkalmazására,

A kanonikus sűrűségfüggvény tulajdonságai Mielőtt áttérünk a kanonikus formalizmus szerinti fundamentális egyenletek termodinamikai alkalmazására, vizsgáljuk meg, milyen a kanonikus sűrűségfüggvény: Korábban láttuk, hogy segítségével kiszámítható az energia várható értéke:

Kanonikus energia szórásnégyzete 1 Vizsgáljuk meg most az energia szórásnégyzetét. Tudjuk: Ennek alapján írhatjuk:

Kanonikus energia szórásnégyzete 1 Vizsgáljuk meg most az energia szórásnégyzetét. Tudjuk: Ennek alapján írhatjuk: A második tag a már ismert várható érték négyzete:

Kanonikus energia szórásnégyzete 2 Mivel alapján , ebből Belátható, hogy ez éppen M(E) =

Kanonikus energia szórásnégyzete 2 Mivel alapján , ebből Belátható, hogy ez éppen M(E) = U deriváltja: Az energia szórásnégyzetére tehát szerencsésen írhatjuk:

Kanonikus energiafluktuáció 1 Gyakran hasznos a szerinti deriválásáról T szerintire áttérni. Ezt a láncszabály

Kanonikus energiafluktuáció 1 Gyakran hasznos a szerinti deriválásáról T szerintire áttérni. Ezt a láncszabály alapján a következőképpen írhatjuk: Ebből: 2 (E) = Az energia „átlagos fluktuációja” éppen a szórásnégyzetgyöke. Az energia várható értékére vonatkoztatott relatív fluktuáció: Mivel U N ( U extenzív), ezért a relatív fluktuáció:

KANONIKUS ENERGIAFLUKTUÁCIÓ 2 Az N részecskeszámmal az U belső energia nő, miközben relatív szórása

KANONIKUS ENERGIAFLUKTUÁCIÓ 2 Az N részecskeszámmal az U belső energia nő, miközben relatív szórása szerint csökken. Ahogy tehát N , Emlékezzünk arra, hogy mikrokanonikus sokaságban az energia állandó. Tehát ha N , akkor a mikrokanonikus formalizmus szerinti termodinamikai mennyiségeket (pl. energiát) kell kapnunk kanonikus formalizmus szerint is. Nézzük, mekkora az energia relatív szórása makroszkópikus mennyiségű anyagra? Ekkor N ~ 1023, és Makroszkopikus méretű testekben tehát az energia szórása abszolút értékének 10– 11 -szerese, azaz mérhetetlenül kicsi.

AZ ENERGIA KANONIKUS SŰRŰSÉGFÜGGVÉNYE Vizsgáljuk meg most az energiaeloszlást kanonikus sokaságon. Tudjuk: E a

AZ ENERGIA KANONIKUS SŰRŰSÉGFÜGGVÉNYE Vizsgáljuk meg most az energiaeloszlást kanonikus sokaságon. Tudjuk: E a rendszer egy makroszkopikus energiája (ez az, ami mérhető). Az E makroszkópikus energiához tartozó mikroállapotok száma (E ). Ez azt jelenti, hogy az E energia degenerációfoka (E ). Ezért Termodinamikai rendszerekre (E ) EN (az S = k ln (E) alapján). Ha tehát a makroállapot E energiájának valószínűségét tekintjük, abban (az 1/ Q konstans szorzó mellett) szerepel még: (E ) állapotsűrűség (az energia nagyon gyorsan növekvő függvénye) e – E Boltzmann-faktor (az energia nagyon gyorsan csökkenő függvénye)

A makroszkopikus energia kanonikus sűrűségfüggvénye 2 P (E ) E N e – E

A makroszkopikus energia kanonikus sűrűségfüggvénye 2 P (E ) E N e – E (E ) M (E ) E Az E sűrűségfüggvénye az e – E és az E N szorzata, ami a gyorsan lecsengő/felfutó függvények miatt egy igen éles csúcs, amely a központi határeloszlás tétel következményeként egy Gauss-függvény. E függvény M (E ) -hez viszonyított relatív szélessége -nel arányos.

Az egyrészecske-energia kanonikus sűrűségfüggvénye 1 Ugyanakkor a (nem degenerált) molekuláris energiaállapotok betöltöttségére érvényes az

Az egyrészecske-energia kanonikus sűrűségfüggvénye 1 Ugyanakkor a (nem degenerált) molekuláris energiaállapotok betöltöttségére érvényes az exponenciálisan lecsengő sűrűségfüggvény: p ( i ) („egyrészecske-eloszlás”) T 1 < T 2 < T 3 T 2 T 3 i Mindig teljesülnie kell a feltételnek.

EKVIPARTÍCIÓ 1 Van még egy fontos tétel, az ún. „ekvipartíció” tétele. Vizsgáljuk meg ehhez

EKVIPARTÍCIÓ 1 Van még egy fontos tétel, az ún. „ekvipartíció” tétele. Vizsgáljuk meg ehhez az egyrészecske-eloszlás alábbi sűrűségfüggvényének felhasználásával a különböző molekuláris módusokhoz tartozó energia várható értékét, amit az alábbiak szerint számíthatunk:

Ekvipartíció 2 Alkalmazzuk először az várhatóérték-számítást a 3 szabadsági fokú transzlációra: A 2 szabadsági

Ekvipartíció 2 Alkalmazzuk először az várhatóérték-számítást a 3 szabadsági fokú transzlációra: A 2 szabadsági fokú rotációra:

Ekvipartíció 3 Az 1 szabadsági fokú rezgésre: Ez csak akkor egyszerűsödik, ha h <<

Ekvipartíció 3 Az 1 szabadsági fokú rezgésre: Ez csak akkor egyszerűsödik, ha h << 1 határesetben sorfejtjük alakra. Ekkor k T >> h , és

Ekvipartíció 4 Hasonlítsuk most össze az energia klasszikus kifejezését a részecskeenergia várható értékével: módus

Ekvipartíció 4 Hasonlítsuk most össze az energia klasszikus kifejezését a részecskeenergia várható értékével: módus klasszikus energia átlagos energia szabadsági fokokkal transzláció rotáció vibráció Általánosítás: minden négyzetes energiatagra (kellően magas hőmérsékleten) 1/2 k T energia jut.

EKVIPARTÍCIÓ TÉTEL Minden négyzetes energiatagra (kellően magas hőmérsékleten) 1/ 2 k T energia jut.

EKVIPARTÍCIÓ TÉTEL Minden négyzetes energiatagra (kellően magas hőmérsékleten) 1/ 2 k T energia jut. A hőkapacitás az energia (várható értékének) hőmérséklet szerinti deriváltja: A táblázatban összefoglaltak alapján a moláris hőkapacitásra írhatjuk: Ennek alapján a moláris hőkapacitásból a molekula szerkezetére lehet következtetni.

Az entrópia mint rendezetlenség 1 Az entrópia általános értelmezéséről Könnyen belátható, hogy az entrópia

Az entrópia mint rendezetlenség 1 Az entrópia általános értelmezéséről Könnyen belátható, hogy az entrópia éppen a sűrűségfüggvény logaritmusának az alakban felírható várható értéke. Mikrokanonikus esetben ugyanis , és ekkor

Az entrópia mint rendezetlenség 2 Kanonikus esetben = , és ekkor =S

Az entrópia mint rendezetlenség 2 Kanonikus esetben = , és ekkor =S

Az entrópia mint rendezetlenség 3 Milyen értelemben nevezhető „rendezetlenségnek” az 1) Ha a pi

Az entrópia mint rendezetlenség 3 Milyen értelemben nevezhető „rendezetlenségnek” az 1) Ha a pi sűrűségfüggvényben valamelyik i-hez 1 tartozik (az összes többi zérus), akkor A rendezetlenség S = – k · 1 · ln 1 = 0 mértéke zérus 2) A rendezetlenség maximuma az egyenletes eloszláshoz tartozik – ha a lehetőségek közül egy sincs kitüntetve. (véletlenszerű bejárás). 3) A rendezetlenség maximuma monoton nő, ha az egyenletes eloszlásban a lehetőségek száma ( ) nő. 4) A rendezetlenség additív a részrendszerek felett: Az entrópiára a „rendezetlenség” kifejezést tehát (csak) a fenti értelemben – azaz körültekintően – lehet alkalmazni.

Az egyensúlyi állandó kanonikus kifejezése 1 Használjuk ki eddig tanult termodinamikai apparátusunkat kémiai reakciók

Az egyensúlyi állandó kanonikus kifejezése 1 Használjuk ki eddig tanult termodinamikai apparátusunkat kémiai reakciók egyensúlyi állandóinak számításához. Fejezzük ki először F-et molekuláris állapotösszegből gázokra: F(0) rögzíti az energiaskálát alkalmazzuk a Stirling formulát Ebből a moláris F (k NA = R felhasználásával):

Az egyensúlyi állandó kanonikus kifejezése 2 T, P = állandó feltételek mellett a kémiai

Az egyensúlyi állandó kanonikus kifejezése 2 T, P = állandó feltételek mellett a kémiai reakció egyensúlyi állandója G -vel fejezhető ki: G = F + PV Csaljunk most egy kicsit, és legyen P V = n R T, azaz korlátozódjunk ideális gázokra: Jelöljük -lal azt a qm-et, amit mellett számíthatunk. Ezzel a jelöléssel:

Az egyensúlyi állandó kanonikus kifejezése 3 Fm(0) -t visszavezethetjük az F = U –TS

Az egyensúlyi állandó kanonikus kifejezése 3 Fm(0) -t visszavezethetjük az F = U –TS segítségével az U(0) -ra: Fm(0) = Um(0) – T(0) S(0) Ha T(0) = 0 K Tudjuk: Fm(0) = Um(0) ; Hasonlítsuk össze ezt a fenti -lal: rövidítés: r Um(0) = r. U 0

Az egyensúlyi állandó kanonikus kifejezése 4 Ezzel (kis csalással) megvan az egyensúlyi állandó statisztikus

Az egyensúlyi állandó kanonikus kifejezése 4 Ezzel (kis csalással) megvan az egyensúlyi állandó statisztikus kifejezése – ami leginkább csak ideális gázelegyekben használható igazán. A kifejezés általában is használható, ha a helyébe az „igazi” Qm, i kanonikus partíciós függvényt helyettesítjük. A kitevőben ilyenkor is a T = 0 hőmérsékleten számított reakcióenergia szerepel. Összefoglalva: A gázreakció egyensúlyi állandója a kanonikus formalizmus alapján: a molekuláris partíciós függvény, r U 0 a reakcióenergia T = 0 K hőmérsékleten

Az egyensúlyi állandó kanonikus kifejezése 5 Írjuk fel ezt két konkrét gázreakcióra molekuláris partíciós

Az egyensúlyi állandó kanonikus kifejezése 5 Írjuk fel ezt két konkrét gázreakcióra molekuláris partíciós függvényekkel: Ha A+B C „bimolekulás reakció” akkor Ha akkor A B „unimolekulás reakció”

Az egyensúlyi állandó kanonikus kifejezése 6 Ugyanez moláris partíciós függvényekkel: Ha A+B C „bimolekulás

Az egyensúlyi állandó kanonikus kifejezése 6 Ugyanez moláris partíciós függvényekkel: Ha A+B C „bimolekulás reakció” akkor Ha akkor A B „unimolekulás reakció”

VÉGE Vég e

VÉGE Vég e