TY - JOUR
T1 - Analytic Modeling of Neural Tissue
T2 - I. A Spherical Bidomain
AU - Schwartz, Benjamin L.
AU - Chauhan, Munish
AU - Sadleir, Rosalind
N1 - Funding Information:
This work was supported by a grant from the NIH (No. R01NS077004) to RJS.
Publisher Copyright:
© 2016, Schwartz et al.
PY - 2016/12/1
Y1 - 2016/12/1
N2 - Presented here is a model of neural tissue in a conductive medium stimulated by externally injected currents. The tissue is described as a conductively isotropic bidomain, i.e. comprised of intra and extracellular regions that occupy the same space, as well as the membrane that divides them, and the injection currents are described as a pair of source and sink points. The problem is solved in three spatial dimensions and defined in spherical coordinates (r, θ, ϕ). The system of coupled partial differential equations is solved by recasting the problem to be in terms of the membrane and a monodomain, interpreted as a weighted average of the intra and extracellular domains. The membrane and monodomain are defined by the scalar Helmholtz and Laplace equations, respectively, which are both separable in spherical coordinates. Product solutions are thus assumed and given through certain transcendental functions. From these electrical potentials, analytic expressions for current density are derived and from those fields the magnetic flux density is calculated. Numerical examples are considered wherein the interstitial conductivity is varied, as well as the limiting case of the problem simplifying to two dimensions due to azimuthal independence. Finally, future modeling work is discussed.
AB - Presented here is a model of neural tissue in a conductive medium stimulated by externally injected currents. The tissue is described as a conductively isotropic bidomain, i.e. comprised of intra and extracellular regions that occupy the same space, as well as the membrane that divides them, and the injection currents are described as a pair of source and sink points. The problem is solved in three spatial dimensions and defined in spherical coordinates (r, θ, ϕ). The system of coupled partial differential equations is solved by recasting the problem to be in terms of the membrane and a monodomain, interpreted as a weighted average of the intra and extracellular domains. The membrane and monodomain are defined by the scalar Helmholtz and Laplace equations, respectively, which are both separable in spherical coordinates. Product solutions are thus assumed and given through certain transcendental functions. From these electrical potentials, analytic expressions for current density are derived and from those fields the magnetic flux density is calculated. Numerical examples are considered wherein the interstitial conductivity is varied, as well as the limiting case of the problem simplifying to two dimensions due to azimuthal independence. Finally, future modeling work is discussed.
KW - Analytic modeling
KW - Bidomain
UR - http://www.scopus.com/inward/record.url?scp=84986630234&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84986630234&partnerID=8YFLogxK
U2 - 10.1186/s13408-016-0041-1
DO - 10.1186/s13408-016-0041-1
M3 - Article
AN - SCOPUS:84986630234
SN - 2190-8567
VL - 6
JO - Journal of Mathematical Neuroscience
JF - Journal of Mathematical Neuroscience
IS - 1
M1 - 9
ER -