Full text
52,479 characters
· extracted from
preprint-html
· click to expand
TMS-EEG reveals causal dynamics of the premotor cortex during musical improvisation | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results TMS-EEG reveals causal dynamics of the premotor cortex during musical improvisation View ORCID Profile Pedro T. Palhares , View ORCID Profile Sasha D’Ambrosio , View ORCID Profile Óscar F. Gonçalves doi: https://doi.org/10.1101/2025.10.04.680222 Pedro T. Palhares 1 Proaction Laboratory, Faculty of Psychology and Educational Sciences, University of Coimbra , Portugal 2 CINEICC, Faculty of Psychology and Educational Sciences, University of Coimbra , Portugal Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Pedro T. Palhares For correspondence: pedro.torres.palhares{at}gmail.com Sasha D’Ambrosio 3 Department of Biomedical and Clinical Sciences, University of Milan , Italy Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Sasha D’Ambrosio Óscar F. Gonçalves 4 RISE-HEALTH, Brainloop Laboratory — CINTESIS.UPT, Universidade Portucalense Infante D. Henrique , Portugal Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Óscar F. Gonçalves Abstract Full Text Info/History Metrics Preview PDF ABSTRACT Musical improvisation illustrates the brain’s capacity for flexible, creative motor control, yet the causal mechanisms underlying this complex behaviour remain poorly understood. We employed transcranial magnetic stimulation combined with electroencephalography (TMS-EEG) to probe state-dependent cortical dynamics in the left dorsal premotor cortex (PMd) of professional jazz pianists (n = 3) during improvisation, sight-reading, and rest. This proof-of-concept study demonstrates the feasibility of combining perturbational neuroscience with ecologically valid musical performance. Multiple convergent analyses revealed distinct cortical signatures during improvisation: reduced local mean field power, decreased phase-locking of evoked responses, and preserved but gain-modulated early components as revealed by Correlated Components Analysis. These findings suggest that improvisation is characterized by attenuated PMd excitability and more variable response timings, while preserving the fundamental architecture of cortical responses. This perturbational signature supports a neural efficiency model of expertise whereby expert musicians achieve creative flexibility through training-induced streamlined, optimized cortical processing. Our results establish TMS-EEG as a powerful approach for investigating the causal dynamics of creative cognition and demonstrate how the brain reconfigures its response properties to support internally driven motor performance. Introduction Musical performance represents one of the most sophisticated forms of human sensorimotor behaviour, requiring the precise coordination of perceptual, cognitive, and motor systems in real time ( Vuust et al., 2022 ; Zatorre et al., 2007 ). Within this domain, musical improvisation offers a particularly compelling model for understanding flexible motor control and creative cognition ( Palhares et al., 2024 ), demanding the spontaneous generation of novel motor sequences while maintaining temporal precision and adhering to complex harmonic constraints ( Pressing, 1988 ). Despite extensive neuroimaging research on the neural correlates of musical improvisation ( Beaty, 2015 ; Loui, 2018 ), the causal dynamics underlying this complex behaviour remain poorly understood, particularly regarding the functional role of premotor regions in orchestrating improvised (internally generated) versus externally cued motor performance. The premotor cortex (PMC) plays a critical role in motor planning, sequence generation, and the transformation of abstract goals into concrete motor commands ( Hoshi & Tanji, 2007 ; Rizzolatti et al., 1996 ). In the context of musical improvisation, neuroimaging studies have shown that the PMC undergoes functional reorganization, displaying altered connectivity patterns with other nodes of the motor network, including increased coupling with prefrontal control regions and sensorimotor areas ( Donnay et al., 2014 ; Pinho et al., 2014 ). Improvisation thus appears to engage a complex interplay between idea generation and evaluation processes, requiring the coordination of multiple brain networks which include auditory-motor regions for perception-action cycles and both medial and lateral prefrontal areas for creative cognition ( Beaty, 2015 ; Loui, 2018 ; Pinho et al., 2015 ). However, the correlational nature of neuroimaging data precludes strong causal inferences about how the PMC actively contributes to these distinct performance states ( Havlík et al., 2023 ; Logothetis, 2008 ; Miniussi & Thut, 2010 ). Transcranial magnetic stimulation combined with electroencephalography (TMS-EEG) offers a unique window into the causal dynamics of cortical function by directly perturbing neural activity and measuring the resulting spatiotemporal response patterns ( Ilmoniemi et al., 1997 ; Miniussi & Thut, 2010 ). Unlike traditional neuroimaging approaches that measure correlates of neural activity, TMS-EEG provides direct access to cortical excitability, connectivity, and state-dependent processing through the analysis of TMS-evoked potentials (TEPs). The morphology, amplitude, and spatiotemporal spread of TEPs reflect intrinsic properties of the stimulated cortical region and its connected networks, with evidence suggesting that these responses are sensitive to ongoing behavioural states ( Nikulin et al., 2003 ) and level of consciousness ( Sarasso et al., 2015 ). The distinction between improvised and sight-read musical performance in contrast with task-free resting state provides an ideal paradigm for investigating task-dependent variations in cortical responses to perturbation. Improvisation requires the integration of multiple sources of information — harmonic knowledge, motor memories, auditory feedback — into a coherent but flexible motor output. In contrast, sight-reading involves a more constrained transformation from visual input to motor output. How these different task demands are reflected in the brain’s response to direct perturbation remains an open question. In this study, we employed TMS-EEG to directly probe the state-dependent dynamics of the left dorsal premotor cortex during musical improvisation, sight-reading, and rest in professional jazz pianists. By analysing multiple complementary features of the TMS-evoked response — including local field power, phase-locking and maximally reproducible component structure — we aimed to characterize the PMC’s perturbational signature across these distinct performance states. More importantly, we aimed to demonstrate the feasibility of conducting TMS-EEG during a complex and ecologically valid motor task, extending the scope of TMS-EEG research beyond highly simplified motor tasks such as finger tapping. Materials & Methods Participants Five right-handed professional jazz pianists (1 female) were recruited for this study. All volunteers were screened for contraindications to TMS ( Rossi et al., 2011 ), including a history of traumatic brain injury, neurological or psychiatric disorders, intracranial metallic implants, drug-release dispensers, metallic tattoos, or cardiac pacemakers. Additional inclusion criteria included: ability to improvise over novel chord progressions on the spot, as depicted in jazz notation; ability to sight-read, as depicted in classical notation; having at least ten years of jazz improvisation experience. Two participants were excluded from the final analysis due to excessive EEG noise. The final sample thus comprised three participants (age: M = 30.0, SD = 3.5 years; see Table 1 ). View this table: View inline View popup Download powerpoint Table 1. Demographic data, musical background and information on experimental sessions for all participants Convenience sampling was used, with recruitment involving outreach to conservatories, university jazz programs, and other relevant institutions. Monetary compensation was not provided. All participants provided written informed consent, and the study was carried out under the ethical standards defined by the Institutional Ethics Committee and the Code of Ethics of the World Medical Association (Declaration of Helsinki). Experimental setup overview This study employed a purpose-designed multimodal TMS-EEG setup tailored for use during live piano performance. The experimental setup integrated four interconnected computers and multiple hardware units (EEG amplifier, TMS and neuronavigation), enabling real-time coordination of stimulation, electrophysiological, and neuronavigational components ( Figure 1 ). Participants performed on a digital piano while neural responses to stimulation and behavioural performance were recorded in parallel. One computer recorded musical performance and generated individualized TMS-click-masking stimuli ( Russo et al., 2022 ), which was mixed with the piano’s audio and delivered via headphones. A second computer acquired EEG data and enabled online monitoring of TMS-evoked potentials ( Casarotto et al., 2022 ). TMS delivery and cross-device synchronization were coordinated by a third computer controlling a central trigger interface. Neuronavigation was managed by a fourth computer, enabling real-time visualization and control of TMS targeting. Download figure Open in new tab Figure 1. Schematic representation of hardware and software integration TEP recording TEPs were recorded using a high-density gel-based cap (64-channel cap for TMS with multitrodes; EASYCAP, Germany) comprising sintered Ag/AgCl electrodes with coaxial cabling and arranged in the 10-20 layout with a forehead ground and reference montage. The cap was connected to a TMS-compatible 64-channel amplifier (eego EE-302; ANT Neuro, Netherlands). Impedances were kept below 10 kΩ, signals were band-pass filtered between 0.1–350 Hz and sampled at 8000 Hz with 24-bit resolution. To minimize contamination of TEPs by auditory responses to the TMS coil click ( ter Braack et al., 2015 ), participants wore noise-cancelling in-ear headphones through which individualized masking noise was continuously delivered. This noise was generated using the TAAC toolbox ( Russo et al., 2022 ), which adaptively reproduces the spectral and temporal features of the specific TMS click and superimposes them with calibrated broadband noise to ensure perceptual masking. TMS targeting and parametrization A focal figure-of-eight coil (mean/outer winding diameter ≈ 50/70 mm; biphasic pulse shape; pulse length ≈ 350 µs) connected to a mobile stimulation unit (Magstim Rapid 2 ; Magstim Co., UK) was used to deliver single-pulse stimulation. Stimulation targeted the dorsal left premotor cortex (BA 6). A Brainsight 3D infrared Tracking Position Sensor Unit (Rogue Research Inc., Canada) was used to co-register the participant’s scalp anatomy and the TMS coil relative to a standardized template MRI from Brainsight software. Stimulation was delivered within a 2 mm spatial error in relation to the ROI while maintaining angular deviations below 1° for both tilt and orientation. The identification of artefact-free stimulation targets within the ROI was guided by real-time visualization of TEPs using procedures adapted from the rt-TEP protocol ( Casarotto et al., 2022 ). While this protocol recommends tuning stimulation intensity based on the peak-to-peak amplitude of the earliest TEP component (typically within 0–30 ms), this approach was not feasible due to a prominent recharge artefact that extended up to ≈ 45 ms post-stimulation. As a result, early components were obscured and could not reliably be evaluated. We therefore adapted the procedure by initiating stimulation at each participant’s resting motor threshold and increasing intensity in steps, prioritizing the amplification of the earliest visible component after the artefact window (typically ≥ 45 ms). For each tested intensity, we explored different coil positions, tilts, and orientations within the ROI, using the rt-TEP interface to minimize muscle activity, eye movements, and discharge-related noise. Once no further improvements could be obtained through repositioning, the intensity was increased—unless doing so resulted in excessive artefact duration, in which case the previous setting was retained. The optimal combination of intensity, coil position, and orientation (i.e., the individualized stimulation “hotspot”) was then fixed in the neuronavigation system for the remainder of the session. All subsequent stimulation blocks for that participant — regardless of condition — were delivered using this exact same spatial and stimulation parameterization to ensure consistency. Stimulation intensity was fixed within each participant and ranged from 60% to 73% of maximum stimulator output (MSO), corresponding to 110% to 115% of their individual resting motor thresholds (RMT). After these preparatory steps, the experimental task began: TMS-evoked responses were recorded in three separate conditions, each comprising at least 200 pulses delivered with a randomly jittered inter-stimulus interval of 2000 to 2300 ms (0.4–0.5 Hz). Experimental task Each participant completed three experimental conditions: jazz improvisation, classical sight-reading, and resting state. All conditions were performed while seated at a digital piano. The first two involved active musical performance, while the third served as a passive baseline. In the improvisation condition ( Improv ), participants were asked to improvise freely on the piano, emulating a solo live performance. Each improvisation was guided by unique 16-bar chord progressions presented in standard jazz notation. These chord sequences were composed specifically for the study, designed to be of comparable difficulty, duration, and structure (4/4 time signature, 120 BPM). In the classical sight-reading condition ( Classical ), participants performed selected pieces from Béla Bartók’s Mikrokosmos . These were presented in classical notation and were novel to the participants, who were instructed to perform each piece as accurately and musically as possible at first sight. In the resting state condition ( Rest ), participants remained seated at the piano and were instructed to remain as still and relaxed as possible without engaging in any overt motor or cognitive task. TEPs were recorded for each condition in separate blocks, and the order of blocks was counterbalanced across participants. The duration of each block varied depending on the time required to deliver at least 200 single TMS pulses. Due to the strict spatial and angular tolerances enforced during stimulation, and the need to pause and readjust the coil when participant movements temporarily disrupted its stability, block durations ranged between 15 and 20 minutes. EEG data preprocessing Data preprocessing was performed on MATLAB R2024b (The MathWorks Inc.) using custom scripts. First, TMS-related artefacts were removed by replacing the signal between –2 and 50 ms relative to the pulse with a clean segment from –9 to –2 ms. This extended window ensured removal of both the pulse artefact and recharge artefact described earlier in the Methods section. Subsequently, trials were segmented from —800 to 800 ms around the stimulus onset and high-pass filtered at 1 Hz. Trials and channels were visually inspected and manually rejected when exhibiting excessive noise. Channels were then re-referenced to the average reference, after which independent component analysis (ICA) was applied to identify and remove residual artefacts related to eye movements, scalp muscle activity, and device interference (EEGLAB runica function; Delorme & Makeig, 2004 ). Interpolation of rejected channels was performed after ICA, followed by downsampling to 1 kHz. TEP analysis Analyses were implemented in MATLAB using custom code operating on EEGLAB-format datasets (EEG structs). We compared Improv, Classical and Rest in either single-participant or group mode. Unless noted, statistics are nonparametric and family-wise error rate (FWER) controlled via permutation tests ( α = .05). Divergence Index (DI) DI provided a preliminary global, time-resolved estimate of how much TMS-evoked EEG responses differed pairwise between conditions by quantifying the proportion of post-stimulus samples showing statistically significant amplitude differences across all channels ( Casarotto et al., 2010 ). After baseline handling via iterative Wilcoxon-based baseline matching ( α = .05), condition means were contrasted sample-by-sample from 50–250 ms post-stimulus. Signals were band-pass filtered at 1–45 Hz and baseline-subtracted using −250–0 ms. FWER across time was controlled with 10,000 label-shuffle permutations (time-wise max-statistic). The DI value is the percentage of post-stimulus samples declared significant. The empirical cutoff of 1.7% ( Casarotto et al., 2010 ) was used as a reference threshold. As an additional control, we computed a subject-specific split-half baseline distribution for Rest (1,000 random half-splits; 1,000 permutations for the threshold), providing a within-subject reference for expected DI under no systematic between-condition differences. Local mean field power (LMFP) LMFP assessed the strength of TMS-elicited currents over the four electrodes directly beneath the coil (FC1, FCz, C1, Cz) (a proxy for local cortical excitability; Fecchio et al., 2017 ). Per-trial signals were band-pass filtered at 1–45 Hz, baseline-corrected with the mean over −400 to −50 ms and transformed as the root-mean-square across the four ROI channels at each time point. Values were summarized in 50–300 ms. Within-condition enhancement relative to baseline was tested via one-sided sign-flip permutations across trials with family-wise error control over time using the maximum statistic ( α = .05; 5,000 permutations). Between-condition differences were tested with label-shuffle permutations (FWER over time via max-statistic; 5,000 permutations). Broadband phase-locking factor (bPLF) To quantify phase consistency of the evoked response, we computed PLF over the same 4-channel ROI in a broadband 8–45 Hz band. Trials were band-pass filtered (4th-order Butterworth), analytic phase was obtained with the Hilbert transform, PLF was computed per channel and averaged across the ROI, and statistics were restricted to 50–300 ms. Within-condition significance was assessed with phase-randomization permutations: on each permutation, every trial’s phase time series was offset by a random constant chosen uniformly from zero to one full cycle (0–360°), PLF was recomputed, and the maximum PLF within 50–300 ms served as the test statistic to control family-wise error over time ( α = .05; 5,000 permutations). Between-condition differences were tested with label-shuffle permutations at the trial level, using the max absolute PLF difference within 50–300 ms as the statistic (FWER by max-statistic; 5,000 permutations). Correlated Components Analysis (CorrCA) We used CorrCA to test whether TMS evokes common, reproducible network-level responses within and across conditions ( Parra et al., 2019 ). CorrCA identifies linear combinations of EEG channels whose time courses are maximally correlated across trials within a post-stimulus window. Mathematically, CorrCA solves a generalized eigendecomposition that maximizes the ratio of between to within-trial covariance, which is equivalent to maximizing inter-trial correlation. CorrCA was run on the four under-coil electrodes (FC1, FCz, C1, Cz) in a 50–150 ms window. Trials were entered into a regularized CorrCA, with the shrinkage parameter ( γ ) chosen by split-half cross-validation. Significance of inter-trial correlations (ISCs) was assessed with circular time-shift permutations (10,000 surrogates) using trimming-aware nulls and family-wise error control by the component-wise max statistic ( α =.05). “Transfer” was evaluated by applying filters learned in a reference condition to target conditions and testing the resulting ISCs against each target’s permutation distribution; evidence is reported as −log 10 ( p max ). Cross-condition correspondence was quantified two ways. First, we matched components one-to-one by absolute map correlation (| r map |) and compared their ERPs via Pearson correlation after sign-alignment. We report r ERP at zero lag and with a small-lag search (± 10 ms). Second, we did strict C1 to C1 comparisons: (i) | r map | for C1 forward maps, (ii) r ERP with max-lag (± 10 ms) after the same sign rule, (iii) the angle between C1 weight vectors, and (iv) a confirmatory waveform check after an orthogonal rotation within the shared subspace (L = 2). For figures, component traces are shown as mean ± SEM with a display-only demean on 48–52 ms; plotting normalization does not affect statistics. Transfer evidence for C1 was summarized in both directions (A to B and B to A) and we report the conservative single value as the minimum of the two −log 10 ( p max ) estimates. Results Divergence Index DI values revealed clear differences in TMS-evoked EEG responses across conditions and participants ( Figure 2 ). All pairwise contrasts showed DI values that substantially exceeded each participant’s split-half baseline threshold and the empirical cutoff of 1.7% ( Casarotto et al., 2010 ), suggesting that the observed differences reflect systematic condition effects rather than intrinsic variability. Download figure Open in new tab Figure 2. Divergence index by condition contrast and participant Local mean field power (LMFP) LMFP was quantified over the four electrodes directly beneath the coil (FC1, FCz, C1, Cz) across the 50–300 ms analysis window ( Figure 3 ). Overall power was highest for Rest ( M = 0.77 µV) and Classical ( M = 0.76 µV), and markedly lower for Improv ( M = 0.45 µV). Significant between-condition contrasts (two-sided FWER < .05), visualised as horizontal black bars beneath the LMFP traces, revealed divergence for Rest vs Improv (41.3% of the window) and Improv vs Classical (31.0%), whereas Rest vs Classical showed no reliable differences. Permutation tests confirmed that these responses were not baseline noise. Within-condition comparisons showed that 98.4% of Rest samples, 86.5% of Classical samples, and 36.5% of Improv samples exceeded their family-wise thresholds (one-sided FWER < .05). Thus, Rest and Classical elicited similarly strong cortical currents, while Improvisation was consistently weaker, driving all significant pairwise effects highlighted in the figure Download figure Open in new tab Figure 3. Local mean field power across conditions Note. The three panels show pairwise contrasts between conditions. Horizontal black bars beneath the traces denote time samples with significant differences between the two conditions in each panel. Local broadband phase-locking factor (bPLF) A broadband (8-45 Hz) PLF was quantified over the four electrodes directly beneath the coil (FC1, FCz, C1, Cz) across the 50–300 ms analysis window. All three conditions exhibited significant broadband phase alignment within 50–300 ms, but to markedly different degrees (see Figure 4 ). Rest showed the most sustained locking, with 43.7% of samples exceeding the permutation threshold (threshold = .097), a peak bPLF of .301, a 110 ms contiguous supra-threshold segment, and the largest area-under-the-curve (AUC = 7.33). Classical showed a similar profile: 34.1% significant (threshold = .102), the highest peak amplitude (.364), a 26 ms duration, and AUC = 7.77, whereas Improv yielded the weakest and least sustained locking (15.9% significant, threshold = .102; peak = .246; duration = 24 ms; AUC = 2.41). Pairwise permutation tests revealed limited inter-condition effects: Rest differed from Improv on 9.5% of the window (|threshold| = .101), Rest and Classical did not differ (0%), and Improvisation differed from Classical on just 4.8% of samples (|threshold| = .116). Overall, bPLF results indicate stronger and more sustained stimulus-locked phase synchrony during Rest and Classical than during Improvisation , with minimal differences between former two. Download figure Open in new tab Figure 4. Local broadband phase-locking factor across conditions Note . The three panels show pairwise contrasts between conditions. Horizontal black bars beneath the traces denote time samples with significant differences between the two conditions in each panel. Correlated Components Analysis (CorrCA) We applied CorrCA over the four electrodes directly beneath the coil (FC1, FCz, C1, Cz) across a focused 50– 150 ms analysis window to avoid later sensory components. CorrCA recovered reliable components in all conditions, with markedly stronger within-condition structure in Rest and Classical than in Improv . The leading inter-trial correlation was comparable for Rest (ISC 1 ≈ .106) and Classical (ISC 1 ≈ .106), but smaller for Improv (ISC 1 ≈ .060). Permutation testing confirmed multiple significant components (max-statistic: Rest = 3, Classical = 2, Improv = 1). A γ sweep showed ISC 1 was stable to regularization, indicating these effects are not driven by the shrinkage choice. Cross-condition transfer for C1 was at the permutation floor for 10,000 permutations (−log 10 ( p max ) = 4.00) for all pairings, indicating very strong generalization. For visualization of Component 1 (C1), baseline-corrected overlays are shown in Figure 5a , and raw-amplitude overlays under a fixed Rest projection ( W ref ) are shown in Figure 5b . Download figure Open in new tab Figure 5. Component 1 (shape-first) — baseline-corrected overlays. (b) Component 1 (amplitude-first) — overlays under fixed Rest projection Note . (a) Each trace represents CorrCA Component 1 per condition, demeaned per trial on 48–52 ms and polarity-aligned to emphasize waveform morphology across conditions, independent of overall amplitude. Displays mean ± SEM in µV. (b) All conditions are projected with the Rest spatial filter ( W ref ) and shown in raw µV with polarity aligned. This display highlights between-condition amplitude differences while preserving a comparable time course. Using Rest as reference, cross-condition matching by forward maps yielded strong topographic correspondence for C1 with Classical (| r map | = 0.97) and substantial correspondence with Improv (| r map | = 0.86). ERP similarity at ± 10 ms was high for Classical ( r ≈ 0.62) and moderate for Improv ( r ≈ 0.37) in the shape-first view ( Figure 5a ), but when projecting all conditions with W ref the C1, ERPs converged for both targets ( r ≈ 0.59–0.64), and a simple L = 2 subspace alignment further equalized them ( r ≈ 0.62; values summarized in Table 2 ). Consistent with this, C1 transfer remained at the permutation floor (−log 10 ( p max ) = 4; see Figure 6 , “Transfer”). Principal-angle analyses further support a rotation account: Classical ’s first weight dimension is nearly parallel to Rest ( θ 1 ≈ 1.8°), whereas Improv shares the first dimension more modestly ( θ 1 ≈ 16.7°) and diverges in the second ( θ 2 ≈ 62.8°), consistent with a low-dimensional rotation rather than a wholesale change. View this table: View inline View popup Download powerpoint Table 2. Summary of CorrCA C1 similarity and generalization across conditions Download figure Open in new tab Figure 6. Component 1 similarity summary (pairwise) Note . Bars summarize (from top to bottom): | r map | between forward maps, r ERP (max-lag, +/-10 ms), and transfer evidence reported as −log 10 ( p max ) using its minimum value. White circles mark best-match overlays (C1 paired to the target’s most similar component by | r map |). Metrics not shown as bars (angle and aligned r ERP ) are reported in Table 2 . A strict pairwise check that fixes the index to “C1 vs C1” highlights this rotation and is summarized in the bar panels ( Figure 6 ): Rest – Improv shows a very low C1-to-C1 map correlation (| r map | = 0.08) even though ERP similarity ( r ≈ 0.61) and bidirectional transfer evidence (−log 10 ( p max ) = 4) remain strong; Rest – Classical shows high agreement on all metrics (| r map | = 0.97; r ≈ 0.62; −log 10 ( p max ) = 4). After within-subspace alignment, the ERP similarity for the “C1 process” recovers to high values ( r ≈ 0.62–0.73) even for pairs involving Improv ( Table 2 ), indicating that Improv ’s C1 is effectively a rotated/reindexed combination of the first few Rest components. In short, the dominant response around 50–150 ms is preserved across conditions. Discrepancies with Improv in strict C1-to-C1 map comparisons reflect component mixing within the shared subspace, not a polarity issue or a loss of the underlying effect. Discussion This study employed TMS-EEG to investigate state-dependent dynamics of the left dorsal premotor cortex during musical improvisation, sight-reading, and rest in professional jazz pianists. Our findings reveal distinct perturbational signatures that differentiate improvised performance from both sight-reading and resting states, providing novel causal insights into how the PMC dynamically reconfigures its response properties to support different modes of motor control. Reduced local excitability during improvisation Our most striking finding was the marked reduction in LMFP during improvisation compared to both rest and classical sight-reading conditions. The approximate 41% reduction relative to rest and classical performance, with no significant difference between rest and classical conditions, suggests that improvisation engages the PMC in a fundamentally distinct excitability state, one in which TMS perturbation elicits substantially smaller local activation. This pattern aligns with a state-dependent gating mechanism whereby improvisation places PMd in a mode that selectively attenuates responses to exogenous, time-locked stimulation while maintaining ongoing, internally guided computations. Previous TMS-EEG studies have demonstrated similar dampening of evoked potentials during active motor states compared to rest, including reduced N100 and late positive components over motor areas during simple movement execution or inhibition ( Nikulin et al., 2003 ; Yamanaka et al., 2013 ), consistent with decreased susceptibility of local circuits when they are behaviourally engaged. Our findings extend this principle to complex creative motor behaviour, suggesting that the task-dependent reduction in cortical responsiveness scales with the cognitive demands and internal focus required for improvisation. The reduced LMFP during improvisation converges with neuroimaging evidence supporting a neural efficiency account of musical improvisation expertise. Expert improvisers show paradoxically lower BOLD activation across frontoparietal executive regions alongside stronger functional connectivity between bilateral DLPFC, PMd, and pre-SMA ( Pinho et al., 2014 , 2015 ). This pattern has been interpreted as reflecting the automation and streamlining of creative motor operations through optimized premotor-prefrontal coupling ( Beaty, 2015 ; Beaty et al., 2016 ). This finding also aligns with Haslinger and colleagues’ (2004) finding of reduced activation of PMd (among other motor network regions) in concert pianists during bimanual coordination, which was argued to be indicative of a training-induced efficiency of cortical systems that could be fundamental to achieve high-level motor skills. Preserved early components with reduced gain Despite the overall reduction in evoked power, CorrCA (which maximizes inter-trial correlation, equivalent to maximizing repeat-reliability of the evoked waveform; Parra et al., 2019 ) revealed that the fundamental structure of the early cortical response (50-150 ms) remains preserved across all conditions. The dominant component showed strong generalization between rest and classical conditions (| r map | = 0.97) and substantial correspondence even with improvisation after accounting for subspace rotation. This dissociation (i.e., preserved component morphology with reduced amplitude) suggests that improvisation does not fundamentally alter the canonical premotor response architecture but rather modulates its gain. This finding reconciles heterogeneity in the literature regarding premotor involvement in improvisation. While some studies report increased PMd activation during improvisation ( de Manzano & Ullén, 2012 ), others show decreased activity with expertise ( Pinho et al., 2014 ). Our results suggest that both observations may reflect different aspects of the same phenomenon: the PMd remains engaged in the same fundamental computations (preserved component structure) but operates in a more efficient, less excitable regime (reduced LMFP). This interpretation aligns with theories proposing that expertise involves not the recruitment of different neural circuits but rather the optimization of existing ones ( Beaty, 2015 ; Beaty et al., 2016 ; Pinho et al., 2014 ). Reduced broadband phase-locking This reduced PLF indicates that while TMS still evokes a reproducible response during improvisation (as shown by CorrCA), the timing of this response is more variable across trials. During rest and classical sight-reading, the PMd maintains more stereotyped response timing to external perturbation, with each TMS pulse triggering a response at a consistent phase. In contrast, during improvisation, the temporal relationship between stimulation and response becomes less deterministic. This finding aligns with the reduced LMFP in suggesting that improvisation engages a cortical state less dominated by stimulus-locked responses. The variable response timing during improvisation may indicate that ongoing endogenous dynamics play a stronger role in shaping when and how the cortex responds to perturbation. This could reflect: (i) stronger ongoing oscillations that interact differently with the TMS pulse on each trial, (ii) a more complex state space where perturbations can follow multiple trajectories, or (iii) active gating mechanisms that modulate response timing based on the current phase of internal processes. Together with our other local measures, this temporal variability further characterizes the distinct cortical state associated with improvisation, one that maintains responsiveness (preserved components) but with reduced amplitude (LMFP) and less predictable timing (PLF). Because inter-trial phase metrics can be influenced by differences in evoked amplitude and latency variability ( van Diepen & Mazaheri, 2018 ), we interpret the PLF attenuation with appropriate caution. Nevertheless, the use of phase-randomization and family-wise error statistical controls, together with the concordant LMFP decrease in the same ROI, makes a purely power-driven account unlikely. Limitations and future directions Several limitations warrant consideration. First, our sample size of three participants limits the generalizability of our findings. While the high internal validity of TMS-EEG studies, arising from the large number of observations per participant per condition, allows us to establish effects with near certainty within our sample, extending these findings to the broader population of jazz musicians warrants caution. The demanding nature of these experimental sessions, which lasted up to five hours and required participants to perform complex musical tasks while maintaining precise coil positioning, made larger-scale data collection exceptionally challenging. This study thus represents a proof-of-concept and methodological innovation rather than a systematic investigation, demonstrating the feasibility of combining TMS-EEG with real musical performance despite the considerable technical and practical obstacles. Second, the presence of recharge artefacts limited our analysis to responses beyond 50 ms. Current-generation TMS systems with customisable coil recharge timing would eliminate this limitation, allowing future investigations to probe these earlier components. Additionally, while our focus on the left PMd was motivated by its established role in internally generated movement, improvisation likely involves distributed networks extending beyond this single region. Future studies employing multiple stimulation sites could map how different nodes of the motor network reconfigure during creative performance. Conclusions This study provides initial causal evidence that musical improvisation fundamentally alters the response properties of the premotor cortex to external perturbation. The convergent findings of reduced local excitability, preserved but gain-modulated components, and decreased phase-locking support a neural efficiency model in which improvisation engages a distinct cortical state optimized for internally driven motor control. More broadly, our findings demonstrate the value of combining perturbational approaches with ecologically valid tasks to understand the neural basis of complex human abilities. By directly probing how cortical circuits respond to external perturbation during real musical performance, we reveal state-dependent dynamics that would remain hidden in traditional activation or correlational studies. This approach opens new avenues for understanding how the brain flexibly reconfigures its computational properties to support the human capacity for creative expression. Research transparency statement Funding This work was supported by FCT – Fundação para a Ciência e Tecnologia, I.P. by project reference and DOI identifier https://doi.org/10.54499/2022.12647.BD . Data availability Datasets, analysis scripts and materials supporting this study will be made available upon publication. Ethics approval statement This study was approved by the Ethics Committee of Faculty of Psychology and Educational Sciences, University of Coimbra (reference code: CEDI/FPCEUC:64/3). Conflict of interest disclosure The authors declare no conflicts of interest. Funder Information Declared Fundação para a Ciência e Tecnologia, https://ror.org/00snfqn58 , 2022.12647.BD Footnotes The only change in this revision is the removal of line numbers, which were interfering with automatic indexing processes. References ↵ Beaty , R. E. ( 2015 ). The neuroscience of musical improvisation . Neuroscience & Biobehavioral Reviews , 51 , 108 – 117 . doi: 10.1016/j.neubiorev.2015.01.004 OpenUrl CrossRef PubMed ↵ Beaty , R. E. , Benedek , M. , Silvia , P. J. , & Schacter , D. L. ( 2016 ). Creative cognition and brain network dynamics . Trends in Cognitive Sciences , 20 ( 2 ), 87 – 95 . doi: 10.1016/j.tics.2015.10.004 OpenUrl CrossRef PubMed ↵ Casarotto , S. , Fecchio , M. , Rosanova , M. , Varone , G. , D’Ambrosio , S. , Sarasso , S. , Pigorini , A. , Russo , S. , Comanducci , A. , Ilmoniemi , R. J. , & Massimini , M. ( 2022 ). The rt-TEP tool: Real-time visualization of TMS-Evoked Potentials to maximize cortical activation and minimize artifacts . Journal of Neuroscience Methods , 370 , 109486 . doi: 10.1016/j.jneumeth.2022.109486 OpenUrl CrossRef ↵ Casarotto , S. , Romero Lauro , L. J. , Bellina , V. , Casali , A. G. , Rosanova , M. , Pigorini , A. , Defendi , S. , Mariotti , M. , & Massimini , M. ( 2010 ). EEG responses to TMS are sensitive to changes in the perturbation parameters and repeatable over time . PLoS ONE , 5 ( 4 ), e10281 . doi: 10.1371/journal.pone.0010281 OpenUrl CrossRef PubMed ↵ de Manzano , Ö. , & Ullén , F. ( 2012 ). Goal-independent mechanisms for free response generation: Creative and pseudorandom performance share neural substrates . NeuroImage , 59 ( 1 ), 772 – 780 . doi: 10.1016/j.neuroimage.2011.07.016 OpenUrl CrossRef PubMed ↵ Delorme , A. , & Makeig , S. ( 2004 ). EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis . Journal of Neuroscience Methods , 134 ( 1 ), 9 – 21 . doi: 10.1016/j.jneumeth.2003.10.009 OpenUrl CrossRef PubMed Web of Science ↵ Donnay , G. F. , Rankin , S. K. , Lopez-Gonzalez , M. , Jiradejvong , P. , & Limb , C. J. ( 2014 ). Neural substrates of interactive musical improvisation: An fMRI study of ‘trading fours’ in jazz . PLOS ONE , 9 ( 2 ), e88665 . doi: 10.1371/journal.pone.0088665 OpenUrl CrossRef PubMed ↵ Fecchio , M. , Pigorini , A. , Comanducci , A. , Sarasso , S. , Casarotto , S. , Premoli , I. , Derchi , C.-C. , Mazza , A. , Russo , S. , Resta , F. , Ferrarelli , F. , Mariotti , M. , Ziemann , U. , Massimini , M. , & Rosanova , M. ( 2017 ). The spectral features of EEG responses to transcranial magnetic stimulation of the primary motor cortex depend on the amplitude of the motor evoked potentials . PLOS ONE , 12 ( 9 ), e0184910 . doi: 10.1371/journal.pone.0184910 OpenUrl CrossRef PubMed ↵ Haslinger , B. , Erhard , P. , Altenmüller , E. , Hennenlotter , A. , Schwaiger , M. , Gräfin von Einsiedel , H. , Rummeny , E. , Conrad , B. , & Ceballos-Baumann , A. O. ( 2004 ). Reduced recruitment of motor association areas during bimanual coordination in concert pianists . Human Brain Mapping , 22 ( 3 ), 206 – 215 . doi: 10.1002/hbm.20028 OpenUrl CrossRef PubMed Web of Science ↵ Havlík , M. , Hlinka , J. , Klírová , M. , Adámek , P. , & Horáček , J. ( 2023 ). Towards causal mechanisms of consciousness through focused transcranial brain stimulation . Neuroscience of Consciousness , 2023 ( 1 ), niad008 . doi: 10.1093/nc/niad008 OpenUrl CrossRef PubMed ↵ Hoshi , E. , & Tanji , J. ( 2007 ). Distinctions between dorsal and ventral premotor areas: Anatomical connectivity and functional properties . Current Opinion in Neurobiology , 17 ( 2 ), 234 – 242 . doi: 10.1016/j.conb.2007.02.003 OpenUrl CrossRef PubMed Web of Science ↵ Ilmoniemi , R. J. , Virtanen , J. , Ruohonen , J. , Karhu , J. , Aronen , H. J. , Näätänen , R. , & Katila , T. ( 1997 ). Neuronal responses to magnetic stimulation reveal cortical reactivity and connectivity . NeuroReport , 8 ( 16 ), 3537 . OpenUrl CrossRef PubMed Web of Science ↵ Logothetis , N. K. ( 2008 ). What we can do and what we cannot do with fMRI . Nature , 453 ( 7197 ), 869 – 878 . doi: 10.1038/nature06976 OpenUrl CrossRef PubMed Web of Science ↵ Loui , P. ( 2018 ). Rapid and flexible creativity in musical improvisation: Review and a model . Annals of the New York Academy of Sciences , 1423 ( 1 ), 138 – 145 . doi: 10.1111/nyas.13628 OpenUrl CrossRef ↵ Miniussi , C. , & Thut , G. ( 2010 ). Combining TMS and EEG offers new prospects in cognitive neuroscience . Brain Topography , 22 ( 4 ), 249 – 256 . doi: 10.1007/s10548-009-0083-8 OpenUrl CrossRef PubMed Web of Science ↵ Nikulin , V. V. , Kičić , D. , Kähkönen , S. , & Ilmoniemi , R. J. ( 2003 ). Modulation of electroencephalographic responses to transcranial magnetic stimulation: Evidence for changes in cortical excitability related to movement . European Journal of Neuroscience , 18 ( 5 ), 1206 – 1212 . doi: 10.1046/j.1460-9568.2003.02858.x OpenUrl CrossRef PubMed Web of Science ↵ Palhares , P. T. , Sas , M. I. , & Gonçalves , Ó. F. ( 2024 ). Music and states of consciousness: A narrative review of the broader significance of music to understanding absorption, mind wandering and creative thought . Neuroscience & Biobehavioral Reviews , 167 , 105920 . doi: 10.1016/j.neubiorev.2024.105920 OpenUrl CrossRef PubMed ↵ Parra , L. C. , Haufe , S. , & Dmochowski , J. P. ( 2019 ). Correlated Components Analysis—Extracting reliable dimensions in multivariate data . doi: 10.51628/001c.7125 OpenUrl CrossRef ↵ Pinho , A. L. , de Manzano , O. , Fransson , P. , Eriksson , H. , & Ullen , F. ( 2014 ). Connecting to create: Expertise in musical improvisation is associated with increased functional connectivity between premotor and prefrontal areas . Journal of Neuroscience , 34 ( 18 ), 6156 – 6163 . doi: 10.1523/JNEUROSCI.4769-13.2014 OpenUrl Abstract / FREE Full Text ↵ Pinho , A. L. , Ullén , F. , Castelo-Branco , M. , Fransson , P. , & de Manzano , Ö. ( 2015 ). Addressing a paradox: Dual strategies for creative performance in introspective and extrospective networks . Cerebral Cortex , 26 ( 7 ), 3052 – 3063 . doi: 10.1093/cercor/bhv130 OpenUrl CrossRef PubMed ↵ J. Sloboda Pressing , J. ( 1988 ). Improvisation: Methods and models . In J. Sloboda , Generative processes in music: The psychology of performance, improvisation, and composition (pp. 129 – 178 ). Oxford University Press . doi: 10.1093/acprof:oso/9780198508465.001.0001 OpenUrl CrossRef ↵ Rizzolatti , G. , Fadiga , L. , Gallese , V. , & Fogassi , L. ( 1996 ). Premotor cortex and the recognition of motor actions . Cognitive Brain Research , 3 ( 2 ), 131 – 141 . doi: 10.1016/0926-6410(95)00038-0 OpenUrl CrossRef PubMed ↵ Rossi , S. , Hallett , M. , Rossini , P. M. , & Pascual-Leone , A. ( 2011 ). Screening questionnaire before TMS: An update . Clinical Neurophysiology , 122 ( 8 ), 1686 . doi: 10.1016/j.clinph.2010.12.037 OpenUrl CrossRef PubMed ↵ Russo , S. , Sarasso , S. , Puglisi , G. E. , Dal Palù , D. , Pigorini , A. , Casarotto , S. , D’Ambrosio , S. , Astolfi , A. , Massimini , M. , Rosanova , M. , & Fecchio , M. ( 2022 ). TAAC - TMS Adaptable Auditory Control: A universal tool to mask TMS clicks . Journal of Neuroscience Methods , 370 , 109491 . doi: 10.1016/j.jneumeth.2022.109491 OpenUrl CrossRef PubMed ↵ Sarasso , S. , Boly , M. , Napolitani , M. , Gosseries , O. , Charland-Verville , V. , Casarotto , S. , Rosanova , M. , Casali , A. G. , Brichant , J.-F. , Boveroux , P. , Rex , S. , Tononi , G. , Laureys , S. , & Massimini , M. ( 2015 ). Consciousness and complexity during unresponsiveness induced by propofol, xenon, and ketamine . Current Biology , 25 ( 23 ), 3099 – 3105 . doi: 10.1016/j.cub.2015.10.014 OpenUrl CrossRef PubMed ↵ ter Braack , E. M. , de Vos , C. C. , & van Putten , M. J. A. M. ( 2015 ). Masking the auditory evoked potential in TMS–EEG: A comparison of various methods . Brain Topography , 28 ( 3 ), 520 – 528 . doi: 10.1007/s10548-013-0312-z OpenUrl CrossRef PubMed ↵ van Diepen , R. M. , & Mazaheri , A. ( 2018 ). The caveats of observing inter-trial phase-coherence in cognitive neuroscience . Scientific Reports , 8 , 2990 . doi: 10.1038/s41598-018-20423-z OpenUrl CrossRef PubMed ↵ Vuust , P. , Heggli , O. A. , Friston , K. J. , & Kringelbach , M. L. ( 2022 ). Music in the brain . Nature Reviews Neuroscience , 23 ( 5 ), 287 – 305 . doi: 10.1038/s41583-022-00578-5 OpenUrl CrossRef PubMed ↵ Yamanaka , K. , Kadota , H. , & Nozaki , D. ( 2013 ). Long-latency TMS-evoked potentials during motor execution and inhibition . Frontiers in Human Neuroscience , 7 , 751 . doi: 10.3389/fnhum.2013.00751 OpenUrl CrossRef ↵ Zatorre , R. J. , Chen , J. L. , & Penhune , V. B. ( 2007 ). When the brain plays music: Auditory–motor interactions in music perception and production . Nature Reviews Neuroscience , 8 ( 7 ), 547 – 558 . doi: 10.1038/nrn2152 OpenUrl CrossRef PubMed Web of Science View the discussion thread. Back to top Previous Next Posted October 05, 2025. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following TMS-EEG reveals causal dynamics of the premotor cortex during musical improvisation Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share TMS-EEG reveals causal dynamics of the premotor cortex during musical improvisation Pedro T. Palhares , Sasha D’Ambrosio , Óscar F. Gonçalves bioRxiv 2025.10.04.680222; doi: https://doi.org/10.1101/2025.10.04.680222 Share This Article: Copy Citation Tools TMS-EEG reveals causal dynamics of the premotor cortex during musical improvisation Pedro T. Palhares , Sasha D’Ambrosio , Óscar F. Gonçalves bioRxiv 2025.10.04.680222; doi: https://doi.org/10.1101/2025.10.04.680222 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Neuroscience Subject Areas All Articles Animal Behavior and Cognition (7617) Biochemistry (17633) Bioengineering (13856) Bioinformatics (41841) Biophysics (21399) Cancer Biology (18529) Cell Biology (25422) Clinical Trials (138) Developmental Biology (13352) Ecology (19860) Epidemiology (2067) Evolutionary Biology (24281) Genetics (15582) Genomics (22461) Immunology (17700) Microbiology (40293) Molecular Biology (17140) Neuroscience (88413) Paleontology (666) Pathology (2823) Pharmacology and Toxicology (4813) Physiology (7632) Plant Biology (15107) Scientific Communication and Education (2042) Synthetic Biology (4284) Systems Biology (9808) Zoology (2267)
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.