Model Reduction of Large Spiking Neurons
| dc.contributor.author | Kellems, Anthony Richard | |
| dc.date.accessioned | 2018-06-19T17:46:05Z | |
| dc.date.available | 2018-06-19T17:46:05Z | |
| dc.date.issued | 2010-04 | |
| dc.date.note | April 2010 | |
| dc.description | This work was also published as a Rice University thesis/dissertation: http://hdl.handle.net/1911/61961 | |
| dc.description.abstract | This thesis introduces and applies model reduction techniques to problems associated with simulation of realistic single neurons. Neurons have complicated dendritic structures and spatially-distributed ionic kinetics that give rise to highly nonlinear dynamics. However, existing model reduction methods compromise the geometry, and thus sacrifice the original input-output relationship. I demonstrate that linear and nonlinear model reduction techniques yield systems that capture the salient dynamics of morphologically accurate neuronal models and preserve the input-output maps while using significantly fewer variables than the full systems. Two main dynamic regimes characterize the voltage response of a neuron, and I demonstrate that different model reduction techniques are well-suited to each regime. Small perturbations from the neuron's rest state fall into the subthreshold regime, which can be accurately described by a linear system. By applying Balanced Truncation (BT), a model reduction technique for general linear systems, I recover subthreshold voltage dynamics, and I provide an efficient Iterative Rational Krylov Algorithm (IRKA), which makes large problems of interest tractable. However, these approximations are not valid once the input to the neuron is sufficient to drive the voltage into the spiking regime, which is characterized by highly nonlinear behavior. To reproduce spiking dynamics, I use a proper orthogonal decomposition (POD) to reduce the number of state variables and a discrete empirical interpolation method (DEIM) to reduce the complexity of the nonlinear terms. The techniques described above are successful, but they inherently assume that the whole neuron is either passive (linear) or active (nonlinear). However, in realistic cells the voltage response at distal locations is nearly linear, while at proximal locations it is very nonlinear. With this observation, I fuse the aforementioned models together to create a reduced coupled model in which each reduction technique is used where it is most advantageous, thereby making it possible to more accurately simulate a larger class of cortical neurons. | |
| dc.format.extent | 160 pp | |
| dc.identifier.citation | Kellems, Anthony Richard. "Model Reduction of Large Spiking Neurons." (2010) <a href="https://hdl.handle.net/1911/102156">https://hdl.handle.net/1911/102156</a>. | |
| dc.identifier.digital | TR10-14 | |
| dc.identifier.uri | https://hdl.handle.net/1911/102156 | |
| dc.language.iso | eng | |
| dc.relation.HasVersion | Kellems, A.R., Chaturantabut, S., Sorensen, D.C., and Cox, S.J. (2010) Morphologically accurate reduced order modeling of spiking neurons. J. Comput. Neurosci. The original publication is available at www.springerlink.com: DOI: 10.1007/s10827-010-0229-4; Kellems, A.R., Roos, D., Xiao, N., and Cox, S.J. (2009) Low-dimensional, morphologically accurate models of subthreshold membrane potential. J. Comput. Neurosci., 27:161-176 The original publication is available at www.springerlink.com: DOI: 10.1007/s10827-008-0134-2 | |
| dc.title | Model Reduction of Large Spiking Neurons | |
| dc.type | Technical report | |
| dc.type.dcmi | Text |
Files
Original bundle
1 - 1 of 1