Healthy VolunteersNeuroimaging & Brain MeasuresPsilocybinPlacebo

EEG microstate dynamics during psilocybin intoxication relate to acute experience and persisting psychological changes

This double-blind, randomised, placebo-controlled crossover study (n=15) in healthy volunteers found that psilocybin altered the timing of EEG brain-state patterns during peak intoxication, with faster switching between states but largely preserved overall state coverage. These changes were linked to the intensity of the acute experience and to self-reported psychological changes 28 days later.

Authors

  • Tomáš Páleníček
  • František Tylš
  • Martin Brunovský

Published

Biorxiv
individual Study

Abstract

Psilocybin and other serotonergic psychedelics show therapeutic promise for psychiatric disorders, yet objective neural correlates linking the acute psychedelic state to persisting psychological outcomes remain limited. Electroencephalography (EEG) microstate analysis characterizes the rapid spatiotemporal organization of large-scale brain activity, offering a millisecond-resolution window into neural dynamics. Here, we examined resting-state EEG microstates in 15 healthy volunteers who participated in a double-blind, randomized, placebo-controlled crossover study of psilocybin, using both data-driven (three-microstate) and canonical (four-microstate) analysis solutions. EEG was recorded at five time points spanning pre-drug baseline, peak intoxication, and recovery. Psilocybin significantly increased the number of global field power (GFP) peaks and reduced microstate lifespan while increasing frequency of occurrence during peak intoxication (50-100 min post-administration), consistent with accelerated transitions between brain states. Notably, microstate coverage was largely preserved, with only a transient difference at peak intoxication in the 2-20Hz band-width, suggesting that access to the repertoire of canonical brain states is broadly maintained despite altered temporal dynamics. Critically, individual differences in microstate dynamics during peak intoxication correlated with both acute subjective experience intensity and self-reported psychological changes measured 28 days post-administration, providing exploratory evidence for a link between acute neural dynamics and longer-term experiential outcomes in healthy volunteers. These findings suggest that psilocybin is associated with altered temporal organization of large-scale brain dynamics with largely preserved microstate coverage, and identify EEG microstates as candidate neural markers for psychedelic-induced alterations in consciousness with potential relevance to therapeutic research.

Unlocked with Blossom Pro

Research Summary of 'EEG microstate dynamics during psilocybin intoxication relate to acute experience and persisting psychological changes'

Editorial

βBlossom's Take

This (pre-print) study adds one more piece to the puzzle of how psychedelics work. Though measures across studies are still quite differentiated, studies like this help fill in some of the details.

Introduction

The authors note that psilocybin and other serotonergic psychedelics have renewed interest as potential treatments for psychiatric disorders, with clinical trials suggesting rapid and sometimes sustained benefits after only one or two administrations. However, they also emphasise that the neural mechanisms linking the acute psychedelic state to later psychological change remain incompletely understood. Earlier research has shown that the subjective quality of the psychedelic experience can predict longer-term outcomes, but relying on self-report alone limits mechanistic inference. At the neural level, prior work has focused more on spectral power, entropy, and large-scale network changes than on the fast spatiotemporal organisation of EEG activity. Against this background, the study set out to examine whether resting-state EEG microstate dynamics during psilocybin intoxication are altered in healthy volunteers, and whether those acute changes relate to both subjective experience and persisting psychological changes 28 days later. The researchers aimed to characterise psilocybin’s effects on microstate temporal organisation and to test whether individual differences in those dynamics track acute phenomenology and later self-reported change. They present this as an exploratory attempt to evaluate EEG microstates as candidate neural markers of psychedelic-induced changes in consciousness with possible therapeutic relevance.

Methods

The study analysed data from a double-blind, randomised, placebo-controlled crossover trial in healthy volunteers. Participants were screened carefully for neurological, psychiatric, substance-use, and major medical disorders; those with personal or first-degree family histories of psychotic illness were excluded. For the present analysis, the authors used the first 20 enrolled participants, but 5 were excluded because of missing sessions or poor-quality EEG, leaving 15 participants in the final dataset. The sessions were separated by at least 28 days. Participants received a weight-adjusted oral dose of pharmaceutical-grade psilocybin of approximately 0.26 mg/kg, or visually identical placebo capsules containing wheat starch only. Each person completed two experimental days, one with psilocybin and one with placebo, in counterbalanced order. During the sessions, they rested in a reclined position under supervision. EEG was recorded at five time points: a 10-minute pre-dose baseline and four 10-minute recordings at 50-60, 90-100, 180-190, and 360-370 minutes after ingestion, covering peak intoxication and recovery. Blood samples for psilocin levels were collected at baseline and at 1, 2, 4, and 6 hours after dosing. Subjective experience was measured with the Altered States of Consciousness questionnaire at the end of each session, acute clinical symptoms were rated repeatedly with the Brief Psychiatric Rating Scale during the session, and persisting effects were assessed 28 days later with the Persisting Effects Questionnaire. EEG was recorded from a 21-channel cap, sampled at 1000 Hz, and preprocessed using filtering, independent component analysis, artefact rejection, average re-referencing, and downsampling. The cleaned data were segmented into 2-second epochs, yielding 40 seconds of clean signal per participant per condition per time point. Microstate analysis was carried out in Python using GFP peaks, which identify time points with the highest signal-to-noise ratio, followed by hierarchical clustering to derive microstate maps. The authors analysed both a three-microstate solution and the canonical four-microstate solution, and repeated all analyses in two frequency bands: 1-40 Hz and 2-20 Hz. The three-microstate 1-40 Hz solution was treated as primary because the cross-validation criterion suggested three states, while the 2-20 Hz four-state solution was kept as a secondary, more traditional comparison. Outcome measures included microstate lifespan, coverage, frequency of occurrence, and transition probabilities, although transition probabilities were calculated for completeness and not analysed statistically. Repeated-measures ANOVA with Benjamini-Hochberg false discovery rate correction was used for GFP peak counts and microstate features, followed by paired t-tests where appropriate. Relationships between microstate measures and experience data were examined using Spearman correlations with false discovery rate correction. To reduce the number of correlations, the researchers aggregated experience variables and microstate features using principal component analysis and correlated the first principal component from each domain.

Results

The psilocybin condition significantly increased GFP peak density compared with placebo, especially during peak intoxication. In the primary 1-40 Hz analysis, condition, time, and their interaction were significant, and the number of GFP peaks was higher at T2, T3, and T4 under psilocybin. In the 2-20 Hz analysis, condition was also significant, and the T3 psilocybin condition showed a higher number of GFP peaks than placebo. Cross-validation indicated that three microstates was the optimal solution across conditions and time points, although the authors also reported the canonical four-state solution. The canonical microstate maps resembled the expected template maps in both frequency bands, with group-level mean correlations above 0.8. For explained variance, the 2-20 Hz four-state solution showed no significant time effect, but in the 1-40 Hz analysis condition was significant: explained variance was slightly higher in placebo than in psilocybin, with an unbiased Cohen’s d of 0.414. The main temporal microstate effects were reductions in lifespan and increases in frequency of occurrence during psilocybin intoxication, most clearly at T2 and T3, when the subjective effects were strongest. These changes were seen in both frequency-band analyses, with strong effect sizes for some microstate-specific comparisons. Coverage was comparatively stable: no significant differences were found in the 1-40 Hz analysis, and in the 2-20 Hz analysis there was a single significant difference at T2. Overall, the authors describe coverage as largely preserved despite altered timing dynamics. Subjectively, all four Altered States of Consciousness dimensions were significantly elevated under psilocybin relative to placebo. Brief Psychiatric Rating Scale scores also increased during intoxication, particularly for thought disturbance and hallucinations, withdrawal and retardation, tension and excitement, and hostile suspiciousness at 70 minutes; anxiety and depression did not differ significantly. The authors also report that these subjective and observer-rated effects aligned with blood psilocin levels and dose, with correlations showing coherent clustering among the psychometric and pharmacokinetic measures. Finally, microstate measures during peak intoxication correlated with experiential outcomes. The aggregated microstate component was associated with acute subjective experience, and at T2 it also correlated with positive persisting effects on the Persisting Effects Questionnaire 28 days later. Negative persisting effects were low overall, with mean scores below 10% of the scale maximum, and were not significantly associated with microstate dynamics.

Discussion

The authors interpret the findings as showing that psilocybin accelerates and destabilises resting-state EEG microstate dynamics during peak intoxication. They argue that the increased number of GFP peaks and the shortened microstate lifespan indicate more rapid switching between global brain states, consistent with earlier work showing that psychedelics increase neural entropy, signal diversity, and large-scale network reorganisation. In their view, microstate peak density may serve as a topography-based complement to complexity metrics such as Lempel-Ziv complexity. They emphasise that the temporal profile of these changes followed psilocybin pharmacokinetics, peaking at T2 and T3 alongside peak plasma psilocin levels and maximal subjective effects. The increased frequency of occurrence is presented as a necessary counterpart to shorter lifespans and as further evidence of faster state transitions. The authors link these observations to broader theories suggesting that psychedelics increase the rate at which the brain explores its state space. By contrast, microstate coverage was largely unchanged, which the authors interpret as suggesting that psilocybin mainly alters the sequencing and stability of brain states rather than the repertoire of states available. They note that this pattern may distinguish the psychedelic state from psychiatric conditions such as schizophrenia and depression, where altered coverage has been reported, although they caution that their small sample limits sensitivity to smaller coverage effects. The authors also see translational interest in the correlations between acute microstate dynamics, subjective intensity, and persisting positive effects 28 days later. They suggest that greater microstate destabilisation during the acute state was associated with stronger later changes in mood, well-being, and personal meaning. However, they stress that these are correlational findings and do not establish causality or therapeutic efficacy. They say that validation in clinical populations is still required before any treatment-related inference can be drawn. In discussing limitations, the authors highlight the small sample size (N=15), which reduces power for time-resolved and smaller effects, particularly in the secondary analysis. They also raise an alternative explanation for the EEG findings: psilocybin-related increases in heart rate, pupil dilation, and sympathetic arousal could potentially contribute to faster EEG fluctuations. Because physiological signals were not simultaneously recorded, they could not rule this out. The authors further state that future work should include concurrent physiological monitoring, clinical samples, multimodal imaging, pharmacological manipulations, and computational modelling to clarify what microstates reflect mechanistically. They conclude that microstate measures are promising pragmatic biomarkers of altered brain dynamics, but their functional interpretation remains incomplete.

Conclusion

The authors conclude that psilocybin markedly reorganised resting-state EEG microstate dynamics, increasing GFP peak density and microstate switching while shortening microstate lifespan, whereas microstate coverage was largely preserved. They state that these findings are consistent with an acceleration and destabilisation of large-scale brain dynamics rather than a major change in the set of accessible brain states. They also conclude that acute microstate differences were associated with both subjective psychedelic experience and persisting positive psychological effects, suggesting potential value as neural markers of psychedelic-associated changes in consciousness, while noting that mechanistic interpretation remains limited.

View full paper sections

SECTION

Alterations in microstate dynamics have been reported across a range of psychiatric and neurological disorders. Schizophrenia and major depression, for example, are associated with reduced microstate duration, altered coverage, and abnormal transition patterns, suggesting that microstates capture clinically meaningful aspects of network dysfunction. More recently, abnormal microstate profiles have been linked to neurodegenerative processes and cognitive decline. These findings position microstates as sensitive markers of large-scale brain dynamics with potential translational relevance. Despite this, the effects of psychedelics on EEG microstate organization have received little systematic investigation. Given that psilocybin profoundly alters conscious experience, destabilizes canonical network structure, and increases neural entropy, we hypothesized that it would also substantially modulate microstate dynamics. We also explored, in a hypothesis-generating fashion, whether inter-individual differences in microstate dynamics during the acute state would relate to self-reported psychological changes assessed 28 days postadministration. To test these hypotheses, we examined resting-state EEG microstates in healthy volunteers undergoing placebo and psilocybin sessions. Our primary aims were (i) to characterize the acute effects of psilocybin on microstate temporal organization, and (ii) to determine whether individual differences in microstate dynamics during the psychedelic state relate to both subjective experience and long-term psychological outcomes previously reported in this cohort. By linking millisecond-scale brain dynamics to both phenomenology and persisting effects, this work seeks to evaluate EEG microstates as candidate neural markers of psychedelic-induced changes in consciousness and their potential therapeutic relevance.

STUDY APPROVAL

The study was approved by the local ethics committee of the Prague Psychiatric Centre/National Institute of Mental Health and by the Czech legal authority, the State Institute for Drug Control. It was approved as a clinical trial registered under the EudraCT No. 2012-004579-37. All volunteers signed an informed consent prior to entering the study.

STUDY DESIGN, PARTICIPANTS, AND RECRUITMENT

Data for the present analysis were collected during the first phase of a larger clinical trial consisting of two arms. In each arm, subjects underwent one psilocybin and one placebo session, using a double-blind, randomized, crossover design. While the first arm focused primarily on EEG data collection, the second arm was designed for fMRI data collection. Sessions were separated by at least a 28-day interval. The comprehensive structure of the full trial is detailed in our previous work. Healthy volunteers were recruited via word-of-mouth and community referral. To ensure the exclusion of neurological, psychiatric, or major medical disorders-as well as current or past substance misuse-all participants underwent a rigorous screening process. This included a medical somatoneurological examination, laboratory testing, the Mini-International Neuropsychiatric Interview (M.I.N.I.), and the Minnesota Multiphasic Personality Inventory-2 (MMPI-2). Candidates with a personal or first-degree family history of psychotic illness were excluded. All participants abstained from alcohol and psychoactive substances for at least 48 hours before each session. For the analyses presented here, we utilized data from the first 20 healthy subjects (8 females; aged 20-35 years). This cohort was recorded using a 21-channel EEG and is identical to the sample described in. Of the 20 enrolled participants, 5 were excluded from the final analyses due to data acquisition problems (N=3 missing sessions; N=2 incomplete EEG recordings due to poor signal quality). All analyses reported here are therefore based on N = 15 participants.

DRUG AND DOSAGE

Participants received a weight-adjusted oral dose of approximately 0.26 mg/kg of pharmaceutical-grade psilocybin encapsulated in gelatin capsules with Amylum tritici (wheat starch) as a binder. Placebo capsules contained Amylum tritici (wheat starch) only and were visually indistinguishable. Capsules containing either 1 or 5 mg of psilocybin were prepared at the IKEM (Institute of Clinical and Experimental Medicine in Prague) pharmacy. The dose of psilocybin was adjusted by combining capsules containing 1 and 5 mg, increasing or decreasing by 1 mg per 5 kg of body weight, with a 75 kg person receiving 20 mg of psilocybin.

EXPERIMENTAL DESIGN

Each participant completed two experimental days (psilocybin, hereafter PSI, and placebo, hereafter PLA) in a counterbalanced order. Throughout the dosing sessions, participants were accompanied by a pair of assisting sitters who served supportive roles and actively participated in data collection. A trained nurse/EEG technician was also present throughout the sessions. Upon arrival, baseline physiological measurements and questionnaires were collected. Subjects received an intravenous cannula and were fitted with the EEG gel base cap, and the session began. After receiving the drug or placebo, participants rested comfortably in a reclined position under the continuous supervision of assisting sitters. Resting-state eyes-closed EEG was recorded at five predefined time points during the whole session: a baseline 10-minute recording acquired before drug ingestion (T1), followed by four 10-minute recordings at 50-60 minutes (T2), 90-100 minutes (T3), 180-190 minutes (T4), and 360-370 minutes (T5) after ingestion. These intervals captured pre-drug baseline, peak and late intoxication, and the post-acute recovery period. No stimulation was presented; participants were instructed to keep their eyes closed and remain still and relaxed. In between the resting-state EEG recordings, subjects listened to music, underwent EEG oddball paradigmsand auditory steady-state responses (ASSR), and were repeatedly examined with the Brief Psychiatric Rating Scale (BPRS). Blood samples for psilocin plasma levels were collected at 5 time points (baseline, 1h, 2h, 4h, and 6h after ingestion) throughout the dosing session.

SUBJECTIVE AND BEHAVIOURAL MEASURES

Subjective and behavioral effects of psilocybin were assessed using a combination of standardized psychometric instruments that capture acute drug effects, clinical symptomatology during intoxication, and persisting subjective effects after the session. Acute subjective drug effects were evaluated using the Altered States of Consciousness (ASC) questionnaire, which participants completed at the end of each experimental session. The ASC comprises 72 items rated on visual analog scales and yields scores across three main dimensions -Oceanic Boundlessness (OSE), Dread of Ego Dissolution (AIA), and Visionary Restructuralization (VUS), as well as a fourth composite score reflecting the overall intensity of the altered state of waking consciousness (Veränderter Wachbewusstseinszustand, VWB), providing a multidimensional characterization of the psychedelic experience. Objective clinical and behavioral changes during the session were assessed using the BPRS. The BPRS was administered repeatedly over the course of each session by assisting sitters, specifically before drug/placebo intake (T=0 min) and at 70 and 180 minutes after ingestion. The scale rates 18 psychiatric symptoms, each on a 7-point scale (0-6, ranging from "not present" to "very severe"), which can be clustered into five symptom domains: FI (anxiety, depression), FII (withdrawal, retardation), FIII (thought disturbance and hallucinations), FIV (tension, excitement), and FV (hostile suspiciousness), as described in. Persisting subjective effects were assessed 28 days after each experimental sessionusing the Persisting Effects Questionnaire (PEQ). The PEQ captures persisting subjective changes attributed to the psychedelic experience, including alterations in mood, attitudes, behavior, well-being, and perceived personal meaning. Items are rated on Likert-type scales and summarized into positive and negative persisting-effect domains, allowing quantification of persisting experiential impact beyond the acute intoxication phase.

PSILOCIN PLASMA LEVELS

Blood samples were drawn via an indwelling cannula inserted before the start of each session. The samples were collected at 5 time points: before ingestion (baseline) and at 1, 2, 4, and 6 hours after ingestion of the capsules. Blood samples were centrifuged at room temperature for 10 min at 4000 rpm. The separated sera were then stored at -20 °C until analyzed. Psilocin in sera was analyzed by gas chromatography/mass spectrometry (GC-MS). The analysis itself was performed on a GC HP model 6890 A with 5973 MSD and capillary HP5-MS. For details on the analyses, see our previous work.

EEG RECORDING

EEG signals were recorded using a 21-channel Ag/AgCl cap arranged according to the extended 10-20 system (original 10-20 nomenclature; ElectroCap International, USA). Four bipolar EOG channels monitored eye movements. Signals were amplified with a BioSDA09 digital EEG amplifier (M&I, Prague, Czech Republic), sampled at 1000 Hz. Electrode impedances were kept below 5 kΩ. All recordings were obtained in an electrically shielded, sound-attenuated room with participants resting in a reclined position.

EEG DATA PREPROCESSING

EEG preprocessing was performed in BrainVision Analyzer 2.1.1 (Brain Products GmbH, Gilching, Germany). Raw data were visually and semi-automatically inspected for artifacts. Signals were filtered using zero-phase infinite impulse response (IIR) Butterworth filters (0.5-100 Hz band-pass, 8th order) together with a 50 Hz notch filter. Independent component analysis (ICA; Fas-tICA, 24 components; 21 scalp EEG channels and four EOG channels used as input) was applied to the continuous data. Artifactual ICA components, primarily those reflecting ocular activity, were identified semi-automatically and removed via inverse ICA before data reconstruction. Following ICA cleaning, EEG was re-referenced to the average of 19 scalp electrodes (Fp1, Fp2, F3, F4, F7, F8, Fz, C3, C4, Cz, T3, T4, T5, T6, P3, P4, Pz, O1, O2). The data were then downsampled to 250 Hz using spline interpolation after anti-alias filtering at 112.5 Hz (24 dB/oct). Continuous EEG recordings were segmented into non-overlapping 2-second epochs, with all previously marked artifact intervals excluded. Only clean epochs were retained for microstate analysis. Preprocessed EEG was exported in ASCII format for subsequent computational processing in Python.

MICROSTATE ANALYSIS

EEG microstate analysis was performed for each participant and condition separately using Python. Within microstate analysis, the multichannel EEG signal is viewed as a series of instantaneous topographies of potential. We extracted the global field power (GFP) peaks in order to account only for temporal points with the highest signal-to-noise ratio. The GFP at each temporal instant is equal to the root mean square across the average-referenced electrodes (i.e. a standard deviation of the signal): with v i (t) being the voltage at electrode i at time instant t, v(t) is the mean voltage across all electrodes in a given time instant t, and n is the number of electrodes. Scalp potential topographies (or maps) at GFP time-series maxima represent the highest field strength and the greatest signal-to-noise ratio. The microstate analysis then includes two steps: first, identifying a set of microstate maps; second, projecting the original multichannel EEG recording onto the basis of microstate maps, thereby converting the EEG signal into a sequence of microstate maps (see Supplementary Figurefor an overview). The maps at the maxima of GFP time-series were submitted into a modified hierarchical clustering algorithm. Briefly, all maps are initially considered as independent clusters and in each iteration, the worst cluster is identified and dissolved; its constituent maps are redistributed to the remaining clusters according to the strongest Pearson product-moment correlation in the form where sums are taken over i electrodes and u and v represent two topographies. These steps are repeated until the desired number of clusters remains (see e.g.). We used the above-described microstate algorithm to cluster original maps from each participant and condition into microstate maps, specific for each subject and condition. Next, in order to obtain conditionspecific maps, we submitted maps from each participant from given condition to another round of clustering separately for each condition, thus obtained mean maps for given condition. All maps (both individual and condition-specific) were labeled class A, B, C, and D based on resemblance (i.e. Pearson product-moment correlation) to normative template maps following Koenig et al.(Supplementary Figure). After finding canonical microstate maps for each subject and each condition, the original multichannel EEG recording was transformed into a temporal sequence of microstate maps by the following procedure: one of the canonical maps is assigned for each GFP peak based on the highest absolute Pearson product-moment correlation (i.e., map polarity is disregarded) and if two subsequent GFP peaks are assigned different maps, we identified the temporal midpoint and created a microstate transition there. No minimum segment duration threshold or temporal smoothing was applied during back-projection. From the sequence of microstate maps, we calculated the following features: Average lifespan of microstates. Lifespan of microstates is calculated as the time during which all successive original topographies were assigned the same microstate class(cf. Supplementary Figure). Microstate coverage. The coverage of a microstate is calculated by taking a ratio of the total time spent in that particular microstate over total recording time. Frequency of occurrence. Microstate frequency is determined by counting unique appearances of each microstate in one second of the recording. Microstate transition probabilities. The number of transitions from each microstate class into other microstate class is counted and normalised to fractions of all between-class transitions, yielding a N × N matrix with zero diagonal of transition probabilities. Transition probabilities were computed for completeness but are not analyzed statistically in the present study, as the primary focus is on temporal (lifespan, frequency of occurrence) and spatial (coverage) microstate parameters. After the EEG data were preprocessed, divided into epochs, and artifacts removed, we were left with 20 epochs of 2 s duration, totaling 40 seconds of the signal per participant, per condition (placebo: PLA vs. psilocybin: PSI ), and time (T1 through T5 ). We then submitted these data to the microstate analysis described above. All our subsequent analyses rely on two different filtering approaches: since participants in our study were administered psilocybin, which shifts relative spectral maximum towards faster frequencies, we opted for a wider bandwidth, that is 1-40 Hz. On the other hand, in order to be able to compare our results with the existing body of literature, we repeated all analyses with a more "traditional" microstate filtering option, that is 2-20 Hz. A fundamental methodological decision concerned the optimal number of microstates to extract. The cross-validation (CV) criterionindicated three microstates as optimal; however, the canonical approach in the microstate literature employs four microstates. As noted by Murray et al., the CV criterion can underestimate the true number of classes in datasets with lower electrode density, and a recent rat study similarly found it to favour three microstates. We therefore analysed both three-and fourmicrostate solutions across both frequency bands, reporting the 1-40 Hz / three-microstate solution as primary (supported by the CV criterion) and the 2-20 Hz / four-microstate solution as secondary (following traditional methodology).

STATISTICAL ANALYSES

To compare the number of GFP peaks across times and conditions, we used repeated-measures ANOVA, followed by pairwise repeated-measures t-test post hoc with the correction for multiple comparison using the Benjamini-Hochberg step-up procedure to control the False Discovery Rate (FDR). Similar procedure 6 (repeated-measures ANOVA + pairwise repeated-measures t-test post hoc) was employed to compare microstate features (lifespan, coverage, and frequency of occurrence). To assess relationships between subject experience data, we used Spearman's rank correlation coefficient ρ to mitigate non-normality, with significance assessed using the Benjamini-Hochberg FDR procedure to control for multiple comparisons across all correlation pairs within each matrix. To minimize the number of correlations between experience data and microstate features and streamline the comparison, we aggregated both the experience data and microstate features using principal component analysis (PCA) and used the first PCA component, explaining the largest proportion of variance in the whole dataset. The first PCA component in the case of microstate features explains approximately 80% of the variance, whereas in the case of experience data it ranges between 54 and 72%.

MICROSTATE ANALYSIS 3.1.1. GLOBAL FIELD POWER AND IDEAL NUMBER OF MICROSTATES

First, we compared the number of peaks in GFP curves between conditions. Our hypothesis was that the number of peaks in the psilocybin condition would be higher for relevant times (i.e. T2, T3, and, possibly, T4 ), since classical psychedelics tend to shift spectral maxima towards higher frequencies, hence the signal itself is faster and therefore we expected the number of peaks to be higher. This was validated using repeated measures ANOVA which identified condition as a significant factor in 2-20 Hz bandwidth, while all (condition, time, and their interaction) factors were regarded as significant in 1-40 Hz bandwidth. Figureillustrates this point: in the right-hand side of the figure (1-40 Hz bandwidth) the number of GFP peaks in psilocybin condition is indeed significantly higher than placebo in T2, T3, and T4 times. The number of GFP peaks for psilocybin T3 condition is significantly higher in 2-20 Hz bandwidth as well (left-hand side of Figure). Figurealso indicates significant differences as revealed by pairwise t-tests with Benjamini-Hochberg procedureemployed in order to control false discovery rate (FDR). Higher number of GFP peaks results in higher number of topographic maps submitted into the microstate clustering algorithm. Traditionally, the microstate studies compute 4 canonical microstates and refer to the seminal early works in this area (e.g. Lehmann et al.or Pascual-Marqui et al.). Given that microstate analysis has received limited systematic investigation in the context of psychedelic EEG, we also tested the ideal number of microstates as obtained from the cross-validation (CV) test proposed by Pascual-Marqui et al.. Cross-validation criterion states that the ideal number of microstates is the number q which minimises modified cross-validation variance σ2 M CV given by the following function: where σ2 µ is the microstate model nonpredictive residual variance for q different microstates and N s is the number of electrodes. We computed the σ2 M CV curve for each data snippet with changing number of states q ∈ {1, . . . , 10} in all times and conditions and then found the minimum within each curve. Figuresummarises the ideal number of states (q that minimises the σ2 M CV curve) in our dataset. Unexpectedly, the ideal number of microstates in each condition and time across the subject (aggregating with median) was 3. Overall, we observe the trend towards lower number of states in wider, 1-40 Hz, bandwidth, however the ideal number of states was 3 for both bandwidths. Based on these results, we decided to continue our analyses within two branches: the main branch is 1-40 Hz bandwidth with decomposition into 3 canonical microstates, while the second is "classical" microstate analysis, that is 2-20 Hz bandwidth with 4 canonical states.

MICROSTATE MAPS

After initial EEG signal analysis with respect to the number of GFP peaks and ideal number of microstates within each bandwidth, we present the actual microstate topographies in Figuresandin supplementary information. The time-condition specific topographies, i.e. the average microstates per condition and time, were obtained firstly by running microstate clustering algorithm on all our data (that is for each participant, each condition, and each time) and all maps from one group were subsequently submitted to another round of clustering. The maps for 1-40 Hz bandwidth, for which the ideal number of microstates was 3 (cf. Figure), are shown in Figurein supplementary information for each of the time-condition combination. Below each map we also report the Pearson product-moment correlation (see eq. 2) between particular map and microstate "template maps " provided by Koenig et al.. These template maps serve as the gold standard in microstate analysis. The correlations indicate (group-level mean Pearson product-moment correlations per bandwidth, per condition, and per time are all C > 0.8; individual maps may fall below this threshold) that all our maps resemble template maps to a high degree and therefore we are able to take the analysis to the next step and compute some microstate statistics using these canonical maps. Similarly as for wider bandwidth, the 2-20 Hz bandwidth maps (see Figurein supplementary information) were compared with template maps from Koenig et al.. We conclude that in 2-20 Hz bandwidth all our maps bear a resemblance to template and the overall correlation is even higher than in the 1-40 Hz case. Consequently, we take also this branch of our analysis to the next step and submit the data to microstate statistics computation. We also tested whether there is a significant difference in explained variance in GFP peaks by our microstate decomposition. In the 2-20 Hz bandwidth with 4 canonical states there was no significant difference between times (F (4, 56) = 2.042, p = 0. microstates, repeated measures ANOVA showed significant dependence on condition with F (1, 14) = 5.675, p = 0.032, where the variance explained in placebo condition was slightly higher (0.707 ± 0.065, mean ± SD) than in the psilocybin condition (0.682 ± 0.055) (unbiased Cohen's d = 0.414). Other factors were not identified as significant (time: F (4, 56) = 2.521, p = 0.074 and interaction condition × time: F (4, 56) = 2.751, p = 0.066).

MICROSTATE FEATURES

After performing base analysis of GFP peaks, cross-validation test for ideal number of canonical maps and visualising the canonical maps themselves for both our filtering paradigms we computed typical microstate features which are the average lifespan of microstates, microstate coverage, frequency of occurrence, and their transition probabilities. Significant differences between placebo and psilocybin condition are marked with vertical lines and their respective p-values as based on pairwise t-test post hoc (preceded by repeated measures ANOVA) and are corrected for multiple comparisons using Benjamini-Hochberg FDR procedure. For tabulated results please refer to Table. The results of average microstate lifespan for both filtering paradigms are shown in Figurewhile the significant differences in average lifespan per conditions, microstate, and time are described in Table. Generally, and in line with our initial hypothesis, the significant differences between placebo and psilocybin 9 conditions are mainly in the T2 and T3 times, when the effect of the intoxication is the strongest with psilocybin condition tending to lower average lifespan, consistent with faster state switching. This is true for both filtering paradigms. Similar results (faster state switching) were observed in interaction with particular microstates: microstates B and C in 1-40 Hz bandwidth. The effect sizes as measured with Cohen's d were strong in this interaction setting (cf. Table; Supplementary Figure). Lastly, even individual microstate factor while disregarding time factor was significantly different in the 1-40 Hz filtering paradigm (microstates A, B, and C). Similarly to the lifespan, coverage results are depicted in Supplementary Figure. In the 2-20 Hz bandwidth, a significant difference in coverage between conditions was found at T2 (pairwise t-test post hoc, p = 0.015, Benjamini-Hochberg FDR corrected). No significant differences were found for any other time or microstate combination, nor in the 1-40 Hz bandwidth. Overall, as the figure shows, the means of coverage were similar across both factors and their interaction, while the interquartile ranges were narrower in the psilocybin condition, particularly at T2 and T3.

PLACEBO VS. PSILOCYBIN DIFFERENCES IN LIFESPAN

Our final measure of interest was the frequency of occurrence. The main results are plotted in Figure. Since frequency of occurrence measures similar underlying dynamics as the average lifespan (the relationship between the two is inverse), we expected similar outcomes as with lifespan. As expected, the main differences between two conditions are observed in T2 and T3 times; moreover, the effect still persists into T4 in the 1-40 Hz filtering paradigm. Similar results were observed in interaction with particular microstates: all three microstates in the 1-40 Hz bandwidth and microstate A in 2-20 Hz bandwidth. The effect sizes as measured with Cohen's d were strong in this interaction setting (cf. Table). Lastly, even individual microstate factor while disregarding time factor was significantly different in the 1-40 Hz filtering paradigm (all three microstates A, B, and C).

RELATIONSHIP BETWEEN MICROSTATES AND PSYCHEDELIC PHENOMENOLOGY

Given the pronounced effects of psilocybin on microstate temporal dynamics, we next examined whether these neural changes relate to subjective experience and clinical ratings collected during the experimental sessions. As expected given the established psychotropic effects of psilocybin, all four ASC dimensions were significantly elevated in the psilocybin condition compared to placebo: Oceanic Boundlessness (OSE: t(Similarly, the Brief Psychiatric Rating Scale revealed significant condition-related differences across most symptom domains during peak intoxication. At baseline (T=0min), no differences were observed between conditions. However, at 70 and 180 minutes post-administration, BPRS factors showed substantial elevations. under psilocybin, particularly in thought disturbance and hallucinations (FIII), withdrawal and retardation (FII), and tension and excitement (FIV); anxiety and depression (FI) showed no significant difference at either time point; hostile suspiciousness (FV) was significantly elevated at T=70 min (p = 0.033, d = 0.861) but not at T=180 min (p = 0.164). Full statistical details are reported in Table. To examine how self-reported subjective experiences (as captured by the ASC scales) relate to observerrated clinical symptoms (as captured by the BPRS) during the psilocybin session, we computed a correlation matrix including psilocin concentration in the blood at 5 different times during the experiment and the dose in mg per participant. The matrix is shown in Figure. As clearly seen from the correlation matrix, individual types of measures (ASC, BPRS at three different times) are all internally consistent, hence they create clusters of high positive correlation in the grand matrix. Moreover, ASC scales correlate with BPRS scales at 70 and 180 minutes into the experiment, and all scales correlate with psilocin concentration in the blood, in particular after 60 and 120 minutes from the start of the experiment. This confirmed that participants experienced a robust altered state of consciousness, as reflected across both self-reported and observer-rated measures and in plasma psilocin concentration. To identify which (if any) microstate feature relates to the changes in psychedelic scales, we plotted Figureas the correlation matrix between aggregated scales (means of ASC, BPRS at relevant times -70 and 180 minutes of the experiment), dose, and psilocin concentration on the one hand and aggregated microstate features on the other. Aggregation via PCA followed the same procedure described in the Statistical Analyses section. To assess whether acute microstate changes relate to persisting psychological effects, we examined correlations between microstate features during peak intoxication (T2, T3) and PEQ scores obtained 28 days post-administration. As shown in Figure, the first PCA component of microstate features correlated with positive persisting effects, particularly for the T2 time point. Notably, this relationship was specific to positive persisting effects; negative persisting effects showed low overall scores (mean negative PEQ score < 10% of scale maximum) and no significant associations with microstate dynamics.

DISCUSSION

The rationale for analyzing both three-and four-microstate solutions across two frequency bands is detailed in the Methods section. Both approaches yielded qualitatively consistent conclusions across the main outcome measures, supporting the robustness of our findings. GFP peak density was significantly elevated at T2 and T3 (and at T4 in the 1-40 Hz analysis; Figure), corresponding to peak subjective effects. Increased GFP peak density may reflect more frequent transitions between moments of high signal-to-noise ratio and, by extension, more rapid changes in global topographic configurations. This finding aligns with a large body of work showing that psychedelics increase neural signal diversity and entropy, as quantified by Lempel-Ziv complexity and related metrics in EEG and MEG. Computational modeling further supports this view, demonstrating that whole-brain models can reproduce entropy increases under psychedelics through decreased synaptic conductance. Multimodal studies further show that increases in electrophysiological entropy covary with alterations in large-scale functional connectivity. In this context, GFP peak density may provide a complementary, topography-based measure of enhanced dynamical richness under psilocybin. The most robust microstate-level effect was a pronounced reduction in microstate lifespan at T2 and T3 (Figure, Table), indicating reduced temporal stability and faster switching between global field configurations. This destabilization is consistent with convergent evidence reviewed in the Introduction, including reductions in alpha and beta power under psychedelicsand disruption of large-scale resting-state networks. Notably, microstate statistics are substantially determined by the spectral and autocorrelation structure of the underlying EEG, suggesting the observed lifespan reductions are tightly linked to psilocybin-induced spectral acceleration rather than representing a fully independent topographic phenomenon. The temporal profile of microstate changes closely followed psilocybin pharmacokinetics. Lifespan reductions peaked at T2 (50-60 min) and T3 (90-100 min), matching peak plasma psilocin levels and maximal subjective effects. This pattern of accelerated state transitions is compatible with previous reports linking psychedelic effects to 5-HT 2A -mediated mechanisms and large-scale network reorganization, and with recent evidence for time-dependent alterations in sensory and cognitive processing. Microstate frequency of occurrence increased at T2 and T3 (Figure, Table), complementing lifespan reductions and consistent with accelerated state transitions. Shorter microstate durations necessarily imply more frequent state onsets within a fixed recording window. Similar acceleration of brain state dynamics has been reported in analyses of metastable topographic sequences and time-resolved EEG under DMT, and spontaneous switching between functional connectivity states has been linked to cognitive performance. Together with complexity-based findings, these results support theoretical models proposing that psychedelics increase the rate at which the brain explores its available state space. In contrast to robust effects on temporal parameters, microstate coverage was largely preserved (Supplementary Figure), with only a transient significant difference at T2 in the 2-20 Hz bandwidth (p = 0.015). Coverage reflects the proportion of time spent in each microstate class and thus the relative dominance of distinct global configurations. Its broad stability-which should be interpreted cautiously given N=15, as the study may be underpowered to detect smaller coverage changes-suggests that psilocybin is primarily associated with changes in how brain states are sequenced and how long they persist, rather than which states are accessible. This dissociation is notable given that altered microstate coverage has been reported in psychiatric conditions such as schizophrenia and depression. The preservation of coverage alongside altered temporal dynamics may therefore distinguish the psychedelic state from pathological microstate profiles, consistent with the notion that psychedelics transiently increase dynamical flexibility without fundamentally disrupting the repertoire of large-scale brain states. Of potential translational interest were associations between microstate dynamics, acute subjective experience, and persisting effects, most prominently at T2. Using principal component analysis to reduce dimensionality, we observed significant correlations between microstate-derived components and experiential components at T2 (Figures). This suggests that inter-individual variation in microstate reorganization tracks variation in subjective phenomenology during the acute psychedelic state. This finding is consistent with prior reports linking electrophysiological measures to experiential intensity, including oscillatory synchronization and signal diversity metrics. Importantly, the same microstate components were also associated with Persisting Effects Questionnaire scores assessed 28 days post-administration, suggesting that the degree of acute microstate destabilization is associated with the magnitude of longer-term psychological change. Greater microstate destabilization at T2 was associated with stronger persisting positive changes in mood, well-being, and personal meaning. However, clinical validation in patient populations is required before any therapeutic inference can be drawn. These results complement emerging evidence that acute neural signatures are associated with persisting outcomes following psychedelic exposure. Within theoretical frameworks such as the entropic brain and REBUS (Relaxed Beliefs Under Psychedelics) models, transient destabilization of hierarchical constraints and increased dynamical flexibility are proposed to facilitate psychological change by enabling revision of maladaptive priors. Neuroplasticity potentially induced by psychedelics may, speculatively, provide a mechanistic substrate for such lasting effects. While our findings are consistent with this view, they remain correlational and do not establish causality. Several limitations warrant consideration. The sample size (N=15), while within the range of psychedelic neuroimaging studies, limits sensitivity to smaller effects and fine-grained analyses. Post-hoc power analysis based on observed effect sizes (α = 0.05, N = 15) indicated adequate power for condition effects in the primary 1-40 Hz analysis: lifespan (η 2 g = 0.063, power = 0.68) and frequency of occurrence (η 2 g = 0.076, power = 0.78). Power was substantially lower for time effects (power ≤ 0.36) and condition × time interactions (power ≤ 0.59) in the primary analysis, and across all factors in the secondary 2-20 Hz analysis (maximum power = 0.42), indicating limited sensitivity to smaller or time-resolved effects. An important alternative interpretation of the GFP and lifespan findings is that psilocybin increases heart rate, pupil dilation, and general sympathetic arousal at peak intoxication, which could drive faster EEG fluctuations independently of psychedelic-specific mechanisms. Physiological signals (heart rate, respiration) were not coregistered in the present study, and their potential contribution to microstate dynamics cannot be ruled out; future studies should include concurrent physiological monitoring and regression of cardiovascular artefacts. Crucially, establishing whether microstate dynamics are associated with therapeutic efficacy will require direct testing in clinical populations undergoing psychedelic-assisted treatment. More generally, integration with multimodal imaging, pharmacological manipulations, and computational modelingwill be essential for clarifying what specific aspects of neural processing microstates reflect. At present, microstate metrics appear well-suited as pragmatic biomarkers of altered brain dynamics, but their functional interpretation remains incomplete, constraining mechanistic insight into psychedelic therapy.

CONCLUSIONS

Psilocybin was associated with a marked reorganization of resting-state EEG microstate dynamics, characterized by an increased number of GFP peaks, reduced microstate lifespan, and higher frequency of occurrence, while microstate coverage was largely preserved. These findings indicate an acceleration and destabilization of large-scale brain dynamics rather than a change in the repertoire of accessible brain states. Importantly, individual differences in microstate dynamics during the acute psychedelic state were associated with both subjective experience and persisting psychological effects, suggesting that microstate dynamics show promise as neural markers of psychedelic-associated alterations in consciousness with potential therapeutic implications, although their mechanistic interpretation remains limited.

SUPPLEMENTARY INFORMATION 7. ACKNOWLEDGMENTS

The authors thank all participants who volunteered for this study, as well as the technical staff at the National Institute of Mental Health (Klecany, Czech Republic) for their assistance with data collection and EEG recordings.

FUNDING

This work was supported by grants from: • Czech Health Research Council (project NU21-04-00307 and NW24-04-00413), • Czech Science Foundation (project 23-07578K), • Ministry of the Interior of the Czech Republic (project VK01010212), • Long-term conceptual development of research organization (RVO 00023752), 1-40 Hz). Significant differences between placebo and psilocybin condition are marked with vertical lines and their respective p-values based on pairwise t-test post hoc (preceded by repeated-measures ANOVA), corrected for multiple comparisons using Benjamini-Hochberg FDR procedure. For tabulated results see Table.

Full Text PDF

Full Paper PDF

Create a free account to open full-text PDFs.

Study Details

References (34)

Papers cited by this study that are also in Blossom

Psychedelic drugs: neurobiology and potential for treatment of psychiatric disorders

Vollenweider, F. X., Preller, K. H. · Nature Reviews Neuroscience (2020)

546 cited
Psilocybin with psychological support for treatment-resistant depression: an open-label feasibility study

Carhart-Harris, R. L., Bolstridge, M., Rucker, J. et al. · Lancet Psychiatry (2016)

1520 cited
2144 cited
Effects of Psilocybin-Assisted Therapy on Major Depressive Disorder

Davis, A. K., Barrett, F. S., May, D. G. et al. · JAMA Psychiatry (2021)

1236 cited
Single-Dose Psilocybin for a Treatment-Resistant Episode of Major Depression

Goodwin, G. M., Aaronson, S. T., Alvarez, O. et al. · New England Journal of Medicine (2022)

1057 cited
The therapeutic potential of psychedelic drugs: past, present, and future

Carhart-Harris, R. L., Goodwin, G. M. · Neuropsychopharmacology (2017)

669 cited
Quality of acute psychedelic experience predicts therapeutic efficacy of psilocybin for treatment-resistant depression

Roseman, L., Nutt, D. J., Carhart-Harris, R. L. · Frontiers in Pharmacology (2018)

805 cited
Psilocybin occasioned mystical-type experiences: immediate and persisting dose-related effects

Griffiths, R. R., Johnson, M. W., Richards, W. A. et al. · Psychopharmacology (2011)

934 cited
Neural correlates of the psychedelic state as determined by fMRI studies with psilocybin

Carhart-Harris, R. L., Erritzoe, D., Williams, T. et al. · PNAS (2012)

1178 cited
Show all 34 references
The entropic brain: a theory of conscious states informed by neuroimaging research with psychedelic drugs

Carhart-Harris, R. L., Leech, R., Shanahan, M. et al. · Frontiers in Human Neuroscience (2014)

1269 cited
Increased spontaneous MEG signal diversity for psychoactive doses of ketamine, LSD and psilocybin

Schartner, M., Carhart-Harris, R. L., Barrett, A. B. et al. · Scientific Reports (2017)

444 cited
Broadband Cortical Desynchronization Underlies the Human Psychedelic State

Muthukumaraswamy, S. D., Carhart-Harris, R. L., Moran, R. J. et al. · Journal of Neuroscience (2013)

498 cited
Psilocybin desynchronizes brain networks

Nicol, G. E. · Nature (2024)

226 cited
Neural correlates of the DMT experience assessed with multivariate EEG

Timmermann, C., Roseman, L., Schartner, M. et al. · Scientific Reports (2019)

321 cited
Psilocybin can occasion mystical-type experiences having substantial and sustained personal meaning and spiritual significance

Griffiths, R. R., Richards, W. A., Mccann, U. et al. · Journal of Psychopharmacology (2006)

1671 cited
Effects of external stimulation on psychedelic state neurodynamics

Mediano, P. A. M., Rosas, F. E., Timmermann, C. et al. · ACS Chemical Neuroscience (2024)

58 cited
Neural complexity is increased after low doses of LSD, but not moderate to high doses of oral THC or methamphetamine

Murray, C., Frohlich, J, Haggarty, C. J., Tare, I. et al. · Neuropsychopharmacology (2024)

29 cited
A whole-brain model of the neural entropy increase elicited by psychedelic drugs

Herzog, R., Mediano, P. A. M., Rosas, F. E. et al. · Scientific Reports (2023)

49 cited
Human brain effects of DMT assessed via EEG-fMRI

Timmermann, C., Roseman, L., Haridas, S. et al. · PNAS (2023)

213 cited
Network control energy reductions under DMT relate to serotonin receptors, signal diversity, and subjective experience

Singleton, S. P., Timmermann, C., Luppi, A. I. et al. · Communications Biology (2025)

5 cited
Default Mode Network Modulation by Psychedelics: A Systematic Review

Gattuso, J. J., Perkins, D., Ruffell, S. G. D. et al. · International Journal of Neuropsychopharmacology (2022)

219 cited
Time-resolved Neural and Experience Dynamics of Medium- and High-dose N,N-Dimethyltryptamine

Lewis-Healey, E., Pallavicini, C., Cavanna, F. et al. · Journal of Cognitive Neuroscience (2025)

3 cited
Neural and subjective effects of inhaled N,N-dimethyltryptamine in natural settings

Pallavicini, C., Cavanna, F., Zamberlan, F. et al. · Journal of Psychopharmacology (2021)

70 cited
REBUS and the Anarchic Brain: Toward a Unified Model of the Brain Action of Psychedelics

Carhart-Harris, R. L., Friston, K. J. · Pharmacological Reviews (2019)

1128 cited
Dynamical exploration of the repertoire of brain networks at rest is modulated by psilocybin

Lord, L. D., Expert, P., Atasoy, S. et al. · NeuroImage (2019)

246 cited
163 cited
Baseline power of theta oscillations predicts mystical-type experiences induced by DMT in a natural setting

Tagliazucchi, E., Zamberlan, F., Cavanna, F. et al. · Frontiers in Psychiatry (2021)

3 cited
Brain dynamics predictive of response to psilocybin for treatment-resistant depression

Vohryzek, J., Cabral, J., Lord, L. D. et al. · Brain Communications (2024)

31 cited