A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

Journal of Computational Physics

Volumn 230, Issue 4, 2011, Pages 1084-1099

  Cáceres, María J ^a   Carrillo, José A ^b   Tao, Louis ^c  

^a UNIVERSITY OF GRANADA   (Spain)

^b UNIVERSITAT AUTÒNOMA DE BARCELONA   (Spain)

^c PEKING UNIVERSITY   (China)

Author keywords

Chang Cooper method; Deterministic simulations; Direct simulation Monte Carlo; Fokker Planck equation; Kinetic equations; Neuronal network; WENO methods

Indexed keywords

COMPUTATION THEORY; DISTRIBUTION FUNCTIONS; INTELLIGENT SYSTEMS; KINETIC ENERGY; KINETIC THEORY; KINETICS; MONTE CARLO METHODS; NEURAL NETWORKS; NEURONS; NONLINEAR EQUATIONS; NUMERICAL METHODS; PROBABILITY DENSITY FUNCTION;

CHANG-COOPER METHOD; DETERMINISTIC SIMULATION; DIRECT SIMULATION MONTE CARLO; FOKKER-PLANCK-EQUATION; KINETIC EQUATIONS; NEURONAL NETWORKS; NONLINEAR FOKKER-PLANCK EQUATIONS; NUMERICAL SOLVERS; STATIONARY STATE; WENO METHOD;

FOKKER PLANCK EQUATION;

EID: 78650564252     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2010.10.027     Document Type: Article

References (58)

1. Abbott L., van Vreeswijk C. Asynchronous states in networks of pulse-coupled oscillators. Phys. Rev. E 1993, 48:1483-1490.
2. Acebrón J.A., Bulsara A.R., Rappel W.-J. The noisy Fitzhugh-Nagumo model: from single elements to globally coupled networks. Phys. Rev. E 2004, 69:026202.
3. Albantakis L., Deco G. The encoding of alternatives in multiple-choice decision making. Proc. Natl. Acad. Sci. USA 2009, 106(25):10308-10313.
4. Anderson J., Lampl I., Gillespie D., Ferster D. The contribution of noise to contrast invariance of orientation tuning in cat visual cortex. Science 2000, 290:1968-1972.
5. Apfaltrer F., Ly C., Tranchina D. Population density methods for stochastic neurons with realistic synaptic kinetics: Firing rate dynamics and fast computational methods. Network: Comp. Neural Syst. 2006, 17:373-418.
6. Barna G., Grobler T., Erdi P. Statistical model of the hippocampal ca3 region, ii. the population framework: model of rhythmic activity in ca3 slice. Biol. Cybern. 1998, 79:309-321.
7. Bonilla L.L., Carrillo J.A., Soler J. Asymptotic behavior of an initial-boundary value problem for the Vlasov-Poisson-Fokker-Planck system. SIAM J. Appl. Math. 1997, 57(5):1343-1372.
8. Brunel N. Dynamics of sparsely connected networks of excitatory and inhibitory spiking networks. J. Comp. Neurosci. 2000, 8:183-208.
9. Brunel N., Hakim V. Fast global oscillations in networks of integrate-and-fire neurons with long fitting rates. Neural Comput. 1999, 11:1621-1671.
10. Brunel N., Wang X.-J. What determines the frequency of fast network oscillations with irregular neural discharge? I. synaptic dynamics and excitation-inhibition balance. J. Neurophysiol. 2003, 90:415-430.
11. Buet C., Cordier S., Dos Santos V. A conservative and entropy scheme for a simplified model of granular media. Transport Theory Statist. Phys. 2004, 33(2):125-155.
12. Cáceres M.J., Carrillo J.A., Gamba I.M., Majorana A., Shu C.-W. Deterministic kinetic solvers for charged particle transport in semiconductor devices. Model. Simul. Sci. Eng. Technol. 2007, 151-171.
13. Cai D., Tao L., McLaughlin D.W. An embedded network approach for scale-up of fluctuation-driven systems with preservation of spike information. PNAS 2004, 101:14288-14293.
14. Cai D., Tao L., Rangan A.V., McLaughlin D.W. Kinetic theory for neuronal network dynamics. Commun. Math. Sci. 2006, 4(1):97-127.
15. Cai D., Tao L., Shelley M.J., McLaughlin D.W. An effective kinetic representation of fluctuation-driven neuronal networks with application to simple and complex cells in visual cortex. PNAS 2004, 101:7757-7762.
16. Carrillo J.A., Gamba I.M., Majorana A., Shu C.-W. A WENO-solver for the transients of Boltzmann-Poisson system for semiconductor devices: performance and comparisons with Monte Carlo methods. J. Comput. Phys. 2003, 184(2):498-525.
17. Carrillo J.A., Gamba I.M., Majorana A., Shu C.-W. 2D semiconductor device simulations by WENO-Boltzmann schemes: efficiency, boundary conditions and comparison to Monte Carlo methods. J. Comput. Phys. 2006, 214(1):55-80.
18. Carrillo J.A., Majorana A., Vecil F. A semi-lagrangian deterministic solver for the semiconductor Boltzmann-Poisson system. Commun. Comput. Phys. 2007, 2:1027-1054.
19. Carrillo J.A., Vecil F. Nonoscillatory interpolation methods applied to Vlasov-based models. SIAM J. Sci. Comput. 2007, 29(3):1179-1206. electronic.
20. Chawanya T., Aoyagi A., Nishikawa T., Okuda K., Kuramoto Y. A model for feature linking via collective oscillations in the primary visual cortex. Biol. Cybern. 1993, 68:483-490.
21. Fourcaud N., Brunel N. Dynamics of the firing probability of noisy integrate-and-fire neurons. Neural Comp. 2002, 14:2057-2110.
22. Gerstner W. Population dynamics of spiking neurons: fast transients, asynchronous states, and locking. Neural Comp. 2000, 12:43-89.
23. Gray C.M., Singer W. Stimulus-specific neuronal oscillations in orientation columns of cat visual cortex. Proc. Natl. Acad. Sci. USA 1989, 86(5):1698-1702.
24. Haskell E., Nykamp D., Tranchina D. Population density methods for large-scale modeling of neuronal networks with realistic synaptic kinetics: cutting the dimension down to size. Network: Compt. Neural Syst. 2001, 12:141-174.
25. Henrie J.A., Shapley R. Lfp power spectra in v1 cortex: the graded effect of stimulus contrast. J. Neurophysiol. 2005, 94(1):479-490.
26. Izhikevich E.M., Edelman G.M. Large-scale model of mammalian thalamocortical systems. Proc. Natl. Acad. Sci. USA 2008, 105(9):593-3598.
27. Jiang G.-S., Shu C.-W. Efficient implementation of weighted ENO schemes. J. Comput. Phys. 1996, 126(1):202-228.
28. Knight B. Dynamics of encoding in a populaton neurons. J. Gen. Physiol. 1972, 59:734-766.
29. Koulakov A., Raghavachari S., Kepecs A., Lisman J. Model for a robust neural integrator. Nat. Neurosci. 2002, 5:775-782.
30. Ly C., Tranchina D. Critical analysis of dimension reduction by a moment closure method in a population density approach to neural network modeling. Neural Comput. 2007, 19(8):2032-2092.
31. Marino J., Schummers J., Lyon D.C., Schwabe L., Beck O., Wiesing P., Obermayer K., Sur M. Invariant computations in local cortical networks with balanced excitation and inhibition. Nature Neurosci. 2005, 8:194-201.
32. Mattia M., Del Giudice P. Population dynamics of interacting spiking neurons. Phys. Rev. E 2002, 66:051917.
33. McLaughlin D., Shapley R., Shelley M., Wielaard J. A neuronal network model of macaque primary visual cortex (V1): orientation selectivity and dynamics in the input layer 4Cα. Proc. Natl. Acad. Sci. USA 2000, 97:8087-8092.
34. Moreno-Bote R., Rinzel J., Rubin N. Noise-induced alternations in an attractor network model of perceptual bistability. J. Neurophysiol. 2007, 98:1125-1139.
35. K. Newhall, G. Kovačič, P. Kramer, A.V. Rangan, D. Cai, Cascade-induced synchrony in stochastically-driven neuronal networks, preprint, 2010.
36. Newhall K., Kovačič G., Kramer P., Zhou D., Rangan A.V., Cai D. Dynamics of current-based, poisson driven, integrate-and-fire neuronal networks. Comm. Math. Sci. 2010, 8:541-600.
37. Nykamp D., Tranchina D. A population density method that facilitates large-scale modeling of neural networks: analysis and application to orientation tuning. J. Comput. Neurosci. 2000, 8:19-50.
38. Nykamp D., Tranchina D. A population density method that facilitates large-scale modeling of neural networks: extension to slow inhibitory synapses. Neural Comput. 2001, 13:511-546.
39. Omurtag A., Kaplan E., Knight B., Sirovich L. A population approach to cortical dynamics with an application to orientation tuning. Network 2000, 11:247-260.
40. Omurtag A., Knight B., Sirovich L. On the simulation of large populations of neurons. J. Comp. Neurosci. 2000, 8:51-63.
41. Pham J., Pakdaman K., Champagnat J., Vibert J. Activity in sparsely connected excitatory neural networks: effect of connectivity neural networks. Neural Networks 1998, 11:415-434.
42. Rangan A.V., Cai D. Maximum-entropy closures for kinetic theories of neuronal network dynamics. Phys. Rev. Lett. 2006, 96:178101.
43. Rangan A.V., Cai D. Fast numerical methods for simulating large-scale integrate-and-fire neuronal networks. J. Comput. Neurosci. 2007, 22(1):81-100.
44. Rangan A.V., Cai D., McLaughlin D.W. Modeling the spatiotemporal cortical activity associated with the line-motion illusion in primary visual cortex. PNAS 2005, 102:18793-18800.
45. Rangan A.V., Cai D., McLaughlin D.W. Quantifying neuronal network dynamics through coarse-grained event trees. PNAS 2008, 105:10990-10995.
46. Rangan A.V., Cai D., Tao L. Numerical methods for solving moment equations in kinetic theory of neuronal network dynamics. J. Comput. Phys. 2007, 221(2):781-798.
47. Rangan A.V., Kovačič G., Cai D. Kinetic theory for neuronal networks with fast and slow excitatory conductances driven by the same spike train. Phys. Rev. E (3) 2008, 77(4):041915.
48. Rudolph M., Destexhe A. Characterization of subthreshold voltage fluctuations in neuronal membranes. Neural Comp. 2003, 15:2577-2618.
49. Shadlen M., Newsome W. The variable discharge of cortical neurons: implications for connectivity, computation and information coding. J. Neurosci. 1998, 18:3870-3896.
50. Shelley M., Tao L. Efficient and accurate time-stepping schemes for integrate-and-fire neuronal networks. J. Comp. Neurosci. 2001, 11:111-119.
51. Shu C.-W., Osher S. Efficient implementation of essentially nonoscillatory shock-capturing schemes. J. Comput. Phys. 1988, 77(2):439-471.
52. Somers D.C., Nelson S.B., Sur M. An emergent model of orientation selectivity in cat visual cortical simple cells. J. Neurosci. 1995, 15(8):5448-5465.
53. Stern E., Kincaid A., Wilson C. Spontaneous subthreshold membrane potential fluctuations and action potential variability of rat corticostriatal and striatal neurons in vivo. J. Neurophysiol. 1997, 77:1697-1715.
54. Tao L., Shelley M., McLaughlin D., Shapley R. An egalitarian network model for the emergence of simple and complex cells in visual cortex. Proc. Natl. Acad. Sci. USA 2004, 101(1):366-371.
55. Treves A. Mean field analysis of neuronal spike dynamics. Network 1993, 4:259-284.
56. Troyer T., Krukowski A., Priebe N., Miller K. Contrast invariant orientation tuning in cat visual cortex with feedforward tuning and correlation based intracortical connectivity. J. Neurosci. 1998, 18:5908-5927.
57. Wang X. Synaptic basis of cortical persistent activity: the importance of NMDA receptors to working memory. J. Neurosci. 1999, 19:9587-9603.
58. Wilbur W., Rinzel J. A theoretical basis for large coefficient of variation and bimodality in neuronal interspike interval distributions. J. Theor. Biol. 1983, 105:345-368.