Towards a mesoscale physical modeling framework for stereotactic-EEG recordings.

Objective. Stereotactic-electroencephalography (SEEG) and scalp EEG recordings can be modeled using mesoscale neural mass population models (NMMs). However, the relationship between those mathematical models and the physics of the measurements is unclear. In addition, it is challenging to represent SEEG data by combining NMMs and volume conductor models due to the intermediate spatial scale represented by these measurements. Approach. We provide a framework combining the multi-compartmental modeling formalism and a detailed geometrical model to simulate the transmembrane currents that appear in layer 3, 5 and 6 pyramidal cells due to a synaptic input. With this approach, it is possible to realistically simulate the current source density (CSD) depth profile inside a cortical patch due to inputs localized into a single cortical layer and the induced voltage measured by two SEEG contacts using a volume conductor model. Based on this approach, we built a framework to connect the activity of a NMM with a volume conductor model and we simulated an example of SEEG signal as a proof of concept. Main results. CSD depends strongly on the distribution of the synaptic inputs onto the different cortical layers and the equivalent current dipole strengths display substantial differences (of up to a factor of four in magnitude in our example). Thus, the inputs coming from different neural populations do not contribute equally to the electrophysiological recordings. A direct consequence of this is that the raw output of NMMs is not a good proxy for electrical recordings. We also show that the simplest CSD model that can accurately reproduce SEEG measurements can be constructed from discrete monopolar sources (one per cortical layer). Significance. Our results highlight the importance of including a physical model in NMMs to represent measurements. We provide a framework connecting microscale neuron models with the neural mass formalism and with physical models of the measurement process that can improve the accuracy of predicted electrophysiological recordings.