Abstract
Electrophysiological signals exhibit both periodic and aperiodic properties. Periodic oscillations have been linked to numerous physiological, cognitive, behavioral and disease states. Emerging evidence demonstrates that the aperiodic component has putative physiological interpretations and that it dynamically changes with age, task demands and cognitive states. Electrophysiological neural activity is typically analyzed using canonically defined frequency bands, without consideration of the aperiodic (1/f-like) component. We show that standard analytic approaches can conflate periodic parameters (center frequency, power, bandwidth) with aperiodic ones (offset, exponent), compromising physiological interpretations. To overcome these limitations, we introduce an algorithm to parameterize neural power spectra as a combination of an aperiodic component and putative periodic oscillatory peaks. This algorithm requires no a priori specification of frequency bands. We validate this algorithm on simulated data, and demonstrate how it can be used in applications ranging from analyzing age-related changes in working memory to large-scale data exploration and analysis.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
All empirical behavioral and physiological data reported and analyzed in this manuscript are secondary uses of data that have previously been published and/or were accessed from openly available data repositories. A copy of the simulated data, as well as the code to regenerate it, is available in the GitHub repository (https://github.com/TomDonoghue/SimFOOOF). EEG data were analyzed from a previously described study33. Open-access MEG data were analyzed from the Human Connectome Project63,64, which is described on the project site (https://www.humanconnectome.org/), and are available through the data portal (https://db.humanconnectome.org/). LFP data were analyzed from rhesus monkeys from a previously described study56. Additional LFP data from rats were accessed from the HC-2 dataset54, which is available from the Collaborative Research in Computational Neuroscience (CRCNS) data sharing portal (https://crcns.org/).
Code availability
Custom code used in this manuscript is predominantly using the Python programming language, v.3.7. In addition, some preprocessing of MEG data was done in MATLAB (R2017a), using the Brainstorm package (https://neuroimage.usc.edu/brainstorm/). The algorithm code is openly available and released under the Apache-2.0 open-source software license. The code for the algorithm is available on GitHub (https://github.com/fooof-tools/fooof), and from PyPi (https://pypi.org/project/fooof/), and includes a dedicated documentation site (https://fooof-tools.github.io/). All of the code used for the analyses is openly available, and indexed on Github (https://github.com/fooof-tools/Paper). This includes all of the code used for the simulations (https://github.com/TomDonoghue/SimFOOOF), the EEG analyses (https://github.com/TomDonoghue/EEGFOOOF) and the MEG analyses (https://github.com/TomDonoghue/MEGFOOOF).
References
Engel, A. K., Fries, P. & Singer, W. Dynamic predictions: oscillations and synchrony in top–down processing. Nat. Rev. Neurosci. 2, 704–716 (2001).
Buzsáki, G. & Draguhn, A. Neuronal oscillations in cortical networks. Science 304, 1926–1929 (2004).
Fries, P. A mechanism for cognitive dynamics: neuronal communication through neuronal coherence. Trends Cogn. Sci. 9, 474–480 (2005).
Voytek, B. et al. Oscillatory dynamics coordinating human frontal networks in support of goal maintenance. Nat. Neurosci. 18, 1318–1324 (2015).
Kopell, N. J., Gritton, H. J., Whittington, M. A. & Kramer, M. A. Beyond the connectome: the dynome. Neuron 83, 1319–1328 (2014).
Voytek, B. & Knight, R. T. Dynamic network communication as a unifying neural basis for cognition, development, aging, and disease. Biol. Psychiatry 77, 1089–1097 (2015).
Buzsáki, G., Logothetis, N. & Singer, W. Scaling brain size, keeping timing: evolutionary preservation of brain rhythms. Neuron 80, 751–764 (2013).
He, B. J. Scale-free brain activity: past, present, and future. Trends Cogn. Sci. 18, 480–487 (2014).
Jasper, H. & Penfield, W. Electrocorticograms in man: effect of voluntary movement upon the electrical activity of the precentral gyrus. Arch. Psychiatr. Nervenkr. 183, 163–174 (1949).
Crone, N. Functional mapping of human sensorimotor cortex with electrocorticographic spectral analysis. I. Alpha and beta event-related desynchronization. Brain 121, 2271–2299 (1998).
Haegens, S., Cousijn, H., Wallis, G., Harrison, P. J. & Nobre, A. C. Inter- and intra-individual variability in alpha peak frequency. NeuroImage 92, 46–55 (2014).
Mierau, A., Klimesch, W. & Lefebvre, J. State-dependent alpha peak frequency shifts: experimental evidence, potential mechanisms and functional implications. Neuroscience 360, 146–154 (2017).
Manning, J. R., Jacobs, J., Fried, I. & Kahana, M. J. Broadband shifts in local field potential power spectra are correlated with single-neuron spiking in humans. J. Neurosci. 29, 13613–13620 (2009).
Miller, K. J. et al. Human motor cortical activity is selectively phase-entrained on underlying rhythms. PLoS Comput. Biol. 8, e1002655 (2012).
Winawer, J. et al. Asynchronous broadband signals are the principal source of the BOLD response in human visual cortex. Curr. Biol. 23, 1145–1153 (2013).
Freeman, W. J. & Zhai, J. Simulated power spectral density (PSD) of background electrocorticogram (ECoG). Cogn. Neurodyn. 3, 97–103 (2009).
Podvalny, E. et al. A unifying principle underlying the extracellular field potential spectral responses in the human cortex. J. Neurophysiol. 114, 505–519 (2015).
Voytek, B. et al. Age-related changes in 1/f neural electrophysiological noise. J. Neurosci. 35, 13257–13265 (2015).
Gao, R., Peterson, E. J. & Voytek, B. Inferring synaptic excitation/inhibition balance from field potentials. NeuroImage 158, 70–78 (2017).
Dustman, R. E., Shearer, D. E. & Emmerson, R. Y. EEG and event-related potentials in normal aging. Prog. Neurobiol. 41, 369–401 (1993).
Cole, S. & Voytek, B. Cycle-by-cycle analysis of neural oscillations. J. Neurophysiol. 122, 849–861 (2019).
Lansbergen, M. M., Arns, M., van Dongen-Boomsma, M., Spronk, D. & Buitelaar, J. K. The increase in theta/beta ratio on resting-state EEG in boys with attention-deficit/hyperactivity disorder is mediated by slow alpha peak frequency. Prog. Neuropsychopharmacol. Biol. Psychiatry 35, 47–52 (2011).
Bullock, T. H., Mcclune, M. C. & Enright, J. T. Are the electroencephalograms mainly rhythmic? Assessment of periodicity in wide-band time series. Neuroscience 121, 233–252 (2003).
Miller, K. J., Sorensen, L. B., Ojemann, J. G. & den Nijs, M. Power-law scaling in the brain surface electric potential. PLoS Comput. Biol. 5, e1000609 (2009).
Groppe, D. M. et al. Dominant frequencies of resting human brain activity as measured by the electrocorticogram. NeuroImage 79, 223–233 (2013).
Buzsáki, G., Anastassiou, C. A. & Koch, C. The origin of extracellular fields and currents—EEG, ECoG, LFP and spikes. Nat. Rev. Neurosci. 13, 407–420 (2012).
He, B. J., Zempel, J. M., Snyder, A. Z. & Raichle, M. E. The temporal structures and functional significance of scale-free brain activity. Neuron 66, 353–369 (2010).
He, W. et al. Co-increasing neuronal noise and beta power in the developing brain. Preprint at bioRxiv https://doi.org/10.1101/839258 (2019).
Robertson, M. M. et al. EEG power spectral slope differs by ADHD status and stimulant medication exposure in early childhood. J. Neurophysiol. 122, 2427–2437 (2019).
Molina, J. L. et al. Memantine effects on electroencephalographic measures of putative excitatory/inhibitory balance in schizophrenia. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 5, 562–568 (2020).
Becker, R., Van de Ville, D. & Kleinschmidt, A. Alpha oscillations reduce temporal long-range dependence in spontaneous human brain activity. J. Neurosci. 38, 755–764 (2018).
Gao, R., van den Brink, R. L., Pfeffer, T. & Voytek, B. Neuronal timescales are functionally dynamic and shaped by cortical microarchitecture. Preprint at bioRxiv https://doi.org/10.1101/2020.05.25.115378 (2020).
Tran, T. T., Hoffner, N. C., LaHue, S. C., Tseng, L. & Voytek, B. Alpha phase dynamics predict age-related visual working memory decline. NeuroImage 143, 196–203 (2016).
Niso, G. et al. OMEGA: the Open MEG Archive. NeuroImage 124, 1182–1187 (2016).
Lim, S. & Goldman, M. S. Balanced cortical microcircuitry for maintaining information in working memory. Nat. Neurosci. 16, 1306–1314 (2013).
Demirtaş, M. et al. Hierarchical heterogeneity across human cortex shapes large-scale neural dynamics. Neuron 101, 1181–1194.e13 (2019).
Frauscher, B. et al. Atlas of the normal intracranial electroencephalogram: neurophysiological awake activity in different cortical areas. Brain 141, 1130–1144 (2018).
Mukamel, R. et al. Coupling between neuronal firing, field potentials, and fMRI in human auditory cortex. Science 309, 951–954 (2005).
Bartoli, E. et al. Functionally distinct gamma range activity revealed by stimulus tuning in human visual cortex. Curr. Biol. 29, 3345–3358.e7 (2019).
Muthukumaraswamy, S. D., Singh, K. D., Swettenham, J. B. & Jones, D. K. Visual gamma oscillations and evoked responses: variability, repeatability and structural MRI correlates. NeuroImage 49, 3349–3357 (2010).
Canolty, R. T. & Knight, R. T. The functional role of cross-frequency coupling. Trends Cogn. Sci. 14, 506–515 (2010).
van der Meij, R., Kahana, M. & Maris, E. Phase-amplitude coupling in human electrocorticography is spatially distributed and phase diverse. J. Neurosci. 32, 111–123 (2012).
Muthukumaraswamy, S. D. & Liley, D. T. J. 1/f electrophysiological spectra in resting and drug-induced states can be explained by the dynamics of multiple oscillatory relaxation processes. NeuroImage 179, 582–595 (2018).
Tran, T. T., Rolle, C. E., Gazzaley, A. & Voytek, B. Linked sources of neural noise contribute to age-related cognitive decline. J. Cogn. Neurosci. 32, 1813–1822 (2020).
Reimann, M. W. et al. A biophysically detailed model of neocortical local field potentials predicts the critical role of active membrane currents. Neuron 79, 375–390 (2013).
Pascual-Marqui, R. D., Valdés-Sosa, P. A. & Alvarez-Amador, A. A parametric model for multichannel EEG spectra. Int. J. Neurosci. 40, 89–99 (1988).
Hughes, A. M., Whitten, T. A., Caplan, J. B. & Dickson, C. T. BOSC: a better oscillation detection method, extracts both sustained and transient rhythms from rat hippocampal recordings. Hippocampus 22, 1417–1428 (2012).
Wen, H. & Liu, Z. Separating fractal and oscillatory components in the power spectrum of neurophysiological signal. Brain Topogr. 29, 13–26 (2016).
Cole, S. R. & Voytek, B. Brain oscillations and the importance of waveform shape. Trends Cogn. Sci. 21, 137–149 (2017).
Caplan, J. B., Bottomley, M., Kang, P. & Dixon, R. A. Distinguishing rhythmic from non-rhythmic brain activity during rest in healthy neurocognitive aging. NeuroImage 112, 341–352 (2015).
Welch, P. The use of fast Fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. IEEE Trans. Audio Electroacoustics 15, 70–73 (1967).
Gao, R. Interpreting the electrophysiological power spectrum. J. Neurophysiol. 115, 628–630 (2016).
Cole, S., Donoghue, T., Gao, R. & Voytek, B. NeuroDSP: a package for neural digital signal processing. J. Open Source Softw. 4, 1272 (2019).
Mizuseki, K., Sirota, A., Pastalkova, E. & Buzsáki, G. Theta oscillations provide temporal windows for local circuit computation in the entorhinal-hippocampal loop. Neuron 64, 267–280 (2009).
Vogel, E. K. & Machizawa, M. G. Neural activity predicts individual differences in visual working memory capacity. Nature 428, 748–751 (2004).
Lara, A. H. & Wallis, J. D. Executive control processes underlying multi-item working memory. Nat. Neurosci. 17, 876–883 (2014).
Haegens, S., Nacher, V., Luna, R., Romo, R. & Jensen, O. Alpha-oscillations in the monkey sensorimotor network influence discrimination performance by rhythmical inhibition of neuronal spiking. Proc. Natl Acad. Sci. USA 108, 19377–19382 (2011).
Gramfort, A. et al. MNE software for processing MEG and EEG data. NeuroImage 86, 446–460 (2014).
Voytek, B. & Knight, R. T. Prefrontal cortex and basal ganglia contributions to visual working memory. Proc. Natl Acad. Sci. USA 107, 18167–18172 (2010).
Bell, A. J. & Sejnowski, T. J. An information-maximisation approach to blind separation and blind deconvolution. Neural Comput. 7, 1129–1159 (1995).
Jas, M., Engemann, D. A., Bekhti, Y., Raimondo, F. & Gramfort, A. Autoreject: automated artifact rejection for MEG and EEG data. NeuroImage 159, 417–429 (2017).
Klimesch, W. EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis. Brain Res. Rev. 29, 169–195 (1999).
Van Essen, D. C. et al. The WU-Minn human connectome project: an overview. NeuroImage 80, 62–79 (2013).
Van Essen, D. C. et al. The Human Connectome Project: a data acquisition perspective. NeuroImage 62, 2222–2231 (2012).
Gross, J. et al. Good practice for conducting and reporting MEG research. NeuroImage 65, 349–363 (2013).
Tadel, F., Baillet, S., Mosher, J. C., Pantazis, D. & Leahy, R. M. Brainstorm: a user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci. 2011, 1–13 (2011).
Nolte, G. & Curio, G. The effect of artifact rejection by signal-space projection on source localization accuracy in MEG measurements. IEEE Trans. Biomed. Eng. 46, 400–408 (1999).
Fischl, B. FreeSurfer. NeuroImage 62, 774–781 (2012).
Huang, M. X., Mosher, J. C. & Leahy, R. M. A sensor-weighted overlapping-sphere head model and exhaustive head model comparison for MEG. Phys. Med. Biol. 44, 423–440 (1999).
Baillet, S. et al. Evaluation of inverse methods and head models for EEG source localization using a human skull phantom. Phys. Med. Biol. 46, 77–96 (2001).
Fonov, V., Evans, A., McKinstry, R., Almli, C. & Collins, D. Unbiased nonlinear average age-appropriate brain templates from birth to adulthood. NeuroImage 47, S102 (2009).
Izhikevich, L., Gao, R., Peterson, E. & Voytek, B. Measuring the average power of neural oscillations. Preprint at bioRxiv https://doi.org/10.1101/441626 (2018).
Acknowledgements
We thank S.R. Cole, B. Postle, R. Hammonds, T. Tran, R. van der Meij and numerous other colleagues on GitHub, bioRxiv and Twitter who contributed to usability testing and provided invaluable comments, discussion and code contributions. M.H. is supported by National Science Foundation (NSF) Graduate Research Fellowship grant no. DGE1106400. P.S. is supported by a UC San Diego Frontiers of Innovation Scholars Program fellowship. R.G. is supported by the Natural Sciences and Engineering Research Council of Canada grant no. NSERC PGS-D, the UC San Diego Kavli Innovative Research Grant, the Frontiers for Innovation Scholars Program fellowship and a Katzin Prize. J.D.W. is supported by National Institute of Mental Health (NIMH) grant no. R01-MH121448 and NIMH grant no. R01-MH117763. R.T.K. is supported by a National Institute of Neurological Disorders and Stroke (NINDS) grant no. R37NS21135, NIMH Conte Center grant no. P50MH109429 and U19 Brain Initiative grant no. U19NS107609. A.S. is supported by NIMH grant no. F32MH75317. B.V. is supported by Sloan Research Fellowship grant no. FG-2015-66057, the Whitehall Foundation grant no. 2017-12-73, the National Science Foundation grant no. BCS-1736028 and the National Institute of General Medical Sciences grant no. R01GM134363-01.
Author information
Authors and Affiliations
Contributions
T.D., M.H., A.S., R.T.K. and B.V. initiated and designed the study. M.H. and B.V. conceived of and coded the initial algorithm. T.D. and E.J.P. designed the final algorithm and coding architecture. T.D., M.H., E.J.P., A.S. and B.V. tested and refined the algorithm. T.D. implemented the final version, wrote the tutorials and maintains the GitHub package. T.D., E.J.P., P.V., R.G. and T.N. contributed code and tested the algorithm against real and simulated data. T.D. and P.S. analyzed the EEG and MEG data. B.V. and A.S. analyzed the LFP data. A.H.L. and J.D.W. contributed LFP data and insights into analyzing them. B.V. and A.S. provided overall supervision and guidance throughout the project. All authors contributed to the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Neuroscience thanks Peter Lakatos and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 False oscillatory power changes and illusory oscillations.
a, Here we took a real neural PSD (blue) and artificially introduced a change in the aperiodic exponent, similar to what is seen in healthy aging18. This PSD was then inverted back to the time domain (right panels). The exponent change manifests as amplitude differences in the time domain. This affects apparent narrowband power when an a priori filter is applied. This is despite the fact that the true oscillatory power relative to the aperiodic component is unaffected. b, Even when no oscillation is present, such as the case with the white and pink (1/f) noise examples here (blue and green, respectively), narrowband filtering gives rise to illusory oscillations where no periodic feature exists in the actual signal, by definition.
Extended Data Fig. 2 Algorithm performance on simulated data across a broader frequency range.
a-c, Power spectra were simulated across the frequency range (1–100 Hz), with two peaks, one in a low range, and one in a high range (see Methods), across five distinct noise levels (1000 spectra per noise level). a, Example power spectra with simulation parameters as aperiodic [offset, knee, exponent] and periodic [center frequency, power, bandwidth]. b, Absolute error of algorithmically identified peak center frequency, separated for the low (3–35 Hz) and high range (50–90 Hz) peaks. c, Absolute error of algorithmically identified aperiodic parameters, offset, knee, and exponent. All violin plots show full distributions, where small white dots represent median values and small box plots show median, first and third quartiles, and ranges. Note that the error axis is log-scaled in b and c.
Extended Data Fig. 3 Algorithm performance on simulated data that violate model assumptions.
a-c, Power spectra were simulated across a broader frequency range (1–100 Hz), with two peaks, one in a low range, and one in a high range (see Methods), across five distinct knees values (1000 spectra per knee value), with a fixed noise level (0.01). Power spectra were parameterized in the ‘fixed’ aperiodic mode (without a knee) to evaluate how sensitive performance is to aperiodic mode. a, Example power spectra with simulation parameters as aperiodic [offset, knee, exponent] and periodic [center frequency, power, bandwidth], showing spectra with knee values of 0 and 150, both fit in the ‘fixed’ aperiodic mode. b, Absolute error of algorithmically identified aperiodic exponent, across spectra with different knee values. Notably, exponent reconstruction is high when spectra with knees are fit without a knee parameter. c, The number of peaks fit by the model, across knee values. Note that all spectra in this group have two peaks, indicating here that the presence of knee’s in ‘fixed’ mode leads to overfitting peaks. d-f, A distinct set of simulations were created in which power spectra were created with asymmetric or skewed peaks (see Methods), across five distinct skew levels (1000 spectra per skew level). d, Example simulated spectra, showing two different skew levels. (e) Absolute error of algorithmically identified peak center frequency, across peak skewness values. f, The number of peaks fit by the model, across peak skewness. Note that all spectra in this set have one peak. g-i, A distinct set of simulations, in which time series were generated with asymmetric oscillations in the time domain, from which power spectra were calculated (see Methods), across five distinct levels of oscillation asymmetry (1000 spectra per asymmetry value). g, Example simulation of an asymmetric oscillation, simulated in the time domain, and the associated power spectrum. Note that the power spectrum displays harmonic peaks. h, Absolute error of algorithmically identified peak center frequency, across oscillation asymmetry values. i, The number of peaks fit by the model, compared across oscillation asymmetry values. Note that these simulations all contained one oscillation in the time domain. All violin plots show full distributions, where small white dots represent median values and small box plots show median, first and third quartiles, and ranges. Note that the error axis is log-scaled in b, e and h.
Extended Data Fig. 4 Oscillation band occurrences and relative powers in the MEG dataset.
a, The proportion of participants for whom an oscillation peak was fit, at each vertex, per band. b, The group level relative power, per band. For each participant, the oscillation power within the band was normalized between 0 and 1, and then averaged across all participants, such that a maximal relative power of 1 would indicate that all participants have the same location of maximal band-specific power. Note that alpha and beta have maximal values approaching 1, reflecting a high level of consistency in location of maximal power, whereas in theta the values are lower, reflecting more variability. The ‘oscillation score’ metric, as presented in Fig. 7, is the result of multiplying the occurrence probability map with the power maps.
Extended Data Fig. 5 Methods comparison of measures of the aperiodic exponent.
a-c, Comparisons of spectral parameterization to a linear fit (in log-log space), as used in BOSC47, on simulated power spectra (1000 per comparison). a, Comparison of a linear fit and spectral parameterization on a low frequency range (2–40 Hz) with one peak. b, Comparison across the same range with multiple (3) peaks. c, Comparison across a broader frequency range (1–150 Hz) with two peaks and an aperiodic knee. In all cases (a–c), spectral parameterization outperforms the linear fit. d-f, Comparisons of spectral parameterization to IRASA48. Groups of simulations mirror those used in (a–c). Note that for these simulations, the data were simulated as time series (1000 per comparison) (see Methods). IRASA and spectral parameterization are comparable for the one peak cases (d), but spectral parameterization is significantly better in the other cases (e,f). Note that as IRASA has both greater absolute error and a systematic estimation bias (see Supplementary Modeling Note). g-i, Example of IRASA and spectral parameterization applied to real data. (g) The IRASA-decomposed aperiodic component, using default settings, in orange, is compared to the original spectrum, in blue. There are still visible non-aperiodic peaks, meaning IRASA did not fully separate out the periodic and aperiodic components. h, The IRASA-decomposed aperiodic component, with increased resampling. This helps remove the peaks, but also increasingly distorts the aperiodic component, especially at the higher frequencies, due to its multi-fractal properties (the presence of a knee). i, The isolated aperiodic component from spectral parameterization (computed as the peak-removed spectrum), showing parameterization can account for concomitant large peaks and knees, providing a better fit to the data. All violin plots show full distributions, where small white dots represent median values and small box plots show median, first and third quartiles, and ranges. Note that the error axis is log scaled in a–f. * indicates a significant difference in between the distributions of errors between methods (paired samples t-tests). Discussion of these methods and results are reported in the Supplementary Modeling Note.
Supplementary information
Supplementary Information
Supplementary Modeling Note and Supplementary Tables 1 and 2.
Rights and permissions
About this article
Cite this article
Donoghue, T., Haller, M., Peterson, E.J. et al. Parameterizing neural power spectra into periodic and aperiodic components. Nat Neurosci 23, 1655–1665 (2020). https://doi.org/10.1038/s41593-020-00744-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41593-020-00744-x
This article is cited by
-
Impaired long-range excitatory time scale predicts abnormal neural oscillations and cognitive deficits in Alzheimer’s disease
Alzheimer's Research & Therapy (2024)
-
A developmental increase of inhibition promotes the emergence of hippocampal ripples
Nature Communications (2024)
-
Multisensory flicker modulates widespread brain networks and reduces interictal epileptiform discharges
Nature Communications (2024)
-
Rhythmicity of neuronal oscillations delineates their cortical and spectral architecture
Communications Biology (2024)
-
Gamma oscillatory complexity conveys behavioral information in hippocampal networks
Nature Communications (2024)