Fast computational E-field dosimetry for transcranial magnetic stimulation using adaptive cross approximation and auxiliary dipole method (ACA-ADM).

Transcranial Magnetic Stimulation (TMS) is a non-invasive brain stimulation technique that uses a coil to induce an electric field (E-field) in the brain and modulate its activity. Many applications of TMS call for the repeated execution of E-field solvers to determine the E-field induced in the brain for different coil placements. However, the usage of solvers for these applications remains impractical because each coil placement requires the solution of a large linear system of equations. We develop a fast E-field solver that enables the rapid evaluation of the E-field distribution for a brain region of interest (ROI) for a large number of coil placements, which is achieved in two stages. First, during the pre-processing stage, the mapping between coil placement and brain ROI E-field distribution is approximated from E-field results for a few coil placements. Specifically, we discretize the mapping into a matrix with each column having the ROI E-field samples for a fixed coil placement. This matrix is approximated from a few of its rows and columns using adaptive cross approximation (ACA). The accuracy, efficiency, and applicability of the new ACA approach are determined by comparing its E-field predictions with analytical and standard solvers in spherical and MRI-derived head models. During the second stage, the E-field distribution in the brain ROI from a specific coil placement is determined by the obtained rows and columns in milliseconds. For many applications, only the E-field distribution for a comparatively small ROI is required. For example, the solver can complete the pre-processing stage in approximately 4 hours and determine the ROI E-field in approximately 40 ms for a 100 mm diameter ROI with less than 2% error enabling its use for neuro-navigation and other applications. Highlight: We developed a fast solver for TMS computational E-field dosimetry, which can determine the ROI E-field in approximately 40 ms for a 100 mm diameter ROI with less than 2% error.

Structural Controllability Predicts Functional Patterns and Brain Stimulation Benefits Associated with Working Memory.

The brain is an inherently dynamic system, and much work has focused on the ability to modify neural activity through both local perturbations and changes in the function of global network ensembles. Network controllability is a recent concept in network neuroscience that purports to predict the influence of individual cortical sites on global network states and state changes, thereby creating a unifying account of local influences on global brain dynamics. While this notion is accepted in engineering science, it is subject to ongoing debates in neuroscience as empirical evidence linking network controllability to brain activity and human behavior remains scarce. Here, we present an integrated set of multimodal brain-behavior relationships derived from fMRI, diffusion tensor imaging, and online repetitive transcranial magnetic stimulation (rTMS) applied during an individually calibrated working memory task performed by individuals of both sexes. The modes describing the structural network system dynamics showed direct relationships to brain activity associated with task difficulty, with difficult-to-reach modes contributing to functional brain states in the hard task condition. Modal controllability (a measure quantifying the contribution of difficult-to-reach modes) at the stimulated site predicted both fMRI activations associated with increasing task difficulty and rTMS benefits on task performance. Furthermore, fMRI explained 64% of the variance between modal controllability and the working memory benefit associated with 5 Hz online rTMS. These results therefore provide evidence toward the functional validity of network control theory, and outline a clear technique for integrating structural network topology and functional activity to predict the influence of stimulation on subsequent behavior.SIGNIFICANCE STATEMENT The network controllability concept proposes that specific cortical nodes are able to steer the brain into certain physiological states. By applying external perturbation to these control nodes, it is theorized that brain stimulation is able to selectively target difficult-to-reach states, potentially aiding processing and improving performance on cognitive tasks. The current study used rTMS and fMRI during a working memory task to test this hypothesis. We demonstrate that network controllability correlates with fMRI modulation because of working memory load and with the behavioral improvements that result from a multivisit intervention using 5 Hz rTMS. This study demonstrates the validity of network controllability and offers a new targeting approach to improve efficacy.

Fast computational optimization of TMS coil placement for individualized electric field targeting.

BACKGROUND: During transcranial magnetic stimulation (TMS) a coil placed on the scalp is used to non-invasively modulate activity of targeted brain networks via a magnetically induced electric field (E-field). Ideally, the E-field induced during TMS is concentrated on a targeted cortical region of interest (ROI). Determination of the coil position and orientation that best achieve this objective presently requires a large computational effort. OBJECTIVE: To improve the accuracy of TMS we have developed a fast computational auxiliary dipole method (ADM) for determining the optimum coil position and orientation. The optimum coil placement maximizes the E-field along a predetermined direction or, alternatively, the overall E-field magnitude in the targeted ROI. Furthermore, ADM can assess E-field uncertainty resulting from precision limitations of TMS coil placement protocols. METHOD: ADM leverages the electromagnetic reciprocity principle to compute rapidly the TMS induced E-field in the ROI by using the E-field generated by a virtual constant current source residing in the ROI. The framework starts by solving for the conduction currents resulting from this ROI current source. Then, it rapidly determines the average E-field induced in the ROI for each coil position by using the conduction currents and a fast-multipole method. To further speed-up the computations, the coil is approximated using auxiliary dipoles enabling it to represent all coil orientations for a given coil position with less than 600 dipoles. RESULTS: Using ADM, the E-fields generated in an MRI-derived head model when the coil is placed at 5900 different scalp positions and 360 coil orientations per position (over 2.1 million unique configurations) can be determined in under 15 min on a standard laptop computer. This enables rapid extraction of the optimum coil position and orientation as well as the E-field variation resulting from coil positioning uncertainty. ADM is implemented in SimNIBS 3.2. CONCLUSION: ADM enables the rapid determination of coil placement that maximizes E-field delivery to a specific brain target. This method can find the optimum coil placement in under 15 min enabling its routine use for TMS. Furthermore, it enables the fast quantification of uncertainty in the induced E-field due to limited precision of TMS coil placement protocols, enabling minimization and statistical analysis of the E-field dose variability.

Computationally optimized ECoG stimulation with local safety constraints.

Direct stimulation of the cortical surface is used clinically for cortical mapping and modulation of local activity. Future applications of cortical modulation and brain-computer interfaces may also use cortical stimulation methods. One common method to deliver current is through electrocorticography (ECoG) stimulation in which a dense array of electrodes are placed subdurally or epidurally to stimulate the cortex. However, proximity to cortical tissue limits the amount of current that can be delivered safely. It may be desirable to deliver higher current to a specific local region of interest (ROI) while limiting current to other local areas more stringently than is guaranteed by global safety limits. Two commonly used global safety constraints bound the total injected current and individual electrode currents. However, these two sets of constraints may not be sufficient to prevent high current density locally (hot-spots). In this work, we propose an efficient approach that prevents current density hot-spots in the entire brain while optimizing ECoG stimulus patterns for targeted stimulation. Specifically, we maximize the current along a particular desired directional field in the ROI while respecting three safety constraints: one on the total injected current, one on individual electrode currents, and the third on the local current density magnitude in the brain. This third set of constraints creates a computational barrier due to the huge number of constraints needed to bound the current density at every point in the entire brain. We overcome this barrier by adopting an efficient two-step approach. In the first step, the proposed method identifies the safe brain region, which cannot contain any hot-spots solely based on the global bounds on total injected current and individual electrode currents. In the second step, the proposed algorithm iteratively adjusts the stimulus pattern to arrive at a solution that exhibits no hot-spots in the remaining brain. We report on simulations on a realistic finite element (FE) head model with five anatomical ROIs and two desired directional fields. We also report on the effect of ROI depth and desired directional field on the focality of the stimulation. Finally, we provide an analysis of optimization runtime as a function of different safety and modeling parameters. Our results suggest that optimized stimulus patterns tend to differ from those used in clinical practice.

Visualizing simulated electrical fields from electroencephalography and transcranial electric brain stimulation: a comparative evaluation.

Electrical activity of neuronal populations is a crucial aspect of brain activity. This activity is not measured directly but recorded as electrical potential changes using head surface electrodes (electroencephalogram - EEG). Head surface electrodes can also be deployed to inject electrical currents in order to modulate brain activity (transcranial electric stimulation techniques) for therapeutic and neuroscientific purposes. In electroencephalography and noninvasive electric brain stimulation, electrical fields mediate between electrical signal sources and regions of interest (ROI). These fields can be very complicated in structure, and are influenced in a complex way by the conductivity profile of the human head. Visualization techniques play a central role to grasp the nature of those fields because such techniques allow for an effective conveyance of complex data and enable quick qualitative and quantitative assessments. The examination of volume conduction effects of particular head model parameterizations (e.g., skull thickness and layering), of brain anomalies (e.g., holes in the skull, tumors), location and extent of active brain areas (e.g., high concentrations of current densities) and around current injecting electrodes can be investigated using visualization. Here, we evaluate a number of widely used visualization techniques, based on either the potential distribution or on the current-flow. In particular, we focus on the extractability of quantitative and qualitative information from the obtained images, their effective integration of anatomical context information, and their interaction. We present illustrative examples from clinically and neuroscientifically relevant cases and discuss the pros and cons of the various visualization techniques.

A semi-automated pipeline for finite element modeling of electric field induced in nonhuman primates by transcranial magnetic stimulation.

BACKGROUND: Transcranial magnetic stimulation (TMS) is used to treat a range of brain disorders by inducing an electric field (E-field) in the brain. However, the precise neural effects of TMS are not well understood. Nonhuman primates (NHPs) are used to model the impact of TMS on neural activity, but a systematic method of quantifying the induced E-field in the cortex of NHPs has not been developed. NEW METHOD: The pipeline uses statistical parametric mapping (SPM) to automatically segment a structural MRI image of a rhesus macaque into five tissue compartments. Manual corrections are necessary around implants. The segmented tissues are tessellated into 3D meshes used in finite element method (FEM) software to compute the TMS induced E-field in the brain. The gray matter can be further segmented into cortical laminae using a volume preserving method for defining layers. RESULTS: Models of three NHPs were generated with TMS coils placed over the precentral gyrus. Two coil configurations, active and sham, were simulated and compared. The results demonstrated a large difference in E-fields at the target. Additionally, the simulations were calculated using two different E-field solvers and were found to not significantly differ. COMPARISON WITH EXISTING METHODS: Current methods segment NHP tissues manually or use automated methods for only the brain tissue. Existing methods also do not stratify the gray matter into layers. CONCLUSION: The pipeline calculates the induced E-field in NHP models by TMS and can be used to plan implant surgeries and determine approximate E-field values around neuron recording sites.

How electroconvulsive therapy works in the treatment of depression: is it the seizure, the electricity, or both?

We have known for nearly a century that triggering seizures can treat serious mental illness, but what we do not know is why. Electroconvulsive Therapy (ECT) works faster and better than conventional pharmacological interventions; however, those benefits come with a burden of side effects, most notably memory loss. Disentangling the mechanisms by which ECT exerts rapid therapeutic benefit from the mechanisms driving adverse effects could enable the development of the next generation of seizure therapies that lack the downside of ECT. The latest research suggests that this goal may be attainable because modifications of ECT technique have already yielded improvements in cognitive outcomes without sacrificing efficacy. These modifications involve changes in how the electricity is administered (both where in the brain, and how much), which in turn impacts the characteristics of the resulting seizure. What we do not completely understand is whether it is the changes in the applied electricity, or in the resulting seizure, or both, that are responsible for improved safety. Answering this question may be key to developing the next generation of seizure therapies that lack these adverse side effects, and ushering in novel interventions that are better, faster, and safer than ECT.

Electric Field Modeling in Personalizing Transcranial Magnetic Stimulation Interventions.

The modeling of transcranial magnetic stimulation (TMS)-induced electric fields (E-fields) is a versatile technique for evaluating and refining brain targeting and dosing strategies, while also providing insights into dose-response relationships in the brain. This review outlines the methodologies employed to derive E-field estimations, covering TMS physics, modeling assumptions, and aspects of subject-specific head tissue and coil modeling. We also summarize various numerical methods for solving the E-field and their suitability for various applications. Modeling methodologies have been optimized to efficiently execute numerous TMS simulations across diverse scalp coil configurations, facilitating the identification of optimal setups or rapid cortical E-field visualization for specific brain targets. These brain targets are extrapolated from neurophysiological measurements and neuroimaging, enabling precise and individualized E-field dosing in experimental and clinical applications. This necessitates the quantification of E-field estimates using metrics that enable the comparison of brain target engagement, functional localization, and TMS intensity adjustments across subjects. The integration of E-field modeling with empirical data has the potential to uncover pivotal insights into the aspects of E-fields responsible for stimulating and modulating brain function and states, enhancing behavioral task performance, and impacting the clinical outcomes of personalized TMS interventions.

Spatio-temporal regularization in linear distributed source reconstruction from EEG/MEG: a critical evaluation.

The high temporal resolution of EEG/MEG data offers a way to improve source reconstruction estimates which provide insight into the spatio-temporal involvement of neuronal sources in the human brain. In this work, we investigated the performance of spatio-temporal regularization (STR) in a current density approach using a systematic comparison to simple ad hoc or post hoc filtering of the data or of the reconstructed current density, respectively. For the used STR approach we implemented a frequency-specific constraint to penalize solutions outside a narrow frequency band of interest. The widely used sLORETA algorithm was adapted for STR and generally used for source reconstruction. STR and filtering approaches were evaluated with respect to spatial localization error and spatial dispersion, as well as to correlation of original and reconstructed source time courses in single source and two source scenarios with fixed source locations and oscillating source waveforms. We used extensive computer simulations and tested all algorithms with different parameter settings (noise levels and regularization parameters) for EEG data. To verify our results, we also used data from MEG phantom measurements. For the investigated scenarios, we did not find any evidence that STR-based methods outperform purely spatial algorithms applied to temporally filtered data. Furthermore, the results show very clearly that the performance of STR depends very much on the choice of regularization parameters.

Real-Time Computation of Brain E-Field for Enhanced Transcranial Magnetic Stimulation Neuronavigation and Optimization.

Transcranial Magnetic Stimulation (TMS) coil placement and pulse waveform current are often chosen to achieve a specified E-field dose on targeted brain regions. TMS neuronavigation could be improved by including real-time accurate distributions of the E-field dose on the cortex. We introduce a method and develop software for computing brain E-field distributions in real-time enabling easy integration into neuronavigation and with the same accuracy as 1st -order finite element method (FEM) solvers. Initially, a spanning basis set (< 400) of E-fields generated by white noise magnetic currents on a surface separating the head and permissible coil placements are orthogonalized to generate the modes. Subsequently, Reciprocity and Huygens' principles are utilized to compute fields induced by the modes on a surface separating the head and coil by FEM, which are used in conjunction with online (real-time) computed primary fields on the separating surface to evaluate the mode expansion. We conducted a comparative analysis of E-fields computed by FEM and in real-time for eight subjects, utilizing two head model types (SimNIBS's 'headreco' and 'mri2mesh' pipeline), three coil types (circular, double-cone, and Figure-8), and 1000 coil placements (48,000 simulations). The real-time computation for any coil placement is within 4 milliseconds (ms), for 400 modes, and requires less than 4 GB of memory on a GPU. Our solver is capable of computing E-fields within 4 ms, making it a practical approach for integrating E-field information into the neuronavigation systems without imposing a significant overhead on frame generation (20 and 50 frames per second within 50 and 20 ms, respectively).

Comparison of coil placement approaches targeting dorsolateral prefrontal cortex in depressed adolescents receiving repetitive transcranial magnetic stimulation: an electric field modeling study.

Background: A promising treatment option for adolescents with treatment-resistant depression is high-frequency repetitive transcranial magnetic stimulation (rTMS) delivered to the left dorsolateral prefrontal cortex (L-DLPFC). Conventional coil placement strategies for rTMS in adults include the 5-cm rule, the Beam F3 method, and the magnetic resonance imaging (MRI) neuronavigation method. The purpose of this study was to compare the three targeting approaches to a computational E-field optimization coil placement method in depressed adolescents. Methods: Ten consenting and assenting depressed adolescents (4 females, age: 15.9 ± 1.1) participated in an open-label rTMS treatment study. Participants were offered MRI-guided rTMS 5 times per week over 6-8 weeks. To compute the induced E-field, a head model was generated based on MRI images, and a figure-8 TMS coil (Neuronetics) was placed over the L-DLPFC using the four targeting approaches. Results: Results show that there was a significant difference in the induced E-field at the L-DLPFC between the four targeting methods ( χ 2 = 24.7, p < 0.001). Post hoc pairwise comparisons show that there was a significant difference between any two of the targeting methods (Holm adjusted p < 0.05), with the 5-cm rule producing the weakest E-field (46.0 ± 17.4 V/m), followed by the F3 method (87.4 ± 35.4 V/m), followed by the MRI-guided (112.1 ± 14.6 V/m), and followed by the computationally optimized method (130.1 ± 18.1 V/m). The Bartlett test of homogeneity of variances show that there was a significant difference in sample variance between the groups ( K 2 = 8.0, p < 0.05), with F3 having the largest variance. In participants who completed the full course of treatment, the median E-field strength in the L-DLPFC was correlated with the change in depression severity ( r = - 0.77, p < 0.05). Conclusions: The E-field models revealed inadequacies of scalp-based targeting methods compared to MRI-guidance. Computational optimization may further enhance E-field dose delivery to the treatment target.

Optimizing TMS Coil Placement Approaches for Targeting the Dorsolateral Prefrontal Cortex in Depressed Adolescents: An Electric Field Modeling Study.

High-frequency repetitive transcranial magnetic stimulation (rTMS) to the left dorsolateral prefrontal cortex (L-DLPFC) shows promise as a treatment for treatment-resistant depression in adolescents. Conventional rTMS coil placement strategies include the 5 cm, the Beam F3, and the magnetic resonance imaging (MRI) neuronavigation methods. The purpose of this study was to use electric field (E-field) models to compare the three targeting approaches to a computational E-field optimization coil placement method in depressed adolescents. Ten depressed adolescents (4 females, age: 15.9±1.1) participated in an open-label rTMS treatment study and were offered MRI-guided rTMS five times per week over 6-8 weeks. Head models were generated based on individual MRI images, and E-fields were simulated for the four targeting approaches. Results showed a significant difference in the induced E-fields at the L-DLPFC between the four targeting methods (χ2=24.7, p<0.001). Post hoc pairwise comparisons showed that there was a significant difference between any two of the targeting methods (Holm adjusted p<0.05), with the 5 cm rule producing the weakest E-field (46.0±17.4V/m), followed by the F3 method (87.4±35.4V/m), followed by MRI-guided (112.1±14.6V/m), and followed by the computational approach (130.1±18.1V/m). Variance analysis showed that there was a significant difference in sample variance between the groups (K2=8.0, p<0.05), with F3 having the largest variance. Participants who completed the full course of treatment had median E-fields correlated with depression symptom improvement (r=-0.77, p<0.05). E-field models revealed limitations of scalp-based methods compared to MRI guidance, suggesting computational optimization could enhance dose delivery to the target.

Multimodal smoking cessation treatment combining transcranial magnetic stimulation, cognitive behavioral therapy, and nicotine replacement therapy in veterans with posttraumatic stress disorder: A feasibility randomized controlled trial protocol.

Tobacco-related deaths exceed those resulting from homicides, suicides, motor vehicle accidence, alcohol consumption, illicit substance use, and acquired immunodeficiency syndrome (AIDS), combined. Amongst U.S. veterans, this trend is particularly concerning given that those suffering from posttraumatic stress disorder (PTSD)-about 11% of those receiving care from the Department of Veterans Affairs (VA)-have triple the risk of developing tobacco use disorder (TUD). The most efficacious strategies being used at the VA for smoking cessation only result in a 23% abstinence rate, and veterans with PTSD only achieve a 4.5% abstinence rate. Therefore, there is a critical need to develop more effective treatments for smoking cessation. Recent studies have revealed the insula as integrally involved in the neurocircuitry of TUD, specifically showing that individuals with brain lesions involving this region had drastically improved quit rates. Some of these studies show a probability of quitting up to 5 times greater compared to non-insula lesioned regions). Altered activity of the insula may be involved in the disruption of the salience network's (SN) connectivity to the executive control network (ECN), which compromises that patient's ability to switch between interoceptive states focused on cravings to executive and cognitive control. Thus, we propose a feasibility phase II randomized controlled trial (RCT) to study a patterned form of repetitive transcranial magnetic stimulation (rTMS), intermittent theta burst stimulation (iTBS), at 90% of the subject's resting motor threshold (rMT) applied over a region in the right post-central gyrus most functionally connected to the right posterior insula. We hypothesize that by increasing functional connectivity between the SN with the ECN to enhance executive control and by decreasing connectivity with the default mode network (DMN) to reduce interoceptive focus on withdrawal symptoms, we will improve smoking cessation outcomes. Fifty eligible veterans with comorbid TUD and PTSD will be randomly assigned to two conditions: active-iTBS + cognitive behavioral therapy (CBT) + nicotine replacement therapy (NRT) (n=25) or sham-iTBS + CBT + NRT (n=25). The primary outcome, feasibility, will be determined by achieving a recruitment of 50 participants and retention rate of 80%. The success of iTBS will be evaluated through self-reported nicotine use, cravings, withdrawal symptoms, and abstinence following quit date (confirmed by bioverification) along with evaluation for target engagement through neuroimaging changes, specifically connectivity differences between the insula and other regions of interest.

Social appropriateness perception of dynamic interactions.

The current study explored the judgment of communicative appropriateness while processing a dialogue between two individuals. All stimuli were presented as audio-visual as well as audio-only vignettes and 24 young adults reported their social impression (appropriateness) of literal, blunt, sarcastic, and teasing statements. On average, teasing statements were rated as more appropriate when processing audio-visual statements compared to the audio-only version of a stimuli, while sarcastic statements were judged as less appropriate with additional visual information. These results indicate a rejection of the Tinge Hypothesis for audio-visual vignettes while confirming it for the reduced, audio-only counterparts. We also analyzed time-frequency EEG data of four frequency bands that have been related to language processing: alpha, beta, theta and low gamma. We found desynchronization in the alpha band literal versus nonliteral items, confirming the assumption that the alpha band reflects stimulus complexity. The analysis also revealed a power increase in the theta, beta and low gamma band, especially when comparing blunt and nonliteral statements in the audio-only condition. The time-frequency results corroborate the prominent role of the alpha and theta bands in language processing and offer new insights into the neural correlates of communicative appropriateness and social aspects of speech perception.

TAP: targeting and analysis pipeline for optimization and verification of coil placement in transcranial magnetic stimulation.

Objective.Transcranial magnetic stimulation (TMS) can modulate brain function via an electric field (E-field) induced in a brain region of interest (ROI). The ROI E-field can be computationally maximized and set to match a specific reference using individualized head models to find the optimal coil placement and stimulus intensity. However, the available software lacks many practical features for prospective planning of TMS interventions and retrospective evaluation of the experimental targeting accuracy.Approach.The TMS targeting and analysis pipeline (TAP) software uses an MRI/fMRI-derived brain target to optimize coil placement considering experimental parameters such as the subject's hair thickness and coil placement restrictions. The coil placement optimization is implemented in SimNIBS 3.2, for which an additional graphical user interface (TargetingNavigator) is provided to visualize/adjust procedural parameters. The coil optimization process also computes the E-field at the target, allowing the selection of the TMS device intensity setting to achieve specific E-field strengths. The optimized coil placement information is prepared for neuronavigation software, which supports targeting during the TMS procedure. The neuronavigation system can record the coil placement during the experiment, and these data can be processed in TAP to quantify the accuracy of the experimental TMS coil placement and induced E-field.Main results.TAP was demonstrated in a study consisting of three repetitive TMS sessions in five subjects. TMS was delivered by an experienced operator under neuronavigation with the computationally optimized coil placement. Analysis of the experimental accuracy from the recorded neuronavigation data indicated coil location and orientation deviations up to about 2 mm and 2°, respectively, resulting in an 8% median decrease in the target E-field magnitude compared to the optimal placement.Significance.TAP supports navigated TMS with a variety of features for rigorous and reproducible stimulation delivery, including planning and evaluation of coil placement and intensity selection for E-field-based dosing.

Modeling of the human skull in EEG source analysis.

We used computer simulations to investigate finite element models of the layered structure of the human skull in EEG source analysis. Local models, where each skull location was modeled differently, and global models, where the skull was assumed to be homogeneous, were compared to a reference model, in which spongy and compact bone were explicitly accounted for. In both cases, isotropic and anisotropic conductivity assumptions were taken into account. We considered sources in the entire brain and determined errors both in the forward calculation and the reconstructed dipole position. Our results show that accounting for the local variations over the skull surface is important, whereas assuming isotropic or anisotropic skull conductivity has little influence. Moreover, we showed that, if using an isotropic and homogeneous skull model, the ratio between skin/brain and skull conductivities should be considerably lower than the commonly used 80:1. For skull modeling, we recommend (1) Local models: if compact and spongy bone can be identified with sufficient accuracy (e.g., from MRI) and their conductivities can be assumed to be known (e.g., from measurements), one should model these explicitly by assigning each voxel to one of the two conductivities, (2) Global models: if the conditions of (1) are not met, one should model the skull as either homogeneous and isotropic, but with considerably higher skull conductivity than the usual 0.0042 S/m, or as homogeneous and anisotropic, but with higher radial skull conductivity than the usual 0.0042 S/m and a considerably lower radial:tangential conductivity anisotropy than the usual 1:10.

Effects of a five-day HD-tDCS application to the right IFG depend on current intensity: A study in children and adolescents with ADHD.

Impaired executive functions in ADHD are associated with hypoactivity of the right inferior frontal gyrus (IFG). This region was targeted via repetitive applications of anodal, high-definition transcranial direct current simulation (HD-tDCS) on five consecutive days in 33 ADHD patients (10-17years) and in a healthy control group (n=13, only sham). Patients received either sham (n=13) or verum tDCS with 0.5mA (n=9) or 0.25mA (n=11) depending on individual cutaneous sensitivity. During stimulation, participants performed a combined working memory and response inhibition paradigm (n-back/nogo). At baseline, post, and a 4-month follow up, electroencephalography was recorded during this task. Moreover, interference control (flanker task) and spatial working memory (spanboard task) were assessed to explore possible transfer effects. Omission errors and reaction time variability in all tasks served as measures of attention. In the 0.25mA group increased nogo commission errors indicated a detrimental tDCS effect on response inhibition. After the 5-day stimulation, attentional improvements in the 0.5mA group were indicated by reduced omission errors and reaction time variability. Variability improvements were still evident at follow up. In all groups, nogo P3 amplitudes were reduced post-stimulation, but in the 0.5mA group this reduction was smaller than in the 0.25mA group. Results of the current study suggest distinct effects of tDCS with different current intensities demonstrating the importance of a deeper understanding on the impact of stimulation parameters and repeated tDCS applications to develop effective tDCS-based therapy approaches in ADHD.

Transcranial alternating current stimulation entrains alpha oscillations by preferential phase synchronization of fast-spiking cortical neurons to stimulation waveform.

Computational modeling and human studies suggest that transcranial alternating current stimulation (tACS) modulates alpha oscillations by entrainment. Yet, a direct examination of how tACS interacts with neuronal spiking activity that gives rise to the alpha oscillation in the thalamo-cortical system has been lacking. Here, we demonstrate how tACS entrains endogenous alpha oscillations in head-fixed awake ferrets. We first show that endogenous alpha oscillations in the posterior parietal cortex drive the primary visual cortex and the higher-order visual thalamus. Spike-field coherence is largest for the alpha frequency band, and presumed fast-spiking inhibitory interneurons exhibit strongest coupling to this oscillation. We then apply alpha-tACS that results in a field strength comparable to what is commonly used in humans (<0.5 mV/mm). Both in these ferret experiments and in a computational model of the thalamo-cortical system, tACS entrains alpha oscillations by following the theoretically predicted Arnold tongue. Intriguingly, the fast-spiking inhibitory interneurons exhibit a stronger entrainment response to tACS in both the ferret experiments and the computational model, likely due to their stronger endogenous coupling to the alpha oscillation. Our findings demonstrate the in vivo mechanism of action for the modulation of the alpha oscillation by tACS.

Network-based rTMS to modulate working memory: The difficult choice of effective parameters for online interventions.

BACKGROUND: Online repetitive transcranialmagnetic stimulation (rTMS) has been shown to modulate working memory (WM) performance in a site-specific manner, with behavioral improvements due to stimulation of the dorsolateral prefrontal cortex (DLPFC), and impairment from stimulation to the lateral parietal cortex (LPC). Neurobehavioral studies have demonstrated that subprocesses of WM allowing for the maintenance and manipulation of information in the mind involve unique cortical networks. Despite promising evidence of modulatory effects of rTMS on WM, no studies have yet demonstrated distinct modulatory control of these two subprocesses. The current study therefore sought to explore this possibility through site-specific stimulation during an online task invoking both skills. METHODS: Twenty-nine subjects completed a 4-day protocol, in which active or sham 5Hz rTMS was applied over the DLPFC and LPC in separate blocks of trials while participants performed tasks that required either maintenance alone, or both maintenance and manipulation (alphabetization) of information. Stimulation targets were defined individually based on fMRI activation and structural network properties. Stimulation amplitude was adjusted using electric field modeling to equate induced current in the target region across participants. RESULTS: Despite the use of advanced techniques, no significant differences or interactions between active and sham stimulation were found. Exploratory analyses testing stimulation amplitude, fMRI activation, and modal controllability showed nonsignificant but interesting trends with rTMS effects. CONCLUSION: While this study did not reveal any significant behavioral changes in WM, the results may point to parameters that contribute to positive effects, such as stimulation amplitude and functional activation.

Conditions for numerically accurate TMS electric field simulation.

BACKGROUND: Computational simulations of the E-field induced by transcranial magnetic stimulation (TMS) are increasingly used to understand its mechanisms and to inform its administration. However, characterization of the accuracy of the simulation methods and the factors that affect it is lacking. OBJECTIVE: To ensure the accuracy of TMS E-field simulations, we systematically quantify their numerical error and provide guidelines for their setup. METHOD: We benchmark the accuracy of computational approaches that are commonly used for TMS E-field simulations, including the finite element method (FEM) with and without superconvergent patch recovery (SPR), boundary element method (BEM), finite difference method (FDM), and coil modeling methods. RESULTS: To achieve cortical E-field error levels below 2%, the commonly used FDM and 1st order FEM require meshes with an average edge length below 0.4 mm, 1st order SPR-FEM requires edge lengths below 0.8 mm, and BEM and 2nd (or higher) order FEM require edge lengths below 2.9 mm. Coil models employing magnetic and current dipoles require at least 200 and 3000 dipoles, respectively. For thick solid-conductor coils and frequencies above 3 kHz, winding eddy currents may have to be modeled. CONCLUSION: BEM, FDM, and FEM all converge to the same solution. Compared to the common FDM and 1st order FEM approaches, BEM and 2nd (or higher) order FEM require significantly lower mesh densities to achieve the same error level. In some cases, coil winding eddy-currents must be modeled. Both electric current dipole and magnetic dipole models of the coil current can be accurate with sufficiently fine discretization.