3.1. Fluktuatsioonid Gaussi lähenduses

Kui vaadeldav süsteem on termodünaamilises tasakaalus, siis selle termodünaamilised karakteristikud omandavad väärtusi, mis on harilikult väga lähedased vastavatele keskväärtustele. Siiski on termodünaamilises tasakaalus olemas nullist erinev tõenäosus selleks, et mingi süsteemi iseloomustava suuruse väärtus antud ajahetkel ei lange kokku tasakaalulise väärtusega. Realselt ostsillerub see suurus ajas juhuslikult väikese amplituudiga oma tasakaalulise väärtuse ümber ning selle juhuslikku kõrvalekaldumist tasakaalulisest väärtusest nimetakse fluktuatsiooniks.

Olgu meil mingi suurus \(a\), mis iseloomustab kas süsteemi tervikuna või siis süsteemi osa. Selle suuruse tasakaalulist väärtust tähistame \(\bar a\). Suuruse \( a\) fluktuatsiooniks on vahe \(\Delta a= a-\bar a\).

3.1.1. Boltzmanni printsiip

Fluktuatsioonide poolfenomenoloogiline kirjeldamine lähtub A. Einsteini valemist fluktuatsiooni tõenäosustiheduse jaoks

(3.1)\[w( a)=Ae^\frac{S( a)}{k_\mathrm{B}},\qquad A^{-1}=\int e^\frac{S( a)}{k_\mathrm{B}}d a.\]

Siin \(S( a)\) on entroopia mittetasakaalulise seisundi jaoks, mis on tekkinud tasakaalulise oleku fluktuatsioonina. Kordaja \(A\) tagab tõenäosustiheduse normeeringu. Valemi kohaselt on \(w( a)d a\) tõenäosus selleks, et parameetri väärtused on vahemikus \([ a, a+d a]\).

Valemi tuletamisel on kasutatud Boltzmanni poolt leitud seos entroopia ja mikrokanoonilise jaotuse vahel, mis oli üldistatud fluktuatsioonidele, \(S( a_2)-S( a_1)=k_\mathrm{B}\ln\frac{w( a_2)}{w( a_1)}\). Üldistuse kohaselt ilmub entroopia ja tõenäosuse vahel põhimõtteline seos. Kui kinnine süsteem asub mittetasakaalulises olekus, siis paratamatult ja iseenesest läheb see relaksatsiooni aja jooksul tasakaaluolekusse, millega kaasneb süsteemi entroopia kasvamine kuni maksimaalse väärtuseni. Samal ajal näeme, et süsteem püüab selles protsessis üle minna vähem tõenäolisest olekust tõenäolisemasse olekusse.

Entroopia seos tõenäosusega muudab ka kategoorilise entroopia kasvu printsiibi tõenäosusliku tähendusega printsiibiks: isoleeritud süsteemis on kõige tõenäolisemad sellised protsessid, mille käigus süsteemi entroopia kasvab või jääb muutumatuks (süsteemi oleku tõenäosus suureneb, või jääb muutumatuks).

3.1.2. Üks fluktueeruv suurus

Suuruse \( a\) tasakaalulise väärtuse korral on süsteemi entroopia maksimaalne \(S( a)|_{ a=\bar a}=\textnormal{max}\). See tähendab, et kehtib funktsiooni \(S\) maksimumi tingimus

\[\frac{\partial S}{\partial a}\bigg|_{ a=\bar a}=0,\qquad \beta\equiv-\frac{\partial^2 S}{\partial a^2}\bigg|_{ a=\bar a}>0.\]

Oletame, et fluktuatsioonid on väikesed ja seega kehtib järgmine \(S( a)\) reaksarendus \(\Delta a\) astmete järgi

\[S( a)=S(\bar a)-\frac{\beta}{2}(\Delta a)^2.\]

Järelikult on fluktuatsiooni tõenäosustihedus (3.1) kujul

\[w( a)=\frac{e^{-\frac{\beta}{2k_\mathrm{B}}( a-\bar a)^2}}{\int e^{-\frac{\beta}{2k_\mathrm{B}}( a-\bar a)^2}d a}.\]

Paneme tähele, et kehtib \(w\sim e^{\Delta S/k_\mathrm{B}}\), kus \(\Delta S=S( a)-S(\bar a)\) on entroopia kõrvalekalle tasakaaluväärtusest seoses \( a\) fluktuatsiooniga. Kuigi me eeldasime, et fluktuatsioonid \(\Delta a\) on väikesed võime me siin integreerida vahemikus \(-\infty\ldots+\infty\), kuna integraalialune funktsioon kahaneb väga kiiresti \(|\Delta a|\) kasvades. Selline integreerimine annab

\[\int e^{-\frac{\beta}{2k_\mathrm{B}}( a-\bar a)^2}d a\approx \int\limits_{-\infty}^\infty e^{-\frac{\beta}{2k_\mathrm{B}}(\Delta a)^2}d\Delta a=\sqrt{\frac{2\pi k_\mathrm{B}}{\beta}},\]

kus viimases võrduses on kasutatud Gaussi integraali valem. Saadud jaotust

\[w( a)=\sqrt{\frac{\beta}{2\pi k_\mathrm{B}}}e^{-\frac{\beta}{2k_\mathrm{B}}( a-\bar a)^2}\]

nimetatakse Gaussi jaotuseks, vt pilti.

../_images/im31_1.svg

Fluktuatsiooni keskväärtuste jaoks saame

(3.2)\[\langle a\rangle=\int_{-\infty}^{\infty} a w( a)d a=\bar a,\qquad \langle(\Delta a)^2\rangle=\int_{-\infty}^{\infty}(\Delta a)^2 w( a)d a=\frac{k_\mathrm{B}}{\beta}.\]

Näeme, et Gaussi jaotus on täielikult määratud, kui on teada kaks keskväärtust \(\langle a\rangle\) ja \(\left\langle(\Delta a)^2\right\rangle\).

3.1.3. Mitu fluktueeruvat suurust

Võtame nüüd vaatluse alla mitme suuruse \( a_1, a_2,\ldots a_n\equiv\{ a\}\) samaaegsed fluktuatsioonid \(\Delta a_1,\Delta a_2,\ldots\Delta a_n\equiv\{\Delta a\}\). Eeldame, et fluktuatsioonid on väikesed. Pidades silmas entroopia maksimaalsust termodünaamilise tasakaalu seisundis

\[\frac{\partial S\big(\{ a\}\big)}{\partial a_i}\bigg|_{\{ a\}=\{\bar a\}}=0,\qquad i=1,2,\ldots,n,\]

esitame fluktuatsioonidest põhjustatud süsteemi entroopia muutuse reaksarendusena fluktuatsioonide astmete järgi

\[S\big(\{ a\}\big)=S\big(\{\bar a\}\big)-\frac{1}{2}\sum_{i,j}\beta_{ij} \Delta a_i\Delta a_j,\]

kus

\[\beta_{ij}=-\frac{\partial^2S\big(\{ a\}\big)}{\partial a_i\partial a_j}\bigg|_{\{ a\}=\{\bar a\}},\qquad \beta_{ij}=\beta_{ji}.\]

Boltzmanni printsiibi (3.1) üldistades meie juhule saame fluktuatsioonide tõenäosustihedust kujul

\[w\big(\{ a\}\big)=\frac{e^{-\frac{1}{2k_\mathrm{B}}\sum_{i,j}\beta_{ij}\Delta a_i\Delta a_j}}{\int_{-\infty}^\infty e^{-\frac{1}{2k_\mathrm{B}}\sum_{i,j}\beta_{ij}\Delta a_i\Delta a_j}d a_1d a_2\ldots d a_n}.\]

Oluline erinevus ühe fluktueeruva suuruse juhuga on see, et eksponent sisaldab segaliikmeid, mis vastavad summa \(i\neq j\) elemendile.

Arvuvame kõigepealt normeeritava integraali. Selleks kasutame lineaarteisendust

\[\Delta a_i=\sum_ju_{ij}\Delta\tilde a_j\]

nii, et see viiks ruutvormi \(\sum_{i,j}\beta_{ij}\Delta a_i\Delta a_j\) diagonaalsele kujule

\[\sum_{i,j}\beta_{ij}\Delta a_i\Delta a_j=\sum_{i,j}\beta_{ij}\sum_{k,l}u_{ik}u_{jl}\Delta\tilde a_k\Delta\tilde a_l\equiv\sum_{k,l}\delta_{kl}\Delta\tilde a_k\Delta\tilde a_l=\sum_k(\Delta\tilde a_k)^2.\]

Selleks peavad teisenduse koefitsiendid rahuldama tingimust

(3.3)\[\sum_{i,j}\beta_{ij}u_{ik}u_{jl}=\delta_{kl}.\]

Seoses \(\Delta a\to\Delta\tilde a\) teisendusega muutub ka diferentsiaalide korrutist

\[d a_1d a_2\ldots d a_n=d\Delta a_1d\Delta a_2\ldots d\Delta a_n= \frac{\partial (\Delta a_i)}{\partial(\Delta\tilde a_j)}d\Delta\tilde a_1d\Delta\tilde a_2\ldots d\Delta\tilde a_n=\textnormal{det}[u_{ij}]d\Delta\tilde a_1d\Delta\tilde a_2\ldots d\Delta\tilde a_n.\]

Jakobiaani \(\textnormal{det}[u_{ij}]\) leidmiseks paneme tähele, et kehtib tingimus (3.3), mille kohaselt

\[\textnormal{det}\Big[\sum_{i,j}\beta_{ij}u_{ik}u_{jl}\Big]=\textnormal{det}[\beta_{ij}]\textnormal{det}[u_{ij}]\textnormal{det}[u_{ij}]=\textnormal{det}[\delta_{ij}]=1,\]

seega \(\textnormal{det}[u_{ij}]=1/\sqrt{\textnormal{det}[\beta_{ij}]}\). Siinjuures tingimus \(\textnormal{det}[\beta_{ij}]>0\) tagab entroopia maksimumi piisava tingimuse (ühe fluktueeruva suuruse korral on vastavaks tingimuseks \(\beta>0\)). Järelikult saame normeeritava integraali jaoks

\[\int_{-\infty}^\infty e^{-\frac{1}{2k_\mathrm{B}}\sum_{i,j}\beta_{ij}\Delta a_i\Delta a_j}d a_1d a_2\ldots d a_n=\frac{1}{\sqrt{\textnormal{det}[\beta_{ij}]}}\prod_{i=1}^n\int_{-\infty}^\infty e^{-\frac{(\Delta\tilde a_i)^2}{2k_\mathrm{B}}}d\Delta\tilde a_i=\sqrt{\frac{(2\pi k_\mathrm{B})^n}{\textnormal{det}[\beta_{ij}]}}\]

Viimase võrduse korral oli jälle kasutatud Gaussi integraali valem.

Tõenäosustihedusel on meil nüüd lõplik kuju

\[w\big(\{ a\}\big)=\sqrt{\frac{\textnormal{det}[\beta_{ij}]}{(2\pi k_\mathrm{B})^n}}e^{-\frac{1}{2k_\mathrm{B}}\sum_{i,j}\beta_{ij}( a_i-\bar a_i)( a_j-\bar a_j)}.\]

Saadud tulemus on Gaussi jaotus mitme fluktueeruva suuruse jaoks, vt. pilti.

../_images/im31_2.svg

Selle jaotuse alusel võib leida, et

(3.4)\[\langle a_i\rangle=\int\limits_{-\infty}^\infty a_iw(\{ a\})d a_1d a_2\ldots d a_n=\bar a_i,\]
(3.5)\[\langle\Delta a_i\Delta a_j\rangle=\int\limits_{-\infty}^\infty \Delta a_i\Delta a_jw(\{ a\})d a_1d a_2\ldots d a_n=k_\mathrm{B}\beta^{-1}_{ij},\]

kus \(\beta^{-1}_{ij}\) on koefitsientidest \(\beta_{ij}\) moodustatud maatriksi pöördmaatriksi elemendid. Kui \(\beta^{-1}_{ij}=0\), siis \(\langle\Delta a_i\Delta a_j\rangle=0\), ehk fluktuatsioonid \(\Delta a_i\) ja \(\Delta a_j\) on statistiliselt sõltumatud.

Toome sisse veel termodünaamilise jõu definitsiooni \(X_i=-\frac{\partial S}{\partial a_i}\). See jõud tekkib süsteemis seoses selle kõrvalekallega tasakaaluasendist. Kasutades Gaussi lähendust saame \(X_i=\sum_j\beta_{ij}\Delta a_j\) ning

\[\langle X_i\rangle=0,\qquad \langle X_i\Delta a_k\rangle=k_\mathrm{B}\delta_{ik},\qquad \langle X_kX_i\rangle= k_\mathrm{B}\beta_{ik}.\]