Discontinuous Galerkin finite element method for solving population density functions of cortical pyramidal and thalamic neuronal populations

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons.

Original languageEnglish
Pages (from-to)150-158
Number of pages9
JournalComputers in Biology and Medicine
Volume57
DOIs
Publication statusPublished - 2015 Feb 1

Fingerprint

Population Density
Probability density function
Neurons
Finite element method
Convergence of numerical methods
Population
Monte Carlo Method
Finite difference method
Partial differential equations
Pyramidal Cells
Monte Carlo methods
Electroencephalography
Brain
Population Dynamics

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Health Informatics

Cite this

@article{118dedcd059e4d9fa33828cddc5c5d61,
title = "Discontinuous Galerkin finite element method for solving population density functions of cortical pyramidal and thalamic neuronal populations",
abstract = "Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons.",
author = "Huang, {Chih Hsu} and Chou-Ching Lin and Ming-Shaung Ju",
year = "2015",
month = "2",
day = "1",
doi = "10.1016/j.compbiomed.2014.12.011",
language = "English",
volume = "57",
pages = "150--158",
journal = "Computers in Biology and Medicine",
issn = "0010-4825",
publisher = "Elsevier Limited",

}

TY - JOUR

T1 - Discontinuous Galerkin finite element method for solving population density functions of cortical pyramidal and thalamic neuronal populations

AU - Huang, Chih Hsu

AU - Lin, Chou-Ching

AU - Ju, Ming-Shaung

PY - 2015/2/1

Y1 - 2015/2/1

N2 - Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons.

AB - Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons.

UR - http://www.scopus.com/inward/record.url?scp=84920278255&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84920278255&partnerID=8YFLogxK

U2 - 10.1016/j.compbiomed.2014.12.011

DO - 10.1016/j.compbiomed.2014.12.011

M3 - Article

C2 - 25557200

AN - SCOPUS:84920278255

VL - 57

SP - 150

EP - 158

JO - Computers in Biology and Medicine

JF - Computers in Biology and Medicine

SN - 0010-4825

ER -