Stationary $1/f^α$ noise in discrete models of the Kardar-Parisi-Zhang class

Dit artikel toont aan dat in discrete modellen van de Kardar-Parisi-Zhang-klasse de hoogtefluctuaties op een vast punt wide-sense stationair zijn en een $1/f^{5/3}$-ruispatroon vertonen, wat de toepasbaarheid van de Wiener-Khinchin-stelling bevestigt.

De kern

Het Geluid van een Groeiende Stad: Waarom Ruimtelijke Ruis een Stabiel Ritme Heeft

Stel je voor dat je naar een drukke bouwplaats kijkt waar duizenden arbeiders tegelijkertijd muren optrekken. Ze gooien bakstenen op elkaar, maar niet perfect netjes. Soms valt een steen scheef, soms wordt hij op een piek gelegd, soms in een kuil. Na verloop van tijd ontstaat er een ruw, hobbelig landschap. In de natuurkunde noemen we dit een "groeigrensvlak" (zoals een laagje verf dat droogt, of bacteriën die een oppervlak bedekken).

De auteurs van dit paper, Rahul en Avinash, hebben gekeken naar de ruis in dit proces. Ze keken niet naar de hele stad, maar naar één specifieke plek op de muur: hoe hoog is die ene baksteen op dit moment, en hoe verandert die hoogte in de loop van de tijd?

Hier is wat ze hebben ontdekt, vertaald naar alledaags taal:

1. Het mysterie van het "1/f"-geluid

In de wereld van geluid en signalen is er een heel bekend fenomeen dat klinkt als een diep, grommend ruisgeluid. Dit wordt "1/f-ruis" genoemd (of "roze ruis"). Je hoort dit in veel dingen: in de muziek van Bach, in het gedrag van sterren, of zelfs in je eigen hartslag. Het is een soort "perfecte chaos" die overal voorkomt.

Vroeger dachten wetenschappers dat dit soort ruis in groeiende systemen (zoals onze bouwplaats) niet stabiel was. Het idee was: "Als je lang genoeg kijkt, wordt het steeds chaotischer en stopt het met een vast ritme." Het was alsof je naar een radio luisterde die steeds harder kraakt en nooit tot rust komt.

2. De grote ontdekking: Het is wel stabiel!

Rahul en Avinash hebben gekeken naar kleine systemen (een kleinere bouwplaats). Ze hebben ontdekt dat als je lang genoeg wacht, deze kleine systemen wel een stabiele toestand bereiken.

• De analogie: Stel je voor dat je een badkamer vol water hebt. Als je de kraan heel langzaam openzet en er zijn geen gaten in de bodem, zal het water uiteindelijk een vast niveau bereiken. Het blijft daar dan "stilstaan" (in een dynamisch evenwicht), zelfs als er nog steeds waterbellen op en neer dansen.
• In hun onderzoek zagen ze dat de hoogte van de muur op één plek, na een lange tijd, begint te gedragen als een stabiel ritme. Het is niet meer een chaotische storm, maar een stabiele zee met golven.

3. De "Vergeten" Regel (Wiener-Khinchin)

Er is een beroemde regel in de fysica, de Wiener-Khinchin-stelling. Deze regel zegt: "Als je geluid stabiel is, dan kun je het geluid (ruis) precies voorspellen door te kijken naar hoe het in de tijd samenhangt."

Vroeger dachten wetenschappers dat deze regel niet werkte voor groeiende systemen, omdat ze dachten dat die systemen nooit stabiel werden.
Maar deze paper zegt: "Wacht even! Voor kleine systemen werkt de regel wél!"
Ze hebben bewezen dat als je de "golven" op je muur meet, die golven een vast patroon hebben. Je kunt dus het geluid van de muur voorspellen door te kijken naar hoe de golven in de tijd op elkaar lijken.

4. De "Snelheidsrem" (De ondergrens)

Een van de coolste dingen die ze zagen, is dat het geluid niet oneindig laag kan worden.

• De analogie: Stel je voor dat je een radio hebt die je kunt afstemmen op heel lage tonen. Bij een normaal radioapparaat zou je bij heel lage tonen alleen maar stilte horen. Maar bij deze "groeiende muur" is er een ondergrens.
• Als je de frequentie (hoe vaak de muur verandert) te laag maakt, stopt het geluid niet met afnemen, maar blijft het op een vaste, constante niveau.
• Dit komt omdat de muur een "geheugen" heeft. De golven op de muur hebben een maximale levensduur. Na die tijd "vergeten" ze elkaar. Dit creëert een soort "muur" in het geluidsspectrum waaronder het geluid niet verder zakt.

5. Waarom is dit belangrijk?

Vroeger dachten we dat dit soort systemen (zoals het groeien van bacteriën of het afzetten van dunne lagen in de chipindustrie) te chaotisch waren om te begrijpen met simpele regels.

Dit paper zegt: "Nee, ze zijn niet chaotisch, ze zijn gewoon even geduldig."
Als je een klein systeem hebt en je wacht lang genoeg, vinden ze een ritme. Ze gedragen zich als een stabiel muziekinstrument in plaats van als een vallende steen.

Samenvattend in één zin:
De auteurs hebben bewezen dat kleine, groeiende oppervlakken (zoals een laagje verf of een baksteenmuur) na verloop van tijd een stabiel ritme vinden, waardoor we de "ruis" die ze maken kunnen begrijpen en voorspellen met bestaande wiskundige regels, zolang we maar niet naar oneindig grote systemen kijken.

Technische samenvatting

Titel: Stationaire $1/f^\alpha$-ruis in discrete modellen van de Kardar-Parisi-Zhang-klasse

Auteurs: Rahul Chhimpa en Avinash Chand Yadav (Banaras Hindu University, India)

---

1. Probleemstelling en Achtergrond

Het artikel adresseert een fundamentele tegenstrijdigheid in de studie van niet-evenwichtssystemen, specifiek binnen de Kardar-Parisi-Zhang (KPZ) universaliteitsklasse. Deze klasse beschrijft de groei van ruwe interfaces (zoals dunne films of bacteriële kolonies) en vertoont kenmerkende ruimtetijdschaling.

• De Kernvraag: Hoe kan men de aard van de tijdsfluctuaties in deze systemen karakteriseren? Bestaande literatuur (o.a. Takeuchi, Henkel et al.) suggereerde dat deze fluctuaties niet-stationair zijn. Dit werd gebaseerd op het ontbreken van een ondergrens (lower cutoff) in het vermogensspectraaldichtheid (PSD) en het vertonen van "aging"-gedrag in de autocorrelatiefunctie voor oneindig grote systemen.
• Het Conflict: Voor stationaire processen geldt de Wiener-Khinchin-stelling, die stelt dat het PSD de Fourier-getransformeerde is van de autocorrelatiefunctie. Echter, als processen niet-stationair zijn (zoals bij unrestricted Brownse beweging), is deze stelling niet direct toepasbaar en ontbreekt er vaak een ondergrensfrequentie in het spectrum.
• Het Doel van het Artikel: De auteurs onderzoeken of het mogelijk is om een stationaire toestand te bereiken in discrete KPZ-modellen binnen een eindige, maar voldoende lange observatietijd, vooral voor kleine systemen. Ze willen bepalen of de Wiener-Khinchin-stelling hier toepasbaar is en of de bekende $1/f^\alpha$-scaling ($\alpha = 5/3$) consistent is met stationair gedrag.

2. Methodologie

De auteurs gebruiken numerieke simulaties van vier discrete celautomaten-modellen die tot de KPZ-universaliteitsklasse behoren:

1. Kim-Kosterlitz (KK) model: Een model met een hellingbeperking.
2. Ballistic Deposition (BD) model: Deeltjes vallen verticaal en hechten aan het eerste contactpunt.
3. Etching Model (EM): Een depositieversie van een etsmodel.
4. Single-Step (SS) model: Een model met groef-initieel voorwaarde.

Simulatie-instellingen:

• Systeem: Eendimensionale roosters met periodieke randvoorwaarden, grootte $L$.
• Tijdsunit: Monolaag-eenheden (1 Monte Carlo-stap = $L$ depositiepogingen).
• Signaal: De fluctuaties in de hoogte op een vaste positie $x$, gedefinieerd als $\delta h_L(t) = h(x,t) - \bar{h}(t)$, waarbij $\bar{h}(t)$ de ruimtelijk gemiddelde hoogte is.
• Observatietijd: De simulatietijd $T$ wordt gekozen als $T \gg L^z$ (waarbij $z$ de dynamische exponent is), om te garanderen dat het systeem de stationaire toestand heeft bereikt en transiënten zijn weggefilterd.

Analyse:

• Berekening van de twee-tijds autocorrelatiefunctie $C(\tau, L) = \langle \delta h_L(t) \delta h_L(t+\tau) \rangle$.
• Berekening van het vermogensspectraaldichtheid (PSD) $S(f, L)$ via Fast Fourier Transform (FFT).
• Toepassing van Finite-Size Scaling (FSS) om de schalingsgedrag te analyseren en kritieke exponenten te bepalen.

3. Belangrijkste Resultaten

A. Autocorrelatie en Stationariteit

• De auteurs vinden dat voor eindige systemen de autocorrelatiefunctie $C(\tau, L)$ na een correlatietijd $\tau_c \sim L^z$ naar nul daalt.
• De functie vertoont niet-exponentieel verval en voldoet aan dynamische schaling:
  $$C(\tau, L) \sim L^{2\chi} F_C(\tau/L^z)$$
  waarbij $\chi$ de ruwheidsexponent is.
• Cruciaal is dat de correlatiefunctie alleen afhangt van het tijdsverschil $\tau = |t_1 - t_2|$, wat aantoont dat het proces wijd-sensie stationair (wide-sense stationary) is binnen de simulatieduur. Dit weerlegt de eerdere aanname dat het proces inherent niet-stationair is voor alle tijden; het is alleen niet-stationair als de observatietijd korter is dan de divergerende correlatietijd van een oneindig systeem.

B. Vermogensspectrum en $1/f^\alpha$-scaling

• Het PSD toont een duidelijke ondergrensfrequentie $f_c \sim L^{-z}$. Onder deze frequentie is het vermogen constant.
• In het niet-triviale frequentiegebied ($L^{-z} \ll f \ll 1/2$) vertoont het spectrum een machtswet-afname:
  $$S(f) \sim f^{-\alpha}$$
• De numerieke analyse levert een spectrale exponent op van $\alpha \approx 5/3$.
• Door FSS toe te passen op zowel de correlatiefunctie als het PSD, vinden de auteurs een perfecte data-collapse. Dit bevestigt dat de schalingsrelaties consistent zijn.

C. Relatie tussen Exponenten

De auteurs leiden een theoretische relatie af tussen de spectrale exponent $\alpha$, de dynamische exponent $z$ en de ruwheidsexponent $\chi$:
$$\alpha = 1 + \frac{2\chi}{z}$$
Voor de KPZ-klasse in 1+1 dimensies zijn de theoretische waarden $\chi = 1/2$ en $z = 3/2$. Dit resulteert in:
$$\alpha = 1 + \frac{2(1/2)}{3/2} = 1 + \frac{2}{3} = \frac{5}{3}$$
Dit komt exact overeen met de numeriek gevonden waarde.

4. Belangrijkste Bijdragen

1. Bevestiging van Stationariteit: Het artikel bewijst dat KPZ-systemen voor eindige systemen wel degelijk een stationaire toestand bereiken binnen een lange observatietijd. Dit is een correctie op eerdere interpretaties die het proces als inherent niet-stationair beschouwden.
2. Toepasbaarheid van Wiener-Khinchin: Omdat het proces stationair is, is de Wiener-Khinchin-stelling geldig. De auteurs tonen aan dat het PSD inderdaad de Fourier-getransformeerde is van de autocorrelatiefunctie, wat de theoretische consistentie van de $1/f^\alpha$-ruis in deze systemen bevestigt.
3. Universeel Gedrag: De resultaten zijn robuust en gelden voor meerdere discrete modellen (KK, BD, EM, SS) binnen dezelfde universaliteitsklasse, wat aantoont dat het gedrag fundamenteel is en niet model-specifiek.
4. Verklaring voor Experimentele Beperkingen: Het artikel legt uit waarom experimenten vaak geen ondergrensfrequentie zien: in grote systemen (of oneindige limieten) divergeert de correlatietijd ($T \sim L^z$), waardoor de ondergrensfrequentie ($f_c \sim L^{-z}$) onbereikbaar laag wordt binnen de beschikbare observatietijd.

5. Conclusie en Significatie

De studie concludeert dat de fluctuaties in de KPZ-universaliteitsklasse wijd-sensie stationair zijn, mits de observatietijd groter is dan de systeemgrootte-afhankelijke correlatietijd. De waargenomen $1/f^{5/3}$-ruis is een direct gevolg van deze stationaire dynamica en de onderliggende schalingswetten.

Deze bevindingen zijn significant omdat ze:

• Een brug slaan tussen numerieke simulaties, theoretische schalingsargumenten en experimentele observaties van ruwe interfaces.
• De validiteit van de Wiener-Khinchin-stelling herstellen voor dit type niet-evenwichtssystemen.
• Een verklaring bieden voor de schijnbare non-stationariteit in eerdere studies: dit is een artefact van het observeren van systemen voordat ze de stationaire toestand hebben bereikt of van het observeren van systemen die te groot zijn voor de beschikbare tijd.

Het werk biedt een helder kader voor het interpreteren van $1/f$-ruis in complexe groeiende systemen en benadrukt het belang van het onderscheiden tussen transiënt en stationair gedrag in niet-evenwichtsfysica.