Fizyka statystyczna/Fluktuacje i ruchy Browna

Z Wikibooks, biblioteki wolnych podręczników.
Fizyka statystyczna
Fizyka statystyczna
Fluktuacje i ruchy Browna

Licencja
Autor: Mirosław Makowiecki
Absolwent UMCS Fizyki Komputerowej Uniwersytetu Marii Curie-Skłodowskiej w Lublinie
Email: miroslaw(kropka)makowiecki(małpa)gmail(kropka)pl
Dotyczy: książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami.
Użytkownika książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami nie zwalnia z odpowiedzialności prawnoautorskiej nieprzeczytanie warunków licencjonowania.
Umowa prawna: Creative Commons: uznanie autorstwa, na tych samych warunkach, z możliwością obowiązywania dodatkowych ograniczeń.
Autor tej książki dołożył wszelką staranność, aby informacje zawarte w książce były poprawne i najwyższej jakości, jednakże nie udzielana jest żadna gwarancja, czy też rękojma. Autor nie jest odpowiedzialny za wykorzystanie informacji zawarte w książce, nawet jeśli wywołaby jakąś szkodę, straty w zyskach, zastoju w prowadzeniu firmy, przedsiębiorstwa lub spółki bądź utraty informacji, niezależnie czy autor (a nawet Wikibooks) został powiadomiony o możliwości wystąpienie szkód. Informacje zawarte w książce mogą być wykorzystane tylko na własną odpowiedzialność.


Poprzedni rozdział: Idealne gazy kwantowe.

Podręcznik: Fizyka statystyczna.

Będziemy się zajmować fluktuacjami energii (zespół kanoniczny i wielki zespół kanoniczny) i fluktuacją ilości cząstek w układzie (wielki zespół kanoniczny).

Definicja fluktuacji bezwzględnej[edytuj]

Fluktuacją nazywamy wartość średnią kwadratową odchylenia wielkości x od jej wartości średniej i będziemy je oznaczać: , zatem z definicji średniego odchylenia możemy napisać:

(26.1)

Zatem ze wzoru (26.1) dostajemy równoważne do poprzedniego wyrażenie:

(26.2)

Na podstawie (26.2) widzimy, że fluktuacja zależy od średniej kwadratu wielkości x i od średniej tegoż wyrażenia, ten wzór (26.2) jest alternatywnym zgodnym wzorem z początkową definicją fluktuacji średniej dla badanego układu.

Definicja fluktuacji względnej[edytuj]

Fluktuacją względną nazywamy wielkość zddefiniowanej za pomocą ilorazu fluktuacji bezwzględnej zdefiniowanej wedle wzoru (26.2) przez wartość statystyczną średnią danego układu, i oznaczmy ją literką "k", omawiana wielkość jest równa:

(26.3)

Fluktuacje w zespole kanonicznym (T,V,N,t)[edytuj]

Średnia energia układu według zespołu kanonicznego, którego wyprowadzenie napisaliśmy w punkcie (17.4) jest ona równoważna do wspomnianego wzoru na tą średnią wzorowi wyrażone poprzez sumę statystyczną "Z", ale nie używają logarytmu naturalnego:

(26.4)

A średnia kwadratu energii układu wyrazimy przez sumę statystyczną naszego zespołu kanonicznego, i którą policzymy z definicji wartości średniej kwadratu energii jakie układ statystyczny może przyjmować:

(26.5)

Średni kwadrat fluktuacji bezwzględnej energii cząstek znajdującej się w układzie zamkniętym, która jak wyprowadzimy jest zależna liniowo od kwadratu temperatury bezwzględnej układu oraz liniowo od ciepła właściwego pod stałą objętością:

(26.6)

W rezultacie na podstawie wzoru (26.6) wzór na kwadrat odchylenia kwadratowego piszemy wedle:

(26.7)

Rząd fluktuacji energii[edytuj]

Jeśli energia wewnętrzna na podstawie gazu doskonałego klasycznego jest wyrażona wzorem (17.11) w układzie zamkniętym, którą dla porządku dziennego powtórzymy jego zapis poniżej:

(26.8)

Ciepło właściwe gazu doskonałego o energii wewnętrznej U, zdefiniowanej przy pomocy wzoru (26.8), dla nieoddziaływających klasycznych cząstek jest wielkością niezależna od temperatury, jest wyrażona:

(26.9)

Względna fluktuacja energii wewnętrznej układu zamkniętego (26.3) przy definicji ciepła właściwego dla tego gazu w tym układzie (26.9) przy wykorzystaniu wzoru na bezwzględną fluktuację energii (26.7), wyraża się:

(26.10)

Oczywiste jest, że według (26.10), mamy końcowy wzór na względna fluktuację liczby cząstek znajdującej się w układzie:

(26.11)

Na podstawie wzoru na względną fluktuację względną (26.11) otrzymamy, że czym większa jest liczba cząstek, to względna fluktuacja energii w danym układzie jest czym mniejsza. Dla liczby cząstek N dążących do nieskończoności względna fluktuacja wynosi zero, czyli nie ma względnych fluktuacji energii.

Fluktuacje w wielkim zespole kanonicznym (T,V,μ,t)[edytuj]

Fluktuacja energii wewnętrznej[edytuj]

Średnia energia układu w wielkim zespołu kanonicznym udowodnionych w punkcie (18.4) jest równa wzorowi wyrażone poprzez wielką sumę statystyczną "Q":

(26.12)

A średnia kwadratu energii układu wyrazimy także przez sumę statystyczną jak w przypadku (26.12) dla naszego wielkiego zespołu kanonicznego jest równa:

(26.13)

Średnia kwadrat fluktuacja energii cząstek znajdującej się w układzie zamkniętym jak się przekonamy, jest ona proporcjonalna liniowo do kwadratu temperatury bezwzględnej i liniowo do ciepła właściwego pod stałym parametrem równym objętości układu:

(26.14)

W rezultacie mamy na podstawie (26.14) równość, którą przepiszemy dla przejrzystości wykładu:

(26.15)

Rząd fluktuacji energii[edytuj]

Jeśli energia wewnętrzna dla wielkiego zespołu kanonicznego dla gazu doskonałego klasycznego jest napisana wedle wzoru (18.43) w układzie zamkniętym oraz pamiętając, że ciepło właściwe pod stałą objętością jest takie jak we wzorze (26.9). Korzystając ze wzoru na energię średnią układu (energia wewnętrzna układu cząstek układu statystycznego), to względna fluktuacja energii jest napisana przez wzór zależnych od średniej liczby cząstek znajdujących się w układzie wedle sposobu:

(26.16)

Na podstawie (26.16), czym większa jest średnia liczba cząstek, to względne fluktuacje energii w danym układzie są czym mniejsze. Dla średniej liczby cząstek układu dążących do nieskończoności nie ma względnych fluktuacji energii.

Fluktuacja liczby cząstek[edytuj]

Średnia liczba cząstek w układzie według wielkiego zespołu kanonicznego, którego wyprowadzenie napisaliśmy w punkcie (18.5) jest ona równoważna do wspomnianego wzoru na tą średnią wyrażonej poprzez sumę statystyczną "Z", ale nie używają logarytmu naturalnego:

(26.17)

A średnia kwadratu ilości cząstek w układzie wyrazimy przez sumę statystyczną naszego zespołu kanonicznego:

(26.18)

A średnia kwadratu liczby cząstek układu wyrazimy przez sumę statystyczna naszego zespołu kanonicznego, i którą policzymy z definicji wartości średniej kwadratu liczby cząstek (26.18) i średniej liczby cząstek (26.17), jakie układ statystyczny może przyjmować:

(26.19)

Średni kwadrat fluktuacji bezwzględnej ilości cząstek znajdującej się w układzie zamkniętym, które jak wyprowadziliśmy jest zależna liniowo od kwadratu temperatury bezwzględnej układu oraz liniowo od pochodnej cząstkowej średniej liczby cząstek względem temperatury bezwzględnej T panujących w naszym układzie. jest napisana:

(26.20)

Opis ruchów Browna przez Einsteina[edytuj]

Weźmy sobie ruchy Browna, którą będziemy obserwowali w jednym wymiarze, dla której znamy rozkład prawdopodobieństwo. Załóżmy, że w przedziale czasowym (t,t+τ) cząstka zmienia swoje położenie x do x+Δ. Niech mamy N cząstek w układzie, dla których dN cząstek przemieściło się z Δ +Δ+dΔ w przedziale czasowym (t, t+τ). Infinitezymalne prawdopodobieństwo, że takie przesunięcie nastąpiło jest wyrażone poprzez:

(26.21)

Prawdopodobieństwow względem małego przemieszczenia jest funkcją parzystą i spełnia warunek fτ(Δ)=fτ(-Δ) Całkowite prawdopodobieństwo określoną względem dowolnego przemieszczenia jest równe jedności i jest wyrażona:

(26.22)

Liczba cząstek, które w przedziale czasowym τ wyszły z przedziału dx piszemy według następującego przestawienia:

(26.23)

Gdy przerwa dx w powyższym wyrażeniu (26.23) dąży do zera, zatem granica całki prawdopodobieństwa staje się zaniedbywalnie mała, wtedy praktycznie 100% cząstek ucieka z przedziału dx. Liczba cząstek znalezione w dx w chwili t+τ, które były tam w przedziale czasowym τ jest zapisane poprzez wzór:

(26.24)

Rozłóżmy sobie w szereg Taylora prawą stronę względem Δ, a jego lewą stronę względem τ i w ten sposób otrzymujemy tożsamość:

(26.25)

Pierwszy wyraz po prawej stronie jest równa jedności, a dla drugiego wyrazu ze względu na parzystość funkcji gęstości prawdopodobieństwa znika ona, i tylko pozostaje nam trzeci wyraz, bo dalsze wyrazy gwałtowanie maleją do zera dla , zatem w ten sposób otrzymujemy równość:

(26.26)

Stałą występujące we wzorze (26.26), którą oznaczymy przez D, i określimy ją względem wielkości infinitezymalnej wielkości τ, którą możemy przepisać według:

(26.27)

Przy definicji stałej D (26.27) dla równości (26.25) otrzymujemy stąd wniosek, która jest związkiem łączącym pierwszą pochodną koncentracji cząstek względem czasu z drugą pochodną koncentracji względem współrzędnej iksowej, który jest opisem ruchów Browna dany przez Einsteina:

(26.28)

Równanie (26.28) możemy uogólnić z jednego do trzech wymiarach zastępując w powyższych x przez , a całka pojedynczą przez całkę potrójną, i w ten sposób równość (26.24) po dokonaniu pewnych operacji według wcześniejszych wspomnień rozkładając go w szereg Taylora otrzymamy równość:


(26.29)

Pierwsza całka występująca po jego prawej stronie jest równa jedności pomnożonej przez n(x,t), druga całka jest równa zero ze względu na parzystość funkcji gęstości prawdopodobieństwa, zatem dla Δxyz=Δ da nam w ostatecznych rozrachunkach:

(26.30)

Wtedy definicja parametru D dla trzech wymiarów możemy przepisać go na podstawie (26.30) przy pomocy definicji definicji gęstości prawdopodobieństwa napisanej dla trzech wymiarów:

(26.31)

Równanie opisujące w trzech wymiarach możemy przestawić, czyli (26.30), przy definicji parametru D (26.31), w postaci:

(26.32)

Związek współczynnika dyfuzji z temperaturą bezwzględną[edytuj]

Równanie ciągłości napisane w trzech wymiarach napisanej dla koncentracji cząstek powiązane z gęstością prądu piszemy:

(26.33)

Pochodną czasową koncentracji cząstek piszemy przy pomocy równania (26.32), które podstawiamy do tożsamości różniczkowej (26.33) i w ten sposób możemy otrzymać tożsamość na gęstość prądu:

(26.34)

Cząstka poruszająca się w płynie osiąga stacjonarną gęstość dryfu proporcjonalne do działającej siły zewnętrznej, który tutaj przestawiamy poprzez współczynnik proporcjonalności η zwanej współczynnikiem ruchliwości:

(26.35)

Gęstość dryfu musi kompensować z gęstością dryfu , który prowadzi do gradientu gęstości n, czyli te prądy muszą być w równowadze, co możemy napisać równość na tą równowagę wzorem:

(26.36)

Jeśli do wzoru (26.36) podstawimy wzór na gęstość prądu płynu (26.34) i wzór na prędkość stacjonarnego dryfu (26.35) i wykorzystując do tego ostatniego wyrażenie siły poprzez gradient potencjału:

(26.37)

Z drugiej jednak strony gęstość cząstek o energii potencjalnej możemy wyrazić poprzez tą energię i wyrazić jego gradient, który podstawimy je do tożsamości (26.37), koncentrację w zależności od energii potencjalnej cząstek i jego gradient możemy przestawić je w jednej linijce, które są wyrażone w zależności od położenia:

(26.38)
(26.39)

Wzór (26.39) podstawiamy do tożsamości (26.37) i w ten sposób otrzymujemy:

(26.40)

Współczynnik dyfuzji D występujący w prawie (26.32), która jest równaniem dyfuzji wyprowadzonego przez Einsteina, przestawiamy w zależności od temperatury T i współczynnika ruchliwości η:

(26.41)

Szumy Nyquista i wyprowadzenie wzoru Nyquista[edytuj]

Termiczny ruch w metalach powoduje tzw. szumy Nyquista, tzn. są to fluktuację od średnich napięcia i natężenia płynącego prądu w metalach, i dla opornika możemy napisać średnią kwadratu odchylenia napięcia od kwadratu wartości średniej:

(26.42)

Wynik (26.41) wiąże fluktuację napięcia z oporem i jest przykładem twierdzenia flutuacyjno-dyssypacyjnego. Weźmy sobie opornik o długości L i oporze R. Na końcach tego przewodnika jest umieszczony opornik R. W tym przewodniku rozchodzą się fale biegnące w przewodniku R, które są pochodzenia termicznego, które są absorbowane całkowicie bez odbić na obu jego końcach. Częstotliwość modu podstawowego, a także częstotliwość o modzie n, dla rozważanego przewodnika o długości L, przestawiamy je w jednej linijce:

(26.43)
(26.44)

Wiedząc, że częstotliwość kołowa ωn jest związana z częstością νn poprzez wzór ωn=2πνn, wtedy energię modu o numerze n możemy określić według rozkładu Bosego-Einsteina sposobem:

(26.45)

Tutaj będziemy stosowali warunki na częstości kołowe dla rozkładu (26.45), dla którego stosując ten warunek, otrzymujemy:

(26.46)
(Rys. 26.1) Schemat elektryczny obwodu elektrycznego służący do wyprowadzenia prawa Nyquista

Widzimy, że energia En nie zależy od numeru modu, zatem jeśli w szerokości Δν na podstawie (26.44) istnieje Δν/ν0 modów, wtedy całkowita energia w naszym pasmie jest iloczynem liczby modów przez kBT, i jeśli jednocześnie zastosujemy wzór (26.43) do wzoru na całkowitą energię, wtedy tą wielkość określamy:

(26.47)

W przewodniku rozchodzą się dwie fale biegnące w przeciwnych kierunkach normalna i odbita, wtedy całkowity czas z jakim jedną z tych fal biegnących przebyła odległość L, czyli przez jeden przewodnik, możemy obliczyć z definicji prędkości:

(26.48)

Całkowita energia absorbowana przez jeden oporniki R, pochodząca od jednej fali biegnących w jednostce czasu określamy jako stosunek (26.47) przez podwojoną prędkość czasu (26.48), jest ona wyrażona jako:

(26.49)

We wzorze (26.49) rozpatrywaliśmy tylko jedną falę biegnącą w jednym kierunku. Z definicji energii traconej przez opornik w jednostce czasu określamy przez iloczyn kwadratu natężenia prądu i oporności odbiornika, co stosując jego definicję do (26.49):

(26.50)

Według rysunku obok napięcie mierzone na obu opornikach R jest wyrażone poprzez U=2IR, co stąd wynika, że I=U/2R, wtedy całkowita moc wydzielana przez oporność jest przepisana poprzez:

(26.51)

Wzór (26.51) możemy podstawić do tożsamości (26.50), i na samym końcu otrzymujemy wzór (26.42), który jest wzorem Nyquista.