A modified Euler-Maruyama method to simulate a one-dimensional sticky diffusion

Technische Samenvatting: Een gemodificeerde Euler-Maruyama-methode voor het simuleren van een eendimensionale plakkerige diffusie

Probleemstelling
Plakkerige diffusies (sticky diffusions) zijn stochastische processen die in staat zijn om een niet-nul tijdsduur door te brengen op een lager-dimensionale grens (bijvoorbeeld een punt in één dimensie). Hoewel ze sinds het werk van Feller in de jaren 1950 wiskundig goed bestudeerd zijn, vormt het simuleren van deze processen een unieke uitdaging: standaard numerieke methoden voor stochastische differentiaalvergelijkingen (SDE's), zoals het Euler-Maruyama (EM) schema, nemen willekeurige stappen in continue ruimte en landen daardoor bijna nooit exact op de grens. Bijgevolg slagen ze er niet in de singulariteit van de maat (de "plakkerigheid") die inherent is aan het proces te vangen. Bestaande alternatieven, zoals benaderingen met continue-tijd Markov-ketens op vooraf gedefinieerde roosters, kampen met schaalbaarheidsproblemen in hogere dimensies en moeilijkheden bij het afhandelen van gekromde grenzen. Boveraan zijn recente methode gebaseerd op tijdverandering (time-change) van gereflecteerde diffusies beperkt tot gevallen waarbij de drift- en diffusiecoëfficiënten identiek zijn in het interieur en op de grens.

Methodologie
De auteurs stellen een gemodificeerd Euler-Maruyama-algoritme voor dat is ontworpen om eendimensionale plakkerige diffusies te simuleren die exact naar en van een plakkerig punt (genomen als $x=0$) kunnen springen. De methode is een minimale modificatie van het standaard gereflecteerde EM-schema, waarbij een probabilistisch mechanisme wordt geïntroduceerd om de grenslaag te behandelen.

Het algoritme verdeelt het domein in twee regio's op basis van de tijdstap $h$:

1. Exterieur ($\Omega_e$): Wanneer het proces ver van de grens is, evolueert het met een standaard gereflecteerde Euler-Maruyama update.
2. Grenslaag ($\Omega_b$): Wanneer het proces dicht bij de grens is (specifiek binnen een afstand die schaalt met $\sqrt{h}$), voert het algoritme een gemengde update uit. Met een waarschijnlijkheid $\lambda(x)$ voert het een gereflecteerde EM-stap uit. Met een waarschijnlijkheid $1 - \lambda(x)$ springt het exact naar het plakkerige punt ($x=0$).

De kern van de methodologie ligt in de afleiding van de sprongkansfunctie $\lambda(x)$. De auteurs leiden $\lambda(x)$ af door te eisen dat de transitieoperator van het discrete schema de transitieoperator van de continue plakkerige diffusie over een enkele tijdstap $h$ benadert. Dit wordt bereikt door:

• De generator van de plakkerige diffusie en het discrete schema te expanderen met behulp van Taylorreeksen.
• De plakkerige randvoorwaarde (een relatie tussen de eerste en tweede ruimtelijke afgeleiden op de grens) te gebruiken om de eerste afgeleide uit te drukken in termen van de tweede.
• $\lambda(x)$ op te lossen zodat de coëfficiënten van de drift- en diffusietermen in de lokale foutexpansie overeenkomen met de continue dynamica tot een gespecificeerde tolerantie.

Voor een Sticky Brownian Motion (SBM) met drift $b=0$ en diffusie $\sigma=1$, is de afgeleide kans:
$$\lambda(x) = \frac{x^2 + 2\kappa x + h(\frac{\kappa}{\sqrt{3h}} + 1)}{x^2 + h + \kappa\sqrt{3h}}$$
waarbij $\kappa$ de plakkerigheidsparameter is. De methode wordt vervolgens gegeneraliseerd naar willekeurige eendimensionale plakkerige diffusies met drift $b(x)$ en diffusie $\sigma(x)$, wat resulteert in een complexere uitdrukking voor $\lambda(x)$ die rekening houdt met de toestand-afhankelijke coëfficiënten.

Belangrijkste Bijdragen

1. Algoritmeontwerp: Het artikel introduceert een eenvoudig, computationeel efficiënt algoritme dat een discreet proces in staat stelt exact op de grens te verblijven, waarmee de beperking van standaard EM-methoden wordt overwonnen.
2. Waarschijnlijkheidsafleiding: De auteurs bieden een rigoureuze afleiding voor de sprongkansfunctie $\lambda(x)$, waardoor wordt gewaarborgd dat het discrete schema consistent is met de generator van de gewenste plakkerige diffusiedynamiek.
3. Convergentieanalyse: Het artikel bewijst dat het voorgestelde schema zwak convergeert met orde 1 (globale fout $O(h)$).
   • Het bewijs ontleedt de globale fout in lokale fouten en het aantal stappen doorgebracht in elke regio.
   • In het exterieur is de lokale fout $O(h^2)$.
   • In de grenslaag is de lokale fout $O(h^{3/2})$.
   • Cruciaal is dat de auteurs bewijzen dat het aantal stappen doorgebracht in de grenslaag schaalt als $O(h^{-1/2})$. Het product van de lokale fout en het aantal stappen levert de globale $O(h)$ convergentie op.
4. Regulariteitsvoorwaarden: Het artikel stelt voldoende voorwaarden vast aan de terminale functie $\phi$ (compatibiliteitsvoorwaarden bij de grenzen) om de regulariteit van de oplossing van de backward equation te waarborgen die vereist is voor de convergentiebewijzen.

Resultaten

• Theoretisch: De auteurs bewijzen de zwakke convergentie van orde 1 voor zowel Sticky Brownian Motion als algemene plakkerige diffusies met toestand-afhankelijke drift en diffusie. De bewijzen steunen op het begrenzen van de lokale fout en het aantonen dat de bezettingstijd van de grenslaag gecontroleerd is.
• Numeriek: Numerieke experimenten valideren deze theoretische bevindingen. Simulaties van SBM en algemene plakkerige diffusies met verschillende drift- en diffusiecoëfficiënten laten empirische convergentiesnelheden zien die dicht bij 1 liggen. De resultaten tonen aan dat de methode effectief de tijd die nabij de plakkerige grens wordt doorgebracht vastlegt, waarbij trajecten langere verblijftijden vertonen naarmate de plakkerigheidsparameter $\kappa$ toeneemt.

Betekenis en Claims
De auteurs positioneren dit werk als een belangrijke stap richting een algemene methode voor het simuleren van willekeurige plakkerige diffusies. Zij stellen dat hun methode voordelig is omdat:

• Het een kleine modificatie is van standaard simulatiemethoden, waardoor het gemakkelijker te integreren is in bestaande kaders (bijv. die die interactieve deeltjessystemen of chemische reacties afhandelen).
• Het niet vereist dat het proces een tijdverandering (time-change) is van een gereflecteerde diffusie, wat het mogelijk maakt om verschillende drift- en diffusiecoëfficiënten op de grens te gebruiken dan in het interieur.
• Het de computationele knelpunten vermijdt van roostergebaseerde Markov-ketenbenaderingen in hogere dimensies.

Het artikel concludeert door op te merken dat hoewel het huidige werk zich op één dimensie richt, de beschreven heuristieken en convergentie-instrumenten een pad bieden voor het uitbreiden van deze ideeën naar meerdimensionale plakkerige diffusies, hoewel het waarborgen van de geldigheid van de schakelwaarschijnlijkheden in hogere dimensies een technische uitdaging blijft. De auteurs suggereren ook dat toekomstig werk hogere-orde schema's en toepassingen op manifolds zou kunnen verkennen.