Translator Disclaimer
Open Access
25 January 2018 Concurrent measurement of cerebral hemodynamics and electroencephalography during transcranial direct current stimulation
Martina Giovannella, David Ibañez, Clara Gregori-Pla, Michal Kacprzak, Guillem Mitjà, Giulio Ruffini, Turgut Durduran
Author Affiliations +
Transcranial direct current stimulation (tDCS) is currently being used for research and treatment of some neurological and neuropsychiatric disorders, as well as for improvement of cognitive functions. In order to better understand cerebral response to the stimulation and to redefine protocols and dosage, its effects must be monitored. To this end, we have used functional diffuse correlation spectroscopy (fDCS) and time-resolved functional near-infrared spectroscopy (TR-fNIRS) together with electroencephalography (EEG) during and after stimulation of the frontal cortex. Twenty subjects participated in two sessions of stimulation with two different polarity montages and twelve also underwent a sham session. Cerebral blood flow and oxyhemoglobin concentration increased during and after active stimulation in the region under the stimulation electrode while deoxyhemoglobin concentration decreased. The EEG spectrum displayed statistically significant power changes across different stimulation sessions in delta (2 to 4 Hz), theta (4 to 8 Hz), and beta (12 to 18 Hz) bands. Results suggest that fDCS and TR-fNIRS can be employed as neuromonitors of the effects of transcranial electrical stimulation and can be used together with EEG.



Transcranial current stimulation (tCS) is a noninvasive form of brain stimulation that applies weak direct electrical currents to the brain through electrodes placed on the scalp. Different types of stimulation are available based on how the current delivered to the brain is modulated.1 The most common approach is to hold the current constant during the stimulation period, this is known as transcranial direct current stimulation (tDCS). This method is commonly referred to as tDCS, but in order to avoid confusion with the optical method, we here use DC-tCS. DC-tCS induces bidirectional, polarity-dependent changes in cerebral excitability of humans2,3 and has been shown to modulate cognition in both healthy subjects (improvement of cognitive function) and in patients (reversing the effects of cognitive deficits after a stroke).4,5 Positive effects have also been shown in treatment of depression, chronic and acute pain, Parkinson’s disease, focal epilepsy, and for improving recovery after stroke.69 A therapeutic application of this technique is appealing due to its relative low cost, to its demonstrated safety,10 and as a substitute for pharmacotherapy, especially for patients with poor drug tolerability.8

In order to define protocols of application and to optimize the individual dosage of stimulation, the effects of the stimulation on the cerebral activity must be monitored in real time.11,12 Availability of a concurrent read-out of the stimulation effects could improve its effectiveness and, ultimately, a personalized stimulation driven by the neuromonitor feedback can be imagined.13 For this purpose, due to the complex nature of cerebral activity, multimodal monitoring is preferable,14,15 such that different monitors integrated in a single set-up can allow to record the effects at different levels, from neuronal activity, to hemodynamics to systemic.16 The neuronal activity is evaluated by electroencephalography (EEG), which measures the synchronous activation of a large population of pyramidal neurons oriented perpendicularly to the scalp. On the other hand, cerebral hemodynamics can be used as a surrogate measure of the cerebral activation assuming that cerebral blood flow (CBF) increases in a region of the brain, where neurons and synapses are activated in order to meet the demand for more energy and that the amount of local increase of blood flow and the changes in oxygen saturation are related to the cellular energy consumption.17

Even if different options are available as neuromonitors of hemodynamics, functional near-infrared spectroscopy (fNIRS) appears to be the best suited for this application.18 The most accessible and commonly employed fNIRS technology uses continuous-wave sources (CW-fNIRS), where constant intensity light sources are employed to measure oxy- (HbO2) and deoxyhemoglobin (Hb) concentration changes as surrogate measures of CBF. fNIRS has advantages that are fundamental for this application. In particular, it allows for continuous noninvasive measurements, does not require the immobility of the subject, and is not necessary to change the environment preferable for the stimulation to integrate CW-fNIRS in the protocol.18 It has indeed already been used as a monitor during DC-tCS, both in animals19 and on humans at rest2022 and during tasks.2224 It has also been used together with EEG to monitor DC-tCS25 and a device integrating the two monitors has been developed and characterized.26,27 Despite of the successes of these studies, it must be noted that CW-fNIRS does not allow a quantification of absolute value of HbO2 and Hb concentration, lacks precision and accuracy,28 and does not give a direct measure of CBF, the primary marker of the neurovascular coupling.

In this work, we propose to use time-resolved fNIRS (TR-fNIRS) and functional diffuse correlation spectroscopy (fDCS) as neuromonitors of the brain stimulation. Using pulsed light of few hundreds of picoseconds in width, TR-fNIRS measures absolute values of Hb and HbO2 concentration by separating the effects of absorption and scattering. The time information can also be used to explore different behaviors in different layers of the tissue.29 On the other hand, fDCS30 allows a direct measure of CBF in the microvasculature of the tissue.

Here, fDCS and TR-fNIRS are used together with EEG to monitor cerebral activity before, during, and after 10 min of DC-tCS over the frontal cortex. This area is the most convenient choice for the diffuse optical monitors because of the lack of hair. Since its stimulation has been proven effective for the improvement of various cognitive functions,3138 for example, working memory enhancement,3942 and for treatment of depression,9,43,44 it is also a relevant area to stimulate and study.




Transcranial Direct Current Stimulation and Electroencephalography Device

DC-tCS was delivered using Neuroelectrics Starstim® (Neuroelectrics, Barcelona, Spain), which is a research-class eight channel transcranial current stimulator that is also capable of measuring EEG. About 12-mm diameter silver chloride electrodes (Ag/AgCl) with conductive saline gel were used for EEG recording and for stimulation. Eight electrode (AF7, AF8, FT7, FT8, TP7, TP8, PO7, and PO8) locations were chosen, as depicted in Fig. 1(a), according to the 10–10 EEG system.45 In particular, the active electrodes during stimulation were AF7 (left frontal lobe) and PO8 (right parieto-occipital lobe). From now on, we will refer to the former (AF7) as the stimulation electrode while to the latter (PO8) as the return electrode. All the eight electrodes functioned as EEG recording electrodes for protocol periods that did not involve stimulation. EEG spectrum was sampled at 500 Hz.

Fig. 1

(a) A schematic of the optical probes, stimulating and return electrodes and EEG recording electrodes superimposed on a standard 10–10 EEG system. (b) Image of the integrated probe that could accommodate optical fibers (both for fDCS and TR-fNIRS) and an electrode. (c) Schematic of the protocol that lasted a total of 50 min divided in five 10-min periods. Stimulation was delivered during the second period.



Diffuse Optical Monitors

Two diffuse optical monitors were employed to monitor cerebral hemodynamics during and after the stimulation. A commercial fDCS device (HemoFloMo, HemoPhotonics S.L., Castelldefels, Spain) with four source (785 nm) and eight parallel detector channels was used to measure CBF in the microvasculature. In addition, a prototype of a TR-fNIRS system (TRS-20, Hamamatsu Photonics K.K., Hamamatsu City, Japan) was used for TR-fNIRS. It works at three wavelengths (760, 800, and 830 nm) and has two independent source–detector channels. Hardware triggers allowed the communication between the two optical neuromonitors in order to measure alternately without interfering with each other. The procedure was fully automatized with fDCS defined as the “master.” This means that fDCS triggered the start of the acquisition of TR-fNIRS and then paused the acquisition until a signal was received from TR-fNIRS. Furthermore, the marker signals were introduced by the operator on the fDCS device, which sent a signal to TR-fNIRS in order to have the markers registered in the TR-fNIRS.

Since the measurement of the two devices was sequential, the temporal resolution of the hemodynamics parameters was the sum of the fDCS averaging time (1.5 s) and the TR-fNIRS integration time. The TR-fNIRS averaging time was adjusted between 3 and 4 s for each subject in order to gather enough photons to have three decades between the background noise and the peak of the collected TR-fNIRS distribution of the time-of-flight (DTOF) of photons.

Since fDCS and TR-fNIRS were both meant to probe the area under the corresponding electrode, an integrated probe holder was designed and 3-D printed with a flexible material. As shown in the image in Fig. 1(b), the flexible probe could accommodate fibers for fDCS (source–detector separation 2.5 cm), for TR-fNIRS (source-detector separation 3.7 cm), and an electrode. Two of these probes were placed on the forehead of the subject, one on the left hemisphere in correspondence with AF7, carrying the stimulation electrode, while the other on the right hemisphere in AF8 position with a purely EEG recording electrode. Two of the four sources available for fDCS were used, one for the left and one for the right hemisphere. For each hemisphere, only one detection position was implemented using a bundle of four single-mode fibers probing the same area, i.e., all the eight detector channels were used. The correlation curves collected by four fibers in the same area were averaged in order to improve the signal-to-noise ratio. On the other hand, the two independent source–detector channels of the TR-fNIRS device were used for the two hemispheres.



Healthy, right-handed adult subjects were recruited for three sessions of combined stimulation and measurement: anodal [i.e., anode in AF7 (stimulation electrode) cathode in PO8 (return electrode)], cathodal (i.e., cathode in AF7, anode in PO8) stimulation sessions, and a sham session. Such a montage creates an electric field that crosses the brain as confirmed by simulation of the electric field distribution in the brain (data not shown). For half of the subjects, the anodal was the first session, while for the other half cathodal was. Sham stimulation was always given during the third session. All sessions were separated by a wash-out period of at least 1 week. The studies were approved by the ethical committee of Hospital Clinic in Barcelona. Each subject signed an informed consent and the study was conducted according to the principles of the Declaration of Helsinki. Exclusion criteria included the history of neurological and psychiatric conditions and the current use of psychoactive medications.

A schematic of the protocol is presented in Fig. 1(c). The protocol consisted of 50-min measurement sessions divided into five periods. After 10 min of baseline (pre-period), 10 min of stimulation was given (during period), then 30 more minutes were recorded (post1, post2, and post3 periods). During the stimulation, a current of 1 mA (resulting in a current density of 0.88  mA/cm2) was delivered to the brain with a ramp up of 30 s at the beginning and a ramp down of 30 s at the end of the stimulation period. Such a duration for ramp up and down was chosen considering that a shorter one could cause an itching and disturbing sensation for the subject. In the sham session, the current was ramped up to 1 mA in 30 s and then brought back to 0 in 30 s. After 8 min, it was again ramped up and down to simulate the real stimulation sessions. This design gave the subjects the same feeling as a real stimulation because the ramp up and down periods are usually the most noticeable moments. However, the stimulation was too short to have significant, durable effect on the brain activity.

For each 10-min period, the subject was asked to close the eyes for the first 2 min (from 0 to 2, from 10 to 12, 20 to 22, 30 to 32, and 40 to 42 min). During the eyes-open condition, subjects were requested to stare at a screen, where a cross was shown to change color or to rotate at random time (but between every 20 or 30 s) and to press a key when the cross configuration changed in order to maintain the subject gaze and attention to a fixed point.


Data Analysis

For TR-fNIRS analysis, three different approaches were used. (1) First, the so-called differential pathlength factor (DPF) analysis was implemented.46 Total count rate and the modified Beer–Lambert’s law were used to calculate the absorption changes and consequently the changes in Hb and HbO2 concentration in order to simulate what a typical CW-fNIRS device would have registered. (2) Second, the DTOF of photons collected after propagation in the tissue was fit to the solution of the photon diffusion equation for a semi-infinite homogeneous medium for the reflectance in the time-resolved regime.47 A measure of optical properties (absorption and reduced scattering coefficient) for each wavelength was derived and an absolute measure of Hb and HbO2 concentrations was calculated. This option gave more accurate results with respect to the DPF analysis because complete information on the tissue optical properties was retrieved. To improve the stability of derived optical properties, the reduced scattering coefficient was kept fixed to the mean of the baseline value (pre-period) during the whole protocol. This approximation assumed that metabolic changes affect only the blood oxygenation but not the tissue scattering properties. (3) The last approach included an analysis of the second central moment (variance) of the DTOF. Changes in the absorption coefficient were derived from changes in the variance.48 Through Monte Carlo simulations of photon propagation into the tissue, it was demonstrated that the variance is more sensitive to changes in the deeper layers compared to the intensity (used for DPF analysis), which is more sensitive to the superficial layer. For this reason, we have compared results from DPF and variance analyses.

DCS measurements were analyzed using the solution for the diffusion equation of the electric field autocorrelation function for the semi-infinite medium30 to derive the blood flow index (BFI). Optical properties derived from TR-fNIRS were used in the analysis. Once derived, the BFI, a measure of CBF, was normalized to the baseline (first 10 min of acquisition) to obtain a measure of relative CBF (rCBF).

Changes in the EEG patterns induced by the brain stimulation were independently analyzed for eyes open and closed conditions. Pre- and post-stimulation EEG recordings were filtered using a finite impulse response filter of 500 coefficients with cut-off frequencies set to 1 and 20 Hz. Filtered recordings were split into 60-s blocks in which power at delta (2 to 4 Hz), theta (4 to 8 Hz), alpha (8 to 13 Hz), and beta (12 to 18 Hz) bands were calculated. Each 60-s block was further split into 2-s windows with a 50% overlap. After individual epoch detrending, those epochs containing samples larger than 75  μV, considered high amplitude artifacts, were rejected. Relative power at the previously defined frequency bands was calculated via trapezoidal power spectral density integration and was normalized with respect to the full band power (2 to 18 Hz). Power at 60-s blocks was calculated as the average power of artifact-free 2-s epochs.


Statistical Data Analysis

Outlier detection was implemented by analyzing the functional depth,49,50 i.e., the centrality of a given curve within a group of trajectories, with each curve representing the response of cerebral hemodynamics or EEG power. The R package51 “fda.usc”52 and the R function “Outliergram”50 were used as tools for this purpose. The outlier detection was run in the two hemispheres independently for the cerebral hemodynamics and for each electrode and band for the EEG response. A subject was defined as an outlier if it was detected by both methods simultaneously and was, thus, removed from the statistical analysis. If a subject was defined an outlier for HbO2 concentration, then it was removed also for Hb concentration, and vice versa.

A linear mixed effects (LME) analysis was performed to test the effect of the stimulation session and time on the cerebral hemodynamics and EEG response using the R package “lme4.”53

For cerebral hemodynamics, we treated the two hemispheres independently building an LME model for each of them. To include time, we selected periods of 4 min during which we averaged the response: minute 14 to 19 as during DC-tCS period, minute 24 to 29, 34 to 39, 44 to 49 as post1, post2, and post3 DC-tCS periods, respectively. Response was normalized to the last 4 min of the pre-period; therefore, this period was not introduced in the analysis.

Proceeding to EEG power change, with respect to initial baseline, each band and electrode was evaluated separately, as well as eyes opened and eyes closed periods. Post-periods were introduced in the LME by averaging the response in each of them.

In all the models, the period and the stimulation type (and their interaction) entered as fixed effects. As random effects, the intercepts for subject and the random slope for period on a subject basis were considered. For EEG, only a random intercept was used since the model would not converge if the random slope had also been added. Residual plots were tested for the inspection of deviations from homoscedasticity or normality. Models were built with forward steps following likelihood ratio tests of the model with the effect in question against the model without the effect in question. The Bayesian information criterion parameter was considered to assess whether the effect was improving the model. Stimulation type was the first to be tested. If found significant, the period was added and a new model was tested. If this step was also successful, the interaction between the two was also tested.

Posthoc analysis was performed with R package “lsmeans”54 to quantify changes in hemodynamics parameters and in EEG power for each stimulation session and period. To this end, least-square means of the linear model were calculated and the different sessions were compared using the Tukey correction for multiple comparisons. For data visualization, a bootstrap sample was calculated for the time series of cerebral hemodynamics parameter. The bootstrap method is based on consecutive and random resampling of the sample distribution55 and was implemented in the R package “fda.usc.”52



A total of 20 healthy subjects (nine female) participated in anodal and cathodal stimulation sessions. Twelve of them came back for a sham session. Technical problems resulted in one subject for TR-fNIRS measurement being discarded for the cathodal session leaving 19 subjects for HbO2 and Hb concentration measurements. In addition, one fDCS measurement was discarded during the sham session leaving 11 subjects for flow measurement during sham stimulation. Technical exclusion criteria for EEG, explained in Sec. 2.4, resulted in 20 subjects for anodal, 16 for cathodal, and 9 for the sham session. As explained above, outliers were detected and excluded from the statistical analysis. The number of subjects surviving the outlier detection for EEG and hemodynamics parameters is presented in Tables 1 and 2, respectively.

Table 1

The number of subjects for EEG that were used for the statistical analysis for each electrode and for each band, after the technical exclusion criteria and outlier detection.


Table 2

The number of subjects for hemodynamics parameters that were used for the statistical analysis after the technical exclusion criteria and outlier detection.

HbO2 and Hb121218191717

In Fig. 2, the time series of rCBF change, averaged over subjects, is shown together with the bootstrap sample. Only results from the left hemisphere (ipsilateral side, under the stimulation electrode) are shown from the three different stimulation sessions.

Fig. 2

rCBF change during the three types of stimulation sessions on the left hemisphere. (a) Sham session, (b) anodal session, and (c) cathodal session. The thick black line is the average, and the gray area represents the bootstrap sample over all the subjects. The black horizontal bar highlights the stimulation period.


When building the LME model from the null one, the stimulation type is found to improve it for blood flow in the left hemisphere, but not in the right one. Therefore, blood flow response depends on the stimulation in the ipsilateral hemisphere, but not in the contralateral one. Period is not an effect that improves the model for either of the hemispheres. Accordingly, we can conclude that blood flow stays constant at the end of the stimulation period (during) and in the post-periods. Posthoc analysis results are summarized in Fig. 6. Blood flow increases significantly from the baseline by 10% [confidence interval (CI): 7% to 14%] for anodal (p<0.001) and 11% (CI: 7% to 15%) for cathodal (p<0.001) stimulation while no increase is seen in the sham session (p=0.07). A comparison between the different sessions revealed a detectable difference between anodal and sham (p=0.008) and cathodal and sham (p=0.003) sessions while no difference is seen between anodal and cathodal sessions.

In Fig. 3, the results for the time evolution of HbO2 concentration from TR-fNIRS measurements using the three different analysis methods are shown for anodal and sham stimulations for the left (ipsilateral) hemisphere. DPF analysis gives an increase of HbO2 concentration starting during the stimulation and persisting after it. Notably, this change is seen both in the anodal and sham sessions. The results obtained with the analysis of variance and by fitting are comparable and show smaller changes than the DPF analysis both in anodal (p<0.001) and sham (p<0.001) sessions.

Fig. 3

Changes in oxyhemoglobin (HbO2) concentration during (a) sham and (b) anodal stimulation, in the left hemisphere, according to the three types of analysis considered for TR-fNIRS: orange for DPF analysis, blue for the variance, and red for the fitting analysis. Continuous line is the average of all the subjects, the lighter colored region represents the bootstrap sample with the same color code. The black horizontal bar highlights the stimulation period.


For this reason, the fit analysis results are used for building the LME model for the changes in HbO2 and Hb concentrations. Furthermore, this is what is most commonly used for TR-fNIRS analysis.29 Figures 4 and 5 show the average time evolution and the bootstrap sample from all the subjects for HbO2 and Hb concentrations, respectively, in the left hemisphere during the three stimulation sessions.

Fig. 4

Oxyhemoglobin (HbO2) concentration during the three types of stimulation sessions on the left hemisphere from the fit analysis. (a) Sham session, (b) anodal session, and (c) cathodal session. The thick line is the average, and the lighter colored region represents the bootstrap sample over all the subjects. The black horizontal bar highlights the stimulation period.


Fig. 5

Deoxyhemoglobin (Hb) concentration during the three types of stimulation sessions on the left hemisphere from the fit analysis. (a) Sham session, (b) anodal session, and (c) cathodal session. The thick line is the average, and the lighter colored region represents the bootstrap sample over all the subjects. The black horizontal bar highlights the stimulation period.


As for the blood flow, among the effects considered in the LME model and listed in Sec. 2, only the stimulation improves the model. For both variables, the changes depend on the stimulation session but a recovery is not seen in the 30 min of post-periods, because the period is not a significant effect for the model.

As summarized in Fig. 6, HbO2 concentration increased by 0.5  μM (CI: 0.1 to 0.8  μM) during the anodal session (p=0.002) and by 0.5  μM (CI: 0.2 to 0.8  μM) during the cathodal (p=0.001) session. As expected, the changes in Hb concentration were opposite to HbO2 concentration. Hb concentration decreased by 0.3  μM (CI: 0.5 to 0.1  μM) during the anodal session (p<0.001) and by 0.3  μM (CI: 0.5 to 0.2  μM) during the cathodal session (p<0.001). No change was seen in the sham session (p=0.1 for HbO2 and p=0.07 for Hb concentration). A difference was detected by our model for HbO2 concentration between anodal and sham (p=0.004) sessions, as well as between cathodal and sham (p=0.003) sessions. On the other hand, a pair comparison for Hb concentration gave a significant difference only between cathodal and sham sessions (p=0.03).

Fig. 6

Changes of cerebral hemodynamics parameters during and after the stimulation sessions, as obtained by the LME model with the 95% confidence interval as the errorbars. N refers to the number of subjects used for the statistical analysis. Red asterisk highlights statistically significant difference between stimulation sessions. Dashed vertical line represents the zero level.


A frequency range from 1 to 20 Hz was considered for the EEG analysis. The lower limit was chosen because artifacts derived from movement and subsequent changes in the skin-electrode impedance affect mostly the low EEG frequencies. Furthermore, a time–frequency analysis of each recording was performed in order to study the signal quality of the dataset. In many recordings, the EEG spectrum was affected by external interferences in frequencies above 22 Hz. It was then decided to choose an upper cutoff frequency of 20 Hz to avoid this noise. As a result, the beta band covers 13 to 20 Hz.

In addition, after the temporal and spectral inspection of the EEG spectrum that was acquired during the stimulation, we have concluded that the quality was not sufficient to perform a reliable EEG analysis. Therefore, only recordings before and after the stimulation were used.

Going through results obtained in EEG power, different changes from resting-state across stimulation modalities were observed only in the eyes-open condition. Figure 7 shows barplots for bands and electrodes where response was affected by the stimulation session according to LME model. Posthoc analysis resulted in a statistically significant power decrease in delta band in TP8 after every stimulation type. Statistically significant differences between anodal and cathodal sessions (p<0.001) and anodal and sham sessions (p=0.02) were observed. A statistically significant power decrease in the theta band after cathodal stimulation in fronto-right region (AF8, FT8) was observed. Statistically significant differences between cathodal and sham sessions (p=0.002 in both electrodes) and cathodal and anodal sessions (p=0.002 in both electrodes) were also observed. A more widespread effect in the beta band was observed. All electrodes in the tempo-parietal (TP7 and TP8) and parieto-occipital (PO7 and PO8) lobes experienced a statistically significant increase after sham, but, not after any active stimulation. In all these electrodes, the difference between anodal and sham sessions (p=0.004 for TP7, p=0.01 for TP8, p<0.001 for PO7, p=0.002 for PO8) and cathodal and sham sessions was significant (p<0.001 for all electrodes).

Fig. 7

Change of EEG power in electrode TP8 in the delta band, electrodes AF8 and FT8 in the theta band, and electrodes TP7, PO7, TP8, PO8 in the beta band caused by the three types of stimulation, as obtained by the LME model, with 95% confidence interval as the errorbars. N refers to the number of subjects used for the statistical analysis. Red asterisk highlights the statistically significant difference between stimulation sessions. Dashed vertical line represents the zero level. Only bands and electrodes with statistical significant changes are shown.




In this work, we have introduced the use of fDCS, together with TR-fNIRS, as neuromonitors of cerebral hemodynamics during and after DC-tCS. The former has allowed the direct measurement of CBF primary biomarker of the neurovascular coupling. The measurement of cerebral hemodynamics was concurrent to EEG, which measures electrical activity. Therefore, information on both hemodynamics and neuronal activity from the same run of the protocol were available.

Before going through the results of the hemodynamics parameters, we note that the fDCS analysis requires the optical properties of the probed tissue as input parameters. In this experiment, the optical properties were measured by TR-fNIRS with fibers with a source–detector separation of 3.7 cm. On the other hand, source–detector separation of fDCS fibers was 2.5 cm, hence, probing a different volume than TR-fNIRS. Since our analysis assumes that the underlying tissue is homogeneous, we can tolerate this difference. It must also be highlighted that since we focus on the measurement of relative CBF, our results are relatively insensitive to small inaccuracies of the bulk optical properties.56

Statistical analysis confirms that CBF increases in the left ipsilateral hemisphere (in the region under the stimulation electrode) during anodal and cathodal sessions but not during the sham session. CBF does not recover during the 30 min following the stimulation. This is consistent with the fact that if the stimulation is applied for a period longer than 9 min, the effect of the stimulation persists long after the stimulation ends.57 Specifically, our statistical model results in a 10% increase during the anodal session and 11% during cathodal and gives no difference between the two stimulation conditions. This amount for CBF change is similar to what was assessed with arterial spin labeling, a magnetic resonance imaging technique,58 but, is much larger than what measured by positron emission tomography (PET).59 In Ref. 59, scans taken during the following hour after the stimulation were averaged, while we have measured only the following 30 min, which may account for some of the difference in our results. The optical probe placed in the right contralateral hemisphere did not detect any change in CBF in this region. The work done with PET59 allowed for whole brain imaging and detected regional increases in different areas of the brain, even far away from the stimulation electrodes, after both anodal and cathodal sessions compared to sham. Although anodal stimulation proved to cause a more widespread increase of regional CBF than cathodal stimulation, the same effect was seen for both stimulation types in the region under the electrode, which is consistent with our results.

It is known that diffuse optical techniques measure a mixture of extra- and intracerebral perfusion.60 This could be a confounding factor since DC-tCS causes a “skin redness” (erythema), which may lead to changes in the skin perfusion of different amounts between active and sham stimulations.10 However, the properties of photon diffusion into the tissue help us in explaining that response detected in this protocol is not mainly due to changes in skin perfusion due to the stimulation induced skin redness. The visitation probability of photons traveling into the tissue has a so-called “banana-shaped” pattern and touches the skin mainly in the region corresponding to the source and detector fibers.61 In our case, the fibers are separated by a distance of 2.5 cm while the electrode, placed in the middle, has a diameter of 1.2 cm. If we consider that various works62,63 have shown that the erythema is concentrated in the area under the electrode, our optical monitors are not particularly sensitive to this redness. Furthermore, a recent study64 measured the temperature under the electrode during the stimulation, detecting a nonsignificant change after 10 min of stimulation (our total time of stimulation). Only after 20 min of application, a significant increase of 1.36 °C was detected. However, we should account for the fact that this work used a larger electrode resulting in a lower current density than ours (0.06 versus 0.88  mA/cm2). Interestingly, there is evidence suggesting that the current intensity and the electrode area are more important parameters leading to increased redness and discomfort due to the stimulation of more cutaneous receptors.62 In the above mentioned Ref. 64, temperature was measured using a current intensity of 2 mA, larger that what was used in our experiment. Therefore, we expect a minimal, nonsignificant temperature increase in our case. Another indication that the observed, long-term changes are not purely due to skin redness is that it is known that erythema disappears within 30 min after the stimulation.65 In our case, the stimulation response did not decrease after 30 min.

Proceeding to TR-fNIRS measurements, results from the different analysis methods considered were compared. To start with, DPF analysis gives us the results that would have been obtained using a typical CW-fNIRS setup for fNIRS, the most commonly used fNIRS technology. DPF analysis shows an increase of about 2  μM in HbO2 concentration in the left hemisphere compatible with what was previously found in the literature.20,22 This increase is seen in active and sham stimulations, where the brain is not actually stimulated which questions its validity. Focusing on the results given by the analysis of variance, the second moment of the DTOF, we retrieve smaller changes in both anodal and sham sessions with respect to DPF analysis. Notably, the variance of TR-fNIRS curves is more sensitive to the deeper layers with respect to the intensity, the parameter used in the DPF analysis.48 This may suggest that results in DPF analysis are contaminated by superficial effects present both during anodal and sham stimulations. Moreover, these results are in accordance with a previous single subject experiment, where TR-fNIRS was used during anodal stimulation.21 In this work, TR-fNIRS measurements were analyzed with a gated analysis.66 This method exploits the depth information encoded in time in the broadened pulse collected after propagation into the tissue and allows to decouple changes in the extracerebral and cerebral layers. Measurements performed on a single subject during anodal stimulation retrieved a larger HbO2 concentration change in the superficial than in the deeper layer, which would be heavily reflected in the DPF analysis. We would like to point out that the time gating analysis of TR-fNIRS was not implemented here because the lack of knowledge of the extracerebral layer thickness on an individual basis compromises the validity of this method.67

The usage of TR-fNIRS has allowed us to improve the separation of the extracranial and intracranial signals for the analysis of the hemoglobin concentration changes. We do not have a comparable time-domain measurement for fDCS at the moment.68,69 However, it has been proved that the relative brain-to-scalp sensitivity is about three times higher for fDCS than for CW-fNIRS.70 Therefore, we believe that our results are reflective of the CBF in a more reliable manner than the CW-fNIRS results.

Since the fit analysis gives results comparable to the variance, and this is the most commonly used method when measuring oxygenation with TR-fNIRS, we have decided to consider the results given by the fit analysis.

In this protocol, HbO2 concentration increases and Hb concentration decreases in the left hemisphere during both anodal and cathodal stimulations, but not during sham. The right frontal lobe (contralateral to the stimulation) does not present changes in the oxygenation as is also true for CBF. Even if changes retrieved were very small, less than 1  μM, HbO2 and Hb concentrations can mark the status of the ongoing stimulation. As a matter of fact, the changes are statistically significant during both anodal and cathodal sessions but not during the sham session. HbO2 concentration is a more robust marker, since response to anodal and to sham sessions is statistically significantly different, as well as the response to the cathodal and to sham sessions. On the other hand, for Hb concentration, the contrast-to-noise ratio is not enough to give a statistically significant difference between anodal and sham sessions, but a decrease is not seen in the sham session.

The directions of Hb and HbO2 concentration changes are consistent with what is expected in the case of a cerebral activation and what was previously found in the literature.20,22 Specifically, due to the neurovascular coupling, when a cerebral region is activated, CBF is recruited to meet the additional demand of oxygen and it increases.17 Concurrent to this, an increase of HbO2 concentration and a decrease of Hb concentration are usually detected.71 This suggests that the brain is over supplied by oxygen by the increase of CBF compared to the oxygen consumption. The details of this process are, in fact, more complex.17 The time series of the hemodynamics parameters in Figs. 2, 4, and 5 suggest a steeper increase of CBF than the HbO2 and Hb concentrations. The above-mentioned mechanism of neurovascular coupling can explain this.71 Right after the stimulus, the supply of CBF does not overcome the need immediately. Only after this happens, a HbO2 concentration increase and Hb decrease is detected. This process was previously verified during in vivo protocols.72 Nonetheless, to make quantitative inferences on the time lag, we should have had a better time resolution than what we have here.

Both the anodal and cathodal stimulation sessions have created similar changes in the hemodynamics parameters. All the reviewed and cited published works58,59 that measure cerebral hemodynamics during DC-tCS show the same pattern for both cathodal and anodal stimulations in the region under the electrode. Our results confirm these findings. In addition, PET studies concurrent to transcranial magnetic stimulation also demonstrate an increase in CBF in the region under the stimulation for both excitatory and inhibitory protocols.73 These results suggest that CBF is not a marker of the polarity of the stimulation, at least in the directly stimulated area, and that regardless of whether the stimulation results in an activation of inhibitory or excitatory synapses, it leads to an increase in blood flow. It has been highlighted that this behavior of CBF reflects the local levels of synaptic activity in intracortical neurons and inputs to those areas rather than the activity of output neurons, which is opposite when the stimulation is excitatory or inhibitory.59,74

The stimulation protocol used in this study influenced EEG rhythmic activity at different bands and electrodes. Comparing the results from previously published studies about DC-tCS on the frontal cortex is not a trivial task since EEG response depends heavily on the electrode montage.75 All of the previous studies used F3 as the stimulation electrode while the return electrode, here referring to the electrode not placed in the frontal cortex with opposite polarity, was placed either on the supraorbital region or on extracerebral regions. By contrast, we have decided to use the AF7 electrode as stimulation electrode because the optical measurement series on the frontal brain region are more reliable due to the absence of hair. The return electrode was chosen to be spatially and functionally distant from the stimulation electrode in order to minimize the interactions and to stimulate the frontal region unilaterally.

Moreover, few studies measured the EEG spectrum at rest while most of them assessed the response to a task and how this can be influenced by the stimulation. It must be noted that a recent systematic review highlighted that results in the EEG spectrum of the same protocol were not consistent within different works.76 They concluded that there is no evidence that the EEG spectrum is affected by the stimulation. It may also be hypothesized that, since it is not trivial to maintain subjects in a well-controlled condition, the variability of the effects is increased.

We note that the protocol is 50 min long, which may have led to a sense of drowsiness and/or fatigue in the subjects, which could alter the cerebral electrophysiology and metabolism. The staring and attention tasks were implemented to minimize this effect. Furthermore, the transition from eyes open and closed and vice versa also helped with this issue. In spite of this, we cannot exclude that some of the subjects may have been affected and we do not have any control over this.

In order to define EEG spectrum results, effects in different frequency bands have been studied since different bands are traditionally related to different brain functions. It is, however, controversial whether the activity at a certain frequency band is related to a single brain function or that complex stimuli are reflected on a single oscillatory activity.77 This may explain the multiplicity of bands affected by our protocol.

The protocol applied in this study influenced brain activity in the delta, theta, and beta bands. No effect was detected in the alpha band, which is consistent with previous works.76

We will focus, from now on, on the bands, where the stimulation has proven to be effective. The changes in the delta band were observed in the right tempo-parietal region after all stimulation sessions but were significantly smaller in the anodal session compared to the other two. On the other hand, changes in the theta band occurred in the right fronto-temporal area exclusively after cathodal stimulation. Interestingly, delta and theta bands are associated with functions linked to the frontal cortex, which was stimulated in our protocol. It is an indication of activity in cognitive and memory performances,78 attention,79 and working memory.80 A response in these bands is consistent with previous work. Specifically, a lower power in delta and beta bands after the stimulation compared to sham has been previously observed,35,81 where pre- and post-stimulation spectral power was not calculated but the comparison between stimulation sessions was done in absolute values. This effect has been found only in the frontal cortex, in contrast to our findings, but both stimulation and return electrodes were placed on the anterior part of the head, while we have an electric field crossing all the brain. A significant and selective decrease of the power of theta band in the right inferior frontal gyrus area has been reported after DC-tCS stimulation along with an improvement in behavioral inhibition.82

Differences in EEG power across different stimulation sessions were also seen in the beta band. Power increased after the sham session in temporal and parietal regions, but not after active stimulation sessions. An unchanged or a decreased level of beta activity is usually observed in tasks where behavior is driven by bottom-up signals, while enhancement occurs if the system has to actively maintain endogenous attention.83 It may be hypothesized that the detected increase in beta activity is linked to an active effort to maintain the cognitive concentration state required for the sustained attentional task along the experimental protocol. This activity was present during sham but not during active stimulation sessions. Several references highlight an influence of DC-tCS when applied on the frontal cortex, on sustained attention.84 In our protocols, the same effects for anodal and cathodal sessions were found in the EEG beta band in the posterior part of the brain, without a differentiation due to the polarity. Intuitively, this is not expected because opposite polarity stimulation should produce an opposite electric field in the brain, i.e., either excitatory or inhibitory. Nonetheless, it must be mentioned that the predominance of excitatory or inhibitory electric field in an area of the brain does not exclude that some “pockets” of the opposite direction are present.85 We can speculate that only excitatory areas influenced the power in the beta band in the posterior area of the brain in our protocol. Unfortunately, this is prohibitively hard to confirm since not only the electrode montage but also the concurrent tasks alter the response to a certain polarity of stimulation,86,87 i.e., it would require a special protocol to address this. For our protocol, the attention task was not the primary goal.

Last, it was verified whether the change in CBF in the left frontal cortex due to stimulation has any correlation with the power change respect to baseline in any of the EEG bands and electrodes. This was done even for the electrodes far away from the optical probes since the source of the EEG signal can lay far from the scalp point, where it is registered and its exact location cannot be well determined due to the inverse problem.88 The Pearson correlation coefficient was calculated and a correlation was defined significant if, by bootstrapping the sample (i.e., repeating the analysis removing one subject at time), the correlation stays significant. After this analysis, we have found a negative correlation in the delta band for electrode FT8 with a Pearson correlation coefficient of R=0.6. This is consistent with a predicted and verified decrease in lower frequency band during activation.89,90 The study of the correlation between the hemodynamics and neuronal activity is a very appealing topic due to the important information on the neurovascular coupling it could provide. It is also a very complicated process since several factors can contribute to the failure of the correlation between EEG activity and hemodynamics.91 A detailed characterization of this is beyond the scope for our paper, which aimed at the introduction of hybrid diffuse optical devices for monitoring the response to the stimulation, and requires a dedicated protocol with more regions probed for hemodynamics.

In summary, we have introduced fDCS as a neuromonitoring technique for following the response to DC-tCS. It was integrated with TR-fNIRS to measure the concentrations of HbO2 and Hb in a more accurate way with respect to CW-fNIRS. We have proved that fDCS and TR-fNIRS can be applied concurrently to EEG measurements, which assess cerebral activity. We were able to have information on both cerebral hemodynamics and neuronal activity in the same run of stimulation.



We have shown that fDCS and TR-fNIRS are suitable tools to monitor cerebral hemodynamics during and after DC-tCS and they can be integrated with EEG, which monitors the neuronal activity. We have found CBF to be a good indicator of the ongoing stimulation since it showed an increase during and after anodal and cathodal stimulation in the region under the stimulation electrode. Results obtained with TR-fNIRS showed small changes of HbO2 and Hb concentrations but the effect of the stimulation was still detectable.

By concurrent EEG recordings, we could follow the modulation of the underlying neural oscillations. DC-tCS over the frontal cortex induced statistically significant power changes across different stimulation sessions in delta, theta, and beta EEG rhythms.


ICFO has equity ownership in the spin-off company HemoPhotonics S.L. Potential financial conflicts of interest and the objectivity of research have been monitored by ICFO’s Knowledge & Technology Transfer Department. No financial conflicts of interest were identified. Neuroelectrics is the manufacturer of the Starstim system.


The study was funded by Fundació CELLEX Barcelona, Ministerio de Economía y Competitividad/FEDER (PHOTODEMENTIA, DPI2015-64358-C2-1-R), Instituto de Salud CarlosIII / FEDER (MEDPHOTAGE, DTS16/00087), the “Severo Ochoa” Programme for Centers of Excellence in R&D (SEV-2015-0522), the Obra social “la Caixa” Foundation (LlumMedBcn), AGAUR-Generalitat (2014SGR-1555) and LASERLAB-EUROPE IV (654148). M. G. is funded by the “Severo Ochoa” Programme for Centers of Excellence in R&D (SEV-2015-0522) and BabyLux (620996 CIP-ICT-PSP-2013-7). We thank Hamamatsu Photonics K.K. and its Spain office for loaning the TRS-20 system and for their collaboration.



G. Ruffini et al., “Transcranial current brain stimulation (tCS): models and technologies,” IEEE Trans. Neural Syst. Rehabil. Eng., 21 (3), 333 –345 (2013). Google Scholar


M. A. Nitsche and W. Paulus, “Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation,” J. Physiol., 527 (3), 633 –639 (2000). JPHYA7 0022-3751 Google Scholar


A. Priori, “Brain polarization in humans: a reappraisal of an old tool for prolonged non-invasive modulation of brain excitability,” Clin. Neurophysiol., 114 (4), 589 –595 (2003). CNEUFU 1388-2457 Google Scholar


V. P. Clark and R. Parasuraman, “Neuroenhancement: enhancing brain and mind in health and in disease,” Neuroimage, 85 889 –894 (2014). NEIMEF 1053-8119 Google Scholar


B. A. Coffman, V. P. Clark and R. Parasuraman, “Battery powered thought: enhancement of attention, learning, and memory in healthy adults using transcranial direct current stimulation,” Neuroimage, 85 895 –908 (2014). NEIMEF 1053-8119 Google Scholar


F. Fregni and A. Pascual-Leone, “Technology Insight: noninvasive brain stimulation in neurology perspectives on the therapeutic potential of rTMS and tDCS,” Nat. Clin. Pract. Neurol., 3 (7), 383 –393 (2007). Google Scholar


M. S. George and G. Aston-Jones, “Noninvasive techniques for probing neurocircuitry and treating illness: vagus nerve stimulation (VNS), transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS),” Neuropsychopharmacology, 35 (1), 301 –316 (2010). NEROEW 0893-133X Google Scholar


A. R. Brunoni et al., “Clinical research with transcranial direct current stimulation (tDCS): challenges and future directions,” Brain Stimul., 5 (3), 175 –195 (2012). Google Scholar


J.-p. Lefaucheur et al., “Evidence-based guidelines on the therapeutic use of transcranial direct current stimulation (tDCS),” Clin. Neurophysiol., 128 (1), 56 –92 (2017). CNEUFU 1388-2457 Google Scholar


A. Antal et al., “Low intensity transcranial electric stimulation: safety, ethical, legal regulatory and application guidelines,” Clin. Neurophysiol., 128 (9), 1774 –1809 (2017). CNEUFU 1388-2457 Google Scholar


H. R. Siebner et al., “Consensus paper: combining transcranial stimulation with neuroimaging,” Brain Stimul., 2 (2), 58 –80 (2009). Google Scholar


T. O. Bergmann et al., “Combining non-invasive transcranial brain stimulation with neuroimaging and electrophysiology: current approaches and future perspectives,” Neuroimage, 140 4 –19 (2016). NEIMEF 1053-8119 Google Scholar


B. Khan et al., “Enhancing motor performance improvement by personalizing non-invasive cortical stimulation with concurrent functional near-infrared spectroscopy and multi-modal motor measurements,” Proc. SPIE, 9305 93051A (2015). PSISDG 0277-786X Google Scholar


M. A. Hunter et al., “Tracking the neuroplastic changes associated with transcranial direct current stimulation: a push for multimodal imaging,” Front. Hum. Neurosci., 7 495 (2013). Google Scholar


M. J. Rosa, J. Daunizeau and K. J. Friston, “EEG-fMRI integration: a critical review of biophysical modeling and data analysis approaches,” J. Integr. Neurosci., 09 (04), 453 –476 (2010). Google Scholar


B. Horwitz and D. Poeppel, “How can EEG/MEG and fMRI/PET data be combined?,” Hum. Brain Mapp., 17 (1), 1 –3 (2002). HBRME7 1065-9471 Google Scholar


J. J. Riera and A. Sumiyoshi, “Brain oscillations: ideal scenery to understand the neurovascular coupling,” Curr. Opin. Neurol., 23 (4), 374 –381 (2010). Google Scholar


R. McKendrick, R. Parasuraman and H. Ayaz, “Wearable functional near infrared spectroscopy (fNIRS) and transcranial direct current stimulation (tDCS): expanding vistas for neurocognitive augmentation,” Front. Syst. Neurosci., 9 1 –14 (2015). Google Scholar


C.-H. Han et al., “Hemodynamic responses in rat brain during transcranial direct current stimulation: a functional near-infrared spectroscopy study,” Biomed. Opt. Express, 5 (6), 1812 (2014). BOEICL 2156-7085 Google Scholar


A. C. Merzagora et al., “Prefrontal hemodynamic changes produced by anodal direct current stimulation,” Neuroimage, 49 (3), 2304 –2310 (2010). NEIMEF 1053-8119 Google Scholar


M. Muthalib et al., “Transcranial direct current stimulation induced modulation of cortical haemodynamics: a comparison between time-domain and continuous-wave functional near-infrared spectroscopy,” Brain Stimul., 8 (2), 392 –393 (2015). Google Scholar


B. Khan et al., “Functional near-infrared spectroscopy maps cortical plasticity underlying altered motor performance induced by transcranial direct current stimulation,” J. Biomed. Opt., 18 (11), 116003 (2013). JBOPFO 1083-3668 Google Scholar


K. Ishikuro et al., “Cerebral functional imaging using near-infrared spectroscopy during repeated performances of motor rehabilitation tasks tested on healthy subjects,” Front. Hum. Neurosci., 8 292 (2014). Google Scholar


K. T. Jones, F. Gözenman and M. E. Berryhill, “The strategy and motivational influences on the beneficial effect of neurostimulation: a tDCS and fNIRS study,” Neuroimage, 105 238 –247 (2015). NEIMEF 1053-8119 Google Scholar


U. Jindal et al., “Development of point of care testing device for neurovascular coupling from simultaneous recording of EEG and NIRS during anodal transcranial direct current stimulation,” IEEE J. Transl. Eng. Health Med., 3 1 –12 (2015). Google Scholar


E. Lareau et al., “Multichannel wearable system dedicated for simultaneous electroencephalography near-infrared spectroscopy real-time data acquisitions,” J. Biomed. Opt., 16 (9), 096014 (2011). JBOPFO 1083-3668 Google Scholar


J. Safaie et al., “Toward a fully integrated wireless wearable EEG-NIRS bimodal acquisition system,” J. Neural Eng., 10 (5), 056001 (2013). 1741-2560 Google Scholar


Y. Hoshi, “Functional near-infrared optical imaging: Utility and limitations in human brain mapping,” Psychophysiology, 40 (4), 511 –520 (2003). PSPHAF 0048-5772 Google Scholar


A. Torricelli et al., “Time domain functional NIRS imaging for human brain mapping,” Neuroimage, 85 28 –50 (2014). NEIMEF 1053-8119 Google Scholar


T. Durduran et al., “Diffuse optics for tissue monitoring and tomography,” Rep. Progr. Phys., 73 076701 (2010). Google Scholar


T. Z. Kincses et al., “Facilitation of probabilistic classification learning by transcranial direct current stimulation of the prefrontal cortex in the human,” Neuropsychologia, 42 (1), 113 –117 (2004). NUPSA6 0028-3932 Google Scholar


R. Ferrucci, F. Mameli and I. Guidi, “Transcranial direct current stimulation improves recognition memory in Alzheimer disease,” Neurology, 71 493 –498 (2008). NEURAI 0028-3878 Google Scholar


C. A. Dockery et al., “Enhancement of planning ability by transcranial direct current stimulation,” J. Neurosci., 29 (22), 7271 –7277 (2009). JNRSDS 0270-6474 Google Scholar


S. Elmer et al., “Direct current induced short-term modulation of the left dorsolateral prefrontal cortex while learning auditory presented nouns,” Behav. Brain Funct., 5 (1), 29 (2009). 1744-9081 Google Scholar


M. Wirth et al., “Effects of transcranial direct current stimulation (tDCS) on behaviour and electrophysiology of language production,” Neuropsychologia, 49 (14), 3989 –3998 (2011). NUPSA6 0028-3932 Google Scholar


V. Fiori et al., “Transcranial direct current stimulation improves word retrieval in healthy and nonfluent aphasic subjects,” J. Cognit. Neurosci., 23 (9), 2309 –2323 (2011). JCONEO 0898-929X Google Scholar


A. H. Javadi and V. Walsh, “Transcranial direct current stimulation (tDCS) of the left dorsolateral prefrontal cortex modulates declarative memory,” Brain Stimul., 5 (3), 231 –241 (2012). Google Scholar


J. T. Nelson et al., “Enhancing vigilance in operators with prefrontal cortex transcranial direct current stimulation (tDCS),” Neuroimage, 85 909 –917 (2014). NEIMEF 1053-8119 Google Scholar


F. Fregni et al., “Anodal transcranial direct current stimulation of prefrontal cortex enhances working memory,” Exp. Brain Res., 166 (1), 23 –30 (2005). EXBRAP 0014-4819 Google Scholar


S. H. Ohn et al., “Time-dependent effect of transcranial direct current stimulation on the enhancement of working memory,” Neuroreport, 19 (1), 43 –47 (2008). NERPEZ 0959-4965 Google Scholar


T. Zaehle et al., “Transcranial direct current stimulation of the prefrontal cortex modulates working memory performance: combined behavioural and electrophysiological evidence,” BMC Neurosci., 12 (1), 2 (2011). 1471-2202 Google Scholar


S. C. Andrews et al., “Improving working memory: the effect of combining cognitive activity and anodal transcranial direct current stimulation to the left dorsolateral prefrontal cortex,” Brain Stimul., 4 (2), 84 –89 (2011). Google Scholar


P. S. Boggio et al., “Go-no-go task performance improvement after anodal transcranial DC stimulation of the left dorsolateral prefrontal cortex in major depression,” J. Affective Disorders, 101 (1–3), 91 –98 (2007). JADID7 Google Scholar


M. A. Nitsche et al., “Treatment of depression with transcranial direct current stimulation (tDCS): a review,” Exp. Neurol., 219 (1), 14 –19 (2009). EXNEAC 0014-4886 Google Scholar


V. Jurcak, D. Tsuzuki and I. Dan, “10/20, 10/10, and 10/5 systems revisited: their validity as relative head-surface-based positioning systems,” Neuroimage, 34 (4), 1600 –1611 (2007). NEIMEF 1053-8119 Google Scholar


D. T. Delpy et al., “Estimation of optical pathlength through tissue from direct time of flight measurement,” Phys. Med. Biol., 33 (12), 1433 –1442 (1988). PHMBA7 0031-9155 Google Scholar


M. Patterson, B. Chance and B. Wilson, “Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties,” Appl. Opt., 28 (12), 2331 –2336 (1989). APOPAI 0003-6935 Google Scholar


A. Liebert et al., “Time-resolved multidistance near-infrared spectroscopy of the adult head: intracerebral and extracerebral absorption changes from moments of distribution of times of flight of photons,” Appl. Opt., 43 (15), 3037 (2004). APOPAI 0003-6935 Google Scholar


M. Febrero, P. Galeano and W. González-Manteiga, “Outlier detection in functional data by depth measures, with application to identify abnormal NOx levels,” Environmetrics, 19 (4), 331 –345 (2008). ENVCEE 1099-095X Google Scholar


A. Arribas-Gil and J. Romo, “Shape outlier detection and visualization for functional data: the outliergram,” Biostatistics, 15 (4), 603 –619 (2014). Google Scholar


R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria (2008). Google Scholar


M. Febrero-Bande and M. Oviedo de la Fuente, “Statistical computing in functional data analysis: the R package fda.usc,” J. Stat. Software, 51 (4), 1 –28 (2012). Google Scholar


D. Bates et al., “Fitting linear mixed-effects models using lme4,” J. Stat. Software, 67 (1), 1 –48 (2015). Google Scholar


R. V. Lenth, “Least-squares means: the R package lsmeans,” J. Stat. Software, 69 (1), 1 –33 (2016). Google Scholar


A. Cuevas, M. Febrero and R. Fraiman, “On the use of the bootstrap for estimating functions with functional data,” Comput. Stat. Data Anal., 51 (2), 1063 –1074 (2006). CSDADW 0167-9473 Google Scholar


T. Durduran and A. G. Yodh, “Diffuse correlation spectroscopy for non-invasive, micro-vascular cerebral blood flow measurement,” Neuroimage, 85 51 –63 (2014). NEIMEF 1053-8119 Google Scholar


N. Lang et al., “Effects of transcranial direct current stimulation over the human motor cortex on corticospinal and transcallosal excitability,” Exp. Brain Res., 156 (4), 439 –443 (2004). EXBRAP 0014-4819 Google Scholar


X. Zheng, D. C. Alsop and G. Schlaug, “Effects of transcranial direct current stimulation (tDCS) on human regional cerebral blood flow,” Neuroimage, 58 (1), 26 –33 (2011). NEIMEF 1053-8119 Google Scholar


N. Lang et al., “How does transcranial DC stimulation of the primary motor cortex alter regional neuronal activity in the human brain?,” Eur. J. Neurosci., 22 (2), 495 –504 (2005). EJONEI 0953-816X Google Scholar


W. B. Baker et al., “Pressure modulation algorithm to separate cerebral hemodynamic signals from extracerebral artifacts,” Neurophotonics, 2 (3), 035004 (2015). Google Scholar


M. S. Patterson et al., “Absorption spectroscopy in tissue-simulating materials: a theoretical and experimental study of photon paths,” Appl. Opt., 34 (1), 22 (1995). APOPAI 0003-6935 Google Scholar


Z. Turi et al., “When size matters: large electrodes induce greater stimulation-related cutaneous discomfort than smaller electrodes at equivalent current density,” Brain Stimul., 7 (3), 460 –467 (2014). Google Scholar


F. Ezquerro et al., “The influence of skin redness on blinding in transcranial direct current stimulation studies: a crossover trial,” Neuromodulation, 20 (3), 248 –255 (2017). Google Scholar


N. Khadka et al., “Minimal heating at the skin surface during transcranial direct current stimulation,” Neuromodulation, (2017). Google Scholar


F. Guarienti et al., “Reducing transcranial direct current stimulation-induced erythema with skin pretreatment: considerations for sham-controlled clinical trials,” Neuromodulation, 18 (4), 261 –265 (2015). Google Scholar


L. Zucchelli et al., “Method for the discrimination of superficial and deep absorption variations by time domain fNIRS,” Biomed. Opt. Express, 4 (12), 2893 –910 (2013). BOEICL 2156-7085 Google Scholar


R. Re et al., “Effect of a thin superficial layer on the estimate of hemodynamic changes in a two-layer medium by time domain NIRS,” Biomed. Opt. Express, 7 (2), 264 (2016). BOEICL 2156-7085 Google Scholar


J. Sutin et al., “Time-domain diffuse correlation spectroscopy,” Optica, 3 (9), 1006 (2016). Google Scholar


M. Pagliazzi et al., “Time domain diffuse correlation spectroscopy with a high coherence pulsed source: in vivo and phantom results,” Biomed. Opt. Express, 8 (11), 5311 (2017). BOEICL 2156-7085 Google Scholar


J. Selb et al., “Sensitivity of near-infrared spectroscopy and diffuse correlation spectroscopy to brain hemodynamics: simulations and experimental findings during hypercapnia,” Neurophotonics, 1 (1), 015005 (2014). Google Scholar


L. Sokoloff, “The physiological and biochemical bases of functional brain imaging,” Cognit. Neurodynamics, 2 (1), 1 –5 (2008). Google Scholar


N. Roche-Labarbe et al., “Somatosensory evoked changes in cerebral oxygen consumption measured non-invasively in premature neonates,” Neuroimage, 85 279 –286 (2014). NEIMEF 1053-8119 Google Scholar


E. Rounis et al., “Frequency specific changes in regional cerebral blood flow and motor system connectivity following rTMS to the primary motor cortex,” Neuroimage, 26 (1), 164 –176 (2005). NEIMEF 1053-8119 Google Scholar


N. K. Logothetis, “The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal,” Philos. Trans. R. Soc. B, 357 (1424), 1003 –1037 (2002). PTRBAE 0962-8436 Google Scholar


N. Accornero et al., “EEG mean frequency changes in healthy subjects during prefrontal transcranial direct current stimulation,” J. Neurophysiol., 112 (6), 1367 –1375 (2014). JONEA4 0022-3077 Google Scholar


J. C. Horvath, J. D. Forte and O. Carter, “Evidence that transcranial direct current stimulation (tDCS) generates little-to-no reliable neurophysiologic effect beyond MEP amplitude modulation in healthy human subjects: a systematic review,” Neuropsychologia, 66 213 –236 (2015). NUPSA6 0028-3932 Google Scholar


E. Basar et al., “Are cognitive processes manifested in event-related gamma, alpha, theta and delta oscillations in the EEG?,” Neurosci. Lett., 259 (3), 165 –168 (1999). NELED5 0304-3940 Google Scholar


W. Klimesch, “EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis,” Brain Res. Rev., 29 (2–3), 169 –195 (1999). BRERD2 0165-0173 Google Scholar


T. Harmony et al., “EEG delta activity: an indicator of attention to internal processing during performance of mental tasks,” Int. J. Psychophysiol., 24 (1–2), 161 –171 (1996). IJPSEE 0167-8760 Google Scholar


P. Sauseng et al., “Control mechanisms in working memory: a possible function of EEG theta oscillations,” Neurosci. Biobehav. Rev., 34 (7), 1015 –1022 (2010). NBREDE 0149-7634 Google Scholar


D. Keeser et al., “Prefrontal direct current stimulation modulates resting EEG and event-related potentials in healthy subjects: a standardized low resolution tomography (sLORETA) study,” Neuroimage, 55 (2), 644 –657 (2011). NEIMEF 1053-8119 Google Scholar


L. Jacobson et al., “Modulating oscillatory brain activity correlates of behavioral inhibition using transcranial direct current stimulation,” Clin. Neurophysiol., 123 (5), 979 –984 (2012). CNEUFU 1388-2457 Google Scholar


A. K. Engel and P. Fries, “Beta-band oscillations-signalling the status quo?,” Curr. Opin. Neurobiol., 20 (2), 156 –165 (2010). COPUEN 0959-4388 Google Scholar


T. E. Gladwin et al., “Enhancement of selective attention by tDCS: interaction with interference in a Sternberg task,” Neurosci. Lett., 512 (1), 33 –37 (2012). NELED5 0304-3940 Google Scholar


I. Laakso et al., “Electric fields of motor and frontal tDCS in a standard brain space: a computer simulation study,” Neuroimage, 137 140 –151 (2016). NEIMEF 1053-8119 Google Scholar


R. Sparing et al., “Bidirectional alterations of interhemispheric parietal balance by non-invasive cortical stimulation,” Brain, 132 (11), 3011 –3020 (2009). BRAIAK 0006-8950 Google Scholar


D. B. Stone and C. D. Tesche, “Transcranial direct current stimulation modulates shifts in global/local attention,” Neuroreport, 20 1115 –1119 (2009). NERPEZ 0959-4965 Google Scholar


R. Grech et al., “Review on solving the inverse problem in EEG source analysis,” J. Neuroeng. Rehabil., 5 (1), 25 (2008). Google Scholar


J. M. Kilner et al., “Hemodynamic correlates of EEG: a heuristic,” Neuroimage, 28 (1), 280 –286 (2005). NEIMEF 1053-8119 Google Scholar


C. Magri et al., “The amplitude and timing of the BOLD signal reflects the relationship between local field potential power at different frequencies,” J. Neurosci., 32 (4), 1395 –1407 (2012). JNRSDS 0270-6474 Google Scholar


P. Ritter and A. Villringer, “Simultaneous EEG-fMRI,” Neurosci. Biobehav. Rev., 30 (6), 823 –838 (2006). NBREDE 0149-7634 Google Scholar


Martina Giovannella is a doctoral student at ICFO—The Institute of Photonic Sciences (Barcelona) in the medical optics group. She received her master’s degree in physics at the University of Pisa (Pisa, Italy). Her research is focused on developing novel diffuse optical instrumentation, algorithms and protocols for brain monitoring for preclinical and clinical studies.

David Ibañez graduated with his Msc in biomedical engineering from the University of Boras (Sweden) in 2009. In 2011, he joined Starlab working as scientific researcher and project manager. His research interests include novel brain–computer interfaces for human–machine interaction, noninvasive brain-stimulation, neuromodulation, biomarkers and machine learning techniques applied to electrophysiological signal processing.

Clara Gregori-Pla is a doctoral student in the medical optics group at ICFO—The Institute of Photonic Sciences (Barcelona). She received her master’s degree in engineering physics in biomedical physics at the Royal Institute of Technology (Stockholm, Sweden). Her research is focused on developing novel diffuse optical instrumentation, algorithms and protocols and applying them for clinical studies.

Michal Kacprzak received his MSc degree in electrical engineering from Warsaw University of Technology in 2002 and his PhD degree from the Institute of Biocybernetics and Biomedical Engineering in 2011, where, from 2012 to 2015, he was a head of the Laboratory of Biomedical Optics. In 2015, he moved as postdoctoral researcher to the Institute of Photonic Sciences in Barcelona. His main area of interest is application of diffuse optics for monitoring and diagnosis of brain tissue.

Guillem Mitjà has a degree in industrial engineering with electronics specialization and obtained a MBA at EAE. After some years of professional experience in the industrial automation, he worked in Neuroelectrics as operations manager with responsibilities on the design and regulatory of electronic medical devices and biosignal sensors and stimulators.

Giulio Ruffini graduated (math/physics) from UC Berkeley and obtained a PhD in physics from UC Davis (1995). In 2000, he cofounded Starlab to transform research into technologies with positive impact. During the FET HIVE project, his team developed multielectrode hybrid EEG-transcranial stimulation devices, and he recently led the first demonstration of noninvasive brain-to-brain communication. In 2011, he cofounded Neuroelectrics to deliver clinical EEG-tCS systems. He collaborates with teams worldwide developing applications of EEG and noninvasive stimulation, as in the Luminous project studying consciousness.

Turgut Durduran was trained at the University of Pennsylvania. In 2009, he moved to ICFO—The Institute of Photonic Sciences, Spain, where he leads the medical optics group. His research interests revolve around the use of diffuse light to noninvasively probe tissue function. The group develops new technologies and algorithms and routinely translates them to preclinical, clinical, and industrial applications.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Martina Giovannella, David Ibañez, Clara Gregori-Pla, Michal Kacprzak, Guillem Mitjà, Giulio Ruffini, and Turgut Durduran "Concurrent measurement of cerebral hemodynamics and electroencephalography during transcranial direct current stimulation," Neurophotonics 5(1), 015001 (25 January 2018).
Received: 27 September 2017; Accepted: 20 December 2017; Published: 25 January 2018

Cited by 11 scholarly publications and 1 patent.




Statistical analysis


Blood circulation

Back to Top