Finite difference/spectral approximations for the fractional cable equation

Authors:
Yumin Lin, Xianjuan Li and Chuanju Xu

Journal:
Math. Comp. **80** (2011), 1369-1396

MSC (2010):
Primary 65M12, 65M06, 65M70, 35S10

DOI:
https://doi.org/10.1090/S0025-5718-2010-02438-X

Published electronically:
December 2, 2010

MathSciNet review:
2785462

Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: The Cable equation has been one of the most fundamental equations for modeling neuronal dynamics. In this paper, we consider the numerical solution of the fractional Cable equation, which is a generalization of the classical Cable equation by taking into account the anomalous diffusion in the movement of the ions in neuronal system. A schema combining a finite difference approach in the time direction and a spectral method in the space direction is proposed and analyzed. The main contribution of this work is threefold: 1) We construct a finite difference/Legendre spectral schema for discretization of the fractional Cable equation. 2) We give a detailed analysis of the proposed schema by providing some stability and error estimates. Based on this analysis, the convergence of the method is rigourously established. We prove that the overall schema is unconditionally stable, and the numerical solution converges to the exact one with order , where and are respectively the time step size, polynomial degree, and regularity in the space variable of the exact solution. and are two different exponents between 0 and 1 involved in the fractional derivatives. 3) Finally, some numerical experiments are carried out to support the theoretical claims.

**1.**F. Amblard, A. C. Maggs, B. Yurke, A. N. Pargellis, and S. Leibler.

Subdiffusion and anomalous local viscoelasticity in actin networks.*Phys. Rev. Lett.*, 77:4470, 1996.**2.**E. Barkai, R. Metzler, and J. Klafter.

From continuous time random walks to the fractional Fokker-Planck equation.*Phys. Rev. E*, 61:132-138, 2000. MR**1736459 (2000j:82035)****3.**C. Bernardi and Y. Maday.*Approximations spectrales de problemès aux limites elliptiques*.

Springer-Verlag, 1992. MR**1208043 (94f:65112)****4.**C. Bernardi and Y. Maday.

Spectral methods.

In P. G. Ciarlet and J. L. Lions, editors,*Handbook of Numerical Analysis*. Elsevier, Amsterdam, 1997. MR**1470226****5.**J. Bouchaud and A. Georges.

Anomalous diffusion in disordered media: Statistical mechanisms: models and physical applications.*Phys. Rep.*, 195:127-293, 1990. MR**1081295****6.**E. Brown, E. Wu, W. Zipfel, and W. Webb.

Measurement of molecular diffusion in solution by multiphoton fluorescence photobleaching recovery.*Biophys. J.*, 77:2837-2849, 1999.**7.**W. H. Deng.

Finite element method for the space and time fractional Fokker-Planck equation.*SIAM J. Numer. Anal.*, 47(1):204-226, 2008. MR**2452858 (2009m:65175)****8.**M. Dentz, A. Cortis, H. Scher, and B. Berkowitz.

Time behaviour of solute in heterogeneous media: Transition from anomalous to normal transport.*Adv. Water Resources*, 27:155-173, 2004.**9.**V. J. Ervin and J. P. Roop.

Variational formulation for the stationary fractional advection dispersion equation.*Numer. Methods Partial Differential Equations*, 22:558-576, 2005. MR**2212226 (2006m:65265)****10.**T. Feder, I. Brust-Mascher, J. Slattery, B. Baird, and W. Webb.

Constrained diffusion or immobile fraction on cell surfaces: A new interpretation.*Biophys. J.*, 70:2767-2773, 1996.**11.**R. Ghosh.

Mobility and clustering of individual low density lipoprotein receptor molecures on the surface of human skin fibroblasts.*Ph.D. thesis*. Cornell University, Ithaca, NY.**12.**R. Ghosh and W. Webb.

Automated detection and tracking of individual and clustered cell surface low density lipoprotein receptor molecures.*Biophys. J.*, 66:1301-1318, 1994.**13.**I. Goychuk, E. Heinsalu, M. Patriarca, G. Schmid, and Pänggi.

Current and universal scaling in anomalous transport.*Phys. Rev. E*, 73:020101, 2006.**14.**B. I. Henry, T. A. M. Langlands, and S. L. Wearne.

Fractional cable models for spiny neuronal dendrites.*Phys. Rev. Lett.*, 100(12):128103, 2008.**15.**A. Hodgkin and A. Huxley.

A quantitative description of membrane current and its application to conduction and excitation in nerve.*J. Physiol.*, 117:500-544, 1952.**16.**J. Jack, D. Noble, and R. Tsien.*Electrical current flow in excitable cells*.

Oxford University Press, Oxford, 1975.**17.**D. Junge.*Nerve and Muscle Excitation (2nd edition ed.)*.

Sinauer Associates, Inc., Sunderland, Massachusetts, 1981.**18.**C. Koch.*Biophysics of Computation, Information Processing in Single neurons, Computational Neuroscience*.

Oxford University Press, New York, 1999.**19.**A. Kusumi, C. Nakada, K. Ritchie, K. Murase, K. Suzuki, H. Murakoshi, R. Kasai, J. Kondo, and T. Fujiwara.

Paradigm shift of the plasma membrane concept from two-dimensional continuum fluid to the partitioned fluid: Highspeed single-molecure tracking of membrane molecures.*Annu. Rev. Biophys. Biomol. Struct.*, 34:351-378, 2005.**20.**T. A. M. Langlands and B. I. Henry.

The accuracy and stability of an implicit solution method for the fractional diffusion equation.*J. Comput. Phys.*, 205(2):719-736, 2005. MR**2135000 (2005m:65174)****21.**E. K. Lenzi, R. S. Mendes, K. S. Fa, and L. C. Malacarne.

Anomalous diffusion: Fractional Fokker-Planck equation and its solutions.*J. Math. Phys.*, 44:2179-2185, 2003. MR**1972771 (2004f:82060)****22.**Xianjuan Li and Chuanju Xu.

A space-time spectral method for the time fractional diffusion equation.*SIAM J. Numerical Analysis*, 47(3):2108-2131, 2009. MR**2519596****23.**Y. M. Lin and C. J. Xu.

Finite difference/spectral approximation for the time fractional diffusion equations.*J. Comput. Phys.*, 2(3):1533-1552, 2007. MR**2349193 (2008i:65200)****24.**F. Mainardi.

Fractional diffusive waves in viscoelastic solids.*Nonlinear Waves in Solids*, pages 93-97, 1995.**25.**R. Metzler and J. Klafter.

The random walk's guide to anomalous diffusion: A fractional dynamics approach.*Phys. Rep.*, 339:1-77, 2000. MR**1809268 (2001k:82082)****26.**R. Metzler, J. Klafter, and I. Sokolov.

Anomalous transport in external fields: Continuous time random walks and fractional diffusion equations extended.*Phys. Rev. E*, 58:1621-1633, 1998.**27.**H. P. Müller, R. Kimmich, and J. Weis.

NMR flow velocity mapping in random percolation model objects: Evidence for a power-law dependence of the volume-averaged velocity on the probe-volume radius.*Phys. Rev. E*, 54:5278-5285, 1996.**28.**I. Podlubny.*Fractional Differential Equations*.

Academic Press, San Diego,, 1999. MR**1658022 (99m:26009)****29.**N. Qian and T. Sejnowski.

An electro-diffusion model for computing membrane potentials and ionic concentrations in branching dendrites, spines and axons.*Biol. Cybern.*, 62:1-15, 1989.**30.**W. Rall.

Branching dendritic trees and motoneuron membrane resistivity.

volume 1, pages 491-527. 1959.**31.**W. Rall.

Core conductor theory and cable properties of neurons.

In R. Poeter, editor,*Handbook of Physiology: The Nervous System, Vol. 1 (Chapter 3)*, pages 39-97. American Physiological Society, Bethesda, MD, 1977.**32.**K. Ritchie.

Detection of non-Browian diffusion in the cell membrane in single molecure tracking.*Biophys. J.*, 88:2266-2277, 2005.**33.**H. Scher and M. Lax.

Stochastic transport in a disordered solid.*Phys. Rev. B*, 7:4491-4502, 1973. MR**0391854 (52:12673)****34.**H. Scher and E. Montroll.

Anomalous transit-time dispersion in amorphous solids.*Phys. Rev. B*, 12:2455-2477, 1975.**35.**I. Segev, J. Fleshman, and R. Burke.

Compartmental models of complex neurons.

In C. Koch and I. Segev, editors,*Methods in Neuronal Modelling*. MIT Press, Cambridge, MA.**36.**R. Simson, B. Yang, S. Moore, P. Doherty, F. Walsh, and K. Jacobson.

Structural mosaicism on the submicron scale in the plasma membrane.*Biophys. J.*, 74:297-308, 1998.**37.**P. Smith, I. Morrison, K. Wilson, N. Fernandez, and R. Cherry.

Constrained diffusion or immobile fraction on cell surfaces: A new interpretation.*Biophys. J.*, 76:3331-3344, 1999.**38.**M. Wachsmuth, T. Weidemann, G. Muller, U. Hoffmann-Rohrer, T.Knoch, W. Waldeck, and J. Langowski.

Analyzing intracellular binding and diffusion with continuous fluorescence photobleaching.*Biophys. J.*, 84:3353-3363, 2003.**39.**E. R. Weeks.

Experimental studies of anomalous diffusion, blocking phenomena, and two-dimensional turbulence.*Ph.D. thesis*. University of Texas at Austin.**40.**W. Wyss.

The fractional diffusion equation.*J. Math. Phys.*, 27(11):2782-2785, 1986. MR**861345 (87m:44008)**

Retrieve articles in *Mathematics of Computation*
with MSC (2010):
65M12,
65M06,
65M70,
35S10

Retrieve articles in all journals with MSC (2010): 65M12, 65M06, 65M70, 35S10

Additional Information

**Yumin Lin**

Affiliation:
School of Mathematical Sciences, Xiamen University, 361005 Xiamen, China

**Xianjuan Li**

Affiliation:
School of Mathematical Sciences, Xiamen University, 361005 Xiamen, China

**Chuanju Xu**

Affiliation:
School of Mathematical Sciences, Xiamen University, 361005 Xiamen, China

Email:
cjxu@xmu.edu.cn

DOI:
https://doi.org/10.1090/S0025-5718-2010-02438-X

Keywords:
Fractional cable equation,
numerical solution,
stability,
convergence.

Received by editor(s):
February 7, 2009

Received by editor(s) in revised form:
April 11, 2010

Published electronically:
December 2, 2010

Additional Notes:
The research of the first author was partially supported by Fujian NSF under Grant S0750017.

The research of the third author was partially supported by National NSF of China under Grant 10531080, the Excellent Young Teachers Program by the Ministry of Education of China, and 973 High Performance Scientific Computation Research Program 2005CB321703

Article copyright:
© Copyright 2010
American Mathematical Society

The copyright for this article reverts to public domain 28 years after publication.