Period-adding bifurcations and chaos in a periodically stimulated excitable neural relaxation oscillator

The response of an excitable neuron to trains of electrical spikes is relevant to the understanding of the neural code. In this paper we study a neurobiologically motivated relaxation oscillator, with appropriately identified fast and slow coordinates, that admits an explicit mathematical analysis. An application of geometric singular perturbation theory shows the existence of an attracting invariant manifold which is used to construct the Fenichel normal form for the system. This facilitates the calculation of the response of the system to pulsatile stimulation and allows the construction of a so-called extended isochronal map. The isochronal map is shown to have a single discontinuity and be of a type that can admit three types of response: mode-locked, quasi-periodic and chaotic. The bifurcation structure of the system is seen to be extremely rich and supports period-adding bifurcations separated by windows of both chaos and periodicity. A bifurcation analysis of the isochronal map is presented in conjunction with a description of the various routes to chaos in this system

Oscillator-based neuronal modeling for seizure progression investigation and seizure control strategy

The coupled oscillator model has previously been used for the simulation of neuronal activities in in vitro rat hippocampal slice seizure data and the evaluation of seizure suppression algorithms. Each model unit can be described as either an oscillator which can generate action potential spike trains without inputs, or a threshold-based unit. With the change of only one parameter, each unit can either be an oscillator or a threshold-based spiking unit. This would eliminate the need for a new set of equations for each type of unit. Previous analysis has suggested that long kernel duration and imbalance of inhibitory feedback can cause the system to intermittently transition into and out of ictal activities. The state transitions of seizure-like events were investigated here; specifically, how the system excitability may change when the system undergoes transitions in the preictal and postictal processes. Analysis showed that the area of the excitation kernel is positively correlated with the mean firing rate of the ictal activity. The kernel duration is also correlated to the amount of ictal activity. The transition into ictal activity involved the escape from the saddle point foci in the state space trajectory identified by using Newton\u27s method. The ability to accurately anticipate and suppress seizures is an important endeavor that has tremendous impact on improving the quality of lives for epileptic patients. The stimulation studies have suggested that an electrical stimulation strategy that uses the intrinsic high complexity dynamics of the biological system may be more effective in reducing the duration of seizure-like activities in the computer model. In this research, we evaluate this strategy on an in vitro rat hippocampal slice magnesium-free model. Simulated postictal field potential data generated by an oscillator-based hippocampal network model was applied to the CA1 region of the rat hippocampal slices through a multi-electrode array (MEA) system. It was found to suppress and delay the onset of future seizures temporarily. The average inter-seizure time was found to be significantly prolonged after postictal stimulation when compared to the negative control trials and bipolar square wave signals. The result suggests that neural signal-based stimulation related to resetting may be suitable for seizure control in the clinical environment

Stochastic neural network dynamics: synchronisation and control

Biological brains exhibit many interesting and complex behaviours. Understanding of the mechanisms behind brain behaviours is critical for continuing advancement in fields of research such as artificial intelligence and medicine. In particular, synchronisation of neuronal firing is associated with both improvements to and degeneration of the brain’s performance; increased synchronisation can lead to enhanced information-processing or neurological disorders such as epilepsy and Parkinson’s disease. As a result, it is desirable to research under which conditions synchronisation arises in neural networks and the possibility of controlling its prevalence. Stochastic ensembles of FitzHugh-Nagumo elements are used to model neural networks for numerical simulations and bifurcation analysis. The FitzHugh-Nagumo model is employed because of its realistic representation of the flow of sodium and potassium ions in addition to its advantageous property of allowing phase plane dynamics to be observed. Network characteristics such as connectivity, configuration and size are explored to determine their influences on global synchronisation generation in their respective systems. Oscillations in the mean-field are used to detect the presence of synchronisation over a range of coupling strength values. To ensure simulation efficiency, coupling strengths between neurons that are identical and fixed with time are investigated initially. Such networks where the interaction strengths are fixed are referred to as homogeneously coupled. The capacity of controlling and altering behaviours produced by homogeneously coupled networks is assessed through the application of weak and strong delayed feedback independently with various time delays. To imitate learning, the coupling strengths later deviate from one another and evolve with time in networks that are referred to as heterogeneously coupled. The intensity of coupling strength fluctuations and the rate at which coupling strengths converge to a desired mean value are studied to determine their impact upon synchronisation performance. The stochastic delay differential equations governing the numerically simulated networks are then converted into a finite set of deterministic cumulant equations by virtue of the Gaussian approximation method. Cumulant equations for maximal and sub-maximal connectivity are used to generate two-parameter bifurcation diagrams on the noise intensity and coupling strength plane, which provides qualitative agreement with numerical simulations. Analysis of artificial brain networks, in respect to biological brain networks, are discussed in light of recent research in sleep theor

Neuronal oscillations: from single-unit activity to emergent dynamics and back

L’objectiu principal d’aquesta tesi és avançar en la comprensió del processament d’informació en xarxes neuronals en presència d’oscil lacions subumbrals. La majoria de neurones propaguen la seva activitat elèctrica a través de sinapsis químiques que són activades, exclusivament, quan el corrent elèctric que les travessa supera un cert llindar. És per aquest motiu que les descàrregues ràpides i intenses produïdes al soma neuronal, els anomenats potencials d’acció, són considerades la unitat bàsica d’informació neuronal, és a dir, el senyal mínim i necessari per a iniciar la comunicació entre dues neurones. El codi neuronal és entès, doncs, com un llenguatge binari que expressa qualsevol missatge (estímul sensorial, memòries, etc.) en un tren de potencials d’acció. Tanmateix, cap funció cognitiva rau en la dinàmica d’una única neurona. Circuits de milers de neurones connectades entre sí donen lloc a determinats ritmes, palesos en registres d’activitat colectiva com els electroencefalogrames (EEG) o els potencials de camp local (LFP). Si els potencials d’acció de cada cèl lula, desencadenats per fluctuacions estocàstiques de les corrents sinàptiques, no assolissin un cert grau de sincronia, no apareixeria aquesta periodicitat a nivell de xarxa. Per tal de poder entendre si aquests ritmes intervenen en el codi neuronal hem estudiat tres situacions. Primer, en el Capítol 2, hem mostrat com una cadena oberta de neurones amb un potencial de membrana intrínsecament oscil latori filtra un senyal periòdic arribant per un dels extrems. La resposta de cada neurona (pulsar o no pulsar) depèn de la seva fase, de forma que cada una d’elles rep un missatge filtrat per la precedent. A més, cada potencial d’acció presinàptic provoca un canvi de fase en la neurona postsinàptica que depèn de la seva posició en l’espai de fases. Els períodes d’entrada capaços de sincronitzar les oscil lacions subumbrals són aquells que mantenen la fase d’arribada dels potencials d’acció fixa al llarg de la cadena. Per tal de què el missatge arribi intacte a la darrera neurona cal, a més a més, que aquesta fase permeti la descàrrega del voltatge transmembrana. En segon cas, hem estudiat una xarxa neuronal amb connexions tant a veïns propers com de llarg abast, on les oscil lacions subumbrals emergeixen de l’activitat col lectiva reflectida en els corrents sinàptics (o equivalentment en el LFP). Les neurones inhibidores aporten un ritme a l’excitabilitat de la xarxa, és a dir, que els episodis en què la inhibició és baixa, la probabilitat d’una descàrrega global de la població neuronal és alta. En el Capítol 3 mostrem com aquest ritme implica l’aparició d’una bretxa en la freqüència de descàrrega de les neurones: o bé polsen espaiadament en el temps o bé en ràfegues d’elevada intensitat. La fase del LFP determina l’estat de la xarxa neuronal codificant l’activitat de la població: els mínims indiquen la descàrrega simultània de moltes neurones que, ocasionalment, han superat el llindar d’excitabilitat degut a un decreixement global de la inhibició, mentre que els màxims indiquen la coexistència de ràfegues en diferents punts de la xarxa degut a decreixements locals de la inhibició en estats globals d’excitació. Aquesta dinàmica és possible gràcies al domini de la inhibició sobre l’excitació. En el Capítol 4 considerem acoblament entre dues xarxes neuronals per tal d’estudiar la interacció entre ritmes diferents. Les oscil lacions indiquen recurrència en la sincronització de l’activitat col lectiva, de manera que durant aquestes finestres temporals una població optimitza el seu impacte en una xarxa diana. Quan el ritme de la població receptora i el de l’emissora difereixen significativament, l’eficiència en la comunicació decreix, ja que les fases de màxima resposta de cada senyal LFP no mantenen una diferència constant entre elles. Finalment, en el Capítol 5 hem estudiat com les oscil lacions col lectives pròpies de l’estat de son donen lloc al fenomen de coherència estocàstica. Per a una intensitat òptima del soroll, modulat per l’excitabilitat de la xarxa, el LFP assoleix una regularitat màxima donant lloc a un període refractari de la població neuronal. En resum, aquesta Tesi mostra escenaris d’interacció entre els potencials d’acció, característics de la dinàmica de neurones individuals, i les oscil lacions subumbrals, fruit de l’acoblament entre les cèl lules i ubiqües en la dinàmica de poblacions neuronals. Els resultats obtinguts aporten funcionalitat a aquests ritmes emergents, agents sincronitzadors i moduladors de les descàrregues neuronals i reguladors de la comunicació entre xarxes neuronals.The main objective of this thesis is to better understand information processing in neuronal networks in the presence of subthreshold oscillations. Most neurons propagate their electrical activity via chemical synapses, which are only activated when the electric current that passes through them surpasses a certain threshold. Therefore, fast and intense discharges produced at the neuronal soma (the action potentials or spikes) are considered the basic unit of neuronal information. The neuronal code is understood, then, as a binary language that expresses any message (sensory stimulus, memories, etc.) in a train of action potentials. Circuits of thousands of interconnected neurons give rise to certain rhythms, revealed in collective activity measures such as electroencephalograms (EEG) and local field potentials (LFP). Synchronization of action potentials of each cell, triggered by stochastic fluctuations of the synaptic currents, cause this periodicity at the network level.To understand whether these rhythms are involved in the neuronal code we studied three situations. First, in Chapter 2, we showed how an open chain of neurons with an intrinsically oscillatory membrane potential filters a periodic signal coming from one of its ends. The response of each neuron (to spike or not) depends on its phase, so that each cell receives a message filtered by the preceding one. Each presynaptic action potential causes a phase change in the postsynaptic neuron, which depends on its position in the phase space. Those incoming periods that are able to synchronize the subthreshold oscillations, keep the phase of arrival of action potentials fixed along the chain. The original message reaches intact the last neuron provided that this phase allows the discharge of the transmembrane voltage.I the second case, we studied a neuronal network with connections to both long range and close neighbors, in which the subthreshold oscillations emerge from the collective activity apparent in the synaptic currents. The inhibitory neurons provide a rhythm to the excitability of the network. When inhibition is low, the likelihood of a global discharge of the neuronal population is high. In Chapter 3 we show how this rhythm causes a gap in the discharge frequency of neurons: either they pulse single spikes or they fire bursts of high intensity. The LFP phase determines the state of the neuronal network, coding the activity of the population: its minima indicate the simultaneous discharge of many neurons, while its maxima indicate the coexistence of bursts due to local decreases of inhibition at global states of excitation. In Chapter 4 we consider coupling between two neural networks in order to study the interaction between different rhythms. The oscillations indicate recurrence in the synchronization of collective activity, so that during these time windows a population optimizes its impact on a target network. When the rhythm of the emitter and receiver population differ significantly, the communication efficiency decreases as the phases of maximum response of each LFP signal do not maintain a constant difference between them.Finally, in Chapter 5 we studied how oscillations typical of the collective sleep state give rise to stochastic coherence. For an optimal noise intensity, modulated by the excitability of the network, the LFP reaches a maximal regularity leading to a refractory period of the neuronal population.In summary, this Thesis shows scenarios of interaction between action potentials, characteristics of the dynamics of individual neurons, and the subthreshold oscillations, outcome of the coupling between the cells and ubiquitous in the dynamics of neuronal populations . The results obtained provide functionality to these emerging rhythms, triggers of synchronization and modulator agents of the neuronal discharges and regulators of the communication between neuronal networks

Electromagnetic Precursors in Complex Layered Media

The dynamical evolution of an ultrawideband electromagnetic pulse as it propagates through a temporally dispersive and attenuative medium is a classical problem in electromagnetic wave theory with considerable practical importance dating back to seminal works conducted in 1914. With the use of modern asymptotic theory and numerical techniques, propagation of canonical pulses into complex (attenuative and dispersive) media have been analyzed and recently extended to nonlinear materials. The materials of interest for this research are modeled after realistic biological tissues. The mathematically rigorous and more accurate physical model of electromagnetic energy transfer into the biological materials modeled will be used as input to the FitzHugh-Nagumo circuit equivalent model for an excitable neuron. This detailed analysis will provide a new point of view to working groups and standardization committees in the eld of non-ionizing radiation safety that is based on so-called athermal effects

Low-dimensional models of single neurons: A review

The classical Hodgkin-Huxley (HH) point-neuron model of action potential generation is four-dimensional. It consists of four ordinary differential equations describing the dynamics of the membrane potential and three gating variables associated to a transient sodium and a delayed-rectifier potassium ionic currents. Conductance-based models of HH type are higher-dimensional extensions of the classical HH model. They include a number of supplementary state variables associated with other ionic current types, and are able to describe additional phenomena such as sub-threshold oscillations, mixed-mode oscillations (subthreshold oscillations interspersed with spikes), clustering and bursting. In this manuscript we discuss biophysically plausible and phenomenological reduced models that preserve the biophysical and/or dynamic description of models of HH type and the ability to produce complex phenomena, but the number of effective dimensions (state variables) is lower. We describe several representative models. We also describe systematic and heuristic methods of deriving reduced models from models of HH type

Nonlinear Dynamics, Synchronisation and Chaos in Coupled FHN Cardiac and Neural Cells

Physiological systems are amongst the most challenging systems to investigate from a mathematically based approach. The eld of mathematical biology is a relatively recent one when compared to physics. In this thesis I present an introduction to the physiological aspects needed to gain access to both cardiac and neural systems for a researcher trained in a mathematically based discipline. By using techniques from nonlinear dynamical systems theory I show a number of results that have implications for both neural and cardiac cells. Examining a reduced model of an excitable biological oscillator I show how rich the dynamical behaviour of such systems can be when coupled together. Quantifying the dynamics of coupled cells in terms of synchronisation measures is treated at length. Most notably it is shown that for cells that themselves cannot admit chaotic solutions, communication between cells be it through electrical coupling or synaptic like coupling, can lead to the emergence of chaotic behaviour. I also show that in the presence of emergent chaos one nds great variability in intervals of activity between the constituent cells. This implies that chaos in both cardiac and neural systems can be a direct result of interactions between the constituent cells rather than intrinsic to the cells themselves. Furthermore the ubiquity of chaotic solutions in the coupled systems may be a means of information production and signaling in neural systems