Chaos is not an invariant of fractional-order systems: the clock shapes orbits, attractors and bifurcations

Technische Samenvatting: "Chaos is not an invariant of fractional-order systems: the clock shapes orbits, attractors and bifurcations"

Probleemstelling
In de kwalitatieve theorie van gewone differentiaalvergelijkingen (systemen met gehele orde) is een fundamenteel principe, vastgesteld sinds Poincaré, dat het vermenigvuldigen van een vectorveld met een positieve toestandsafhankelijke functie (een "klok") alleen de snelheid verandert waarmee banen worden doorlopen, maar niet de geometrie van de banen zelf. Deze herparametrisatie-invariantie impliceert dat fasenportretten, limietcycli, attractoren en bifurcatiestructuren eigenschappen zijn van het vectorveld alleen, onafhankelijk van de tijdswet.

Fractionele systemen (specifiek gedefinieerd via de Caputo-afgeleide) bezitten echter een niet-lokale geheugenkern die het verleden van de geschiedenis van de toestand weegt. Het onderzoek onderzoekt of deze klassieke herparametrisatie-invariantie standhoudt voor deze systemen. Specifiek wordt de vraag gesteld: als een fractaal systeem $D^\alpha x = f(x)$ wordt gewijzigd naar $D^\alpha x = g(x)f(x)$ (waarbij $g(x)$ een positieve toestandsafhankelijke functie is), beschrijven de oplossingen dan dezelfde geometrische banen, waarbij ze slechts met verschillende snelheden worden doorlopen? De auteur stelt dat voor fractale orden $\alpha \in (0, 1)$ deze invariantie faalt.

Methodologie
De studie maakt gebruik van een combinatie van analytische expansie en numerieke verificatie:

1. Analytische Afleiding: De auteur analyseert het kortetermijngedrag van oplossingen aan de Caputo-differentiaalvergelijking $D^\alpha x = F(x)$ met behulp van de equivalente Volterra-integraalvergelijking. Er wordt een Taylor-achtige expansie uitgevoerd in machten van $t^\alpha$ voor zowel het referentiedisysteem ($F=f$) als het geklokte systeem ($F=gf$), waarbij de oplossing $x(t)$ en de gereparametriseerde referentie-oplossing met de getransformeerde tijd $s(t) = \int_0^t g(x(\tau)) d\tau$ worden vergeleken, om de discrepantietermen te isoleren.
2. Universele Functies: De analyse levert twee universele functies van de orde $\alpha$ op, aangeduid als $c(\alpha)$ en $d(\alpha)$, die de afwijking van de invariantie kwantificeren.
3. Numerieke Verificatie: De theoretische voorspellingen worden getest op vier verschillende systemen:
   • Een nietlineaire oscillator (om coëfficiënten $c(\alpha)$ en $d(\alpha)$ te verifiëren).
   • De fractale van der Pol-oscillator (om de deformatie van de limietcyclus te onderzoeken).
   • Het fractale Rössler-systeem (om periodedubbelingscascades en bifurcatiedrempels te onderzoeken).
   • Het fractale Chen-systeem (om de stabiliteit van chaotische regimes en periodieke vensters te testen).
   • Twee numerieke schema's werden gebruikt: de fractale Adams-Bashforth-Moulton predictor-corrector methode en een diffusieve representatiemethode (superpositie van exponentiële modi) om robuustheid te garanderen.

Belangrijkste Bijdragen en Resultaten

1. Falen van Klok-covariantie: Het artikel bewijst dat voor $0 < \alpha < 1$, de banen van een fractaal systeem niet invariant zijn onder een toestandsafhankelijke klok.

   • Tijdwet Discrepantie ($t^{2\alpha}$): Op de eerste niet-verdwijnende orde is het verschil tussen de geklokte oplossing en de gereparametriseerde referentie-oplossing tangent aan de baan. Deze discrepantie wordt beheerst door de universele functie:
     $$c(\alpha) = \frac{1}{\Gamma(1+2\alpha)} - \frac{\alpha}{\Gamma(1+\alpha)\Gamma(2+\alpha)}$$
     Deze functie verdwijnt alleen bij $\alpha=1$ (waardoor de integer-orde stelling wordt hersteld) en nadert 1 als $\alpha \to 0$.
   • Geometrische Discrepantie ($t^{3\alpha}$): Op de volgende orde ontstaat een transversale component, wat bewijst dat de banen geometrisch verschillende curven zijn. Dit wordt beheerst door een tweede universele functie $d(\alpha)$, die voor alle $\alpha \in (0, 1)$ niet nul is.
   • Numerieke Bevestiging: De coëfficiënten $c(\alpha)$ en $d(\alpha)$ zijn met hoge precisie numeriek geverifieerd (fout < 0,1% voor $c$ en < 1% voor $d$) tegen de theoretische voorspellingen.

2. Invariantie van Lokale Stabiliteit: De auteur toont aan dat lokale eigenschappen invariant blijven. De evenwichtspunten, hun lineaire stabiliteit en de Matignon-kritische orde (de drempel $\alpha^*$ waarbij stabiliteit verloren gaat, bepaald door het argument van de eigenwaarden) blijven onveranderd door de klok. Dit komt omdat het Matignon-criterium afhangt van de argumenten van de eigenwaarden, die niet worden beïnvloed door vermenigvuldiging met een positieve scalaire $g(x^*)$.

3. Niet-invariantie van Globale Dynamica: Hoewel lokale stabiliteit behouden blijft, worden globale structuren vervormd.

   • Limietcycli: In de fractale van der Pol-oscillator verandert het wijzigen van de klok de vorm en amplitude van de limietcyclus met tot wel 16% voor intermediaire orden, een fenomeen dat in integer-orde onmogelijk is.
   • Bifurcaties: In het fractale Rössler-systeem verschuift de klok de waarden van lokale maxima, maar verandert het niet de sequentie van periodedubbelingsbifurcaties of de kritische orde voor het ontstaan van chaos.
   • Chaos versus Periodiciteit: In het fractale Chen-systeem kan een toestandsafhankelijke Arrhenius-klok (afhankelijk van een temperatuur-achtige variabele $z$) een periodiek regime transformeren in een chaotisch regime (of vice versa) bij dezelfde orde $\alpha$ en vanuit dezelfde beginconditie. Specifiek, bij $\alpha=0,94$ is het systeem periodiek met een triviale klok en chaotisch met een specifieke Arrhenius-klok.

4. Algemene Regel: Het artikel stelt een regel vast voor de verschuiving van bifurcaties:

   • Lokale bifurcaties (evenwichtspuntsstabiliteit, Hopf-bifurcaties) bewegen niet.
   • Globale bifurcaties (periodedubbelingen, crisis, geboorte/dood van periodieke vensters) verschuiven lineair met de sterkte van de klok.
   • De "fragiliteit" van een systeem voor een klokverandering wordt bepaald door de nabijheid van een globale bifurcatie.

Betekenis en Claims
Het artikel claimt dat voor fractale systemen de "klok" niet louter een wiskundig gemak is, maar een fysiek onderdeel van het model.

• Fysieke Interpretatie: In integer-orde modellen is het vermenigvuldigen van een veld met een toestandsafhankelijke factor (bijv. een temperatuurafhankelijke snelheid) gelijk aan het herdefiniëren van de tijdseenheid. In fractale modellen is deze keuze een andere fysieke hypothese:

  • $D^\alpha x = g(x)f(x)$ impliceert dat geheugen zich ophoopt in de fysieke tijd.
  • $D^\alpha X = f(X)$ in een intrinsieke tijd $s$ impliceert dat geheugen zich ophoopt in de intrinsieke/materiële tijd (een concept gebruikt in visco-elasticiteit en glasachtige dynamica).
    Het artikel stelt dat dit twee verschillende modellen zijn met verschillende voorspellingen, en dat het verschil wordt gekwantificeerd door de universele functies $c(\alpha)$ en $d(\alpha)$.

• Implicaties voor Chaostheorie: De auteur concludeert dat "chaos niet invariant is" van fractale systemen. Minimale orden voor chaos, bifurcatiediagrammen en de vormen van vreemde attractoren zijn eigenschappen van het specifieke veld-klok-paar, en niet invarianten van het systeem zelf. Bijgevolg is de noodzakelijke voorwaarde voor chaos op basis van Matignon's criterium robuust, maar het daadwerkelijke optreden van chaos hangt af van de specifieke formulering van de kinetische factoren.

• Praktisch Gevolg: Onderzoekers die fractale modellen gebruiken, moeten expliciet verklaren of het geheugen zich ophoopt in de fysieke tijd of in een intrinsieke materiële tijd. Een fractale attractor getekend in de toestandsruimte is een eigenschap van het vectorveld en de klok die wordt gebruikt om het te draaien; een andere klok zou een andere attractor opleveren, mogelijk met een andere aard (periodiek versus chaotisch).