Modeling Brain Resonance Phenomena Using a Neural Mass Model
Identifieur interne : 001D53 ( Ncbi/Merge ); précédent : 001D52; suivant : 001D54Modeling Brain Resonance Phenomena Using a Neural Mass Model
Auteurs : Andreas Spiegler [Allemagne] ; Thomas R. Knösche [Allemagne] ; Karin Schwab [Allemagne] ; Jens Haueisen [Allemagne] ; Fatihcan M. Atay [Allemagne]Source :
- PLoS Computational Biology [ 1553-734X ] ; 2011.
Abstract
Stimulation with rhythmic light flicker (photic driving) plays an important role in the diagnosis of schizophrenia, mood disorder, migraine, and epilepsy. In particular, the adjustment of spontaneous brain rhythms to the stimulus frequency (entrainment) is used to assess the functional flexibility of the brain. We aim to gain deeper understanding of the mechanisms underlying this technique and to predict the effects of stimulus frequency and intensity. For this purpose, a modified Jansen and Rit neural mass model (NMM) of a cortical circuit is used. This mean field model has been designed to strike a balance between mathematical simplicity and biological plausibility. We reproduced the entrainment phenomenon observed in EEG during a photic driving experiment. More generally, we demonstrate that such a single area model can already yield very complex dynamics, including chaos, for biologically plausible parameter ranges. We chart the entire parameter space by means of characteristic Lyapunov spectra and Kaplan-Yorke dimension as well as time series and power spectra. Rhythmic and chaotic brain states were found virtually next to each other, such that small parameter changes can give rise to switching from one to another. Strikingly, this characteristic pattern of unpredictability generated by the model was matched to the experimental data with reasonable accuracy. These findings confirm that the NMM is a useful model of brain dynamics during photic driving. In this context, it can be used to study the mechanisms of, for example, perception and epileptic seizure generation. In particular, it enabled us to make predictions regarding the stimulus amplitude in further experiments for improving the entrainment effect.
Url:
DOI: 10.1371/journal.pcbi.1002298
PubMed: 22215992
PubMed Central: 3245303
Links toward previous steps (curation, corpus...)
- to stream Pmc, to step Corpus: 002183
- to stream Pmc, to step Curation: 002183
- to stream Pmc, to step Checkpoint: 001A86
Links to Exploration step
PMC:3245303Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Modeling Brain Resonance Phenomena Using a Neural Mass Model</title>
<author><name sortKey="Spiegler, Andreas" sort="Spiegler, Andreas" uniqKey="Spiegler A" first="Andreas" last="Spiegler">Andreas Spiegler</name>
<affiliation wicri:level="3"><nlm:aff id="aff1"><addr-line>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig</wicri:regionArea>
<placeName><region type="land" nuts="1">Saxe (Land)</region>
<region type="district" nuts="2">District de Leipzig</region>
<settlement type="city">Leipzig</settlement>
</placeName>
</affiliation>
<affiliation wicri:level="1"><nlm:aff id="aff2"><addr-line>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau</wicri:regionArea>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Knosche, Thomas R" sort="Knosche, Thomas R" uniqKey="Knosche T" first="Thomas R." last="Knösche">Thomas R. Knösche</name>
<affiliation wicri:level="3"><nlm:aff id="aff1"><addr-line>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig</wicri:regionArea>
<placeName><region type="land" nuts="1">Saxe (Land)</region>
<region type="district" nuts="2">District de Leipzig</region>
<settlement type="city">Leipzig</settlement>
</placeName>
</affiliation>
</author>
<author><name sortKey="Schwab, Karin" sort="Schwab, Karin" uniqKey="Schwab K" first="Karin" last="Schwab">Karin Schwab</name>
<affiliation wicri:level="1"><nlm:aff id="aff3"><addr-line>Bernstein Group for Computational Neuroscience Jena, Institute of Medical Statistics, Computer Sciences and Documentation, Jena University Hospital, Friedrich Schiller University Jena, Jena, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Bernstein Group for Computational Neuroscience Jena, Institute of Medical Statistics, Computer Sciences and Documentation, Jena University Hospital, Friedrich Schiller University Jena, Jena</wicri:regionArea>
<wicri:noRegion>Jena</wicri:noRegion>
<wicri:noRegion>Jena</wicri:noRegion>
<wicri:noRegion>Jena</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Haueisen, Jens" sort="Haueisen, Jens" uniqKey="Haueisen J" first="Jens" last="Haueisen">Jens Haueisen</name>
<affiliation wicri:level="1"><nlm:aff id="aff2"><addr-line>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau</wicri:regionArea>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Atay, Fatihcan M" sort="Atay, Fatihcan M" uniqKey="Atay F" first="Fatihcan M." last="Atay">Fatihcan M. Atay</name>
<affiliation wicri:level="3"><nlm:aff id="aff4"><addr-line>Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Max Planck Institute for Mathematics in the Sciences, Leipzig</wicri:regionArea>
<placeName><region type="land" nuts="1">Saxe (Land)</region>
<region type="district" nuts="2">District de Leipzig</region>
<settlement type="city">Leipzig</settlement>
</placeName>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">22215992</idno>
<idno type="pmc">3245303</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3245303</idno>
<idno type="RBID">PMC:3245303</idno>
<idno type="doi">10.1371/journal.pcbi.1002298</idno>
<date when="2011">2011</date>
<idno type="wicri:Area/Pmc/Corpus">002183</idno>
<idno type="wicri:Area/Pmc/Curation">002183</idno>
<idno type="wicri:Area/Pmc/Checkpoint">001A86</idno>
<idno type="wicri:Area/Ncbi/Merge">001D53</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Modeling Brain Resonance Phenomena Using a Neural Mass Model</title>
<author><name sortKey="Spiegler, Andreas" sort="Spiegler, Andreas" uniqKey="Spiegler A" first="Andreas" last="Spiegler">Andreas Spiegler</name>
<affiliation wicri:level="3"><nlm:aff id="aff1"><addr-line>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig</wicri:regionArea>
<placeName><region type="land" nuts="1">Saxe (Land)</region>
<region type="district" nuts="2">District de Leipzig</region>
<settlement type="city">Leipzig</settlement>
</placeName>
</affiliation>
<affiliation wicri:level="1"><nlm:aff id="aff2"><addr-line>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau</wicri:regionArea>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Knosche, Thomas R" sort="Knosche, Thomas R" uniqKey="Knosche T" first="Thomas R." last="Knösche">Thomas R. Knösche</name>
<affiliation wicri:level="3"><nlm:aff id="aff1"><addr-line>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig</wicri:regionArea>
<placeName><region type="land" nuts="1">Saxe (Land)</region>
<region type="district" nuts="2">District de Leipzig</region>
<settlement type="city">Leipzig</settlement>
</placeName>
</affiliation>
</author>
<author><name sortKey="Schwab, Karin" sort="Schwab, Karin" uniqKey="Schwab K" first="Karin" last="Schwab">Karin Schwab</name>
<affiliation wicri:level="1"><nlm:aff id="aff3"><addr-line>Bernstein Group for Computational Neuroscience Jena, Institute of Medical Statistics, Computer Sciences and Documentation, Jena University Hospital, Friedrich Schiller University Jena, Jena, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Bernstein Group for Computational Neuroscience Jena, Institute of Medical Statistics, Computer Sciences and Documentation, Jena University Hospital, Friedrich Schiller University Jena, Jena</wicri:regionArea>
<wicri:noRegion>Jena</wicri:noRegion>
<wicri:noRegion>Jena</wicri:noRegion>
<wicri:noRegion>Jena</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Haueisen, Jens" sort="Haueisen, Jens" uniqKey="Haueisen J" first="Jens" last="Haueisen">Jens Haueisen</name>
<affiliation wicri:level="1"><nlm:aff id="aff2"><addr-line>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau</wicri:regionArea>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
<wicri:noRegion>Ilmenau</wicri:noRegion>
</affiliation>
</author>
<author><name sortKey="Atay, Fatihcan M" sort="Atay, Fatihcan M" uniqKey="Atay F" first="Fatihcan M." last="Atay">Fatihcan M. Atay</name>
<affiliation wicri:level="3"><nlm:aff id="aff4"><addr-line>Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany</addr-line>
</nlm:aff>
<country xml:lang="fr">Allemagne</country>
<wicri:regionArea>Max Planck Institute for Mathematics in the Sciences, Leipzig</wicri:regionArea>
<placeName><region type="land" nuts="1">Saxe (Land)</region>
<region type="district" nuts="2">District de Leipzig</region>
<settlement type="city">Leipzig</settlement>
</placeName>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS Computational Biology</title>
<idno type="ISSN">1553-734X</idno>
<idno type="eISSN">1553-7358</idno>
<imprint><date when="2011">2011</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Stimulation with rhythmic light flicker (photic driving) plays an important role in the diagnosis of schizophrenia, mood disorder, migraine, and epilepsy. In particular, the adjustment of spontaneous brain rhythms to the stimulus frequency (entrainment) is used to assess the functional flexibility of the brain. We aim to gain deeper understanding of the mechanisms underlying this technique and to predict the effects of stimulus frequency and intensity. For this purpose, a modified Jansen and Rit neural mass model (NMM) of a cortical circuit is used. This mean field model has been designed to strike a balance between mathematical simplicity and biological plausibility. We reproduced the entrainment phenomenon observed in EEG during a photic driving experiment. More generally, we demonstrate that such a single area model can already yield very complex dynamics, including chaos, for biologically plausible parameter ranges. We chart the entire parameter space by means of characteristic Lyapunov spectra and Kaplan-Yorke dimension as well as time series and power spectra. Rhythmic and chaotic brain states were found virtually next to each other, such that small parameter changes can give rise to switching from one to another. Strikingly, this characteristic pattern of unpredictability generated by the model was matched to the experimental data with reasonable accuracy. These findings confirm that the NMM is a useful model of brain dynamics during photic driving. In this context, it can be used to study the mechanisms of, for example, perception and epileptic seizure generation. In particular, it enabled us to make predictions regarding the stimulus amplitude in further experiments for improving the entrainment effect.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Gebber, Gl" uniqKey="Gebber G">GL Gebber</name>
</author>
<author><name sortKey="Zhong, S" uniqKey="Zhong S">S Zhong</name>
</author>
<author><name sortKey="Lewis, C" uniqKey="Lewis C">C Lewis</name>
</author>
<author><name sortKey="Barman, Sm" uniqKey="Barman S">SM Barman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Herrmann, Cs" uniqKey="Herrmann C">CS Herrmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Silberstein, Rb" uniqKey="Silberstein R">RB Silberstein</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Niedermeyer, E" uniqKey="Niedermeyer E">E Niedermeyer</name>
</author>
<author><name sortKey="Lopes Da Silva, Fh" uniqKey="Lopes Da Silva F">FH Lopes da Silva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cobb, S" uniqKey="Cobb S">S Cobb</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jin, Y" uniqKey="Jin Y">Y Jin</name>
</author>
<author><name sortKey="Potkin, Sg" uniqKey="Potkin S">SG Potkin</name>
</author>
<author><name sortKey="Sandman, Ca" uniqKey="Sandman C">CA Sandman</name>
</author>
<author><name sortKey="Bunney, We" uniqKey="Bunney W">WE Bunney</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lazarev, Vv" uniqKey="Lazarev V">VV Lazarev</name>
</author>
<author><name sortKey="Simpson, Dm" uniqKey="Simpson D">DM Simpson</name>
</author>
<author><name sortKey="Schubsky, Bm" uniqKey="Schubsky B">BM Schubsky</name>
</author>
<author><name sortKey="Deazevedo, Lc" uniqKey="Deazevedo L">LC Deazevedo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Takahashi, T" uniqKey="Takahashi T">T Takahashi</name>
</author>
<author><name sortKey="Matsuoka, H" uniqKey="Matsuoka H">H Matsuoka</name>
</author>
<author><name sortKey="Okuma, T" uniqKey="Okuma T">T Okuma</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lazarev, Vv" uniqKey="Lazarev V">VV Lazarev</name>
</author>
<author><name sortKey="Infantosi, Af" uniqKey="Infantosi A">AF Infantosi</name>
</author>
<author><name sortKey="Valencio De Campos, D" uniqKey="Valencio De Campos D">D Valencio-de-Campos</name>
</author>
<author><name sortKey="Deazevedo, Lc" uniqKey="Deazevedo L">LC deAzevedo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miranda De Sa, Am" uniqKey="Miranda De Sa A">AM Miranda de Sa</name>
</author>
<author><name sortKey="Infantosi, Af" uniqKey="Infantosi A">AF Infantosi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fedotchev, Ai" uniqKey="Fedotchev A">AI Fedotchev</name>
</author>
<author><name sortKey="Bondar, At" uniqKey="Bondar A">AT Bondar</name>
</author>
<author><name sortKey="Konovalov, Vf" uniqKey="Konovalov V">VF Konovalov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kawaguchi, T" uniqKey="Kawaguchi T">T Kawaguchi</name>
</author>
<author><name sortKey="Jijiwa, H" uniqKey="Jijiwa H">H Jijiwa</name>
</author>
<author><name sortKey="Watanabe, S" uniqKey="Watanabe S">S Watanabe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rappelsberger, P" uniqKey="Rappelsberger P">P Rappelsberger</name>
</author>
<author><name sortKey="Muller Paschinger, Ib" uniqKey="Muller Paschinger I">IB Muller-Paschinger</name>
</author>
<author><name sortKey="Petsche, H" uniqKey="Petsche H">H Petsche</name>
</author>
<author><name sortKey="Pockberger, H" uniqKey="Pockberger H">H Pockberger</name>
</author>
<author><name sortKey="Prohaska, O" uniqKey="Prohaska O">O Prohaska</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sakamoto, H" uniqKey="Sakamoto H">H Sakamoto</name>
</author>
<author><name sortKey="Inouye, T" uniqKey="Inouye T">T Inouye</name>
</author>
<author><name sortKey="Shinosaki, K" uniqKey="Shinosaki K">K Shinosaki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kalitzin, S" uniqKey="Kalitzin S">S Kalitzin</name>
</author>
<author><name sortKey="Parra, J" uniqKey="Parra J">J Parra</name>
</author>
<author><name sortKey="Velis, Dn" uniqKey="Velis D">DN Velis</name>
</author>
<author><name sortKey="Lopes Da Silva, Fh" uniqKey="Lopes Da Silva F">FH Lopes da Silva</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parra, J" uniqKey="Parra J">J Parra</name>
</author>
<author><name sortKey="Kalitzin, Sn" uniqKey="Kalitzin S">SN Kalitzin</name>
</author>
<author><name sortKey="Iriarte, J" uniqKey="Iriarte J">J Iriarte</name>
</author>
<author><name sortKey="Blanes, W" uniqKey="Blanes W">W Blanes</name>
</author>
<author><name sortKey="Velis, Dn" uniqKey="Velis D">DN Velis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lazarev, Vv" uniqKey="Lazarev V">VV Lazarev</name>
</author>
<author><name sortKey="Pontes, A" uniqKey="Pontes A">A Pontes</name>
</author>
<author><name sortKey="Mitrofanov, Aa" uniqKey="Mitrofanov A">AA Mitrofanov</name>
</author>
<author><name sortKey="Deazevedo, Lc" uniqKey="Deazevedo L">LC deAzevedo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Tommaso, M" uniqKey="De Tommaso M">M de Tommaso</name>
</author>
<author><name sortKey="Stramaglia, S" uniqKey="Stramaglia S">S Stramaglia</name>
</author>
<author><name sortKey="Schoffelen, Jm" uniqKey="Schoffelen J">JM Schoffelen</name>
</author>
<author><name sortKey="Guido, M" uniqKey="Guido M">M Guido</name>
</author>
<author><name sortKey="Libro, G" uniqKey="Libro G">G Libro</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jin, Y" uniqKey="Jin Y">Y Jin</name>
</author>
<author><name sortKey="Potkin, Sg" uniqKey="Potkin S">SG Potkin</name>
</author>
<author><name sortKey="Sandman, Ca" uniqKey="Sandman C">CA Sandman</name>
</author>
<author><name sortKey="Bunney, We" uniqKey="Bunney W">WE Bunney</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bojak, I" uniqKey="Bojak I">I Bojak</name>
</author>
<author><name sortKey="Liley, Dt" uniqKey="Liley D">DT Liley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Breakspear, M" uniqKey="Breakspear M">M Breakspear</name>
</author>
<author><name sortKey="Roberts, Ja" uniqKey="Roberts J">JA Roberts</name>
</author>
<author><name sortKey="Terry, Jr" uniqKey="Terry J">JR Terry</name>
</author>
<author><name sortKey="Rodrigues, S" uniqKey="Rodrigues S">S Rodrigues</name>
</author>
<author><name sortKey="Mahant, N" uniqKey="Mahant N">N Mahant</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Freeman, Wj" uniqKey="Freeman W">WJ Freeman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jansen, Bh" uniqKey="Jansen B">BH Jansen</name>
</author>
<author><name sortKey="Rit, Vg" uniqKey="Rit V">VG Rit</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lopes Da Silva, Fh" uniqKey="Lopes Da Silva F">FH Lopes da Silva</name>
</author>
<author><name sortKey="Hoeks, A" uniqKey="Hoeks A">A Hoeks</name>
</author>
<author><name sortKey="Smits, H" uniqKey="Smits H">H Smits</name>
</author>
<author><name sortKey="Zetterberg, Lh" uniqKey="Zetterberg L">LH Zetterberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Markounikau, V" uniqKey="Markounikau V">V Markounikau</name>
</author>
<author><name sortKey="Igel, C" uniqKey="Igel C">C Igel</name>
</author>
<author><name sortKey="Grinvald, A" uniqKey="Grinvald A">A Grinvald</name>
</author>
<author><name sortKey="Jancke, D" uniqKey="Jancke D">D Jancke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nunez, P" uniqKey="Nunez P">P Nunez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S Coombes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
<author><name sortKey="Robinson, Pa" uniqKey="Robinson P">PA Robinson</name>
</author>
<author><name sortKey="Breakspear, M" uniqKey="Breakspear M">M Breakspear</name>
</author>
<author><name sortKey="Friston, K" uniqKey="Friston K">K Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Abbott, Lf" uniqKey="Abbott L">LF Abbott</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hindmarsh, Jl" uniqKey="Hindmarsh J">JL Hindmarsh</name>
</author>
<author><name sortKey="Rose, Rm" uniqKey="Rose R">RM Rose</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hodgkin, Al" uniqKey="Hodgkin A">AL Hodgkin</name>
</author>
<author><name sortKey="Huxley, Af" uniqKey="Huxley A">AF Huxley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Finger, S" uniqKey="Finger S">S Finger</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shepherd, Gm" uniqKey="Shepherd G">GM Shepherd</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Basar, E" uniqKey="Basar E">E Basar</name>
</author>
<author><name sortKey="Gonder, A" uniqKey="Gonder A">A Gonder</name>
</author>
<author><name sortKey="Ozesmi, C" uniqKey="Ozesmi C">C Ozesmi</name>
</author>
<author><name sortKey="Ungan, P" uniqKey="Ungan P">P Ungan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Basar, E" uniqKey="Basar E">E Basar</name>
</author>
<author><name sortKey="Gonder, A" uniqKey="Gonder A">A Gonder</name>
</author>
<author><name sortKey="Ozesmi, C" uniqKey="Ozesmi C">C Ozesmi</name>
</author>
<author><name sortKey="Ungan, P" uniqKey="Ungan P">P Ungan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bressler, Sl" uniqKey="Bressler S">SL Bressler</name>
</author>
<author><name sortKey="Kelso, Jas" uniqKey="Kelso J">JAS Kelso</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="John, Er" uniqKey="John E">ER John</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Livanov, Mn" uniqKey="Livanov M">MN Livanov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nunez, Pl" uniqKey="Nunez P">PL Nunez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Waters, J" uniqKey="Waters J">J Waters</name>
</author>
<author><name sortKey="Schaefer, A" uniqKey="Schaefer A">A Schaefer</name>
</author>
<author><name sortKey="Sakmann, B" uniqKey="Sakmann B">B Sakmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Contreras, D" uniqKey="Contreras D">D Contreras</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Spruston, N" uniqKey="Spruston N">N Spruston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Frank, Td" uniqKey="Frank T">TD Frank</name>
</author>
<author><name sortKey="Daffertshofer, A" uniqKey="Daffertshofer A">A Daffertshofer</name>
</author>
<author><name sortKey="Beek, Pj" uniqKey="Beek P">PJ Beek</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Harrison, Lm" uniqKey="Harrison L">LM Harrison</name>
</author>
<author><name sortKey="David, O" uniqKey="David O">O David</name>
</author>
<author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Knight, Bw" uniqKey="Knight B">BW Knight</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
<author><name sortKey="Haken, H" uniqKey="Haken H">H Haken</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wilson, Hr" uniqKey="Wilson H">HR Wilson</name>
</author>
<author><name sortKey="Cowan, Jd" uniqKey="Cowan J">JD Cowan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atay, Fm" uniqKey="Atay F">FM Atay</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atay, Fm" uniqKey="Atay F">FM Atay</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jansen, Bh" uniqKey="Jansen B">BH Jansen</name>
</author>
<author><name sortKey="Zouridakis, G" uniqKey="Zouridakis G">G Zouridakis</name>
</author>
<author><name sortKey="Brandt, Me" uniqKey="Brandt M">ME Brandt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lopes Da Silva, Fh" uniqKey="Lopes Da Silva F">FH Lopes da Silva</name>
</author>
<author><name sortKey="Van Rotterdam, A" uniqKey="Van Rotterdam A">A van Rotterdam</name>
</author>
<author><name sortKey="Barts, P" uniqKey="Barts P">P Barts</name>
</author>
<author><name sortKey="Van Heusden, E" uniqKey="Van Heusden E">E van Heusden</name>
</author>
<author><name sortKey="Burr, W" uniqKey="Burr W">W Burr</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zetterberg, Lh" uniqKey="Zetterberg L">LH Zetterberg</name>
</author>
<author><name sortKey="Kristiansson, L" uniqKey="Kristiansson L">L Kristiansson</name>
</author>
<author><name sortKey="Mossberg, K" uniqKey="Mossberg K">K Mossberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thomson, Am" uniqKey="Thomson A">AM Thomson</name>
</author>
<author><name sortKey="Bannister, Ap" uniqKey="Bannister A">AP Bannister</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Haeusler, S" uniqKey="Haeusler S">S Haeusler</name>
</author>
<author><name sortKey="Maass, W" uniqKey="Maass W">W Maass</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="David, O" uniqKey="David O">O David</name>
</author>
<author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wendling, F" uniqKey="Wendling F">F Wendling</name>
</author>
<author><name sortKey="Bartolomei, F" uniqKey="Bartolomei F">F Bartolomei</name>
</author>
<author><name sortKey="Bellanger, Jj" uniqKey="Bellanger J">JJ Bellanger</name>
</author>
<author><name sortKey="Chauvel, P" uniqKey="Chauvel P">P Chauvel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wendling, F" uniqKey="Wendling F">F Wendling</name>
</author>
<author><name sortKey="Bellanger, Jj" uniqKey="Bellanger J">JJ Bellanger</name>
</author>
<author><name sortKey="Bartolomei, F" uniqKey="Bartolomei F">F Bartolomei</name>
</author>
<author><name sortKey="Chauvel, P" uniqKey="Chauvel P">P Chauvel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Grimbert, F" uniqKey="Grimbert F">F Grimbert</name>
</author>
<author><name sortKey="Faugeras, O" uniqKey="Faugeras O">O Faugeras</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Touboul, J" uniqKey="Touboul J">J Touboul</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Spiegler, A" uniqKey="Spiegler A">A Spiegler</name>
</author>
<author><name sortKey="Kiebel, Sj" uniqKey="Kiebel S">SJ Kiebel</name>
</author>
<author><name sortKey="Atay, Fm" uniqKey="Atay F">FM Atay</name>
</author>
<author><name sortKey="Knosche, Tr" uniqKey="Knosche T">TR Knösche</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="David, O" uniqKey="David O">O David</name>
</author>
<author><name sortKey="Harrison, L" uniqKey="Harrison L">L Harrison</name>
</author>
<author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="David, O" uniqKey="David O">O David</name>
</author>
<author><name sortKey="Kiebel, Sj" uniqKey="Kiebel S">SJ Kiebel</name>
</author>
<author><name sortKey="Harrison, Lm" uniqKey="Harrison L">LM Harrison</name>
</author>
<author><name sortKey="Mattout, J" uniqKey="Mattout J">J Mattout</name>
</author>
<author><name sortKey="Kilner, Jm" uniqKey="Kilner J">JM Kilner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sotero, Rc" uniqKey="Sotero R">RC Sotero</name>
</author>
<author><name sortKey="Trujillo Barreto, Nj" uniqKey="Trujillo Barreto N">NJ Trujillo-Barreto</name>
</author>
<author><name sortKey="Iturria Medina, Y" uniqKey="Iturria Medina Y">Y Iturria-Medina</name>
</author>
<author><name sortKey="Carbonell, F" uniqKey="Carbonell F">F Carbonell</name>
</author>
<author><name sortKey="Jimenez, Jc" uniqKey="Jimenez J">JC Jimenez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiebel, Sj" uniqKey="Kiebel S">SJ Kiebel</name>
</author>
<author><name sortKey="Garrido, Mi" uniqKey="Garrido M">MI Garrido</name>
</author>
<author><name sortKey="Moran, R" uniqKey="Moran R">R Moran</name>
</author>
<author><name sortKey="Chen, Cc" uniqKey="Chen C">CC Chen</name>
</author>
<author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiebel, Sj" uniqKey="Kiebel S">SJ Kiebel</name>
</author>
<author><name sortKey="David, O" uniqKey="David O">O David</name>
</author>
<author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiebel, Sj" uniqKey="Kiebel S">SJ Kiebel</name>
</author>
<author><name sortKey="Garrido, Mi" uniqKey="Garrido M">MI Garrido</name>
</author>
<author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Moran, Rj" uniqKey="Moran R">RJ Moran</name>
</author>
<author><name sortKey="Stephan, Ke" uniqKey="Stephan K">KE Stephan</name>
</author>
<author><name sortKey="Seidenbecher, T" uniqKey="Seidenbecher T">T Seidenbecher</name>
</author>
<author><name sortKey="Pape, Hc" uniqKey="Pape H">HC Pape</name>
</author>
<author><name sortKey="Dolan, Rj" uniqKey="Dolan R">RJ Dolan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buzsaki, G" uniqKey="Buzsaki G">G Buzsáki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buzsaki, G" uniqKey="Buzsaki G">G Buzsáki</name>
</author>
<author><name sortKey="Draguhn, A" uniqKey="Draguhn A">A Draguhn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Daneshyari, M" uniqKey="Daneshyari M">M Daneshyari</name>
</author>
<author><name sortKey="Kamkar, Ll" uniqKey="Kamkar L">LL Kamkar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ghosh Dastidar, S" uniqKey="Ghosh Dastidar S">S Ghosh-Dastidar</name>
</author>
<author><name sortKey="Adeli, H" uniqKey="Adeli H">H Adeli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Raiesdana, S" uniqKey="Raiesdana S">S Raiesdana</name>
</author>
<author><name sortKey="Golpayegani, S" uniqKey="Golpayegani S">S Golpayegani</name>
</author>
<author><name sortKey="Firoozabadi, Smp" uniqKey="Firoozabadi S">SMP Firoozabadi</name>
</author>
<author><name sortKey="Habibabadi, Jm" uniqKey="Habibabadi J">JM Habibabadi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buchli, J" uniqKey="Buchli J">J Buchli</name>
</author>
<author><name sortKey="Righetti, L" uniqKey="Righetti L">L Righetti</name>
</author>
<author><name sortKey="Ijspeert, Aj" uniqKey="Ijspeert A">AJ Ijspeert</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pikovsky, A" uniqKey="Pikovsky A">A Pikovsky</name>
</author>
<author><name sortKey="Rosenblum, M" uniqKey="Rosenblum M">M Rosenblum</name>
</author>
<author><name sortKey="Kurths, J" uniqKey="Kurths J">J Kurths</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shil Nikov, Lp" uniqKey="Shil Nikov L">LP Shil'nikov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Beer, Rd" uniqKey="Beer R">RD Beer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Veen, L" uniqKey="Van Veen L">L van Veen</name>
</author>
<author><name sortKey="Liley, Dt" uniqKey="Liley D">DT Liley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rossler, Oe" uniqKey="Rossler O">OE Rössler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schwab, K" uniqKey="Schwab K">K Schwab</name>
</author>
<author><name sortKey="Ligges, C" uniqKey="Ligges C">C Ligges</name>
</author>
<author><name sortKey="Jungmann, T" uniqKey="Jungmann T">T Jungmann</name>
</author>
<author><name sortKey="Hilgenfeld, B" uniqKey="Hilgenfeld B">B Hilgenfeld</name>
</author>
<author><name sortKey="Haueisen, J" uniqKey="Haueisen J">J Haueisen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S Coombes</name>
</author>
<author><name sortKey="Osbaldestin, Ah" uniqKey="Osbaldestin A">AH Osbaldestin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Levi, M" uniqKey="Levi M">M Levi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lofaro, T" uniqKey="Lofaro T">T LoFaro</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Osipov, Gv" uniqKey="Osipov G">GV Osipov</name>
</author>
<author><name sortKey="Kurths, J" uniqKey="Kurths J">J Kurths</name>
</author>
<author><name sortKey="Zhou, C" uniqKey="Zhou C">C Zhou</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vadivasova, Te" uniqKey="Vadivasova T">TE Vadivasova</name>
</author>
<author><name sortKey="Balanov, Ag" uniqKey="Balanov A">AG Balanov</name>
</author>
<author><name sortKey="Sosnovtseva, Ov" uniqKey="Sosnovtseva O">OV Sosnovtseva</name>
</author>
<author><name sortKey="Postnov, De" uniqKey="Postnov D">DE Postnov</name>
</author>
<author><name sortKey="Mosekilde, E" uniqKey="Mosekilde E">E Mosekilde</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brons, M" uniqKey="Brons M">M Brons</name>
</author>
<author><name sortKey="Gross, P" uniqKey="Gross P">P Gross</name>
</author>
<author><name sortKey="Bar Eli, K" uniqKey="Bar Eli K">K Bar-Eli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Meron, E" uniqKey="Meron E">E Meron</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nayfeh, Ah" uniqKey="Nayfeh A">AH Nayfeh</name>
</author>
<author><name sortKey="Sanchez, Ne" uniqKey="Sanchez N">NE Sanchez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parlitz, U" uniqKey="Parlitz U">U Parlitz</name>
</author>
<author><name sortKey="Lauterborn, W" uniqKey="Lauterborn W">W Lauterborn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parlitz, U" uniqKey="Parlitz U">U Parlitz</name>
</author>
<author><name sortKey="Lauterborn, W" uniqKey="Lauterborn W">W Lauterborn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Park, Eh" uniqKey="Park E">EH Park</name>
</author>
<author><name sortKey="Zaks, Ma" uniqKey="Zaks M">MA Zaks</name>
</author>
<author><name sortKey="Kurths, J" uniqKey="Kurths J">J Kurths</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zaks, Ma" uniqKey="Zaks M">MA Zaks</name>
</author>
<author><name sortKey="Park, Eh" uniqKey="Park E">EH Park</name>
</author>
<author><name sortKey="Rosenblum, Mg" uniqKey="Rosenblum M">MG Rosenblum</name>
</author>
<author><name sortKey="Kurths, J" uniqKey="Kurths J">J Kurths</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aihara, K" uniqKey="Aihara K">K Aihara</name>
</author>
<author><name sortKey="Matsumoto, G" uniqKey="Matsumoto G">G Matsumoto</name>
</author>
<author><name sortKey="Ikegaya, Y" uniqKey="Ikegaya Y">Y Ikegaya</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Flaherty, Je" uniqKey="Flaherty J">JE Flaherty</name>
</author>
<author><name sortKey="Hoppensteadt, Fc" uniqKey="Hoppensteadt F">FC Hoppensteadt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sanchez, Ne" uniqKey="Sanchez N">NE Sanchez</name>
</author>
<author><name sortKey="Nayfeh, Ah" uniqKey="Nayfeh A">AH Nayfeh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aizawa, Y" uniqKey="Aizawa Y">Y Aizawa</name>
</author>
<author><name sortKey="Uezu, T" uniqKey="Uezu T">T Uezu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhusubaliyev, Zt" uniqKey="Zhusubaliyev Z">ZT Zhusubaliyev</name>
</author>
<author><name sortKey="Laugesen, Jl" uniqKey="Laugesen J">JL Laugesen</name>
</author>
<author><name sortKey="Mosekilde, E" uniqKey="Mosekilde E">E Mosekilde</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, G R" uniqKey="Wang G">G-R Wang</name>
</author>
<author><name sortKey="Chen, S G" uniqKey="Chen S">S-G Chen</name>
</author>
<author><name sortKey="Hao, B L" uniqKey="Hao B">B-L Hao</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gilbert, T" uniqKey="Gilbert T">T Gilbert</name>
</author>
<author><name sortKey="Gammon, Rw" uniqKey="Gammon R">RW Gammon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mettin, R" uniqKey="Mettin R">R Mettin</name>
</author>
<author><name sortKey="Parlitz, U" uniqKey="Parlitz U">U Parlitz</name>
</author>
<author><name sortKey="Lauterborn, W" uniqKey="Lauterborn W">W Lauterborn</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pikovsky, As" uniqKey="Pikovsky A">AS Pikovsky</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Crevier, Dw" uniqKey="Crevier D">DW Crevier</name>
</author>
<author><name sortKey="Meister, M" uniqKey="Meister M">M Meister</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aihara, K" uniqKey="Aihara K">K Aihara</name>
</author>
<author><name sortKey="Numajiri, T" uniqKey="Numajiri T">T Numajiri</name>
</author>
<author><name sortKey="Matsumoto, G" uniqKey="Matsumoto G">G Matsumoto</name>
</author>
<author><name sortKey="Kotani, M" uniqKey="Kotani M">M Kotani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Engel, Ak" uniqKey="Engel A">AK Engel</name>
</author>
<author><name sortKey="Singer, W" uniqKey="Singer W">W Singer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fries, P" uniqKey="Fries P">P Fries</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schnitzler, A" uniqKey="Schnitzler A">A Schnitzler</name>
</author>
<author><name sortKey="Gross, J" uniqKey="Gross J">J Gross</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tallon Baudry, C" uniqKey="Tallon Baudry C">C Tallon-Baudry</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Varela, F" uniqKey="Varela F">F Varela</name>
</author>
<author><name sortKey="Lachaux, Jp" uniqKey="Lachaux J">JP Lachaux</name>
</author>
<author><name sortKey="Rodriguez, E" uniqKey="Rodriguez E">E Rodriguez</name>
</author>
<author><name sortKey="Martinerie, J" uniqKey="Martinerie J">J Martinerie</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Visani, E" uniqKey="Visani E">E Visani</name>
</author>
<author><name sortKey="Varotto, G" uniqKey="Varotto G">G Varotto</name>
</author>
<author><name sortKey="Binelli, S" uniqKey="Binelli S">S Binelli</name>
</author>
<author><name sortKey="Fratello, L" uniqKey="Fratello L">L Fratello</name>
</author>
<author><name sortKey="Franceschetti, S" uniqKey="Franceschetti S">S Franceschetti</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Siniatchkin, M" uniqKey="Siniatchkin M">M Siniatchkin</name>
</author>
<author><name sortKey="Groppa, S" uniqKey="Groppa S">S Groppa</name>
</author>
<author><name sortKey="Jerosch, B" uniqKey="Jerosch B">B Jerosch</name>
</author>
<author><name sortKey="Muhle, H" uniqKey="Muhle H">H Muhle</name>
</author>
<author><name sortKey="Kurth, C" uniqKey="Kurth C">C Kurth</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Robinson, Pa" uniqKey="Robinson P">PA Robinson</name>
</author>
<author><name sortKey="Rennie, Cj" uniqKey="Rennie C">CJ Rennie</name>
</author>
<author><name sortKey="Rowe, Dl" uniqKey="Rowe D">DL Rowe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dafilis, Mp" uniqKey="Dafilis M">MP Dafilis</name>
</author>
<author><name sortKey="Frascoli, F" uniqKey="Frascoli F">F Frascoli</name>
</author>
<author><name sortKey="Cadusch, Pj" uniqKey="Cadusch P">PJ Cadusch</name>
</author>
<author><name sortKey="Liley, Dtj" uniqKey="Liley D">DTJ Liley</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dafilis, Mp" uniqKey="Dafilis M">MP Dafilis</name>
</author>
<author><name sortKey="Liley, Dtj" uniqKey="Liley D">DTJ Liley</name>
</author>
<author><name sortKey="Cadusch, Pj" uniqKey="Cadusch P">PJ Cadusch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Freeman, Wj" uniqKey="Freeman W">WJ Freeman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Freeman, Wj" uniqKey="Freeman W">WJ Freeman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Freeman, Wj" uniqKey="Freeman W">WJ Freeman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Korn, H" uniqKey="Korn H">H Korn</name>
</author>
<author><name sortKey="Faure, P" uniqKey="Faure P">P Faure</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Makeig, S" uniqKey="Makeig S">S Makeig</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miller, Kj" uniqKey="Miller K">KJ Miller</name>
</author>
<author><name sortKey="Leuthardt, Ec" uniqKey="Leuthardt E">EC Leuthardt</name>
</author>
<author><name sortKey="Schalk, G" uniqKey="Schalk G">G Schalk</name>
</author>
<author><name sortKey="Rao, Rp" uniqKey="Rao R">RP Rao</name>
</author>
<author><name sortKey="Anderson, Nr" uniqKey="Anderson N">NR Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y Zhang</name>
</author>
<author><name sortKey="Wang, X" uniqKey="Wang X">X Wang</name>
</author>
<author><name sortKey="Bressler, Sl" uniqKey="Bressler S">SL Bressler</name>
</author>
<author><name sortKey="Chen, Y" uniqKey="Chen Y">Y Chen</name>
</author>
<author><name sortKey="Ding, M" uniqKey="Ding M">M Ding</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Doesburg, Sm" uniqKey="Doesburg S">SM Doesburg</name>
</author>
<author><name sortKey="Kitajo, K" uniqKey="Kitajo K">K Kitajo</name>
</author>
<author><name sortKey="Ward, Lm" uniqKey="Ward L">LM Ward</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Babloyantz, A" uniqKey="Babloyantz A">A Babloyantz</name>
</author>
<author><name sortKey="Salazar, Jm" uniqKey="Salazar J">JM Salazar</name>
</author>
<author><name sortKey="Nicolis, C" uniqKey="Nicolis C">C Nicolis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gallez, D" uniqKey="Gallez D">D Gallez</name>
</author>
<author><name sortKey="Babloyantz, A" uniqKey="Babloyantz A">A Babloyantz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Di Prisco, Gv" uniqKey="Di Prisco G">GV Di Prisco</name>
</author>
<author><name sortKey="Freeman, Wj" uniqKey="Freeman W">WJ Freeman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Freeman, Wj" uniqKey="Freeman W">WJ Freeman</name>
</author>
<author><name sortKey="Di Prisco, Gv" uniqKey="Di Prisco G">GV Di Prisco</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Babloyantz, A" uniqKey="Babloyantz A">A Babloyantz</name>
</author>
<author><name sortKey="Destexhe, A" uniqKey="Destexhe A">A Destexhe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Iasemidis, Ld" uniqKey="Iasemidis L">LD Iasemidis</name>
</author>
<author><name sortKey="Sackellares, Jc" uniqKey="Sackellares J">JC Sackellares</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Iasemidis, Ld" uniqKey="Iasemidis L">LD Iasemidis</name>
</author>
<author><name sortKey="Sackellares, Jc" uniqKey="Sackellares J">JC Sackellares</name>
</author>
<author><name sortKey="Zaveri, Hp" uniqKey="Zaveri H">HP Zaveri</name>
</author>
<author><name sortKey="Williams, Wj" uniqKey="Williams W">WJ Williams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lehnertz, K" uniqKey="Lehnertz K">K Lehnertz</name>
</author>
<author><name sortKey="Elger, Ce" uniqKey="Elger C">CE Elger</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Steyn Ross, Da" uniqKey="Steyn Ross D">DA Steyn-Ross</name>
</author>
<author><name sortKey="Steyn Ross, Ml" uniqKey="Steyn Ross M">ML Steyn-Ross</name>
</author>
<author><name sortKey="Sleigh, Jw" uniqKey="Sleigh J">JW Sleigh</name>
</author>
<author><name sortKey="Wilson, Mt" uniqKey="Wilson M">MT Wilson</name>
</author>
<author><name sortKey="Gillies, Ip" uniqKey="Gillies I">IP Gillies</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Spiegler, A" uniqKey="Spiegler A">A Spiegler</name>
</author>
<author><name sortKey="Graimann, B" uniqKey="Graimann B">B Graimann</name>
</author>
<author><name sortKey="Pfurtscheller, G" uniqKey="Pfurtscheller G">G Pfurtscheller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Blinowska, Kj" uniqKey="Blinowska K">KJ Blinowska</name>
</author>
<author><name sortKey="Kus, R" uniqKey="Kus R">R Kus</name>
</author>
<author><name sortKey="Kaminski, M" uniqKey="Kaminski M">M Kaminski</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kaminski, M" uniqKey="Kaminski M">M Kaminski</name>
</author>
<author><name sortKey="Ding, M" uniqKey="Ding M">M Ding</name>
</author>
<author><name sortKey="Truccolo, Wa" uniqKey="Truccolo W">WA Truccolo</name>
</author>
<author><name sortKey="Bressler, Sl" uniqKey="Bressler S">SL Bressler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Breakspear, M" uniqKey="Breakspear M">M Breakspear</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Marten, F" uniqKey="Marten F">F Marten</name>
</author>
<author><name sortKey="Rodrigues, S" uniqKey="Rodrigues S">S Rodrigues</name>
</author>
<author><name sortKey="Benjamin, O" uniqKey="Benjamin O">O Benjamin</name>
</author>
<author><name sortKey="Richardson, Mp" uniqKey="Richardson M">MP Richardson</name>
</author>
<author><name sortKey="Terry, Jr" uniqKey="Terry J">JR Terry</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rodrigues, S" uniqKey="Rodrigues S">S Rodrigues</name>
</author>
<author><name sortKey="Barton, D" uniqKey="Barton D">D Barton</name>
</author>
<author><name sortKey="Szalai, R" uniqKey="Szalai R">R Szalai</name>
</author>
<author><name sortKey="Benjamin, O" uniqKey="Benjamin O">O Benjamin</name>
</author>
<author><name sortKey="Richardson, Mp" uniqKey="Richardson M">MP Richardson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watson, Ab" uniqKey="Watson A">AB Watson</name>
</author>
<author><name sortKey="Nachmias, J" uniqKey="Nachmias J">J Nachmias</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wilke, M" uniqKey="Wilke M">M Wilke</name>
</author>
<author><name sortKey="Mueller, Km" uniqKey="Mueller K">KM Mueller</name>
</author>
<author><name sortKey="Leopold, Da" uniqKey="Leopold D">DA Leopold</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Steriade, M" uniqKey="Steriade M">M Steriade</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gabernet, L" uniqKey="Gabernet L">L Gabernet</name>
</author>
<author><name sortKey="Jadhav, Sp" uniqKey="Jadhav S">SP Jadhav</name>
</author>
<author><name sortKey="Feldman, De" uniqKey="Feldman D">DE Feldman</name>
</author>
<author><name sortKey="Carandini, M" uniqKey="Carandini M">M Carandini</name>
</author>
<author><name sortKey="Scanziani, M" uniqKey="Scanziani M">M Scanziani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kimura, F" uniqKey="Kimura F">F Kimura</name>
</author>
<author><name sortKey="Itami, C" uniqKey="Itami C">C Itami</name>
</author>
<author><name sortKey="Ikezoe, K" uniqKey="Ikezoe K">K Ikezoe</name>
</author>
<author><name sortKey="Tamura, H" uniqKey="Tamura H">H Tamura</name>
</author>
<author><name sortKey="Fujita, I" uniqKey="Fujita I">I Fujita</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wehr, M" uniqKey="Wehr M">M Wehr</name>
</author>
<author><name sortKey="Zador, Am" uniqKey="Zador A">AM Zador</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Miller, Kd" uniqKey="Miller K">KD Miller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Priebe, Nj" uniqKey="Priebe N">NJ Priebe</name>
</author>
<author><name sortKey="Ferster, D" uniqKey="Ferster D">D Ferster</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sun, Qq" uniqKey="Sun Q">QQ Sun</name>
</author>
<author><name sortKey="Huguenard, Jr" uniqKey="Huguenard J">JR Huguenard</name>
</author>
<author><name sortKey="Prince, Da" uniqKey="Prince D">DA Prince</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Otnes, Rk" uniqKey="Otnes R">RK Otnes</name>
</author>
<author><name sortKey="Enochson, Ld" uniqKey="Enochson L">LD Enochson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chen, Zm" uniqKey="Chen Z">ZM Chen</name>
</author>
<author><name sortKey="Djidjeli, K" uniqKey="Djidjeli K">K Djidjeli</name>
</author>
<author><name sortKey="Price, Wg" uniqKey="Price W">WG Price</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Baier, G" uniqKey="Baier G">G Baier</name>
</author>
<author><name sortKey="Klein, M" uniqKey="Klein M">M Klein</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kantz, H" uniqKey="Kantz H">H Kantz</name>
</author>
<author><name sortKey="Schreiber, T" uniqKey="Schreiber T">T Schreiber</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jampolsky, A" uniqKey="Jampolsky A">A Jampolsky</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Leijten, Fs" uniqKey="Leijten F">FS Leijten</name>
</author>
<author><name sortKey="Dekker, E" uniqKey="Dekker E">E Dekker</name>
</author>
<author><name sortKey="Spekreijse, H" uniqKey="Spekreijse H">H Spekreijse</name>
</author>
<author><name sortKey="Kasteleijn Nolst Trenite, Dg" uniqKey="Kasteleijn Nolst Trenite D">DG Kasteleijn-Nolst Trenite</name>
</author>
<author><name sortKey="Van Emde Boas, W" uniqKey="Van Emde Boas W">W Van Emde Boas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wolf, A" uniqKey="Wolf A">A Wolf</name>
</author>
<author><name sortKey="Swift, Jb" uniqKey="Swift J">JB Swift</name>
</author>
<author><name sortKey="Swinney, Hl" uniqKey="Swinney H">HL Swinney</name>
</author>
<author><name sortKey="Vastano, Ja" uniqKey="Vastano J">JA Vastano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atay, Fm" uniqKey="Atay F">FM Atay</name>
</author>
<author><name sortKey="Altintas, Y" uniqKey="Altintas Y">Y Altintas</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS Comput Biol</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">ploscomp</journal-id>
<journal-title-group><journal-title>PLoS Computational Biology</journal-title>
</journal-title-group>
<issn pub-type="ppub">1553-734X</issn>
<issn pub-type="epub">1553-7358</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">22215992</article-id>
<article-id pub-id-type="pmc">3245303</article-id>
<article-id pub-id-type="publisher-id">PCOMPBIOL-D-11-00101</article-id>
<article-id pub-id-type="doi">10.1371/journal.pcbi.1002298</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Biology</subject>
<subj-group><subject>Computational Biology</subject>
<subj-group><subject>Computational Neuroscience</subject>
<subj-group><subject>Circuit Models</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group><subject>Neuroscience</subject>
<subj-group><subject>Computational Neuroscience</subject>
<subj-group><subject>Circuit Models</subject>
</subj-group>
</subj-group>
<subj-group><subject>Sensory Systems</subject>
<subj-group><subject>Visual System</subject>
</subj-group>
</subj-group>
<subj-group><subject>Cognitive Neuroscience</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Mathematics</subject>
<subj-group><subject>Nonlinear Dynamics</subject>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Medicine</subject>
<subj-group><subject>Diagnostic Medicine</subject>
<subj-group><subject>Clinical Neurophysiology</subject>
<subj-group><subject>Electroencephalography</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>Modeling Brain Resonance Phenomena Using a Neural Mass Model</article-title>
<alt-title alt-title-type="running-head">Modeling Brain Resonance Phenomena</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Spiegler</surname>
<given-names>Andreas</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Knösche</surname>
<given-names>Thomas R.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Schwab</surname>
<given-names>Karin</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Haueisen</surname>
<given-names>Jens</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Atay</surname>
<given-names>Fatihcan M.</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Institute for Biomedical Engineering and Informatics, Technical University Ilmenau, Ilmenau, Germany</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>Bernstein Group for Computational Neuroscience Jena, Institute of Medical Statistics, Computer Sciences and Documentation, Jena University Hospital, Friedrich Schiller University Jena, Jena, Germany</addr-line>
</aff>
<aff id="aff4"><label>4</label>
<addr-line>Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Sporns</surname>
<given-names>Olaf</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1">Indiana University, United States of America</aff>
<author-notes><corresp id="cor1">* E-mail: <email>spiegler@cbs.mpg.de</email>
</corresp>
<fn fn-type="con"><p>Conceived and designed the experiments: AS TRK JH. Performed the experiments: AS KS. Analyzed the data: AS KS FMA. Wrote the paper: AS TRK KS JH FMA.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><month>12</month>
<year>2011</year>
</pub-date>
<pmc-comment> Fake ppub added to accomodate plos workflow change from 03/2008 and 03/2009 </pmc-comment>
<pub-date pub-type="ppub"><month>12</month>
<year>2011</year>
</pub-date>
<pub-date pub-type="epub"><day>22</day>
<month>12</month>
<year>2011</year>
</pub-date>
<volume>7</volume>
<issue>12</issue>
<elocation-id>e1002298</elocation-id>
<history><date date-type="received"><day>21</day>
<month>1</month>
<year>2011</year>
</date>
<date date-type="accepted"><day>25</day>
<month>10</month>
<year>2011</year>
</date>
</history>
<permissions><copyright-statement>Spiegler et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</copyright-statement>
<copyright-year>2011</copyright-year>
</permissions>
<abstract><p>Stimulation with rhythmic light flicker (photic driving) plays an important role in the diagnosis of schizophrenia, mood disorder, migraine, and epilepsy. In particular, the adjustment of spontaneous brain rhythms to the stimulus frequency (entrainment) is used to assess the functional flexibility of the brain. We aim to gain deeper understanding of the mechanisms underlying this technique and to predict the effects of stimulus frequency and intensity. For this purpose, a modified Jansen and Rit neural mass model (NMM) of a cortical circuit is used. This mean field model has been designed to strike a balance between mathematical simplicity and biological plausibility. We reproduced the entrainment phenomenon observed in EEG during a photic driving experiment. More generally, we demonstrate that such a single area model can already yield very complex dynamics, including chaos, for biologically plausible parameter ranges. We chart the entire parameter space by means of characteristic Lyapunov spectra and Kaplan-Yorke dimension as well as time series and power spectra. Rhythmic and chaotic brain states were found virtually next to each other, such that small parameter changes can give rise to switching from one to another. Strikingly, this characteristic pattern of unpredictability generated by the model was matched to the experimental data with reasonable accuracy. These findings confirm that the NMM is a useful model of brain dynamics during photic driving. In this context, it can be used to study the mechanisms of, for example, perception and epileptic seizure generation. In particular, it enabled us to make predictions regarding the stimulus amplitude in further experiments for improving the entrainment effect.</p>
</abstract>
<abstract abstract-type="summary"><title>Author Summary</title>
<p>Neuroscience aims to understand the enormously complex function of the normal and diseased brain. This, in turn, is the key to explaining human behavior and to developing novel diagnostic and therapeutic procedures. We develop and use models of mean activity in a single brain area, which provide a balance between tractability and plausibility. We use such a model to explain the resonance phenomenon in a photic driving experiment, which is routinely applied in the diagnosis of various diseases including epilepsy, migraine, schizophrenia and depression. Based on the model, we make predictions on the outcome of similar resonance experiments with periodic stimulation of the patients or participants. Our results are important for researchers and clinicians analyzing brain or behavioral data following periodic input.</p>
</abstract>
<counts><page-count count="17"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>Electrophysiological measurements such as <italic>magneto-</italic>
and <italic>electroencephalography</italic>
(M/EEG), <italic>local field potentials</italic>
(LFP) or single unit recordings contain rich information on brain function, which may be related to specific cognitive processes, to general brain states, or to certain pathological conditions. For example, it is known that stimulation by repetitive light flashes entrains the intrinsic alpha EEG rhythm (i.e., frequency entrainment). Neurons in the human visual cortex synchronize their firing to the frequency of flickering light (at rates of about 5 to 30 Hz), causing the EEG alpha frequency to change toward the stimulation frequency <xref ref-type="bibr" rid="pcbi.1002298-Gebber1">[1]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Herrmann1">[2]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Silberstein1">[3]</xref>
. Clinically, this resonance effect is called photic driving. The effect of photic stimulation of the human EEG was first studied in the 1930s and 40s <xref ref-type="bibr" rid="pcbi.1002298-Niedermeyer1">[4]</xref>
. As early as 1947, photic driving was reported in three cases as a potential cause for epileptic activity in patients <xref ref-type="bibr" rid="pcbi.1002298-Cobb1">[5]</xref>
. A review of the clinical routine can be found in Niedermeyer et al. <xref ref-type="bibr" rid="pcbi.1002298-Niedermeyer1">[4]</xref>
. The occurrence of this effect is often interpreted as an indicator for the functional flexibility of the cortex and thus as a sign of healthiness. Today, photic driving is widely used as an activation method in clinical practice, for instance, in epilepsy, migraine, schizophrenia or depression <xref ref-type="bibr" rid="pcbi.1002298-Jin1">[6]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Lazarev1">[7]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Takahashi1">[8]</xref>
. Note, however, that only 50 to 80% of healthy volunteers show a response in the alpha range of EEG <xref ref-type="bibr" rid="pcbi.1002298-Lazarev2">[9]</xref>
, .</p>
<p>Basic properties of the alpha rhythm during photic driving have been investigated by electroencephalographic methods <xref ref-type="bibr" rid="pcbi.1002298-Fedotchev1">[11]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kawaguchi1">[12]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Rappelsberger1">[13]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Sakamoto1">[14]</xref>
. A closer examination of electroencephalographic photic driving effects was given by Herrmann <xref ref-type="bibr" rid="pcbi.1002298-Herrmann1">[2]</xref>
. In that investigation, a flicker stimulus from 1 to 100 Hz in 1-Hz steps was presented. Miranda de Sa and Infantosi <xref ref-type="bibr" rid="pcbi.1002298-MirandadeSa1">[10]</xref>
stimulated at 4, 5, 6, 8, 10, and 12 Hz and showed that stimulation close to the alpha peak was much more effective. The quantification of photic driving from EEG as well as MEG recordings was carried out for the first time by Kalitzin and Parra <xref ref-type="bibr" rid="pcbi.1002298-Kalitzin1">[15]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Parra1">[16]</xref>
. They estimated the phase clustering index of harmonically related frequency components in the EEG and MEG of normal controls and epileptic patients during light stimulations with 10, 15 and 20 Hz. Topographic effects of encephalographic photic driving in the case of children and adolescents were described by Lazarev et al. <xref ref-type="bibr" rid="pcbi.1002298-Lazarev2">[9]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Lazarev3">[17]</xref>
, for patients with migraine by de Tommaso et al.<xref ref-type="bibr" rid="pcbi.1002298-deTommaso1">[18]</xref>
, and for patients with schizophrenia by Jin et al. <xref ref-type="bibr" rid="pcbi.1002298-Jin2">[19]</xref>
.</p>
<p>In order to gain further insight into mechanisms underlying such brain resonance effects and their relevance to brain function and pathology, as well as to make predictions concerning the stimulation parameters, generative models can be used. Such models are called biologically plausible if their state variables and parameters are biophysically meaningful. By fitting the model parameters to measurements, one can test hypotheses on the implementation of brain function. To ensure that this inversion is mathematically tractable and at the same time physically meaningful, the model must strike a balance between mathematical simplicity and biological realism. One class of models designed to meet these criteria is referred to as <italic>neural mass models</italic>
(NMMs) (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Bojak1">[20]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Breakspear1">[21]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Freeman1">[22]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-LopesdaSilva1">[24]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Markounikau1">[25]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Nunez1">[26]</xref>
). NMMs describe neural function at a mesoscopic level <xref ref-type="bibr" rid="pcbi.1002298-Coombes1">[27]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Deco1">[28]</xref>
, in contrast to single neuron models such as simple integrate-and-fire models <xref ref-type="bibr" rid="pcbi.1002298-Abbott1">[29]</xref>
and the more elaborate Hodgkin and Huxley type of models (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Hindmarsh1">[30]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Hodgkin1">[31]</xref>
). NMMs quantify the mean firing rates and mean postsynaptic potentials (PSPs) of neuronal populations, the <italic>neural masses</italic>
(NMs). Although, at the microscopic level, single neurons are considered the primary computational units of the brain's architecture <xref ref-type="bibr" rid="pcbi.1002298-Finger1">[32]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Shepherd1">[33]</xref>
, it is also widely accepted that relevant information processing underlying brain function in both healthy and diseased states can be carried out by ensembles of interacting neurons at the mesoscopic level (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Coombes1">[27]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Deco1">[28]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Basar1">[34]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Basar2">[35]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Bressler1">[36]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jirsa1">[37]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-McIntosh1">[38]</xref>
). In other words, NMMs describe brain activity on a scale that is highly relevant to brain function <xref ref-type="bibr" rid="pcbi.1002298-John1">[39]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Livanov1">[40]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Nunez2">[41]</xref>
. Moreover, when EEG or MEG data are used, NMMs have also the advantage that they predict exactly what is measured by these modalities, namely coherent activity in entire populations of neurons.</p>
<p>However, this type of modeling also involves a number of simplifications that may lead to limitations. First of all, it is based on a simplified notion of the function of a neuron, namely the firing rate model: The neuron convolves the rate of incoming spikes with an alpha-shaped function and thereby generates a change in membrane potential (PSP), and produces an output spike rate that is a non-linear (e.g., sigmoid) function of the membrane potential. These are the most important aspects of neuronal function. However, in the brain, things are usually more complicated. For examples, , modeling is made difficult due to feedback influence of action potentials on the dendritic membrane potentials (back propagation) <xref ref-type="bibr" rid="pcbi.1002298-Waters1">[42]</xref>
, specific intrinsic firing patterns (e.g., bursting) <xref ref-type="bibr" rid="pcbi.1002298-Contreras1">[43]</xref>
and dendritic hierarchies <xref ref-type="bibr" rid="pcbi.1002298-Spruston1">[44]</xref>
. It remains to be investigated if and to what extent such physiological details affect the properties of NMs at the mesoscopic level. A second simplification that leads to limitations is that spike time dependent effects will be missed since the model relies on firing rates rather than on actual spikes. Third, as the distributions of the neural states are simply described by their means, the impact of higher statistical moments is ignored. In order to capture the variability within a NM, one may use the Fokker-Planck formalism <xref ref-type="bibr" rid="pcbi.1002298-Frank1">[45]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Harrison1">[46]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Knight1">[47]</xref>
. Finally, NMMs approximate the spatial scale of neuronal populations to be point-like <xref ref-type="bibr" rid="pcbi.1002298-Freeman1">[22]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-LopesdaSilva1">[24]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Deco1">[28]</xref>
, ignoring the domain of spatial dynamics. In that line, the approach can be generalized, leading to neural field models <xref ref-type="bibr" rid="pcbi.1002298-Bojak1">[20]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Markounikau1">[25]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Nunez1">[26]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Coombes1">[27]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Deco1">[28]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jirsa2">[48]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Wilson1">[49]</xref>
, which take into account the spatial extent of neural circuitry by dealing with aggregated activities in the vicinity of a given location. This puts field theories somewhere between neural mass theories and discrete neuronal networks, allowing them to address, for instance, distance-dependent delays. A quantitative analysis of neural field models can be found in Atay and Hutt <xref ref-type="bibr" rid="pcbi.1002298-Atay1">[50]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Atay2">[51]</xref>
, for example.</p>
<p>In this work, we use a particular local network of NMs first described by Jansen and Rit <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen2">[52]</xref>
, based on earlier works of Lopes da Silva et al. <xref ref-type="bibr" rid="pcbi.1002298-LopesdaSilva1">[24]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-LopesdaSilva2">[53]</xref>
and Zetterberg et al. <xref ref-type="bibr" rid="pcbi.1002298-Zetterberg1">[54]</xref>
. This NMM comprises an elementary circuit of three interconnected NMs (i.e., pyramidal cells and excitatory and inhibitory interneurons) meant to account for a cortical area, such as the primary visual cortex in our photic driving experiment. Although local neuronal circuits can be very complex <xref ref-type="bibr" rid="pcbi.1002298-Thomson1">[55]</xref>
and may be modeled using more than three NMs (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Haeusler1">[56]</xref>
), the circuit used here is the most reduced representation of the features that are relevant for the temporal dynamics, that is, positive and negative feedback loops. The Jansen and Rit structure has been shown to account for both oscillatory <xref ref-type="bibr" rid="pcbi.1002298-David1">[57]</xref>
and seizure-like EEG recordings <xref ref-type="bibr" rid="pcbi.1002298-Wendling1">[58]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Wendling2">[59]</xref>
. Its dynamic behavior, in terms of stabilities and bifurcations, was first characterized by Grimbert and Faugeras <xref ref-type="bibr" rid="pcbi.1002298-Grimbert1">[60]</xref>
and, more generally, by Touboul <xref ref-type="bibr" rid="pcbi.1002298-Touboul1">[61]</xref>
and Spiegler et al. <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. Several such NMMs can be combined to describe networks of coupled cortical areas and account for more complex transient and oscillatory behaviors <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-David1">[57]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Wendling2">[59]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-David2">[63]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-David3">[64]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Sotero1">[65]</xref>
. The Bayesian inversion of such network NMMs given M/EEG data (referred to as <italic>dynamic causal modeling</italic>
(DCM) <xref ref-type="bibr" rid="pcbi.1002298-David3">[64]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kiebel1">[66]</xref>
) has been successfully used for the analysis of event-related <xref ref-type="bibr" rid="pcbi.1002298-David3">[64]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kiebel2">[67]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kiebel3">[68]</xref>
and steady-state responses <xref ref-type="bibr" rid="pcbi.1002298-Moran1">[69]</xref>
.</p>
<p>To date, the dynamics of this system has been systematically investigated only under the assumption of constant extrinsic input levels, thereby allowing the system to settle in a stable state (e.g., fixed point or limit cycle) <xref ref-type="bibr" rid="pcbi.1002298-Grimbert1">[60]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Touboul1">[61]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. However, in a photic driving experiment, one has to consider rhythmic input. Moreover, the model's response to such input is also of great importance in many other settings, since, in the brain, such local neural circuits are embedded in global brain networks and may experience high amplitude time-varying input from other parts of the brain. Because neuronal ensembles tend to oscillate intrinsically, such input is very often periodic, as evidenced by the widespread occurrence of rhythmic activity in both extracranial and intracranial recordings <xref ref-type="bibr" rid="pcbi.1002298-Buzski1">[70]</xref>
.</p>
<p>In this paper, we use a continuous-time periodic function as model input approximating a periodic train of pulses. In this continuous function, each single pulse is similar (but not equal) to the single event used by Jansen and Rit for eliciting visual evoked potentials <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen2">[52]</xref>
, or used in dynamic causal modeling (e.g., <xref ref-type="bibr" rid="pcbi.1002298-David2">[63]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kiebel1">[66]</xref>
). We systematically vary both amplitude (intensity) and frequency of the stimulation within the effective ranges provided by Spiegler et al. <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. We find the frequency entrainment effect spreading over broader stimulus frequencies for higher stimulus intensities, while away from the entrainment ranges, we find complex behavior, including periodic, quasi-periodic, and chaotic dynamics. The latter behavior, in particular, provides continuous spectra. Networks of such chaotic NMMs (incorporating network variability, for example, by different characteristic constants of time and potential) can be used to describe colored noise sources that produce continuous portions in the spectra, such as 1/<italic>f</italic>
-characteristics, that are commonly observed in M/EEG or LFP data <xref ref-type="bibr" rid="pcbi.1002298-Buzski2">[71]</xref>
. Finally, we fit the output of the periodically forced NMM to data from the photic driving experiment in terms of the largest Lyapunov exponent and frequency detuning. The largest Lyapunov exponent measures the exponential separation or convergence of nearby trajectories. It thereby quantifies the predictability or, at the other extreme, the chaoticity of the behavior of the system and has been demonstrated to be an important marker for pathologically altered brain dynamics, especially in epilepsy <xref ref-type="bibr" rid="pcbi.1002298-Daneshyari1">[72]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-GhoshDastidar1">[73]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Raiesdana1">[74]</xref>
. In this way, we show that the NMM is a suitable model for the dynamics of brain resonance phenomena at the cortical level and demonstrate that useful predictions concerning the parameter choice of entrainment experiments can be derived.</p>
<p>To our knowledge, this is the first study to investigate a photic driving experiment using a NMM. We demonstrate that with this NMM, one can explain effects of complex behavior in such an experiment. The results also indicate that a relatively simple model of a local neural circuit is capable of producing surprisingly complex and diverse phenomena, which are observable in brain data and relevant to the explanation of brain function.</p>
</sec>
<sec id="s2"><title>Results</title>
<p>In our previous work on the extended Jansen and Rit neural mass model (NMM) for a cortical area <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
, we found a self-sustained oscillation due to a stable limit cycle with a certain intrinsic frequency for constant input. Forcing such a limit cycle with periodic input to the NMM causes accelerations and/or decelerations of the oscillation, depending on the timing. If their cumulative effect is non-zero, entrainment occurs. For more details on the precise mechanism of entrainment effects, see <xref ref-type="bibr" rid="pcbi.1002298-Buchli1">[75]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Pikovsky1">[76]</xref>
.</p>
<p>Indeed, we observe frequency entrainment, that is, the cortical area responds with the stimulus frequency instead of the intrinsic frequency, thus forming a plateau in the frequency-detuning curves (see colored ranges in <xref ref-type="fig" rid="pcbi-1002298-g001">Figure 1</xref>
). The detuning curve shows intrinsic behavior that is characterized by a typical repetitive s-shape. With increasing frequency, this s-shape becomes more pronounced. For the model with the stimulus amplitude that fits the experimental data best, this s-shape pattern of the detuning-curve is frequently interrupted near the stimulus hitting the intrinsic frequency 0.5≤<italic>η</italic>
/<italic>η</italic>
<sub>int</sub>
≤1.5 by complex behavior such as chaos (see below, as well as <xref ref-type="fig" rid="pcbi-1002298-g001">Figure 1</xref>
and <xref ref-type="fig" rid="pcbi-1002298-g002">Figure 2</xref>
). Apart from the interruption of the repetitive s-shape pattern by irregularities around the intrinsic frequency <italic>η</italic>
<sub>int</sub>
, the general trend of frequency detuning seems to be shifted towards the intrinsic frequency <italic>η</italic>
<sub>int</sub>
by the response frequencies <italic>η</italic>
<sub>resp</sub>
, which explains the experimental data (see B in <xref ref-type="fig" rid="pcbi-1002298-g002">Figure 2</xref>
). Moreover, stimulating near the intrinsic frequency the response matches the stimulation frequency and entrainment occurs (see <xref ref-type="fig" rid="pcbi-1002298-g001">Figure 1</xref>
and <xref ref-type="fig" rid="pcbi-1002298-g002">Figure 2</xref>
).</p>
<fig id="pcbi-1002298-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g001</object-id>
<label>Figure 1</label>
<caption><title>Frequency entrainment effects in a periodically forced neural mass model of a cortical area.</title>
<p>A frequency-detuning curve refers to the ratio of stimulus to characteristic mean response frequency plotted against the ratio of stimulus to intrinsic frequency, for the normalized stimulus amplitude of 1.5. The entrainment ranges around the intrinsic frequency and its subharmonics are shown in red and in blue/green. Such an entrainment around the intrinsic frequency can be found in photic driving experiments (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Schwab1">[81]</xref>
).</p>
</caption>
<graphic xlink:href="pcbi.1002298.g001"></graphic>
</fig>
<fig id="pcbi-1002298-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g002</object-id>
<label>Figure 2</label>
<caption><title>Entrainment effect found in the experimental data and model.</title>
<p>For the experiment, the mean over subjects is shown as white lines and the region between the 5% and 95% quantiles is covered by red areas. For the model (black dots), the amplitude configuration that best fits the largest Lyapunov exponents of the experiments is used (see Comparison section in <xref ref-type="sec" rid="s2">Materials and Methods</xref>
). The entrainment effect is shown for <bold>A</bold>
the stimulus frequency range of the model and <bold>B</bold>
the stimulus range used in the experiments (green area in A).</p>
</caption>
<graphic xlink:href="pcbi.1002298.g002"></graphic>
</fig>
<p>Outside the entrainment ranges, more complex dynamics, including periodic, quasi-periodic and chaotic behavior, are observed (see <xref ref-type="fig" rid="pcbi-1002298-g003">Figure 3</xref>
). Periodic and quasi-periodic behavior are associated with discrete power spectra with frequency peaks <italic>η</italic>
<sub>i</sub>
that are commensurable (i.e., ∑<italic>k</italic>
<sub>i</sub>
<italic>η</italic>
<sub>i</sub>
= 0 for some non-zero integers <italic>k</italic>
<sub>i</sub>
) for the periodic state and incommensurable (i.e., ∑<italic>k</italic>
<sub>i</sub>
<italic>η</italic>
<sub>i</sub>
≠0 for any set of non-zero integers <italic>k</italic>
<sub>i</sub>
) for the quasi-periodic case. Chaotic behavior is indicated by non-closed bounded trajectories in state space, broadband continuous spectra and positive Lyapunov exponents (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
). Here, chaotic regimes arise by traversing a homoclinic Shil'nikov bifurcation (see Figure 13 in <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
) for non-rational ratios between the frequencies of the stimulation and the intrinsic model kinetics. This route to chaos <xref ref-type="bibr" rid="pcbi.1002298-Shilnikov1">[77]</xref>
has also been identified in more theoretical neural models (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Beer1">[78]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-vanVeen1">[79]</xref>
). Arnol'd tongues or mode-locking structures (i.e., entrainment regions in the parameter space <xref ref-type="bibr" rid="pcbi.1002298-Pikovsky1">[76]</xref>
) are apparent in <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
as a result of negative largest Lyapunov exponents. At low amplitudes, we observe several distinct ranges of such mode locking, which seem to merge or overlap at higher amplitudes. Note that chaotic “islands” occur at incommensurable ratios between stimulation frequencies and intrinsic limit cycles and interrupt frequency locking. For example, at a stimulus amplitude of <italic>ζ</italic>
= 0.8, entrainment occurs for stimulus frequencies between 0<<italic>η</italic>
≤0.06075, 0.06831≤<italic>η</italic>
≤0.07403 and 0.09474≤<italic>η</italic>
≤0.1206 interrupted by chaotic regimes between 0.004835≤<italic>η</italic>
≤0.03464. At a stimulus amplitude of <italic>ζ</italic>
= 2.4, entrainment occurs for stimulus frequencies between 0<<italic>η</italic>
≤0.1447, 0.1545≤<italic>η</italic>
≤0.1608 and <italic>η</italic>
>0.1749 interrupted by chaotic regimes between 0.0365≤<italic>η</italic>
≤0.04397, 0.05237≤<italic>η</italic>
≤0.06818 and 0.08906≤<italic>η</italic>
≤0.09547. Note that these entrainment ranges are rough estimates due to the finite sampling of the parameter space and due to the occurring “islands” of chaos. The chaotic regimes that are present in the parameter space feature a single positive largest Lyapunov exponent that is equal to the entropy of the attracting set (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
).</p>
<fig id="pcbi-1002298-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g003</object-id>
<label>Figure 3</label>
<caption><title>Complex behavior occurring in the periodically forced neural mass model of a single area.</title>
<p>Orbits, time series, and power spectra (columns) are shown for three configurations (rows) displaying (top-down) periodic (normalized input amplitude; normalized input frequency: 3.6301; 9.33·10<sup>−2</sup>
), quasi-periodic (1.5; 7.59·10<sup>−2</sup>
) and chaotic behavior (3.6301; 7.05·10<sup>−2</sup>
). The orbits are in the state space of normalized postsynaptic potentials of pyramidal cells caused by both interneurons (<italic>x</italic>
<sub>30</sub>
) as well as at both excitatory and inhibitory interneurons caused by pyramidal cells (<italic>x</italic>
<sub>31</sub>
and <italic>x</italic>
<sub>32</sub>
). The red circle represents the stable limit cycle (i.e., harmonic oscillation) arising from Andronov-Hopf bifurcations performed by the unperturbed system. The time series and the power spectra are shown for the normalized postsynaptic potentials of pyramidal cells (which are related to M/EEG). Periodic behavior is characterized by a closed orbit (limit cycle) and a discrete power spectrum with peaks at commensurable frequencies. Quasi-periodic behavior is characterized by trajectories forming an invariant <italic>n</italic>
-dimensional torus and discrete power spectra with peaks at incommensurable frequencies. Chaotic behavior is indicated by a strange attractor, that is, a bounded attracting set in which all trajectories are unstable and nearby trajectories diverge locally from each other exponentially, and a broadband power spectrum.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g003"></graphic>
</fig>
<fig id="pcbi-1002298-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g004</object-id>
<label>Figure 4</label>
<caption><title>Largest Lyapunov exponent in parameter space.</title>
<p>The map shows the largest Lyapunov exponent <italic>λ</italic>
<sub>1</sub>
as a function of stimulus amplitude and frequency, indicating the sensitivity of the periodically forced neural mass model to initial conditions. Positive exponents (magenta to black) reflect diverging trajectories irrespective of how close they are, and thus chaos in the system. They also measure the entropy of an attracting set for all cases because the maximum number of positive Lyapunov exponents for any parameter configuration is one. Zero exponents (white) indicate neutral stability, and negative exponents (cyan to yellow) reflect frequency locking. Arnol'd tongue structures (i.e., resonance zones) are indicated by negative largest Lyapunov exponents due to the phase locking between system kinetics and the stimulus. The red arrow indicates the amplitude for which the experimental data fits best (see <xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
and Comparison section in <xref ref-type="sec" rid="s2">Materials and Methods</xref>
). Several parameter regions with scattered, presumably fractal patterns of chaotic regime are shown at a finer resolution of normalized stimulus amplitude and frequency.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g004"></graphic>
</fig>
<p>By studying the Lyapunov spectra, configurations are discovered where the system has two zero Lyapunov exponents and evolves on a two-dimensional invariant torus, indicating quasi- and bi-periodicity (see <xref ref-type="fig" rid="pcbi-1002298-g005">Figure 5</xref>
). In general, the model is dissipative (i.e., the sum of Lyapunov exponents is negative) and does not exhibit hyperchaos, which is a higher form of chaos with at least two directions of hyperbolic instability on the attractor <xref ref-type="bibr" rid="pcbi.1002298-Rssler1">[80]</xref>
(see <xref ref-type="sec" rid="s2">Materials and Methods</xref>
for further explanation), as seen from the observation that the second largest Lyapunov exponent is non-positive and the Kaplan-Yorke dimension (see <xref ref-type="sec" rid="s2">Materials and Methods</xref>
) never reaches or exceeds the value of two (see <xref ref-type="fig" rid="pcbi-1002298-g006">Figure 6</xref>
). This means that the dynamics are low dimensional, not only for the periodic, but also for the chaotic regimes, as compared to the dimensionality of the system (which is six plus one dimension for the force). In general, the maximum Kaplan-Yorke dimension is a non-integer because of the complex geometry of the attractor. The periodic forcing seems to work mostly in the direction of entrainment, and although there are occasional “islands” of chaotic regimes, the regular forcing does not let the dynamics become exceedingly chaotic.</p>
<fig id="pcbi-1002298-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g005</object-id>
<label>Figure 5</label>
<caption><title>Occurrence of quasi-periodic behavior forming a two-torus surface in state space.</title>
<p>Two-dimensional tori are indicated by two zero Lyapunov exponents (shown in black dots) in the parameter space of stimulus amplitude and frequency. The red arrow indicates the amplitude for which the experimental data fits best (see <xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
and Comparison section in <xref ref-type="sec" rid="s2">Materials and Methods</xref>
). The parameter regions of chaotic patterns that were selected for recomputing at a finer resolution (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
and <xref ref-type="fig" rid="pcbi-1002298-g006">Figure 6</xref>
) appear here mostly as white areas.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g005"></graphic>
</fig>
<fig id="pcbi-1002298-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g006</object-id>
<label>Figure 6</label>
<caption><title>Kaplan-Yorke dimension of the periodically forced neural mass model in parameter space.</title>
<p>The Kaplan-Yorke dimension <italic>D</italic>
<sub>KY</sub>
given by Equations (8) and (9) never goes above 1.7, thus hyperchaos does not exist in the model. The red arrow indicates the amplitude for which the experimental data fits best (see <xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
and Comparison section in <xref ref-type="sec" rid="s2">Materials and Methods</xref>
). Several parameter regions with scattered, presumably fractal, patterns of chaotic regime were selected for recomputing at a finer resolution of normalized stimulus amplitude and frequency.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g006"></graphic>
</fig>
<p>Furthermore, we find that the model is indeed able to explain frequency entrainment that is observable during a photic driving experiment (see also <xref ref-type="bibr" rid="pcbi.1002298-Gebber1">[1]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Herrmann1">[2]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Schwab1">[81]</xref>
). Note that Figure 5 in Schwab et al. <xref ref-type="bibr" rid="pcbi.1002298-Schwab1">[81]</xref>
contains an error in the labeling of the <italic>y</italic>
-axes. Each graph in this figure correctly plots the ratio of stimulus to response frequency (<italic>y</italic>
-axis) against the ratio of stimulus to alpha frequency (<italic>x</italic>
-axis). In this case, a horizontal line indicates an entrainment effect, while absence of entrainment would result in a diagonal line. We estimated the largest Lyapunov exponents from the data (see Methods section). In order to probe the stability of this estimate, we repeated it with the same data after adding various levels of Gaussian noise. The pattern of the Lyapunov exponents as function of stimulus frequency appears to be quite stable except for very low signal-to-noise-ratios SNR≤3 dB (see supplementary <xref ref-type="supplementary-material" rid="pcbi.1002298.s001">Figure S1</xref>
). We compare our model outcome with these experimental Lyapunov exponents (see the section Experimental data in the <xref ref-type="sec" rid="s2">Materials and Methods</xref>
section and <xref ref-type="fig" rid="pcbi-1002298-g007">Figure 7</xref>
for the experimental paradigm) and find a particular stimulus amplitude for which, for all ten subjects, the model predicts Lyaponov exponents that are in close agreement with those estimated from the experimental data (<xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
and <xref ref-type="table" rid="pcbi-1002298-t001">Table 1</xref>
), with the amplitude being close to <italic>ζ</italic>
= 3.6692 for all ten subjects. In seven of the subjects, the correlation between model prediction and measurement over stimulus frequencies was significant (<italic>p</italic>
<0.05, corrected). A bootstrap test yielded a probability of error (significance) for the mean over subjects of 6.2% (see also <xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8 (B)</xref>
and <xref ref-type="table" rid="pcbi-1002298-t001">Table 1</xref>
). For three subjects (numbers 3, 6 and 7), the individual fit was not significant (see <xref ref-type="table" rid="pcbi-1002298-t001">Table 1</xref>
). Interestingly, this is reflected in the means and the standard deviations of the shift-and-scale parameters <italic>u</italic>
and <italic>v</italic>
(see the section Comparison in the <xref ref-type="sec" rid="s2">Materials and Methods</xref>
section and supplementary <xref ref-type="supplementary-material" rid="pcbi.1002298.s002">Figure S2</xref>
). For the corresponding model configuration, we present a compact representation in <xref ref-type="fig" rid="pcbi-1002298-g009">Figure 9</xref>
and describe the system states qualitatively in <xref ref-type="table" rid="pcbi-1002298-t002">Table 2</xref>
. In the range of stimulus frequency <italic>η</italic>
between 0 and 0.0534, the system performs limit cycles and appears to undergo a cascade of period-adding bifurcations <xref ref-type="bibr" rid="pcbi.1002298-Coombes2">[82]</xref>
with descending stimulus frequency. This local bifurcation consists of saddle-node bifurcations in which a (<italic>n</italic>
+1)-periodic orbit arises out of a <italic>n</italic>
-periodic orbit for <italic>n</italic>
∈ N<sub>1</sub>
<xref ref-type="bibr" rid="pcbi.1002298-Levi1">[83]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-LoFaro1">[84]</xref>
. For stimulus frequencies <italic>η</italic>
above 0.0534, ranges of chaotic, periodic and quasi-periodic behavior occur. Due to the high dimensionality of the system, an instructive presentation in the form of a video is available, comprising orbits (PSPs), time series, and power spectra (see <xref ref-type="supplementary-material" rid="pcbi.1002298.s003">Video S1</xref>
).</p>
<fig id="pcbi-1002298-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g007</object-id>
<label>Figure 7</label>
<caption><title>Experimental design of the flicker stimulation study.</title>
<p>The LEDs were powered for half of each period. The raise and decay time for the LEDs was measured to be 100 µs.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g007"></graphic>
</fig>
<fig id="pcbi-1002298-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g008</object-id>
<label>Figure 8</label>
<caption><title>Comparison of model and data from a photic driving experiment.</title>
<p>The largest normalized Lyapunov exponents calculated from the model show very good agreement with those obtained from experimental time series. The largest Lyapunov exponents for the average over all subjects and for the nearest neighbors of the model with the stimulus amplitude <italic>ζ</italic>
that fits best (<italic>ζ</italic>
= 3.6301) are plotted against the ratio of stimulus to intrinsic alpha frequency in (A). The largest Lyapunov exponents for the average over all subjects are normalized to the same range as for the model. The green area covers the standard deviation of the largest normalized Lyapunov exponents over all subjects. The comparison based on the minimization of the average of the minimum relative error (i.e., distance/maximum distance) between the normalized largest Lyapunov exponents of the model (comprising the four nearest neighbors) and of the experiment over stimulus frequencies is calculated for all stimulus amplitudes of the model. The mean error <italic>ε</italic>
of model and average over subjects is shown in (B), where the red area covers the 5% and 95% quantiles. For the significant amplitude range, the quantiles are drawn in green. Significant amplitudes are consistent over subjects as shown in <xref ref-type="table" rid="pcbi-1002298-t001">Table 1</xref>
. For more details, see Comparison section in <xref ref-type="sec" rid="s2">Materials and Methods</xref>
.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g008"></graphic>
</fig>
<fig id="pcbi-1002298-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.g009</object-id>
<label>Figure 9</label>
<caption><title>Bifurcation diagram for normalized stimulus amplitude <italic>ζ</italic>
= 3.6301.</title>
<p>The vertical axis is the normalized postsynaptic potential <italic>x</italic>
<sub>32</sub>
on pyramidal cells caused by inhibitory interneurons, that is, the coordinate of the intersection points (black dots) of trajectories with the Poincaré hyperplane after neglecting initial transients. The horizontal axis is the normalized stimulus frequency. The regimes are color-coded and indicated by the horizontal line. Periodic regimes (red) exist, for instance, for frequencies ranging from 0 to 5.34·10<sup>−2</sup>
. In this range, the system appears to undergo a period-adding bifurcation cascade by decreasing the normalized stimulus frequency. Chaotic (blue) and quasi-periodic regimes (green) occur, for example, for frequencies ranging from 5.3·10<sup>−2</sup>
to 6.23·10<sup>−2</sup>
(scattered dots) and between 17.21·10<sup>−2</sup>
and 18.72·10<sup>−2</sup>
. The classification can be also taken from <xref ref-type="table" rid="pcbi-1002298-t002">Table 2</xref>
.</p>
</caption>
<graphic xlink:href="pcbi.1002298.g009"></graphic>
</fig>
<table-wrap id="pcbi-1002298-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.t001</object-id>
<label>Table 1</label>
<caption><title>Single subject comparison of largest Lyapunov exponents.</title>
</caption>
<alternatives><graphic id="pcbi-1002298-t001-1" xlink:href="pcbi.1002298.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Participant</td>
<td align="left" rowspan="1" colspan="1">Normalized amplitude <italic>ζ</italic>
</td>
<td align="left" rowspan="1" colspan="1">Mean errors <italic>ε</italic>
in %</td>
<td align="left" rowspan="1" colspan="1">Significance (<italic>p</italic>
<.05)</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">19</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">2</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">2.89</td>
<td align="left" rowspan="1" colspan="1">30</td>
<td align="left" rowspan="1" colspan="1">no</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">4</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">17</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">16</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">3.47</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">no</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">2.34</td>
<td align="left" rowspan="1" colspan="1">26</td>
<td align="left" rowspan="1" colspan="1">no</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">23</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">9</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">20</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">10</td>
<td align="left" rowspan="1" colspan="1">3.67</td>
<td align="left" rowspan="1" colspan="1">22</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Mean</td>
<td align="left" rowspan="1" colspan="1">3.6301</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Medial</td>
<td align="left" rowspan="1" colspan="1">3.6301</td>
<td align="left" rowspan="1" colspan="1">13</td>
<td align="left" rowspan="1" colspan="1">yes</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt101"><label></label>
<p>Single subject comparisons fit the model at normalized stimulus amplitude around of approximately <italic>ζ</italic>
= 3.67 with mean errors <italic>ε</italic>
over all normalized stimulus amplitudes around 20%. The amplitudes that best fit the data are significant for the mean and median as well as for seven out of the ten subjects. The calculation is described in the Comparison section in <xref ref-type="sec" rid="s2">Materials and Methods</xref>
.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<table-wrap id="pcbi-1002298-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1002298.t002</object-id>
<label>Table 2</label>
<caption><title>Dynamic regimes occurring for normalized stimulus amplitude <italic>ζ</italic>
= 3.6301.</title>
</caption>
<alternatives><graphic id="pcbi-1002298-t002-2" xlink:href="pcbi.1002298.t002"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1">Orbit</td>
<td align="left" rowspan="1" colspan="1">Limit cycle</td>
<td align="left" rowspan="1" colspan="1">Two-torus</td>
<td align="left" rowspan="1" colspan="1">Strange attractor</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">Occurrence in frequency range (·10<sup>−2</sup>
)</td>
<td align="left" rowspan="1" colspan="1">>0 to 18.72</td>
<td align="left" rowspan="1" colspan="1">17.21 to 18.72</td>
<td align="left" rowspan="1" colspan="1">5.34 to 6.23; 7 to 7.18; 7.52 to 8.04; 8.92 to 9.29; 10.33 to 11; 11.96 to 17.21</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt102"><label></label>
<p>A limit cycle appearing as periodic oscillations is a closed orbit in state space. An invariant torus indicates quasi-periodic oscillations that manifest themselves in the power spectra with peaks at incommensurate frequencies. A strange attractor is a bounded attracting set in which all trajectories are unstable and nearby trajectories locally diverge from each other exponentially, as evidenced by a positive Lyapunov exponent, as well as a broadband power spectrum.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3"><title>Discussion</title>
<p>In this study, we analyzed the behavior of the periodically forced extended Jansen and Rit neural mass model (NMM) as a function of amplitude and frequency of the stimulus within biologically plausible ranges. The system investigated exhibits interesting and complex dynamics, including chaos. As an important result, the model was able to account for the EEG dynamics of a photic driving experiment. Photic driving paradigms are of great importance in clinical practice <xref ref-type="bibr" rid="pcbi.1002298-Jin1">[6]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Lazarev1">[7]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Takahashi1">[8]</xref>
. In this type of experiment, the dominant brain rhythm during rest, the alpha rhythm (around 10 Hz), is entrained by a periodic visual stimulus.</p>
<sec id="s3a"><title>Relation to previous results</title>
<p>It should be pointed out that many aspects of our results are in close agreement with previous studies of other types of periodically driven oscillators (see, for example, <xref ref-type="bibr" rid="pcbi.1002298-Buchli1">[75]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Pikovsky1">[76]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Osipov1">[85]</xref>
). Frequency entrainment effects have been described in, for example, the Rössler system <xref ref-type="bibr" rid="pcbi.1002298-Vadivasova1">[86]</xref>
, the Oregonator model <xref ref-type="bibr" rid="pcbi.1002298-Brons1">[87]</xref>
describing chemical oscillators such as the Belousov-Zhabotinsky reaction (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Meron1">[88]</xref>
), the Duffing oscillator describing mechanical pendulums with flexible elements <xref ref-type="bibr" rid="pcbi.1002298-Nayfeh1">[89]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Parlitz1">[90]</xref>
, the van der Pol oscillator modeling electrical triode circuits <xref ref-type="bibr" rid="pcbi.1002298-Parlitz2">[91]</xref>
, the Lorenz system describing turbulent convection in hydrodynamics <xref ref-type="bibr" rid="pcbi.1002298-Park1">[92]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Zaks1">[93]</xref>
, and the Hodgkin-Huxley model of a neuron <xref ref-type="bibr" rid="pcbi.1002298-Aihara1">[94]</xref>
. Overlapping or merging mode-locking regions in the parameter space were also discovered in a periodically driven van der Pol oscillator <xref ref-type="bibr" rid="pcbi.1002298-Flaherty1">[95]</xref>
.</p>
<p>While reverse <italic>periodic-adding</italic>
cascades appear to be the route to chaos in this study of a periodically forced Jansen and Rit model, in a number of previously investigated systems, cascades of period-doubling led to chaos, for example, in the Duffing oscillator <xref ref-type="bibr" rid="pcbi.1002298-Nayfeh1">[89]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Parlitz1">[90]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Sanchez1">[96]</xref>
, the Lorenz system <xref ref-type="bibr" rid="pcbi.1002298-Zaks1">[93]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Aizawa1">[97]</xref>
, the Rössler system <xref ref-type="bibr" rid="pcbi.1002298-Zhusubaliyev1">[98]</xref>
, the Brusselator <xref ref-type="bibr" rid="pcbi.1002298-Wang1">[99]</xref>
and the Oregonator <xref ref-type="bibr" rid="pcbi.1002298-Brons1">[87]</xref>
. In the van der Pol oscillator, both routes – period-adding <xref ref-type="bibr" rid="pcbi.1002298-Levi1">[83]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Gilbert1">[100]</xref>
and period-doubling cascades – occur <xref ref-type="bibr" rid="pcbi.1002298-Parlitz2">[91]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Mettin1">[101]</xref>
. On the other hand, our results concerning the route to chaos are in line with findings in a periodically stimulated excitable neural relaxation oscillator <xref ref-type="bibr" rid="pcbi.1002298-Coombes1">[27]</xref>
and a simple model of the Belousov-Zhabotinsky reaction <xref ref-type="bibr" rid="pcbi.1002298-Pikovsky2">[102]</xref>
. Crevier and Meister <xref ref-type="bibr" rid="pcbi.1002298-Crevier1">[103]</xref>
describe retinal (<italic>electroretinogram</italic>
, ERG) and cortical responses (LFP and <italic>visual evoked potential</italic>
, VEP) to periodic flashes of light in salamander and humans. They also found complex behavior such as frequency entrainment in experimental data as well as in their model. In contrast to our findings, they found a cascade of period-doubling bifurcations (in both data and model) that leads to chaotic regimes in their model. Finally, quasi-periodic solutions have also been reported for various systems, such as the van der Pol oscillator <xref ref-type="bibr" rid="pcbi.1002298-Parlitz2">[91]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Mettin1">[101]</xref>
, the Oregonator <xref ref-type="bibr" rid="pcbi.1002298-Brons1">[87]</xref>
, the Rössler system <xref ref-type="bibr" rid="pcbi.1002298-Zhusubaliyev1">[98]</xref>
and the Hodgkin-Huxley model <xref ref-type="bibr" rid="pcbi.1002298-Aihara1">[94]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Aihara2">[104]</xref>
.</p>
</sec>
<sec id="s3b"><title>Modeling the dynamics in photic driving data</title>
<p>We applied the concept of a periodically forced oscillator to model brain resonance effects. In the brain, such periodic input might stem from rhythmic stimulation of the brain, such as in the photic driving paradigm, or from the output of other oscillating brain areas. Such coupling between (oscillating) processes inside and outside the brain has been discussed as important for the processing of information (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Engel1">[105]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Fries1">[106]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Schnitzler1">[107]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-TallonBaudry1">[108]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Varela1">[109]</xref>
). We described the dominant intrinsic brain rhythm using the NMM performing a self-sustained oscillation, generated by an Andronov-Hopf bifurcation <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. Generally, resonance phenomena such as frequency entrainment in photic driving experiments can be explained by the concept of a periodically forced oscillator. Applying periodic input to an oscillatory system will change the current phase of the oscillation and frequency entrainment (i.e., phase locking) occurs if the sum of phase changes is nonzero over time <xref ref-type="bibr" rid="pcbi.1002298-Buchli1">[75]</xref>
(see <xref ref-type="sec" rid="s2">Results</xref>
for more details). It is expected that the dynamics of the system depend on timing, that is, the ratio between stimulus and intrinsic frequencies, as well as the intensity of stimulation. While the general trend of the frequency-detuning curves is similar for our model and experimental data, there are numerous deviations (see <xref ref-type="fig" rid="pcbi-1002298-g002">Figure 2</xref>
). These might be explained by the simplicity of the model. In the brain, many neuronal circuits are likely to be concurrently active and deviating behavior might be canceled out.</p>
<p>In our simulations, we found that the dynamics of the periodically forced extended Jansen and Rit NMM feature a rich mosaic of complex behavior (<xref ref-type="fig" rid="pcbi-1002298-g003">Figure 3</xref>
). From the parameter space analysis presented in <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
, it can be seen that both flicker intensity and frequency are critical parameters. As expected based on theory <xref ref-type="bibr" rid="pcbi.1002298-Buchli1">[75]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Pikovsky1">[76]</xref>
, the state space analysis reveals that the system is entrained by the stimulus frequency (see <xref ref-type="fig" rid="pcbi-1002298-g001">Figure 1</xref>
) where the entrainment regions (i.e., plateaus in the frequency-detuning curve) around the intrinsic frequency become wider with increasing stimulus intensity (results not shown). Also, stimulus frequencies below the intrinsic frequency lead to decelerations of the intrinsic rhythm of the modeled cortex and vice versa. This phenomenon is reflected by the ratio of stimulus to response frequencies <italic>η</italic>
/<italic>η</italic>
<sub>resp</sub>
above and below the diagonal in <xref ref-type="fig" rid="pcbi-1002298-g001">Figure 1</xref>
for stimulus frequencies below or above the intrinsic frequency (i.e., <italic>η</italic>
/<italic>η</italic>
<sub>int</sub>
<1 and <italic>η</italic>
/<italic>η</italic>
<sub>int</sub>
>1), respectively.</p>
<p>In regions of the parameter space without entrainment, complicated interaction between stimulus and intrinsic kinetics leads to periodic, quasi-periodic, and chaotic behavior, as indexed by the largest Lyapunov exponents and the Kaplan-Yorke dimension. Areas with different dynamic behavior form fractal structures in parameter space (<xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
to <xref ref-type="fig" rid="pcbi-1002298-g005"></xref>
<xref ref-type="fig" rid="pcbi-1002298-g006">Figure 6</xref>
) so that rhythmic and chaotic brain states are found virtually next to each other and even small parameter changes can give rise to a switch from one to another. For these parameter configurations, different forms of the extrinsic periodic input would affect the specific pattern of chaotic regimes in the parameter space, but not the qualitative behavior if the parameter of the stimulus shape <italic>δ</italic>
ranges between 109 and 130, as was found by additional simulations with different stimulus shape parameters (results not shown). On this account, the shape of the extrinsic input is an important model parameter for investigating the occurrence of complex regimes that needs to be investigated in the future.</p>
<p>It should, however, be pointed out that this result has been obtained from a purely deterministic model without any added or modulating noise. If noise is added to the input, this would cause jitter in its amplitude and frequency, and thereby impose a blur on the pattern depicted in <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
. However, the gross patterns are expected to survive; that is, areas with a high density of “chaotic” configurations (e.g., around 6 Hz and amplitudes between 4 and 5 mV) will feature a lower degree of predictability than areas without any configurations with positive largest Lyapunov exponents (e.g., around 11 Hz, same amplitude range). Strikingly, this is corroborated by the fact that the characteristic patterns of unpredictability generated by the model were also found with reasonable accuracy in the noisy experimental data (<xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
). We identified a particular stimulus amplitude, where, for all subjects, the Lyapunov exponents are in close agreement between experiment and model (<xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
and <xref ref-type="table" rid="pcbi-1002298-t001">Table 1</xref>
). We found that the profile of the characteristic Lyapunov spectra for the stimulus amplitude that best fits the data is preserved when noise is added to the input for a signal-to-noise-ratio (SNR) up to 10 dB (for more details, see Model in the <xref ref-type="sec" rid="s2">Materials and Methods</xref>
section). The intensity that best fits our experimental data is located in the upper portion of the effective range for exciting inhibitory interneurons. Since the largest Lyapunov exponent reflects fundamental properties of the current dynamic regime of a system (as evidenced, for example, by its sensitivity to pathological states of the brain, such as epilepsy, see <xref ref-type="bibr" rid="pcbi.1002298-Daneshyari1">[72]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-GhoshDastidar1">[73]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Raiesdana1">[74]</xref>
), the fact that our model predicts its dependence on the most important stimulus parameter (frequency) corroborates the validity of the model.</p>
<p>Consequently, we predict that a decrease in stimulus intensity in photic driving experiments would shrink and an increase would broaden the ranges of frequency entrainment (i.e., the plateaus in the frequency-detuning curve). Our model also predicts that saturation effects become important starting with approximately 1.3 times the currently applied stimulus intensity and for intensities close to zero. A stimulus increase between 1 and 1.3 times the current intensity could lead to an improved entrainment effect (i.e., broadened range). Such broadening of the entrainment range is particularly important because in clinical practice, the individual alpha frequency is usually unknown. It is important to know how great an increase in the stimulus intensity still improves the entrainment effect and hence makes the photic driving more reliable.</p>
<p>Although the effect of photic driving has long been known, and standard examination in neurology includes intermittent photic stimulation in patients with suspected photosensitive epilepsy, the exact pathomechanism is not well understood. It is known that the <italic>photoparoxysmal response</italic>
(PPR) is inheritable. In terms of electrophysiology, photosensitive epilepsy seems to be associated with changes in oscillatory activity. For example, Parra et al. <xref ref-type="bibr" rid="pcbi.1002298-Parra1">[16]</xref>
found enhanced gamma band synchrony and hypothesize that “ … some sort of recruitment or dynamic capture of neurons into larger assemblies appears to precede the epileptic chain reaction (ictal cascade) that ends in a paroxysmal oscillation, the PPR.” Likewise, Visani et al. <xref ref-type="bibr" rid="pcbi.1002298-Visani1">[110]</xref>
confirmed the potential importance of gamma band activity and found alpha band activity relevant to the PPR. Using transcranial magnetic stimulation, Siniatchkin et al. <xref ref-type="bibr" rid="pcbi.1002298-Siniatchkin1">[111]</xref>
found evidence that an increased excitability of the occipital but not the motor cortex might be associated with the PPR. The above studies indicate that a model including more than one area might be needed to further elucidate the pathomechanism of the PPR. Our model can be extended to give such experimental predictions or explanations for experimental findings. However, concrete simulations with, for instance, increased excitability in the occipital cortex and regular excitability in a second region are beyond the scope of this paper.</p>
<p>In short, we show that our model is capable of accounting for major aspects of the photic driving paradigm. This sets the scene for future work that will explore the predictions of the model in health and disease in more detail based on additional experimental data. Furthermore, a systematic exploration of the parameter space of the model with respect to brain resonance is needed. All this requires substantial efforts and is beyond the scope of the current proof-of-principle paper.</p>
<p>A principal limitation of our study is the modeling of the thalamus as independent signal generator, neglecting the cortico-thalamic feedback loop. However, we tested a model of the thalamus according to Robinson et al. <xref ref-type="bibr" rid="pcbi.1002298-Robinson1">[112]</xref>
and found that, at least for the parameters of the cortical model used in this work, the simple signal generator approach yields a good approximation. Future work will include measurements and explicit modeling of the thalamo-cortical loops.</p>
</sec>
<sec id="s3c"><title>Chaos in the brain?</title>
<p>Another issue which must be discussed is whether and to what extent our results support the idea of chaotic dynamics in the brain. The model investigated here describes complex, partially chaotic, dynamics at the mesoscopic spatial scale, which captures mass action of neural ensembles <xref ref-type="bibr" rid="pcbi.1002298-Deco1">[28]</xref>
. Chaotic dynamic regimes have been shown before in mesoscopic models of the cortex <xref ref-type="bibr" rid="pcbi.1002298-Dafilis1">[113]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Dafilis2">[114]</xref>
and of the olfactory bulb (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Freeman2">[115]</xref>
). Concerning the brain, there is evidence for chaotic behavior at different hierarchical levels, from single neurons to entire neural ensembles <xref ref-type="bibr" rid="pcbi.1002298-Freeman3">[116]</xref>
. A suitable means to experimentally access neural activity at the mesoscopic level is provided by M/EEG, which records the summed activity of 10<sup>5</sup>
to 10<sup>9</sup>
, mainly cortical, neurons <xref ref-type="bibr" rid="pcbi.1002298-Livanov1">[40]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Nunez2">[41]</xref>
. M/EEG data describe high-dimensional, noisy, nonlinear, and non-autonomous processes <xref ref-type="bibr" rid="pcbi.1002298-Freeman4">[117]</xref>
, which render it difficult to distinguish between stochastic and complex deterministic processes like deterministic chaos. Accordingly, although there is some evidence for chaos in such data (e.g., in epilepsy), the issue remains controversial (for a discussion, see <xref ref-type="bibr" rid="pcbi.1002298-Korn1">[118]</xref>
and the references cited therein). However, irrespective of whether the complexity of M/EEG fulfils the exact mathematical criteria of deterministic chaos, the parsimonious NMM, as shown here, helps to better describe the dynamics of such data and the underlying brain processes.</p>
<p>Apart from brain rhythms in characteristic frequency bands (e.g., the alpha rhythm), complex behavior with noise-like characteristics causes the continuous spectral components in these data. This can be interpreted as filtered noise (e.g., stochastic sensory input) or described by nonlinear deterministic processes. We have shown that periodically driven NMMs may explain the continuous spectral components of M/EEG without having to consider noisy input processes. Other NMM studies often apply a stochastic input process with the effect that the spectra are more realistically widened around an intrinsic frequency of interest (e.g., alpha band) (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-David1">[57]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Wendling1">[58]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Sotero1">[65]</xref>
). It is, however, an advantage, if these continuous spectral components can be modeled and controlled as intrinsic phenomena of the neural circuits, because there is evidence that broad spectral components are also modulated by cognitive processes and hence their generative processes play a role in information processing (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Makeig1">[119]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Miller1">[120]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Zhang1">[121]</xref>
). This is corroborated by the postulated prominent role of chaos in information processing (see, for example, <xref ref-type="bibr" rid="pcbi.1002298-Freeman2">[115]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Freeman3">[116]</xref>
).</p>
<p>Furthermore, the complex behavior of the NMM for certain parameter sets or ranges could be used to explain ordered sequences of dynamic regimes and multi-stability in M/EEG data by producing a temporal hierarchy <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. Such ordered sequences have been observed in, for example, perception (e.g., mono- and binocular rivalry <xref ref-type="bibr" rid="pcbi.1002298-Doesburg1">[122]</xref>
, Necker-cube illusion), stages of sleep <xref ref-type="bibr" rid="pcbi.1002298-Babloyantz1">[123]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Gallez1">[124]</xref>
, changes in attention or vigilance, learning and training such as odor recognition <xref ref-type="bibr" rid="pcbi.1002298-DiPrisco1">[125]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Freeman5">[126]</xref>
, progression of disease such as epilepsy <xref ref-type="bibr" rid="pcbi.1002298-Babloyantz2">[127]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Iasemidis1">[128]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Iasemidis2">[129]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Lehnertz1">[130]</xref>
, and effects of medication. State transitions or multi-stabilities appear since the brain is subjected to multiple high-dimensional stimuli from both exogenous (e.g., vision or haptic) and endogenous processes (e.g., endocrine or circulatory system), and is highly dependent to the current on the current individual state (e.g., vigilance, sleep or attention). For example, one can interpret the quasi-periodic behavior in <xref ref-type="fig" rid="pcbi-1002298-g003">Figure 3</xref>
as multi-stability. However, the orbits are sensitive to noise, albeit in terms of fine structure and the associated sequences, rather than the overall structure.</p>
<p>One way to achieve ordered sequences of dynamic regimes that are sufficiently robust against noise is to adequately change the state space through parameter changes that are slower than the state dynamics producing a temporal hierarchy <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. For example, one could incorporate a second model with kinetics slower than the NMM (e.g., representing metabolic processes or the neuroendocrine system) that controls a subset of the NMM parameters (e.g., couplings in terms of synaptic plasticity). In this way, the dynamic behavior of the NMM may change qualitatively through passing bifurcations and thus sequences of the complex regimes will be occur. A fine example of this approach is provided by Steyn-Ross et al. <xref ref-type="bibr" rid="pcbi.1002298-SteynRoss1">[131]</xref>
, who modeled the succession of slow wave and REM sleep phases in humans using a mean field model. The parameters of the model were controlled by the states of a low-kinetics model describing the levels of acetylcholine and somnogens (such as adenosine).</p>
<p>Note that the directions of parameter changes play an important role in parameter ranges of the system where a hysteresis occurs (see, for example, in <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
<xref ref-type="fig" rid="pcbi-1002298-g002">Figure 2</xref>
and <xref ref-type="fig" rid="pcbi-1002298-g004">Figures 4</xref>
to <xref ref-type="fig" rid="pcbi-1002298-g005"></xref>
<xref ref-type="fig" rid="pcbi-1002298-g006">6</xref>
: branch type-I A and B, and type-II AB to CC). This previous study <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
provides a catalogue of regimes that is potentially helpful to prevent the system from hysteretic behavior or, quite the reverse, to perform hysteresis.</p>
<p>The present study is the first to find complex types of behavior like entrainment, chaos, and periodic and quasi-periodic motion in a periodically forced Jansen and Rit NMM for a single cortical area for biologically plausible parameter ranges without considering noise processes. Such dynamics are observable in brain data and relevant to the explanation of brain function. We demonstrate that with the NMM, one can explain brain resonance phenomena like frequency entrainment in a clinically relevant photic driving experiment. It should be pointed out that, at this stage, the aim of our model has not been to directly improve the diagnostics of mental illnesses, but rather to allow deeper understanding of the mechanisms underlying a diagnostic tool and thereby pave the way for future new treatments and diagnosis techniques. As a logical next step, the model should be applied to pathological cases in order to specify what disease-specific inferences can be made.</p>
</sec>
<sec id="s3d"><title>Simplifying assumptions in the model</title>
<p>As any model, our model features a number of simplifications with respect to reality. The mean-field model studied embodies structural (e.g., local neural circuitry) as well as functional approximations (e.g., mean postsynaptic potential (PSP), mean firing rates and its conversions) of neural circuits to describe brain dynamics at the mesoscopic and the macroscopic levels, which are accessible, for instance, to LFP and M/EEG. Simplifications appear at all levels of modeling: the description of single cell behaviors, the modeling of neural masses (NMs) based on a single cell description (i.e., firing rate neuron), the description of the local neural circuitry and the description of networks of brain areas.</p>
<p>On the single cell level, we consider the firing rate instead of individual action potentials. Moreover, only two types of synaptic kinetics are modeled, which leads to two types of neurons that either excite or inhibit other neurons. In the brain, there is a great diversity of electrophysiological neuron types that differ in their specific input and output operations <xref ref-type="bibr" rid="pcbi.1002298-Contreras1">[43]</xref>
.</p>
<p>On the population level, the distribution of states (i.e., PSPs and firing rates) is described by their means, while variances and higher-order statistics are left out of consideration.</p>
<p>The local neural circuitry of the cortex is characterized by a wealth of distinguishable populations and their interconnections (see, for example, <xref ref-type="bibr" rid="pcbi.1002298-Thomson1">[55]</xref>
). In the model structure used here, this is simplified by simply considering pyramidal cells and two feedback loops established by inhibitory and excitatory interneurons.</p>
<p>Finally, in this work, we deal with a cortical area mean-field model, without considering projections to the rest of the brain. In particular, the thalamus, which might play a role here, is modeled only in terms of its output.</p>
<p>Moreover, since the retina is fully illuminated by the flicker that drives much of the visual cortex (see <xref ref-type="sec" rid="s2">Materials and Methods</xref>
for further explanation), we consider the entire primary visual cortex as a single source using a simple NMM for a single cortical area. Of course, the primary visual cortex is much more complicated than a single Jansen and Rit circuit, not to speak of its incorporation in brain-wide networks. Hence, our model can only represent a subset of the dynamics of the entire system. Nonetheless, our results show that the model can account for the main phenomena in the photic driving paradigm. However, one must be aware that the future availability of new or more detailed data might necessitate an extension of the model.</p>
</sec>
<sec id="s3e"><title>Link to extended brain networks</title>
<p>In this study, we show that a simple local cortical area model is already capable of performing relevant complex dynamics, particularly in response to periodic inhibitory feed-forward stimulations. Based on the fact that such a local cortical circuitry of neural populations is embedded in large-scale networks that can span the whole brain and also include subcortical structures such as the thalamus, the question arises to what extent network interactions might contribute to the complexity of brain signals such as M/EEG.</p>
<p>The present work might also contribute to the understanding of large-scale networks. In particular, the present results can be applied to inhibitory feed-forward interactions in networks between two local area models, where one model periodically performs spikes that drive the other model. The frequency entrainment or locking phenomena we found here can thus be interpreted as an effect of network interactions, which might have an impact on functional or effective brain connectivity measures such as the phase correlation (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Varela1">[109]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Spiegler2">[132]</xref>
) or the Granger causality (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Blinowska1">[133]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kaminski1">[134]</xref>
).</p>
<p>From the modeling perspective, one can obtain a dynamic regime of a local cortical area such as quasi-periodic behavior within a network by frequency locking through feed-forward inhibition from another local area by considering the following steps: (i) tuning the driving cortical area so that it performs (spiky) rhythms, (ii) selecting the dynamic regime depending on stimulus amplitude and frequency (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
to <xref ref-type="fig" rid="pcbi-1002298-g005"></xref>
<xref ref-type="fig" rid="pcbi-1002298-g006">Figure 6</xref>
), (iii) adjusting the characteristic time constant of the driver to tune the frequency of the driver to the required stimulus frequency, and (iv) selecting the characteristic potential and/or the coupling parameter so that amplitude of the driver fits. The parameters for the driving cortical area can be taken from the catalogue of regimes presented in our previous work <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
. Moreover, in order to best reproduce a specific phenomenon, for instance in M/EEG data, this catalogue helps to balance the complexity of a network, in particular, whether a single area model is sufficient or not. The effective extrinsic input ranges of a cortical area model <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
can be used to determine the coupling parameters between areas in order to prevent a network or individual cortical areas from saturating.</p>
<p>The use of these approaches to control or set up a network depends on the complexity of the graph. For instance, several bidirectional connections or feedback loops usually make a setting more difficult. In such complex graphs, one can expect more complicated behavior than for a single local area model, such as hyperchaos or phase locking of several (chaotic) regimes. However, one might to have to perform a separate analysis for the network.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="s4"><title>Materials and Methods</title>
<sec id="s4a"><title>Model</title>
<p>A generative model for brain measurements such as M/EEG can be specified by two separate systems: the state system <italic>f</italic>
explaining the usually hidden neuronal states <bold>x</bold>
(e.g., the mean postsynaptic potentials (PSPs) of neuronal populations that potentially generate M/EEG), and the observer system <italic>g</italic>
relating the neuronal states to the measurements <bold>z</bold>
:<disp-formula><graphic xlink:href="pcbi.1002298.e001"></graphic>
<label>(1)</label>
</disp-formula>
and<disp-formula><graphic xlink:href="pcbi.1002298.e002"></graphic>
<label>(2)</label>
</disp-formula>
where <bold>L</bold>
(<italic>∂</italic>
/<italic>∂t</italic>
) is a temporal differentiation operator, <bold>p</bold>
denotes the extrinsic inputs, and <bold>s<sub>x</sub>
</bold>
and <bold>s<sub>z</sub>
</bold>
parameterize state and observer system, respectively. For the state system <italic>f</italic>
, we use a neural mass model (NMM) of a cortical area. For the observer system <italic>g</italic>
, we use a simple linear relationship, as we simply consider one area (i.e., source), because the retina is fully illuminated during the photic-driving experiment that presumably drives much of primary visual cortex in parallel, and therefore no elaborate forward modeling is needed. The state system will be explained in more detail in the following paragraphs.</p>
<p>The NMM of Jansen and Rit <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen2">[52]</xref>
describes a local network representing a cortical area. This basic circuit, consisting of three interacting neural masses (NMs), namely <italic>pyramidal cells</italic>
(PCs: NM 3) with feedback loops mediated by <italic>excitatory</italic>
and <italic>inhibitory interneurons</italic>
(EINs and IINs: NMs 1 and 2), has been described in a number of previous studies (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Spruston1">[44]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen2">[52]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-LopesdaSilva2">[53]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-David1">[57]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Wendling1">[58]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Wendling2">[59]</xref>
).</p>
<p>Note that the feedback loops may also be modeled dynamically (see, for example, <xref ref-type="bibr" rid="pcbi.1002298-Nunez1">[26]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jirsa2">[48]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Robinson1">[112]</xref>
when also considering propagation delays, see, for example, <xref ref-type="bibr" rid="pcbi.1002298-Breakspear2">[135]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Marten1">[136]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Rodrigues1">[137]</xref>
). However, here we assume connections within a single area, resulting in transmission times which are shorter than the characteristic (dendritic) time constant <italic>τ</italic>
= 10 ms. Therefore, it is sufficient to describe the feedback connection by a gain constant.</p>
<p>With this NMM, the mean neuronal states can be described by a system of six nonlinearly coupled first-order ordinary differential equations:</p>
<list list-type="simple"><list-item><p>pyramidal cells (3) to (excitatory and inhibitory) interneurons (1 and 2, combined to 0)<disp-formula><graphic xlink:href="pcbi.1002298.e003"></graphic>
<label>(3)</label>
</disp-formula>
</p>
</list-item>
<list-item><p>excitatory interneurons (1) to pyramidal cells (3)<disp-formula><graphic xlink:href="pcbi.1002298.e004"></graphic>
<label>(4)</label>
</disp-formula>
</p>
</list-item>
<list-item><p>inhibitory interneurons (2) to pyramidal cells (3)<disp-formula><graphic xlink:href="pcbi.1002298.e005"></graphic>
<label>(5)</label>
</disp-formula>
</p>
</list-item>
</list>
<p>where the state vector <bold>x</bold>
= (<italic>x</italic>
<sub>03</sub>
, <italic>x</italic>
<sub>31</sub>
, <italic>x</italic>
<sub>32</sub>
, <italic>y</italic>
<sub>30</sub>
, <italic>y</italic>
<sub>31</sub>
, <italic>y</italic>
<sub>32</sub>
)<sup>T</sup>
contains the normalized mean PSPs <italic>x<sub>ba</sub>
</italic>
and currents <italic>y<sub>ba</sub>
</italic>
at NM <italic>b</italic>
caused by NM <italic>a</italic>
. The extrinsic afferents T projected to NM <italic>b</italic>
are denoted by <italic>x<sub>b</sub>
</italic>
<sub>T</sub>
. The average synaptic gains or the average numbers of synaptic contacts established between the two NMs <italic>a</italic>
and <italic>b</italic>
are represented by the constants <italic>α<sub>ba</sub>
</italic>
. Furthermore, <italic>β</italic>
is the ratio of <italic>e</italic>
xcitatory to <italic>i</italic>
nhibitory dendritic time constant <italic>β</italic>
= <italic>τ</italic>
<sub>e</sub>
/<italic>τ</italic>
<sub>i</sub>
and, in the formulas (3) to (5), the dot indicates the derivatives with respect to the normalized time <italic>κ</italic>
= <italic>t</italic>
/<italic>τ</italic>
, where <italic>τ</italic>
is the characteristic time scale. The transfer function O(<italic>x<sub>b</sub>
</italic>
) that converts the normalized mean PSP <italic>x<sub>b</sub>
</italic>
= Σ<italic><sub>a</sub>
x<sub>ba</sub>
</italic>
(i.e., the normalized potential at the axonal hillock) to the normalized mean firing rate is taken to have a sigmoidal shape O(<italic>x<sub>b</sub>
</italic>
) = 1/(1+<italic>γ</italic>
exp(−<italic>x<sub>b</sub>
</italic>
)), where <italic>γ</italic>
represents the distribution of firing thresholds within a NM. The normalized and generalized Equations (3) to (5) correspond to the Jansen and Rit model <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
with the characteristic time constant <italic>τ</italic>
= <italic>τ</italic>
<sub>e</sub>
, the coupling parameter <italic>α<sub>ba</sub>
</italic>
= 2<italic>e</italic>
<sub>0</sub>
<italic>r c<sub>ba</sub>
H</italic>
<sub>e,i</sub>
<italic>τ</italic>
<sup>2</sup>
/<italic>τ</italic>
<sub>e,i</sub>
, the sigmoid parameter <italic>γ</italic>
= exp(<italic>υ</italic>
<sub>0</sub>
<italic>r</italic>
), and the states <italic>x<sub>ba</sub>
</italic>
(<italic>κ</italic>
) = <italic>r<sub>b</sub>
υ<sub>ba</sub>
</italic>
(<italic>τκ</italic>
), with the following parameters: maximum firing rate 2<italic>e</italic>
<sub>0</sub>
, the slope of the sigmoid <italic>r</italic>
, the mean number of synaptic contacts <italic>c<sub>ba</sub>
</italic>
, and the <italic>e</italic>
xcitatory and <italic>i</italic>
nhibitory synaptic gains <italic>H</italic>
<sub>e,i</sub>
(for more details, see <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
). Note that we use normalized parameters and variables in the rest of this paper without further indicating this.</p>
<p>In this work, we explore the dynamics of the single-area model as a function of amplitude and frequency of a periodic input. This input consist of brief pulses similar to ones used by Jansen and Rit for eliciting visual evoked potentials <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Jansen2">[52]</xref>
, or used in dynamic causal modeling (e.g., <xref ref-type="bibr" rid="pcbi.1002298-David2">[63]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kiebel1">[66]</xref>
). These pulses are meant to represent the impulse response of the visual pathway, which has been investigated experimentally by a number of researchers (see <xref ref-type="bibr" rid="pcbi.1002298-Jansen2">[52]</xref>
, and the references cited therein) and described analytically by Watson and Nachmias <xref ref-type="bibr" rid="pcbi.1002298-Watson1">[138]</xref>
.</p>
<p>In the following, we will specify the parameter space to be investigated. The system described by Equations (3) to (5) has nine parameters, namely couplings <italic>α<sub>ba</sub>
</italic>
with <italic>ba</italic>
= {13, 23, 31, 32}, kinetic ratio <italic>β</italic>
, sigmoid parameter <italic>γ</italic>
, and extrinsic inputs <italic>x<sub>b</sub>
</italic>
<sub>T</sub>
with <italic>b</italic>
= {1, 2, 3}.</p>
<p>Jansen and Rit <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
proposed a specific parameter set for the NMM of a cortical area, based on a thorough discussion of the literature. The normalization of time and potentials by Jansen's excitatory dendritic time constant and sigmoid slope (<italic>τ</italic>
= 10 ms and <italic>r</italic>
= 0.56 mV<sup>−1</sup>
, respectively) leads to the following dimensionless parameters in our model: couplings <italic>α</italic>
<sub>13</sub>
= 12.285, <italic>α</italic>
<sub>23</sub>
= <italic>α</italic>
<sub>13</sub>
/4, <italic>α</italic>
<sub>31</sub>
= 4<italic>α</italic>
<sub>13</sub>
/5, <italic>α</italic>
<sub>32</sub>
= −11<italic>α</italic>
<sub>13</sub>
/13, kinetic ratio <italic>β</italic>
= 0.5 and sigmoid parameter <italic>γ</italic>
= 28.7892. The extrinsic inputs on the three NMs are taken to be constant for EINs <italic>x</italic>
<sub>1T</sub>
= 0 and PCs <italic>x</italic>
<sub>3T</sub>
= 3.36, and time-variant for IINs in the form of periodic pulses <italic>x</italic>
<sub>2T</sub>
= <italic>ζ</italic>
exp(−2<italic>δ</italic>
cos<sup>2</sup>
(<italic>θ</italic>
)), with the angle <italic>θ</italic>
<disp-formula><graphic xlink:href="pcbi.1002298.e006"></graphic>
<label>(6)</label>
</disp-formula>
specified by stimulus amplitude <italic>ζ</italic>
and stimulus frequency <italic>η</italic>
(<italic>δ</italic>
controls the shape and is set to <italic>δ</italic>
= 110). Such a peaky waveform has been found in the lateral geniculate nucleus of the thalamus in response to square visual stimuli <xref ref-type="bibr" rid="pcbi.1002298-Wilke1">[139]</xref>
. Interestingly, a very similar waveform can be generated using a NMM of the thalamus, as proposed by Robinson et al. <xref ref-type="bibr" rid="pcbi.1002298-Robinson1">[112]</xref>
, which takes into account the intra-thalamic and thalamo-cortical feedback loops (e.g.,<xref ref-type="bibr" rid="pcbi.1002298-Steriade1">[140]</xref>
). In this model, a strong inhibitory influence of the reticular nucleus on the thalamic relay cells during the relaying of external sensory stimulation, such as an on/off waveform of flickering lights, sharpens the cortical input to render it pulse-like. The time-variant input to the IINs may represent thalamic feed-forward input. This type of disynaptic <italic>feed-forward inhibition</italic>
has been described as crucial for bottom-up processing in the somatosensory (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Gabernet1">[141]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Kimura1">[142]</xref>
), auditory (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Wehr1">[143]</xref>
), and visual (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Miller2">[144]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Priebe1">[145]</xref>
) systems of rodents. Moreover, the literature provides evidence that feed-forward inhibition (e.g., from layer IV IINs driven by thalamus) dominates excitation (from thalamus) (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Miller2">[144]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Sun1">[146]</xref>
). Also, our previous model analysis of the Jansen and Rit circuit reveals the importance of input on IINs for controlling cortical behavior <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
.</p>
<p>In the absence of stimulation (i.e., <italic>x</italic>
<sub>2T</sub>
= 0), the system intrinsically performs limit cycle oscillations arising from Andronov-Hopf bifurcations, appearing as harmonic oscillations with a frequency of approximately <italic>η</italic>
<sub>intr</sub>
= 0.108 (see bifurcation diagram and phase portraits, Figure 2 and Figure 3 in <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
). Applying the characteristic dendritic time constant <italic>τ</italic>
= 10 ms as defined above, this corresponds to the parameter set proposed by Jansen and Rit <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
and an actual oscillation frequency of <italic>f</italic>
= 10.8 Hz, and can be used to describe alpha rhythms in brain signals. This characteristic time constant is used in all results reported in this work. Note, however, that varying the characteristic time constant <italic>τ</italic>
only scales the neuronal states <bold>x</bold>
(<italic>κ</italic>
) in time <italic>t</italic>
= <italic>τ κ</italic>
and thus the frequency <italic>f</italic>
= <italic>τ</italic>
<sup>−1</sup>
<italic>η</italic>
while the states <bold>x</bold>
, the form of time signals and the underlying mechanisms such as bifurcations remain unaffected. Hence, the frequency depends on the choice of the characteristic time constant <italic>τ</italic>
and thus the normalization embraces all cases of <italic>τ</italic>
. In order to study the system with periodic stimulation around the intrinsic frequency (<italic>η</italic>
<sub>intr</sub>
= 0.108), the stimulus frequency <italic>η</italic>
is taken to range from 0 to 0.19. The stimulus frequency is nonlinearly sampled ensuring (<italic>η</italic>
Δ<italic>κ</italic>
)<sup>−1</sup>
∈N<sub>1</sub>
with the sampling interval Δ<italic>κ</italic>
, so that the pulses are well sampled. The stimulus amplitude <italic>ζ</italic>
is linearly sampled from 0 to 4.1 to cover the effective range of excitatory inputs on IINs within the limit cycle which exists when extrinsic input on IINs is constant (see Figure 8 in <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
). Since the specification of the effective extrinsic input ranges is based on an analysis of the invariant transfer function (sigmoid function) <xref ref-type="bibr" rid="pcbi.1002298-Spiegler1">[62]</xref>
, these ranges are valid for any type of input, no matter whether it is constant or time-variant.</p>
<p>In summary, for analysis, we consider a system of seven first-order ordinary differential equations (Equations (3) to (6)) describing the (neuronal) states <bold>x<sup>*</sup>
</bold>
= (<italic>x</italic>
<sub>03</sub>
, <italic>x</italic>
<sub>31</sub>
, <italic>x</italic>
<sub>32</sub>
, <italic>y</italic>
<sub>30</sub>
, <italic>y</italic>
<sub>31</sub>
, <italic>y</italic>
<sub>32</sub>
, <italic>θ</italic>
)<sup>T</sup>
specified by two parameters <bold>p</bold>
= (<italic>ζ</italic>
, <italic>η</italic>
)<sup>T</sup>
.</p>
<p>We study the differential equations (3) to (6) numerically using the fourth-fifth order Runge-Kutta method over <italic>κ</italic>
= 30<bold>·</bold>
10<sup>3</sup>
in time (which equals 5 minutes for <italic>τ</italic>
<sub>JR</sub>
= 10 ms, according to Jansen and Rit <xref ref-type="bibr" rid="pcbi.1002298-Jansen1">[23]</xref>
) with a relative tolerance of 10<sup>−11</sup>
, and then linearly sampled with an interval Δ<italic>κ</italic>
= 10<sup>−2</sup>
for further analysis. From the last 6·10<sup>3</sup>
samples (last minute if <italic>τ</italic>
<sub>JR</sub>
= 10 ms), the histograms of each state were computed using the optimal number of bins <xref ref-type="bibr" rid="pcbi.1002298-Otnes1">[147]</xref>
. Using the state equations, we also compute the characteristic mean frequency of each attractor <xref ref-type="bibr" rid="pcbi.1002298-Osipov1">[85]</xref>
. Characteristic mean frequency is the time average of a trajectory over the angle velocity at points along an <italic>n</italic>
-dimensional curvature forming an attractor in state space. To study the complex behavior, we computed the power density spectra of the time series (last 6·10<sup>3</sup>
samples) using the fast Fourier transform, especially for the time series of the PSPs of the PCs, which are reflected in M/EEG.</p>
<p>We also compute the characteristic Lyapunov spectra, that is, all six Lyapunov exponents <italic>λ</italic>
<sub>1</sub>
><italic>λ</italic>
<sub>2</sub>
…><italic>λ</italic>
<sub>6</sub>
directly from the differential equations (3) to (6), using the Fortran algorithm by Chen et al. <xref ref-type="bibr" rid="pcbi.1002298-Chen1">[148]</xref>
, integrated for <italic>κ</italic>
= 1073,742.00 using a constant sample interval Δ<italic>κ</italic>
= 10<sup>−3</sup>
. Chen et al. <xref ref-type="bibr" rid="pcbi.1002298-Chen1">[148]</xref>
used a constant time-step fourth-order Adams-Bashforth integration method and a QR-reorthoginalization that also preserves orthogonality for higher-dimensional systems. The time interval is sufficiently long to stably estimate the characteristic Lyapunov spectra (error<10<sup>−6</sup>
). The Lyapunov spectrum gives a quantitative measure of the sensitivity of the states of the system dependent on the initial conditions, or, more precisely, the average rate of divergence or convergence of two neighboring trajectories in the state space. Furthermore, the whole Lyapunov spectrum enables statements of hyperchaos.</p>
<p>Hyperchaos is a higher form of chaos with at least two rather than one directions of hyperbolic instability on the attractor <xref ref-type="bibr" rid="pcbi.1002298-Rssler1">[80]</xref>
indicated by two or more positive Lyapunov exponents and by a Kaplan-Yorke dimension larger than two. Such a hyperchaotic attractor appears as a ‘folded-towel’ structure through a continuous stretching and folding in, at least, two independent directions of the state space <xref ref-type="bibr" rid="pcbi.1002298-Baier1">[149]</xref>
. Such behavior was first reported by Rössler in 1979 <xref ref-type="bibr" rid="pcbi.1002298-Rssler1">[80]</xref>
. Generally, a system that performs hyperchaos must be of at least four dimensions.</p>
<p>Due to the computational effort required, only the largest exponent or a few of the largest ones are calculated in most of the existing literature. Here, we compute the whole characteristic Lyapunov spectra running on a massive parallel cluster system of the advanced computing unit at the Computer Center, Ilmenau University of Technology. We also select several regions from the parameter space with scattered, presumably fractal, patterns of chaotic regimes (i.e., positive largest Lyapunov exponents) for recomputing at a finer stimulus amplitude and frequency resolution (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
through <xref ref-type="fig" rid="pcbi-1002298-g005"></xref>
<xref ref-type="fig" rid="pcbi-1002298-g006">Figure 6</xref>
).</p>
<p>We probe the stability of the characteristic Lyapunov spectra by adding a stochastic term to the stimulus <italic>x</italic>
<sub>2T</sub>
. Although the Gaussian noise process that we used is not autocorrelated and could lead to errors due to the constant integration step size of the Adam-Bashforth method, the estimation of the characteristic Lyapunov spectra (for the stimulus amplitude that fits the experimental data best) is stable up to a signal-to-noise-ratio (SNR) of 10 dB, especially for the 1∶1 entrainment region (i.e., <italic>η</italic>
≈<italic>η</italic>
<sub>int</sub>
). However, the stochastic term changes the characteristic Lyapunov spectra specifically, for instance, at stimulus frequencies <italic>η</italic>
around 2/3 of the intrinsic frequency <italic>η</italic>
<sub>int</sub>
. For this stimulus frequency range (0.5817<<italic>η</italic>
/<italic>η</italic>
<sub>int</sub>
<0.7632), the profile is qualitatively preserved for mild noise with a SNR up to 17 dB. We determined the SNR as the ratio of the variances of the deterministic and stochastic portions of the stimulus. The variance of the deterministic terms <italic>σ</italic>
<sup>2</sup>
(<italic>x</italic>
<sub>2T</sub>
) (i.e., periodic pulses) is given as follows<disp-formula><graphic xlink:href="pcbi.1002298.e007"></graphic>
<label>(7)</label>
</disp-formula>
where <italic>I</italic>
<sub>0</sub>
is the modified Bessel function of the first kind, <italic>δ</italic>
is the shape parameter and <italic>ζ</italic>
is the amplitude of the stimulus <italic>x</italic>
<sub>T2</sub>
.</p>
<p>The knowledge of the whole spectrum enables us to derive the Kaplan-Yorke dimension <xref ref-type="bibr" rid="pcbi.1002298-Kantz1">[150]</xref>
given by<disp-formula><graphic xlink:href="pcbi.1002298.e008"></graphic>
<label>(8)</label>
</disp-formula>
where <italic>k</italic>
is such that<disp-formula><graphic xlink:href="pcbi.1002298.e009"></graphic>
<label>(9)</label>
</disp-formula>
</p>
<p>The Kaplan-Yorke dimension measures the upper bound of the Hausdorff dimension and is similar to the information dimension (entropy) or correlation dimension of an attractor. The Hausdorff dimension quantifies the complexity of the geometry of the attractor. For example, the Hausdorff dimension of a point is zero, of a line is one, of a plane is two, but irregular sets, such as fractals or the attractors found in this work, can feature non-integer Hausdorff dimensions. We divide the state space by classifying the behavior of the system qualitatively. To that end, we specify a Poincaré map <italic>P</italic>
by choosing a suitable hyperplane transverse to the limit cycle of the unperturbed system. A Poincaré map <italic>P</italic>
considers the intersections of a trajectory existing in the <italic>d</italic>
-dimensional state space with a hyperplane of dimension <italic>d</italic>
−1. The resulting discrete series of intersection points allow the characterization of the dynamics near periodic solutions. Finally, to study the relationship between system perturbation and system response in terms of synchronization and frequency entrainment, we compute the frequency-detuning curves <xref ref-type="bibr" rid="pcbi.1002298-Pikovsky1">[76]</xref>
; that is, the difference of the response frequency (characteristic frequency or largest peak in the spectrum) and the stimulus frequency plotted against the stimulus frequency.</p>
</sec>
<sec id="s4b"><title>Experimental data</title>
<p>Experimental data were obtained by performing a photic driving experiment that was adapted to the individual alpha frequency of the subjects. Data were previously published by Schwab et al. <xref ref-type="bibr" rid="pcbi.1002298-Schwab1">[81]</xref>
. The aim of this former study was the quantification of frequency entrainment in the alpha rhythms that was most effective in the region around individual alpha and half alpha. Ten healthy participants (22 to 40 years of age, 5 male; 5 female) were stimulated by an intermittent flickering light, while EEG (32 channels, enhanced 10-20 system with a 10-10 system over the occipital region, Compumedics Neuroscan, El Paso, USA) was recorded. EEG was sampled at 1000 Hz and hardware-filtered between 0.1 Hz and 300 Hz. An initial resting condition of 60 seconds was recorded to define the individual alpha rhythm of each participant. The individual alpha frequency measured ranged from 9.5 Hz to 11.8 Hz. After this period, flicker stimulations were conducted for 15 fixed frequencies with an alpha ratio (stimulus/individual alpha frequency) ranging from 0.4 to 1.6 in each participant (randomized order of presentation). The flicker stimuli were generated by two LEDs outside the measurement chamber and were delivered via optical fibers to about 9 cm in front of the closed eyes of the subjects in order to ensure relatively stable luminance over subjects and a fully illuminated retina. The closed eyelid diffuses the flickering light from the optical fiber (with its viewing angle) so that the whole retina is illuminated (e.g., <xref ref-type="bibr" rid="pcbi.1002298-Jampolsky1">[151]</xref>
, <xref ref-type="bibr" rid="pcbi.1002298-Leijten1">[152]</xref>
). Each stimulation frequency was presented in a sequence of 20 trains. A single train contained 40 flashes and was followed by a resting period (4 s). The complete experimental design is summarized in <xref ref-type="fig" rid="pcbi-1002298-g007">Figure 7</xref>
.</p>
<p>One EEG channel located in the occipital region (O<sub>1</sub>
) was examined per participant. Data were filtered and down-sampled to 200 Hz. For each participant, periods of 62.5 s (<italic>n</italic>
= 12500 data points) were analyzed for the 15 flicker frequencies presented (the shortest available data length of the individual flicker blocks F1 through F15 is 62.5 s over all participants investigated).</p>
<p>The estimation of the largest Lyapunov exponent was based on the approach of Wolf et al. <xref ref-type="bibr" rid="pcbi.1002298-Wolf1">[153]</xref>
. An embedding dimension of 16, a time delay of 9 (≈50 ms) and an evolving time of 5 (≈25 ms) was used for the investigation of flicker stimulations. Embedding parameters were defined according to Atay and Altintas <xref ref-type="bibr" rid="pcbi.1002298-Atay3">[154]</xref>
.</p>
</sec>
<sec id="s4c"><title>Comparison</title>
<p>Our periodically driven deterministic model exhibits chaos or otherwise complex behavior in certain parameter ranges (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
). A chaotic regime can be considered as a source of noise. Empirical data as from the photic driving experiment (see Experimental data) generally represent a highly noisy (nonlinear) signal. The sources of this noise are diverse and range from technical noise (e.g., Johnson-Nyquist noise of sensors, 50/60 Hz powerline interferences), via non-brain biological noise (such as transpiration, muscle activities of the heart or eye movements), and unrelated brain activity, to the chaotic behavior of the actual stimulus processing. Thus, mapping experimental data to any biologically motivated model with reasonable accuracy is an extremely challenging task. In particular, comparing identified system variables from model and data with respect to their absolute values needs to be done delicately. Consequently, we chose to compare the Lyapunov exponents as scalar measures of the variations of the regimes of the system with respect to a well-defined “external” parameter, here the stimulus frequency (available for both model and experiment). In contrast to the experimental design, our model has the stimulus amplitude <italic>ζ</italic>
as a parameter in addition to the stimulus frequency <italic>η</italic>
. For this reason, we searched for the stimulus amplitude where the model best fits the data.</p>
<p>The calculated 15 largest Lyapunov exponents from our experimental data are all positive due to the background noise. For our model, this is the case only if chaos arises (see <xref ref-type="fig" rid="pcbi-1002298-g004">Figure 4</xref>
). The absolute values of the Lyaponov exponents can therefore not be compared directly. If, however, we assume that the unpredictability of the experimental data is partially due to background noise (which does not depend on the stimulus frequency) and partially due to the intrinsic dynamics of the modeled system, it makes sense to compare the pattern of dependency of the Lyapunov exponent from the stimulus frequency instead. For this reason, we compared the largest Lyapunov exponents as computed from the model to the largest Lyapunov exponent computed from the data, normalized to the same range as the model-based exponent (−9.6972·10<sup>−2</sup>
≤<italic>λ</italic>
<sub>1, Model</sub>
≤−9.4435·10<sup>−5</sup>
), by a shift-and-scale transformation <italic>u</italic>
+<italic>v</italic>
·<italic>λ</italic>
<sub>1</sub>
. For the means and standard deviations of <italic>u</italic>
and <italic>v</italic>
, please refer to supplementary <xref ref-type="supplementary-material" rid="pcbi.1002298.s002">Figure S2</xref>
. Interestingly, the three subjects for whom the individual fits were not significant (number 3, 6 and 7; see <xref ref-type="sec" rid="s2">Results</xref>
section) are clearly noticeable here in terms of means and standard deviations of <italic>u</italic>
and <italic>v</italic>
. The offset of the experimental Lyaponov exponents <italic>u</italic>
can be regarded as a multiplicative process <italic>R</italic>
of divergence <italic>S</italic>
<sub>1</sub>
(<italic>κ</italic>
), because |<italic>S</italic>
<sub>1</sub>
(<italic>κ</italic>
)|≈<italic>R</italic>
<sub>1</sub>
•exp(<italic>v κ·λ</italic>
<sub>1</sub>
) with <italic>R</italic>
<sub>1</sub>
= exp(<italic>u</italic>
). Background activity and more general unspecified processes may be included in <italic>R</italic>
.</p>
<p>Since the ratio between the sampling rates on the frequency axis is 4.6 between the model (69 samples) and the experimental data (15 samples), we compare an experimental data point with the four nearest neighbors in the model. The comparison and detection of the model configuration that fits the experimental data best comprise seven steps: (i) select the four nearest neighbors along the stimulus frequency axis of a amplitude configuration of our model to a query point in the experiment (i.e., the response to an experimentally applied ratio of stimulus to intrinsic alpha frequency), (ii) calculate the Euclidean distances in the plane spanned by frequency and largest normalized Lyapunov exponent between each of the nearest neighbors and the experimental data point (this way an agreement between model and data Lyaponov exponents is weighted according to the agreement between the frequencies they belong to), (iii) determine the maximum Euclidean distance between the nearest neighbors and the experimental data point, where the largest normalized Lyapunov exponent of a nearest neighbor is set to the lower (min(<italic>λ</italic>
<sub>1, Model</sub>
) = −9.6972·10<sup>−2</sup>
) or to the upper bound (max(<italic>λ</italic>
<sub>1, Model</sub>
) = −9.4435·10<sup>−5</sup>
) of the model-based largest Lyapunov exponent if the query point in the experiment is greater than −4.8533·10<sup>−2</sup>
(i.e., min(<italic>λ</italic>
<sub>1, Model</sub>
)/2+max(<italic>λ</italic>
<sub>1, Model</sub>
)/2) or not, respectively, (iv) calculate the relative error as the ratio of the distance of a nearest neighbor to its maximum distance, (v) detect the nearest neighbor with the minimum relative error for each experimental data point and average these errors over all 15 data points (for the different experimental frequencies), (vi) repeat steps (i) to (iv) for each amplitude of the model stimulation, and (vii) find the model configuration with the stimulus amplitude that fits the data best by detecting the minimum of the averaged minimum relative errors (i.e., mean error <italic>ε</italic>
in graph B of <xref ref-type="fig" rid="pcbi-1002298-g008">Figure 8</xref>
).</p>
<p>In order to test the significance of the comparison results, we (i) compute Pearson's (linear) correlation coefficient for each stimulus amplitude <italic>ζ</italic>
as well as for each subject (and also for the average over subjects) between the largest normalized Lyapunov exponents of model and data as function of the ratios of stimulus to intrinsic alpha frequency, and (ii) test its significance by applying a Student's t-test. Due to the multiple comparisons of the 106 amplitudes and each four nearest neighbors, we used the Bonferroni correction for the significance level <italic>p</italic>
= 0.05 corrected by <italic>p</italic>
′ = <italic>p</italic>
/(4<italic>·</italic>
106).</p>
<p>Moreover, in order to further substantiate our findings, we performed a bootstrap test. We randomized the sequence of experimental Lyapunov exponents along the frequency axis so that their distribution remained the same but any putative frequency-dependence was destroyed. Then we applied our fit method and recorded the fit error. We repeated this 5000 times and obtained an estimate of the error distribution. Counting the occurrences of errors that are below the one obtained with the true sequence of frequencies, we obtained an estimate of the probability that such an error could have been achieved by chance.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pcbi.1002298.s001"><label>Figure S1</label>
<caption><p><bold>Largest Lyapunov exponents of the empirical data under the presence of noise.</bold>
The exponents are characterized by mean (white lines) and standard deviation (red area) over subjects with: <bold>A</bold>
no noise and noise; <bold>B</bold>
SNR = 13 dB; <bold>C</bold>
SNR = 10 dB; <bold>D</bold>
SNR = 7 dB; <bold>E</bold>
SNR = 3 dB; and <bold>F</bold>
SNR = 0 dB.</p>
<p>(EPS)</p>
</caption>
<media xlink:href="pcbi.1002298.s001.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1002298.s002"><label>Figure S2</label>
<caption><p><bold>Means and standard deviations of shifting and scaling parameters </bold>
<bold><italic>u</italic>
</bold>
<bold> and </bold>
<bold><italic>v</italic>
</bold>
<bold> for all ten subjects.</bold>
</p>
<p>(EPS)</p>
</caption>
<media xlink:href="pcbi.1002298.s002.eps" mimetype="application" mime-subtype="postscript"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pcbi.1002298.s003"><label>Video S1</label>
<caption><p><bold>Impact of stimulus frequency for normalized stimulus amplitude </bold>
<bold><italic>ζ</italic>
</bold>
<bold> = 3.6301.</bold>
The large diagram in the middle shows the attractor as a function of the postsynaptic potentials of the three neural masses: the normalized potential that inhibitory interneurons (IINs) cause on pyramidal cells (PCs), the normalized potential that excitatory interneurons (EINs) cause on PCs, and the normalized potential that PCs cause on both interneurons (INs). The small diagrams in the upper left- and right-hand corner show the corresponding time series and the spectra (amplitude and time is normalized) respectively. The heading provides information about the current normalized stimulus frequency that increases with time. The video shows the unperturbed system for the first six seconds with the limit cycle (large diagram) producing a sinusoidal rhythm (left diagram) that appears as a narrow peak at <italic>η</italic>
= 0.108 in the spectrum (right diagram). All three curves (all in red) remain visible for the rest of the video for comparison. The three curves corresponding to a current perturbation are shown in green.</p>
<p>(AVI)</p>
</caption>
<media xlink:href="pcbi.1002298.s003.avi" mimetype="video" mime-subtype="x-msvideo"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>We would like to thank Henning Schwanbeck (Computer Center, Ilmenau University of Technology), who graciously provided timely computational support, and Stefan J. Kiebel for very helpful comments on an earlier version of this manuscript.</p>
</ack>
<fn-group><fn fn-type="conflict"><p>The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="financial-disclosure"><p>This work was supported by the Max Planck Society and the Deutsche Forschungsgemeinschaft (grant HA2899/7-1, KN588/3 and WI1166/9-1). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p>
</fn>
</fn-group>
<ref-list><title>References</title>
<ref id="pcbi.1002298-Gebber1"><label>1</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Gebber</surname>
<given-names>GL</given-names>
</name>
<name><surname>Zhong</surname>
<given-names>S</given-names>
</name>
<name><surname>Lewis</surname>
<given-names>C</given-names>
</name>
<name><surname>Barman</surname>
<given-names>SM</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Human brain alpha rhythm: nonlinear oscillation or filtered noise?</article-title>
<source>Brain Res</source>
<volume>818</volume>
<fpage>556</fpage>
<lpage>560</lpage>
<pub-id pub-id-type="pmid">10082847</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Herrmann1"><label>2</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Herrmann</surname>
<given-names>CS</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Human EEG responses to 1–100 Hz flicker: resonance phenomena in visual cortex and their potential correlation to cognitive phenomena.</article-title>
<source>Exp Brain Res</source>
<volume>137</volume>
<fpage>346</fpage>
<lpage>353</lpage>
<pub-id pub-id-type="pmid">11355381</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Silberstein1"><label>3</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Silberstein</surname>
<given-names>RB</given-names>
</name>
</person-group>
<year>1995</year>
<article-title>Steady-state visually evoked potentials, brain resonances, and cognitive processes.</article-title>
<person-group person-group-type="editor"><name><surname>Nunez</surname>
<given-names>PL</given-names>
</name>
</person-group>
<source>Neocortical dynamics and human EEG rhythms</source>
<publisher-loc>Oxford</publisher-loc>
<publisher-name>Oxford University Press</publisher-name>
<fpage>272</fpage>
<lpage>303</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Niedermeyer1"><label>4</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Niedermeyer</surname>
<given-names>E</given-names>
</name>
<name><surname>Lopes da Silva</surname>
<given-names>FH</given-names>
</name>
</person-group>
<year>1987</year>
<article-title>Electroencephalography.</article-title>
<source>Basic principles, clinical applications and related fields</source>
<publisher-loc>Baltimore-Munich</publisher-loc>
<publisher-name>Urban & Schwarzenberg</publisher-name>
</element-citation>
</ref>
<ref id="pcbi.1002298-Cobb1"><label>5</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Cobb</surname>
<given-names>S</given-names>
</name>
</person-group>
<year>1947</year>
<article-title>Photic driving as a cause of clinical seizures in epileptic patients.</article-title>
<source>Arch Neurol Psychiatry</source>
<volume>58</volume>
<fpage>70</fpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jin1"><label>6</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jin</surname>
<given-names>Y</given-names>
</name>
<name><surname>Potkin</surname>
<given-names>SG</given-names>
</name>
<name><surname>Sandman</surname>
<given-names>CA</given-names>
</name>
<name><surname>Bunney</surname>
<given-names>WE</given-names>
<suffix>Jr</suffix>
</name>
</person-group>
<year>1997</year>
<article-title>Electroencephalographic photic driving in patients with schizophrenia and depression.</article-title>
<source>Biol Psychiatry</source>
<volume>41</volume>
<fpage>496</fpage>
<lpage>499</lpage>
<pub-id pub-id-type="pmid">9034545</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Lazarev1"><label>7</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Lazarev</surname>
<given-names>VV</given-names>
</name>
<name><surname>Simpson</surname>
<given-names>DM</given-names>
</name>
<name><surname>Schubsky</surname>
<given-names>BM</given-names>
</name>
<name><surname>Deazevedo</surname>
<given-names>LC</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Photic driving in the electroencephalogram of children and adolescents: harmonic structure and relation to the resting state.</article-title>
<source>Braz J Med Biol Res</source>
<volume>34</volume>
<fpage>1573</fpage>
<lpage>1584</lpage>
<pub-id pub-id-type="pmid">11717711</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Takahashi1"><label>8</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Takahashi</surname>
<given-names>T</given-names>
</name>
<name><surname>Matsuoka</surname>
<given-names>H</given-names>
</name>
<name><surname>Okuma</surname>
<given-names>T</given-names>
</name>
</person-group>
<year>1979</year>
<article-title>High amplitude photic driving evoked by flickering-pattern in neuropsychiatric patients–with special reference to epileptics.</article-title>
<source>Folia Psychiatr Neurol Jpn</source>
<volume>33</volume>
<fpage>305</fpage>
<lpage>306</lpage>
<pub-id pub-id-type="pmid">520944</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Lazarev2"><label>9</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Lazarev</surname>
<given-names>VV</given-names>
</name>
<name><surname>Infantosi</surname>
<given-names>AF</given-names>
</name>
<name><surname>Valencio-de-Campos</surname>
<given-names>D</given-names>
</name>
<name><surname>deAzevedo</surname>
<given-names>LC</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Topographic aspects of photic driving in the electroencephalogram of children and adolescents.</article-title>
<source>Braz J Med Biol Res</source>
<volume>37</volume>
<fpage>879</fpage>
<lpage>891</lpage>
<pub-id pub-id-type="pmid">15264032</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-MirandadeSa1"><label>10</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Miranda de Sa</surname>
<given-names>AM</given-names>
</name>
<name><surname>Infantosi</surname>
<given-names>AF</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Evaluating the entrainment of the alpha rhythm during stroboscopic flash stimulation by means of coherence analysis.</article-title>
<source>Med Eng Phys</source>
<volume>27</volume>
<fpage>167</fpage>
<lpage>173</lpage>
<pub-id pub-id-type="pmid">15642512</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Fedotchev1"><label>11</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Fedotchev</surname>
<given-names>AI</given-names>
</name>
<name><surname>Bondar</surname>
<given-names>AT</given-names>
</name>
<name><surname>Konovalov</surname>
<given-names>VF</given-names>
</name>
</person-group>
<year>1990</year>
<article-title>Stability of resonance EEG reactions to flickering light in humans.</article-title>
<source>Int J Psychophysiol</source>
<volume>9</volume>
<fpage>189</fpage>
<lpage>193</lpage>
<pub-id pub-id-type="pmid">2228753</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kawaguchi1"><label>12</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kawaguchi</surname>
<given-names>T</given-names>
</name>
<name><surname>Jijiwa</surname>
<given-names>H</given-names>
</name>
<name><surname>Watanabe</surname>
<given-names>S</given-names>
</name>
</person-group>
<year>1993</year>
<article-title>The dynamics of phase relationships of alpha waves during photic driving.</article-title>
<source>Electroencephalogr Clin Neurophysiol</source>
<volume>87</volume>
<fpage>88</fpage>
<lpage>96</lpage>
<pub-id pub-id-type="pmid">7691545</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Rappelsberger1"><label>13</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Rappelsberger</surname>
<given-names>P</given-names>
</name>
<name><surname>Muller-Paschinger</surname>
<given-names>IB</given-names>
</name>
<name><surname>Petsche</surname>
<given-names>H</given-names>
</name>
<name><surname>Pockberger</surname>
<given-names>H</given-names>
</name>
<name><surname>Prohaska</surname>
<given-names>O</given-names>
</name>
<etal></etal>
</person-group>
<year>1979</year>
<article-title>About the intracortical genesis of spontaneous activity and photic driving: EEG-histological correlations in the visual cortex in rabbits EEG EMG Z.</article-title>
<source>Elektroenzephalogr Elektromyogr Verwandte Geb</source>
<volume>10</volume>
<fpage>175</fpage>
<lpage>183</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Sakamoto1"><label>14</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Sakamoto</surname>
<given-names>H</given-names>
</name>
<name><surname>Inouye</surname>
<given-names>T</given-names>
</name>
<name><surname>Shinosaki</surname>
<given-names>K</given-names>
</name>
</person-group>
<year>1993</year>
<article-title>Preservation of alpha rhythm shortly after photic driving.</article-title>
<source>Int J Neurosci</source>
<volume>73</volume>
<fpage>227</fpage>
<lpage>233</lpage>
<pub-id pub-id-type="pmid">8169058</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kalitzin1"><label>15</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kalitzin</surname>
<given-names>S</given-names>
</name>
<name><surname>Parra</surname>
<given-names>J</given-names>
</name>
<name><surname>Velis</surname>
<given-names>DN</given-names>
</name>
<name><surname>Lopes da Silva</surname>
<given-names>FH</given-names>
</name>
</person-group>
<year>2002</year>
<article-title>Enhancement of phase clustering in the EEG/MEG gamma frequency band anticipates transitions to paroxysmal epileptiform activity in epileptic patients with known visual sensitivity.</article-title>
<source>IEEE Trans Biomed Eng</source>
<volume>49</volume>
<fpage>1279</fpage>
<lpage>1286</lpage>
<pub-id pub-id-type="pmid">12450358</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Parra1"><label>16</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Parra</surname>
<given-names>J</given-names>
</name>
<name><surname>Kalitzin</surname>
<given-names>SN</given-names>
</name>
<name><surname>Iriarte</surname>
<given-names>J</given-names>
</name>
<name><surname>Blanes</surname>
<given-names>W</given-names>
</name>
<name><surname>Velis</surname>
<given-names>DN</given-names>
</name>
<etal></etal>
</person-group>
<year>2003</year>
<article-title>Gamma-band phase clustering and photosensitivity: is there an underlying mechanism common to photosensitive epilepsy and visual perception?</article-title>
<source>Brain</source>
<volume>126</volume>
<fpage>1164</fpage>
<lpage>1172</lpage>
<pub-id pub-id-type="pmid">12690055</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Lazarev3"><label>17</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Lazarev</surname>
<given-names>VV</given-names>
</name>
<name><surname>Pontes</surname>
<given-names>A</given-names>
</name>
<name><surname>Mitrofanov</surname>
<given-names>AA</given-names>
</name>
<name><surname>deAzevedo</surname>
<given-names>LC</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>Interhemispheric asymmetry in EEG photic driving coherence in childhood autism.</article-title>
<source>Clin Neurophysiol</source>
<volume>121</volume>
<fpage>145</fpage>
<lpage>152</lpage>
<pub-id pub-id-type="pmid">19951847</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-deTommaso1"><label>18</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>de Tommaso</surname>
<given-names>M</given-names>
</name>
<name><surname>Stramaglia</surname>
<given-names>S</given-names>
</name>
<name><surname>Schoffelen</surname>
<given-names>JM</given-names>
</name>
<name><surname>Guido</surname>
<given-names>M</given-names>
</name>
<name><surname>Libro</surname>
<given-names>G</given-names>
</name>
<etal></etal>
</person-group>
<year>2003</year>
<article-title>Steady-state visual evoked potentials in the low frequency range in migraine: a study of habituation and variability phenomena.</article-title>
<source>Int J Psychophysiol</source>
<volume>49</volume>
<fpage>165</fpage>
<lpage>174</lpage>
<pub-id pub-id-type="pmid">12919718</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jin2"><label>19</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jin</surname>
<given-names>Y</given-names>
</name>
<name><surname>Potkin</surname>
<given-names>SG</given-names>
</name>
<name><surname>Sandman</surname>
<given-names>CA</given-names>
</name>
<name><surname>Bunney</surname>
<given-names>WE</given-names>
</name>
</person-group>
<year>1998</year>
<article-title>Topographic analysis of EEG photic driving in patients with schizophrenia following clozapine treatment.</article-title>
<source>Clin Electroencephalogr</source>
<volume>29</volume>
<fpage>73</fpage>
<lpage>78</lpage>
<pub-id pub-id-type="pmid">9571294</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Bojak1"><label>20</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Bojak</surname>
<given-names>I</given-names>
</name>
<name><surname>Liley</surname>
<given-names>DT</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>Axonal velocity distributions in neural field equations.</article-title>
<source>PLoS Comput Biol</source>
<volume>6</volume>
<fpage>e1000653</fpage>
<pub-id pub-id-type="pmid">20126532</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Breakspear1"><label>21</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Breakspear</surname>
<given-names>M</given-names>
</name>
<name><surname>Roberts</surname>
<given-names>JA</given-names>
</name>
<name><surname>Terry</surname>
<given-names>JR</given-names>
</name>
<name><surname>Rodrigues</surname>
<given-names>S</given-names>
</name>
<name><surname>Mahant</surname>
<given-names>N</given-names>
</name>
<etal></etal>
</person-group>
<year>2006</year>
<article-title>A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis.</article-title>
<source>Cereb Cortex</source>
<volume>16</volume>
<fpage>1296</fpage>
<lpage>1313</lpage>
<pub-id pub-id-type="pmid">16280462</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Freeman1"><label>22</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Freeman</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<year>1975</year>
<source>Mass action in the nervous system : examination of the neurophysiological basis of adaptive behavior through the EEG</source>
<publisher-loc>New York</publisher-loc>
<publisher-name>Academic Press</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jansen1"><label>23</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jansen</surname>
<given-names>BH</given-names>
</name>
<name><surname>Rit</surname>
<given-names>VG</given-names>
</name>
</person-group>
<year>1995</year>
<article-title>Electroencephalogram and visual evoked potential generation in a mathematical model of coupled columns.</article-title>
<source>Biol Cybern</source>
<volume>73</volume>
<fpage>357</fpage>
<lpage>366</lpage>
<pub-id pub-id-type="pmid">7578475</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-LopesdaSilva1"><label>24</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Lopes da Silva</surname>
<given-names>FH</given-names>
</name>
<name><surname>Hoeks</surname>
<given-names>A</given-names>
</name>
<name><surname>Smits</surname>
<given-names>H</given-names>
</name>
<name><surname>Zetterberg</surname>
<given-names>LH</given-names>
</name>
</person-group>
<year>1974</year>
<article-title>Model of brain rhythmic activity. The alpha-rhythm of the thalamus.</article-title>
<source>Kybernetik</source>
<volume>15</volume>
<fpage>27</fpage>
<lpage>37</lpage>
<pub-id pub-id-type="pmid">4853232</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Markounikau1"><label>25</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Markounikau</surname>
<given-names>V</given-names>
</name>
<name><surname>Igel</surname>
<given-names>C</given-names>
</name>
<name><surname>Grinvald</surname>
<given-names>A</given-names>
</name>
<name><surname>Jancke</surname>
<given-names>D</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>A Dynamic Neural Field Model of Mesoscopic Cortical Activity Captured with Voltage-Sensitive Dye Imaging.</article-title>
<source>PLoS Comput Biol</source>
<volume>6</volume>
<fpage>e1000919</fpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Nunez1"><label>26</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Nunez</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>1974</year>
<article-title>The brain wave equation: a model for the EEG.</article-title>
<source>Math Biosci</source>
<volume>21</volume>
<fpage>279</fpage>
<lpage>297</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Coombes1"><label>27</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Coombes</surname>
<given-names>S</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>Large-scale neural dynamics: simple and complex.</article-title>
<source>Neuroimage</source>
<volume>52</volume>
<fpage>731</fpage>
<lpage>739</lpage>
<pub-id pub-id-type="pmid">20096791</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Deco1"><label>28</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>Robinson</surname>
<given-names>PA</given-names>
</name>
<name><surname>Breakspear</surname>
<given-names>M</given-names>
</name>
<name><surname>Friston</surname>
<given-names>K</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>The Dynamic Brain: From Spiking Neurons to Neural Masses and Cortical Fields.</article-title>
<source>PLoS Comput Biol</source>
<volume>4</volume>
<fpage>e1000092</fpage>
<pub-id pub-id-type="pmid">18769680</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Abbott1"><label>29</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Abbott</surname>
<given-names>LF</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Lapicque's introduction of the integrate-and-fire model neuron (1907).</article-title>
<source>Brain Res Bull</source>
<volume>50</volume>
<fpage>303</fpage>
<lpage>304</lpage>
<pub-id pub-id-type="pmid">10643408</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Hindmarsh1"><label>30</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Hindmarsh</surname>
<given-names>JL</given-names>
</name>
<name><surname>Rose</surname>
<given-names>RM</given-names>
</name>
</person-group>
<year>1984</year>
<article-title>A model of neuronal bursting using three coupled first order differential equations.</article-title>
<source>Proc R Soc London, Ser B</source>
<volume>221</volume>
<fpage>87</fpage>
<lpage>102</lpage>
<pub-id pub-id-type="pmid">6144106</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Hodgkin1"><label>31</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Hodgkin</surname>
<given-names>AL</given-names>
</name>
<name><surname>Huxley</surname>
<given-names>AF</given-names>
</name>
</person-group>
<year>1952</year>
<article-title>A Quantitative Description of Membrane Current and its Application to Conduction and Excitation in Nerve.</article-title>
<source>J Physiol</source>
<volume>117</volume>
<fpage>500</fpage>
<lpage>544</lpage>
<pub-id pub-id-type="pmid">12991237</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Finger1"><label>32</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Finger</surname>
<given-names>S</given-names>
</name>
</person-group>
<year>2001</year>
<source>Origins of Neuroscience : A History of Explorations into Brain Function</source>
<publisher-loc>Oxford; New York</publisher-loc>
<publisher-name>Oxford University Press</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Shepherd1"><label>33</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Shepherd</surname>
<given-names>GM</given-names>
</name>
</person-group>
<year>1991</year>
<source>Foundations of the neuron doctrine</source>
<publisher-loc>New York; Oxford</publisher-loc>
<publisher-name>Oxford University Press</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Basar1"><label>34</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Basar</surname>
<given-names>E</given-names>
</name>
<name><surname>Gonder</surname>
<given-names>A</given-names>
</name>
<name><surname>Ozesmi</surname>
<given-names>C</given-names>
</name>
<name><surname>Ungan</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>1975</year>
<article-title>Dynamics of brain rhythmic and evoked potentials. III. Studies in the auditory pathway, reticular formation, and hippocampus during sleep.</article-title>
<source>Biol Cybern</source>
<volume>20</volume>
<fpage>161</fpage>
<lpage>169</lpage>
<pub-id pub-id-type="pmid">1203341</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Basar2"><label>35</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Basar</surname>
<given-names>E</given-names>
</name>
<name><surname>Gonder</surname>
<given-names>A</given-names>
</name>
<name><surname>Ozesmi</surname>
<given-names>C</given-names>
</name>
<name><surname>Ungan</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>1975</year>
<article-title>Dynamics of brain rhythmic and evoked potentials. II. Studies in the auditory pathway, reticular formation, and hippocampus during the waking stage.</article-title>
<source>Biol Cybern</source>
<volume>20</volume>
<fpage>145</fpage>
<lpage>160</lpage>
<pub-id pub-id-type="pmid">1203340</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Bressler1"><label>36</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Bressler</surname>
<given-names>SL</given-names>
</name>
<name><surname>Kelso</surname>
<given-names>JAS</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Cortical coordination dynamics and cognition.</article-title>
<source>Trends Cogn Sci</source>
<volume>5</volume>
<fpage>26</fpage>
<lpage>36</lpage>
<pub-id pub-id-type="pmid">11164733</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jirsa1"><label>37</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Connectivity and dynamics of neural information processing.</article-title>
<source>Neuroinformatics</source>
<volume>2</volume>
<fpage>183</fpage>
<lpage>204</lpage>
<pub-id pub-id-type="pmid">15319516</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-McIntosh1"><label>38</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Towards a network theory of cognition.</article-title>
<source>Neural Netw</source>
<volume>13</volume>
<fpage>861</fpage>
<lpage>870</lpage>
<pub-id pub-id-type="pmid">11156197</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-John1"><label>39</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>John</surname>
<given-names>ER</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>A field theory of consciousness.</article-title>
<source>Conscious Cogn</source>
<volume>10</volume>
<fpage>184</fpage>
<lpage>213</lpage>
<pub-id pub-id-type="pmid">11414714</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Livanov1"><label>40</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Livanov</surname>
<given-names>MN</given-names>
</name>
</person-group>
<year>1977</year>
<source>Spatial organization of cerebral processes</source>
<publisher-loc>New York</publisher-loc>
<publisher-name>Wiley</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Nunez2"><label>41</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Nunez</surname>
<given-names>PL</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Toward a quantitative description of large-scale neocortical dynamic function and EEG.</article-title>
<source>Behav Brain Sci</source>
<volume>23</volume>
<fpage>371</fpage>
<lpage>398; discussion 399–437</lpage>
<pub-id pub-id-type="pmid">11301576</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Waters1"><label>42</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Waters</surname>
<given-names>J</given-names>
</name>
<name><surname>Schaefer</surname>
<given-names>A</given-names>
</name>
<name><surname>Sakmann</surname>
<given-names>B</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Backpropagating action potentials in neurones: measurement, mechanisms and potential functions.</article-title>
<source>Prog Biophys Mol Biol</source>
<volume>87</volume>
<fpage>145</fpage>
<lpage>170</lpage>
<pub-id pub-id-type="pmid">15471594</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Contreras1"><label>43</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Contreras</surname>
<given-names>D</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Electrophysiological classes of neocortical neurons.</article-title>
<source>Neural Netw</source>
<volume>17</volume>
<fpage>633</fpage>
<lpage>646</lpage>
<pub-id pub-id-type="pmid">15288889</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Spruston1"><label>44</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Spruston</surname>
<given-names>N</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Pyramidal neurons: dendritic structure and synaptic integration.</article-title>
<source>Nat Rev Neurosci</source>
<volume>9</volume>
<fpage>206</fpage>
<lpage>221</lpage>
<pub-id pub-id-type="pmid">18270515</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Frank1"><label>45</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Frank</surname>
<given-names>TD</given-names>
</name>
<name><surname>Daffertshofer</surname>
<given-names>A</given-names>
</name>
<name><surname>Beek</surname>
<given-names>PJ</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Multivariate Ornstein-Uhlenbeck processes with mean-field dependent coefficients: application to postural sway.</article-title>
<source>Phys Rev E Stat Nonlin Soft Matter Phys</source>
<volume>63</volume>
<fpage>011905</fpage>
<pub-id pub-id-type="pmid">11304285</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Harrison1"><label>46</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Harrison</surname>
<given-names>LM</given-names>
</name>
<name><surname>David</surname>
<given-names>O</given-names>
</name>
<name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Stochastic models of neuronal dynamics.</article-title>
<source>Philos Trans R Soc Lond B Biol Sci</source>
<volume>360</volume>
<fpage>1075</fpage>
<lpage>1091</lpage>
<pub-id pub-id-type="pmid">16087449</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Knight1"><label>47</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Knight</surname>
<given-names>BW</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Dynamics of encoding in neuron populations: some general mathematical features.</article-title>
<source>Neural Comput</source>
<volume>12</volume>
<fpage>473</fpage>
<lpage>518</lpage>
<pub-id pub-id-type="pmid">10769319</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jirsa2"><label>48</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>Haken</surname>
<given-names>H</given-names>
</name>
</person-group>
<year>1996</year>
<article-title>Field theory of electromagnetic brain activity.</article-title>
<source>Phys Rev Lett</source>
<volume>77</volume>
<fpage>960</fpage>
<lpage>963</lpage>
<pub-id pub-id-type="pmid">10062950</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wilson1"><label>49</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wilson</surname>
<given-names>HR</given-names>
</name>
<name><surname>Cowan</surname>
<given-names>JD</given-names>
</name>
</person-group>
<year>1973</year>
<article-title>A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue.</article-title>
<source>Kybernetik</source>
<volume>13</volume>
<fpage>55</fpage>
<lpage>80</lpage>
<pub-id pub-id-type="pmid">4767470</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Atay1"><label>50</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Atay</surname>
<given-names>FM</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Stability and bifurcations in neural fields with finite propagation speed and general connectivity.</article-title>
<source>SIAM J Appl Math</source>
<volume>65</volume>
<fpage>644</fpage>
<lpage>666</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Atay2"><label>51</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Atay</surname>
<given-names>FM</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Neural fields with distributed transmission speeds and long-range feedback delays.</article-title>
<source>SIAM J Appl Dynam Syst</source>
<volume>5</volume>
<fpage>670</fpage>
<lpage>698</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jansen2"><label>52</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jansen</surname>
<given-names>BH</given-names>
</name>
<name><surname>Zouridakis</surname>
<given-names>G</given-names>
</name>
<name><surname>Brandt</surname>
<given-names>ME</given-names>
</name>
</person-group>
<year>1993</year>
<article-title>A neurophysiologically-based mathematical model of flash visual evoked potentials.</article-title>
<source>Biol Cybern</source>
<volume>68</volume>
<fpage>275</fpage>
<lpage>283</lpage>
<pub-id pub-id-type="pmid">8452897</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-LopesdaSilva2"><label>53</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Lopes da Silva</surname>
<given-names>FH</given-names>
</name>
<name><surname>van Rotterdam</surname>
<given-names>A</given-names>
</name>
<name><surname>Barts</surname>
<given-names>P</given-names>
</name>
<name><surname>van Heusden</surname>
<given-names>E</given-names>
</name>
<name><surname>Burr</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>1976</year>
<article-title>Models of neuronal populations: the basic mechanisms of rhythmicity.</article-title>
<source>Prog Brain Res</source>
<volume>45</volume>
<fpage>281</fpage>
<lpage>308</lpage>
<pub-id pub-id-type="pmid">1013341</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Zetterberg1"><label>54</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Zetterberg</surname>
<given-names>LH</given-names>
</name>
<name><surname>Kristiansson</surname>
<given-names>L</given-names>
</name>
<name><surname>Mossberg</surname>
<given-names>K</given-names>
</name>
</person-group>
<year>1978</year>
<article-title>Performance of a model for a local neuron population.</article-title>
<source>Biol Cybern</source>
<volume>31</volume>
<fpage>15</fpage>
<lpage>26</lpage>
<pub-id pub-id-type="pmid">728488</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Thomson1"><label>55</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Thomson</surname>
<given-names>AM</given-names>
</name>
<name><surname>Bannister</surname>
<given-names>AP</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Interlaminar connections in the neocortex.</article-title>
<source>Cereb Cortex</source>
<volume>13</volume>
<fpage>5</fpage>
<lpage>14</lpage>
<pub-id pub-id-type="pmid">12466210</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Haeusler1"><label>56</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Haeusler</surname>
<given-names>S</given-names>
</name>
<name><surname>Maass</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>2007</year>
<article-title>A statistical analysis of information-processing properties of lamina-specific cortical microcircuit models.</article-title>
<source>Cereb Cortex</source>
<volume>17</volume>
<fpage>149</fpage>
<lpage>162</lpage>
<pub-id pub-id-type="pmid">16481565</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-David1"><label>57</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>David</surname>
<given-names>O</given-names>
</name>
<name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>A neural mass model for MEG/EEG: coupling and neuronal dynamics.</article-title>
<source>Neuroimage</source>
<volume>20</volume>
<fpage>1743</fpage>
<lpage>1755</lpage>
<pub-id pub-id-type="pmid">14642484</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wendling1"><label>58</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wendling</surname>
<given-names>F</given-names>
</name>
<name><surname>Bartolomei</surname>
<given-names>F</given-names>
</name>
<name><surname>Bellanger</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Chauvel</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>2002</year>
<article-title>Epileptic fast activity can be explained by a model of impaired GABAergic dendritic inhibition.</article-title>
<source>Eur J Neurosci</source>
<volume>15</volume>
<fpage>1499</fpage>
<lpage>1508</lpage>
<pub-id pub-id-type="pmid">12028360</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wendling2"><label>59</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wendling</surname>
<given-names>F</given-names>
</name>
<name><surname>Bellanger</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Bartolomei</surname>
<given-names>F</given-names>
</name>
<name><surname>Chauvel</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals.</article-title>
<source>Biol Cybern</source>
<volume>83</volume>
<fpage>367</fpage>
<lpage>378</lpage>
<pub-id pub-id-type="pmid">11039701</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Grimbert1"><label>60</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Grimbert</surname>
<given-names>F</given-names>
</name>
<name><surname>Faugeras</surname>
<given-names>O</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Bifurcation analysis of Jansen's neural mass model.</article-title>
<source>Neural Comput</source>
<volume>18</volume>
<fpage>3052</fpage>
<lpage>3068</lpage>
<pub-id pub-id-type="pmid">17052158</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Touboul1"><label>61</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Touboul</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>2009</year>
<source>Nonlinear and Stochastic Models in Neurosciences [PhD Thesis]</source>
<publisher-loc>Nice - Sophia Antipolis</publisher-loc>
<publisher-name>École Polytechnique</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Spiegler1"><label>62</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Spiegler</surname>
<given-names>A</given-names>
</name>
<name><surname>Kiebel</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Atay</surname>
<given-names>FM</given-names>
</name>
<name><surname>Knösche</surname>
<given-names>TR</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>Bifurcation analysis of neural mass models: Impact of extrinsic inputs and dendritic time constants.</article-title>
<source>Neuroimage</source>
<volume>52</volume>
<fpage>1041</fpage>
<lpage>1058</lpage>
<pub-id pub-id-type="pmid">20045068</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-David2"><label>63</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>David</surname>
<given-names>O</given-names>
</name>
<name><surname>Harrison</surname>
<given-names>L</given-names>
</name>
<name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Modelling event-related responses in the brain.</article-title>
<source>Neuroimage</source>
<volume>25</volume>
<fpage>756</fpage>
<lpage>770</lpage>
<pub-id pub-id-type="pmid">15808977</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-David3"><label>64</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>David</surname>
<given-names>O</given-names>
</name>
<name><surname>Kiebel</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Harrison</surname>
<given-names>LM</given-names>
</name>
<name><surname>Mattout</surname>
<given-names>J</given-names>
</name>
<name><surname>Kilner</surname>
<given-names>JM</given-names>
</name>
<etal></etal>
</person-group>
<year>2006</year>
<article-title>Dynamic causal modeling of evoked responses in EEG and MEG.</article-title>
<source>Neuroimage</source>
<volume>30</volume>
<fpage>1255</fpage>
<lpage>1272</lpage>
<pub-id pub-id-type="pmid">16473023</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Sotero1"><label>65</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Sotero</surname>
<given-names>RC</given-names>
</name>
<name><surname>Trujillo-Barreto</surname>
<given-names>NJ</given-names>
</name>
<name><surname>Iturria-Medina</surname>
<given-names>Y</given-names>
</name>
<name><surname>Carbonell</surname>
<given-names>F</given-names>
</name>
<name><surname>Jimenez</surname>
<given-names>JC</given-names>
</name>
</person-group>
<year>2007</year>
<article-title>Realistically coupled neural mass models can generate EEG rhythms.</article-title>
<source>Neural Comput</source>
<volume>19</volume>
<fpage>478</fpage>
<lpage>512</lpage>
<pub-id pub-id-type="pmid">17206872</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kiebel1"><label>66</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kiebel</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Garrido</surname>
<given-names>MI</given-names>
</name>
<name><surname>Moran</surname>
<given-names>R</given-names>
</name>
<name><surname>Chen</surname>
<given-names>CC</given-names>
</name>
<name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Dynamic Causal Modeling for EEG and MEG.</article-title>
<source>Hum Brain Mapp</source>
<volume>30</volume>
<fpage>1866</fpage>
<lpage>1876</lpage>
<pub-id pub-id-type="pmid">19360734</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kiebel2"><label>67</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kiebel</surname>
<given-names>SJ</given-names>
</name>
<name><surname>David</surname>
<given-names>O</given-names>
</name>
<name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Dynamic causal modelling of evoked responses in EEG/MEG with lead field parameterization.</article-title>
<source>Neuroimage</source>
<volume>30</volume>
<fpage>1273</fpage>
<lpage>1284</lpage>
<pub-id pub-id-type="pmid">16490364</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kiebel3"><label>68</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kiebel</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Garrido</surname>
<given-names>MI</given-names>
</name>
<name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<year>2007</year>
<article-title>Dynamic causal modelling of evoked responses: The role of intrinsic connections.</article-title>
<source>Neuroimage</source>
<volume>36</volume>
<fpage>332</fpage>
<lpage>345</lpage>
<pub-id pub-id-type="pmid">17462916</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Moran1"><label>69</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Moran</surname>
<given-names>RJ</given-names>
</name>
<name><surname>Stephan</surname>
<given-names>KE</given-names>
</name>
<name><surname>Seidenbecher</surname>
<given-names>T</given-names>
</name>
<name><surname>Pape</surname>
<given-names>HC</given-names>
</name>
<name><surname>Dolan</surname>
<given-names>RJ</given-names>
</name>
<etal></etal>
</person-group>
<year>2009</year>
<article-title>Dynamic causal models of steady-state responses.</article-title>
<source>Neuroimage</source>
<volume>44</volume>
<fpage>796</fpage>
<lpage>811</lpage>
<pub-id pub-id-type="pmid">19000769</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Buzski1"><label>70</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Buzsáki</surname>
<given-names>G</given-names>
</name>
</person-group>
<year>2006</year>
<source>Rhythms of the brain</source>
<publisher-loc>Oxford; New York</publisher-loc>
<publisher-name>Oxford University Press</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Buzski2"><label>71</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Buzsáki</surname>
<given-names>G</given-names>
</name>
<name><surname>Draguhn</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Neuronal Oscillations in Cortical Networks.</article-title>
<source>Science</source>
<volume>304</volume>
<fpage>1926</fpage>
<lpage>1929</lpage>
<pub-id pub-id-type="pmid">15218136</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Daneshyari1"><label>72</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Daneshyari</surname>
<given-names>M</given-names>
</name>
<name><surname>Kamkar</surname>
<given-names>LL</given-names>
</name>
</person-group>
<year>2010</year>
<source>Epileptic EEG: A Comprehensive Study of Nonlinear Behavior. Advances in Computational Biology</source>
<publisher-loc>Berlin</publisher-loc>
<publisher-name>Springer-Verlag Berlin</publisher-name>
<fpage>677</fpage>
<lpage>683</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-GhoshDastidar1"><label>73</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Ghosh-Dastidar</surname>
<given-names>S</given-names>
</name>
<name><surname>Adeli</surname>
<given-names>H</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Chaos in the Brain: Novel Methodologies for Epilepsy Diagnosis and Seizure Detection.</article-title>
<person-group person-group-type="editor"><name><surname>Skiadas</surname>
<given-names>CH</given-names>
</name>
<name><surname>Dimotikalis</surname>
<given-names>I</given-names>
</name>
<name><surname>Skiadas</surname>
<given-names>C</given-names>
</name>
</person-group>
<publisher-loc>Singapore</publisher-loc>
<publisher-name>World Scientific Publ Co Pte Ltd</publisher-name>
<fpage>138</fpage>
<lpage>148</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Raiesdana1"><label>74</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Raiesdana</surname>
<given-names>S</given-names>
</name>
<name><surname>Golpayegani</surname>
<given-names>S</given-names>
</name>
<name><surname>Firoozabadi</surname>
<given-names>SMP</given-names>
</name>
<name><surname>Habibabadi</surname>
<given-names>JM</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>On the discrimination of patho-physiological states in epilepsy by means of dynamical measures.</article-title>
<source>Comput Biol Med</source>
<volume>39</volume>
<fpage>1073</fpage>
<lpage>1082</lpage>
<pub-id pub-id-type="pmid">19880102</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Buchli1"><label>75</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Buchli</surname>
<given-names>J</given-names>
</name>
<name><surname>Righetti</surname>
<given-names>L</given-names>
</name>
<name><surname>Ijspeert</surname>
<given-names>AJ</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Engineering entrainment and adaptation in limit cycle systems - From biological inspiration to applications in robotics.</article-title>
<source>Biol Cybern</source>
<volume>95</volume>
<fpage>645</fpage>
<lpage>664</lpage>
<pub-id pub-id-type="pmid">17146662</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Pikovsky1"><label>76</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Pikovsky</surname>
<given-names>A</given-names>
</name>
<name><surname>Rosenblum</surname>
<given-names>M</given-names>
</name>
<name><surname>Kurths</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>2003</year>
<source>Synchronization: a universal concept in nonlinear sciences</source>
<publisher-loc>Cambridge</publisher-loc>
<publisher-name>Cambridge University Press</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Shilnikov1"><label>77</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Shil'nikov</surname>
<given-names>LP</given-names>
</name>
</person-group>
<year>1969</year>
<article-title>On a new type of bifurcation in multidimensional dynamical systems.</article-title>
<source>Sov Math Dokl</source>
<volume>10</volume>
<fpage>1368</fpage>
<lpage>1371</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Beer1"><label>78</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Beer</surname>
<given-names>RD</given-names>
</name>
</person-group>
<year>1995</year>
<article-title>On the Dynamics of Small Continuous-Time Recurrent Neural Networks.</article-title>
<source>Adapt Behav</source>
<volume>3</volume>
<fpage>469</fpage>
<lpage>509</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-vanVeen1"><label>79</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>van Veen</surname>
<given-names>L</given-names>
</name>
<name><surname>Liley</surname>
<given-names>DT</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Chaos via Shilnikov's saddle-node bifurcation in a theory of the electroencephalogram.</article-title>
<source>Phys Rev Lett</source>
<volume>97</volume>
<fpage>208101</fpage>
<pub-id pub-id-type="pmid">17155719</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Rssler1"><label>80</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Rössler</surname>
<given-names>OE</given-names>
</name>
</person-group>
<year>1979</year>
<article-title>An Equation for Hyperchaos.</article-title>
<source>Phys Lett A</source>
<volume>71</volume>
<fpage>155</fpage>
<lpage>157</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Schwab1"><label>81</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Schwab</surname>
<given-names>K</given-names>
</name>
<name><surname>Ligges</surname>
<given-names>C</given-names>
</name>
<name><surname>Jungmann</surname>
<given-names>T</given-names>
</name>
<name><surname>Hilgenfeld</surname>
<given-names>B</given-names>
</name>
<name><surname>Haueisen</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<year>2006</year>
<article-title>Alpha entrainment in human electroencephalogram and magnetoencephalogram recordings.</article-title>
<source>Neuroreport</source>
<volume>17</volume>
<fpage>1829</fpage>
<lpage>1833</lpage>
<pub-id pub-id-type="pmid">17164673</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Coombes2"><label>82</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Coombes</surname>
<given-names>S</given-names>
</name>
<name><surname>Osbaldestin</surname>
<given-names>AH</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Period-adding bifurcations and chaos in a periodically stimulated excitable neural relaxation oscillator.</article-title>
<source>Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics</source>
<volume>62</volume>
<fpage>4057</fpage>
<lpage>4066</lpage>
<pub-id pub-id-type="pmid">11088930</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Levi1"><label>83</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Levi</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1990</year>
<article-title>A Period-Adding Phenomenon.</article-title>
<source>SIAM J Appl Math</source>
<volume>50</volume>
<fpage>943</fpage>
<lpage>955</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-LoFaro1"><label>84</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>LoFaro</surname>
<given-names>T</given-names>
</name>
</person-group>
<year>1996</year>
<article-title>Period-adding bifurcations in a one parameter family of interval maps.</article-title>
<source>Math Comput Model</source>
<volume>24</volume>
<fpage>27</fpage>
<lpage>41</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Osipov1"><label>85</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Osipov</surname>
<given-names>GV</given-names>
</name>
<name><surname>Kurths</surname>
<given-names>J</given-names>
</name>
<name><surname>Zhou</surname>
<given-names>C</given-names>
</name>
</person-group>
<year>2007</year>
<source>Synchronization in oscillatory networks</source>
<publisher-loc>Berlin; New York</publisher-loc>
<publisher-name>Springer</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Vadivasova1"><label>86</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Vadivasova</surname>
<given-names>TE</given-names>
</name>
<name><surname>Balanov</surname>
<given-names>AG</given-names>
</name>
<name><surname>Sosnovtseva</surname>
<given-names>OV</given-names>
</name>
<name><surname>Postnov</surname>
<given-names>DE</given-names>
</name>
<name><surname>Mosekilde</surname>
<given-names>E</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Synchronization in driven chaotic systems: Diagnostics and bifurcations.</article-title>
<source>Phys Lett A</source>
<volume>253</volume>
<fpage>66</fpage>
<lpage>74</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Brons1"><label>87</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Brons</surname>
<given-names>M</given-names>
</name>
<name><surname>Gross</surname>
<given-names>P</given-names>
</name>
<name><surname>Bar-Eli</surname>
<given-names>K</given-names>
</name>
</person-group>
<year>1997</year>
<article-title>Circle maps and the Devil's staircase in a periodically perturbed Oregonator.</article-title>
<source>Int J Bifurcat Chaos</source>
<volume>7</volume>
<fpage>2621</fpage>
<lpage>2628</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Meron1"><label>88</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Meron</surname>
<given-names>E</given-names>
</name>
</person-group>
<year>1992</year>
<article-title>Pattern-Formation in Excitable Media.</article-title>
<source>Phys Rep</source>
<volume>218</volume>
<fpage>1</fpage>
<lpage>66</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Nayfeh1"><label>89</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Nayfeh</surname>
<given-names>AH</given-names>
</name>
<name><surname>Sanchez</surname>
<given-names>NE</given-names>
</name>
</person-group>
<year>1989</year>
<article-title>Bifurcations in a Forced Softening Duffing Oscillator.</article-title>
<source>Int J Nonlinear Mech</source>
<volume>24</volume>
<fpage>483</fpage>
<lpage>497</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Parlitz1"><label>90</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Parlitz</surname>
<given-names>U</given-names>
</name>
<name><surname>Lauterborn</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>1985</year>
<article-title>Superstructure in the Bifurcation Set of the Duffing Equation X+Dx+X+X3 = Fcos(W)T).</article-title>
<source>Phys Lett A</source>
<volume>107</volume>
<fpage>351</fpage>
<lpage>355</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Parlitz2"><label>91</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Parlitz</surname>
<given-names>U</given-names>
</name>
<name><surname>Lauterborn</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>1987</year>
<article-title>Period-Doubling Cascades and Devils Staircases of the Driven Van Der Pol Oscillator.</article-title>
<source>Phys Rev A</source>
<volume>36</volume>
<fpage>1428</fpage>
<lpage>1434</lpage>
<pub-id pub-id-type="pmid">9899002</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Park1"><label>92</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Park</surname>
<given-names>EH</given-names>
</name>
<name><surname>Zaks</surname>
<given-names>MA</given-names>
</name>
<name><surname>Kurths</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Phase synchronization in the forced Lorenz system.</article-title>
<source>Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics</source>
<volume>60</volume>
<fpage>6627</fpage>
<lpage>6638</lpage>
<pub-id pub-id-type="pmid">11970583</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Zaks1"><label>93</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Zaks</surname>
<given-names>MA</given-names>
</name>
<name><surname>Park</surname>
<given-names>EH</given-names>
</name>
<name><surname>Rosenblum</surname>
<given-names>MG</given-names>
</name>
<name><surname>Kurths</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Alternating locking ratios in imperfect phase synchronization.</article-title>
<source>Phys Rev Lett</source>
<volume>82</volume>
<fpage>4228</fpage>
<lpage>4231</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Aihara1"><label>94</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Aihara</surname>
<given-names>K</given-names>
</name>
<name><surname>Matsumoto</surname>
<given-names>G</given-names>
</name>
<name><surname>Ikegaya</surname>
<given-names>Y</given-names>
</name>
</person-group>
<year>1984</year>
<article-title>Periodic and Non-Periodic Responses of a Periodically Forced Hodgkin-Huxley Oscillator.</article-title>
<source>J Theor Biol</source>
<volume>109</volume>
<fpage>249</fpage>
<lpage>269</lpage>
<pub-id pub-id-type="pmid">6482467</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Flaherty1"><label>95</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Flaherty</surname>
<given-names>JE</given-names>
</name>
<name><surname>Hoppensteadt</surname>
<given-names>FC</given-names>
</name>
</person-group>
<year>1978</year>
<article-title>Frequency Entrainment of a Forced Van Der Pol Oscillator.</article-title>
<source>Studies in Appl Math</source>
<volume>58</volume>
<fpage>5</fpage>
<lpage>15</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Sanchez1"><label>96</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Sanchez</surname>
<given-names>NE</given-names>
</name>
<name><surname>Nayfeh</surname>
<given-names>AH</given-names>
</name>
</person-group>
<year>1990</year>
<article-title>Prediction of Bifurcations in a Parametrically Excited Duffing Oscillator.</article-title>
<source>Int J Nonlinear Mech</source>
<volume>25</volume>
<fpage>163</fpage>
<lpage>176</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Aizawa1"><label>97</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Aizawa</surname>
<given-names>Y</given-names>
</name>
<name><surname>Uezu</surname>
<given-names>T</given-names>
</name>
</person-group>
<year>1982</year>
<article-title>Topological Aspects in Chaos and in 2k-Period Doubling Cascade.</article-title>
<source>Prog Theor Phys</source>
<volume>67</volume>
<fpage>982</fpage>
<lpage>985</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Zhusubaliyev1"><label>98</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Zhusubaliyev</surname>
<given-names>ZT</given-names>
</name>
<name><surname>Laugesen</surname>
<given-names>JL</given-names>
</name>
<name><surname>Mosekilde</surname>
<given-names>E</given-names>
</name>
</person-group>
<year>2010</year>
<article-title>From multi-layered resonance tori to period-doubled ergodic tori.</article-title>
<source>Phys Lett A</source>
<volume>374</volume>
<fpage>2534</fpage>
<lpage>2538</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wang1"><label>99</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>G-R</given-names>
</name>
<name><surname>Chen</surname>
<given-names>S-G</given-names>
</name>
<name><surname>Hao</surname>
<given-names>B-L</given-names>
</name>
</person-group>
<year>1984</year>
<article-title>Intermittent chaos in the forced Brusselator.</article-title>
<source>Chinese Phys</source>
<volume>4</volume>
<fpage>284</fpage>
<lpage>293</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Gilbert1"><label>100</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Gilbert</surname>
<given-names>T</given-names>
</name>
<name><surname>Gammon</surname>
<given-names>RW</given-names>
</name>
</person-group>
<year>2000</year>
<article-title>Stable oscillations and Devil's staircase in the Van der Pol oscillator.</article-title>
<source>Int J Bifurcat Chaos</source>
<volume>10</volume>
<fpage>155</fpage>
<lpage>164</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Mettin1"><label>101</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Mettin</surname>
<given-names>R</given-names>
</name>
<name><surname>Parlitz</surname>
<given-names>U</given-names>
</name>
<name><surname>Lauterborn</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>1993</year>
<article-title>Bifurcation structure of the driven van der Pol oscillator.</article-title>
<source>Int J Bifurcat Chaos</source>
<volume>3</volume>
<fpage>1529</fpage>
<lpage>1555</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Pikovsky2"><label>102</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Pikovsky</surname>
<given-names>AS</given-names>
</name>
</person-group>
<year>1981</year>
<article-title>A Dynamical Model for Periodic and Chaotic Oscillations in the Belousov-Zhabotinsky Reaction.</article-title>
<source>Phys Lett A</source>
<volume>85</volume>
<fpage>13</fpage>
<lpage>16</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Crevier1"><label>103</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Crevier</surname>
<given-names>DW</given-names>
</name>
<name><surname>Meister</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1998</year>
<article-title>Synchronous period-doubling in flicker vision of salamander and man.</article-title>
<source>J Neurophysiol</source>
<volume>79</volume>
<fpage>1869</fpage>
<lpage>1878</lpage>
<pub-id pub-id-type="pmid">9535954</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Aihara2"><label>104</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Aihara</surname>
<given-names>K</given-names>
</name>
<name><surname>Numajiri</surname>
<given-names>T</given-names>
</name>
<name><surname>Matsumoto</surname>
<given-names>G</given-names>
</name>
<name><surname>Kotani</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1986</year>
<article-title>Structures of Attractors in Periodically Forced Neural Oscillators.</article-title>
<source>Phys Lett A</source>
<volume>116</volume>
<fpage>313</fpage>
<lpage>317</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Engel1"><label>105</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Engel</surname>
<given-names>AK</given-names>
</name>
<name><surname>Singer</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Temporal binding and the neural correlates of sensory awareness.</article-title>
<source>Trends Cogn Sci</source>
<volume>5</volume>
<fpage>16</fpage>
<lpage>25</lpage>
<pub-id pub-id-type="pmid">11164732</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Fries1"><label>106</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Fries</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>A mechanism for cognitive dynamics: neuronal communication through neuronal coherence.</article-title>
<source>Trends Cogn Sci</source>
<volume>9</volume>
<fpage>474</fpage>
<lpage>480</lpage>
<pub-id pub-id-type="pmid">16150631</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Schnitzler1"><label>107</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Schnitzler</surname>
<given-names>A</given-names>
</name>
<name><surname>Gross</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Normal and pathological oscillatory communication in the brain.</article-title>
<source>Nat Rev Neurosci</source>
<volume>6</volume>
<fpage>285</fpage>
<lpage>296</lpage>
<pub-id pub-id-type="pmid">15803160</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-TallonBaudry1"><label>108</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Tallon-Baudry</surname>
<given-names>C</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Attention and awareness in synchrony.</article-title>
<source>Trends Cogn Sci</source>
<volume>8</volume>
<fpage>523</fpage>
<lpage>525</lpage>
<pub-id pub-id-type="pmid">15556020</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Varela1"><label>109</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Varela</surname>
<given-names>F</given-names>
</name>
<name><surname>Lachaux</surname>
<given-names>JP</given-names>
</name>
<name><surname>Rodriguez</surname>
<given-names>E</given-names>
</name>
<name><surname>Martinerie</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>The brainweb: phase synchronization and large-scale integration.</article-title>
<source>Nat Rev Neurosci</source>
<volume>2</volume>
<fpage>229</fpage>
<lpage>239</lpage>
<pub-id pub-id-type="pmid">11283746</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Visani1"><label>110</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Visani</surname>
<given-names>E</given-names>
</name>
<name><surname>Varotto</surname>
<given-names>G</given-names>
</name>
<name><surname>Binelli</surname>
<given-names>S</given-names>
</name>
<name><surname>Fratello</surname>
<given-names>L</given-names>
</name>
<name><surname>Franceschetti</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<year>2010</year>
<article-title>Photosensitive epilepsy: spectral and coherence analyses of EEG using 14 Hz intermittent photic stimulation.</article-title>
<source>Clin Neurophysiol</source>
<volume>121</volume>
<fpage>318</fpage>
<lpage>324</lpage>
<pub-id pub-id-type="pmid">20071221</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Siniatchkin1"><label>111</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Siniatchkin</surname>
<given-names>M</given-names>
</name>
<name><surname>Groppa</surname>
<given-names>S</given-names>
</name>
<name><surname>Jerosch</surname>
<given-names>B</given-names>
</name>
<name><surname>Muhle</surname>
<given-names>H</given-names>
</name>
<name><surname>Kurth</surname>
<given-names>C</given-names>
</name>
<etal></etal>
</person-group>
<year>2007</year>
<article-title>Spreading photoparoxysmal EEG response is associated with an abnormal cortical excitability pattern.</article-title>
<source>Brain</source>
<volume>130</volume>
<fpage>78</fpage>
<lpage>87</lpage>
<pub-id pub-id-type="pmid">17121743</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Robinson1"><label>112</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Robinson</surname>
<given-names>PA</given-names>
</name>
<name><surname>Rennie</surname>
<given-names>CJ</given-names>
</name>
<name><surname>Rowe</surname>
<given-names>DL</given-names>
</name>
</person-group>
<year>2002</year>
<article-title>Dynamics of large-scale brain activity in normal arousal states and epileptic seizures.</article-title>
<source>Phys Rev E Stat Nonlin Soft Matter Phys</source>
<volume>65</volume>
<fpage>041924</fpage>
<pub-id pub-id-type="pmid">12005890</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Dafilis1"><label>113</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Dafilis</surname>
<given-names>MP</given-names>
</name>
<name><surname>Frascoli</surname>
<given-names>F</given-names>
</name>
<name><surname>Cadusch</surname>
<given-names>PJ</given-names>
</name>
<name><surname>Liley</surname>
<given-names>DTJ</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Chaos and generalised multistability in a mesoscopic model of the electroencephalogram.</article-title>
<source>Physica D</source>
<volume>238</volume>
<fpage>1056</fpage>
<lpage>1060</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Dafilis2"><label>114</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Dafilis</surname>
<given-names>MP</given-names>
</name>
<name><surname>Liley</surname>
<given-names>DTJ</given-names>
</name>
<name><surname>Cadusch</surname>
<given-names>PJ</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Robust chaos in a model of the electroencephalogram: Implications for brain dynamics.</article-title>
<source>Chaos</source>
<volume>11</volume>
<fpage>474</fpage>
<lpage>478</lpage>
<pub-id pub-id-type="pmid">12779484</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Freeman2"><label>115</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Freeman</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<year>1987</year>
<article-title>Simulation of chaotic EEG patterns with a dynamic model of the olfactory system.</article-title>
<source>Biol Cybern</source>
<volume>56</volume>
<fpage>139</fpage>
<lpage>150</lpage>
<pub-id pub-id-type="pmid">3593783</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Freeman3"><label>116</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Freeman</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<year>1992</year>
<article-title>Tutorial in Neurobiology: From Single Neurons to Brain Chaos.</article-title>
<source>Int J Bifurcat Chaos</source>
<volume>2</volume>
<fpage>451</fpage>
<lpage>482</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Freeman4"><label>117</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Freeman</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Evidence from human scalp electroencephalograms of global chaotic itinerancy.</article-title>
<source>Chaos</source>
<volume>13</volume>
<fpage>1067</fpage>
<lpage>1077</lpage>
<pub-id pub-id-type="pmid">12946200</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Korn1"><label>118</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Korn</surname>
<given-names>H</given-names>
</name>
<name><surname>Faure</surname>
<given-names>P</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Is there chaos in the brain? II. Experimental evidence and related models.</article-title>
<source>Comptes Rendus Biologies</source>
<volume>326</volume>
<fpage>787</fpage>
<lpage>840</lpage>
<pub-id pub-id-type="pmid">14694754</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Makeig1"><label>119</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Makeig</surname>
<given-names>S</given-names>
</name>
</person-group>
<year>1993</year>
<article-title>Auditory event-related dynamics of the EEG spectrum and effects of exposure to tones.</article-title>
<source>Electroencephalogr Clin Neurophysiol</source>
<volume>86</volume>
<fpage>283</fpage>
<lpage>293</lpage>
<pub-id pub-id-type="pmid">7682932</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Miller1"><label>120</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Miller</surname>
<given-names>KJ</given-names>
</name>
<name><surname>Leuthardt</surname>
<given-names>EC</given-names>
</name>
<name><surname>Schalk</surname>
<given-names>G</given-names>
</name>
<name><surname>Rao</surname>
<given-names>RP</given-names>
</name>
<name><surname>Anderson</surname>
<given-names>NR</given-names>
</name>
<etal></etal>
</person-group>
<year>2007</year>
<article-title>Spectral changes in cortical surface potentials during motor movement.</article-title>
<source>J Neurosci</source>
<volume>27</volume>
<fpage>2424</fpage>
<lpage>2432</lpage>
<pub-id pub-id-type="pmid">17329441</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Zhang1"><label>121</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Zhang</surname>
<given-names>Y</given-names>
</name>
<name><surname>Wang</surname>
<given-names>X</given-names>
</name>
<name><surname>Bressler</surname>
<given-names>SL</given-names>
</name>
<name><surname>Chen</surname>
<given-names>Y</given-names>
</name>
<name><surname>Ding</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Prestimulus cortical activity is correlated with speed of visuomotor processing.</article-title>
<source>J Cogn Neurosci</source>
<volume>20</volume>
<fpage>1915</fpage>
<lpage>1925</lpage>
<pub-id pub-id-type="pmid">18370597</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Doesburg1"><label>122</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Doesburg</surname>
<given-names>SM</given-names>
</name>
<name><surname>Kitajo</surname>
<given-names>K</given-names>
</name>
<name><surname>Ward</surname>
<given-names>LM</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Increased gamma-band synchrony precedes switching of conscious perceptual objects in binocular rivalry.</article-title>
<source>Neuroreport</source>
<volume>16</volume>
<fpage>1139</fpage>
<lpage>1142</lpage>
<pub-id pub-id-type="pmid">16012336</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Babloyantz1"><label>123</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Babloyantz</surname>
<given-names>A</given-names>
</name>
<name><surname>Salazar</surname>
<given-names>JM</given-names>
</name>
<name><surname>Nicolis</surname>
<given-names>C</given-names>
</name>
</person-group>
<year>1985</year>
<article-title>Evidence of chaotic dynamics of brain activity during the sleep cycle.</article-title>
<source>Phys Lett A</source>
<volume>111</volume>
<fpage>152</fpage>
<lpage>156</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Gallez1"><label>124</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Gallez</surname>
<given-names>D</given-names>
</name>
<name><surname>Babloyantz</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>1991</year>
<article-title>Predictability of human EEG: a dynamical approach.</article-title>
<source>Biol Cybern</source>
<volume>64</volume>
<fpage>381</fpage>
<lpage>391</lpage>
<pub-id pub-id-type="pmid">2049414</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-DiPrisco1"><label>125</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Di Prisco</surname>
<given-names>GV</given-names>
</name>
<name><surname>Freeman</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<year>1985</year>
<article-title>Odor-related bulbar EEG spatial pattern analysis during appetitive conditioning in rabbits.</article-title>
<source>Behav Neurosci</source>
<volume>99</volume>
<fpage>964</fpage>
<lpage>978</lpage>
<pub-id pub-id-type="pmid">3843313</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Freeman5"><label>126</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Freeman</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Di Prisco</surname>
<given-names>GV</given-names>
</name>
</person-group>
<year>1986</year>
<article-title>Spatial patterns differences with discriminated odors manifest chaotic and limit cycles attractors in olfactory bulb of rabbits.</article-title>
<person-group person-group-type="editor"><name><surname>Palm</surname>
<given-names>G</given-names>
</name>
<name><surname>Aartsen</surname>
<given-names>A</given-names>
</name>
</person-group>
<source>Brain Theory</source>
<publisher-loc>Berlin</publisher-loc>
<publisher-name>Springer</publisher-name>
<fpage>97</fpage>
<lpage>119</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Babloyantz2"><label>127</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Babloyantz</surname>
<given-names>A</given-names>
</name>
<name><surname>Destexhe</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>1986</year>
<article-title>Low-dimensional chaos in an instance of epilepsy.</article-title>
<source>Proc Natl Acad Sci U S A</source>
<volume>83</volume>
<fpage>3513</fpage>
<lpage>3517</lpage>
<pub-id pub-id-type="pmid">3085091</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Iasemidis1"><label>128</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Iasemidis</surname>
<given-names>LD</given-names>
</name>
<name><surname>Sackellares</surname>
<given-names>JC</given-names>
</name>
</person-group>
<year>1991</year>
<article-title>The evolution with time of the spatial distribution of the largest Lyapunov exponent on the human epileptic cortex.</article-title>
<person-group person-group-type="editor"><name><surname>Duke</surname>
<given-names>D</given-names>
</name>
<name><surname>Pritchard</surname>
<given-names>W</given-names>
</name>
</person-group>
<source>Measuring Chaos in the Human Brain</source>
<publisher-loc>Singapore</publisher-loc>
<publisher-name>World Scientific</publisher-name>
<fpage>49</fpage>
<lpage>82</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Iasemidis2"><label>129</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Iasemidis</surname>
<given-names>LD</given-names>
</name>
<name><surname>Sackellares</surname>
<given-names>JC</given-names>
</name>
<name><surname>Zaveri</surname>
<given-names>HP</given-names>
</name>
<name><surname>Williams</surname>
<given-names>WJ</given-names>
</name>
</person-group>
<year>1990</year>
<article-title>Phase space topography and the Lyapunov exponent of electrocorticograms in partial seizures.</article-title>
<source>Brain Topogr</source>
<volume>2</volume>
<fpage>187</fpage>
<lpage>201</lpage>
<pub-id pub-id-type="pmid">2116818</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Lehnertz1"><label>130</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Lehnertz</surname>
<given-names>K</given-names>
</name>
<name><surname>Elger</surname>
<given-names>CE</given-names>
</name>
</person-group>
<year>1998</year>
<article-title>Can Epileptic Seizures be Predicted? Evidence from Nonlinear Time Series Analysis of Brain Electrical Activity.</article-title>
<source>Phys Rev Lett</source>
<volume>80</volume>
<fpage>5019</fpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-SteynRoss1"><label>131</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Steyn-Ross</surname>
<given-names>DA</given-names>
</name>
<name><surname>Steyn-Ross</surname>
<given-names>ML</given-names>
</name>
<name><surname>Sleigh</surname>
<given-names>JW</given-names>
</name>
<name><surname>Wilson</surname>
<given-names>MT</given-names>
</name>
<name><surname>Gillies</surname>
<given-names>IP</given-names>
</name>
<etal></etal>
</person-group>
<year>2005</year>
<article-title>The sleep cycle modelled as a cortical phase transition.</article-title>
<source>J Biol Phys</source>
<volume>31</volume>
<fpage>547</fpage>
<lpage>569</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Spiegler2"><label>132</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Spiegler</surname>
<given-names>A</given-names>
</name>
<name><surname>Graimann</surname>
<given-names>B</given-names>
</name>
<name><surname>Pfurtscheller</surname>
<given-names>G</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Phase coupling between different motor areas during tongue-movement imagery.</article-title>
<source>Neurosci Lett</source>
<volume>369</volume>
<fpage>50</fpage>
<lpage>54</lpage>
<pub-id pub-id-type="pmid">15380306</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Blinowska1"><label>133</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Blinowska</surname>
<given-names>KJ</given-names>
</name>
<name><surname>Kus</surname>
<given-names>R</given-names>
</name>
<name><surname>Kaminski</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>2004</year>
<article-title>Granger causality and information flow in multivariate processes.</article-title>
<source>Phys Rev E Stat Nonlin Soft Matter Phys</source>
<volume>70</volume>
<fpage>050902</fpage>
<pub-id pub-id-type="pmid">15600583</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kaminski1"><label>134</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kaminski</surname>
<given-names>M</given-names>
</name>
<name><surname>Ding</surname>
<given-names>M</given-names>
</name>
<name><surname>Truccolo</surname>
<given-names>WA</given-names>
</name>
<name><surname>Bressler</surname>
<given-names>SL</given-names>
</name>
</person-group>
<year>2001</year>
<article-title>Evaluating causal relations in neural systems: granger causality, directed transfer function and statistical assessment of significance.</article-title>
<source>Biol Cybern</source>
<volume>85</volume>
<fpage>145</fpage>
<lpage>157</lpage>
<pub-id pub-id-type="pmid">11508777</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Breakspear2"><label>135</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Breakspear</surname>
<given-names>M</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<year>2007</year>
<article-title>Neuronal dynamics and brain connectivity.</article-title>
<person-group person-group-type="editor"><name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
</person-group>
<source>Handbook of Brain Connectivity</source>
<publisher-loc>Berlin, Heidelberg</publisher-loc>
<publisher-name>Springer</publisher-name>
<fpage>3</fpage>
<lpage>64</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Marten1"><label>136</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Marten</surname>
<given-names>F</given-names>
</name>
<name><surname>Rodrigues</surname>
<given-names>S</given-names>
</name>
<name><surname>Benjamin</surname>
<given-names>O</given-names>
</name>
<name><surname>Richardson</surname>
<given-names>MP</given-names>
</name>
<name><surname>Terry</surname>
<given-names>JR</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Onset of polyspike complexes in a mean-field model of human electroencephalography and its application to absence epilepsy.</article-title>
<source>Philos Transact A Math Phys Eng Sci</source>
<volume>367</volume>
<fpage>1145</fpage>
<lpage>1161</lpage>
<pub-id pub-id-type="pmid">19218156</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Rodrigues1"><label>137</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Rodrigues</surname>
<given-names>S</given-names>
</name>
<name><surname>Barton</surname>
<given-names>D</given-names>
</name>
<name><surname>Szalai</surname>
<given-names>R</given-names>
</name>
<name><surname>Benjamin</surname>
<given-names>O</given-names>
</name>
<name><surname>Richardson</surname>
<given-names>MP</given-names>
</name>
<etal></etal>
</person-group>
<year>2009</year>
<article-title>Transitions to spike-wave oscillations and epileptic dynamics in a human cortico-thalamic mean-field model.</article-title>
<source>J Comput Neurosci</source>
<volume>27</volume>
<fpage>507</fpage>
<lpage>526</lpage>
<pub-id pub-id-type="pmid">19499316</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Watson1"><label>138</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Watson</surname>
<given-names>AB</given-names>
</name>
<name><surname>Nachmias</surname>
<given-names>J</given-names>
</name>
</person-group>
<year>1977</year>
<article-title>Patterns of temporal interaction in the detection of gratings.</article-title>
<source>Vision Res</source>
<volume>17</volume>
<fpage>893</fpage>
<lpage>902</lpage>
<pub-id pub-id-type="pmid">595394</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wilke1"><label>139</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wilke</surname>
<given-names>M</given-names>
</name>
<name><surname>Mueller</surname>
<given-names>KM</given-names>
</name>
<name><surname>Leopold</surname>
<given-names>DA</given-names>
</name>
</person-group>
<year>2009</year>
<article-title>Neural activity in the visual thalamus reflects perceptual suppression.</article-title>
<source>Proc Natl Acad Sci U S A</source>
<volume>106</volume>
<fpage>9465</fpage>
<lpage>9470</lpage>
<pub-id pub-id-type="pmid">19458249</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Steriade1"><label>140</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Steriade</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Coherent oscillations and short-term plasticity in corticothalamic networks.</article-title>
<source>Trends Cogn Sci</source>
<volume>22</volume>
<fpage>337</fpage>
<lpage>345</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Gabernet1"><label>141</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Gabernet</surname>
<given-names>L</given-names>
</name>
<name><surname>Jadhav</surname>
<given-names>SP</given-names>
</name>
<name><surname>Feldman</surname>
<given-names>DE</given-names>
</name>
<name><surname>Carandini</surname>
<given-names>M</given-names>
</name>
<name><surname>Scanziani</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>2005</year>
<article-title>Somatosensory integration controlled by dynamic thalamocortical feed-forward inhibition.</article-title>
<source>Neuron</source>
<volume>48</volume>
<fpage>315</fpage>
<lpage>327</lpage>
<pub-id pub-id-type="pmid">16242411</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kimura1"><label>142</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Kimura</surname>
<given-names>F</given-names>
</name>
<name><surname>Itami</surname>
<given-names>C</given-names>
</name>
<name><surname>Ikezoe</surname>
<given-names>K</given-names>
</name>
<name><surname>Tamura</surname>
<given-names>H</given-names>
</name>
<name><surname>Fujita</surname>
<given-names>I</given-names>
</name>
<etal></etal>
</person-group>
<year>2010</year>
<article-title>Fast activation of feedforward inhibitory neurons from thalamic input and its relevance to the regulation of spike sequences in the barrel cortex.</article-title>
<source>J Physiol</source>
<volume>588</volume>
<fpage>2769</fpage>
<lpage>2787</lpage>
<pub-id pub-id-type="pmid">20530116</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wehr1"><label>143</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wehr</surname>
<given-names>M</given-names>
</name>
<name><surname>Zador</surname>
<given-names>AM</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Balanced inhibition underlies tuning and sharpens spike timing in auditory cortex.</article-title>
<source>Nature</source>
<volume>426</volume>
<fpage>442</fpage>
<lpage>446</lpage>
<pub-id pub-id-type="pmid">14647382</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Miller2"><label>144</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Miller</surname>
<given-names>KD</given-names>
</name>
</person-group>
<year>2003</year>
<article-title>Understanding layer 4 of the cortical circuit: a model based on cat V1.</article-title>
<source>Cereb Cortex</source>
<volume>13</volume>
<fpage>73</fpage>
<lpage>82</lpage>
<pub-id pub-id-type="pmid">12466218</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Priebe1"><label>145</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Priebe</surname>
<given-names>NJ</given-names>
</name>
<name><surname>Ferster</surname>
<given-names>D</given-names>
</name>
</person-group>
<year>2008</year>
<article-title>Inhibition, spike threshold, and stimulus selectivity in primary visual cortex.</article-title>
<source>Neuron</source>
<volume>57</volume>
<fpage>482</fpage>
<lpage>497</lpage>
<pub-id pub-id-type="pmid">18304479</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Sun1"><label>146</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Sun</surname>
<given-names>QQ</given-names>
</name>
<name><surname>Huguenard</surname>
<given-names>JR</given-names>
</name>
<name><surname>Prince</surname>
<given-names>DA</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Barrel cortex microcircuits: thalamocortical feedforward inhibition in spiny stellate cells is mediated by a small number of fast-spiking interneurons.</article-title>
<source>J Neurosci</source>
<volume>26</volume>
<fpage>1219</fpage>
<lpage>1230</lpage>
<pub-id pub-id-type="pmid">16436609</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Otnes1"><label>147</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Otnes</surname>
<given-names>RK</given-names>
</name>
<name><surname>Enochson</surname>
<given-names>LD</given-names>
</name>
</person-group>
<year>1972</year>
<source>Digital time series analysis</source>
<publisher-loc>New York</publisher-loc>
<publisher-name>Wiley</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Chen1"><label>148</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Chen</surname>
<given-names>ZM</given-names>
</name>
<name><surname>Djidjeli</surname>
<given-names>K</given-names>
</name>
<name><surname>Price</surname>
<given-names>WG</given-names>
</name>
</person-group>
<year>2006</year>
<article-title>Computing Lyapunov exponents based on the solution expression of the variational system.</article-title>
<source>Appl Math Comput</source>
<volume>174</volume>
<fpage>982</fpage>
<lpage>996</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Baier1"><label>149</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Baier</surname>
<given-names>G</given-names>
</name>
<name><surname>Klein</surname>
<given-names>M</given-names>
</name>
</person-group>
<year>1991</year>
<source>A Chaotic hierarchy</source>
<publisher-loc>Singapore; Teaneck, NJ</publisher-loc>
<publisher-name>World Scientific</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Kantz1"><label>150</label>
<element-citation publication-type="book"> <person-group person-group-type="author"><name><surname>Kantz</surname>
<given-names>H</given-names>
</name>
<name><surname>Schreiber</surname>
<given-names>T</given-names>
</name>
</person-group>
<year>2004</year>
<source>Nonlinear time series analysis</source>
<publisher-loc>Cambridge; New York</publisher-loc>
<publisher-name>Cambridge University Press</publisher-name>
<size units="page"></size>
</element-citation>
</ref>
<ref id="pcbi.1002298-Jampolsky1"><label>151</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Jampolsky</surname>
<given-names>A</given-names>
</name>
</person-group>
<year>1994</year>
<article-title>Consequences of retinal image clarity versus occlusion (absent) versus diffusion.</article-title>
<source>Trans Am Ophthalmol Soc</source>
<volume>92</volume>
<fpage>349</fpage>
<lpage>373; discussion 373–346</lpage>
<pub-id pub-id-type="pmid">7886872</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Leijten1"><label>152</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Leijten</surname>
<given-names>FS</given-names>
</name>
<name><surname>Dekker</surname>
<given-names>E</given-names>
</name>
<name><surname>Spekreijse</surname>
<given-names>H</given-names>
</name>
<name><surname>Kasteleijn-Nolst Trenite</surname>
<given-names>DG</given-names>
</name>
<name><surname>Van Emde Boas</surname>
<given-names>W</given-names>
</name>
</person-group>
<year>1998</year>
<article-title>Light diffusion in photosensitive epilepsy.</article-title>
<source>Electroencephalogr Clin Neurophysiol</source>
<volume>106</volume>
<fpage>387</fpage>
<lpage>391</lpage>
<pub-id pub-id-type="pmid">9680149</pub-id>
</element-citation>
</ref>
<ref id="pcbi.1002298-Wolf1"><label>153</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Wolf</surname>
<given-names>A</given-names>
</name>
<name><surname>Swift</surname>
<given-names>JB</given-names>
</name>
<name><surname>Swinney</surname>
<given-names>HL</given-names>
</name>
<name><surname>Vastano</surname>
<given-names>JA</given-names>
</name>
</person-group>
<year>1985</year>
<article-title>Determining Lyapunov Exponents from a Time-Series.</article-title>
<source>Physica D</source>
<volume>16</volume>
<fpage>285</fpage>
<lpage>317</lpage>
</element-citation>
</ref>
<ref id="pcbi.1002298-Atay3"><label>154</label>
<element-citation publication-type="journal"> <person-group person-group-type="author"><name><surname>Atay</surname>
<given-names>FM</given-names>
</name>
<name><surname>Altintas</surname>
<given-names>Y</given-names>
</name>
</person-group>
<year>1999</year>
<article-title>Recovering smooth dynamics from time series with the aid of recurrence plots.</article-title>
<source>Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics</source>
<volume>59</volume>
<fpage>6593</fpage>
<lpage>6598</lpage>
<pub-id pub-id-type="pmid">11969647</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
<affiliations><list><country><li>Allemagne</li>
</country>
<region><li>District de Leipzig</li>
<li>Saxe (Land)</li>
</region>
<settlement><li>Leipzig</li>
</settlement>
</list>
<tree><country name="Allemagne"><region name="Saxe (Land)"><name sortKey="Spiegler, Andreas" sort="Spiegler, Andreas" uniqKey="Spiegler A" first="Andreas" last="Spiegler">Andreas Spiegler</name>
</region>
<name sortKey="Atay, Fatihcan M" sort="Atay, Fatihcan M" uniqKey="Atay F" first="Fatihcan M." last="Atay">Fatihcan M. Atay</name>
<name sortKey="Haueisen, Jens" sort="Haueisen, Jens" uniqKey="Haueisen J" first="Jens" last="Haueisen">Jens Haueisen</name>
<name sortKey="Knosche, Thomas R" sort="Knosche, Thomas R" uniqKey="Knosche T" first="Thomas R." last="Knösche">Thomas R. Knösche</name>
<name sortKey="Schwab, Karin" sort="Schwab, Karin" uniqKey="Schwab K" first="Karin" last="Schwab">Karin Schwab</name>
<name sortKey="Spiegler, Andreas" sort="Spiegler, Andreas" uniqKey="Spiegler A" first="Andreas" last="Spiegler">Andreas Spiegler</name>
</country>
</tree>
</affiliations>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/HapticV1/Data/Ncbi/Merge
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001D53 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Ncbi/Merge/biblio.hfd -nk 001D53 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Ticri/CIDE |area= HapticV1 |flux= Ncbi |étape= Merge |type= RBID |clé= PMC:3245303 |texte= Modeling Brain Resonance Phenomena Using a Neural Mass Model }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Ncbi/Merge/RBID.i -Sk "pubmed:22215992" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Ncbi/Merge/biblio.hfd \ | NlmPubMed2Wicri -a HapticV1
![]() | This area was generated with Dilib version V0.6.23. | ![]() |