Designing a planar vector field to investigate the role of a slow variable in an enhanced mean-field model during general anesthesia.

International audienceLocal mean-field models (MFMs) describe regional brain activities by some connected differential equations. In an overall view, constituting variables of these differential equations can be divided to very fast, fast and slow variables. In this article we propose a method that can be used to determine role of a slow variable in behavior of MFMs. Very fast variables can be adiabatically removed from the equations. Isoclines of fast and slow variables and their corresponding vector field can provide valuable information about model behavior and role of the slow variable in it. The vector field of our interested MFM that is an enhanced MFM designed specially for general anesthesia, is a 3D field (one slow and two fast variables) and it is not so convenient for visually inspecting the role of the slow variable in this model. To afford this problem we design a 2D (planar) vector filed that only considers the slow variable and one of the fast variables

Removing Ocular Artifacts from EEG Signals Using Adaptive Filtering and ARMAX Modeling

EEG signal is one of the oldest measures of brain activity that has been used vastly for clinical diagnoses and biomedical researches. However, EEG signals are highly contaminated with various artifacts, both from the subject and from equipment interferences. Among these various kinds of artifacts, ocular noise is the most important one. Since many applications such as BCI require online and real-time processing of EEG signal, it is ideal if the removal of artifacts is performed in an online fashion. Recently, some methods for online ocular artifact removing have been proposed. One of these methods is ARMAX modeling of EEG signal. This method assumes that the recorded EEG signal is a combination of EOG artifacts and the background EEG. Then the background EEG is estimated via estimation of ARMAX parameters. The other recently proposed method is based on adaptive filtering. This method uses EOG signal as the reference input and subtracts EOG artifacts from recorded EEG signals. In this paper we investigate the efficiency of each method for removing of EOG artifacts. A comparison is made between these two methods. Our undertaken conclusion from this comparison is that adaptive filtering method has better results compared with the results achieved by ARMAX modeling

Person Identification by Using AR Model for EEG Signals

A direct connection between ElectroEncephaloGram (EEG) and the genetic information of individuals has been investigated by neurophysiologists and psychiatrists since 1960’s; and it opens a new research area in the science. This paper focuses on the person identification based on feature extracted from the EEG which can show a direct connection between EEG and the genetic information of subjects. In this work the full EO EEG signal of healthy individuals are estimated by an autoregressive (AR) model and the AR parameters are extracted as features. Here for feature vector constitution, two methods have been proposed; in the first method the extracted parameters of each channel are used as a feature vector in the classification step which employs a competitive neural network and in the second method a combination of different channel parameters are used as a feature vector. Correct classification scores at the range of 80% to 100% reveal the potential of our approach for person classification/identification and are in agreement to the previous researches showing evidence that the EEG signal carries genetic information. The novelty of this work is in the combination of AR parameters and the network type (competitive network) that we have used. A comparison between the first and the second approach imply preference of the second one

From oscillatory transcranial current stimulation to scalp EEG changes: a biophysical and physiological modeling study.

International audienceBoth biophysical and neurophysiological aspects need to be considered to assess the impact of electric fields induced by transcranial current stimulation (tCS) on the cerebral cortex and the subsequent effects occurring on scalp EEG. The objective of this work was to elaborate a global model allowing for the simulation of scalp EEG signals under tCS. In our integrated modeling approach, realistic meshes of the head tissues and of the stimulation electrodes were first built to map the generated electric field distribution on the cortical surface. Secondly, source activities at various cortical macro-regions were generated by means of a computational model of neuronal populations. The model parameters were adjusted so that populations generated an oscillating activity around 10 Hz resembling typical EEG alpha activity. In order to account for tCS effects and following current biophysical models, the calculated component of the electric field normal to the cortex was used to locally influence the activity of neuronal populations. Lastly, EEG under both spontaneous and tACS-stimulated (transcranial sinunoidal tCS from 4 to 16 Hz) brain activity was simulated at the level of scalp electrodes by solving the forward problem in the aforementioned realistic head model. Under the 10 Hz-tACS condition, a significant increase in alpha power occurred in simulated scalp EEG signals as compared to the no-stimulation condition. This increase involved most channels bilaterally, was more pronounced on posterior electrodes and was only significant for tACS frequencies from 8 to 12 Hz. The immediate effects of tACS in the model agreed with the post-tACS results previously reported in real subjects. Moreover, additional information was also brought by the model at other electrode positions or stimulation frequency. This suggests that our modeling approach can be used to compare, interpret and predict changes occurring on EEG with respect to parameters used in specific stimulation configurations

Effects of transcranial Direct Current Stimulation (tDCS) on cortical activity: A computational modeling study.

International audienceAlthough it is well-admitted that transcranial Direct Current Stimulation (tDCS) allows for interacting with brain endogenous rhythms, the exact mechanisms by which externally-applied fields modulate the activity of neurons remain elusive. In this study a novel computational model (a neural mass model including subpopulations of pyramidal cells and inhibitory interneurons mediating synaptic currents with either slow or fast kinetics) of the cerebral cortex was elaborated to investigate the local effects of tDCS on neuronal populations based on an in-vivo experimental study. Model parameters were adjusted to reproduce evoked potentials (EPs) recorded from the somatosensory cortex of the rabbit in response to air-puffs applied on the whiskers. EPs were simulated under control condition (no tDCS) as well as under anodal and cathodal tDCS fields. Results first revealed that a feed-forward inhibition mechanism must be included in the model for accurate simulation of actual EPs (peaks and latencies). Interestingly, results revealed that externally-applied fields are also likely to affect interneurons. Indeed, when interneurons get polarized then the characteristics of simulated EPs become closer to those of real EPs. In particular, under anodal tDCS condition, more realistic EPs could be obtained when pyramidal cells were depolarized and, simultaneously, slow (resp. fast) interneurons became de- (resp. hyper-) polarized. Geometrical characteristics of interneurons might provide some explanations for this effect

A neural population model of the bi-phasic EEG-power spectrum during general anaesthesia

International audienceThe neuronal mechanisms of general anaesthesia are still poorly understood, though the induction of analgesia, amnesia, immobility and loss of consciousness by anaesthetic agents is well-established in hospital practice. To shed some light onto these mysterious effects, the chapter analyzes mathematically a neural field model describing the neural population dynamics by an integro-differential equation. The power spectrum is derived and compared to experimental results

Population based models of cortical drug response: insights from anaesthesia

A great explanatory gap lies between the molecular pharmacology of psychoactive agents and the neurophysiological changes they induce, as recorded by neuroimaging modalities. Causally relating the cellular actions of psychoactive compounds to their influence on population activity is experimentally challenging. Recent developments in the dynamical modelling of neural tissue have attempted to span this explanatory gap between microscopic targets and their macroscopic neurophysiological effects via a range of biologically plausible dynamical models of cortical tissue. Such theoretical models allow exploration of neural dynamics, in particular their modification by drug action. The ability to theoretically bridge scales is due to a biologically plausible averaging of cortical tissue properties. In the resulting macroscopic neural field, individual neurons need not be explicitly represented (as in neural networks). The following paper aims to provide a non-technical introduction to the mean field population modelling of drug action and its recent successes in modelling anaesthesia

Scaling Effects and Spatio-Temporal Multilevel Dynamics in Epileptic Seizures

Epileptic seizures are one of the most well-known dysfunctions of the nervous system. During a seizure, a highly synchronized behavior of neural activity is observed that can cause symptoms ranging from mild sensual malfunctions to the complete loss of body control. In this paper, we aim to contribute towards a better understanding of the dynamical systems phenomena that cause seizures. Based on data analysis and modelling, seizure dynamics can be identified to possess multiple spatial scales and on each spatial scale also multiple time scales. At each scale, we reach several novel insights. On the smallest spatial scale we consider single model neurons and investigate early-warning signs of spiking. This introduces the theory of critical transitions to excitable systems. For clusters of neurons (or neuronal regions) we use patient data and find oscillatory behavior and new scaling laws near the seizure onset. These scalings lead to substantiate the conjecture obtained from mean-field models that a Hopf bifurcation could be involved near seizure onset. On the largest spatial scale we introduce a measure based on phase-locking intervals and wavelets into seizure modelling. It is used to resolve synchronization between different regions in the brain and identifies time-shifted scaling laws at different wavelet scales. We also compare our wavelet-based multiscale approach with maximum linear cross-correlation and mean-phase coherence measures

Characterization of K-Complexes and Slow Wave Activity in a Neural Mass Model

NREM sleep is characterized by two hallmarks, namely K-complexes (KCs) during sleep stage N2 and cortical slow oscillations (SOs) during sleep stage N3. While the underlying dynamics on the neuronal level is well known and can be easily measured, the resulting behavior on the macroscopic population level remains unclear. On the basis of an extended neural mass model of the cortex, we suggest a new interpretation of the mechanisms responsible for the generation of KCs and SOs. As the cortex transitions from wake to deep sleep, in our model it approaches an oscillatory regime via a Hopf bifurcation. Importantly, there is a canard phenomenon arising from a homoclinic bifurcation, whose orbit determines the shape of large amplitude SOs. A KC corresponds to a single excursion along the homoclinic orbit, while SOs are noise-driven oscillations around a stable focus. The model generates both time series and spectra that strikingly resemble real electroencephalogram data and points out possible differences between the different stages of natural sleep

Investigating the electrode-electrolyte interface modelling in cochlear implants

Objective. Proposing a good electrode-electrolyte interface (EEI) model and properly identifying relevant parameters may help designing safer and more optimized auditory nerve fiber stimulation and recording in cochlear implants (CI). However, in literature, EEI model parameter values exhibit large variability. We aim to explain some root causes of this variability using the Cole model and its simpler form, the Basic RC model. Approach. We use temporal and spectral methods and fit the models to stimulation pulse voltage response (SPVR) and electrochemical impedance spectroscopy (EIS) data. Main Results. Temporal fittings show that there are multiple sets of model parameters that provide a good fit to the SPVR data. Therefore, small methodological differences in literature may result in different model fits. While these models share similar characteristics at high frequencies >500 Hz, the SPVR fitting is blind to low frequencies, thus it cannot correctly estimate the Faradaic resistor. Similarly, the polarization capacitor and its fractional order are not estimated robustly (capacitor variations in the nano- to micro-farad range) due to limited observation of mid-range frequencies. EIS provides a good model fit down to & SIM;3Hz, and thus robust estimation for the polarization capacitor. At lower frequencies charge mechanisms may modify the EEI, requiring multi-compartment Cole model fitting to EIS to improve the estimation of Faradaic characteristics. Our EIS data measurements down to 0.05Hz show that a two-compartment Cole model is sufficient to explain the data. Significance. Our study describes the scope and limitation of SPVR and EIS fitting methods, by which literature variability is explained among CI EEI models. The estimation of mid-to-low-frequency characteristics of the CI EEI is not in the scope of the SPVR method. EIS provides a better fit; however, its results should not be extrapolated to unobserved frequencies where new charge transfer mechanisms may emerge at the EEI