pith. sign in

arxiv: 2104.07445 · v3 · submitted 2021-04-15 · 🧬 q-bio.NC · math.DS

Simulations Approaching Data: Cortical Slow Waves in Inferred Models of the Whole Hemisphere of Mouse

Cristiano Capone , Chiara De Luca , Giulia De Bonis , Robin Gutzen , Irene Bernava , Elena Pastorelli , Francesco Simula , Cosimo Lupo
show 6 more authors
Leonardo Tonielli Anna Letizia Allegra Mascaro Francesco Resta Francesco Pavone Micheal Denker Pier Stanislao Paolucci
This is my paper

Pith reviewed 2026-05-24 14:05 UTC · model grok-4.3

classification 🧬 q-bio.NC math.DS
keywords cortical slow wavesmean-field modelwide-field calcium imagingdata-driven inferencemouse cortexneuromodulationtraveling waveswhole-hemisphere dynamics
0
0 comments X

The pith

A two-loop inference method lets a mean-field model match the traveling slow waves recorded across the mouse cortex.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents a data-driven approach that infers parameters for a mean-field model of cortical activity from wide-field calcium imaging of an entire mouse hemisphere. An inner optimization loop maximizes the likelihood of the model given the data, while an outer loop tunes a periodic neuromodulation term by directly comparing observables that quantify the speed, direction, and non-stationarity of slow waves. The resulting simulations reproduce most of the nonlinear and time-varying features seen in the experimental recordings. This matters because it supplies a practical route from high-resolution imaging data to mechanistic simulations that can explore how brain states are shaped by neuromodulation.

Core claim

The two-loop inference procedure yields a mean-field model whose simulated activity reproduces the spatio-temporal features of cortical slow waves, including their non-stationary and nonlinear dynamics, as observed in whole-hemisphere wide-field calcium imaging of the mouse brain.

What carries the argument

Two-loop inference on a mean-field model: inner loop maximizes likelihood of the model parameters, outer loop optimizes periodic neuromodulation by matching wave-characterization observables.

If this is right

  • The inferred model supplies a quantitative description of how periodic neuromodulation controls the propagation of slow waves.
  • Direct comparison of observables allows systematic testing of whether a given mean-field description is sufficient at the hemispheric scale.
  • The same inference pipeline can be used to generate families of models that track transitions between different brain states.
  • Simulations built this way become testable predictions for how external perturbations should alter wave patterns.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could be applied to other large-scale recording techniques such as voltage-sensitive dyes or multi-electrode arrays.
  • If the approach generalizes, it suggests that mean-field reductions are adequate for capturing the dominant slow-wave phenomenology without full spiking-network detail.
  • One could check whether the inferred neuromodulation waveform aligns with independent measurements of neuromodulator concentrations.

Load-bearing premise

The chosen wave observables together with the mean-field approximation capture the essential mechanisms that generate the observed slow-wave dynamics without needing extra unmodeled biological factors.

What would settle it

New calcium-imaging recordings from a different animal or brain state in which the optimized model, even after re-tuning the neuromodulation parameters, fails to reproduce the measured wave speeds, directions, or non-stationary statistics.

Figures

Figures reproduced from arXiv: 2104.07445 by Anna Letizia Allegra Mascaro, Chiara De Luca, Cosimo Lupo, Cristiano Capone, Elena Pastorelli, Francesco Pavone, Francesco Resta, Francesco Simula, Giulia De Bonis, Irene Bernava, Leonardo Tonielli, Micheal Denker, Pier Stanislao Paolucci, Robin Gutzen.

Figure 1
Figure 1. Figure 1: Experimental setup and sketch of the model. A. The cortical activity of the right hemisphere of a mouse (under Ketamine/Xylazine anesthesia) is recorded through wide-field calcium imaging (GCaMP6f indicator). B. The calcium indicator response is slow if compared with the time scales of neuron activity. Applying to the original signal (upper) a proper deconvolution function fitted from single-spike response… view at source ↗
Figure 2
Figure 2. Figure 2: Summary of wave propagation properties in experimental trials for two different mice. Measures for the quantitative characterization and comparison of waves accumulated across all pixels and all trials in mouse 1 (upper row) and 2 (lower row), respectively. A. Local wave velocity distributions measured over the original dataset at a spatial resolution of 0.1 mm (black thick line) and the one downsampled at… view at source ↗
Figure 3
Figure 3. Figure 3: Inner Loop. A. The inner loop: likelihood maximization. B,C,D,E,F,G. Spatial maps of estimated parameters after 700 iterations of iRprop [35], that we further improved by setting priors for the expected spatial decay law of lateral connectivity. Parameters k0, e, λ, and a are the connectivity parameters (weight, eccentricity, spatial scale, and anisotropy), as defined in Methods; Iext and b are the local e… view at source ↗
Figure 4
Figure 4. Figure 4: Outer Loop. A. The second step of the inference method consists in comparing simulations and data, to optimize neuro-modulation amplitude and period hyper-parameters. Here, we have a bi-directional flow of information from data to model and the other way round. B. Direct comparison of the distributions of velocity (left panel), waves direction (central panel), and IWI (right panel) between data and simulat… view at source ↗
Figure 5
Figure 5. Figure 5: Dynamics comparison. Raster plot of the detected waves in experimental data (top), in simulated data without neuro-modulation (output of the inner loop, center), and in simulated data with optimal neuro-modulation (output of the outer loop, bottom). For visualization purposes, only 30 s out of the 40 s of the recording (mouse 1, trial 1) are reported. All the waves detected in the three scenarios are class… view at source ↗
Figure 6
Figure 6. Figure 6: Validation of propagation modes. Analysis over the full set of 6 experimental trials and corresponding optimal simulations (for each of the two mice). Optimal simulations are selected trial-wise following the grid search approach, as described in section The Outer Loop. A and E. Direct comparison between experimental and simulated distributions of wave velocity (left panel), direction (central panel), and … view at source ↗
read the original abstract

Thanks to novel, powerful brain activity recording techniques, we can create data-driven models from thousands of recording channels and large portions of the cortex, which can improve our understanding of brain-states neuromodulation and the related richness of traveling waves dynamics. We investigate the inference of data-driven models and the comparison among experiments and simulations, through the characterization of the spatio-temporal features of cortical waves in experimental recordings and simulations. Inference is built in two steps: the inner loop that optimizes by likelihood maximization a mean-field model, and the outer loop that optimizes a periodic neuro-modulation by relying on direct comparison of observables apt for the characterization of cortical slow waves. The model is capable to reproduce most of the features of the non-stationary and non-linear dynamics displayed by the high-resolution recording of the in-vivo mouse brain obtained by wide-field calcium imaging techniques. The proposed approach is of interest for both experimental and computational neuroscientists.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper describes a two-loop inference procedure for a mean-field model of cortical slow-wave dynamics across the mouse hemisphere. The inner loop performs likelihood maximization to fit mean-field parameters to wide-field calcium imaging data; the outer loop then tunes periodic neuromodulation parameters by direct comparison to chosen observables that characterize traveling waves. The central claim is that the resulting model reproduces most of the non-stationary and non-linear spatio-temporal features observed in the experimental recordings.

Significance. If the reproduction can be shown to be robust and not an artifact of fitting the same observables used for parameter tuning, the two-loop framework would offer a practical route for constructing data-driven mean-field models that incorporate neuromodulatory effects on wave propagation. This could be useful for both experimentalists seeking to interpret imaging data and modelers interested in linking global modulation to emergent cortical dynamics.

major comments (3)
  1. [Abstract / outer-loop Methods] Abstract and Methods (outer-loop description): the claim that the model 'reproduces most of the features' of the non-stationary/non-linear dynamics is not accompanied by quantitative metrics (e.g., error on wave-speed histograms, power spectra, or higher-order spatio-temporal correlations), cross-validation details, or an explicit count of tested observables. Without these, it is impossible to judge whether the agreement is substantive or limited to the summary statistics used in the outer loop.
  2. [Methods (outer loop)] Methods (two-loop procedure): the outer loop optimizes the periodic neuromodulation parameters by direct comparison to observables extracted from the identical experimental dataset used for the inner-loop likelihood fit. This renders the reported reproduction at least partly by construction rather than an independent test of predictive power; an explicit statement of which observables were held out, or a comparison against a null model without the outer loop, is needed to assess whether the mean-field plus global modulation captures the essential mechanisms.
  3. [Results (wave characterization)] Results (comparison of waves): the mean-field approximation averages out local heterogeneity and microscopic connectivity that can sustain or shape traveling waves. If the chosen observables omit perturbation responses or higher-order correlations, agreement on summary statistics may be superficial; the manuscript should report at least one test (e.g., response to focal perturbation or spatial correlation functions) that would falsify the mean-field description.
minor comments (2)
  1. [Methods] Notation for the mean-field variables and the periodic modulation term should be introduced with explicit equations rather than descriptive text only.
  2. [Figures] Figure captions should state the number of experimental sessions/animals and the exact observables used for outer-loop optimization.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below, agreeing where revisions are needed to clarify the strength of the reproduction claims and the validation of the two-loop procedure.

read point-by-point responses

  1. Referee: [Abstract / outer-loop Methods] Abstract and Methods (outer-loop description): the claim that the model 'reproduces most of the features' of the non-stationary/non-linear dynamics is not accompanied by quantitative metrics (e.g., error on wave-speed histograms, power spectra, or higher-order spatio-temporal correlations), cross-validation details, or an explicit count of tested observables. Without these, it is impossible to judge whether the agreement is substantive or limited to the summary statistics used in the outer loop.

    Authors: We agree that the manuscript would benefit from explicit quantitative metrics to support the reproduction claim. In the revised version we will add error measures (e.g., Wasserstein distance or MSE on wave-speed histograms and power spectra), report the exact number of observables used in the outer loop, and include a brief description of any cross-validation performed during the inner-loop likelihood maximization. revision: yes

  2. Referee: [Methods (outer loop)] Methods (two-loop procedure): the outer loop optimizes the periodic neuromodulation parameters by direct comparison to observables extracted from the identical experimental dataset used for the inner-loop likelihood fit. This renders the reported reproduction at least partly by construction rather than an independent test of predictive power; an explicit statement of which observables were held out, or a comparison against a null model without the outer loop, is needed to assess whether the mean-field plus global modulation captures the essential mechanisms.

    Authors: The referee is correct that the outer-loop observables are drawn from the same recordings. The separation between loops is that the inner loop maximizes likelihood on the raw calcium time series while the outer loop matches derived wave statistics; nevertheless this remains a form of in-sample tuning. We will revise the Methods to state this limitation explicitly and add a null-model comparison (mean-field without outer-loop neuromodulation) to quantify the added explanatory power of the periodic modulation. revision: yes

  3. Referee: [Results (wave characterization)] Results (comparison of waves): the mean-field approximation averages out local heterogeneity and microscopic connectivity that can sustain or shape traveling waves. If the chosen observables omit perturbation responses or higher-order correlations, agreement on summary statistics may be superficial; the manuscript should report at least one test (e.g., response to focal perturbation or spatial correlation functions) that would falsify the mean-field description.

    Authors: We acknowledge that the mean-field description necessarily averages local heterogeneity. The current observables target global traveling-wave statistics that the model is designed to reproduce. We will add a short discussion of this limitation and, where the existing data permit, include at least one additional diagnostic (spatial correlation functions) that could potentially falsify the mean-field approximation. revision: partial

Circularity Check

0 steps flagged

No significant circularity: explicit data-driven fitting procedure

full rationale

The paper describes an explicit two-step inference procedure (inner-loop likelihood maximization on a mean-field model; outer-loop optimization of periodic neuromodulation via direct observable comparison) applied to the same experimental dataset. The claim that the resulting model reproduces the observed features is the intended and direct outcome of this fitting process, not a first-principles derivation or independent prediction that reduces to its inputs by construction. No self-citation load-bearing steps, self-definitional equations, ansatz smuggling, or renaming of known results appear in the provided abstract or method description. This is a standard model-fitting workflow whose reproduction of fitted observables is expected and transparent.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract supplies no explicit parameter counts or invented entities; the mean-field model and optimization loops rest on standard domain assumptions whose details are not visible.

free parameters (1)
  • periodic neuromodulation parameters
    Optimized in the outer loop by direct comparison to wave observables extracted from the recording data.
axioms (1)
  • domain assumption Mean-field model is an adequate approximation for the cortical population dynamics underlying slow waves
    Invoked as the inner model whose parameters are optimized by likelihood maximization.

pith-pipeline@v0.9.0 · 5753 in / 1268 out tokens · 22618 ms · 2026-05-24T14:05:32.970289+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

78 extracted references · 78 canonical work pages · 1 internal anchor

  1. [1]

    Lin, M. Z. & Schnitzer, M. J. Genetically encoded indicators of neuronal activity. Nature neuroscience 19, 1142–1153 (2016)

    work page 2016

  2. [2]

    Sabatini, B. L. & Tian, L. Imaging neurotransmitter and neuromodulator dynamics in vivo with genet- ically encoded indicators. Neuron 108, 17–32 (2020)

    work page 2020

  3. [3]

    Mohajerani, M. H. et al. Spontaneous cortical activity alternates between motifs defined by regional axonal projections. Nature neuroscience 16, 1426 (2013)

    work page 2013

  4. [4]

    K., Dickson, C

    Greenberg, A., Abadchi, J. K., Dickson, C. T. & Mohajerani, M. H. New waves: Rhythmic electrical field stimulation systematically alters spontaneous slow dynamics across mouse neocortex. Neuroimage 174, 328–339 (2018)

    work page 2018

  5. [5]

    Akemann, W. et al. Imaging neural circuit dynamics with a voltage-sensitive fluorescent protein.Journal of neurophysiology 108, 2323–2337 (2012)

    work page 2012

  6. [6]

    & Konnerth, A

    Grienberger, C. & Konnerth, A. Imaging calcium in neurons. Neuron 73, 862–885 (2012)

    work page 2012

  7. [7]

    Chen, T.-W. et al. Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499, 295– 300 (2013)

    work page 2013

  8. [8]

    & Kn¨ opfel, T

    Shimaoka, D., Song, C. & Kn¨ opfel, T. State-dependent modulation of slow wave motifs towards awak- ening. Frontiers in cellular neuroscience 11, 108 (2017)

    work page 2017

  9. [9]

    Scott, G. et al. Voltage imaging of waking mouse cortex reveals emergence of critical neuronal dynamics. Journal of Neuroscience 34, 16611–16620 (2014)

    work page 2014

  10. [10]

    Fagerholm, E. D. et al. Cortical entropy, mutual information and scale-free dynamics in waking mice. Cerebral cortex 26, 3945–3952 (2016)

    work page 2016

  11. [11]

    Wright, P. W. et al. Functional connectivity structure of cortical calcium dynamics in anesthetized and awake mice. PloS one 12, e0185759 (2017). 21

    work page 2017

  12. [12]

    P., Chan, A

    Vanni, M. P., Chan, A. W., Balbi, M., Silasi, G. & Murphy, T. H. Mesoscale mapping of mouse cor- tex reveals frequency-dependent cycling between distinct macroscale functional modules. Journal of Neuroscience 37, 7513–7533 (2017)

    work page 2017

  13. [13]

    Xie, Y. et al. Resolution of high-frequency mesoscale intracortical maps using the genetically encoded glutamate sensor iGluSnFR. Journal of Neuroscience 36, 1261–1272 (2016)

    work page 2016

  14. [14]

    & Komiyama, T

    Ren, C. & Komiyama, T. Characterizing Cortex-Wide Dynamics with Wide-Field Calcium Imaging. Journal of Neuroscience 41, 4160–4168. issn: 0270-6474. eprint: https : / / www . jneurosci . org / content/41/19/4160.full.pdf. https://www.jneurosci.org/content/41/19/4160 (2021)

    work page 2021

  15. [15]

    Celotto, M. et al. Analysis and model of cortical slow waves acquired with optical techniques. Methods and protocols 3, 14 (2020)

    work page 2020

  16. [16]

    Brier, L. M. et al. Separability of calcium slow waves and functional connectivity during wake, sleep, and anesthesia. Neurophotonics 6, 035002 (2019)

    work page 2019

  17. [17]

    & Robinson, P

    Van Albada, S., Kerr, C., Chiang, A., Rennie, C. & Robinson, P. Neurophysiological changes with age probed by inverse modeling of EEG spectra. Clinical neurophysiology 121, 21–38 (2010)

    work page 2010

  18. [18]

    Jirsa, V. K. et al. The virtual epileptic patient: individualized whole-brain models of epilepsy spread. Neuroimage 145, 377–388 (2017)

    work page 2017

  19. [19]

    Karoly, P. J. et al. Seizure pathways: A model-based investigation. PLoS computational biology 14, e1006403 (2018)

    work page 2018

  20. [20]

    Aqil, M., Atasoy, S., Kringelbach, M. L. & Hindriks, R. Graph neural fields: A framework for spa- tiotemporal dynamical models on the human connectome. PLOS Computational Biology 17, 1–29. https://doi.org/10.1371/journal.pcbi.1008310 (Jan. 2021)

    work page doi:10.1371/journal.pcbi.1008310 2021

  21. [21]

    & Del Giudice, P

    Capone, C., Gigante, G. & Del Giudice, P. Spontaneous activity emerging from an inferred network model captures complex spatio-temporal dynamics of spike data. Scientific reports 8, 17056 (2018)

    work page 2018

  22. [22]

    J., Segev, R

    Schneidman, E., Berry, M. J., Segev, R. & Bialek, W. Weak pairwise correlations imply strongly corre- lated network states in a neural population. Nature 440, 1007–1012 (2006)

    work page 2006

  23. [23]

    & Del Giudice, P

    Capone, C., Filosa, C., Gigante, G., Ricci-Tersenghi, F. & Del Giudice, P. Inferring synaptic structure in presence of neural interaction time scales. PloS one 10, e0118412 (2015)

    work page 2015

  24. [24]

    P., Gr¨ un, S

    Rostami, V., Mana, P. P., Gr¨ un, S. & Helias, M. Bistability, non-ergodicity, and inhibition in pairwise maximum-entropy models. PLoS computational biology 13, e1005762 (2017)

    work page 2017

  25. [25]

    & Boucsein, C

    Schnepel, P., Kumar, A., Zohar, M., Aertsen, A. & Boucsein, C. Physiology and Impact of Horizontal Connections in Rat Neocortex. Cerebral Cortex 25, 3818–3835. issn: 1047-3211. eprint: https : / / academic.oup.com/cercor/article-pdf/25/10/3818/14101924/bhu265.pdf . https://doi.org/ 10.1093/cercor/bhu265 (Nov. 2014)

    work page doi:10.1093/cercor/bhu265 2014

  26. [26]

    Olcese, U. et al. Spike-Based Functional Connectivity in Cerebral Cortex and Hippocampus: Loss of Global Connectivity Is Coupled to Preservation of Local Connectivity During Non-REM Sleep. Journal of Neuroscience 36, 7676–7692. issn: 0270-6474. eprint: https://www.jneurosci.org/content/36/ 29/7676.full.pdf. https://www.jneurosci.org/content/36/29/7676 (2016)

    work page 2016

  27. [27]

    & Mattia, M

    Capone, C. & Mattia, M. Speed hysteresis and noise shaping of traveling fronts in neural fields: role of local circuitry and nonlocal connectivity. Scientific reports 7, 1–10 (2017)

    work page 2017

  28. [28]

    Waves, bumps, and patterns in neural field theories

    Coombes, S. Waves, bumps, and patterns in neural field theories. Biological cybernetics 93, 91–108 (2005)

    work page 2005

  29. [29]

    A., Rennie, C

    Robinson, P. A., Rennie, C. J. & Wright, J. J. Propagation and stability of waves of electrical activity in the cerebral cortex. Physical Review E 56, 826 (1997)

    work page 1997

  30. [30]

    Li, P. et al. Measuring sharp waves and oscillatory population activity with the genetically encoded calcium indicator GCaMP6f. Frontiers in cellular neuroscience 13, 274 (2019)

    work page 2019

  31. [32]

    Pedregosa, F. et al. Scikit-learn: Machine Learning in Python. Journal of Machine Learning Research 12, 2825–2830 (2011)

    work page 2011

  32. [33]

    Pazienti, A., Galluzzi, A., Dasilva, M., Sanchez-Vives, M. V. & Mattia, M. Slow waves form expanding, memory-rich mesostates steered by local excitability in fading anesthesia. Iscience 25, 103918 (2022)

    work page 2022

  33. [34]

    Gutzen, R. et al. Comparing apples to apples – Using a modular and adaptable analysis pipeline to compare slow cerebral rhythms across heterogeneous datasets 2022. https://arxiv.org/abs/2211. 08527

    work page 2022

  34. [35]

    & H¨ usken, M

    Igel, C. & H¨ usken, M. Improving the Rprop learning algorithm in Proceedings of the second international ICSC symposium on neural computation (NC 2000) 2000 (2000), 115–121

    work page 2000

  35. [36]

    & Destexhe, A

    El Boustani, S. & Destexhe, A. A master equation formalism for macroscopic modeling of asynchronous irregular activity states. Neural computation 21, 46–100 (2009)

    work page 2009

  36. [37]

    & Del Giudice, P

    Gigante, G., Mattia, M. & Del Giudice, P. Diverse population-bursting modes of adapting spiking neurons. Physical Review Letters 98, 148101 (2007)

    work page 2007

  37. [38]

    & Paolucci, P

    Capone, C., Pastorelli, E., Golosio, B. & Paolucci, P. S. Sleep-like slow oscillations improve visual clas- sification through synaptic homeostasis and memory association in a thalamo-cortical model. Scientific Reports 9, 8990 (2019)

    work page 2019

  38. [39]

    & Destexhe, A

    Di Volo, M., Romagnoni, A., Capone, C. & Destexhe, A. Biologically realistic mean-field models of conductance-based networks of spiking neurons with adaptation. Neural computation 31, 653–680 (2019)

    work page 2019

  39. [40]

    & Destexhe, A

    Capone, C., Di Volo, M., Romagnoni, A., Mattia, M. & Destexhe, A. State-dependent mean-field for- malism to model different activity states in conductance-based networks of spiking neurons. Physical Review E 100, 062413 (2019)

    work page 2019

  40. [41]

    & H¨ usken, M

    Igel, C. & H¨ usken, M. Empirical evaluation of the improved Rprop learning algorithms.Neurocomputing 50, 105–123 (2003)

    work page 2003

  41. [42]

    & Hertz, J

    Roudi, Y., Dunn, B. & Hertz, J. Multi-neuronal activity and functional connectivity in cell assemblies. Current opinion in neurobiology 32, 38–44 (2015)

    work page 2015

  42. [43]

    & Hertz, J

    Tyrcha, J., Roudi, Y., Marsili, M. & Hertz, J. The effect of nonstationarity on models inferred from neural data. Journal of Statistical Mechanics: Theory and Experiment 2013, P03005 (2013)

    work page 2013

  43. [44]

    & Ferrari, U

    Nghiem, T.-A., Telenczuk, B., Marre, O., Destexhe, A. & Ferrari, U. Maximum-entropy models reveal the excitatory and inhibitory correlation structures in cortical neuronal activity. Physical Review E 98, 012402 (2018)

    work page 2018

  44. [45]

    Pillow, J. W. et al. Spatio-temporal correlations and visual signalling in a complete neuronal population. Nature 454, 995–999 (2008)

    work page 2008

  45. [46]

    M., Meister, M

    Park, I. M., Meister, M. L., Huk, A. C. & Pillow, J. W. Encoding and decoding in parietal cortex during sensorimotor decision-making. Nature neuroscience 17, 1395–1403 (2014)

    work page 2014

  46. [47]

    Weber, A. I. & Pillow, J. W. Capturing the dynamical repertoire of single neurons with generalized linear models. Neural computation 29, 3260–3289 (2017)

    work page 2017

  47. [48]

    Capone, C. et al. Slow waves in cortical slices: how spontaneous activity is shaped by laminar structure. Cerebral cortex 29, 319–335 (2019)

    work page 2019

  48. [49]

    & Edelman, G

    Tononi, G., Sporns, O. & Edelman, G. M. Measures of degeneracy and redundancy in biological networks. Proceedings of the National Academy of Sciences 96, 3257–3262. issn: 0027-8424. eprint: https://www. pnas.org/content/96/6/3257.full.pdf. https://www.pnas.org/content/96/6/3257 (1999)

    work page 1999

  49. [50]

    Melozzi, F. et al. Individual structural features constrain the mouse functional connectome. Proceedings of the National Academy of Sciences 116, 26961–26969. issn: 0027-8424. eprint: https://www.pnas. org/content/116/52/26961.full.pdf. https://www.pnas.org/content/116/52/26961 (2019)

    work page 2019

  50. [51]

    & Rudrauf, D

    Marrelec, G., Mess´ e, A., Giron, A. & Rudrauf, D. Functional Connectivity’s Degenerate View of Brain Computation. PLoS Computational Biology 12(10): e1005031. issn: 1553-7358. https://journals. plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1005031 (2016). 23

    work page doi:10.1371/journal.pcbi.1005031 2016

  51. [52]

    Barson, D. et al. Simultaneous mesoscopic and two-photon imaging of neuronal activity in cortical circuits. Nature methods 17, 107–113 (2020)

    work page 2020

  52. [53]

    & Talukdar, P

    Saxena, A., Tripathi, A. & Talukdar, P. Improving multi-hop question answering over knowledge graphs using knowledge base embeddings in Proceedings of the 58th annual meeting of the association for com- putational linguistics (2020), 4498–4507

    work page 2020

  53. [54]

    & Komiyama, T

    Ren, C. & Komiyama, T. Wide-field calcium imaging of cortex-wide activity in awake, head-fixed mice. STAR Protocols 2, 100973. issn: 2666-1667. https://www.sciencedirect.com/science/article/ pii/S2666166721006791 (2021)

    work page 2021

  54. [55]

    Cardin, J. A. Functional flexibility in cortical circuits. Current opinion in neurobiology 58, 175–180 (2019)

    work page 2019

  55. [56]

    Makino, H. et al. Transformation of cortex-wide emergent properties during motor learning. Neuron 94, 880–890 (2017)

    work page 2017

  56. [57]

    Li, T. et al. Earthquakes Induced by Wastewater Disposal near Musreau Lake, Alberta, 2018–2020. Seismological Society of America 93, 727–738 (2022)

    work page 2018

  57. [58]

    Pastorelli, E. et al. Scaling of a Large-Scale Simulation of Synchronous Slow-Wave and Asynchronous Awake-Like Activity of a Cortical Model With Long-Range Interconnections. Frontiers in Systems Neuroscience 13, 33. issn: 1662-5137. https://www.frontiersin.org/article/10.3389/fnsys. 2019.00033 (2019)

    work page doi:10.3389/fnsys 2019

  58. [59]

    Accessed: 2021-04-14

    Slow Wave Analysis Pipeline (SWAP): Integrating multi-scale data and the output of simulations in a reproducible and adaptable pipeline https://wiki.ebrains.eu/bin/view/Collabs/slow- wave- analysis-pipeline. Accessed: 2021-04-14

    work page 2021

  59. [60]

    & Sejnowski, T

    Muller, L., Chavane, F., Reynolds, J. & Sejnowski, T. J. Cortical travelling waves: mechanisms and computational principles. Nature Reviews Neuroscience 19, 255–268 (2018)

    work page 2018

  60. [61]

    Golosio, B. et al. Thalamo-cortical spiking model of incremental learning combining perception, context and NREM-sleep. PLoS Computational Biology 17, e1009045 (2021)

    work page 2021

  61. [62]

    Oliver, P. A. et al. Clinical effectiveness of intravenous racemic ketamine infusions in a large commu- nity sample of patients with treatment-resistant depression, suicidal ideation, and generalized anxiety symptoms: a retrospective chart review. The Journal of Clinical Psychiatry 83, 42811 (2022)

    work page 2022

  62. [63]

    Tort-Colet, N., Capone, C., Sanchez-Vives, M. V. & Mattia, M. Attractor competition enriches cortical dynamics during awakening from anesthesia. Cell Reports 35, 109270 (2021)

    work page 2021

  63. [64]

    Terzano, M. et al. The cyclic alternating pattern as a physiologic component of normal NREM sleep. Sleep 8, 137–145 (1985)

    work page 1985

  64. [65]

    Akaike information criterion

    Hu, S. Akaike information criterion. Center for Research in Scientific Computation 93 (2007)

    work page 2007

  65. [66]

    & Destexhe, A

    Muller, L., Reynaud, A., Chavane, F. & Destexhe, A. The stimulus-evoked population response in visual cortex of awake monkey is a propagating wave. Nature communications 5, 1–14 (2014)

    work page 2014

  66. [67]

    R., Silberstein, R

    Burkitt, G. R., Silberstein, R. B., Cadusch, P. J. & Wood, A. W. Steady-state visual evoked potentials and travelling waves. Clinical Neurophysiology 111, 246–258 (2000)

    work page 2000

  67. [68]

    & Paolucci, P

    Muratore, P., Capone, C. & Paolucci, P. S. Target spike patterns enable efficient and biologically plau- sible learning for complex temporal tasks. PloS one 16, e0247014 (2021)

    work page 2021

  68. [69]

    J., Rajan, K., Escola, G

    DePasquale, B., Cueva, C. J., Rajan, K., Escola, G. S. & Abbott, L. full-FORCE: A target-based method for training recurrent networks. PloS one 13, e0191527 (2018)

    work page 2018

  69. [70]

    & Paolucci, P

    Capone, C., Muratore, P. & Paolucci, P. S. Error-based or target-based? A unified framework for learning in recurrent spiking networks. PLoS computational biology 18, e1010221 (2022)

    work page 2022

  70. [71]

    N., Jin, B., Chow, C

    Chong, M. N., Jin, B., Chow, C. W. & Saint, C. Recent developments in photocatalytic water treatment technology: a review. Water research 44, 2997–3027 (2010)

    work page 2010

  71. [72]

    Tuckwell, H. C. Introduction to theoretical neurobiology: volume 2, nonlinear and stochastic theories (Cambridge University Press, 1988). 24

    work page 1988

  72. [73]

    & Amit, D

    Mascaro, M. & Amit, D. J. Effective neural response function for collective population states. Network: Computation in Neural Systems 10, 351–373 (1999)

    work page 1999

  73. [74]

    Reynolds, D. A. Gaussian mixture models. Encyclopedia of biometrics 741 (2009)

    work page 2009

  74. [75]

    Virtanen, P. et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nature Methods 17, 261–272 (2020)

    work page 2020

  75. [76]

    On Wasserstein Two Sample Testing and Related Families of Nonparametric Tests

    Ramdas, A., Garcia, N. & Cuturi, M. On Wasserstein Two Sample Testing and Related Families of Nonparametric Tests 2015. https://arxiv.org/abs/1509.02237

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  76. [77]

    De Bonis, G. et al. Analysis Pipeline for Extracting Features of Cortical Slow Oscillations. Frontiers in Systems Neuroscience 13, 70. issn: 1662-5137. https://www.frontiersin.org/article/10.3389/ fnsys.2019.00070 (2019)

    work page arXiv 2019

  77. [78]

    Resta, F., Allegra Mascaro, A. L. & Pavone, F. Study of Slow Waves (SWs) propagation through wide- field calcium imaging of the right cortical hemisphere of GCaMP6f mice (v2) 2021. https : / / kg . ebrains.eu/search/instances/Dataset/28e65cf1-ce13-4c12-92dc-743b0cb66862

    work page 2021

  78. [79]

    external

    Claeskens, G. & Hjort, N. L. Model Selection and Model Averaging https://EconPapers.repec.org/ RePEc:cup:cbooks:9780521852258 (Cambridge University Press, 2008). Supporting information Variability across trials To show the variability of the observed phenomena across different trials from the same mouse, in Fig. S1 we provided the experimental distribution...

    work page 2008