Review on biophysical modelling and simulation studies for transcranial magnetic stimulation

Stimulus effects of extremely low-frequency electric field exposure on calcium oscillations in a human cortical spheroid

High-intensity, low-frequency (1 Hz to 100 kHz) electric and magnetic fields (EF and MF) cause electrical excitation of the nervous system via an induced EF (iEF) in living tissue. However, the biological properties and thresholds of stimulus effects on synchronized activity in a three-dimensional (3D) neuronal network remain uncertain. In this study, we evaluated changes in neuronal network activity during extremely low-frequency EF (ELF-EF) exposure by measuring intracellular calcium ([Ca2+]i) oscillations, which reflect neuronal network activity. For ELF-EF exposure experiments, we used a human cortical spheroid (hCS), a 3D-cultured neuronal network generated from human induced pluripotent stem cell (hiPSC)-derived cortical neurons. A 50 Hz sinusoidal ELF-EF exposure modulated [Ca2+]i oscillations with dependencies on exposure intensity and duration. Based on the experimental setup and results, the iEF distribution inside the hCS was estimated using high-resolution numerical dosimetry. The numerical estimation revealed threshold values ranging between 255-510 V/m (peak) and 131-261 V/m (average). This indicates that thresholds of neuronal excitation in the hCS were equivalent to those of a thin nerve fiber.

Anatomical measurements and field modeling to assess transcranial magnetic stimulation motor and non-motor effects

OBJECTIVE

Explore how anatomical measurements and field modeling can be leveraged to improve investigations of transcranial magnetic stimulation (TMS) effects on both motor and non-motor TMS targets.

METHODS

TMS motor effects (targeting the primary motor cortex [M1]) were evaluated using the resting motor threshold (rMT), while TMS non-motor effects (targeting the superior temporal gyrus [STG]) were assessed using a pain memory task. Anatomical measurements included scalp-cortex distance (SCD) and cortical thickness (CT), whereas field modeling encompassed the magnitude of the electric field (E) induced by TMS.

RESULTS

Anatomical measurements and field modeling values differed significantly between M1 and STG. For TMS motor effects, rMT was correlated with SCD, CT, and E values at M1 (p < 0.05). No correlations were found between these metrics for the STG and TMS non-motor effects (pain memory; all p-values > 0.05).

CONCLUSION

Although anatomical measurements and field modeling are closely related to TMS motor effects, their relationship to non-motor effects - such as pain memory - appear to be much more tenuous and complex, highlighting the need for further advancement in our use of TMS and virtual lesion paradigms.

Dosimetry for repetitive transcranial magnetic stimulation: a translational study from Alzheimer's disease patients to controlledin vitroinvestigations

Objective.Recent studies have indicated that repetitive transcranial magnetic stimulation (rTMS) could enhance cognition in Alzheimer's Disease (AD) patients, but to now the molecular-level interaction mechanisms driving this effect remain poorly understood. While cognitive scores have been the primary measure of rTMS effectiveness, employing molecular-based approaches could offer more precise treatment predictions and prognoses. To reach this goal, it is fundamental to assess the electric field (E-field) and the induced current densities (J) within the stimulated brain areas and to translate these values toin vitrosystems specifically devoted in investigating molecular-based interactions of this stimulation.Approach.This paper offers a methodological procedure to guide dosimetric assessment to translate the E-field induced in humans (in a specific pilot study) intoin vitrosettings. Electromagnetic simulations on patients' head models and cellular holders were conducted to characterize exposure conditions and determine necessary adjustments forin vitroreplication of the same dose delivered in humans using the same stimulating coil.Main results.Our study highlighted the levels of E-field andJinduced in the target brain region and showed that the computed E-field andJwere different among patients that underwent the treatment, so to replicate the exposure to thein vitrosystem, we have to consider a range of electric quantities as reference. To match the E-field to the levels calculated in patients' brains, an increase of at least the 25% in the coil feeding current is necessary whenin vitrostimulations are performed. Conversely, to equalize current densities, modifications in the cells culture medium conductivity have to be implemented reducing it to one fifth of its value.Significance.This dosimetric assessment and subsequent experimental adjustments are essential to achieve controlledin vitroexperiments to better understand rTMS effects on AD cognition. Dosimetry is a fundamental step for comparing the cognitive effects with those obtained by stimulating a cellular model at an equal dose rigorously evaluated.

Prediction of hippocampal electric field in time series induced by TI-DMS with temporal convolutional network

Temporal interference deep-brain magnetic stimulation (TI-DMS) induces rhythmic electric field (EF) in the hippocampus to normalize cognitive function. The rhythmic time series of the hippocampal EF is essential for the assessment of TI-DMS. However, the finite element method (FEM) takes several hours to obtain the time series of EF. In order to reduce the time cost, the temporal convolutional network (TCN) model is adopted to predict the time series of hippocampal EF induced by TI-DMS. It takes coil configuration and loaded current as input and predicts the time series of maximum and mean values of the left and right hippocampal EF. The prediction takes only a few seconds. The model parameter combination of kernel size and layers is selected optimally by cross-validation method. The experimental results for multiple subjects show that the R2 of all the time series predicted by the model exceed 0.98. And the prediction accuracy is even higher as the input parameters approach the training set. These results demonstrate that the adopted model can quickly predict the time series of hippocampal EF induced by TI-DMS with relatively high accuracy, which is beneficial for future clinical applications.

The Effect of Stimulation Intensity, Sampling Frequency, and Sample Synchronization in TMS-EEG on the TMS Pulse Artifact Amplitude and Duration

Transcranial magnetic stimulation (TMS) coupled with electroencephalography (EEG) possesses diagnostic and therapeutic benefits. However, TMS provokes a large pulse artifact that momentarily obscures the cortical response, presenting a significant challenge for EEG data interpretation. We examined how stimulation intensity (SI), EEG sampling frequency (Fs) and synchronization of stimulation with EEG sampling influence the amplitude and duration of the pulse artifact. In eight healthy subjects, single-pulse TMS was administered to the primary motor cortex, due to its well-documented responsiveness to TMS. We applied two different SIs (90% and 120% of resting motor threshold, representing the commonly used subthreshold and suprathreshold levels) and Fs (conventional 5 kHz and high frequency 20 kHz) both with TMS synchronized with the EEG sampling and the conventional non-synchronized setting. Aside from removal of the DC-offset and epoching, no preprocessing was performed to the data. Using a random forest regression model, we identified that Fs had the largest impact on both the amplitude and duration of the pulse artifact, with median variable importance values of 1.444 and 1.327, respectively, followed by SI (0.964 and 1.083) and sampling synchronization (0.223 and 0.248). This indicated that Fs and SI are crucial for minimizing prediction error and thus play a pivotal role in accurately characterizing the pulse artifact. The results of this study enable focusing some of the study design parameters to minimize TMS pulse artifact, which is essential for both enhancing the reliability of clinical TMS-EEG applications and improving the overall integrity and interpretability of TMS-EEG data.

Comparison of Transcranial Magnetic Stimulation Dosimetry between Structured and Unstructured Grids Using Different Solvers

In recent years, the interest in transcranial magnetic stimulation (TMS) has surged, necessitating deeper understanding, development, and use of low-frequency (LF) numerical dosimetry for TMS studies. While various ad hoc dosimetric models exist, commercial software tools like SimNIBS v4.0 and Sim4Life v7.2.4 are preferred for their user-friendliness and versatility. SimNIBS utilizes unstructured tetrahedral mesh models, while Sim4Life employs voxel-based models on a structured grid, both evaluating induced electric fields using the finite element method (FEM) with different numerical solvers. Past studies primarily focused on uniform exposures and voxelized models, lacking realism. Our study compares these LF solvers across simplified and realistic anatomical models to assess their accuracy in evaluating induced electric fields. We examined three scenarios: a single-shell sphere, a sphere with an orthogonal slab, and a MRI-derived head model. The comparison revealed small discrepancies in induced electric fields, mainly in regions of low field intensity. Overall, the differences were contained (below 2% for spherical models and below 12% for the head model), showcasing the potential of computational tools in advancing exposure assessment required for TMS protocols in different bio-medical applications.

Efficacy of personalized rTMS to enhance upper limb function in subacute stroke patients: a protocol for a multi-center, randomized controlled study

Background

Repetitive transcranial magnetic stimulation (rTMS) is widely used therapy to enhance motor deficit in stroke patients. To date, rTMS protocols used in stroke patients are relatively unified. However, as the pathophysiology of stroke is diverse and individual functional deficits are distinctive, more precise application of rTMS is warranted. Therefore, the objective of this study was to determine the effects of personalized protocols of rTMS therapy based on the functional reserve of each stroke patient in subacute phase.

Methods

This study will recruit 120 patients with stroke in subacute phase suffering from the upper extremity motor impairment, from five different hospitals in Korea. The participants will be allocated into three different study conditions based on the functional reserve of each participant, measured by the results of TMS-induced motor evoked potentials (MEPs), and brain MRI with diffusion tensor imaging (DTI) evaluations. The participants of the intervention-group in the three study conditions will receive different protocols of rTMS intervention, a total of 10 sessions for 2 weeks: high-frequency rTMS on ipsilesional primary motor cortex (M1), high-frequency rTMS on ipsilesional ventral premotor cortex, and high-frequency rTMS on contralesional M1. The participants of the control-group in all three study conditions will receive the same rTMS protocol: low-frequency rTMS on contralesional M1. For outcome measures, the following assessments will be performed at baseline (T0), during-intervention (T1), post-intervention (T2), and follow-up (T3) periods: Fugl-Meyer Assessment (FMA), Box-and-block test, Action Research Arm Test, Jebsen-Taylor hand function test, hand grip strength, Functional Ambulatory Category, fractional anisotropy measured by the DTI, and brain network connectivity obtained from MRI. The primary outcome will be the difference of upper limb function, as measured by FMA from T0 to T2. The secondary outcomes will be the differences of other assessments.

Discussion

This study will determine the effects of applying different protocols of rTMS therapy based on the functional reserve of each patient. In addition, this methodology may prove to be more efficient than conventional rTMS protocols. Therefore, effective personalized application of rTMS to stroke patients can be achieved based on their severity, predicted mechanism of motor recovery, or functional reserves.

Clinical trial registration

https://clinicaltrials.gov/, identifier NCT06270238.

Individual Prediction of Electric Field Induced by Deep-Brain Magnetic Stimulation With CNN-Transformer

Deep-brain Magnetic Stimulation (DMS) can improve the symptoms caused by Alzheimer's disease by inducing rhythmic electric field in the deep brain, and the induced electric field is rhythm-dependent. However, calculating the induced electric field requires building a voxel model of the brain for the stimulated object, which usually takes several hours. In order to obtain the rhythm-dependent electric field induced by DMS in real time, we adopt a CNN-Transformer model to predict it. A data set with a sample size of 7350 is established for the training and testing of the model. 10-fold cross validation is used to determine the optimal hyperparameters for training CNN-Transformer. The combination of 5-layer CNN and 6-layer Transformer is verified as the optimal combination of CNN-Transformer model. The experimental results show that the CNN-Transformer model can complete the prediction in 0.731s (CPU) or 0.042s (GPU), and the overall performance metrics of prediction can reach: MAE =0.0269, RMSE =0.0420, MAPE =4.61% and R2=0.9627. The prediction performance of the CNN-Transformer model for the hippocampal electric field is better than that of the brain grey matter electric field, and the stimulation rhythm has less influence on the model performance than the coil configuration. Taking the same dataset to train and test the separate CNN model and Transformer model, it is found that CNN-Transformer has better prediction performance than the separate CNN model and Transformer model in the task of predicting electric field induced by DMS.

Predicting the hotspot location and motor threshold prior to transcranial magnetic stimulation using electric field modelling

Objective.To investigate whether the motor threshold (MT) and the location of the motor hotspot in transcranial magnetic stimulation (TMS) can be predicted with computational models of the induced electric field.Approach.Individualized computational models were constructed from structural magnetic resonance images of ten healthy participants, and the induced electric fields were determined with the finite element method. The models were used to optimize the location and direction of the TMS coil on the scalp to produce the largest electric field at a predetermined cortical target location. The models were also used to predict how the MT changes as the magnetic coil is moved to various locations over the scalp. To validate the model predictions, the motor evoked potentials were measured from the first dorsal interosseous (FDI) muscle with TMS in the ten participants. Both computational and experimental methods were preregistered prior to the experiments.Main results.Computationally optimized hotspot locations were nearly as accurate as those obtained using manual hotspot search procedures. The mean Euclidean distance between the predicted and the measured hotspot locations was approximately 1.3 cm with a 0.8 cm bias towards the anterior direction. Exploratory analyses showed that the bias could be removed by changing the cortical target location that was used for the prediction. The results also indicated a statistically significant relationship (p< 0.001) between the calculated electric field and the MT measured at several locations on the scalp.Significance.The results show that the individual TMS hotspot can be located using computational analysis without stimulating the subject or patient even once. Adapting computational modelling would save time and effort in research and clinical use of TMS.

A computational analysis of transcranial magnetic stimulation in patients with cranial defects and skull plate implants

We utilized computational analysis to investigate the impact of skull defects and skull implants on the TMS-induced EF. Our findings revealed a noteworthy alteration in the induced EF when acute skull defects were present. When high-conductivity titanium plates were used, we observed a pronounced increase in the peak EF, accompanied by a shift in the induced EF from the center towards both ends of the implant. These findings underscore the importance of carefully considering skull defects and implant materials during TMS.

A Review of Formulations, Boundary Value Problems and Solutions for Numerical Computation of Transcranial Magnetic Stimulation Fields

Since the inception of the transcranial magnetic stimulation (TMS) technique, it has become imperative to numerically compute the distribution of the electric field induced in the brain. Various models of the coil-brain system have been proposed for this purpose. These models yield a set of formulations and boundary conditions that can be employed to calculate the induced electric field. However, the literature on TMS simulation presents several of these formulations, leading to potential confusion regarding the interpretation and contribution of each source of electric field. The present study undertakes an extensive compilation of widely utilized formulations, boundary value problems and numerical solutions employed in TMS fields simulations, analyzing the advantages and disadvantages associated with each used formulation and numerical method. Additionally, it explores the implementation strategies employed for their numerical computation. Furthermore, this work provides numerical expressions that can be utilized for the numerical computation of TMS fields using the finite difference and finite element methods. Notably, some of these expressions are deduced within the present study. Finally, an overview of some of the most significant results obtained from numerical computation of TMS fields is presented. The aim of this work is to serve as a guide for future research endeavors concerning the numerical simulation of TMS.

Revisiting the focality of non-invasive brain stimulation - Implications for studies of human cognition

Non-invasive brain stimulation techniques are popular tools to investigate brain function in health and disease. Although transcranial magnetic stimulation (TMS) is widely used in cognitive neuroscience research to probe causal structure-function relationships, studies often yield inconclusive results. To improve the effectiveness of TMS studies, we argue that the cognitive neuroscience community needs to revise the stimulation focality principle - the spatial resolution with which TMS can differentially stimulate cortical regions. In the motor domain, TMS can differentiate between cortical muscle representations of adjacent fingers. However, this high degree of spatial specificity cannot be obtained in all cortical regions due to the influences of cortical folding patterns on the TMS-induced electric field. The region-dependent focality of TMS should be assessed a priori to estimate the experimental feasibility. Post-hoc simulations allow modeling of the relationship between cortical stimulation exposure and behavioral modulation by integrating data across stimulation sites or subjects.

Mapping Brain Motor Functions Using Transcranial Magnetic Stimulation with a Volume Conductor Model and Electrophysiological Experiments

Transcranial magnetic stimulation (TMS) activates brain cells in a noninvasive manner and can be used for mapping brain motor functions. However, the complexity of the brain anatomy prevents the determination of the exact location of the stimulated sites, resulting in the limitation of the spatial resolution of multiple targets. The aim of this study is to map two neighboring muscles in cortical motor areas accurately and quickly. Multiple stimuli were applied to the subject using a TMS stimulator to measure the motor-evoked potentials (MEPs) in the corresponding muscles. For each stimulation condition (coil location and angle), the induced electric field (EF) in the brain was computed using a volume conductor model for an individualized head model of the subject constructed from magnetic resonance images. A post-processing method was implemented to determine a TMS hotspot using EF corresponding to multiple stimuli, considering the amplitude of the measured MEPs. The dependence of the computationally estimated hotspot distribution on two target muscles was evaluated (n = 11). The center of gravity of the first dorsal interosseous cortical representation was lateral to the abductor digiti minimi by a minimum of 2 mm. The localizations were consistent with the putative sites obtained from previous EF-based studies and fMRI studies. The simultaneous cortical mapping of two finger muscles was achieved with only several stimuli, which is one or two orders of magnitude smaller than that in previous studies. Our proposal would be useful in the preoperative mapping of motor or speech areas to plan brain surgery interventions.

Computation of transcranial magnetic stimulation electric fields using self-supervised deep learning

Electric fields (E-fields) induced by transcranial magnetic stimulation (TMS) can be modeled using partial differential equations (PDEs). Using state-of-the-art finite-element methods (FEM), it often takes tens of seconds to solve the PDEs for computing a high-resolution E-field, hampering the wide application of the E-field modeling in practice and research. To improve the E-field modeling's computational efficiency, we developed a self-supervised deep learning (DL) method to compute precise TMS E-fields. Given a head model and the primary E-field generated by TMS coils, a DL model was built to generate a E-field by minimizing a loss function that measures how well the generated E-field fits the governing PDE. The DL model was trained in a self-supervised manner, which does not require any external supervision. We evaluated the DL model using both a simulated sphere head model and realistic head models of 125 individuals and compared the accuracy and computational speed of the DL model with a state-of-the-art FEM. In realistic head models, the DL model obtained accurate E-fields that were significantly correlated with the FEM solutions. The DL model could obtain precise E-fields within seconds for whole head models at a high spatial resolution, faster than the FEM. The DL model built for the simulated sphere head model also obtained an accurate E-field whose average difference from the analytical E-fields was 0.0054, comparable to the FEM solution. These results demonstrated that the self-supervised DL method could obtain precise E-fields comparable to the FEM solutions with improved computational speed.

Transcranial magnetic stimulation treatment in Alzheimer's disease: a meta-analysis of its efficacy as a function of protocol characteristics and degree of personalization

Alzheimer's disease (AD) represents the most common type of neurodegenerative disorder. Although our knowledge on the causes of AD remains limited and no curative treatments are available, several interventions have been proposed in trying to improve patients' symptomatology. Among those, transcranial magnetic stimulation (TMS) has been shown a promising, safe and noninvasive intervention to improve global cognitive functioning. Nevertheless, we currently lack agreement between research studies on the optimal stimulation protocol yielding the highest efficacy in these patients. To answer this query, we conducted a systematic literature search in PubMed, PsycINFO and Scopus databases and meta-analysis of studies published in the last 10 years (2010-2021) according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Differently from prior published meta-analytic work, we investigated whether protocols that considered participants-specific neuroimaging scans for the selection of individualized stimulation targets held more successful outcomes compared to those relying on a generalized targeting selection criteria. We then compared the effect sizes of subsets of studies based on additional protocol characteristics (frequency, duration of intervention, number of stimulation sites, use of concomitant cognitive training and patients' educational level). Our results confirm TMS efficacy in improving global cognitive functioning in mild-to-moderate AD patients, but also highlight the flaws of current protocols characteristics, including a possible lack of sufficient personalization in stimulation protocols.

Boosting psychological change: Combining non-invasive brain stimulation with psychotherapy

Mental health disorders and substance use disorders are a leading cause of morbidity and mortality worldwide, and one of the most important challenges for public health systems. While evidence-based psychotherapy is generally pursued to address mental health challenges, psychological change is often hampered by non-adherence to treatments, relapses, and practical barriers (e.g., time, cost). In recent decades, Non-invasive brain stimulation (NIBS) techniques have emerged as promising tools to directly target dysfunctional neural circuitry and promote long-lasting plastic changes. While the therapeutic efficacy of NIBS protocols for mental illnesses has been established, neuromodulatory interventions might also be employed to support the processes activated by psychotherapy. Indeed, combining psychotherapy with NIBS might help tailor the treatment to the patient's unique characteristics and therapeutic goal, and would allow more direct control of the neuronal changes induced by therapy. Herein, we overview emerging evidence on the use of NIBS to enhance the psychotherapeutic effect, while highlighting the next steps in advancing clinical and research methods toward personalized intervention approaches.

The Relation between Induced Electric Field and TMS-Evoked Potentials: A Deep TMS-EEG Study

Transcranial magnetic stimulation (TMS) in humans induces electric fields (E-fields, EF) that perturb and modulate the brain’s endogenous neuronal activity and result in the generation of TMS-evoked potentials (TEPs). The exact relation of the characteristics of the induced E-field and the intensity of the brains’ response, as measured by electroencephalography (EEG), is presently unclear. In this pilot study, conducted on three healthy subjects and two patients with generalized epilepsy (total: 3 males, 2 females, mean age of 26 years; healthy: 2 males, 1 female, mean age of 25.7 years; patients: 1 male, 1 female, mean age of 26.5 years), we investigated the temporal and spatial relations of the E-field, induced by single-pulse stimuli, and the brain’s response to TMS. Brain stimulation was performed with a deep TMS device (BrainsWay Ltd., Jerusalem, Israel) and an H7 coil placed over the central area. The induced EF was computed on personalized anatomical models of the subjects through magneto quasi-static simulations. We identified specific time instances and brain regions that exhibit high positive or negative associations of the E-field with brain activity. In addition, we identified significant correlations of the brain’s response intensity with the strength of the induced E-field and finally prove that TEPs are better correlated with E-field characteristics than with the stimulator’s output. These observations provide further insight in the relation between E-field and the ensuing cortical activation, validate in a clinically relevant manner the results of E-field modeling and reinforce the view that personalized approaches should be adopted in the field of non-invasive brain stimulation.

From the Group to the Individual in Schizophrenia Spectrum Disorders: Biomarkers of Social Cognitive Impairments and Therapeutic Translation

People with schizophrenia spectrum disorders (SSDs) often experience persistent social cognitive impairments, associated with poor functional outcome. There are currently no approved treatment options for these debilitating symptoms, highlighting the need for novel therapeutic strategies. Work to date has elucidated differential social processes and underlying neural circuitry affected in SSDs, which may be amenable to modulation using neurostimulation. Further, advances in functional connectivity mapping and electric field modeling may be used to identify individualized treatment targets to maximize the impact of brain stimulation on social cognitive networks. Here, we review literature supporting a roadmap for translating functional connectivity biomarker discovery to individualized treatment development for social cognitive impairments in SSDs. First, we outline the relevance of social cognitive impairments in SSDs. We review machine learning approaches for dimensional brain-behavior biomarker discovery, emphasizing the importance of individual differences. We synthesize research showing that brain stimulation techniques, such as repetitive transcranial magnetic stimulation, can be used to target relevant networks. Further, functional connectivity-based individualized targeting may enhance treatment response. We then outline recent approaches to account for neuroanatomical variability and optimize coil positioning to individually maximize target engagement. Overall, the synthesized literature provides support for the utility and feasibility of this translational approach to precision treatment. The proposed roadmap to translate biomarkers of social cognitive impairments to individualized treatment is currently under evaluation in precision-guided trials. Such a translational approach may also be applicable across conditions and generalizable for the development of individualized neurostimulation targeting other behavioral deficits.

Multi-scale neural decoding and analysis

Objective.

Complex spatiotemporal neural activity encodes rich information related to behavior and cognition. Conventional research has focused on neural activity acquired using one of many different measurement modalities, each of which provides useful but incomplete assessment of the neural code. Multi-modal techniques can overcome tradeoffs in the spatial and temporal resolution of a single modality to reveal deeper and more comprehensive understanding of system-level neural mechanisms. Uncovering multi-scale dynamics is essential for a mechanistic understanding of brain function and for harnessing neuroscientific insights to develop more effective clinical treatment.

Approach.

We discuss conventional methodologies used for characterizing neural activity at different scales and review contemporary examples of how these approaches have been combined. Then we present our case for integrating activity across multiple scales to benefit from the combined strengths of each approach and elucidate a more holistic understanding of neural processes.

Main results.

We examine various combinations of neural activity at different scales and analytical techniques that can be used to integrate or illuminate information across scales, as well the technologies that enable such exciting studies. We conclude with challenges facing future multi-scale studies, and a discussion of the power and potential of these approaches.

Significance.

This roadmap will lead the readers toward a broad range of multi-scale neural decoding techniques and their benefits over single-modality analyses. This Review article highlights the importance of multi-scale analyses for systematically interrogating complex spatiotemporal mechanisms underlying cognition and behavior.

Computational exploration of epidural cortical stimulation using a realistic head model

Motor cortex stimulation, either non-invasively or with implanted electrodes, has been applied worldwide as a treatment for intractable neuropathic pain syndromes. Although computer simulations of non-invasive brain stimulation have been investigated largely to optimize protocols and improve our understanding of underlying mechanisms using a realistic head model, computational studies of invasive cortical stimulation are rare and limited to very simplified cortical models. In this paper, we present an anatomically realistic head model for epidural cortical stimulation that includes the most sophisticated epidural electrodes with an insulating paddle. The head model predicted the stimulus-induced field strengths according to two different stimulation techniques, bipolar and monopolar stimulations. We found that the stimulus-induced field focused on the precentral and postcentral gyri because of the epidural lead's invasiveness. Different stimulation configurations influenced the shape of the field markedly, and complex patterns of inward and outward directions of the radial field were observed in bipolar stimulation compared to those in monopolar stimulation. The spatial distributions of field strength showed that the optimal stimulation varied according to the target areas. In conclusion, we proposed an anatomically realistic head model and a sophisticated epidural lead to simulate epidural cortical stimulation-induced field strengths and identified the importance of such detailed modeling for epidural cortical stimulation because of the current's shunting through the cerebrospinal fluid.

Influence of segmentation accuracy in structural MR head scans on electric field computation for TMS and tES

In several diagnosis and therapy procedures based on electrostimulation effect, the internal physical quantity related to the stimulation is the induced electric field. To estimate the induced electric field in an individual human model, the segmentation of anatomical imaging, such as magnetic resonance image (MRI) scans, of the corresponding body parts into tissues is required. Then, electrical properties associated with different annotated tissues are assigned to the digital model to generate a volume conductor. However, the segmentation of different tissues is a tedious task with several associated challenges specially with tissues appear in limited regions and/or low-contrast in anatomical images. An open question is how segmentation accuracy of different tissues would influence the distribution of the induced electric field. In this study, we applied parametric segmentation of different tissues to exploit the segmentation of available MRI to generate different quality of head models using deep learning neural network architecture, named ForkNet. Then, the induced electric field are compared to assess the effect of model segmentation variations. Computational results indicate that the influence of segmentation error is tissue-dependent. In brain, sensitivity to segmentation accuracy is relatively high in cerebrospinal fluid (CSF), moderate in gray matter (GM) and low in white matter for transcranial magnetic stimulation (TMS) and transcranial electrical stimulation (tES). A CSF segmentation accuracy reduction of 10% in terms of Dice coefficient (DC) lead to decrease up to 4% in normalized induced electric field in both applications. However, a GM segmentation accuracy reduction of 5.6% DC leads to increase of normalized induced electric field up to 6%. Opposite trend of electric field variation was found between CSF and GM for both TMS and tES. The finding obtained here would be useful to quantify potential uncertainty of computational results.