
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)13458-5
10.1016/j.heliyon.2024.e37427
e37427
Research Article
Physiology-inspired bifocal fronto-parietal tACS for working memory enhancement
Pupíková Monika a
Maceira-Elvira Pablo bc1
Harquel Sylvain b1
Šimko Patrik a
Popa Traian bc
Gajdoš Martin def
Lamoš Martin def
Nencha Umberto bg
Mitterová Kristína ad
Šimo Adam f
Hummel Friedhelm C. friedhelm.hummel@epfl.ch
bcg2⁎
Rektorová Irena adf2
a Applied Neuroscience Research Group, Central European Institute of Technology, Masaryk University, Brno, Czech Republic
b Neuro-X Institute (INX), École Polytechnique Fédérale de Lausanne (EPFL), Chemin des Mines 9, 1202, CH, Geneva, Switzerland
c Neuro-X Institute (INX), EPFL Valais, Clinique Romande de Réadaptation Sion, Switzerland
d International Clinical Research Center, Faculty of Medicine and St. Anne's University Hospital, Masaryk University, Brno, Czech Republic
e Brain and Mind Research, Central European Institute of Technology, Masaryk University, Brno, Czech Republic
f First Department of Neurology, Faculty of Medicine and St. Anne's University Hospital, Masaryk University, Brno, Czech Republic
g Clinical Neuroscience, University of Geneva Medical School, Geneva, Switzerland
⁎ Corresponding author. Clinical Neuroengineering, Neuro-X Institute (INX), École Polytechnique Fédérale de Lausanne (EPFL), Chemin des Mines 9, 1202, CH, Geneva, Switzerland. friedhelm.hummel@epfl.ch
1 Authors contributed equally.

2 Authors contibuted equally

06 9 2024
30 9 2024
06 9 2024
10 18 e374276 5 2024
14 8 2024
3 9 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Aging populations face significant cognitive challenges, particularly in working memory (WM). Transcranial alternating current stimulation (tACS) offer promising avenues for cognitive enhancement, especially when inspired by brain physiology. This study (NCT04986787) explores the effect of multifocal tACS on WM performance in healthy older adults, focusing on fronto-parietal network modulation. Individualized physiology-inspired tACS applied to the fronto-parietal network was investigated in two blinded cross-over experiments. The first experiment involved monofocal/bifocal theta-tACS to the fronto-parietal network, while in the second experiment cross-frequency theta-gamma interactions between these regions were explored. Participants have done online WM tasks under the stimulation conditions. Network connectivity was assessed via rs-fMRI and multichannel electroencephalography. Prefrontal monofocal theta tACS modestly improved WM accuracy over sham (d = 0.30). Fronto-parietal stimulation enhanced WM task processing speed, with the strongest effects for bifocal in-phase theta tACS (d = 0.41). Cross-frequency stimulations modestly boosted processing speed with or without impairing task accuracy depending on the stimulation protocol. This research adds to the understanding of physiology-inspired brain stimulation for cognitive enhancement in older subjects.

Highlights

• Individualized frontal monofocal theta tACS enhanced high load working memory (WM) task accuracy in healthy older subjects.

• Fronto-parietal cross-frequency theta-gamma tACS enhanced WM processing speed at the expense of task accuracy in most cases.

• Fronto-parietal theta in-phase tACS enhanced WM processing speed without any detrimental impact on task accuracy.

• Magnitude and symmetry of induced electric fields in bi-focal theta tACS was related to tACS-induced speed enhancement.

Keywords

tACS
Multifocal
Orchestrated brain stimulation
Healthy aging
Cognition
Systems neuroscience
Electric field modelling
Neuroimaging
Working memory
Abbreviations

EEG electroencephalography

EF electric fields

fMRI functional magnetic resonance imaging

FPCN fronto-parietal control network

NIBS non-invasive brain stimulation

PAC phase-amplitude coupling

rs-fc resting state functional connectivity

rIPL right inferior parietal lobule

rMFG right middle frontal gyrus

rs-fMRI resting state fMRI

RT reaction times

tACS transcranial alternating current stimulation

WM working memory
==== Body
pmc1 Introduction

Global demographics indicate a rising proportion of older adults. Despite prolonged life expectancy, the duration of healthy cognitive years has remained unchanged [1]. Working memory (WM) represents a vital system, maintaining and manipulating transient information necessary for diverse tasks, which deteriorates with aging [2]. The age-compromised efficacy of WM can be attributed to various neurobiological factors, with notable changes in the functioning of fronto-parietal networks. Animal and human research has consistently emphasized integrative interactions between the prefrontal and parietal cortices for WM [[3], [4], [5], [6]].

The coordinated interplay within and between brain regions [7] orchestrated by brain oscillations and their intra- and cross-frequency interactions have been shown to be central to WM processes [8,9]. The fascinating multi-tiered organization of these oscillations has been described, revealing how they interconnect functional brain systems across diverse spatiotemporal scales, providing the brain with a sophisticated solution for handling cognitive tasks [10,11]. Theta oscillations (4–8 Hz) in particular coordinate neural assemblies, ensuring the temporal consistency of items retained in WM [12]. A specific form of cross-frequency coupling – theta (4–8 Hz) gamma (>30 Hz) phase-amplitude coupling (PAC) – where the amplitude of faster gamma rhythms aligns with the phase of the slower theta rhythms, was shown to temporally organize information in WM. Nested gamma cycles may code for individual items held in WM [12,13].

Based on the current understanding of key nodes of the fronto-parietal cognitive network, non-invasive brain stimulation (NIBS) interventions have been employed to modulate their activity and enhance our understanding of their roles in WM processes [[14], [15], [16], [17]]. Most studies targeted single nodes in the dorsolateral prefrontal cortex, inferior frontal gyrus, inferior parietal lobule, and superior temporal gyrus using either repetitive transcranial magnetic stimulation or transcranial direct current stimulation [[18], [19], [20], [21], [22], [23]]. These studies did not demonstrate any conclusive results in terms of which brain region should be targeted for enhancing distinct WM functions in specific age/disease groups. Based on the understanding of the functional role of oscillatory regional and interregional interactions, a different concept of interventional strategies based on transcranial alternating-current stimulation (tACS) has been developed and evaluated [3,24,25]. tACS has shown potential in selectively modulating large-scale cortical networks in a causal, frequency and phase specific manner with enduring behavioral effects [3,17,[24], [25], [26]]. Growing evidence from tACS studies implies the feasibility of modulating WM processing and capacity by targeting cortical regions with theta-band tACS [3,27]. Notably, stimulation of the left prefrontal cortex with a complex theta-gamma waveform during a WM task led to a noticeable modulation of WM performance in healthy young people [3]. This opens an interesting interventional window to enhance WM in healthy older people.

Our current study assesses how the principles of intra-areal and inter-areal synchronization and cross-frequency coupling impact WM performance in healthy older subjects following two tACS experiments. In the first experiment, monofocal theta tACS and synchronized theta tACS between two distant regions were tested in the same subjects. The latter stimulation protocol is motivated by the idea of increasing interregional synchronization and connectivity within a network [24,26,28,29], with the hypothesis that a reduced phase lag (in-phase stimulation) between sites would promote an optimal inter-areal coupling and thus, optimal communication and efficiency [30,31]. The second theoretical framework examined here relates to the concept of cross-frequency theta-gamma PAC during WM. In an effort to synchronize cross-regional neural activities, we employed cross-frequency theta-gamma entrainment of interregional oscillatory interactions within the fronto-parietal network via bifocal tACS.

We established spatially and temporally individualized protocols for both experiments. Our primary outcome was WM task performance (normalized accuracies and response times) as measured by an n-back task. As little is known about factors causing interindividual differences in stimulation responses, we explored EEG and electric field properties and their relations to stimulation outcome. Given the known association of PAC with WM performance [25], we examined whether the baseline PAC is associated with tACS responsiveness. Changes in the offline resting-state functional connectivity (rs-FC) within the stimulated network were also studied using fMRI as a potential readout of stimulation effects.

2 Methods

2.1 Study design and procedure

This is a multicentric study consisting of two experiments on two cohorts investigating the effects of neuromodulation of the fronto-parietal network by means of mono- or bifocal theta-theta/theta-gamma tACS. An individualized approach in terms of spatial targets and stimulation frequency was adopted to maximize the personalization of the stimulation [24,25,[32], [33], [34]]. Each cohort underwent a randomized cross-over study with different tACS stimulation protocols, including sham stimulations, to determine the effects on WM-task performance (see Fig. 1). The primary outcome was the performance (normalized task accuracy and response times) during an n-back task conducted online during the stimulation (Fig. 2 and Supplementary materials for task description). The trial was registered (2021-04-19) in ClinicalTrial.gov under NCT04986787. See Supplementary materials for details.Fig. 1 A. Design of the study. Prior to the stimulation sessions, participants underwent baseline multichannel scalp EEG during the task and task-fMRI measurements to personalize the stimulation parameters in terms of individual frequency and target location. Left: For the MUNI cohort, the stimulation protocols aimed to entrain oscillations within the theta range and synchronize fronto-parietal regions. Participants received monofocal stimulation in frontal or parietal areas and bifocal frontal and parietal stimulation in and out of phase. In the MUNI cohort only, resting-state fMRI (rs-fMRI) scans were acquired before and immediately after each stimulation session. The tACS and n-back tasks were conducted in the NIBS laboratory adjacent to the MRI scanner, with subjects easily moved between the two labs with a maximum of 5 min between the examinations. An alternating current was delivered using pairs of concentric rubber electrodes (Outer electrode: Ø 75 mm with a hole Ø 40 mm, inner electrode Ø 20 mm) Right: For the EPFL cohort, the stimulation protocols aimed to affect cross-frequency coupling between theta and gamma oscillations. Individualized gamma and theta cross-frequency tACS stimulation over the prefrontal and parietal cortex was applied. Continuous theta stimulation was applied while gamma stimulation was delivered either continuously or in bursts aligned to the theta peak. There was no before-and-after fMRI. B. Stimulation montage for bifocal stimulation. For monofocal stimulation, the electrode placement was the same, but only one current source was active. Image created with the use of SimNIBS software. Head is derived from a standard template provided. C. Stimulation protocols for inter-areal synchronization in experiment 1 (MUNI) D. Stimulation protocols for inter-areal cross-frequency stimulation in experiment 2 (EPFL).

Fig. 1

Fig. 2 A. Schematic diagram of a study timeline. B. Main behavioral outcome – n-back task modified from Ref. [35]; task performed online during each stimulation condition. Participants completed two levels of difficulty during the stimulation. Participants indicate whether the current stimulus matches the one from n steps (2-back/3-back) earlier in the sequence. The task was performed also during fMRI to spatially individualize targets based on each participants activation.

Fig. 2

2.2 Participants

In total, two cohorts of healthy participants were recruited, one in each research center (MUNI cohort: n = 20, EPFL cohort: n = 20), with intact cognition and no serious neuropsychiatric conditions (see Table 1 and Supplementary material). None of the participants had contraindications for either MRI or tACS. Each subject signed the informed consent form in accordance with the ethics codes and relevant regulations approved by the ethics committee of Masaryk University (EKV-2020-019) and from Swiss Ethics (SNCTP000003928 | BASEC2020-00127).Table 1 Mean demographic, neuropsychological screening and EEG data., MoCA = Montreal Cognitive Assessment, SD = Standard deviation.

Table 1		Cohort	
MUNI	EPFL	
n = 20	n = 20	
Age	Mean	64.49	67.6	
SD	(6.98)	(4.70)	
Men/Women Ratio		11/9	8/12	
MoCA	Mean	27.21	26.35	
(SD)	(2.27)	(2.30)	
Individual Theta (Hz)	Mean	4.51	4.64	
(Min-Max)	(4.17–4.98)	(4.32–5.28)	
Individual Gamma (Hz)	Mean	–	34.04	
(Min-Max)	–	(30.06–37.04)	

2.3 Transcranial alternating current stimulation

tACS was performed through two battery-driven stimulators attached to concentric electrodes placed over the task-fMRI individualized (see Supplementary materials) areas of the right middle frontal gyrus (rMFG) and right inferior parietal lobule (rIPL) for 20 min. Each concentric electrode pair was connected to a separate stimulator (current source, see Fig. 1B). A pair of concentric electrodes consists of smaller inner electrode, which is a target electrode (inner electrode Ø 20 mm) and a larger reference electrode around the target electrode (outer electrode: Ø 75 mm with a hole Ø 40 mm). The current flows between the inner and outer electrode. We used T1 MRI scan-based frameless stereotactic neuro-navigation to specify the exact location of the inner electrode center in each individual. The waveform of the stimulation was sinusoidal without DC offset. Individualized peak frequencies within the theta (4–7 Hz, both cohorts) and gamma (>30 Hz, EPFL cohort) range were calculated for each individual based on the pre-stimulation task EEG (see Table 1). The sham stimulation was applied with the same settings, but the stimulator was turned off after 30 s. For details see Fig. 1 and Supplementary materials. The impedance was kept below 20 kΩ.

2.3.1 Inter-areal synchronization tACS – MUNI cohort

A continuous, sinusoidal stimulation within a theta range was applied at an intensity of 3 mA peak-to-peak. Protocols included monofocal theta stimulation over the frontal site only and over the parietal site only, and bifocal theta stimulation of both sites with a relative phase shift of 0° (in-phase) and 180°(out-of-phase). For monofocal stimulation protocols, electrodes were placed over both stimulation sites, despite only one of stimulators being active.

2.3.2 Cross-frequency inter-areal tACS – EPFL cohort

This stimulation consisted of two components: continuous slow-wave theta frequency applied over the frontal or parietal cortices with concurrent fast gamma over the other pair of electrodes applied both continuously or in bursts. Gamma-burst stimulation was synchronized with the peaks within the continuous theta frequency and lasted 1000 ms. This results in four different stimulation protocols: slow theta over the frontal/parietal areas with continuous gamma (θ Pre - γ Par; γ Pre - θ Par resp., 2 mA peak-to-peak) or gamma applied in bursts (θ Pre - γ Burst, γ Burst - θ Par resp., 4 mA peak-to-peak) over the other area.

2.4 MRI data acquisition and pre-processing

The MRI data was acquired via 3.0 T Magnetom Siemens Prisma at both centers. For the structural MRI data, the T1 MPRAGE sequence was used. A gradient-echo, T2 echo-planar imaging sequence was used for the rs-fMRI and task-fMRI. The fMRI data were preprocessed and analyzed with SPM12 in MATLAB R2019a. For the data preprocessing pipeline and control for spatial abnormalities, see Supplementary materials.

2.5 Pre-stimulation task EEG data acquisition and baseline theta-gamma phase-amplitude coupling estimation

Participants underwent EEG recording concurrently with the pre-stimulation n-back task, in a shielded room using a 256-channel (MUNI) or 64 active channel (EPFL) system with task stimuli projected on a 150 cm screen. Task recording lasted approximately 20 min. The EEG sampling rate was 1 kHz with Cz (MUNI) or FCz (EPFL) as a reference electrode. EEG preprocessing involved epoching (between −7 and −10 to + 1 s relative to stimulus presentation for 2-back and 3-back tasks respectively), data filtering, artifact removal via ICA, and average re-referencing. PAC was computed in the −2.5 s–0 s time window using the Mean Vector Length method [36] on phase frequencies between 4 and 7 Hz (1 Hz step) and amplitude frequencies between 18 and 45 Hz (1 Hz step). Maximum electrode-level theta-gamma PAC was computed to assess WM load variations, contrasted between 2-back and 3-back tasks (contrast referred as Δ PAC scores) and analyzed for their relationship to performance. This approach aimed to better capture the PAC related to WM due to the variation in cognitive load across the tasks. The contrast was extracted from the frontal and parietal stimulation targets (frontal and parietal ROI). See Supplementary materials for more details.

2.6 Pre- and post-stimulation resting-state fMRI data analysis in MUNI cohort

After preprocessing, the resting-state functional connectivity (rs-FC) was assessed between the fronto-parietal control network (FPCN) major seeds using nine coordinates as described previously [37], see Supplementary material for details. We extracted fMRI signals from the region of interest (spheres, radius 6 mm) around seed positions and used means as representative signals. To analyze connectivity, we used the Pearson correlation between representative signals, converted to z values using the Fisher r-to-z transformation.

2.7 Montage current simulations

We simulated the electric fields (EF) induced at each stimulation site in all participants a posteriori, in order to estimate the influence of their magnitude on performance. Using SimNIBS [38], we constructed head models for each participant based on T1 and T2 MRI images. We extracted the magnitude of the EF induced within a 10 mm radius of each target. We calculated additional measures for the bifocal stimulation conditions, namely symmetry (i.e., relative magnitude of the field at one site with respect to the second site) and joint magnitude (i.e., the vectoral combination of the magnitude at both sites).

2.8 Statistical analysis

Data analyses were performed using R software (lme4 package). For the behavioral data analysis, our primary outcomes were normalized accuracy and normalized response times, which refer to the relative change in accuracy/response times as compared to the first block of each session. By normalizing the data to the first block of each session, we accounted for variability between sessions, as we did not measure pre-stimulation task performance on each stimulation sessions (see Fig. 2A). For response time, we considered the response times of the correct trials only. The behavioral data were checked for outliers (more than 3 SD from the group mean). We conducted a statistical analysis using linear mixed models (LMM) to investigate the interaction effect between the stimulation protocol and cognitive load on measures of normalized accuracy and response time, as well as accuracy and response time. The initial linear mixed model, referred to as the first-level model, included the stimulation protocol and cognitive load as fixed effects, along with their interaction. We then expanded the analysis to second-level models, which incorporated additional fixed factors such as training blocks, and random intercepts/factors of subject and days (for detailed model choice process, please refer to the reports on statistics). To ensure the appropriateness of our models, we assessed their fit by comparing the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). These criteria take into account model overfitting by penalizing the number of parameters included in the model. In the MUNI cohort, the final LMM for normalized accuracy included also random intercept for subjects and days, and the LMM for normalized speed included random intercept for subjects. In the EPFL cohort, the final LMM for normalized accuracy included fixed effect for blocks with no additional random intercept/factors, and the LMM for normalized speed included a random intercept for days and subjects. For the resting-state fMRI analysis, LMM for each network pair were calculated. The differences of the z values of the connectivity changes (FC post–FC pre) were used as a dependent variable. For the final model, the variable stimulation protocol was used as an independent variable imputed as fixed effect, with no random intercepts or factors, as they did not improve the model fit. Post-hoc analyses were performed for the significant connectivity seed pairs. We further correlated the differences of the z values of the connectivity changes (FC post–FC pre) with within-session learning curves of accuracy and response time separately for each stimulation protocol. For montage current simulations, the LMM were computed to assess the interaction between stimulation protocol, cognitive load, and magnitude at each site (2 variables)/symmetry/joint magnitude which were imputed as fixed effects on measures of normalized accuracy and response time, as well as accuracy and response time (RT). For pre-stimulation Δ PAC scores, behavioral outcomes (normalized accuracy and RT) were taken as dependent variable with stimulation condition, Δ PAC Frontal and Δ PAC parietal were taken as fixed factors with no additional random intercepts/factors. Statistical significance (threshold p < 0.05) was assessed using the ANOVA function with Satterthwaite's approximations in the lmerTest package. Post-hoc pairwise comparisons were conducted by computing the estimated marginal means using the emmeans package. Effect sizes are reported as partial eta-squared (pη2) for F tests and Cohen's d for pairwise comparisons and were obtained using the effectsize package.

3 Results

3.1 Experiment 1. effect of monofocal and bifocal theta tACS on working memory (MUNI)

The participants in the MUNI cohort received either monofocal tACS to frontal/parietal areas or bifocal tACS to fronto-parietal regions (i.e., in-phase and out-of-phase). The accuracy was significantly different among the n-back task conditions (F[1,2567] = 240.37, η2 = 0.09, p < 0.0001). When examined separately, we did not find a significant difference among the interventions in either the 2-back (F[4,1243] = 1.16, η2 = 0.0037, p = 0.32) or the 3-back tasks (F[4,1268] = 0.62, η2 = 0.0019, p = 0.64; Fig. S1 Left). When comparing the accuracy normalized to the first block, we found a significant interaction between the stimulation and the tasks (F[4,2567] = 8.30, η2 = 0.01, p < 0.0001). Participants had a significantly larger accuracy increase in the 3-back task when receiving monofocal frontal stimulation compared to the sham condition (T[1273] = 3.48, d = 0.30, p = 0.004; see Fig. 3 left).Fig. 3 Normalized accuracy in the n-back tasks for the MUNI (experiment 1) and EPFL (experiment 2) cohorts. *sig. p < 0.05.

Fig. 3

We found response times to differ significantly among tasks (F[1,2571] = 89.83, η2 = 0.03, p < 0.0001). However, they did not vary significantly between the stimulation protocols in either the 2-back (F[4,1276] = 1.60, η2 = 0.005, p = 0.16) or the 3-back (F[4,1276] = 1.88, η2 = 0.005, p = 0.11) tasks (see Fig. S2 Left). In terms of the normalized speed, we found the change to be different among the tasks (F[1,2567] = 48.55, η2 = 0.02, p < 0.0001). Post-hoc comparisons revealed that participants experienced significantly different changes in speed in the 3-back task when receiving monofocal frontal (d = 0.27), monofocal parietal (d = 0.286), and in-phase bifocal (d = 0.41) stimulations as compared to the sham stimulation, without compromising task accuracy (see Fig. 4 left).Fig. 4 Normalized speed in the n-back tasks for the MUNI (experiment 1) and EPFL (experiment 2) cohorts. *sig. p < 0.05.

Fig. 4

3.2 Experiment 2. effects of cross-frequency tACS on working memory. (EPFL)

The participants in the EPFL cohort received cross-frequency tACS in individualized gamma and theta frequencies applied over the prefrontal and parietal cortices, with theta stimulation kept constant and gamma stimulation applied either continuously or in bursts. In regard to accuracy, a significant interaction between the tasks and the blocks of each difficulty was found (F[1,2454] = 13.18, η2 = 0.0053, p = 0.0002), suggesting performance changing at different rates. No significant difference was found between the stimulation protocols in either the 2-back (F [4,39] = 0.72, η2 = 0.06, p = 0.57) or the 3-back (F[4,83] = 0.32, η2 = 0.02, p = 0.86) tasks compared to sham (see Fig. S1 Right). For normalized accuracy, significant effects of the stimulation were observed in both the 2-back (F [4] = 14.29, η2 = 0.04, p < 0.0001) and the 3-back (F [4] = 12.29, η2 = 0.04, p < 0.0001) tasks. In the 2-back task, the participants experienced a significantly smaller improvement in accuracy when receiving parietal-theta-prefrontal-gamma-burst tACS (T[1238] = 5.52, d = 0.49, p < 0.0001); in the 3-back tasks, participants experienced a significantly smaller improvement in accuracy in all interventions except for prefrontal-theta-parietal-gamma-continuous, when compared to sham stimulation (see Fig. 3 right). In other words, most of the active stimulation protocols led to normalized accuracy impairment compared to sham.

There was a significant difference in response time among the tasks (F[1,2452] = 4.42, η2 = 0.001, p = 0.03). In the 3-back task, there was a significant difference among the stimulation types (F [4,40] = 4.29, η2 = 0.24, p = 0.004), with participants being significantly faster with the parietal-theta-prefrontal-gamma-burst stimulation compared to the sham (T [4] = 4.6, d = 0.74, p = 0.03; see Fig. S2 right). The normalized speed differed among the tasks (F [1] = 14.22, η2 = 0.005, p = 0.0001) with a significant effect of stimulation in both the 2-back (F [4] = 10.45, η2 = 0.03, p < 0.0001) and the 3-back (F [4] = 7.91, η2 = 0.02, p < 0.0001) tasks. In the 2-back, the change in speed was significantly larger in both the prefrontal-gamma-continuous and theta-parietal (T[1225] = 4.83, d = 0.42, p < 0.0001) and the prefrontal-theta and parietal-gamma-burst (T[1225] = 4.46, d = 0.40, p = 0.0001). In the 3-back task, participants experienced a significantly larger change in speed when receiving the prefrontal-theta and parietal-gamma-continuous (T[1237] = 2.96, d = 0.24, p = 0.02) and the parietal-theta and prefrontal-gamma-burst (T[1237] = 3.04, d = 0.27, p = 0.02) tACS stimulation (see Fig. 4 right).

3.3 Effect of monofocal and bifocal theta tACS upon rs-fMRI connectivity

Within the fronto-parietal network only, the connectivity change between the right anterior prefrontal cortex (raPFC) and right dorsolateral prefrontal cortex (rdlPFC) (F [4] = 2.66, p = 0.04) was affected by the stimulation condition and was significantly different from the sham stimulation condition (t [41] = 2.33, p = 0.022). This difference resulted from a decrease in raPFC-rdlPFC connectivity following stimulation compared to pre-stimulation levels (t [38] = 2.068, p = 0.046, see Fig. 5). Connectivity changes were not correlated with any behavioral changes caused by the stimulation. For full results see Supplementary material and the Report on Statistics.Fig. 5 A Visualization of FPCN network seeds used for the rs-fMRI data analysis. raPFC-rdlPFC connection is highlighted by a line. B Connectivity between the raPFC-rdlPFC. Legend: l/raPFC = left/right anterior prefrontal cortex, ACC – anterior cingulate cortex, l/raIPL = left/right anterior intraparietal lobule, l/rdlPFC = left/right dorsolateral prefrontal cortex, l/rINS = left/right insula.

Fig. 5

3.4 PAC scores and effect of stimulation

We found a relationship between the Δ PAC scores at the parietal target and the normalized accuracy under different stimulation conditions (F [3] = 2.78, η2 = 0.12, p = 0.04). Post-hoc comparisons revealed that the relationship of parietal Δ PAC scores and normalized accuracy differed between the frontal theta condition and bifocal theta stimulation conditions (frontal vs. in-phase p = 0.033; frontal vs. out-of-phase p = 0.046). Notably, the correlations between parietal Δ PAC scores and the relevant mean normalized accuracies were not significant. For all results, see Supplementary materials.

3.5 Influence of the magnitude of the induced electric fields (EFs) on WM enhancement: experiment 1

We found effects of individually induced EFs on normalized response times, but not on normalized accuracies in the bifocal theta-theta tACS conditions. We found a significant triple interaction between the stimulation and the magnitude of the EF induced in the prefrontal target and that of the parietal target (F [1] = 5.05, η2 = 0.07, p = 0.02). To explore the interaction, we introduced a feature representing the symmetry of EF at the prefrontal and parietal targets, defined by their ratio, and the vector magnitude of both fields to measure the total injected energy (see Table S2). These metrics were integrated into the linear mixed models, replacing the frontal/parietal magnitude variables. Our analysis confirmed significant differences in response time changes across stimulation conditions for symmetry of EFs (F [1] = 5.68, η2 = 0.1, p = 0.02), but not for their joint magnitude (F [3] = 1.17, p = 0.33).

The stimulation-related differences appear to be most pronounced when both EF magnitudes were similar at both stimulation sites, in both the in-phase and the out-of-phase conditions, with the opposite direction of changes (see Fig. 6). In the in-phase condition, higher EF magnitudes predicted more pronounced reaction time enhancement. Conversely, out-of-phase stimulation with higher EF magnitudes led to more pronounced speed decrements, especially in the 3-back task.Fig. 6 Visualization of triple interactions among the bifocal stimulation conditions and the magnitude of the EFs induced at both stimulation sites (i.e., prefrontal and parietal cortices) acting on the speed in the 2-back and the 3-back task. The dots represent individual datapoints; the planes reflect the statistical model adjusted to these points. In the in-phase condition, the model describes steeper improvements in speed predicted for individuals exposed to EFs of higher magnitude. In contrast, the model describes speed decreases in individuals receiving out-of-phase stimulation who were predicted to have speed decreases when exposed to EFs of higher magnitude. This detrimental effect is only apparent in the 2-back and not the 3-back task. Note that the vertical axis has the same scale in all four panels. Even though there is no numerical grading along the axes, all plots have the same scale.

Fig. 6

3.6 Influence of the magnitude of the induced EF on WM enhancement: experiment 2

We found effects of the individually induced EF on normalized accuracies, but not on normalized response times in the cross-frequency tACS protocols. We observed a significant triple interaction between the stimulation and the magnitude of electrical fields at both stimulation sites (F [3] = 3.04, η2 = 0.7, p = 0.03), which was particularly significant in the 3-back task (F [3] = 4.08, η2 = 0.17, p = 0.01). The accuracy change was modulated differently across conditions (see Fig. 7). In continuous conditions, with gamma tACS on the prefrontal cortex and theta-tACS on the parietal cortex, accuracy improvement increased with higher EF magnitudes, particularly when fields at both sites were similar. In contrast, when theta tACS was applied to the prefrontal cortex and gamma tACS to the parietal cortex, the greatest improvement occurred with lateralized fields, decreasing with symmetric and intense fields. Additional modeling to examine these interactions revealed no significant effects for either symmetry or joint magnitude of the EFs.Fig. 7 Visualization of triple interactions among the stimulation and the magnitude of the EFs induced at both stimulation sites (i.e., prefrontal and parietal cortices) acting on the accuracy in the 3-back task. The dots represent individual datapoints; the planes reflect the statistical model adjusted to these points. Specifically, when applying tACS within the gamma band to the prefrontal cortex and theta-tACS to the parietal cortex (stimulation condition 1), the rate of improvement in accuracy within a training session became steeper as a function of the magnitude of the induced EF, and this effect was more pronounced when the magnitude was similar at both stimulation sites. In contrast, when applying theta to the prefrontal cortex and gamma tACS to the parietal cortex (stimulation condition 2), the rate of improvement was steepest when the EFs were lateralized, and least pronounced when they were both symmetric and of high intensity. These interactions seem to be more attenuated in both conditions including gamma bursts (i.e., stimulation conditions 3 and 4). Note that even though there is no numerical grading along the axes, all the vertical axes in all plots have the same scale.

Fig. 7

4 Discussion

An effective exchange of information across neural structures necessitates the spatiotemporal coordination of oscillatory activity as outlined by oscillatory hierarchy [9,42]. Slower frequencies such as theta oscillations foster this inter-areal, cross-frequency neural communication. Specifically, increased fronto-parietal theta synchrony is associated with the maintenance and manipulation of information in WM [43,44]. Older populations were shown to be predisposed toward long-range theta desynchronization in comparison to younger adults [25,[39], [45], [46]]. We explored the online cognitive effects of individualized, monofocal, and multifocal physiology-based tACS in healthy older adults, focusing on enhancing intra-areal and inter-regional cognitive processing within the fronto-parietal network. To achieve these goals, we studied: (1) the entrainment of individualized theta-band neural oscillations at frontal and parietal sites separately as well as synchronously, and (2) inspired by the role of physiological theta-gamma coupling in WM [3,4], we examined the effects of cross-frequency theta and gamma tACS applied to the prefrontal and parietal cortices.

We found that the most pronounced effects were observed in the n-back task with higher cognitive loads. The monofocal prefrontal stimulation was the only stimulation condition that led to modest yet significant enhancements in WM task accuracy, surpassing other fronto-parietal stimulation protocols. On the other hand, fronto-parietal stimulation tended to modulate WM task processing speed with or without concomitant impairment of the task accuracy. Positive mild effects on response times without impacting the task accuracy were observed with the bifocal fronto-parietal theta in-phase tACS, followed by monofocal theta stimulation at parietal and frontal targets, and cross-frequency prefrontal theta and parietal gamma-continuous tACS.

4.1 Effects of intra-areal and inter-areal theta synchronization

Previous research has shown that tACS-induced phase synchronization of frontal and parietal regions improves WM performance [3,24,26,43], while desynchronization impairs it [47]. In our study, only the monofocal frontal stimulation provided modest accuracy modulation of the 3-back task as compared to sham stimulation. While no significant effect on accuracy was observed, the fronto-parietal in-phase condition significantly enhanced processing speed, particularly under higher WM loads (3-back), consistent with prior findings of faster reaction times using in-phase fronto-parietal stimulation [3,24,26,43,47,48]. Notably, this study extends previous works by examining monofocal conditions; both frontal and parietal monofocal stimulations led to slightly faster reaction times as compared to sham. There are some inconsistencies between prior monofocal stimulation studies and our findings. While the efficacy of frontal monofocal theta upon WM performance has been reported [49,50], a recent review identified the parietal node of the fronto-parietal network as superior for WM enhancement, both in reaction times and accuracies [51]. Discrepancies between our findings and the results of others may arise from using the n-back task, which strains the frontal lobe and its connections due to the higher emphasis of updating demands [52]. This contrasts with other studies that employed tasks focusing on pure WM capacity or simpler match-to-sample tasks (as listed in Nissim et al. [51]). Additionally, age might have been a major contributing factor. Unlike most prior studies involving young adults, we included older participants. Using resting state EEG, Perinelli et al. [53] reported decreased frontal intra-areal connectivity and lower modularity (less segregated network) in older versus younger adults; fronto-parietal and fronto-temporal inter-areal connectivity was rather increased [53,54]. Notably, connectivity differences between older and younger adults were mostly contributed by the theta (and alpha) band. Therefore, modulating frontal oscillation by monofocal theta tACS might bring significant (albeit modest) WM task improvement specifically in older subjects.

Interestingly, the behavioral effects of frontal monofocal tACS were accompanied by significant pre-post changes in resting state fMRI connectivity within the stimulated fronto-parietal network. Specifically, we observed a reduction in connectivity strength between the right anterior prefrontal cortex (aPFC) and the right dorsolateral prefrontal cortex, as opposed to unchanged connectivity in the sham condition. While Violante et al. [43] reported enhanced online connectivity within stimulated FPCN during theta in-phase tACS, Abellaneda-Perez et al. [55] showed, similarly to our study result, that prefrontal theta tACS led to post-versus pre-stimulation rs-connectivity decrease between discrete seeds of FPCN. While both the anterior prefrontal cortex and DLPFC are core components of the FPCN, they appear to have different roles in WM processing. The MFG is part of a subnetwork critical to WM performance, and aPFC seems to be part of a subnetwork more aligned with introspective processes [56], multitasking [40] or delayed response tasks [57]. The changes in offline rs-connectivity were not directly associated with the online behavioral changes in this study, pointing out that the offline rs-connectivity change is a coarse and indirect measure of the online tACS-induced network effects.

4.2 Effects of fronto-parietal theta-gamma cross-frequency tACS

Enhanced cross-frequency coupling, particularly gamma oscillations nested within theta cycles in prefrontal areas, have been linked to better WM performance [12]. This study expands upon this by exploring physiologically inspired multifocal stimulation and its effects on inter-areal interactions. Stimulation protocols involved the application of continuous slow-wave theta frequency over either the frontal or parietal cortices, coupled with concurrent fast gamma frequency stimulation over the other electrode pair, delivered either continuously or in bursts.

Unexpectedly, our stimulation conditions led to enhanced speed processing at the cost of task accuracy; this effect was more prominent in the WM task with a higher cognitive load. In fact, prefrontal-theta and parietal gamma-continuous tACS was the only configuration that led to a significant though modest increase in processing speed without a concomitant speed-accuracy trade-off. Our observation is in agreement with a recent meta-modeling study by Wischnewski et al. [58], which synthesizes findings from 28 placebo-controlled tACS studies. They showed dose-dependent improvement in working memory task performance after prefrontal theta and parietal gamma stimulation. Further, our finding that parietal-theta when paired with prefrontal-gamma stimulation negatively impacted WM performance aligns well with their reported detrimental effects on accuracy for prefrontal gamma modulation. Negative effects of prefrontal gamma stimulation were specific to accuracy; no detrimental effects were reported for response times. Speed-accuracy trade-off can manifest through the use of heuristics and mental shortcuts to ease the cognitive load associated with decision-making [59,60], potentially strengthened by our stimulation. Previous studies hypothesized engagement of the superior medial prefrontal cortices and the hyper-direct pathway connecting the prefrontal cortex to the subthalamic nucleus (STN) [[61], [62], [63]]; abnormal connectivity of the pathway has been implicated in impulse control disorders and might subserve the speed-accuracy trade-off [[64], [65], [66]]. Alternatively, the overall effectiveness of inter-areal cross-frequency coupling scenarios in the current study could have been compromised due to the misalignment between the externally induced fast gamma component and the inherent local theta rhythms leading to distorted theta-phase locking [67,68]; for visualization see Fig. 8), specifically in gamma-burst scenarios. Continuous gamma, delivered uniformly across theta phases, might be less affected allowing local theta to modulate gamma bursts through variations in local excitability [30]. Notably, inter-regional oscillatory mechanisms are dynamic, as recent studies show high-frequency gamma activity nested within various phases of the slow theta wave, depending on cognitive demand. This enables dynamic coupling or decoupling of the fronto-parietal control network to posterior areas, adjusted to cognitive effort [69].Fig. 8 Visualization of the proposed misalignment between the externally induced fast gamma component and the inherent local theta rhythms in different stimulation conditions. Consequently, the theta-gamma coupling might be disproportionately skewed, which in turn would distort specific theta phase-locking present between prefrontal and distant posterior areas.

Fig. 8

Future studies should examine the integration of naturally occurring phase-lags in inter-areal theta activity and consider systems that superimpose both fast and slow rhythms in a closed-loop fashion. Specifically closed-loop systems that measure brain activity and apply targeted stimulation based on these measurements might be well suited in respect to optimal phase lag.

4.3 Individualization of stimulation paradigms for WM enhancement

In our study, we adopted an individualized approach, using each participant's activation maps and determining their peak frequencies, not to compare individualized versus standard methods but to maximize stimulation benefits based on current knowledge [32]. The substantial variation in targets within the large MFG region underscored the need for individualization. On the other hand, our study's mean frequency of 4.58 Hz aligns closely with the 4.5 Hz average reported in another tACS study [24], suggesting limited frequency variability. This narrow range of frequency variability raises questions about the necessity of frequency individualization, at least in older healthy subjects. A standardized frequency might suffice, thereby simplifying the tACS protocol without compromising efficacy.

The electrical field modelling data offer another possibility for individualization. Prior research indicates that individual anatomical differences lead to variations in electric field magnitude [70], affecting stimulation outcomes. In the present study, various stimulation condition outcomes depended on the magnitude and particularly on the symmetry of electrical fields induced at the two stimulated sites. This was most evident in the theta in-phase condition, where more balanced fields led to better behavioral effects, possibly more efficiently enhancing the fronto-parietal theta coherence crucial for demanding WM processes [44]. Thus, our results emphasize the need for pre-stimulation optimization by current modelling [71,41].

4.4 Limitations

This study has several limitations that require careful consideration. Notably, while our objective was to affect behavior by modulating the oscillatory activity within the fronto-parietal network in a physiology-inspired way by tACS, we only assessed behavioral changes without monitoring and adjusting the tACS to the individual “background” ongoing oscillatory brain activity by e.g. online EEG. Notably, EEG recording with concurrent tACS is currently limited due to tACS artifacts and problems with their elimination [72,73]. Another limitation is the fact that only one session of stimulation was studied. Further studies are needed to explore the effects of multiple tACS sessions. The fields induced beneath the electrodes were rather low in our study and could have led to small effects sizes [74]. This might stem from our focus on focality which is often achieved at the expense of field intensity [75,76]. Notably, we found that the symmetry of the induced electrical field magnitudes influenced the task performance during various stimulation conditions.

4.5 Conclusion

We used individualized, monofocal and multifocal, physiology-inspired tACS to modulate the fronto-parietal network with the aim of enhancing WM performance in healthy older individuals. The main effects were found for the 3-back task (with a higher cognitive load). Only prefrontal monofocal theta stimulation produced a modest yet significant improvement in WM task accuracy as compared to the sham. Fronto-parietal bifocal stimulation conditions enhanced WM task processing speed; however, this enhancement was at the expense of task accuracy in most cases, eventually leading to WM impairment. Fronto-parietal theta in-phase stimulation produced mild to moderate speed acceleration without any detrimental effect on task accuracy. The speed enhancement of this stimulation protocol was related to the magnitude and particularly to the symmetry of the induced electric fields at both stimulation targets, thus stressing the importance of current modelling for stimulation protocol optimization. Our research contributes to developing tACS-based intervention strategies to enhance WM in healthy older adults; further work is warranted in neurological conditions with WM deficits.

Ethics statement

Each subject signed the informed consent form in accordance with the ethics codes and relevant regulations approved by the ethics committee of Masaryk University (EKV-2020-019) and from Swiss Ethics (SNCTP000003928 | BASEC2020-00127).

Data availability statement

The datasets generated and/or analyzed during the current study are not publicly available due accordance with informed consents signed by study participants but are available from the corresponding author on a request.

CRediT authorship contribution statement

Monika Pupíková: Writing – review & editing, Writing – original draft, Visualization, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. Pablo Maceira-Elvira: Writing – review & editing, Visualization, Formal analysis, Data curation. Sylvain Harquel: Writing – review & editing, Visualization, Formal analysis, Data curation. Patrik Šimko: Writing – review & editing, Project administration, Investigation, Data curation, Conceptualization. Traian Popa: Writing – review & editing, Validation, Project administration, Methodology, Investigation, Data curation, Conceptualization. Martin Gajdoš: Validation, Software, Formal analysis, Data curation. Martin Lamoš: Software, Formal analysis, Data curation. Umberto Nencha: Investigation, Data curation. Kristína Mitterová: Investigation, Formal analysis. Adam Šimo: Investigation. Friedhelm C. Hummel: Writing – review & editing, Supervision, Methodology, Funding acquisition, Conceptualization. Irena Rektorová: Writing – review & editing, Supervision, Methodology, Funding acquisition, Conceptualization.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A Supplementary data

The following are the Supplementary data to this article.Multimedia component 1

Multimedia component 1

Multimedia component 2

Multimedia component 2

Multimedia component 3

Multimedia component 3

Multimedia component 4

Multimedia component 4

Acknowledgements

This project has received funding from La-Grant of the 10.13039/100000001 Swiss National Science Foundation (SNSF) in cooperation with the 10.13039/501100001824 Czech Science Foundation (GAČR) GF21-13462L (Novel individualized brain stimulation network-based approaches to improve cognition in healthy seniors and patients with MCI), from the Czech Ministry of Health NU21-04-00652 (Implication of cognitive reserve in non-pharmacological intervention outcomes), and from the grant LX22NPO5107 (MEYS): Financed by 10.13039/501100000780 European Union – Next Generation EU. We acknowledge the contribution of the core facility MAFIL supported by MEYS CR (LM2023050 Czech-BioImaging), part of the Euro-BioImaging (www.eurobioimaging.eu) ALM and Medical Imaging Node (Brno, CZ). This project was further supported by the 10.13039/501100001711 Swiss National Science Foundation (SNSF , 320030L_197899 , NiBS-iCog) to F.H.; the 10.13039/501100023369 Defitech Foundation (Morges, CH) to F.H.; the 10.13039/100015473 Wyss Center for Bio and Neuroengineering to F.H.; Personalized Health and Related Technologies (PHRT 2017-205) of the ETH Domain to F.H. We thank Anne Johnson for the English editing.

Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.heliyon.2024.e37427.
==== Refs
References

1 GBD 2019 Demographics Collaborators Global age-sex-specific fertility, mortality, healthy life expectancy (HALE), and population estimates in 204 countries and territories, 1950-2019: a comprehensive demographic analysis for the Global Burden of Disease Study 2019 Lancet 396 2020 1160 1203 https://doi:10.1016/S0140-6736(20)30977-6 33069325
2 Park D.C. Reuter-Lorenz P. The adaptive brain: aging and neurocognitive scaffolding Annu. Rev. Psychol. 60 2009 173 196 10.1146/annurev.psych.59.103006.093656 19035823
3 Alekseichuk I. Turi Z. Amador de Lara G. Antal A. Paulus W. Spatial working memory in humans depends on theta and high gamma synchronization in the prefrontal cortex Curr. Biol. 26 2016 1513 1521 10.1016/J.CUB.2016.04.035 27238283
4 Axmacher N. Henseler M.M. Jensen O. Weinreich I. Elger C.E. Fell J. Cross-frequency coupling supports multi-item working memory in the human hippocampus Proc. Natl. Acad. Sci. U.S.A. 107 2010 3228 3233 10.1073/PNAS.0911531107/SUPPL_FILE/PNAS.200911531SI.PDF 20133762
5 Cowan N. What are the differences between long-term, short-term, and working memory? Prog. Brain Res. 2008 323 338 10.1016/S0079-6123(07)00020-9
6 Sederberg P.B. Schulze-Bonhage A. Madsen J.R. Bromfield E.B. Litt B. Brandt A. Kahana M.J. Gamma oscillations distinguish true from false memories Psychol. Sci. 18 2007 927 932 10.1111/j.1467-9280.2007.02003.x 17958703
7 Cohen J.D. Perlstein W.M. Braver T.S. Nystrom L.E. Noll D.C. Jonides J. Smith E.E. Temporal dynamics of brain activation during a working memory task Nature 386 1997 604 608 10.1038/386604a0 9121583
8 Johnson E.L. Dewar C.D. Solbakk A.-K. Endestad T. Meling T.R. Knight R.T. Bidirectional frontoparietal oscillatory systems support working memory 10.1016/j.cub.2017.05.046 2017
9 Roux F. Uhlhaas P.J. Working memory and neural oscillations: alpha–gamma versus theta–gamma codes for distinct WM information? Trends Cognit. Sci. 18 2014 16 25 10.1016/J.TICS.2013.10.010 24268290
10 Buzsáki G. Logothetis N. Singer W. Size Scaling Brain Timing Keeping Evolutionary preservation of brain rhythms Neuron 80 2013 751 764 10.1016/j.neuron.2013.10.002 24183025
11 Constantinidis C. Klingberg T. The neuroscience of working memory capacity and training Nat. Rev. Neurosci. 17 2016 438 449 10.1038/nrn.2016.43 27225070
12 Lisman J.E. Jensen O. The theta-gamma neural code Neuron 77 2013 1002 1016 10.1016/j.neuron.2013.03.007 23522038
13 Lisman J.E. Idiart M.A. Storage of 7 +/- 2 short-term memories in oscillatory subcycles Science 267 1995 1512 1515 10.1126/science.7878473 7878473
14 Šimko P. Kent J.A. Rektorova I. Is non-invasive brain stimulation effective for cognitive enhancement in Alzheimer's disease? An updated meta-analysis Clin. Neurophysiol. 144 2022 23 40 10.1016/j.clinph.2022.09.010 36215904
15 Pievani M. Pini L. Ferrari C. Pizzini F.B. Boscolo Galazzo I. Cobelli C. Cotelli M. Manenti R. Frisoni G.B. Coordinate-based meta-analysis of the default mode and salience network for target identification in non-invasive brain stimulation of Alzheimer's disease and behavioral variant Frontotemporal dementia networks J. Alzheim. Dis. 57 2017 825 843 10.3233/JAD-161105
16 Grover S. Wen W. Viswanathan V. Gill C.T. Reinhart R.M.G. Long-lasting, dissociable improvements in working memory and long-term memory in older adults with repetitive neuromodulation Nat. Neurosci. 25 2022 1237 1246 10.1038/s41593-022-01132-3 35995877
17 Reinhart R.M.G. Disruption and rescue of interareal theta phase coupling and adaptive behavior Proc. Natl. Acad. Sci. USA 114 2017 11542 11547 10.1073/pnas.1710257114 29073084
18 Pupíková M. Šimko P. Gajdoš M. Rektorová I. Modulation of working memory and resting-state fMRI by tDCS of the right frontoparietal network Neural Plast. 2021 2021 10.1155/2021/5594305
19 Di Lazzaro V. Bella R. Benussi A. Bologna M. Borroni B. Capone F. Chen K.-H.S. Chen R. V Chistyakov A. Classen J. Kiernan M.C. Koch G. Lanza G. Lefaucheur J.-P. Matsumoto H. Nguyen J.-P. Orth M. Pascual-Leone A. Rektorova I. Simko P. Taylor J.-P. Tremblay S. Ugawa Y. Dubbioso Af R. Ranieri F. Diagnostic contribution and therapeutic perspectives of transcranial magnetic stimulation in dementia Clin. Neurophysiol. 132 2021 2568 2607 10.1016/j.clinph.2021.05.035 34482205
20 Eliasova I. Anderkova L. Marecek R. Rektorova I. Non-invasive brain stimulation of the right inferior frontal gyrus may improve attention in early Alzheimer's disease: a pilot study J. Neurol. Sci. 346 2014 318 322 10.1016/j.jns.2014.08.036 25216556
21 Anderkova L. Eliasova I. Marecek R. Janousova E. Rektorova I. Distinct pattern of gray matter atrophy in mild Alzheimer's disease impacts on cognitive outcomes of noninvasive brain stimulation J. Alzheim. Dis. 48 2015 251 260 10.3233/JAD-150067
22 Cappon D. Jahanshahi M. Bisiacchi P. Value and efficacy of transcranial direct current stimulation in the cognitive rehabilitation: a critical review since 2000 Front. Neurosci. 10 2016 157 10.3389/fnins.2016.00157 27147949
23 Elder G.J. Taylor J.-P. Transcranial magnetic stimulation and transcranial direct current stimulation: treatments for cognitive and neuropsychiatric symptoms in the neurodegenerative dementias? Alzheimer's Res. Ther. 6 2014 74 10.1186/s13195-014-0074-1 25478032
24 Draaisma L.R. Wessel M.J. Moyne M. Morishita T. Hummel F.C. Targeting the frontoparietal network using bifocal transcranial alternating current stimulation during a motor sequence learning task in healthy older adults Brain Stimul. 15 2022 968 979 10.1016/j.brs.2022.06.012 35792318
25 Reinhart R.M.G. Nguyen J.A. Working memory revived in older adults by synchronizing rhythmic brain circuits Nat. Neurosci. 22 2019 820 827 10.1038/s41593-019-0371-x 30962628
26 Polanía R. Nitsche M.A. Korman C. Batsikadze G. Paulus W. The importance of timing in segregated theta phase-coupling for cognitive performance Curr. Biol. 22 2012 1314 1318 10.1016/j.cub.2012.05.021 22683259
27 Wolinski N. Cooper N.R. Sauseng P. Romei V. The speed of parietal theta frequency drives visuospatial working memory capacity PLoS Biol. 16 2018 1 17 10.1371/journal.pbio.2005348
28 Schwab B.C. Misselhorn J. Engel A.K. Modulation of large-scale cortical coupling by transcranial alternating current stimulation Brain Stimul. 12 2019 1187 1196 10.1016/j.brs.2019.04.013 31101568
29 Vieira P.G. Krause M.R. Pack C.C. tACS entrains neural activity while somatosensory input is blocked PLoS Biol. 18 2020 e3000834 10.1371/journal.pbio.3000834
30 Fries P. Neuron perspective rhythms for cognition: communication through coherence Neuron 88 2015 220 235 10.1016/j.neuron.2015.09.034 26447583
31 Deco G. Buehlmann A. Masquelier T. Hugues E. The role of rhythmic neural synchronization in rest and task conditions Front. Hum. Neurosci. 5 2011 1 6 10.3389/fnhum.2011.00004 21283556
32 Woods A.J. Antal A. Bikson M. Boggio P.S. Brunoni A.R. Celnik P. Cohen L.G. Fregni F. Herrmann C.S. Kappenman E.S. Knotkova H. Liebetanz D. Miniussi C. Miranda P.C. Paulus W. Priori A. Reato D. Stagg C. Wenderoth N. Nitsche M.A. A technical guide to tDCS, and related non-invasive brain stimulation tools Clin. Neurophysiol. 127 2016 1031 1048 10.1016/j.clinph.2015.11.012 26652115
33 Aktürk T. de Graaf T.A. Güntekin B. Hanoğlu L. Sack A.T. Enhancing memory capacity by experimentally slowing theta frequency oscillations using combined EEG-tACS Sci. Rep. 12 2022 14199 10.1038/s41598-022-18665-z
34 Živanović M. Bjekić J. Konstantinović U. Filipović S.R. Effects of online parietal transcranial electric stimulation on associative memory: a direct comparison between tDCS, theta tACS, and theta-oscillatory tDCS Sci. Rep. 12 2022 14091 10.1038/s41598-022-18376-5
35 Thielscher A. Antunes A. Saturnino G.B. Field modeling for transcranial magnetic stimulation: a useful tool to understand the physiological effects of TMS? 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) 2015 IEEE 222 225 10.1109/EMBC.2015.7318340
36 Canolty R.T. Edwards E. Dalal S.S. Soltani M. Nagarajan S.S. Kirsch H.E. Berger M.S. Barbare N.M. Knight R.T. High gamma power is phase-locked to theta oscillations in human neocortex Science 313 2006 1626 1628 10.1126/SCIENCE.1128115/SUPPL_FILE/CANOLTY_SOM.PDF 1979 16973878
37 Gao W. Lin W. Frontal parietal control network regulates the anti-correlated default and dorsal attention networks Hum. Brain Mapp. 2012 192 202 10.1002/hbm.21204 21391263
38 Thielscher A. Antunes A. Saturnino G.B. Field modeling for transcranial magnetic stimulation: a useful tool to understand the physiological effects of TMS? Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS 2015-Novem 2015 222 225 10.1109/EMBC.2015.7318340
39 Li L. Zhao D. Age-related inter-region EEG coupling changes during the control of bottom-Up and top-Down attention Front. Aging Neurosci. 7 2015 171929 10.3389/FNAGI.2015.00223/BIBTEX
40 Koechlin E. Hyafil A. Anterior prefrontal function and the limits of human decision-making Science 318 2007 594 598 10.1126/science.1142995 17962551
41 Alekseichuk I. Falchier A.Y. Linn G. Xu T. Milham M.P. Schroeder C.E. Opitz A. Electric field dynamics in the brain during multi-electrode transcranial electric stimulation Nat. Commun. 10 2019 2573 10.1038/s41467-019-10581-7 31189931
42 Clayton M.S. Yeung N. Cohen Kadosh R. The roles of cortical oscillations in sustained attention Trends Cognit. Sci. 19 2015 188 195 10.1016/J.TICS.2015.02.004 25765608
43 Violante I.R. Li L.M. Carmichael D.W. Lorenz R. Leech R. Hampshire A. Rothwell J.C. Sharp D.J. Externally induced frontoparietal synchronization modulates network dynamics and enhances working memory performance Elife 6 2017 1 22 10.7554/eLife.22001
44 Sauseng P. Klimesch W. Schabus M. Doppelmayr M. Fronto-parietal EEG coherence in theta and upper alpha reflect central executive functions of working memory Int. J. Psychophysiol. 57 2005 97 103 10.1016/J.IJPSYCHO.2005.03.018 15967528
45 Courtney S.M. Hinault T. When the time is right: temporal dynamics of brain activity in healthy aging and dementia Prog. Neurobiol. 203 2021 102076 10.1016/j.pneurobio.2021.102076
46 Tóth B. Kardos Z. File B. Boha R. Stam C.J. Molnár M. Frontal midline theta connectivity is related to efficiency of WM maintenance and is affected by aging Neurobiol. Learn. Mem. 114 2014 58 69 10.1016/J.NLM.2014.04.009 24758899
47 Alekseichuk I. Pabel S.C. Antal A. Paulus W. Intrahemispheric theta rhythm desynchronization impairs working memory Restor. Neurol. Neurosci. 35 2017 147 158 10.3233/RNN-160714 28059806
48 Röhner F. Breitling C. Rufener K.S. Heinze H.-J. Hinrichs H. Krauel K. Sweeney-Reed C.M. Modulation of working memory using transcranial electrical stimulation: a direct comparison between tacs and tdcs Front. Neurosci. 12 2018 10.3389/fnins.2018.00761
49 Meiron O. Lavidor M. Prefrontal oscillatory stimulation modulates access to cognitive control references in retrospective metacognitive commentary Clin. Neurophysiol. 125 2014 77 82 10.1016/j.clinph.2013.06.013 23831184
50 Hoy K.E. Bailey N. Arnold S. Windsor K. John J. Daskalakis Z.J. Fitzgerald P.B. The effect of γ-tACS on working memory performance in healthy controls Brain Cognit. 101 2015 51 56 10.1016/j.bandc.2015.11.002 26580743
51 Nissim N.R. McAfee D.C. Edwards S. Prato A. Lin J.X. Lu Z. Coslett H.B. Hamilton R.H. Efficacy of transcranial alternating current stimulation in the enhancement of working memory performance in healthy adults: a systematic meta-analysis Neuromodulation: Technology at the Neural Interface 26 2023 728 737 10.1016/j.neurom.2022.12.014 36759231
52 Zanto T.P. Rubens M.T. Thangavel A. Gazzaley A. Causal role of the prefrontal cortex in top-down modulation of visual processing and working memory Nat. Neurosci. 14 2011 656 661 10.1038/nn.2773 21441920
53 Perinelli A. Assecondi S. Tagliabue C.F. Mazza V. Power shift and connectivity changes in healthy aging during resting-state EEG Neuroimage 256 2022 119247 10.1016/j.neuroimage.2022.119247
54 Zangrossi A. Zanzotto G. Lorenzoni F. Indelicato G. Cannas Aghedu F. Cermelli P. Bisiacchi P.S. Resting-state functional brain connectivity predicts cognitive performance: an exploratory study on a time-based prospective memory task Behav. Brain Res. 402 2021 113130 10.1016/J.BBR.2021.113130
55 Abellaneda-Pérez K. Vaqué-Alcázar L. Perellón-Alfonso R. Bargalló N. Kuo M.F. Pascual-Leone A. Nitsche M.A. Bartrés-Faz D. Differential tDCS and tACS effects on working memory-related neural activity and resting-state connectivity Front. Neurosci. 13 2020 10.3389/fnins.2019.01440
56 Murphy A.C. Bertolero M.A. Papadopoulos L. Lydon-Staley D.M. Bassett D.S. Multimodal network dynamics underpinning working memory Nat. Commun. 11 1 11 2020 1 13 10.1038/s41467-020-15541-0 2020 31911652
57 Reynolds J.R. West R. Braver T. Distinct neural circuits support transient and sustained processes in prospective memory and working memory Cerebr. Cortex 19 2009 1208 1221 10.1093/CERCOR/BHN164
58 Wischnewski M. Berger T.A. Opitz A. Alekseichuk I. Causal functional maps of brain rhythms in working memory Proc. Natl. Acad. Sci. U.S.A. 121 2024 10.1073/PNAS.2318528121/-/DCSUPPLEMENTAL
59 Tversky A. Kahneman D. Judgment under uncertainty: heuristics and biases Science 185 1974 1124 1131 10.1126/science.185.4157.1124 17835457
60 Gigerenzer G. Czerlinski J. Martignon L. How good are fast and frugal heuristics? Decision Science and Technology 1999 81 103 10.1007/978-1-4615-5089-1_6
61 Reppert T.R. Heitz R.P. Schall J.D. Neural mechanisms for executive control of speed-accuracy trade-off Cell Rep. 42 2023 113422 10.1016/j.celrep.2023.113422
62 Standage D. Blohm G. Dorris M.C. On the neural implementation of the speed-accuracy trade-off Front. Neurosci. 8 2014 88630 10.3389/FNINS.2014.00236/BIBTEX
63 Forstmann B.U. Anwander A. Schäfer A. Neumann J. Brown S. Wagenmakers E.-J. Bogacz R. Turner R. Cortico-striatal connections predict control over speed and accuracy in perceptual decision making Proc. Natl. Acad. Sci. USA 107 2010 15916 15920 10.1073/pnas.1004932107 20733082
64 Weigard A. Beltz A. Reddy S.N. Wilson S.J. Characterizing the role of the pre‐SMA in the control of speed/accuracy trade‐off with directed functional connectivity mapping and multiple solution reduction Hum. Brain Mapp. 40 2019 1829 1843 10.1002/hbm.24493 30569619
65 van Veen V. Krug M.K. Carter C.S. The neural and computational basis of controlled speed-accuracy tradeoff during task performance J. Cognit. Neurosci. 20 2008 1952 1965 10.1162/jocn.2008.20146 18416686
66 Filmer H.L. Ballard T. Don K. Amarasekera R. Sewell D.K. Dux P.E. The causal role of the prefrontal and superior medial frontal cortices in the incidental manipulation of decision strategies Neuropsychologia 179 2023 108466 10.1016/j.neuropsychologia.2022.108466
67 Jones M.W. Wilson M.A. Theta rhythms coordinate hippocampal–prefrontal interactions in a spatial memory task PLoS Biol. 3 2005 e402 10.1371/JOURNAL.PBIO.0030402 16279838
68 Rutishauser U. Ross I.B. Mamelak A.N. Schuman E.M. Human memory strength is predicted by theta-frequency phase-locking of single neurons Nature 464 2010 903 907 10.1038/nature08860 20336071
69 Berger B. Griesmayr B. Minarik T. Biel A.L. Pinal D. Sterr A. Sauseng P. Dynamic regulation of interregional cortical communication by slow brain oscillations during working memory Nat. Commun. 10 2019 4242 10.1038/s41467-019-12057-0 31534123
70 Kasten F.H. Duecker K. Maack M.C. Meiser A. Herrmann C.S. Integrating electric field modeling and neuroimaging to explain inter-individual variability of tACS effects Nat. Commun. 10 2019 5427 10.1038/s41467-019-13417-6 31780668
71 Saturnino G.B. Madsen K.H. Siebner H.R. Thielscher A. How to target inter-regional phase synchronization with dual-site Transcranial Alternating Current Stimulation Neuroimage 163 2017 68 80 10.1016/J.NEUROIMAGE.2017.09.024 28919407
72 Neuling T. Ruhnau P. Weisz N. Herrmann C.S. Demarchi G. Faith and oscillations recovered: on analyzing EEG/MEG signals during tACS Neuroimage 147 2017 960 963 10.1016/j.neuroimage.2016.11.022 27888060
73 Holzmann R. Koppehele-Gossel J. Voss U. Klimke A. Investigating nuisance effects induced in EEG during tACS application Front. Hum. Neurosci. 15 2021 1 14 10.3389/fnhum.2021.637080
74 Alekseichuk I. Wischnewski M. Opitz A. A minimum effective dose for (transcranial) alternating current stimulation Brain Stimul. 15 2022 1221 1222 10.1016/j.brs.2022.08.018 36044976
75 Alam M. Truong D.Q. Khadka N. Bikson M. Spatial and polarity precision of concentric high-definition transcranial direct current stimulation (HD-tDCS) Phys. Med. Biol. 61 2016 4506 10.1088/0031-9155/61/12/4506 27223853
76 Mikkonen M. Laakso I. Tanaka S. Hirata A. Cost of focality in TDCS: interindividual variability in electric fields Brain Stimul. 13 2020 117 124 10.1016/J.BRS.2019.09.017 31606449
