Dataset of concurrent EEG, ECG, and behavior with multiple doses of transcranial electrical stimulation

We present a dataset combining human-participant high-density electroencephalography (EEG) with physiological and continuous behavioral metrics during transcranial electrical stimulation (tES). Data include within participant application of nine High-Definition tES (HD-tES) types, targeting three cortical regions (frontal, motor, parietal) with three stimulation waveforms (DC, 5 Hz, 30 Hz); more than 783 total stimulation trials over 62 sessions with EEG, physiological (ECG, EOG), and continuous behavioral vigilance/alertness metrics. Experiment 1 and 2 consisted of participants performing a continuous vigilance/alertness task over three 70-minute and two 70.5-minute sessions, respectively. Demographic data were collected, as well as self-reported wellness questionnaires before and after each session. Participants received all 9 stimulation types in Experiment 1, with each session including three stimulation types, with 4 trials per type. Participants received two stimulation types in Experiment 2, with 20 trials of a given stimulation type per session. Within-participant reliability was tested by repeating select sessions. This unique dataset supports a range of hypothesis testing including interactions of tDCS/tACS location and frequency, brain-state, physiology, fatigue, and cognitive performance.


Participants.
A total of twenty neurologically typical individuals (7 females, 13 males) between the ages of 19 − 43 (median age: 30; mean age: 29.10 ± 6.75) were recruited from the New York metropolitan area. Before enrolment, participants underwent a tES eligibility screening procedure where they were asked questions regarding their medical history, psychiatric history, and prior drug use (non/illicit) including history of seizures, depression, surgeries, pain, ear trauma, heart conditions, skin allergies, implantable medical devices (non/metallic), alcohol dependency, brain lesions, loss of consciousness, sleep disorders, and fainting at the sight of blood. Participants who met the eligibility criteria were then enrolled (see Fig. 1a,b). After passing the eligibility screening, informed consent was obtained for each participant. Experimental procedures were reviewed and approved by the Western Institutional Review Board and all procedures were conducted in accordance with the ethical guidelines set forth by the Declaration of Helsinki in 1964 and its later amendments. All participants were financially compensated for their participation.
Prior to each experimental session participants were asked to maintain their normal sleep-wake cycle, refrain from consuming alcohol or overly caffeinated foods/ beverages at least 4 h before their session, and refrain from eating any heavy meals at least 3 h before their visit. In order to facilitate expeditious EEG set-up, participants were also asked to refrain from wearing any heavy make-up or facial moisturizer, and refrain from wearing any hair products or scalp treatments 76 .
One participant was excluded from the experiments (participant 17) due to inability to follow instructions in performing the experiment's behavioral task (see Table 1). Four participants (participant number : 12, 15, 21, 22) were invited back to repeat both their experimental sessions for Experiment 2 (see Experimental Overview). Upon returning to repeat Experiment 2, returning participants were assigned a new participant number based on the incremental sequence of participant number assignment (i.e. 1 st participant enrolled assigned to participant number 01, 2 nd participant enrolled assigned to participant number 02, etc.). Participant 12 returned once and was assigned an updated participant number: 19; participant 15 returned once and was assigned an updated participant number: 18; participant 21 returned twice and was assigned updated participant numbers: 25 and 26; and participant 22 returned twice and was assigned updated participant numbers: 23 and 24 (see Table 2).
Questionnaires. Upon enrollment participants completed a demographic questionnaire to assess their age, gender, height, weight, years of education, handedness, English proficiency, typical exercise time each week, history of electrical stimulation, and sleep quality (quantified with the Pittsburgh Sleep Quality Index; PSQI 77 ). These questionnaires were only administered once upon participant's enrollment on their initial session (hence empty cells in GX_Subject Info & Behavioral Data.xlsx for participants who repeated Experiment 2; see Demographics, PSQI, and Pre Post Questionnaires).
At the beginning of each session participants completed a prequestionnaire consisting of questions regarding their activities in the past 24 h (see GX_Demo_PreQuest_Scales.pdf, page 2), in addition, they also completed assessment scales immediately before and after each session (Fig. 1c-h). These scales included assessments of sleepiness (Karolinska Sleepiness Scale; a 9-point scale ranging from 1-Extremely alert to 9-Very sleepy, great effort to keep awake, fighting sleep), discomfort scale (ranging from 1-No discomfort to 9-Extreme discomfort), pain (numeric pain scale; ranging from 0-No pain to 10-Worst possible pain), mood (ranging 1-Saddest to 9-Happiest), anxiety (ranging 1-Most relaxed to 9-Most tense), and energy level (ranging 1-Most tired to 9-Most energetic). Wellness scales were modeled after the standardized English version of the Karolinka Sleepiness Scale (KSS) [78][79][80] , which range from 1 to 9; and the pain scale utilized was modeled after standard Numeric Rating Scales (NRS) or Numeric Pain Scales (NPS) that conventionally range from 0 to 10 81,82 . ear plugs to aid in sound attenuation. Participants were seated ~57 cm in front of a 17-inch LCD monitor (Model: E173FPb; Dell Corp., Texas, USA; set at a 60 Hz refresh rate) and were asked to continuously perform a Compensatory Tracking Task (CTT) over the course of each session. On-screen task instructions presented to participants stated: Your goal in this task is to keep the ball near the target using your mouse or pointing device. The ball moves on its own, but is impacted by how you move the mouse. If possible, keep the ball in the center of the ring, or as close as possible to the inner circle as you can. The goal of the 2-dimensional task was to keep a moving circle constantly constrained near an annulus at the center of the screen (Fig. 2a). The circle was endowed with several kinematic properties including oscillatory and dampening forces 47 . The CTT was sourced and administered through the open source PEBL software version 2.1, an experimental psychology task program built upon C++ 83 . Participants provided task input through their dominant hand via a trackball (Logitech TrackMan Marble Mouse; Model: 904360-0403; Logitech International S.A., Lausanne, Switzerland) and were left uninterrupted throughout the duration of the task (even in cases where participants were in visible hypnagogic states; consistent with prior CTT designs). Before the start of each session participants performed a short practice session that lasted between 1-3 mins. The duration of the task was set to 70 mins for both Experiment 1 and 2 (extended for an additional 0.5 mins). The CTT ran continuously and uninterrupted over the duration of both experiments and participants were blinded to each experiment's block design and stimulation type. EEG data acquisition. EEG data were acquired using a wired Waveguard cap containing 32 Ag/AgCl recording channels (ANT Neuro, Hengelo, The Netherlands) and 29 interleaved plastic HD-holders (Soterix Medical Inc., New York, USA). The HD-holders were used to house the stimulation materials including conductive gel and stimulation electrodes (see HD-tES). For both stimulation and recording, conductive gel (SignaGel, Parker Laboratories Inc., New Jersey, USA) was placed between electrodes and the scalp using blunt tip (15-gauge; Cortech Solutions Inc., North Carolina, USA) syringes.
EEG recording electrodes were located at standard locations following the 10/10 international placement system (Fig. 2b). Signals were sampled at 2 kHz, referenced relative to CPz, online grounded at AFz, and amplified with an eego sport amplifier (ANT Neuro, Hengelo, The Netherlands). The amplifier bandwidth was set between 0-520 Hz and the acquisition voltages range was set to 1 V peak-to-peak. Prior to each recording, scalp impedances were monitored to ensure impedance levels were below 20 kΩ. To time-lock the CTT with the concurrent www.nature.com/scientificdata www.nature.com/scientificdata/ EEG a trigger was sent to the EEG amplifier at the start of the CTT (see Data Records). Following data acquisition, data were exported to a .cnt file.
Physiological monitoring. Cardiac activity was acquired with a lead I electrocardiogram (ECG) configuration and ocular motor activity with horizontal electrooculogram (EOG) configuration. Both were acquired with bipolar snap electrodes, concurrently with EEG and HD-tES. For lead I ECG, bipolar snap electrodes were placed on participants' chest approximately 5 cm below the center of their left (anode) and right (cathode) clavicle bone. A ground electrode was placed on participants' left hip, in close proximity to participants' iliac crest. For horizontal EOG, electrodes were placed at the outer canthus of the participants' left (anode) and right (cathode) eye. Prior to all peripheral electrode placement, the application sites on the skin were gently swabbed with alcohol pads (BD Alcohol Swabs, Becton Dickinson, New Jersey, USA) to increase electrode adhesion and contact quality.  www.nature.com/scientificdata www.nature.com/scientificdata/ A total of 9 Ag/AgCl sintered ring stimulation electrodes (Soterix Medical Inc., New York, USA) were placed at standard EEG 10/10 locations in manner forming three possible Nx1 HD-tES 84 ; where N = 3, 4, 4 for frontal, motor, and parietal, respectively. For each montage, N electrodes were selected for the outer (surround) ring electrode, and one electrode was selected as the center electrode. Note, some electrode locations were shared across montages with varied "ring" or "center" assignments. In this way, a single 9 electrode position HD-tES set-up was prepared for all experiments. For frontal stimulation, outer or surround electrodes (N = 3) were placed at AF3, FT7, FC3 and the return or center electrode was placed at F5 (Fig. 2b,c). For motor stimulation outer or surround electrodes (N = 4) were placed at FT7, FC3, CP3, TP7, and the return or center electrode was placed at C5 (Fig. 2b,c). For parietal stimulation outer or surround electrodes (N = 4) were placed at C5, C1, P1, TP7 and the return or center electrode was placed at CP3 (Fig. 2b,c). To visualize stimulation electrode placement in realistic 3D space, an MRI-derived head model was used together with the ROAST toolbox 27,85 in MATLAB (2018b and 2019b; MathWorks, Massachusetts, USA). Each montage shared stimulation sites and no two stimulation montages were applied at the same time.

HD-tES.
In terms of stimulation waveforms, a monophasic DC (0 Hz) or biphasic sinusoidal waveform (5 or 30 Hz) was applied. For DC, for each cephalic location (frontal, motor, parietal) the respective center electrode was used as the cathode with the surround electrodes as anodes. For biphasic sinusoidal stimulation all electrodes switched between being an anode and cathode at a rate set by the frequency of stimulation (i.e. 5 Hz or 30 Hz).
Using conventional procedures 86 , stimulation electrodes were placed in plastic HD-holders (Soterix Medical Inc., New York, USA), which were embedded in the EEG cap (interleaved with EEG electrode positions). To  www.nature.com/scientificdata www.nature.com/scientificdata/ prepare for electrode placement, participants' hair was parted with a blunt Q-tip or blunt tip syringe at each stimulation location through the annulus of the plastic HD-holders, then the HD-holders were slowly backfilled with ~15 mL of HD gel (Soterix Medical Inc., New York, USA). The stimulation electrodes were then placed in the HD-holders and were fully submerged in electrolyte gel, then the plastic holder was covered with the associated plastic covering to secure the electrode in place.
Stimulation was delivered using the MxN 9-channel high-definition transcranial electrical stimulation stimulator (Soterix Medical Inc., New York, USA). The multichannel stimulation device supplied a current controlled current source where the current amplitude of each individual stimulation channel could be adjusted. Stimulation channel impedances were measured relative to F5, prior to the commencement of stimulation. During each session, the stimulator was triggered using a pulse generated in MATLAB and relayed to a National Instruments (NI) -DAQmx device (National Instruments Inc., Texas, USA). A corresponding trigger was sent to the EEG amplifier to time-lock stimulation and EEG (Fig. 2d). For details on stimulation settings for each stimulation dose see Table 3.  www.nature.com/scientificdata www.nature.com/scientificdata/ Prior to the start of each session participants underwent a pre stimulation test to determine their tolerance and comfort with each stimulation montage. Stimulation was delivered for approximately 20 secs (5-sec ramp-up/ down) at an initial net current of 1 mA (peak-to-peak) and 0.5 mA (peak-to-peak) for Experiment 1 and 2, respectively. Participants were then asked to rate the level of pain/discomfort on a numeric pain scale (0-No pain/ discomfort to 10-Worst pain you've ever felt). For Experiment 1, if participants rated a 5 or above, the stimulation intensity was halved (0.5 mA minimum amplitude); whereas for Experiment 2 if participants rated below 5, the stimulation intensity was increased to 1 mA. In both experiments the adjusted stimulation intensity was 0.5 mA at minimum and was retested to ensure that it was well tolerated. All participants in Experiment 2 received 1 mA of stimulation (see Table 2) whereas stimulation intensities for Experiment 1 were either 0.5 or 1 mA, based on pre stimulation testing (see Table 1). For Experiment 1, each session (of the 3 sessions per participant) consisted of three stimulation enabled blocks (10 mins each) interleaved with four stimulation off blocks (10 mins each). For each session, three of the nine total stimulation conditions were pseudo-randomly assigned to stimulation enabled blocks -to ensure participants received each of the nine stimulation conditions across the nine total stimulation enabled blocks in Experiment 1. During each stimulation enabled block the assigned stimulation condition was administered as 4 consecutive stimulation trials (see Fig. 2e,f). Each stimulation trial consisted of a 5-sec ramp-up period, 30-sec stimulation period, and a 5-sec ramp down period.

Experimental overview. In
For Experiment 2, each of the two sessions (one session for F30 and one session for M30) consisted of one stimulation off blocks (20 mins), which was followed by five consecutive stimulation enabled blocks (10 mins each; Fig. 2g,h). For each session, only one of the two stimulation conditions (F30 or M30) were pseudo-randomly assigned to be administered during the five stimulation enabled blocks. The assigned stimulation condition was administered as 4 consecutive trials during each stimulation enabled blocks. Thus, a total of 20 trials of stimulation (30 secs each plus 5-sec ramp up and down) were delivered for each session in Experiment 2.
A precise experimental sequence with EEG (channel Cz), behavioral timeseries, and associated triggers are illustrated for Experiment 1 and Experiment 2 (Fig. 3) for exemplary participants 10 and 24, respectively. This demonstrates participant-wise implementation of the block design in Fig. 2e-h from the data collection perspective. In all experiments, participants continuously performed a task, with associated performance, EEG, and physiology recorded.
At their first session of Experiment 1 (Fig. 3b), participant 10 was pseudo-randomly assigned to receive stimulation conditions M0, P0, and M5 (indicated by insets). Each stimulation condition was repeated over 4   www.nature.com/scientificdata www.nature.com/scientificdata/ consecutive trials in a block (Stim Enabled); for each trial current was ramped up (over 5 secs), sustained for 30 secs at a maximum designated current intensity (0.5 or 1 mA; Max Stim Current), then ramped back down. Session 2 of Experiment 1, for participant 10, continued similarly with F0, P30, and F30 stimulation applied during the stimulation enabled blocks (Fig. 3c). Session 3 of Experiment 1, for participant 10, also continued similarly with F5, M30, and P5 stimulation applied during the stimulation enabled blocks (Fig. 3d).
In session 1 of Experiment 2 (Fig. 3f), participant 24 received the M30 stimulation condition (indicated by insets). The M30 stimulation condition was repeated over 20 consecutive trials divided across five consecutive stimulation enabled blocks (Stim Enabled); for each trial current was ramped up (over 5 secs), sustained for 30 secs at a maximum designated current intensity (0.5 or 1 mA; Max Stim Current), then ramped back down. Session 2 of Experiment 2, for participant 24, continued similarly with the F30 stimulation condition applied during the stimulation enabled blocks (Fig. 3g).
The band-pass filtering applied in Fig. 3 for display purposes, produces large voltage artifacts in the EEG during ramp-up and ramp-down, as well as removing DC voltage artifacts, which is present in unfiltered EEG-stimulation data 21 .

Data Records
All raw data described in the text 87 can be accessed directly at: https://doi.org/10.5281/zenodo.3837212. The main data repository contains raw data for Experiment 1 and 2, in.cnt format (see EEG and Behavioral Data) 87 ; in addition, down sampled data for Experiment 1 and 2 are indexed in linked repositories 88,89 in.mat format, which is compatible with MATLAB and Python. The dataset is also provided according to the Brain Imaging Data Structure (BIDS) 90  All experimental sessions and applied stimulation conditions are summarized in Table 1 and Table 2 for Within the Demographic page of the GX_Subject Info & Behavioral Data.xlsx file the columns indicated each participants' (rows) participant number (Sub#), date of data collection (month/date/year), age (in years), gender (Male or Female), height (in cm), weight (in kg), ethnicity, race, years of education, handedness (left or right handed), whether English was their first language, their proficiency in English, whether they exercise or not, how many hours they exercise, if they've ever had electrical stimulation before, and how recently have they had electrical stimulation.
Within the PSQI page of the GX_Subject Info & Behavioral Data.xlsx file the data were arranged by participant (rows) and PSQI questions (columns). Each column other than participant number (Subj#) and date, are labeled according to their correspondence to each PSQI question. For example, PSQI-5-a corresponds to PSQI question number 5, part a, asking: During the past month, how often have you had trouble sleeping because you cannot get to sleep within 30 minutes? Notes were made when participants indicated responses that were outside the scope of the questionnaire or if decisions were made to average or quantize participants' responses, when appropriate. The full list of PSQI questions as well as scoring the PSQI can be found here 94 .
Within the All Behavioral page of the GX_Subject Info & Behavioral Data.xlsx file data were arranged according to participant number (Sub#), session label, session number, file number (File Num), a count of all records, date of data collection, presumed start and end times of each session, stimulation amplitude used for each session and each montage (one number indicates that a single stimulation amplitude was applied for all stimulation montages within a session), session type/montage order (Arms), responses to all questions in the Pre-Questionnaire www.nature.com/scientificdata www.nature.com/scientificdata/ participants' assigned number and that participant's experimental session number. For example, for participant 6, session 2 the corresponding folder will be 0602.
Within each session folder is the raw EEG, ECG, and EOG data. These data were exported to.cnt format and accompanying.evt file, without any applied filters or data augmentation. Each.cnt file is labeled with a project pseudonym (GX) followed by the participant number, year of data acquisition, month of acquisition, date of acquisition and end time of acquisition in 24-hour clock notation. For example, GX_07_2019-11-15_20-06-07. cnt indicates that this data file was collected from participant 7, on November (11) 15 th , 2019 and the recording ended at 20:06:07 hrs (08:06:07 PM). The session folder also contains a.mat file and a.txt file. These files contain session timers as well as the montage/s entered for each session. Within the MATLABfilestream < participant number > .mat files relevant variables include Montages, and stim. For Experiment 1 the Montages variable contains the three stimulation conditions that were applied for the defined session. For example, for participant 8 the Montages file for session 2 (0802) indicates 'P30' , 'F0' , 'P0' . This indicates that the P30 stimulation was applied in stimulation block 1, whereas F0 was applied during stimulation block 2, and P0 was applied during stimulation block 3 (see Fig. 3a Stim Enabled time periods). For Experiment 2 the Montages variable contains the single stimulation condition that was applied for the whole session, the stimulation condition was repeated three times within the variable for code execution dependencies. For example, for participant 19 the Montages file for session 2 (1902) indicates 'M30' , 'M30' , 'M30' . This indicates that the M30 stimulation condition was applied for all stimulation trials (20 trials total) during session 2 (see Fig. 3e Stim Enabled time periods). The stim variable within the MATLABfilestream < participant number > .mat files on the other hand kept track of parameters that were passed to PEBL and variables defined within MATLAB for file execution. The amp variable defined the stimulation amplitude that was intended to be applied during each session (see Table 1, 2 for stimulation amplitudes applied). The blocklenmins and the blocklen variables defined the length of each block as outlined in Fig. 2e,g in units of mins and secs, respectively. The variables timeon, rampup, rampdown, and repeat defined the length of time stimulation was applied for, the ramp up and ramp down times, and the amount of times stimulation was repeated within a block (Fig. 2f,h), respectively. Other variables within the stim variable and within the MATLABfilestream < participant number > .mat files contained parameters that tracked parallel experimental timelines. The text files, MATLABfilestream < participant number > .txt were used as additional time parameter trackers and were not pertinent to any study outcome. Within each session folder is also a folder with identical session naming. This subfolder contains the exported CTT data as well as a CTT summary in.csv and.txt formats, respectively. The.cnt and.evt can be used in both MATLAB and python using free and open access accompanying libraries 95 .
The CTT data are located within each session's folder in a subfolder with an identical name. For example, for participant 19 session 2, their CTT data can be found by navigating through 1902 > 1902 > ptracker-1902. csv. The ptracker-< participant number + session number > .csv file contains metrics collected from the trackball. Some of these useful metrics include the participant number (subnum), the time stamp (time), the cursor's x-position (posX), the cursor's y-position (posY), the change in time steps in ms (timeDelta), and the cursor's radial deviation from the center of the screen (deviation). For additional details on CTT output see PEBL's ptracker documentation.
The physiological (EEG, ECG, EOG) and the behavioral (CTT) data can be time aligned by using triggers sent to the EEG amplifier. The start of the CTT is indicated by the first Block Start trigger (code:02), which can be found in the EEG data once it has been imported for post processing. The start of each stimulation trial is indicated with a stimulation start trigger (code:16). After the stimulation start trigger was delivered, to the stimulator and EEG amplifier, the stimulation current was ramped up over the course of 5 secs until it reached the desired intensity. After 30 secs of stimulation at the desired intensity a stimulation stop trigger (code:32) was delivered to the stimulator and EEG amplifier, and the stimulation intensity was ramped down over the course of 5 secs. This process repeated 4 times for each stimulation enabled block.

technical Validation
Our technical validation consisted of consulting computational current flow modeling, quality control of behavioral metrics, and quality control of EEG and ECG voltage changes.
A high-resolution finite element model (FEM) was generated based on an MR image (1 mm 3 voxel size) in order to predict current flow pattern in three different montages used in the study (i.e. frontal, motor, parietal). . The total stimulating current in each montage (i.e. frontal, motor, parietal) was applied at the center electrode and the four (three for frontal) surrounding electrodes were assigned as ground. The FEM model predicted the expected magnitude of current (electric fields) reaching underlying brain tissue and that electrode montages encapsulated the 3 brain regions of interest (frontal, motor, parietal; Fig. 4a-f). Mean EEG recordings of scalp voltages during stimulation, were in accord with the magnitude and spatial distribution of FEM model predictions of scalp voltage for each montage (Fig. 4d-f). To examine the designated frequency of stimulation applied during each stimulation trial, the Welch power spectral density (PSD; Fig. 4g); during stimulation topoplots; and time-frequency spectrograms (Fig. 4k-p) for EEG data were computed for pre, during, and post stimulation, accordingly, for all participants. Together these metrics aided in corroborating the applied stimulation frequency (with the PSD and spectrograms) as well as the stimulation's spatial location (with the topoplots). The stimulation voltage artifacts in the trial-wise timeseries data (Fig. 4k,n) also aided in stimulation montage and frequency corroboration. These participant and trial-wise data are fully indexed for EEG PSDs 91      www.nature.com/scientificdata www.nature.com/scientificdata/ In addition, we extracted and tabulated the PSD at the most dominant frequency during stimulation for all participants and all trials, across both experiments (see Online-only Table 1 for Experiment 1 and Online-only Table 2 for Experiment 2). All dominant frequencies matched the frequency of programmed stimulation (i.e. 0 Hz for F0, M0, P0; 5 Hz for F5, M5, P5; and 30 Hz for F30, M30, P30). These are tabulated together with behavioral data (CTT deviation change) for each trial of all participants across both experiments. These tabulated data also contain trial-wise means and standard deviations, whereas group-wise sample sizes, means, standard deviations, and confidence intervals can be found in Table 4. Metrics for participants who repeated the experiment were averaged across repeats before being added to group-wise metrics.
CTT data were examined after each experimental session to ensure data were a continuous stream that approximated the time of the EEG, ECG and EOG recordings. Since CTT data were sampled at the frequency of the experimental monitor's screen refresh rate, data were sometimes nonuniformly sampled at 60 Hz. This nonuniformity was eliminated post hoc with uniform sampling using the resample function in MATLAB. The CTT data were then examined 30 secs before, during, and after stimulation. The mean across trials for each time period was computed for each stimulation montage across participants for both Experiment 1 and 2 (Fig. 4h,i). For Experiment 2, we examined the average percent change in deviation (taken over 20 trials) for each participant and each stimulation type (F30 and M30; Fig. 4i). We then expanded our examination to look at each participants' performance (percent change in deviation) on a trial-to-trial level for each stimulation type (F30 and M30; Fig. 4j). Four of the participants who repeated the experiments were denoted with unique participant numbers including participant number 12,19; participant number 15, 18; participant number 21, 25, 26; and participant number 22, 23, 24 ( Fig. 4i,j). Two trials for session 1 and 2 for participant 24 are examined in detail from different perspectives including as a timeseries signal plotted with concurrent CTT (Fig. 4k,n); as a time-frequency breakdown (spectrogram) of the EEG voltage from C3 (Fig. 4l,o), and as behavioral comparison of CTT data pre, during, and post stimulation (Fig. 4m,p). Behavioral CTT data for all three repeats for participant number 22, 23, 24 were examined in detail for F30 and M30 stimulation montages (Fig. 4q,r).
During EEG, ECG, and EOG data collection; data were monitored in real time to ensure continuous data streaming. For participant 01 session 01 and session 02, data collection was halted at the end of the sessions (0101 and 0102) due to a technical error. Session 0101 contained 1 trial of F0 for ~2 mins, whereas session 0102 contained 1 trial of F0 for 30 sec. For participant 02 session 01 (data record 0201), the first trial of F30 contained a ramp-up/down time of 30 secs rather than 5 secs. Subsequent trials for 0201 contained a 5 sec ramp-up/down time. EEG amplitudes over time were examined to confirm expected stimulation-generated scalp voltage (voltage artifact see 21,22 ) during stimulation and concurrency with ECG and CTT data (Fig. 4k,n). EEG spectra during stimulation were examined to ensure that they contained significant power at the frequencies of stimulation (0 Hz, 5 Hz, 30 Hz; Fig. 4l,o). The CTT deviations, for the aforementioned examinations, were computed and its distributions were examined on a trial-by-trial basis (Fig. 4q,r).

Code availability
The latest version of all accompanying code for this dataset can be acquired within this repository: https://github. com/ngebodh/GX_tES_EEG_Physio_Behavior. MATLAB, version 2018b and 2019b were utilized with functions from EEGlab 96 , Raincloud plots toolbox 97 , and ANT neuro's import functions 95 .   Table 4. Experiment 1 and 2 group-wise CTT deviation change and EEG PSD. The EEG PSD (dB/Hz) at the most dominant frequency during stimulation at electrode FC5 for frontal stimulation, C3 for motor stimulation, and CP5 for parietal stimulation; as well as the CTT deviation change (%) is tabulated across participants and all trials for each stimulation condition. The sample size, mean, standard deviation (Std), and 95% confidence interval (95% CI (+/−)) for each metric is also tabulated.