1 / 50

Scale-up of Cortical Representations in Fluctuation Driven Settings

Scale-up of Cortical Representations in Fluctuation Driven Settings. David W. McLaughlin Courant Institute & Center for Neural Science New York University http://www.cims.nyu.edu/faculty/dmac/ IAS – May‘03. I. Background. Our group has been modeling a local patch (1mm 2 ) of a

Download Presentation

Scale-up of Cortical Representations in Fluctuation Driven Settings

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Scale-up of Cortical Representationsin Fluctuation Driven Settings David W. McLaughlin Courant Institute & Center for Neural Science New York University http://www.cims.nyu.edu/faculty/dmac/ IAS – May‘03

  2. I. Background Our group has been modeling a local patch (1mm2) of a single layer of Primary Visual Cortex Now, we want to “scale-up”   David Cai David Lorentz Robert Shapley Michael Shelley   Louis Tao Jacob Wielaard

  3. Local Patch 1mm2

  4. Lateral Connections and Orientation -- Tree Shrew Bosking, Zhang, Schofield & Fitzpatrick J. Neuroscience, 1997

  5. From one “input layer” to several layers Why scale-up?

  6. II. Our Large-Scale Model • A detailed, fine scale model of a local patch of input layer of Primary Visual Cortex; • Realistically constrained by experimental data; • Integrate & Fire, point neuron model (16,000 neurons per sq mm). Refs: --- PNAS (July, 2000) --- J Neural Science (July, 2001) --- J Comp Neural Sci (2002) --- J Comp Neural Sci (2002) --- http://www.cims.nyu.edu/faculty/dmac/

  7. III. Cortical Networks Have Very Noisy Dynamics Strong temporal fluctuations On synaptic timescale Fluctuation driven spiking

  8. Fluctuation-driven spiking (fluctuations on the synaptic time scale) Solid: average ( over 72 cycles) Dashed: 10 temporal trajectories

  9. IV. Coarse-Grained Asymptotic Representations For “scale-up” in fluctuation driven systems

  10. A Regular Cortical Map ----  500   ----

  11. For the rest of the talk, consider one coarse-grained cell, containing several hundred neurons

  12. We’ll replace the 200 neurons in this CG cell by an effective pdf representation • For convenience of presentation, I’ll sketch the derivation of the reduction for 200 excitatory Integrate & Fire neurons • Later, I’ll show results with inhibition included – as well as “simple” & “complex” cells.

  13. N excitatory neurons (within one CG cell) • all toall coupling; • AMPA synapses (with time scale ) t vi = -(v – VR) – gi (v-VE)  t gi = - gi + l f (t – tl) + (Sa/N) l,k (t – tlk) (g,v,t)  N-1 i=1,N E{[v – vi(t)] [g – gi(t)]}, Expectation “E” taken over modulated incoming (from LGN neurons) Poisson spike train.

  14. t vi = -(v – VR) – gi (v-VE)  t gi = - gi + l f (t – tl) + (Sa/N) l,k (t – tlk) Evolution of pdf -- (g,v,t): t  = -1v {[(v – VR) + g (v-VE)] } + g {(g/) } + 0(t) [(v, g-f/, t) - (v,g,t)] + N m(t) [(v, g-Sa/N, t) - (v,g,t)],where 0(t) = modulated rate of incoming Poisson spike train; m(t) = average firing rate of the neurons in the CG cell =  J(v)(v,g; )|(v= 1) dg where J(v)(v,g; ) = -{[(v – VR) + g (v-VE)] }

  15. t  = -1v {[(v – VR) + g (v-VE)] } + g {(g/) } + 0(t) [(v, g-f/, t) - (v,g,t)] + N m(t) [(v, g-Sa/N, t) - (v,g,t)], N>>1; f << 1; 0 f = O(1); t  = -1v {[(v – VR) + g (v-VE)] } + g {[g – G(t)]/) } + g2/ gg  + … where g2 = 0(t) f2 /(2) + m(t) (Sa)2 /(2N) G(t) = 0(t) f + m(t) Sa

  16. Moments • (g,v,t) • (g)(g,t) =  (g,v,t) dv • (v)(v,t) =  (g,v,t) dg • 1(v)(v,t) =  g (g,tv) dg where (g,v,t) = (g,tv) (v)(v,t) Integrating (g,v,t) eq over v yields:  t (g) =g {[g – G(t)]) (g)} + g2 gg (g)

  17. Fluctuations in g are Gaussian  t (g) =g {[g – G(t)]) (g)} + g2 gg (g)

  18. Integrating (g,v,t) eq over g yields: t (v) = -1v [(v – VR) (v) + 1(v)(v-VE) (v)] Integrating [g (g,v,t)] eq over g yields an equation for 1(v)(v,t) =  g (g,tv) dg, where (g,v,t) = (g,tv) (v)(v,t)

  19. t 1(v) = - -1[1(v) – G(t)] + -1{[(v – VR) + 1(v)(v-VE)] v 1(v)} + 2(v)/ ((v)) v [(v-VE) (v)] + -1(v-VE) v2(v) where 2(v) = 2(v) – (1(v))2 . Closure: (i) v2(v) = 0; (ii) 2(v) = g2

  20. Thus, eqs for (v)(v,t) & 1(v)(v,t): t (v) = -1v [(v – VR) (v)+ 1(v)(v-VE) (v)] t 1(v) = - -1[1(v) – G(t)] + -1{[(v – VR) + 1(v)(v-VE)] v 1(v)} + g2 / ((v)) v [(v-VE) (v)] Together with a diffusion eq for (g)(g,t):  t (g) =g {[g – G(t)]) (g)} + g2 gg (g)

  21. But we can go further for AMPA (  0): t 1(v) = - [1(v) – G(t)] + -1{[(v – VR) + 1(v)(v-VE)] v 1(v)} + g2/ ((v)) v [(v-VE) (v)] Recall g2 = f2/(2) 0(t) + m(t) (Sa)2 /(2N); Thus, as   0, g2 = O(1). Let   0: Algebraically solve for 1(v) : 1(v) = G(t) + g2/ ((v)) v [(v-VE) (v)]

  22. Result: A Fokker-Planck eq for (v)(v,t):  t (v) = v{[(1 + G(t) + g2/ ) v – (VR + VE (G(t) + g2/ ))](v) + g2/ (v- VE)2v(v)} g2/ -- Fluctuations in g Seek steady state solutions – ODE in v, which will be good for scale-up.

  23. Remarks: (i) Boundary Conditions; (ii) Inhibition, spatial coupling of CG cells, simple & complex cells have been added; (iii) N   yields “mean field” representation. Next, use one CG cell to (i) Check accuracy of this pdf representation; (ii) Get insight about mechanisms in fluctuation driven systems.

  24. Bistability and Hysteresis • Network of Simple and Complex — Excitatory only N=16! N=16 Mean­Driven: Fluctuation­Driven: Relatively Strong Cortical Coupling:

  25. Bistability and Hysteresis • Network of Simple and Complex — Excitatory only N=16! Mean­Driven: Relatively Strong Cortical Coupling:

  26. Three Levels of Cortical Amplification : • 1) Weak Cortical Amplification • No Bistability/Hysteresis • 2) Near Critical Cortical Amplification • 3) Strong Cortical Amplification • Bistability/Hysteresis • (2) (1) • (3) • I&F • Excitatory Cells Shown • Possible Mechanism • for Orientation Tuning of Complex Cells • Regime 2 for far-field/well-tuned Complex Cells • Regime 1 for near-pinwheel/less-tuned With inhib, simple & cmplx (2) (1)

  27. New Pdf for fluctuation driven systems -- accurate and efficient • (g,v,t) -- (i) Evolution eq, with jumps from incoming spikes; (ii) Jumps smoothed to diffusion in g by a “large N expansion” • (g)(g,t) =  (g,v,t) dv -- diffuses as a Gaussian • (v)(v,t) =  (g,v,t) dg ; 1(v)(v,t) =  g (g,tv) dg • Coupled (moment) eqs for (v)(v,t) & 1(v)(v,t) , which are not closed [but depend upon 2(v)(v,t) ] • Closure -- (i) v2(v) = 0; (ii) 2(v) = g2 , where 2(v) = 2(v) – (1(v))2 . •   0  eq for 1(v)(v,t) solved algebraically in terms of (v)(v,t), resulting in a Fokker-Planck eq for (v)(v,t)

  28. Scale-up & Dynamical Issuesfor Cortical Modeling of V1 • Temporal emergence of visual perception • Role of spatial & temporal feedback -- within and between cortical layers and regions • Synchrony & asynchrony • Presence (or absence) and role of oscillations • Spike-timing vs firing rate codes • Very noisy, fluctuation driven system • Emergence of an activity dependent, separation of time scales • But often no (or little) temporal scale separation

  29. Fluctuation-driven spiking (very noisy dynamics, on the synaptic time scale) Solid: average ( over 72 cycles) Dashed: 10 temporal trajectories

  30. Simple Complex Expt. Results Shapley’s Lab (unpublished)  4B # of cells  4C CV (Orientation Selectivity)

  31. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Complex Excitatory Cells • Mean-Driven

  32. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Complex Excitatory Cells • Mean-Driven

  33. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Complex Excitatory Cells • Mean-Driven

  34. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Complex Excitatory Cells • Mean-Driven

  35. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Simple Excitatory Cells • Mean-Driven

  36. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Simple Excitatory Cells • Mean-Driven

  37. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Simple Excitatory Cells • Mean-Driven

  38. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Simple Excitatory Cells • Mean-Driven

  39. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Simple Excitatory Cells • Mean-Driven

  40. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Complex Excitatory Cells • Fluctuation-Driven

  41. Incorporation of Inhibitory Cells • 4 Population Dynamics • Simple: • Excitatory • Inhibitory • Complex: • Excitatory • Inhibitory • Simple Excitatory Cells • Fluctuation-Driven

More Related