Przeskocz do treści

Delta mi!

Twierdzenie o powracaniu i pewne zagadki nierównowagowej mechaniki statystycznej (I)

Krzysztof Turzyński i Henryk Żołądek

o artykule ...

  • Publikacja w Delcie: sierpień 2015
  • Publikacja elektroniczna: 31-07-2015
  • Autor: Henryk Żołądek
    Afiliacja: Instytut Matematyki, Wydział Matematyki, Informatyki i Mechaniki, Uniwersytet Warszawski
  • Wersja do druku [application/pdf]: (767 KB)

Wymienione w tytule tematy uchodzą za trudne. Tym ważniejsze jest zacząć od rzeczy prostych i powszechnie znanych...

obrazek

Weźmy pewien zbiór o skończonej liczbie elementów i jego różnowartościowe przekształcenie na siebie, czyli permutację. Dla lepszego efektu wizualnego będziemy sobie wyobrażać zbiór czarnych lub białych kwadratów wypełniających pola kwadratowej tablicy 32 ×32 oraz permutację, dla której w i-tym wierszu i | j-tej kolumnie tablicy umieszczamy kwadrat poprzednio znajdujący się w wierszu o numerze |(i + j)mod 32 oraz kolumnie o numerze (2i + j) mod 32. Wyniki wielokrotnego złożenia takiej permutacji dla pewnego wyjściowego ułożenia kwadratów przedstawione są na rysunku 1: po wykonaniu 32 iteracji otrzymujemy taki sam wzór jak na początku.

obrazek

Rys. 1

Rys. 1

Czy podobne efekty - powracanie zbioru do stanu wyjściowego po wielokrotnym wykonaniu pewnego przekształcenia - możemy zaobserwować także dla zbiorów nieskończonych? Aby odpowiedzieć na to, intuicyjnie zrozumiałe, ale nie do końca precyzyjne pytanie, przypatrzmy się jeszcze trzem dość prostym przykładom.

Obrót na okręgu

obrazek

Rys. 2 Obrót na okręgu

Rys. 2 Obrót na okręgu

Niezbyt skomplikowanym przykładem zbioru nieskończonego jest odcinek |[0, 1]. Żeby uniknąć kłopotów związanych z tym, iż odcinek ma końce, możemy utożsamić jego lewy koniec z prawym - otrzymamy wtedy okrąg długości 1, którego punkty ponumerowane są liczbami rzeczywistymi od 0 do 1 (tyle że jedynkę skleiliśmy z zerem). Dla takiego okręgu oraz α ∈(0, 1) określimy przekształcenie R α następującym wzorem:

Rα(x) = x +α − ⌊x +α ⌋.

Symbol części całkowitej ⌊⋅⌋, występujący w tej definicji, zapewnia, że gdy liczba x + α przekracza 1, przekształcenie Rα umieszcza ją z powrotem na odcinku [0,1] w takiej samej odległości od 0, w jakiej liczba |x + α znajduje się "na prawo" od 1. Z tego względu przekształcenie |Rα możemy wyobrażać sobie jako przesuwanie punktów naszego okręgu o odcinek α (Rys. 2). Jeżeli α jest liczbą wymierną, którą możemy zapisać w postaci ułamka nieskracalnego |pq, nietrudno będzie przekonać się, że po q-krotnym wykonaniu przekształcenia R α każdy punkt okręgu wróci na swoje wyjściowe miejsce. Możemy tę własność powracania wyrazić w mądry sposób, mówiąc, że dla dowolnego punktu x0, leżącego na okręgu, jego dodatnia orbita

𝒪+(x ) = {R (x ),R2(x ),R3(x ),...} 0 α 0 α 0 α 0

przechodzi przez punkt |x0 (oznaczamy R2α(x0) = Rα(R α(x0)) = R2α(x0), |R3(x ) = R (R (x )) = R (x ) α 0 α 2α 0 3α 0 itd.). Dla α niewymiernego tak sformułowana własność powracania nie zachodzi, nie sposób się jednak oprzeć wrażeniu, że w wyniku wielokrotnego przekształcania punkt |x0 powróci "w pobliże" swego położenia wyjściowego. Możemy to wyrazić następującym stwierdzeniem

Twierdzenie (Własność powracania). Dla dowolnego punktu początkowego x0 i dowolnego jego otoczenia U orbita |𝒪+(x0) przecina |U.

Używając mądrych słów, moglibyśmy powiedzieć, że R | α ma własność topologicznej tranzytywności oraz że jest przekształceniem prawie okresowym.

Skośny obrót torusa

obrazek

Rys. 3 Sklejenie przeciwległych boków kwadratu zgodnie z kierunkiem strzałek daje torus.

Rys. 3 Sklejenie przeciwległych boków kwadratu zgodnie z kierunkiem strzałek daje torus.

Innym prostym przykładem zbioru nieskończonego jest kwadrat o boku 1. Żeby nie kłopotać się jego brzegiem, wykonamy operację podobną do sklejania końców odcinka, tyle że teraz należy skleić przeciwległe boki tak, by kierunki strzałek na rysunku 3 się zgadzały. W ten sposób otrzymujemy torus. Punkty tego torusa możemy opisać dwiema współrzędnymi |(x,y) odpowiednich punktów kwadratu, przy czym |0⩽ x,y < 1. Na torusie możemy zdefiniować następujące przekształcenie:

T(x,y) = (x,x + y− ⌊x +y ⌋)

przyporządkowujące danemu punktowi torusa pewien inny punkt. Ponieważ współrzędna x nie zmienia się przy tym przekształceniu, działanie |T sprowadza się do przesuwania punktów odcinka x = const o |x "w górę" (lub, równoważnie, obrotu odpowiedniego okręgu o x), a zatem przekształcenia rozważanego w poprzednim przykładzie. Wynika stąd, że sformułowana powyżej własność powracania jest spełniona także dla przekształcenia |T.

Siodło-węzeł na okręgu

obrazek

Rys. 4 Wykres funkcji

T x x -1sin2 πx 8

oraz graficzny sposób znajdowania wartości kolejnych iteracji tego przekształcenia.

Rys. 4 Wykres funkcji

 1 T x x - sin2 πx 8

oraz graficzny sposób znajdowania wartości kolejnych iteracji tego przekształcenia.

Po tej rozgrzewce możemy przystąpić do analizy nieco trudniejszego przykładu. Rozważmy przekształcenie T [0,1) [0,1) określone wzorem:

T (x) = x + 1-sin2(π x). 8

Możemy uważać je za różnowartościowe przekształcenie okręgu na siebie. Jak wyglądają orbity tego przekształcenia? Najprościej przekonać się o tym, rysując wykres T (x) i pamiętając, że |Tn(x) = T(T n−1(x)). Wynik iterowania przekształcenia T dla pewnego dowolnie wybranego punktu x0 przedstawiony jest na rysunku 4. Widzimy, że dla dowolnego punktu początkowego x0 ≠ 0 orbita 𝒪+(x0) składa się z wyrazów ciągu |xn = T n(x0), który jest zbieżny do |... No właśnie, powiedzielibyśmy odruchowo, że jest zbieżny do 1, ale utożsamiliśmy jedynkę z zerem. Widzimy zatem, że punkt 0 ma ciekawą własność: odpycha punkty leżące odeń na prawo, a przyciąga te leżące na lewo. Przez analogię do sytuacji, z jakimi spotykamy się przy przekształceniach dwuwymiarowych, możemy nazwać go punktem typu siodło-węzeł.

Łatwo zauważyć, że jedynie jednopunktowa orbita 𝒪+(0) = {0} punktu 0 ma rozważaną poprzednio własność powracania. Dla x ≠ 0 0 orbita  + |𝒪 (x0) rozmija się z dostatecznie małym otoczeniem punktu |x0.

Miara

Na podstawie trzech przykładów doszliśmy do wniosku, że jedne przekształcenia mają własność powracania, a inne jej nie mają, przynajmniej w wersji, jaką podaliśmy. Czy nie da się lepiej? Czy można sformułować jakieś ogólniejsze twierdzenie dotyczące powracania? Okazuje się, że można. Niezbędne będzie do tego pojęcie miary zbioru, a dokładniej miary probabilistycznej.

Sama nazwa wskazuje na to, że pojęcie miary wprowadzono po to, by określać, jak "duży" jest zbiór. Ścisła definicja miary jest dość techniczna, choć intuicyjna, ograniczymy się zatem do podania własności miary. Otóż miara na zbiorze X jest funkcją określoną na pewnej rodzinie podzbiorów |X o odpowiednich własnościach. Funkcja ta, którą będziemy oznaczać przez |µ, spełnia:

  • µ (∅) = 0,
  • dla |A zachodzi µ(B ∖ A)
  • jeśli |A1, są parami rozłączne, to |µ( Ai)

Miarę nazywamy probabilistyczną, jeśli )=1. µ(X

Nikt nie zdziwi się chyba, że dobrą miarą dla odcinka jest jego długość, miarą dla (otwartego) podzbioru płaszczyzny - jego pole. Rozszerzenie tej intuicji na bardziej złożone przypadki prowadzi do pojęcia miary Lebesgue'a, w które nie będziemy się bardziej szczegółowo zagłębiać. Istnieje wszakże jeszcze jeden prosty przykład miary. Każda gospodyni prowadząca tzw. salon wie, że o jego randze decyduje nie ogólna liczba gości, lecz liczba gości znamienitych. Najprostszym odpowiednikiem takiego pomysłu jest wyróżnienie w zbiorze | X pewnego punktu | x0 i powiedzenie, że |µ(A) wtedy i tylko wtedy, gdy x0 ∈A, a w przeciwnym przypadku |µ(A) Taką miarę nazywamy miarą Diraca i oznaczamy δ x0 ; jej uogólnienie na dowolną skończoną liczbę wyróżnionych punktów jest oczywiste.

Dla przekształcenia X T X miara µ jest niezmiennicza, jeśli  −1 |µ(T (A)) Dla obrotu na okręgu, Rα, oczywistym przykładem miary niezmienniczej jest miara Lebesgue'a, w szczególności długość odcinka nie zmienia się przy równomiernym przesunięciu wszystkich punktów. Ale dla obrotu postaci |R p~q są jeszcze inne miary niezmiennicze - postaci

 1- µ = q (δx0 + δx0+1q + ...+ δx0+ qq 1). (*)

Jeśli obrót jest niewymierny, jedyną miarą niezmienniczą jest miara Lebesgue'a. Spełnia ona wówczas następujące twierdzenie

Twierdzenie (Własność ergodyczności). Jeśli |T−1(A) to µ (A) lub µ (A) ; jedyne podzbiory niezmiennicze - zbiór pusty i cały okrąg - są trywialne w sensie miarowym.

W przypadku skośnego obrotu torusa naturalną miarą niezmienniczą dla T jest miara Lebesgue'a: prostokąt |[a,b]× [c,d] jest przekształcany na równoległobok o wysokości |b− a i długości podstawy d − c. Ale istnieją też miary postaci

µ([a,b] ×[c,d]) = δ ([a,b]) ⋅µ ([c,d]), x0 1

gdzie µ 1 jest miarą Lebesgue'a na okręgu (miary "skupione" na okręgu danym przez x = x 0 ) lub, dla |x = p/q, 0 miarą daną wzorem |(∗).

Interesują nas często "typowe" własności, spełnione dla prawie każdego obiektu. Istnieją dwa sposoby rozumienia tego zwrotu: "prawie każdy" oznacza każdy z wyjątkiem elementów skończonego zbioru punktów lub każdy z wyjątkiem elementów pewnego zbioru miary zero. Będziemy się dalej posługiwać drugim znaczeniem tego zwrotu. Możemy teraz sformułować pewne bardzo ważne twierdzenie.

obrazek

Twierdzenie (Poincarégo o powracaniu). Niech przekształcenie XT X zachowuje probabilistyczną miarę µ . Wtedy dla każdego zbioru A o mierze (A) µ prawie każdy punkt |x0∈ A powraca do A, tzn. dla pewnego n | > 0 zachodzi |Tn(x0) ∈A.

Dowód. Niech B będzie zbiorem tych punktów A, które nigdy nie powracają do |A. W szczególności, dla n | > 0 i x ∈ B mamy |Tn(x) ∉B, czyli |x∉ T −n(B), a więc |T−n(B) ∩ B = ∅. Stąd

T−1(B) ∩T −1− n(B) = T −1(T −n(B) ∩ B) = T −1(∅ ) =∅ .

Analogicznie, dla dowolnego |k∈ N zachodzi |T−k(B) ∩ T−k−n(B) = ∅. Zatem zbiory postaci |T−n(B) są parami rozłączne oraz mają taką samą miarę. Gdyby zachodziło µ(B) > 0, to z własności c) w definicji miary mielibyśmy |µ( T −n(B)) = Pn∞ 0µ(B) = ∞ , co byłoby sprzeczne z tym, że | T−n(B) ⊂ X oraz )=1. µ(X Musi zatem zachodzić µ(B) = 0.


Teraz możemy na nowo rozważyć siodło-węzeł na okręgu. Zero musi się znajdować w każdym podzbiorze A o dodatniej mierze niezmienniczej. Ale wtedy zwrot prawie każdy punkt x0 ∈ A oznacza, że x0 | = 0. Ten punkt nie opuszcza zbioru A.

Przekształcenie piekarza

obrazek

Rys. 5 Przekształcenie piekarza

Rys. 5 Przekształcenie piekarza

Omówione na wstępie przykłady ilustrują dość szczególne mechanizmy powracania. Warto uzupełnić tę listę jeszcze jednym przykładem. Rozważmy przekształcenie kwadratu jednostkowego dane wzorem

 ⎧ y 1 T(x,y) = ⎪⎪⎨(2x, 2) dla x < 2, ⎪⎪(2x − 1, y+1) dla x ⩾ 1. ⎩ 2 2

Działanie tego przekształcenia pokazano na rysunku 5: spłaszcza i rozciąga ono wyjściowy kwadrat do prostokąta o podstawie 2 i wysokości 1 2, następnie rozcina uzyskany prostokąt na dwa prostokąty o podstawie 1, a w końcu układa te prostokąty jeden na drugim, tak że znowu powstaje kwadrat jednostkowy. Z uwagi na to rozciąganie i ściskanie (operacje doskonale znane każdemu, kto kiedykolwiek wyrabiał ciasto drożdżowe) przekształcenie piekarza zachowuje (dwuwymiarową) miarę Lebesgue'a. Okazuje się, że przekształcenie to ma następującą własność

Twierdzenie (Własność mieszania).

 −n lni m∞ µ(T (A)

czyli po dużej liczbie iteracji przekształcenia T obraz |Tn(A) dowolnego zbioru A | rozkłada się prawie równomiernie po całym kwadracie.

Własność mieszania jest silniejsza od własności ergodyczności; ta druga wynika z pierwszej. W szczególności obrót niewymierny na okręgu nie jest mieszający, chociaż jest ergodyczny.

Podsumowanie

Chociaż dowód twierdzenia Poincarégo o powracaniu wydaje się oczywisty i mało subtelny, twierdzenie to wywoływało (i najwidoczniej nadal wywołuje) dyskusje natury filozoficznej. Są one najczęściej związane z próbami zastosowania go do układów złożonych z dużej liczby cząstek - standardowym przykładem jest gaz skupiony w połówce zbiornika, z którego usuwamy przegrodę. Tymi tematami zajmiemy się w drugiej części artykułu.