Enhanced tES and tDCS computational models by meninges emulation

Perspectives on Optimized Transcranial Electrical Stimulation Based on Spatial Electric Field Modeling in Humans

Background: Transcranial electrical stimulation (tES) generates an electric field (or current density) in the brain through surface electrodes attached to the scalp. Clinical significance has been demonstrated, although with moderate and heterogeneous results partly due to a lack of control of the delivered electric currents. In the last decade, computational electric field analysis has allowed the estimation and optimization of the electric field using accurate anatomical head models. This review examines recent tES computational studies, providing a comprehensive background on the technical aspects of adopting computational electric field analysis as a standardized procedure in medical applications. Methods: Specific search strategies were designed to retrieve papers from the Web of Science database. The papers were initially screened based on the soundness of the title and abstract and then on their full contents, resulting in a total of 57 studies. Results: Recent trends were identified in individual- and population-level analysis of the electric field, including head models from non-neurotypical individuals. Advanced optimization techniques that allow a high degree of control with the required focality and direction of the electric field were also summarized. There is also growing evidence of a correlation between the computationally estimated electric field and the observed responses in real experiments. Conclusions: Computational pipelines and optimization algorithms have reached a degree of maturity that provides a rationale to improve tES experimental design and a posteriori analysis of the responses for supporting clinical studies.

Quasistatic approximation in neuromodulation

We define and explain the quasistatic approximation (QSA) as applied to field modeling for electrical and magnetic stimulation. Neuromodulation analysis pipelines include discrete stages, and QSA is applied specifically when calculating the electric and magnetic fields generated in tissues by a given stimulation dose. QSA simplifies the modeling equations to support tractable analysis, enhanced understanding, and computational efficiency. The application of QSA in neuro-modulation is based on four underlying assumptions: (A1) no wave propagation or self-induction in tissue, (A2) linear tissue properties, (A3) purely resistive tissue, and (A4) non-dispersive tissue. As a consequence of these assumptions, each tissue is assigned a fixed conductivity, and the simplified equations (e.g., Laplace's equation) are solved for the spatial distribution of the field, which is separated from the field's temporal waveform. Recognizing that electrical tissue properties may be more complex, we explain how QSA can be embedded in parallel or iterative pipelines to model frequency dependence or nonlinearity of conductivity. We survey the history and validity of QSA across specific applications, such as microstimulation, deep brain stimulation, spinal cord stimulation, transcranial electrical stimulation, and transcranial magnetic stimulation. The precise definition and explanation of QSA in neuromodulation are essential for rigor when using QSA models or testing their limits.

An adaptive h-refinement method for the boundary element fast multipole method for quasi-static electromagnetic modeling

Objective.In our recent work pertinent to modeling of brain stimulation and neurophysiological recordings, substantial modeling errors in the computed electric field and potential have sometimes been observed for standard multi-compartment head models. The goal of this study is to quantify those errors and, further, eliminate them through an adaptive mesh refinement (AMR) algorithm. The study concentrates on transcranial magnetic stimulation (TMS), transcranial electrical stimulation (TES), and electroencephalography (EEG) forward problems.Approach.We propose, describe, and systematically investigate an AMR method using the boundary element method with fast multipole acceleration (BEM-FMM) as the base numerical solver. The goal is to efficiently allocate additional unknowns to critical areas of the model, where they will best improve solution accuracy. The implemented AMR method's accuracy improvement is measured on head models constructed from 16 Human Connectome Project subjects under problem classes of TES, TMS, and EEG. Errors are computed between three solutions: an initial non-adaptive solution, a solution found after applying AMR with a conservative refinement rate, and a 'silver-standard' solution found by subsequent 4:1 global refinement of the adaptively-refined model.Main results.Excellent agreement is shown between the adaptively-refined and silver-standard solutions for standard head models. AMR is found to be vital for accurate modeling of TES and EEG forward problems for standard models: an increase of less than 25% (on average) in number of mesh elements for these problems, efficiently allocated by AMR, exposes electric field/potential errors exceeding 60% (on average) in the solution for the unrefined models.Significance.This error has especially important implications for TES dosing prediction-where the stimulation strength plays a central role-and for EEG lead fields. Though the specific form of the AMR method described here is implemented for the BEM-FMM, we expect that AMR is applicable and even required for accurate electromagnetic simulations by other numerical modeling packages as well.

Transcranial electric stimulation modulates firing rate at clinically relevant intensities

Notwithstanding advances with low-intensity transcranial electrical stimulation (TES), there remain questions about the efficacy of clinically realistic electric fields on neuronal function. We used Neuropixels 2.0 probe with 384 channels in an in-vivo rat model of TES to detect effects of weak fields on neuronal firing rate. High-density field mapping and computational models verified field intensity (1 V/m in hippocampus per 50 μA of applied skull currents). We demonstrate that electric fields below 0.5 V/m acutely modulate firing rate in 5% of neurons recorded in the hippocampus. At these intensities, average firing rate effects increased monotonically with electric field intensity at a rate of 7 % per V/m. For the majority of excitatory neurons, firing increased for cathodal stimulation and diminished for anodal stimulation. While more diverse, the response of inhibitory neurons followed a similar pattern on average, likely as a result of excitatory drive. Our results indicate that responses to TES at clinically relevant intensities are driven by a fraction of high-responder excitatory neurons, with polarity-specific effects. We conclude that transcranial electric stimulation is an effective neuromodulator at clinically realistic intensities.

Outcome measures for electric field modeling in tES and TMS: A systematic review and large-scale modeling study

BACKGROUND

Electric field (E-field) modeling is a potent tool to estimate the amount of transcranial magnetic and electrical stimulation (TMS and tES, respectively) that reaches the cortex and to address the variable behavioral effects observed in the field. However, outcome measures used to quantify E-fields vary considerably and a thorough comparison is missing.

OBJECTIVES

This two-part study aimed to examine the different outcome measures used to report on tES and TMS induced E-fields, including volume- and surface-level gray matter, region of interest (ROI), whole brain, geometrical, structural, and percentile-based approaches. The study aimed to guide future research in informed selection of appropriate outcome measures.

METHODS

Three electronic databases were searched for tES and/or TMS studies quantifying E-fields. The identified outcome measures were compared across volume- and surface-level E-field data in ten tES and TMS modalities targeting two common targets in 100 healthy individuals.

RESULTS

In the systematic review, we extracted 308 outcome measures from 202 studies that adopted either a gray matter volume-level (n = 197) or surface-level (n = 111) approach. Volume-level results focused on E-field magnitude, while surface-level data encompassed E-field magnitude (n = 64) and normal/tangential E-field components (n = 47). E-fields were extracted in ROIs, such as brain structures and shapes (spheres, hexahedra and cylinders), or the whole brain. Percentiles or mean values were mostly used to quantify E-fields. Our modeling study, which involved 1,000 E-field models and > 1,000,000 extracted E-field values, revealed that different outcome measures yielded distinct E-field values, analyzed different brain regions, and did not always exhibit strong correlations in the same within-subject E-field model.

CONCLUSIONS

Outcome measure selection significantly impacts the locations and intensities of extracted E-field data in both tES and TMS E-field models. The suitability of different outcome measures depends on the target region, TMS/tES modality, individual anatomy, the analyzed E-field component and the research question. To enhance the quality, rigor, and reproducibility in the E-field modeling domain, we suggest standard reporting practices across studies and provide four recommendations.

Dose-response in modulating brain function with transcranial direct current stimulation: From local to network levels

Understanding the dose-response relationship is crucial in studying the effects of brain stimulation techniques, such as transcranial direct current stimulation (tDCS). The dose-response relationship refers to the relationship between the received stimulation dose and the resulting response, which can be described as a function of the dose at various levels, including single/multiple neurons, clusters, regions, or networks. Here, we are focused on the received stimulation dose obtained from computational head models and brain responses which are quantified by functional magnetic resonance imaging (fMRI) data. In this randomized, triple-blind, sham-controlled clinical trial, we recruited sixty participants with methamphetamine use disorders (MUDs) as a sample clinical population who were randomly assigned to receive either sham or active tDCS. Structural and functional MRI data, including high-resolution T1 and T2-weighted MRI, resting-state functional MRI, and a methamphetamine cue-reactivity task fMRI, were acquired before and after tDCS. Individual head models were generated using the T1 and T2-weighted MRI data to simulate electric fields. In a linear approach, we investigated the associations between electric fields (received dose) and changes in brain function (response) at four different levels: voxel level, regional level (using atlas-based parcellation), cluster level (identifying active clusters), and network level (task-based functional connectivity). At the voxel level, regional level, and cluster level, no FDR-corrected significant correlation was observed between changes in functional activity and electric fields. However, at the network level, a significant positive correlation was found between frontoparietal connectivity and the electric field at the frontopolar stimulation site (r = 0.42, p corrected = 0.02; medium effect size). Our proposed pipeline offers a methodological framework for analyzing tDCS effects by exploring dose-response relationships at different levels, enabling a direct link between electric field variability and the neural response to tDCS. The results indicate that network-based analysis provides valuable insights into the dependency of tDCS neuromodulatory effects on the individual's regional current dose. Integration of dose-response relationships can inform dose optimization, customization, or the extraction of predictive/treatment-response biomarkers in future brain stimulation studies.

The Rehabilitation Potential of Neurostimulation for Mild Traumatic Brain Injury in Animal and Human Studies

Neurostimulation carries high therapeutic potential, accompanied by an excellent safety profile. In this review, we argue that an arena in which these tools could provide breakthrough benefits is traumatic brain injury (TBI). TBI is a major health problem worldwide, with the majority of cases identified as mild TBI (mTBI). MTBI is of concern because it is a modifiable risk factor for dementia. A major challenge in studying mTBI is its inherent heterogeneity across a large feature space (e.g., etiology, age of injury, sex, treatment, initial health status, etc.). Parallel lines of research in human and rodent mTBI can be collated to take advantage of the full suite of neuroscience tools, from neuroimaging (electroencephalography: EEG; functional magnetic resonance imaging: fMRI; diffusion tensor imaging: DTI) to biochemical assays. Despite these attractive components and the need for effective treatments, there are at least two major challenges to implementation. First, there is insufficient understanding of how neurostimulation alters neural mechanisms. Second, there is insufficient understanding of how mTBI alters neural function. The goal of this review is to assemble interrelated but disparate areas of research to identify important gaps in knowledge impeding the implementation of neurostimulation.

An Adaptive H-Refinement Method for the Boundary Element Fast Multipole Method for Quasi-static Electromagnetic Modeling

Objective

In our recent work pertinent to modeling of brain stimulation and neurophysiological recordings, substantial modeling errors in the computed electric field and potential have sometimes been observed for standard multi-compartment head models. The goal of this study is to quantify those errors and, further, eliminate them through an adaptive mesh refinement (AMR) algorithm. The study concentrates on transcranial magnetic stimulation (TMS), transcranial electrical stimulation (TES), and electroencephalography (EEG) forward problems.

Approach

We propose, describe, and systematically investigate an AMR method using the Boundary Element Method with Fast Multipole Acceleration (BEM-FMM) as the base numerical solver. The goal is to efficiently allocate additional unknowns to critical areas of the model, where they will best improve solution accuracy.The implemented AMR method's accuracy improvement is measured on head models constructed from 16 Human Connectome Project subjects under problem classes of TES, TMS, and EEG. Errors are computed between three solutions: an initial non-adaptive solution, a solution found after applying AMR with a conservative refinement rate, and a "silver-standard" solution found by subsequent 4:1 global refinement of the adaptively-refined model.

Main Results

Excellent agreement is shown between the adaptively-refined and silver-standard solutions for standard head models. AMR is found to be vital for accurate modeling of TES and EEG forward problems for standard models: an increase of less than 25% (on average) in number of mesh elements for these problems, efficiently allocated by AMR, exposes electric field/potential errors exceeding 60% (on average) in the solution for the unrefined models.

Significance

This error has especially important implications for TES dosing prediction - where the stimulation strength plays a central role - and for EEG lead fields. Though the specific form of the AMR method described here is implemented for the BEM-FMM, we expect that AMR is applicable and even required for accurate electromagnetic simulations by other numerical modeling packages as well.

Multi-scale multi-physics model of brain interstitial water flux by transcranial Direct Current Stimulation

Objective. Transcranial direct current stimulation (tDCS) generates sustained electric fields in the brain, that may be amplified when crossing capillary walls (across blood-brain barrier, BBB). Electric fields across the BBB may generate fluid flow by electroosmosis. We consider that tDCS may thus enhance interstitial fluid flow.Approach. We developed a modeling pipeline novel in both (1) spanning the mm (head),μm (capillary network), and then nm (down to BBB tight junction (TJ)) scales; and (2) coupling electric current flow to fluid current flow across these scales. Electroosmotic coupling was parametrized based on prior measures of fluid flow across isolated BBB layers. Electric field amplification across the BBB in a realistic capillary network was converted to volumetric fluid exchange.Main results. The ultrastructure of the BBB results in peak electric fields (per mA of applied current) of 32-63Vm-1across capillary wall and >1150Vm-1in TJs (contrasted with 0.3Vm-1in parenchyma). Based on an electroosmotic coupling of 1.0 × 10-9- 5.6 × 10-10m3s-1m2perVm-1, peak water fluxes across the BBB are 2.44 × 10-10- 6.94 × 10-10m3s-1m2, with a peak 1.5 × 10-4- 5.6 × 10-4m3min-1m3interstitial water exchange (per mA).Significance. Using this pipeline, the fluid exchange rate per each brain voxel can be predicted for any tDCS dose (electrode montage, current) or anatomy. Under experimentally constrained tissue properties, we predicted tDCS produces a fluid exchange rate comparable to endogenous flow, so doubling fluid exchange with further local flow rate hot spots ('jets'). The validation and implication of such tDCS brain 'flushing' is important to establish.

Quasi-static pipeline in electroconvulsive therapy computational modeling

BACKGROUND

Computational models of current flow during Electroconvulsive Therapy (ECT) rely on the quasi-static assumption, yet tissue impedance during ECT may be frequency specific and change adaptively to local electric field intensity.

OBJECTIVES

We systematically consider the application of the quasi-static pipeline to ECT under conditions where 1) static impedance is measured before ECT and 2) during ECT when dynamic impedance is measured. We propose an update to ECT modeling accounting for frequency-dependent impedance.

METHODS

The frequency content on an ECT device output is analyzed. The ECT electrode-body impedance under low-current conditions is measured with an impedance analyzer. A framework for ECT modeling under quasi-static conditions based on a single device-specific frequency (e.g., 1 kHz) is proposed.

RESULTS

Impedance using ECT electrodes under low-current is frequency dependent and subject specific, and can be approximated at >100 Hz with a subject-specific lumped parameter circuit model but at <100 Hz increased non-linearly. The ECT device uses a 2 μA 800 Hz test signal and reports a static impedance that approximate 1 kHz impedance. Combined with prior evidence suggesting that conductivity does not vary significantly across ECT output frequencies at high-currents (800-900 mA), we update the adaptive pipeline for ECT modeling centered at 1 kHz frequency. Based on individual MRI and adaptive skin properties, models match static impedance (at 2 μA) and dynamic impedance (at 900 mA) of four ECT subjects.

CONCLUSIONS

By considering ECT modeling at a single representative frequency, ECT adaptive and non-adaptive modeling can be rationalized under a quasi-static pipeline.

UncertainSCI: Uncertainty quantification for computational models in biomedicine and bioengineering

BACKGROUND

Computational biomedical simulations frequently contain parameters that model physical features, material coefficients, and physiological effects, whose values are typically assumed known a priori. Understanding the effect of variability in those assumed values is currently a topic of great interest. A general-purpose software tool that quantifies how variation in these parameters affects model outputs is not broadly available in biomedicine. For this reason, we developed the 'UncertainSCI' uncertainty quantification software suite to facilitate analysis of uncertainty due to parametric variability.

METHODS

We developed and distributed a new open-source Python-based software tool, UncertainSCI, which employs advanced parameter sampling techniques to build polynomial chaos (PC) emulators that can be used to predict model outputs for general parameter values. Uncertainty of model outputs is studied by modeling parameters as random variables, and model output statistics and sensitivities are then easily computed from the emulator. Our approaches utilize modern, near-optimal techniques for sampling and PC construction based on weighted Fekete points constructed by subsampling from a suitably randomized candidate set.

RESULTS

Concentrating on two test cases-modeling bioelectric potentials in the heart and electric stimulation in the brain-we illustrate the use of UncertainSCI to estimate variability, statistics, and sensitivities associated with multiple parameters in these models.

CONCLUSION

UncertainSCI is a powerful yet lightweight tool enabling sophisticated probing of parametric variability and uncertainty in biomedical simulations. Its non-intrusive pipeline allows users to leverage existing software libraries and suites to accurately ascertain parametric uncertainty in a variety of applications.

Neurocapillary-Modulation

OBJECTIVES

We consider two consequences of brain capillary ultrastructure in neuromodulation. First, blood-brain barrier (BBB) polarization as a consequence of current crossing between interstitial space and the blood. Second, interstitial current flow distortion around capillaries impacting neuronal stimulation.

MATERIALS AND METHODS

We developed computational models of BBB ultrastructure morphologies to first assess electric field amplification at the BBB (principle 1) and neuron polarization amplification by the presence of capillaries (principle 2). We adapt neuron cable theory to develop an analytical solution for maximum BBB polarization sensitivity.

RESULTS

Electrical current crosses between the brain parenchyma (interstitial space) and capillaries, producing BBB electric fields (EBBB) that are >400x of the average parenchyma electric field (ĒBRAIN), which in turn modulates transport across the BBB. Specifically, for a BBB space constant (λBBB) and wall thickness (dth-BBB), the analytical solution for maximal BBB electric field (EABBB) is given as: (ĒBRAIN × λBBB)/dth-BBB. Electrical current in the brain parenchyma is distorted around brain capillaries, amplifying neuronal polarization. Specifically, capillary ultrastructure produces ∼50% modulation of the ĒBRAIN over the ∼40 μm inter-capillary distance. The divergence of EBRAIN (Activating function) is thus ∼100 kV/m2 per unit ĒBRAIN.

CONCLUSIONS

BBB stimulation by principle 1 suggests novel therapeutic strategies such as boosting metabolic capacity or interstitial fluid clearance. Whereas the spatial profile of EBRAIN is traditionally assumed to depend only on macroscopic anatomy, principle 2 suggests a central role for local capillary ultrastructure-which impact forms of neuromodulation including deep brain stimulation (DBS), spinal cord stimulation (SCS), transcranial magnetic stimulation (TMS), electroconvulsive therapy (ECT), and transcranial electrical stimulation (tES)/transcranial direct current stimulation (tDCS).

Optimized APPS-tDCS electrode position, size, and distance doubles the on-target stimulation magnitude in 3000 electric field models

Transcranial direct current stimulation (tDCS) is a widely used noninvasive brain stimulation technique with mixed results to date. A potential solution is to apply more efficient stimulation to ensure that each participant receives sufficient cortical activation. In this four-part study, we used electric field (E-field) modeling to systematically investigate the cortical effects of conventional and novel tDCS electrode montages, with the goal of creating a new easily adoptable form of tDCS that induces higher and more focal E-fields. We computed 3000 anatomically accurate, MRI-based E-field models using 2 mA tDCS to target the left primary motor cortex in 200 Human Connectome Project (HCP) participants and tested the effects of: 1. Novel Electrode Position, 2. Electrode Size, and 3. Inter-Electrode Distance on E-field magnitude and focality. In particular, we examined the effects of placing electrodes surrounding the corticomotor target in the anterior and posterior direction (anterior posterior pad surround tDCS; APPS-tDCS). We found that electrode position, electrode size, and inter-electrode distance all significantly impact the cortical E-field magnitude and focality of stimulation (all p < 0.0001). At the same 2 mA scalp stimulation intensity, APPS-tDCS with smaller than conventional 1 × 1 cm electrodes surrounding the neural target deliver more than double the on-target cortical E-field (APPS-tDCS: average of 0.55 V/m from 2 mA; M1-SO and bilateral M1: both 0.27 V/m from 2 mA) while stimulating only a fraction of the off-target brain regions; 2 mA optimized APPS-tDCS produces 4.08 mA-like cortical E-fields. In sum, this new optimized APPS-tDCS method produces much stronger cortical stimulation intensities at the same 2 mA scalp intensity. APPS-tDCS also more focally stimulates the cortex at the intended target, using simple EEG coordinate locations and without MRI scans. This APPS-tDCS method is adoptable to any existing, commercially available tDCS device and can be used to ensure sufficient cortical activation in each person. Future directions include testing whether APPS-tDCS produces larger and more consistent therapeutic tDCS effects.

Personalized tDCS for Focal Epilepsy—A Narrative Review: A Data-Driven Workflow Based on Imaging and EEG Data

Conventional transcranial electric stimulation(tES) using standard anatomical positions for the electrodes and standard stimulation currents is frequently not sufficiently selective in targeting and reaching specific brain locations, leading to suboptimal application of electric fields. Recent advancements in in vivo electric field characterization may enable clinical researchers to derive better relationships between the electric field strength and the clinical results. Subject-specific electric field simulations could lead to improved electrode placement and more efficient treatments. Through this narrative review, we present a processing workflow to personalize tES for focal epilepsy, for which there is a clear cortical target to stimulate. The workflow utilizes clinical imaging and electroencephalography data and enables us to relate the simulated fields to clinical outcomes. We review and analyze the relevant literature for the processing steps in the workflow, which are the following: tissue segmentation, source localization, and stimulation optimization. In addition, we identify shortcomings and ongoing trends with regard to, for example, segmentation quality and tissue conductivity measurements. The presented processing steps result in personalized tES based on metrics like focality and field strength, which allow for correlation with clinical outcomes.

The effect of meninges on the electric fields in TES and TMS. Numerical modeling with adaptive mesh refinement

BACKGROUND

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.

OBJECTIVE

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.

METHOD

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.

RESULTS

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.

CONCLUSION

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.

Individually optimized multi-channel tDCS for targeting somatosensory cortex

OBJECTIVE

Transcranial direct current stimulation (tDCS) is a non-invasive neuro-modulation technique that delivers current through the scalp by a pair of patch electrodes (2-Patch). This study proposes a new multi-channel tDCS (mc-tDCS) optimization method, the distributed constrained maximum intensity (D-CMI) approach. For targeting the P20/N20 somatosensory source at Brodmann area 3b, an integrated combined magnetoencephalography (MEG) and electroencephalography (EEG) source analysis is used with individualized skull conductivity calibrated realistic head modeling.

METHODS

Simulated electric fields (EF) for our new D-CMI method and the already known maximum intensity (MI), alternating direction method of multipliers (ADMM) and 2-Patch methods were produced and compared for the individualized P20/N20 somatosensory target for 10 subjects.

RESULTS

D-CMI and MI showed highest intensities parallel to the P20/N20 target compared to ADMM and 2-Patch, with ADMM achieving highest focality. D-CMI showed a slight reduction in intensity compared to MI while reducing side effects and skin level sensations by current distribution over multiple stimulation electrodes.

CONCLUSION

Individualized D-CMI montages are preferred for our follow up somatosensory experiment to provide a good balance between high current intensities at the target and reduced side effects and skin sensations.

SIGNIFICANCE

An integrated combined MEG and EEG source analysis with D-CMI montages for mc-tDCS stimulation potentially can improve control, reproducibility and reduce sensitivity differences between sham and real stimulations.

Standard Non-Personalized Electric Field Modeling of Twenty Typical tDCS Electrode Configurations via the Computational Finite Element Method: Contributions and Limitations of Two Different Approaches

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation procedure to modulate cortical excitability and related brain functions. tDCS can effectively alter multiple brain functions in healthy humans and is suggested as a therapeutic tool in several neurological and psychiatric diseases. However, variability of results is an important limitation of this method. This variability may be due to multiple factors, including age, head and brain anatomy (including skull, skin, CSF and meninges), cognitive reserve and baseline performance level, specific task demands, as well as comorbidities in clinical settings. Different electrode montages are a further source of variability between tDCS studies. A procedure to estimate the electric field generated by specific tDCS electrode configurations, which can be helpful to adapt stimulation protocols, is the computational finite element method. This approach is useful to provide a priori modeling of the current spread and electric field intensity that will be generated according to the implemented electrode montage. Here, we present standard, non-personalized model-based electric field simulations for motor, dorsolateral prefrontal, and posterior parietal cortex stimulation according to twenty typical tDCS electrode configurations using two different current flow modeling software packages. The resulting simulated maximum intensity of the electric field, focality, and current spread were similar, but not identical, between models. The advantages and limitations of both mathematical simulations of the electric field are presented and discussed systematically, including aspects that, at present, prevent more widespread application of respective simulation approaches in the field of non-invasive brain stimulation.

Does participant's age impact on tDCS induced fields? Insights from computational simulations

Objective: Understanding the induced current flow from transcranial direct current stimulation (tDCS) is essential for determining the optimal dose and treatment. Head tissue conductivities play a key role in the resulting electromagnetic fields. However, there exists a complicated relationship between skull conductivity and participant age, that remains unclear. We explored how variations in skull electrical conductivities, particularly as a suggested function of age, affected tDCS induced electric fields.Approach: Simulations were employed to compare tDCS outcomes for different intensities across head atlases of varying age. Three databases were chosen to demonstrate differing variability in skull conductivity with age and how this may affect induced fields. Differences in tDCS electric fields due to proposed age-dependent skull conductivity variation, as well as deviations in grey matter, white matter and scalp, were compared and the most influential tissues determined.Main results: tDCS induced peak electric fields significantly negatively correlated with age, exacerbated by employing proposed age-appropriate skull conductivity (according to all three datasets). Uncertainty in skull conductivity was the most sensitive to changes in peak fields with increasing age. These results were revealed to be directly due to changing skull conductivity, rather than head geometry alone. There was no correlation between tDCS focality and age.Significance: Accurate and individualised head anatomy andin vivoskull conductivity measurements are essential for modelling tDCS induced fields. In particular, age should be taken into account when considering stimulation dose to precisely predict outcomes.

Effect of skull thickness and conductivity on current propagation for noninvasively injected currents

Objective.When currents are injected into the scalp, e.g. during transcranial current stimulation, the resulting currents generated in the brain are substantially affected by the changes in conductivity and geometry of intermediate tissue. In this work, we introduce the concept of 'skull-transparent' currents, for which the changing conductivity does not significantly alter the field while propagating through the head.Approach.We establish transfer functions relating scalp currents to head potentials in accepted simplified models of the head, and find approximations for which skull-transparency holds. The current fields resulting from specified current patterns are calculated in multiple head models, including MRI heads and compared with homogeneous heads to characterize the transparency. Experimental validation is performed by measuring the current field in head phantoms.Main results.The main theoretical result is derived from observing that at high spatial frequencies, in the transfer function relating currents injected into the scalp to potential generated inside the head, the conductivity terms form a multiplicative factor and do not otherwise influence the transfer function. This observation is utilized to design injected current waveforms that maintain nearly identical focusing patterns independently of the changes in skull conductivity and thickness for a wide range of conductivity and thickness values in an idealized spherical head model as well as in a realistic MRI-based head model. Experimental measurements of the current field in an agar-based head phantom confirm the transparency of these patterns.Significance.Our results suggest the possibility that well-chosen patterns of current injection result in precise focusing inside the brain even withouta prioriknowledge of exact conductivities of intermediate layers.

Neurovascular-modulation: A review of primary vascular responses to transcranial electrical stimulation as a mechanism of action

BACKGROUND

The ubiquitous vascular response to transcranial electrical stimulation (tES) has been attributed to the secondary effect of neuronal activity forming the classic neurovascular coupling. However, the current density delivered transcranially concentrates in: A) the cerebrospinal fluid of subarachnoid space where cerebral vasculature resides after reaching the dural and pial surfaces and B) across the blood-brain-barrier after reaching the brain parenchyma. Therefore, it is anticipated that tES has a primary vascular influence.

OBJECTIVES

Focused review of studies that demonstrated the direct vascular response to electrical stimulation and studies demonstrating evidence for tES-induced vascular effect in coupled neurovascular systems.

RESULTS

tES induces both primary and secondary vascular phenomena originating from four cellular elements; the first two mediating a primary vascular phenomenon mainly in the form of an immediate vasodilatory response and the latter two leading to secondary vascular effects and as parts of classic neurovascular coupling: 1) The perivascular nerves of more superficially located dural and pial arteries and medium-sized arterioles with multilayered smooth muscle cells; and 2) The endothelial lining of all vessels including microvasculature of blood-brain barrier; 3) Astrocytes; and 4) Neurons of neurovascular units.

CONCLUSION

A primary vascular effect of tES is highly suggested based on various preclinical and clinical studies. We explain how the nature of vascular response can depend on vessel anatomy (size) and physiology and be controlled by stimulation waveform. Further studies are warranted to investigate the mechanisms underlying the vascular response and its contribution to neural activity in both healthy brain and pathological conditions - recognizing many brain diseases are associated with alteration of cerebral hemodynamics and decoupling of neurovascular units.

Development of a Non-invasive Deep Brain Stimulator With Precise Positioning and Real-Time Monitoring of Bioimpedance

Methods by which to achieve non-invasive deep brain stimulation via temporally interfering with electric fields have been proposed, but the precision of the positioning of the stimulation and the reliability and stability of the outputs require improvement. In this study, a temporally interfering electrical stimulator was developed based on a neuromodulation technique using the interference modulation waveform produced by several high-frequency electrical stimuli to treat neurodegenerative diseases. The device and auxiliary software constitute a non-invasive neuromodulation system. The technical problems related to the multichannel high-precision output of the device were solved by an analog phase accumulator and a special driving circuit to reduce crosstalk. The function of measuring bioimpedance in real time was integrated into the stimulator to improve effectiveness. Finite element simulation and phantom measurements were performed to find the functional relations among the target coordinates, current ratio, and electrode position in the simplified model. Then, an appropriate approach was proposed to find electrode configurations for desired target locations in a detailed and realistic mouse model. A mouse validation experiment was carried out under the guidance of a simulation, and the reliability and positioning accuracy of temporally interfering electric stimulators were verified. Stimulator improvement and precision positioning solutions promise opportunities for further studies of temporally interfering electrical stimulation.

Accurate and robust whole-head segmentation from magnetic resonance images for individualized head modeling

Transcranial brain stimulation (TBS) has been established as a method for modulating and mapping the function of the human brain, and as a potential treatment tool in several brain disorders. Typically, the stimulation is applied using a one-size-fits-all approach with predetermined locations for the electrodes, in electric stimulation (TES), or the coil, in magnetic stimulation (TMS), which disregards anatomical variability between individuals. However, the induced electric field distribution in the head largely depends on anatomical features implying the need for individually tailored stimulation protocols for focal dosing. This requires detailed models of the individual head anatomy, combined with electric field simulations, to find an optimal stimulation protocol for a given cortical target. Considering the anatomical and functional complexity of different brain disorders and pathologies, it is crucial to account for the anatomical variability in order to translate TBS from a research tool into a viable option for treatment. In this article we present a new method, called CHARM, for automated segmentation of fifteen different head tissues from magnetic resonance (MR) scans. The new method compares favorably to two freely available software tools on a five-tissue segmentation task, while obtaining reasonable segmentation accuracy over all fifteen tissues. The method automatically adapts to variability in the input scans and can thus be directly applied to clinical or research scans acquired with different scanners, sequences or settings. We show that an increase in automated segmentation accuracy results in a lower relative error in electric field simulations when compared to anatomical head models constructed from reference segmentations. However, also the improved segmentations and, by implication, the electric field simulations are affected by systematic artifacts in the input MR scans. As long as the artifacts are unaccounted for, this can lead to local simulation differences up to 30% of the peak field strength on reference simulations. Finally, we exemplarily demonstrate the effect of including all fifteen tissue classes in the field simulations against the standard approach of using only five tissue classes and show that for specific stimulation configurations the local differences can reach 10% of the peak field strength.