Towards mesoscopic modeling of firing neurons: a feasibility study
Independent thesis Advanced level (professional degree), 20 credits / 30 HE creditsStudent thesis
Ion channel models are related to non-equilibrium statistical physics, fluid mechanics and electromagnetism. Some classes of ordinary differential equations that model ion channels can be seen as a limit of finite state-space continuous-time Markov chains. The purpose of this thesis is to qualitatively investigate the numerical results of systems of equations that incorporate ion channels modeled by such Markov chains and an electrical circuit model of a single neuron with isopotential extracellular space. This may be useful for making more detailed micro-physical simulations of neurons.
A subset of the Rallpack benchmarks is conducted in order to evaluate the accuracy of the electrical circuit model of the transmembrane voltage propagation. In order to test the tau-leap method employed to simulate the Markov-chain based ion channel models a cylindrical geometry is implemented. Convergence properties are presented in terms of mean interspike intervals of the transmembrane voltages for different time- and spatial discretisations. Accuracy of the tau-leap method is presented in relation to the deterministic versions of the ion channel models.
The results show that the method used to simulate the transmembrane voltages is accurate and that while the tau-leap method is convergent in the mean interspike interval sense, it is not conclusive how accurate it is compared to the corresponding ordinary differential equations or how efficient it is.
Place, publisher, year, edition, pages
UPTEC F, ISSN 1401-5757 ; 14011
tau-leaping, ion channels
IdentifiersURN: urn:nbn:se:uu:diva-222876OAI: oai:DiVA.org:uu-222876DiVA: diva2:712486
Master Programme in Engineering Physics
Nyberg, TomasLötstedt, Per