I ett quasi-integrabelt Hamiltonskt system utan störande krafter finns bevarade integraler som varaktigt bestämmer systemets rörelse. När hysteretiska krafter, viskoelastiska krafter, dämpning med fraktionell derivata eller tidsfördröjda krafter introduceras, bryts dessa idealiska bevarandeegenskaper delvis. Systemet är inte längre fullständigt integrabelt, men kan ofta behandlas som “nära” integrabelt med små störningar – och då blir metoder för stokastisk medelvärdesbildning centrala.

Hysteretiska krafter: dessa är krafter vars respons inte bara beror på den aktuella läget eller fördröjningen utan också på hur systemet tidigare belastats. De kan ge upphov till cykler där banan i fasutrymme uppvisar minne och “loopar”. Jämlikgöring av hysteretiska krafter innebär att man söker en ekvivalent kraftmodell som på medellång tid reproducerar samma energi- eller effektöverföring som den hysteretiska komponenten – detta möjliggör en enklare form för medelvärdesanalys.

Viscoelastiska krafter kombinerar elasticitet och viskositet: svaren är beroende inte bara på momentan töjning utan på dess tidsderivata. Dessa krafter kan ofta representeras med standardmodeller såsom Kelvin–Voigt eller Maxwell, och i quasi-integrabla system behandlas de som små störningar på de integrabla Hamiltonska delarna.

Dämpning med fraktionell derivata: när dämpningen inte är proportionell mot första tidsderivatan utan snarare involverar derivator av icke-heltalsordning, får systemet minneseffekter som decays långsammare och ibland icke-lokala i tiden. Detta inför komplexitet vad gäller stabilitet och medelvärdesvärden, men kan ofta approximativt införlivas i en stokastisk framework för att bearbeta långtidseffekter.

Tidsfördröjda krafter: om kraften beror på tidigare lägen eller tidigare hastigheter uppstår fördröjda återkopplingar. Dessa kan leda till resonansfenomen eller instabiliteter om fördröjningen är betydande i förhållande till tidsskalan hos systemets naturliga vibrationer.

Stokastisk medelvärdesbildning: ett huvudverktyg för att analysera hur små slumpmässiga påfrestningar eller brus (white noise, färgat brus, etc.) påverkar sådana quasi-integrabla system. Man “förenklar” genom att ersätta de snabba variationerna med deras medelvärde över tid eller ensemble, vilket leder till effektiv dynamik för de långsamma variablerna – ofta amplitud och fas. Systemets resonans eller icke-resonansförhållanden, interna resonanser, påverkar starkt hur medelvärdesprocessen gestaltar sig.

Det är också väsentligt att förstå hur dessa modeller används i naturvetenskapliga och tekniska sammanhang: predator‑bytesmodeller, ekologiska system med konkurrens och mättnad, rörelse hos aktiva Brown‑partiklar, DNA‑molekyler, och strukturer utsatta för vibrationer.

Viktigt: systemets stabilitetkonfigurationer (jämvikt, periodicitet, kaos) kan förändras radikalt när någon av dessa störningskrafter tillåts vara tillräckligt stor. Små parametrar kan leda till fenomen som bifurkationer eller språng i stabilitet. Beräkning av stationära sannolikhetsfördelningar, tid till första passage, effektiva reaktionshastigheter (reaction rate) kräver ofta analysera både deterministiska och stokastiska delar och hur de samverkar.

Hur intern resonans och stokastiska metoder påverkar Hamiltoniansystem

I moderna tillämpningar av stokastiska metoder inom fysik och matematik används den så kallade Itô-stokastiska differentialekvationen för att modellera system som är både stokastiska och Hamiltoniansystem. Dessa metoder används för att förstå dynamiken i system med externa störningar eller med fluktuationer som inte kan beskrivas genom deterministiska ekvationer. När det gäller generella Hamiltoniansystem, där interna resonanser kan förekomma, spelar dessa metoder en central roll i att förenkla och approximera komplexa system.

För att analysera dessa system införs ett begrepp om den genomsnittliga lösningen. Ett vanligt sätt att approximera stokastiska system är att ersätta tidsmedelvärde med ett rumsmedelvärde i relation till de parametrar som styr systemet, som till exempel q1, q2, och så vidare. Detta kan leda till en förenklad version av den ursprungliga FPK-ekvationen som inte bara förutser systemets fördelning utan också ger en bild av hur systemet rör sig i en typ av genomsnittlig potential.

När systemet är resonant, vilket innebär att vissa variabler är kopplade i en periodisk eller nära periodisk relation, kan den ursprungliga Hamiltoniansystemet betraktas som "komplett integrerbart". I dessa fall kan de dynamiska ekvationerna förenklas ytterligare genom att reducera antalet variabler och därmed få en mer hanterbar modell. Denna process innebär att man inför en uppsättning angulära variabler som beskriver systemets rörelse på en torus, vilket är ett vanligt sätt att behandla integrerbara system. Dessa variabler utnyttjas sedan för att skriva om de stokastiska differentialekvationerna, vilket ger en kraftig reduktion i komplexiteten av de ursprungliga ekvationerna.

En sådan approximation innebär att systemet kan beskrivas som en Markov-diffusionsprocess, vilket gör det möjligt att använda teorier som Khasminskiis teorem. Enligt detta teorem konvergerar systemet till ett n+α+M-dimensionellt system, där α representerar antalet interna resonanser och M antalet långsamt varierande parametrar i systemet. Denna konvergens innebär att vi kan ersätta tidsmedelvärden med fasmedelvärden och på så sätt förenkla lösningen av systemet. Genom att använda dessa metoder, tillsammans med stokastiska differentialekvationer, kan vi lösa systemet för att få en stationär lösning som beskriver sannolikhetsfördelningen för de olika variablerna.

Vid lösning av dessa ekvationer kan vi dessutom inkludera både initialvillkor och randvillkor, som till exempel periodiska randvillkor, vilket gör det möjligt att få en fullständig bild av systemets långsiktiga beteende. När den stationära lösningen väl har erhållits, kan den använda för att approximera sannolikhetsfördelningen för systemets ursprungliga parametrar.

För system där de generella Hamiltoniansystemen är separerbara, vilket innebär att de kan delas upp i oberoende delar, kan en ytterligare förenkling göras. Om subsystemen är både fullständigt integrerbara och resonanta, kan vi använda den separerbara formen för att förenkla både drift- och diffusionskoefficienterna i de stokastiska differentialekvationerna. Detta leder till en ännu mer precis approximation av systemets beteende, där effekterna av interna resonanser och externa störningar kan hanteras effektivt.

För att förstå hela komplexiteten i sådana system är det viktigt att inte bara fokusera på de initiala och stationära lösningarna utan också att uppmärksamma hur dessa lösningar interagerar med varandra över tid. Medan den stationära lösningen ger en bra approximation av systemets fördelning vid ett givet tillfälle, är det dynamiska beteendet hos systemet också av stor betydelse för att fullt förstå hur dessa stokastiska processer utvecklas.

När man arbetar med sådana system bör läsaren även vara medveten om att även små störningar kan ha stora effekter på systemets långsiktiga beteende. Dessa system kan uppvisa olika typer av fenomen som inte kan förutses genom enkla analytiska lösningar, vilket gör den stokastiska metodiken ovärderlig för en mer exakt förståelse av deras dynamik.

Hur stochastiska genomsnittsmetoden tillämpas på aktiva Browniska partiklar

De aktiva Browniska partiklarna som rör sig genom ett system beskrivet av Hamiltonska ekvationer kan analyseras med hjälp av stochastiska genomsnittsmetoder, vilket gör det möjligt att hantera komplexa resonansförhållanden och icke-linjära dynamiska system. En sådan metod tillämpas på nästan-integrerbara och resonanta Hamiltonsystem, som beskrivs i sektion 5.2.2. Här analyserar vi rörelsen hos dessa partiklar, där både interna resonansförhållanden och de stokastiska differentialekvationerna spelar en central roll.

För att formulera det matematiskt kan skillnaden i fasvinkel representeras som ψ=θ2θ1\psi = \theta_2 - \theta_1, där θ1\theta_1 och θ2\theta_2 är fasvinklarna för partiklarna. Stokastiska differentialekvationer som beskriver dessa system kan härledas genom att använda genomsnittsmetoden, vilket resulterar i en komplex ekvation där dämpning, amplituder och olika parametrar som γ1\gamma_1, γ2\gamma_2, och DD samverkar.

När man tillämpar metoden på systemet där de naturliga frekvenserna är lika, ω\omega, uppstår ett förhållande av intern resonans enligt ekvation (5.75) i volym 1. Detta gör att systemets rörelse kan beskrivas genom stokastiska genomsnitt och att partikelrörelserna kan analyseras som en Markov-process med diffusionsbeteende. För att lösa denna process används den förenklade Fokker-Planck-ekvationen (FPK), som modellerar systemets stationära sannolikhetsfördelning.

De moment som härleds ur dessa ekvationer, som till exempel a1a_1, a2a_2 och aψa_\psi, beskriver förändringarna i partiklarna över tid och är relaterade till dämpningseffekterna samt de stokastiska processerna i systemet. Genom att använda dessa moment kan man härleda lösningar på systemets stationära sannolikhetsdichtefunktion (PDF), vilket gör det möjligt att beräkna sannolikheter för partiklarna att befinna sig i olika positioner eller hastigheter.

Det är viktigt att förstå att den stokastiska genomsnittsmetoden inte bara tillämpas på system med Rayleigh-dämpning, utan även på andra typer av dämpningsmekanismer. Genom att byta ut Rayleigh-dämpning mot en annan form, som till exempel Erdmann-dämpning, kan man få en annan energihanteringsmekanism för systemet. Detta gör det möjligt att anpassa metoden till ett bredare spektrum av fysiska system, inklusive biologiska system som aktiva Browniska partiklar i cellrörelse.

För att lösa den reducerade FPK-ekvationen, där sannolikheten för systemets stationära tillstånd beskrivs av p(h1,h2,ψ)p(h_1, h_2, \psi), kan den analytiska lösningen härledas. Lösningen resulterar i en formel som relaterar sannolikhetsdichtefunktionen till de parametrar som definierar systemet, såsom dämpningskoefficienterna γ1\gamma_1, γ2\gamma_2 och DD. När lösningen väl är härledd kan man använda den för att förstå dynamiken hos aktiva partiklar, som till exempel celler i biologiska system.

För att verkligen förstå de fysiska implikationerna av dessa modeller är det avgörande att relatera teoretiska resultat med experimentella data. Genom att jämföra simuleringar och experimentella resultat för aktiva Browniska partiklar, som exempelvis granulocytrörelse i biologiska system, kan man validera och förbättra modellen. Det är här de statistiska metoderna och den stochastiska genomsnittsmetoden verkligen visar sin nytta, eftersom de gör det möjligt att extrahera viktiga fysikaliska parametrar från verkliga data.

I tillämpningar som involverar biologiska system eller andra komplexa dynamiska system, där många olika faktorer och parametrar spelar en roll, ger stochastiska metoder en robust ram för att analysera och förstå dessa system. Förutom att ge insikt i rörelse och interaktioner hos aktiva partiklar kan dessa metoder också användas för att modellera mer komplexa situationer, såsom cellrörelser i vävnader eller komplexa fysikaliska processer i materialvetenskap.

För att verkligen bemästra denna metod är det också viktigt att förstå de grundläggande antagandena om systemet, inklusive vad som sker vid gränserna för systemets stabilitet, och hur resonansförhållanden kan påverka systemets långsiktiga beteende. Stokastiska processer, särskilt i närvaro av resonans, kan leda till överraskande och icke-intuitiva beteenden, vilket gör att noggranna simuleringar och experimentella kontroller är nödvändiga för att säkerställa riktigheten hos de teoretiska modellerna.