Estimarea disipației de pe traiectorii staționare unice

Sursa: Edgar Roldan, Juan M.R. Parrondo (2010) Estimating dissipation from single stationary trajectories, https://doi.org/10.48550/arXiv.1004.2831 (Rezumat)

Articolul analizează o problemă centrală din termodinamică și mecanica statistică: cum poate fi estimată disiparea, adică producerea de entropie, pornind doar de la observarea unei singure traiectorii staționare a unui sistem. Autorii pornesc de la ideea că ireversibilitatea unui proces fizic este legată direct de producerea de entropie. Cu alte cuvinte, dacă o evoluție temporală poate fi deosebită de aceeași evoluție rulată invers în timp, această asimetrie conține informații despre energia disipată și despre faptul că sistemul se află sau nu într-o stare de neechilibru.

Instrumentul matematic principal folosit este entropia relativă, sau distanța Kullback-Leibler, care măsoară cât de diferite sunt două distribuții de probabilitate. În contextul articolului, cele două distribuții comparate sunt probabilitatea de apariție a unei traiectorii înainte în timp și probabilitatea traiectoriei inverse. Dacă aceste probabilități sunt identice, procesul este reversibil din punct de vedere statistic; dacă diferă, procesul produce entropie. Autorii arată că, pentru stări staționare de neechilibru, rata producerii de entropie poate fi legată de entropia relativă dintre traiectoriile directe și cele inversate. Chiar și atunci când observatorul nu cunoaște mecanismul fizic intern al sistemului, simpla analiză statistică a seriilor temporale poate oferi o limită inferioară pentru disipare.

Un aspect important al articolului este faptul că informația completă despre sistem nu este întotdeauna necesară. Dacă se observă doar o parte a variabilelor relevante, entropia relativă calculată din datele disponibile nu mai oferă neapărat valoarea exactă a producerii de entropie, ci o limită inferioară. Aceasta este totuși foarte valoroasă, deoarece permite detectarea neechilibrului chiar și în situații în care fluxurile obișnuite măsurabile lipsesc sau sunt nule. Autorii subliniază astfel o legătură profundă între fizica disipării și statistica seriilor temporale: o mărime fizică, entropia produsă, poate fi estimată din proprietăți pur statistice ale datelor.

Pentru a face metoda aplicabilă, articolul discută estimatori practici ai entropiei relative. Prima metodă constă în numărarea directă a șirurilor de lungime finită dintr-o traiectorie discretă și compararea frecvenței lor cu frecvența șirurilor inverse. Pentru procese Markoviene, această metodă este eficientă, deoarece este suficientă analiza tranzițiilor dintre perechi de stări. Pentru procese mai generale, autorii definesc estimatori bazați pe blocuri de lungime mai mare și folosesc o extrapolare pentru a aproxima limita la blocuri foarte lungi. Ei menționează și o strategie bazată pe parsarea șirurilor, inspirată de algoritmul de compresie Lempel-Ziv, dar se concentrează în principal pe metoda distribuțiilor empirice.

Ca studiu de caz, autorii folosesc un „ratchet” discret cu potențial intermitent, un model de particulă care se deplasează într-o rețea unidimensională sub acțiunea unui potențial periodic asimetric care se activează și se dezactivează. Schema modelului este prezentată în figura 1, unde particula poate face salturi între mai multe stări, iar potențialul comută cu o anumită rată. În cazul în care sunt observate toate variabilele sistemului — atât poziția particulei, cât și starea potențialului — dinamica este Markoviană, iar entropia relativă calculată din tranziții coincide cu disiparea medie pe pas. Figura 2 arată această concordanță pentru potențiale slabe sau moderate, dar evidențiază și dificultăți de eșantionare pentru potențiale mari, deoarece unele salturi rare, dar importante, pot lipsi dintr-o traiectorie finită.

Situația devine mai interesantă atunci când se observă doar poziția particulei, nu și starea potențialului. În acest caz, procesul observat nu mai este în general Markovian, iar estimarea disipării este mai dificilă. Figura 3 arată că estimările bazate pe blocuri mai lungi de observații reușesc să distingă clar cazul de echilibru de starea staționară de neechilibru și reproduc corect ordinul de mărime al disipării, deși o subestimează. Această subestimare apare deoarece observarea parțială ascunde o parte din ireversibilitatea reală a sistemului.

Un rezultat deosebit de relevant apare în analiza cu forță externă, prezentată în figura 4. În apropierea forței de oprire, fluxul spațial mediu dispare, astfel încât metodele clasice bazate doar pe curenți observabili ar sugera greșit absența disipării. Totuși, entropia relativă calculată din blocuri de lungime mai mare detectează în continuare caracterul de neechilibru al procesului. Aceasta demonstrează forța metodei propuse: ea poate identifica disiparea chiar și atunci când nu există fluxuri măsurabile evidente.

În concluzie, articolul arată că asimetria temporală a unei serii staționare conține informații directe despre producerea de entropie. Metoda propusă nu necesită cunoașterea detaliată a mecanismului fizic care generează datele și poate funcționa chiar cu informație parțială, oferind cel puțin o limită inferioară pentru disipare. Autorii sugerează aplicații importante în sisteme biologice, unde metoda ar putea ajuta la distingerea proceselor active de cele pasive și chiar la estimarea consumului de ATP. În sens mai larg, lucrarea extinde ideea principiului lui Landauer, arătând că ireversibilitatea logică sau statistică observată într-un comportament măsurabil impune un cost termodinamic minim.