Causal functional maps of brain rhythms in working memory.

Human working memory is a key cognitive process that engages multiple functional anatomical nodes across the brain. Despite a plethora of correlative neuroimaging evidence regarding the working memory architecture, our understanding of critical hubs causally controlling overall performance is incomplete. Causal interpretation requires cognitive testing following safe, temporal, and controllable neuromodulation of specific functional anatomical nodes. Such experiments became available in healthy humans with the advance of transcranial alternating current stimulation (tACS). Here, we synthesize findings of 28 placebo-controlled studies (in total, 1,057 participants) that applied frequency-specific noninvasive stimulation of neural oscillations and examined working memory performance in neurotypical adults. We use a computational meta-modeling method to simulate each intervention in realistic virtual brains and test reported behavioral outcomes against the stimulation-induced electric fields in different brain nodes. Our results show that stimulating anterior frontal and medial temporal theta oscillations and occipitoparietal gamma rhythms leads to significant dose-dependent improvement in working memory task performance. Conversely, prefrontal gamma modulation is detrimental to performance. Moreover, we found distinct spatial expression of theta subbands, where working memory changes followed orbitofrontal high-theta modulation and medial temporal low-theta modulation. Finally, all these results are driven by changes in working memory accuracy rather than processing time measures. These findings provide a fresh view of the working memory mechanisms, complementary to neuroimaging research, and propose hypothesis-driven targets for the clinical treatment of working memory deficits.

Dissociation of Centrally and Peripherally Induced Transcranial Magnetic Stimulation Effects in Nonhuman Primates.

Transcranial magnetic stimulation (TMS) is a noninvasive brain stimulation method that is rapidly growing in popularity for studying causal brain-behavior relationships. However, its dose-dependent centrally induced neural mechanisms and peripherally induced sensory costimulation effects remain debated. Understanding how TMS stimulation parameters affect brain responses is vital for the rational design of TMS protocols. Studying these mechanisms in humans is challenging because of the limited spatiotemporal resolution of available noninvasive neuroimaging methods. Here, we leverage invasive recordings of local field potentials in a male and a female nonhuman primate (rhesus macaque) to study TMS mesoscale responses. We demonstrate that early TMS-evoked potentials show a sigmoidal dose-response curve with stimulation intensity. We further show that stimulation responses are spatially specific. We use several control conditions to dissociate centrally induced neural responses from auditory and somatosensory coactivation. These results provide crucial evidence regarding TMS neural effects at the brain circuit level. Our findings are highly relevant for interpreting human TMS studies and biomarker developments for TMS target engagement in clinical applications.SIGNIFICANCE STATEMENT Transcranial magnetic stimulation (TMS) is a widely used noninvasive brain stimulation method to stimulate the human brain. To advance its utility for clinical applications, a clear understanding of its underlying physiological mechanisms is crucial. Here, we perform invasive electrophysiological recordings in the nonhuman primate brain during TMS, achieving a spatiotemporal precision not available in human EEG experiments. We find that evoked potentials are dose dependent and spatially specific, and can be separated from peripheral stimulation effects. This means that TMS-evoked responses can indicate a direct physiological stimulation response. Our work has important implications for the interpretation of human TMS-EEG recordings and biomarker development.

Anatomical details affect electric field predictions for non-invasive brain stimulation in non-human primates.

Non-human primates (NHPs) have become key for translational research in noninvasive brain stimulation (NIBS). However, in order to create comparable stimulation conditions for humans it is vital to study the accuracy of current modeling practices across species. Numerical models to simulate electric fields are an important tool for experimental planning in NHPs and translation to human studies. It is thus essential whether and to what extent the anatomical details of NHP models agree with current modeling practices when calculating NIBS electric fields. Here, we create highly accurate head models of two non-human primates (NHP) MR data. We evaluate how muscle tissue and head field of view (depending on MRI parameters) affect simulation results in transcranial electric and magnetic stimulation (TES and TMS). Our findings indicate that the inclusion of anisotropic muscle can affect TES electric field strength up to 22% while TMS is largely unaffected. Additionally, comparing a full head model to a cropped head model illustrates the impact of head field of view on electric fields for both TES and TMS. We find opposing effects between TES and TMS with an increase up to 24.8% for TES and a decrease up to 24.6% for TMS for the cropped head model compared to the full head model. Our results provide important insights into the level of anatomical detail needed for NHP head models and can inform future translational efforts for NIBS studies.

Are we really targeting and stimulating DLPFC by placing transcranial electrical stimulation (tES) electrodes over F3/F4?

In many clinical trials involving transcranial electrical stimulation (tES), target electrodes are typically placed over DLPFC with the assumption that this will primarily stimulate the underlying brain region. However, our study aimed to evaluate the electric fields (EF) that are actually delivered and identify prefrontal regions that may be inadvertently targeted in DLPFC tES. Head models were generated from the Human Connectome Project database's T1 + T2-weighted MRIs of 80 healthy adults. Two common DLPFC montages were simulated; symmetric-F4/F3, and asymmetric-F4/Fp1. Averaged EF was extracted from (1) the center of the target electrode (F4), and (2) the top 1% of voxels showing the strongest EF in individualized EF maps. Interindividual variabilities were quantified with the standard deviation of EF peak location/value. Similar steps were repeated with 66 participants with methamphetamine use disorder (MUDs) as an independent clinical population. In healthy adults, the group-level location of EF peaks was situated in the medial-frontopolar, and the individualized EF peaks were positioned in a cube with a volume of 29 cm3 /46 cm3 (symmetric/asymmetric montages). EFs in the frontopolar area were significantly higher than EF "under" the target electrode in both symmetric (peak: 0.41 ± 0.06, F4:0.22 ± 0.04) and asymmetric (peak: 0.38 ± 0.04, F4:0.2 ± 0.04) montages (Heges'g > 0.7). Similar results with slight between-group differences were found in MUDs. We highlighted that in common DLPFC tES montages, in addition to interindividual/intergroup variability, the frontopolar received the highest EFs rather than DLPFC as the main target. We specifically recommended considering the potential involvement of the frontopolar area as a mechanism underlying the effectiveness of DLPFC tES protocols.

Effects of transcranial alternating current stimulation on spiking activity in computational models of single neocortical neurons.

Neural oscillations are a key mechanism for information transfer in brain circuits. Rhythmic fluctuations of local field potentials control spike timing through cyclic membrane de- and hyperpolarization. Transcranial alternating current stimulation (tACS) is a non-invasive neuromodulation method which can directly interact with brain oscillatory activity by imposing an oscillating electric field on neurons. Despite its increasing use, the basic mechanisms of tACS are still not fully understood. Here, we investigate in a computational study the effects of tACS on morphologically realistic neurons with ongoing spiking activity. We characterize the membrane polarization as a function of electric field strength and subsequent effects on spiking activity in a set of 25 neurons from different neocortical layers. We find that tACS does not affect the firing rate of investigated neurons for electric field strengths applicable to human studies. However, we find that the applied electric fields entrain the spiking activity of large pyramidal neurons and large basket neurons at < 1 mV/mm field strengths. Our model results are in line with recent experimental studies and can provide a mechanistic framework to understand the effects of oscillating electric fields on single neuron activity. They highlight the importance of neuron morphology and biophysics in responsiveness to electrical stimulation.

In situ electrochemical observation of anisotropic lattice contraction of La0.6Sr0.4FeO3-δ electrodes during pulsed laser deposition.

La0.6Sr0.4FeO3-δ (LSF) electrodes were grown on different electrolyte substrates by pulsed laser deposition (PLD) and their oxygen exchange reaction (OER) resistance was tracked in real-time by in situ PLD impedance spectroscopy (i-PLD) inside the PLD chamber. This enables measurements on pristine surfaces free from any contaminations and the direct observation of thickness dependent properties. As substrates, yttria-stabilized zirconia single crystals (YSZ) were used for polycrystalline LSF growth and La0.95Sr0.05Ga0.95Mg0.05O3-δ (LSGM) single crystals or YSZ single crystals with a 5 nm buffer-layer of Gd0.2Ce0.8O2-δ for epitaxial LSF film growth. While polycrystalline LSF electrodes show a constant OER resistance in a broad thickness range, epitaxially grown LSF electrodes exhibit a continuous and strong increase of the OER resistance with film thickness until ≈60 nm. In addition, the activation energy of the OER resistance increases by 0.23 eV compared to polycrystalline LSF. High resolution transmission electron microscopy (HRTEM) and X-ray diffraction (XRD) measurements reveal an increasing contraction of the out-of-plane lattice parameter in the epitaxial LSF electrodes over electrode thickness. Defect thermodynamic simulations suggest that the decrease of the LSF unit cell volume is accompanied by a lowering of the oxygen vacancy concentration, explaining both the resistive increase and the increased activation energy.

Increasing propensity to mind-wander by transcranial direct current stimulation? A registered report.

Transcranial direct current stimulation (tDCS) has been proposed to be able to modulate different cognitive functions. However, recent meta-analyses conclude that its efficacy is still in question. Recently, an increase in subjects' propensity to mind-wander has been reported as a consequence of anodal stimulation of the left dorsolateral prefrontal cortex (Axelrod et al., Proceedings of the National Academy of Sciences of the United States of America, 112, 2015). In addition, an independent group found a decrease in mind wandering after cathodal stimulation of the same region. These findings seem to indicate that high-level cognitive processes such as mind wandering can reliably be influenced by non-invasive brain stimulation. However, these previous studies used low sample sizes and are as such subject to concerns regarding the replicability of their findings. In this registered report, we implement a high-powered replication of Axelrod et al. (2015) finding that mind-wandering propensity can be increased by anodal tDCS. We used Bayesian statistics and a preregistered sequential-sampling design resulting in a total sample size of N = 192 participants collected across three different laboratories. Our findings show support against a stimulation effect on self-reported mind-wandering scores. The effect was small, in the opposite direction as predicted and not reliably different from zero. Using a Bayes Factor specifically designed to test for replication success, we found strong evidence against a successful replication of the original study. Finally, even when combining data from both the original and replication studies, we could not find evidence for an effect of anodal stimulation. Our results underline the importance of designing studies with sufficient power to detect evidence for or against behavioural effects of non-invasive brain stimulation techniques, preferentially using robust Bayesian statistics in preregistered reports.

Traumatic axonal injury influences the cognitive effect of non-invasive brain stimulation.

Non-invasive brain stimulation has been widely investigated as a potential treatment for a range of neurological and psychiatric conditions, including brain injury. However, the behavioural effects of brain stimulation are variable, for reasons that are poorly understood. This is a particular challenge for traumatic brain injury, where patterns of damage and their clinical effects are heterogeneous. Here we test the hypothesis that the response to transcranial direct current stimulation following traumatic brain injury is dependent on white matter damage within the stimulated network. We used a novel simultaneous stimulation-MRI protocol applying anodal, cathodal and sham stimulation to 24 healthy control subjects and 35 patients with moderate/severe traumatic brain injury. Stimulation was applied to the right inferior frontal gyrus/anterior insula node of the salience network, which was targeted because our previous work had shown its importance to executive function. Stimulation was applied during performance of the Stop Signal Task, which assesses response inhibition, a key component of executive function. Structural MRI was used to assess the extent of brain injury, including diffusion MRI assessment of post-traumatic axonal injury. Functional MRI, which was simultaneously acquired to delivery of stimulation, assessed the effects of stimulation on cognitive network function. Anodal stimulation improved response inhibition in control participants, an effect that was not observed in the patient group. The extent of traumatic axonal injury within the salience network strongly influenced the behavioural response to stimulation. Increasing damage to the tract connecting the stimulated right inferior frontal gyrus/anterior insula to the rest of the salience network was associated with reduced beneficial effects of stimulation. In addition, anodal stimulation normalized default mode network activation in patients with poor response inhibition, suggesting that stimulation modulates communication between the networks involved in supporting cognitive control. These results demonstrate an important principle: that white matter structure of the connections within a stimulated brain network influences the behavioural response to stimulation. This suggests that a personalized approach to non-invasive brain stimulation is likely to be necessary, with structural integrity of the targeted brain networks an important criterion for patient selection and an individualized approach to the selection of stimulation parameters.

Limitations of ex vivo measurements for in vivo neuroscience.

A long history of postmortem studies has provided significant insight into human brain structure and organization. Cadavers have also proven instrumental for the measurement of artifacts and nonneural effects in functional imaging, and more recently, the study of biophysical properties critical to brain stimulation. However, death produces significant changes in the biophysical properties of brain tissues, making an ex vivo to in vivo comparison complex, and even questionable. This study directly compares biophysical properties of electric fields arising from transcranial electric stimulation (TES) in a nonhuman primate brain pre- and postmortem. We show that pre- vs. postmortem, TES-induced intracranial electric fields differ significantly in both strength and frequency response dynamics, even while controlling for confounding factors such as body temperature. Our results clearly indicate that ex vivo cadaver and in vivo measurements are not easily equitable. In vivo examinations remain essential to establishing an adequate understanding of even basic biophysical phenomena in vivo.

Comparative modeling of transcranial magnetic and electric stimulation in mouse, monkey, and human.

Transcranial magnetic stimulation (TMS) and transcranial electric stimulation (TES) are increasingly popular methods to noninvasively affect brain activity. However, their mechanism of action and dose-response characteristics remain under active investigation. Translational studies in animals play a pivotal role in these efforts due to a larger neuroscientific toolset enabled by invasive recordings. In order to translate knowledge gained in animal studies to humans, it is crucial to generate comparable stimulation conditions with respect to the induced electric field in the brain. Here, we conduct a finite element method (FEM) modeling study of TMS and TES electric fields in a mouse, capuchin and macaque monkeys, and a human model. We systematically evaluate the induced electric fields and analyze their relationship to head and brain anatomy. We find that with increasing head size, TMS-induced electric field strength first increases and then decreases according to a two-term exponential function. TES-induced electric field strength strongly decreases from smaller to larger specimen with up to 100x fold differences across species. Our results can serve as a basis to compare and match stimulation parameters across studies in animals and humans.

Cognitive enhancement with Salience Network electrical stimulation is influenced by network structural connectivity.

The Salience Network (SN) and its interactions are important for cognitive control. We have previously shown that structural damage to the SN is associated with abnormal functional connectivity between the SN and Default Mode Network (DMN), abnormal DMN deactivation, and impaired response inhibition, which is an important aspect of cognitive control. This suggests that stimulating the SN might enhance cognitive control. Here, we tested whether non-invasive transcranial direct current stimulation (TDCS) could be used to modulate activity within the SN and enhance cognitive control. TDCS was applied to the right inferior frontal gyrus/anterior insula cortex during performance of the Stop Signal Task (SST) and concurrent functional (f)MRI. Anodal TDCS improved response inhibition. Furthermore, stratification of participants based on SN structural connectivity showed that it was an important influence on both behavioural and physiological responses to anodal TDCS. Participants with high fractional anisotropy within the SN showed improved SST performance and increased activation of the SN with anodal TDCS, whilst those with low fractional anisotropy within the SN did not. Cathodal stimulation of the SN produced activation of the right caudate, an effect which was not modulated by SN structural connectivity. Our results show that stimulation targeted to the SN can improve response inhibition, supporting the causal influence of this network on cognitive control and confirming it as a target to produce cognitive enhancement. Our results also highlight the importance of structural connectivity as a modulator of network to TDCS, which should guide the design and interpretation of future stimulation studies.

Blinding is compromised for transcranial direct current stimulation at 1 mA for 20 min in young healthy adults.

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation method that is frequently used to study cortical excitability changes and their impact on cognitive functions in humans. While most stimulators are capable of operating in double-blind mode, the amount of discomfort experienced during tDCS may break blinding. Therefore, specifically designed sham stimulation protocols are being used. The "fade-in, short-stimulation, fade-out" (FSF) protocol has been used in hundreds of studies and is commonly believed to be indistinguishable from real stimulation applied at 1 mA for 20 min. We analysed subjective reports of 192 volunteers, who either received real tDCS (n = 96) or FSF tDCS (n = 96). Participants reported more discomfort for real tDCS and correctly guessed the condition above chance-level. These findings indicate that FSF does not ensure complete blinding and that better active sham protocols are needed.

Brain state and polarity dependent modulation of brain networks by transcranial direct current stimulation.

Despite its widespread use in cognitive studies, there is still limited understanding of whether and how transcranial direct current stimulation (tDCS) modulates brain network function. To clarify its physiological effects, we assessed brain network function using functional magnetic resonance imaging (fMRI) simultaneously acquired during tDCS stimulation. Cognitive state was manipulated by having subjects perform a Choice Reaction Task or being at "rest." A novel factorial design was used to assess the effects of brain state and polarity. Anodal and cathodal tDCS were applied to the right inferior frontal gyrus (rIFG), a region involved in controlling activity large-scale intrinsic connectivity networks during switches of cognitive state. tDCS produced widespread modulation of brain activity in a polarity and brain state dependent manner. In the absence of task, the main effect of tDCS was to accentuate default mode network (DMN) activation and salience network (SN) deactivation. In contrast, during task performance, tDCS increased SN activation. In the absence of task, the main effect of anodal tDCS was more pronounced, whereas cathodal tDCS had a greater effect during task performance. Cathodal tDCS also accentuated the within-DMN connectivity associated with task performance. There were minimal main effects of stimulation on network connectivity. These results demonstrate that rIFG tDCS can modulate the activity and functional connectivity of large-scale brain networks involved in cognitive function, in a brain state and polarity dependent manner. This study provides an important insight into mechanisms by which tDCS may modulate cognitive function, and also has implications for the design of future stimulation studies.

On the importance of precise electrode placement for targeted transcranial electric stimulation.

Transcranial electric stimulation (TES) is an increasingly popular method for non-invasive modulation of brain activity and a potential treatment for neuropsychiatric disorders. However, there are concerns about the reliability of its application because of variability in TES-induced intracranial electric fields across individuals. While realistic computational models offer can help to alleviate these concerns, their direct empirical validation is sparse, and their practical implications are not always clear. In this study, we combine direct intracranial measurements of electric fields generated by TES in surgical epilepsy patients with computational modeling. First, we directly validate the computational models and identify key parameters needed for accurate model predictions. Second, we derive practical guidelines for a reliable application of TES in terms of the precision of electrode placement needed to achieve a desired electric field distribution. Based on our results, we recommend electrode placement accuracy to be < 1 cm for a reliable application of TES across sessions.

Real-time impedance monitoring of oxygen reduction during surface modification of thin film cathodes.

Improvement of solid oxide fuel cells strongly relies on the development of cathode materials with high catalytic activity for the oxygen reduction reaction. Excellent activity was found for perovskite-type oxides such as La1-xSrxCoO3-δ (LSC), but performance degradation, probably caused by surface composition changes, hinders exploitation of the full potential of LSC. This study reveals that the potentially very high activity of the LSC surface can be traced back to few very active sites. Already tiny amounts of SrO, for example, 4% of a monolayer, deposited on an LSC surface, lead to severe deactivation. Co, on the other hand, causes (re-)activation, suggesting that active sites are strongly related to Co being present at the surface. These insights could be gained by a novel method to measure changes of the electrochemical performance of thin film electrodes in situ, while modifying their surface: impedance spectroscopy measurements during deposition of well-defined fractions of monolayers of Sr-, Co- and La-oxides by single laser pulses in a pulsed laser deposition chamber.

CO2 Reduction on the Pre-reduced Mixed Ionic-Electronic Conducting Perovskites La0.6 Sr-0.4 FeO3-δ and SrTi0.7 Fe0.3 O3-δ.

The activity of the pre-reduced perovskites La0.6 Sr0.4 FeO3-δ (LSF64) and SrTi0.7 Fe0.3 O3-δ (STF73) for the CO2 reduction to CO was investigated with special focus on the reactivity of oxide-dissolved hydrogen. This is of particular interest in hydrogen solid-oxide electrolysis cell (H-SOEC) technology, where proton-conducting ceramics are used and the reaction 2e- +2H+ +CO2 →CO+H2 O is of central importance. To clarify if hydrogen dissolved in LSF64 and STF73 partakes in the CO2 reduction, temperature-programmed reduction (TPR) in H2 , followed by temperature-programmed reoxidation (TPO) in CO2 and, moreover, temperature-programmed desorption (TPD) of ad- and absorbed species were utilized. The experiments reveal that 50 mol % of the CO2 is converted by hydrogen dissolved in STF73 and reacts quantitatively. On the other hand, LSF64 converts less than 20 mol % of CO2 via dissolved hydrogen and a residual of bulk OH is still detectable after CO2 -TPO.

Assessing Variations in Areal Organization for the Intrinsic Brain: From Fingerprints to Reliability.

Resting state fMRI (R-fMRI) is a powerful in-vivo tool for examining the functional architecture of the human brain. Recent studies have demonstrated the ability to characterize transitions between functionally distinct cortical areas through the mapping of gradients in intrinsic functional connectivity (iFC) profiles. To date, this novel approach has primarily been applied to iFC profiles averaged across groups of individuals, or in one case, a single individual scanned multiple times. Here, we used a publically available R-fMRI dataset, in which 30 healthy participants were scanned 10 times (10 min per session), to investigate differences in full-brain transition profiles (i.e., gradient maps, edge maps) across individuals, and their reliability. 10-min R-fMRI scans were sufficient to achieve high accuracies in efforts to "fingerprint" individuals based upon full-brain transition profiles. Regarding test-retest reliability, the image-wise intraclass correlation coefficient (ICC) was moderate, and vertex-level ICC varied depending on region; larger durations of data yielded higher reliability scores universally. Initial application of gradient-based methodologies to a recently published dataset obtained from twins suggested inter-individual variation in areal profiles might have genetic and familial origins. Overall, these results illustrate the utility of gradient-based iFC approaches for studying inter-individual variation in brain function.

Transcranial alternating current stimulation modulates spontaneous low frequency fluctuations as measured with fMRI.

Transcranial alternating current stimulation (tACS) is a promising tool for modulating brain oscillations. Combining tACS with functional magnetic resonance imaging (fMRI), we recently showed that tACS applied over the occipital cortex did not exert its strongest effect on regions below the electrodes, but mainly on more distant fronto-parietal regions. Theoretically, this effect could be explained by tACS-induced modulation of functional connectivity between directly stimulated areas and more distant but anatomically and functionally connected regions. In the present study, we aimed to characterize the effect of tACS on low frequency fMRI signal fluctuations. We employed simultaneous fMRI-tACS in 20 subjects during resting state (eyes open with central fixation for ~8min). Subjects received tACS at different frequencies (10, 16, 40Hz) and with different electrode montages (Cz-Oz, P5-P6) previously used in behavioral studies. Electric field simulations showed that tACS over Cz-Oz directly stimulates occipital cortex, while tACS over P5-P6 primarily targets parietal cortices. Group-level simulation-based functional connectivity maps for Cz-Oz and P5-P6 resembled the visual and fronto-parietal control resting-state networks, respectively. The effects of tACS were frequency and partly electrode montage dependent. In regions where frequency-dependent effects of tACS were observed, 10 and 40Hz tACS generally induced opposite effects. Most tACS effects on functional connectivity were observed between, as opposed to within, resting-state networks. The left fronto-parietal control network showed the most extensive frequency-dependent modulation in functional connectivity, mainly with occipito-parietal regions, where 10Hz tACS increased and 40Hz tACS decreased correlation values. Taken together, our results show that tACS modulates local spontaneous low frequency fluctuations and their correlations with more distant regions, which should be taken into account when interpreting tACS effects on brain function.

An integrated framework for targeting functional networks via transcranial magnetic stimulation.

Transcranial magnetic stimulation (TMS) is a powerful investigational tool for in vivo manipulation of regional or network activity, with a growing number of potential clinical applications. Unfortunately, the vast majority of targeting strategies remain limited by their reliance on non-realistic brain models and assumptions that anatomo-functional relationships are 1:1. Here, we present an integrated framework that combines anatomically realistic finite element models of the human head with resting functional MRI to predict functional networks targeted via TMS at a given coil location and orientation. Using data from the Human Connectome Project, we provide an example implementation focused on dorsolateral prefrontal cortex (DLPFC). Three distinct DLPFC stimulation zones were identified, differing with respect to the network to be affected (default, frontoparietal) and sensitivity to coil orientation. Network profiles generated for DLPFC targets previously published for treating depression revealed substantial variability across studies, highlighting a potentially critical technical issue.

Determinants of the electric field during transcranial direct current stimulation.

Transcranial direct current stimulation (tDCS) causes a complex spatial distribution of the electric current flow in the head which hampers the accurate localization of the stimulated brain areas. In this study we show how various anatomical features systematically shape the electric field distribution in the brain during tDCS. We constructed anatomically realistic finite element (FEM) models of two individual heads including conductivity anisotropy and different skull layers. We simulated a widely employed electrode montage to induce motor cortex plasticity and moved the stimulating electrode over the motor cortex in small steps to examine the resulting changes of the electric field distribution in the underlying cortex. We examined the effect of skull thickness and composition on the passing currents showing that thinner skull regions lead to higher electric field strengths. This effect is counteracted by a larger proportion of higher conducting spongy bone in thicker regions leading to a more homogenous current over the skull. Using a multiple regression model we could identify key factors that determine the field distribution to a significant extent, namely the thicknesses of the cerebrospinal fluid and the skull, the gyral depth and the distance to the anode and cathode. These factors account for up to 50% of the spatial variation of the electric field strength. Further, we demonstrate that individual anatomical factors can lead to stimulation "hotspots" which are partly resistant to electrode positioning. Our results give valuable novel insights in the biophysical foundation of tDCS and highlight the importance to account for individual anatomical factors when choosing an electrode montage.