Despite the widespread use of transcranial magnetic stimulation (TMS), the precise cortical locations underlying the resulting physiological and behavioral effects are still only coarsely known. To date, mapping strategies have relied on projection approaches (often termed "center of gravity" approaches) or maximum electric field value evaluation, and therefore localize the stimulated cortical site only approximately and indirectly. Focusing on the motor cortex, we present and validate a novel method to reliably determine the effectively stimulated cortical site at the individual subject level. The approach combines measurements of motor evoked potentials (MEPs) at different coil positions and orientations with numerical modeling of induced electric fields. We identify sharply bounded cortical areas, around the gyral crowns and rims of the motor hand area, as the origin of MEPs and show that the magnitude of the tangential component and the overall magnitude of the electric field are most relevant for the observed effect. To validate our approach, we identified the coil location and orientation that produces the maximal electric field at the predicted stimulation site, and then experimentally show that this location produces MEPs more efficiently than other tested locations/orientations. Moreover, we used extensive uncertainty and sensitivity analyses to verify the robustness of the method and identify the most critical model parameters. Our generic approach improves the localization of the cortical area effectively stimulated by TMS and may be transferred to other modalities such as language mapping.

Uncertainty surrounding ohmic tissue conductivity impedes accurate calculation of the electric fields generated by non-invasive brain stimulation. We present an efficient and generic technique for uncertainty and sensitivity analyses, which quantifies the reliability of field estimates and identifies the most influential parameters. For this purpose, we employ a non-intrusive generalized polynomial chaos expansion to compactly approximate the multidimensional dependency of the field on the conductivities. We demonstrate that the proposed pipeline yields detailed insight into the uncertainty of field estimates for transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS), identifies the most relevant tissue conductivities, and highlights characteristic differences between stimulation methods. Specifically, we test the influence of conductivity variations on (i) the magnitude of the electric field generated at each gray matter location, (ii) its normal component relative to the cortical sheet, (iii) its overall magnitude (indexed by the 98th percentile), and (iv) its overall spatial distribution. We show that TMS fields are generally less affected by conductivity variations than tDCS fields. For both TMS and tDCS, conductivity uncertainty causes much higher uncertainty in the magnitude as compared to the direction and overall spatial distribution of the electric field. Whereas the TMS fields were predominantly influenced by gray and white matter conductivity, the tDCS fields were additionally dependent on skull and scalp conductivities. Comprehensive uncertainty analyses of complex systems achieved by the proposed technique are not possible with classical methods, such as Monte Carlo sampling, without extreme computational effort. In addition, our method has the advantages of directly yielding interpretable and intuitive output metrics and of being easily adaptable to new problems.

Conceptual knowledge is central to human cognition. The left posterior inferior parietal lobe (pIPL) is implicated by neuroimaging studies as a multimodal hub representing conceptual knowledge related to various perceptual-motor modalities. However, the causal role of left pIPL in conceptual processing remains unclear. Here, we transiently disrupted left pIPL function with transcranial magnetic stimulation (TMS) to probe its causal relevance for the retrieval of action and sound knowledge. We compared effective TMS over left pIPL with sham TMS, while healthy participants performed three different tasks-lexical decision, action judgment, and sound judgment-on words with a high or low association to actions and sounds. We found that pIPL-TMS selectively impaired action judgments on low sound-low action words. For the first time, we directly related computational simulations of the TMS-induced electrical field to behavioral performance, which revealed that stronger stimulation of left pIPL is associated with worse performance for action but not sound judgments. These results indicate that left pIPL causally supports conceptual processing when action knowledge is task-relevant and cannot be compensated by sound knowledge. Our findings suggest that left pIPL is specialized for the retrieval of action knowledge, challenging the view of left pIPL as a multimodal conceptual hub.

When modeling transcranial electrical stimulation (TES) and transcranial magnetic stimulation (TMS) in the brain, the meninges - dura, arachnoid, and pia mater - are often neglected due to high computational costs. We investigate the impact of the meningeal layers on the cortical electric field in TES and TMS while considering the headreco segmentation as the base model. We use T1/T2 MRI data from 16 subjects and apply the boundary element fast multipole method with adaptive mesh refinement, which enables us to accurately solve this problem and establish method convergence at reasonable computational cost. We compare electric fields in the presence and absence of various meninges for two brain areas (M1_HAND and DLPFC) and for several distinct TES and TMS setups. Maximum electric fields in the cortex for focal TES consistently increase by approximately 30% on average when the meninges are present in the CSF volume. Their effect on the maximum field can be emulated by reducing the CSF conductivity from 1.65 S/m to approximately 0.85 S/m. In stark contrast to that, the TMS electric fields in the cortex are only weakly affected by the meningeal layers and slightly (∼6%) decrease on average when the meninges are included. Our results quantify the influence of the meninges on the cortical TES and TMS electric fields. Both focal TES and TMS results are very consistent. The focal TES results are also in a good agreement with a prior relevant study. The solver and the mesh generator for the meningeal layers (compatible with SimNIBS) are available online.

We derived computationally efficient average response models of different types of cortical neurons, which are subject to external electric fields from Transcranial Magnetic Stimulation. We used 24 reconstructions of pyramidal cells (PC) from layer 2/3, 245 small, nested, and large basket cells from layer 4, and 30 PC from layer 5 with different morphologies for deriving average models. With these models, it is possible to efficiently estimate the stimulation thresholds depending on the underlying electric field distribution in the brain, without having to implement and compute complex neuron compartment models. The stimulation thresholds were determined by exposing the neurons to TMS-induced electric fields with different angles, intensities, pulse waveforms, and field decays along the somato-dendritic axis. The derived average response models were verified by reference simulations using a high-resolution realistic head model containing several million neurons. The relative errors of the estimated thresholds between the average model and the reference model ranged between -3% and 3.7% in 98% of the cases, while the computation time was only a fraction of a second compared to several weeks. Finally, we compared the model behavior to TMS experiments and observed high correspondence to the orientation sensitivity of motor evoked potentials. The derived models were compared to the classical cortical column cosine model and to simplified ball-and-stick neurons. It was shown that both models oversimplify the complex interplay between the electric field and the neurons and do not adequately represent the directional sensitivity of the different cell types. The derived models are simple to apply and only require the TMS-induced electric field in the brain as input variable. The models and code are available to the general public in open-source repositories for integration into TMS studies to estimate the expected stimulation thresholds for an improved dosing and treatment planning in the future.