Using reciprocity for relating the simulation of transcranial current stimulation to the EEG forward problem.

To explore the relationship between transcranial current stimulation (tCS) and the electroencephalography (EEG) forward problem, we investigate and compare accuracy and efficiency of a reciprocal and a direct EEG forward approach for dipolar primary current sources both based on the finite element method (FEM), namely the adjoint approach (AA) and the partial integration approach in conjunction with a transfer matrix concept (PI). By analyzing numerical results, comparing to analytically derived EEG forward potentials and estimating computational complexity in spherical shell models, AA turns out to be essentially identical to PI. It is then proven that AA and PI are also algebraically identical even for general head models. This relation offers a direct link between the EEG forward problem and tCS. We then demonstrate how the quasi-analytical EEG forward solutions in sphere models can be used to validate the numerical accuracies of FEM-based tCS simulation approaches. These approaches differ with respect to the ease with which they can be employed for realistic head modeling based on MRI-derived segmentations. We show that while the accuracy of the most easy to realize approach based on regular hexahedral elements is already quite high, it can be significantly improved if a geometry-adaptation of the elements is employed in conjunction with an isoparametric FEM approach. While the latter approach does not involve any additional difficulties for the user, it reaches the high accuracies of surface-segmentation based tetrahedral FEM, which is considerably more difficult to implement and topologically less flexible in practice. Finally, in a highly realistic head volume conductor model and when compared to the regular alternative, the geometry-adapted hexahedral FEM is shown to result in significant changes in tCS current flow orientation and magnitude up to 45° and a factor of 1.66, respectively.

The role of blood vessels in high-resolution volume conductor head modeling of EEG.

Reconstruction of the electrical sources of human EEG activity at high spatio-temporal accuracy is an important aim in neuroscience and neurological diagnostics. Over the last decades, numerous studies have demonstrated that realistic modeling of head anatomy improves the accuracy of source reconstruction of EEG signals. For example, including a cerebro-spinal fluid compartment and the anisotropy of white matter electrical conductivity were both shown to significantly reduce modeling errors. Here, we for the first time quantify the role of detailed reconstructions of the cerebral blood vessels in volume conductor head modeling for EEG. To study the role of the highly arborized cerebral blood vessels, we created a submillimeter head model based on ultra-high-field-strength (7T) structural MRI datasets. Blood vessels (arteries and emissary/intraosseous veins) were segmented using Frangi multi-scale vesselness filtering. The final head model consisted of a geometry-adapted cubic mesh with over 17×10(6) nodes. We solved the forward model using a finite-element-method (FEM) transfer matrix approach, which allowed reducing computation times substantially and quantified the importance of the blood vessel compartment by computing forward and inverse errors resulting from ignoring the blood vessels. Our results show that ignoring emissary veins piercing the skull leads to focal localization errors of approx. 5 to 15mm. Large errors (>2cm) were observed due to the carotid arteries and the dense arterial vasculature in areas such as in the insula or in the medial temporal lobe. Thus, in such predisposed areas, errors caused by neglecting blood vessels can reach similar magnitudes as those previously reported for neglecting white matter anisotropy, the CSF or the dura - structures which are generally considered important components of realistic EEG head models. Our findings thus imply that including a realistic blood vessel compartment in EEG head models will be helpful to improve the accuracy of EEG source analyses particularly when high accuracies in brain areas with dense vasculature are required.

The influence of sulcus width on simulated electric fields induced by transcranial magnetic stimulation.

Volume conduction models can help in acquiring knowledge about the distribution of the electric field induced by transcranial magnetic stimulation. One aspect of a detailed model is an accurate description of the cortical surface geometry. Since its estimation is difficult, it is important to know how accurate the geometry has to be represented. Previous studies only looked at the differences caused by neglecting the complete boundary between cerebrospinal fluid (CSF) and grey matter (Thielscher et al 2011 NeuroImage 54 234-43, Bijsterbosch et al 2012 Med. Biol. Eng. Comput. 50 671-81), or by resizing the whole brain (Wagner et al 2008 Exp. Brain Res. 186 539-50). However, due to the high conductive properties of the CSF, it can be expected that alterations in sulcus width can already have a significant effect on the distribution of the electric field. To answer this question, the sulcus width of a highly realistic head model, based on T1-, T2- and diffusion-weighted magnetic resonance images, was altered systematically. This study shows that alterations in the sulcus width do not cause large differences in the majority of the electric field values. However, considerable overestimation of sulcus width produces an overestimation of the calculated field strength, also at locations distant from the target location.

Complete electrode model in EEG: relationship and differences to the point electrode model.

In electroencephalography (EEG) source analysis, a primary current density generated by the neural activity of the brain is reconstructed from external electrode voltage measurements. This paper focuses on accurate and effective simulations of EEG through the complete electrode model (CEM). The CEM allows for the incorporation of the electrode size, shape and effective contact impedance into the forward simulation. Both neural currents in the brain and shunting currents between the electrodes and the skin can affect the measured voltages in the CEM. The goal of this study was to investigate the CEM by comparing it with the point electrode model (PEM), which is the current standard electrode model for EEG. We used a three-dimensional, realistic and high-resolution finite element head model as the reference computational domain in the comparison. The PEM could be formulated as a limit of the CEM, in which the effective impedance of each electrode goes to infinity and the size tends to zero. Numerical results concerning the forward and inverse errors and electrode voltage strengths with different impedances and electrode sizes are presented. Based on the results obtained, limits for extremely high and low impedance values of the shunting currents are suggested.