Nonlinear dynamics and neural systems: synchronization and modeling by Mark Alan Kramer BA (Oberlin College) 2001 A dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Applied Science and Technology in the GRADUATE DIVISION of the UNIVERSITY OF CALIFORNIA, BERKELEY Committee in charge: Professor Andrew J. Szeri, Chair Professor Edgar Knobloch Professor Robert T. Knight Fall 2005 UMI Number: 3210316 Copyright 2005 by Kramer, Mark Alan All rights reserved. UMI UMI Microform 3210316 Copyright 2006 by ProQuest Information and Learning Company.
All rights reserved. This microform edition is protected against unauthorized copying under Title 17, United States Code. ProQuest Information and Learning Company 300 North Zeeb Road P. Box 1346 Ann Arbor, MI 48106-1346 The dissertation of Mark Alan Kramer is approved: Chair Date Date Date University of California, Berkeley Fall 2005 Nonlinear dynamics and neural systems: synchronization and modeling Copyright 2005 by Mark Alan Kramer Abstract Nonlinear dynamics and neural systems: synchronization and modeling by Mark Alan Kramer Doctor of Philosophy in Applied Science and Technology University of California, Berkeley Professor Andrew J.
Szeri, Chair We study the electrical activity of the human cortex in two ways. First, we state seven cou- pling measures to analyze electroencephalogram and electrocortiogram time series. We apply these measures to simulated and observed data, and we use the measures to deduce changes in coupling induced by auditory stimuli and produced by dementia. Second, we define a mathematical model of the spatially averaged, mean-field cortical electrical activity recorded by the electroencephalograph and electrocortiograph.
We compare the model results with ictal electrocortical data collected from four human subjects, and we show that the observed and simulated results agree in two important ways. We use the model to develop three methods for controlling seizures through electrical stim- ulation and to suggest the physiological mechanisms — and points of leverage for therapies — of epilepsy. Szeri Dissertation Committee Chair Contents List of Figures List of Tables 1 Introduction 2 Coupling measures 2. Q Q Q Q Q Q HH HH kg va 2.ẶẶẶ Example: Henonmap.
00002 2 eee eee eee Example: coupled Réssler oscillators .- Example: Oscillatory Bursts .1 Bursting data versus noise. ee Applicaton: auditory ECoG ERPdata. Application: discrimination between healthy and demented subjects .1 Clinical Diagnosis and Data Collecion. ee iii xviii 3 Model 3.
Observational Data: ECoG SeizureRecordings. Model: Dimensionless SPDEs .4 Simulations: Dimensionless ODEs .1 Example: Dimensionless ODES at P„= 11.2 Example: Dimensionless ODEs at Đ„„= 548. Q Q Q Q Q Q Q Q Q Q g n g kg kg k ki kg 3.7 Bifurcation control of the selzing cor@xX. Ặ TQ SH Ra 3.
Control ofstochastic partial differentialequalons .8 Additional routes toseizure 2. ee 4 Conclusion Bibliography ii 71 71 73 76 81 86 90 93 100 103 106 114 130 135 139 153 158 161 163 169 180 184 List of Figures 2.2 (a) The electric potential of one ensemble member recorded by one electrode in the ECoG ERP experiment we discuss in detail in Section 2. The stimulus occurs at t= 0 ms. Note the oscillatory burst between 30 ms and 160 ms.
The asterisk marks a point on the oscillatory burst. We chose d = 3 and t = 10 for illustrative purposes and project the 3-dimensional embedding onto the plane of the page. The asterisk in this figure corresponds to the asterisk in (a). The five nearest neighbors to this point are marked by triangles.
The nearest neighbors are temporally proximal to the fiducial point because the data set is short. (a) Ten ensemble members embedded using d = 3 and t = 10. The thickest curve is the ensemble member shown in Figures 2. The other nine en- semble members are difficult to distinguish.
The point x*[n*] is marked with an asterisk. The thickest curve is the the trajectory of x‘[n]. The point x“[n*] is marked with an asterisk. The thin curves are trajectories of nine other ensemble members.
The nearest ensemble neighbors are marked with triangles.5 (a) The first ensemble members of s*[n] (solid line) and r*[n] (dashed line) gener- ated from the unidirectionally coupled non-identical Henon map for 50s <n < 200 s. The coupling for 100 s <n < 150 s is not obvious. We plot the center time of each window along the horizontal axis, the lag time along the vertical axis, and the value of the cross-correlation in linear greyscale with values greater than 0.8 in black, less than —0.8 in white, and near 0. All WCC figures follow this color scheme unless otherwise indicated.
The WCC reveals a strong correlation between s*[n] and r* [n] at zero lag for 100 s <n < 150s. We plot the center time of each window along the horizontal axis, the frequency along the vertical axis, and the value of the coherence in linear greyscale with values greater than 0.8 in black and near 0. The coher- ence between s“[n] and r“[n] is strong for all frequencies when 100 s <n < 150 Soe ee Computation of the embedding parameters for the time series generated from the unidirectionally coupled non-identical Henon map. (a) The average mutual in- formation (AMI) of the concatenated s*[n] as a function of time lag.
No relative minimum exists for any lag. (b) The percentage of false nearest neighbors of the concatenated s*[m] as a function of the embedding dimension. The value asymp- totes to a small, positive number for dimension 4 and greater. This is due to the random initial values of s*[n] for each ensemble member.
Synchronization measures applied to the unidirectionally coupled non-identical Henon map. All of the measures are smoothed over a window of size 11 at each time point. All three measures increase during the interval of nonlinear coupling (100 s <n < 150 s) between the chaotic time series. (b) The time shifted synchronization measure T(x[n,n]|y) smoothed over a two-dimensional window of size 11 at each time point.
Note that the horizontal and vertical axes show time along ensembles r*[n] and s‘[n], respectively. In the contour plot, there are five evenly spaced contour levels, ranging from 0. Unless defined otherwise, all T(x[{n,7]|y) figures follow this grey-scale scheme. The diag- onal line in the figure corresponds to the location of zero time lag.
The contour plot shows synchronization occurs with time shift 7 = 0 during the time interval 100 s <n< 150 s. (c) The windowed phase synchronization. We plot the center time of each window along the horizontal axis, the phase (in radians) along the vertical axis, and the value of the phase synchronization in linear greyscale, with values greater than 0.1 in black and near 0. A region of strong phase synchro- nization occurs at angles near 0.0 (or equivalently near 27) for 100 s <n < 150 Soo iv 24 27 30 2.9 Example data and linear analysis for time series generated from the coupled Rössler oscillators system.
The coupling for 40 s < ¢ < 60 s is not apparent. The color scheme is the same as that used to create Figure 2. The color scheme is the same as that used to create Figure 2. ee Computation of the embedding parameters for the time series generated from the coupled Réssler oscillators.
(a) The average mutual information (AMD) of the con- catenated s‘[n] as a function of time lag. The first relative minimum occurs at a lag of 1 s. (b) The percentage of false nearest neighbors of the concatenated s*[n] as a function of the embedding dimension. The value asymptotes to a small, positive number for dimensions greater than5.
Synchronization measures applied to the coupled Rössler system. All of the measures are smoothed over a window of size 11 at each time point. Two of the measures H(x{n]|y) and N(x[n]|y) decrease during the known interval of moderate coupling between s*[n] and r*[n]. (b) The time shifted synchronization measure T(x[n,1]|y) smoothed over a two-dimensional window of size 11 at each time point.
Note that the horizontal and vertical axes show time along ensembles r*[n] and s*[n], respectively. The color scheme is the same as in Figure 2. The diagonal line in the figure corresponds to the location of zero time lag. The region of maximum T(x[n,n]|y) begins approximately 4 s above the diagonal.
The color scheme is the same as that used to create Figure 2. An interval of weak phase synchroniza- tion occurs at angles less than 1.0 radian and near 27 radians for 40 s < t < 60 Soo The ensembles of bursting data and noisy data, and the linear coupling measures. The weak oscillatory response in s*[n] be- tween 100 ms and 150 ms is hidden in the noise. The oscillatory response in ensemble s for 100 ms < n < 150 ms is apparent.
The color scheme is the same as that used to create Figure 2. The color scheme is the same as that used to create Figure 2. We find the WC is near zero for all values of frequency and time.13 Synchronization measures applied to the burst versus noise system. All three measures are smoothed over a window of size 11 ms at each time point.
Both H(x{n]|y) and N(x[n]|y) fluctuate between 0.45 over the entire time interval and suggest no obvious synchronization be- tween the ensembles, as expected. S(x{n]|y) increases during the interval 100 ms < n < 140 ms, and therefore suggests an increased synchronization between the en- sembles during this interval. This incorrect interpretation is a consequence of the increase in R(x*[n]) during the oscillatory burst, as explained in the text. The plotting and color scheme are the same as that used to create Figure 2.
The T(x[n,n]ly) result reveals no coupling between the ensembles. (c) The windowed phase synchronization. The plotting and color scheme are the same as that used to create Figure 2. This The ensembles of bursting data and the linear coupling measures.
The weak oscillatory responses of both time series are mostly hidden in the noise. The oscillatory bursts, hidden in the single ensemble member pair of (a), are revealed here in the ensemble averaged ERPs. The color scheme is the same as that used to create Figure 2. The WCC detects moderate cross-correlation for 130 ms < t < 140 ms and a time lag near —6 ms.
The color scheme is the same as that used to create Figure 2. The WC detect strong coherence between s*[n] and r*[n] for 100 ms <¢< 150ms. eee Synchronization measures applied to the ensembles of bursting data. (a) Two syn- chronization measures: H(x{n]|y) (dotted line), and N(x[n]|y) (solid line), smoothed over a window of size 11 ms at each time point.
Neither H(x[n]|y) nor N(x[n]|y) ac- curately captures the synchronization between the two ensemble for 100 ms <n < 160 ms. (b) The time shifted synchronization T(x[n,n]|y) applied to the ensemble of oscillatory bursting data of Figure 2. The plotting and color scheme are the same as that used to create Figure 2. This measure detects the synchro- nization between the two ensembles for 100 ms <n < 135 ms in s*[n] (along the vertical axis) and for 100 ms <n < 135 ms in r* [n] (along the horizontal axis).
(c) The windowed phase synchronization. The plotting and color scheme are the same as that used to create Figure 2. This measure also reveals no coupling between the ensembles. Ặ Ặ Q Q SH Ha The synchronization measure 7'(x{n,\||y) applied to the bursting ensembles of data with weak noise (SNR % 100).
The plotting and color scheme are the same as that used to create Figure 2. A region of strong synchronization occurs for 90 ms <n < 150 ms in r*[n] along the horizontal axis and 100 ms <n < 160 ms in s‘[n] along the vertical axis. Ặ Ặ Q TQ HQ HQ HH va vi 43 45 2.