Mathematical and numerical modeling of coupled oxygen dynamics and neuronal electrophysiology

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Daniele, Francesco, Saglio, Caterina B. Leimer, Pagani, Stefano, Antonietti, Paola F.
Natura: Preprint
Pubblicazione: 2025
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866911526236979200
author Daniele, Francesco
Saglio, Caterina B. Leimer
Pagani, Stefano
Antonietti, Paola F.
author_facet Daniele, Francesco
Saglio, Caterina B. Leimer
Pagani, Stefano
Antonietti, Paola F.
contents Modeling and simulating how oxygen supply shapes neuronal excitability is crucial for advancing the understanding of brain function in pathological scenarios, such as ischemia. This condition is caused by a reduced blood supply, leading to the deprivation of oxygen and other metabolites; this energy deficit impairs ionic pumps and causes cellular swelling. In this work, this phenomenon is modeled through a volumetric variation law that links cell swelling to local oxygen concentration and the percentage of blood flow reduction. The swelling law links volume changes to local oxygen and the degree of blood-flow depression, providing a simple mechanistic pathway from hypoxia to tortuosity-driven transport impairment. The interplay between oxygen supply and excitability in brain tissue is described by coupling the monodomain model with specific neuronal ionic and metabolic models that characterize ion and metabolite concentration dynamics. The numerical approximation of this coupled multiscale problem is particularly challenging, owing to the presence of sharp and fast-propagating wavefronts and complex geometrical domains. To address these challenges, suitable space- and time-adaptive schemes are employed to capture the action potential dynamics accurately. This multiscale model is discretized in space with the high-order p-adaptive discontinuous Galerkin method on polygonal and polyhedral grids (PolyDG) and integrated in time with a Crank-Nicolson scheme. We numerically investigate different pathological scenarios on a two-dimensional idealized square domain and on a realistic geometry, both discretized with a polygonal grid, analyzing how subclinical and severe ischemia can affect brain electrophysiology and metabolic concentrations.
format Preprint
id arxiv_https___arxiv_org_abs_2509_22863
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Mathematical and numerical modeling of coupled oxygen dynamics and neuronal electrophysiology
Daniele, Francesco
Saglio, Caterina B. Leimer
Pagani, Stefano
Antonietti, Paola F.
Numerical Analysis
Modeling and simulating how oxygen supply shapes neuronal excitability is crucial for advancing the understanding of brain function in pathological scenarios, such as ischemia. This condition is caused by a reduced blood supply, leading to the deprivation of oxygen and other metabolites; this energy deficit impairs ionic pumps and causes cellular swelling. In this work, this phenomenon is modeled through a volumetric variation law that links cell swelling to local oxygen concentration and the percentage of blood flow reduction. The swelling law links volume changes to local oxygen and the degree of blood-flow depression, providing a simple mechanistic pathway from hypoxia to tortuosity-driven transport impairment. The interplay between oxygen supply and excitability in brain tissue is described by coupling the monodomain model with specific neuronal ionic and metabolic models that characterize ion and metabolite concentration dynamics. The numerical approximation of this coupled multiscale problem is particularly challenging, owing to the presence of sharp and fast-propagating wavefronts and complex geometrical domains. To address these challenges, suitable space- and time-adaptive schemes are employed to capture the action potential dynamics accurately. This multiscale model is discretized in space with the high-order p-adaptive discontinuous Galerkin method on polygonal and polyhedral grids (PolyDG) and integrated in time with a Crank-Nicolson scheme. We numerically investigate different pathological scenarios on a two-dimensional idealized square domain and on a realistic geometry, both discretized with a polygonal grid, analyzing how subclinical and severe ischemia can affect brain electrophysiology and metabolic concentrations.
title Mathematical and numerical modeling of coupled oxygen dynamics and neuronal electrophysiology
topic Numerical Analysis
url https://arxiv.org/abs/2509.22863