“…For each EEG derivation, power spectral density (PSD) estimates were computed using the Welch's method with 2-s Hanning windows and 50% overlap (applying the pwelch function in the MATLAB signal processing toolbox). For each patient, average PSD was calculated over the same four channels used for TMS/EEG analysis (contralesional, perilesional): average PSD across segments was computed for each frequency bin and then further averaged across bins pertaining to the classical frequency ranges: delta (0.5-4.5Hz), theta (4.5-8 Hz), alpha (8)(9)(10)(11)(12), beta (12)(13)(14)(15)(16)(17)(18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29)(30) and gamma (30-40 Hz) ( Figure S1A). In parallel to the PSD analysis, filtered continuous EEG data were re-referenced according to a longitudinal bipolar montage based on the 10-20 system (78) including the following EEG derivations: 'Fp1_F5', 'F7_T3', 'T3_T5', 'T5_O1' 'Fp1_F3', 'F3_C3', 'C3_P3', 'P3_O1', 'Fz_Cz', 'Cz_Pz' ,'Fp2_F6', 'F6_T4', 'T4_T6', 'T6_O2', 'Fp2_F4', 'F4_C4', 'C4_P4', and 'P4_O2'.…”