Splitting methods for stochastic Hodgkin-Huxley type systems and a localized fundamental mean-square convergence theorem
par
Pierre Etoré(LJK (Université Grenoble Alpes))
→
Europe/Paris
Amphi Schwartz
Amphi Schwartz
Description
Motivated by the study of splitting schemes for the Hodgkin-Huxley model, we introduce a novel localized version of the fundamental mean-square convergence theorem for numerical schemes for SDEs with locally Lipschitz coefficients. Indeed existing fundamental theorems for mean-square convergence for schemes for stochastic differential equations (SDEs) usually require globally or one-sided Lipschitz continuous coefficients. Those assumptions are not satisfied by the Hodgkin-Huxley model. Specifically, we show that if the coefficients of the SDE are locally Lipschitz and a numerical scheme is locally consistent in the mean-square sense of order q>1 then it is locally mean-square convergent with rate q-1. Building on this result,we further prove that global mean-square convergence follows, provided that both the exact solution and its numerical approximation admit bounded 2pth moments for some p>1.
Having at hand these new convergence results we perform fully the study of convergence of splitting schemes for the Hodgkin-Huxley model. This model describes the neural activity of a giant squid and is characterized by a conditionally linear drift structure. We perform innovative proofs of local consistency and boundedness of moments. In addition, we establish key structure-preserving properties of the splitting methods, in particular state-space preservation. Numerical experiments support the theoretical results and demonstrate that the proposed splitting schemes significantly outperform Euler-Maruyama type methods in preserving the qualitative features of the model.