Abstract
Neurons in the brain are submerged into oscillating extracellular potential produced by synchronized synaptic currents. The dynamics of these oscillations is one of the principal characteristics of neurophysiological activity, broadly studied in basic neuroscience and used in applications. However, our interpretation of the brain waves' structure and hence our understanding of their functions depend on the mathematical and computational approaches used for data analysis. The oscillatory nature of the wave dynamics favors Fourier methods, which have dominated the field for several decades and currently constitute the only systematic approach to brain rhythms. In the following study, we outline an alternative framework for analyzing waves of local field potentials (LFPs) and discuss a set of new structures that it uncovers: a discrete set of frequency-modulated oscillatory processes—the brain wave oscillons and their transient spectral dynamics.
1. Introduction
1.1. Motivation
Brain waves are manifestations of synchronized neuronal currents widely used for describing neurophysiological activity (Fries, 2005; Buzsáki, 2011; Thut et al., 2012; Cannon et al., 2014). However, our understanding of these phenomena depends fundamentally on mathematical and computational tools used for analyzing the recorded Local Field Potentials (LFPs). Most computational methods are based on breaking the signal into a combination of basic components suggested by the study specifics, e.g., wavelet analysis is most appropriate for studying time-localized events, such as ripples or spindles (Battaglia et al., 2004; Bosnyakova et al., 2006; Sitnikova et al., 2009; Luijtelaar et al., 2011), whereas Fourier decomposition is used for describing the oscillatory patterns of LFPs (Roopun et al., 2008; Aru et al., 2015; Lozano-Soldevilla et al., 2016; Cole and Voytek, 2017). Since most techniques are backed up by a completeness theorem, it may appear that selecting a specific decomposition is only a matter of convenience. This, however, is not the case: given that physiological mechanisms of the LFP oscillations and their functions are not yet fully understood, the task of establishing a physically adequate description of the signal's structure is not idle (Kopell et al., 2010; Buzsáki at al., 2012; Sreenivasan and D'Esposito, 2019). One may draw here a historical parallel with the use of the Ptolemaic system, in which every movement of a celestial object could be decomposed into a sufficient system of epicycles (Hanson, 1960; Van der Waerden, 1974, 1982; Babb, 1977). However, it was the discovery of the heliocentric system that eventually revealed the physical laws governing planetary motion (Gallavotti, 2001).
1.2. Approach
Discrete Fourier Transform (DFT) converts data series into a superposition of discrete harmonics with fixed frequencies, proportional to a certain base frequency ω0 (Brigham, 1988). This built-in rigidity of the Fourier spectra leads to a well-known conflict between the temporal and the frequency resolutions, manifested in many fields, from biology to Quantum Mechanics, which limits the method's resolution (Folland and Sitaram, 1997; Grünbaum, 2003). In the following, we use an alternative technique—Discrete Padé Transform (DPT)1, which also converts data points into a superposition of harmonics. However, the DPT harmonics are free to change frequencies independently, adapting their values on a moment to moment basis through the Padé Approximation Theory algorithms (Baker and Graves-Morris, 1996; Bessis, 1996; Bessis and Perotti, 2009; Perotti et al., 2013, 2019; DeVito and Dabaghian, 2014; Perotti and Wojtylak, 2018).
The spectrograms of the LFPs recorded in the CA1 area of the rat's hippocampus, built using a “sliding window” version of DPT (refer to section 4), reveal patterns that open a novel perspective on the analyses of extracellular field dynamics. First, there appear to be two types of reconstructed frequencies. The first kind changes regularly across time, leaving distinct traces in the spectrogram—the spectral waves (Figure 1A). The frequencies of the second kind assume sporadic values from moment to moment and correspond to instantaneous “irregular” harmonics with much lower amplitudes. The nature of these two classes of harmonics can be explained based on several subtle theorems of Complex Analysis (Steinhaus, 1929; Froissart, 1973; Gilewicz and Pindor, 1997; Gilewicz and Kryakin, 2003). In essence, it turns out that the irregular harmonics represent the signal's noise component, ξ(t), whereas the regular, stable frequencies define its genuine oscillatory part, r(t) (Bessis, 1996; Bessis and Perotti, 2009; Perotti et al., 2013, 2019; DeVito and Dabaghian, 2014; Perotti and Wojtylak, 2018). Interestingly, the superposition of the regular harmonics, which typically constitute only 1 − 5% of the full set, captures the shape of the original signal with over 90 − 95% precision. Correspondingly, the contribution of the remaining 95 − 99% harmonics is small, typically less than 5 − 10% of the signal's amplitude (Figure 1B). Thus, according to the DPT, the brain waves consist of a few phase-modulated waves embedded into a weak noise background.
The individual oscillatory terms in (1), , are referred below as brain wave oscillons (DeVito and Dabaghian, 2014; Perotti et al., 2019).
Figure 1
Since the decomposition (1) emerges through empirical analyses, with no a priori assumptions or ansatzs, the oscillons may capture the physical organization of synchronized neuronal activity and help link empirical observations to theoretical models (Berger, 1933; Hoppensteadt and Izhikevich, 1997; Boashash, 2003; Vugt et al., 2007; Colgin, 2016).
Second, higher temporal resolutions reveal a quasiperiodic pattern of the reconstructed frequencies,
where ωq,0 is the mean frequency, and ωq,i are the magnitudes of the embedded undulations with frequencies Ωq,i and phases φq,i (Figure 1C). Importantly, the mean frequencies of the spectral waves dovetail with the mean frequencies of the traditional (i.e., Fourier-defined) rhythms. For example, the mean frequency of the lowest spectral wave (about 8 Hz) matches the mean θ-frequency and the mean frequency of the next spectral wave (about 32 Hz) aligns with the characteristic slow-γ frequency. Furthermore, the spectral undulation magnitudes are consistent with the widths of the corresponding Fourier bands (Senior et al., 2008; Carr et al., 2012; Colgin, 2015), which allows using the standard nomenclature, e.g., ωθ(t) for the spectral θ-wave, ωγ1(t) for the spectral slow-γ wave, and to write the oscillon decomposition (1) in the form
The analysis of the spectral waves carried in Perotti et al. (2019) was motivated by the assumption that, in a given physiological state, the magnitudes ωq,i and the embedded frequencies Ωq,i in the expansion (2) are relatively stable and extractable through Fourier-based analyses, such as Welch's transform (Welch, 1967; Proakis and Manolakis, 1996). Indeed, the power profiles of approximately 1 s long segments of spectral waves exhibit consistent series of isolated peaks, suggesting that the hippocampal oscillons are driven by a discrete and comparatively scarce set of spectral harmonics (Figure 1D). However, further analyses revealed that the spectral dynamics are substantially more complex, as discussed below.
2. Results
Since each instantaneous set of DPT frequencies is computed independently based on a finite number of data points, the resulting frequency patterns exhibit gaps and irregularities (Figure 2A). To capture the underlying continuous spectral dynamics (2), we reconstructed the contiguous pattern of frequencies and amplitudes by interpolating2 the “raw,” intermittent point traces over the full set of sampled times (Figure 2B). We then used the Welch's method (Welch, 1967; Proakis and Manolakis, 1996) to estimate the spectral density of the lowest spectral wave, ωθ(t). Specifically, about 12 s long LFP trace was split into Δt ≈ 2.5 s long, highly overlapping segments, , centered at times 𝔗 ≡ [t1, t2, …, tn],
with ti+1 − ti ≈ 1 ms or less, and then Welch's procedure was applied to each segment. Arranging the resulting power profiles along the discrete time axis yields a three-dimensional (3D) spectrogram shown in Figure 2C, which demonstrates several curious features.
Figure 2
First, the lateral sections of the spectrogram—the instantaneous power profiles—exhibit a series of peaks, commonly situated within discrete frequency ranges (at ΔΩθ1 ≈ 4 − 7 Hz, then at ΔΩθ2 ≈ 10 − 14 Hz, then at ΔΩθ3 ≈ 16 − 19 Hz, etc.), separated by “valleys” in which peaks are rare (first extending over Hz, second over Hz, etc). This pattern was previously observed through static spectral power profiles such as the one shown in Figure 1D or in Perotti et al. (2019). The second surprising feature of the spectrogram is that most peaks are localized not only in frequency but also in time: a typical peak grows and abates over a few 100 ms periods. Third, many peaks are recurrent, appearing and disappearing repeatedly at about the same frequency Ωθ,i. Overall, the pattern illustrated in Figure 2C suggests that the oscillons' spectra are perturbed by a series of pulses that sporadically activate and wear off, as the animal navigates.
Simulated Data
To validate the qualitative conclusion drawn from Figure 2C, we simulated a superposition of two oscillons with the spectral wave parameters derived from the recorded data. For example, the waves illustrated in Figures 2A,B were generated for the mean θ and γ frequencies, ωθ,0 ≈ 8 Hz and ωγ,0 ≈ 32 Hz, along with various specific sets of the reconstructed embedded frequencies, e.g., Ωθ,* ≈ {1.9, 4.4, 6.2, 8.2, …} Hz and Ωγ,* = {1.4, 3.9, 6.9, 9.2, …} Hz. These values were used to build “synthetic” θ and γ oscillons with spectral frequencies
which were then processed using DPT algorithms.
As expected, the stable frequencies reconstructed through DPT procedures produce clear point traces across the spectrogram, and the reconstructed contiguous spectral waves match the input data (Figures 3A,B). Correspondingly, the peaks representing the embedded frequencies appear in the Welch's spectrogram in correct positions and remain steady, nearly unchanged over the entire duration of the signal (Figure 3C). Furthermore, numerous computational experiments with synthetic oscillons produced no spurious peaks or other artifacts suggestive of the patterns visible in Figure 2C.
Figure 3
The ostensible difference between the spectrogram produced by the simulated oscillons with constant spectral waves (Figure 3C) and the ones reconstructed from the recorded LFP data (Figure 2C) suggests that the hippocampal extracellular field dynamics may not be described by quasiperiodic series (4a) with steady coefficients. The time-localized peaks visible in Figure 2C suggest that the hippocampal frequency spectra are disturbed by rapid, transient processes that appear for a short time and rapidly disappear. To verify this possibility, we applied the DPT analyses to a numerically generated signal in which the spectral waves with constant coefficients (4) were replaced by a superposition of harmonics with time-localized spectral magnitudes,
where are narrow ( ms) Gaussian pulses localized at a few discrete moments (Figure 4A). These “spectral kicks” are clearly manifested on Welch's spectrogram computed directly for the simulated spectral θ- and slow-γ waves (Figure 4B), but they do not significantly alter the reconstructed stable frequency traces (Figure 4C). Applying DPT analyses to the corresponding oscillons produces Welch's spectrograms that bear an uncanny resemblance to the spectrograms obtained for hippocampal LFP spectral dynamics (Figure 4D). The latter result suggests that the hippocampal oscillons may exhibit elaborate behaviors that include rapid, nonstationary spectral modulations that may be due to the extracellular field's endogenous dynamics or to inputs from parahippocampal or cortical networks.
Figure 4
3. Discussion
Discrete Fourier Transform techniques currently provide the most commonly used semantics and the main framework for interpreting the structure and physiological functions of the brain waves (Roopun et al., 2008; Buzsáki, 2011; Colgin, 2016). DPT offers an alternative, high-resolution technique that leads to a novel perspective on the LFP's oscillatory component, extracted from its “noise shell.” Specifically, DPT analyses indicate that the conventional, i.e., Fourier-defined θ, γ, and other brain waves conceal elaborate, frequency-modulated oscillatory processes—the oscillons, that may reflect physical dynamics of the extracellular fields.
The term “oscillons” is currently used in several fields, to designate, e.g., quasi-stable solutions of dynamic equations in field theory and cosmology (Gleiser, 1994; Copeland et al., 1995; Kasuya et al., 2003; Amin and Shirokoff, 2010) (also refer to Bogolubsky and Makhankov, 1976) or quasi-stationary undulations in granular media (Umbanhowar et al., 1996; Cerda et al., 1997). In this context, the physical origins of the brain wave oscillons require additional studies. Some properties of the oscillons' dynamics dovetail with predictions of theoretical models that aim to explain the coherent dynamics of extracellular fields through synchronization of neuronal activity in excitatory and inhibitory networks (Hoppensteadt and Izhikevich, 1998; Izhikevich, 1999a,b, 2000; Neda et al., 2000a,b). For example, the Kuramoto model of emergent synchronization (Strogatz, 2000; Arenas et al., 2008) describes networks of weakly interacting phasors with close natural frequencies ωi,
As the coupling strengths λij between the oscillators increase, the network transitions from a disordered to a partially synchronous and then to a globally synchronized state with a net phase
The form of the Equations (6) suggests that the expansion (2) should provide a natural ansatz for describing the functional form of the synchronized phase (7). Correspondingly, the initial analyses of oscillons (Perotti et al., 2019) were carried out under the assumption that spectral waves behave as almost-periodic functions with slowly varying coefficients, given that gradual changes of θ and γ bandwidths and their means, coupled to the animal's speed and acceleration, are well documented (Richard et al., 2013; Lu et al., 2020; Kropff et al., 2021). However, the current study suggests that oscillon dynamics involve not only slow but also rapid changes. In particular, it turns out that rapid dynamics affect not only the bandwidths and mean frequencies, but also the embedded frequencies, yielding time-localized “spectral pulses” that may reflect external inputs into the hippocampal CA1 area from other brain parts, e.g., from the hippocampal CA3 area or the medial entorhinal cortex (Brun et al., 2002; Kesner, 2007; Langston et al., 2010; Yamamoto and Tonegawa, 2017).
4. Methods
4.1. Discrete Fourier and Padé Transforms
Discrete Fourier Transform is produced by convolving the data values, s1, s2, …, sN, with a discrete set of harmonics with fixed frequencies,
arranged uniformly over the unit circle in the complex plane. The closer is the discrete frequency ωl = 2πl/N to the frequency of the signal's constituent waves, , the bigger is the contribution of the corresponding harmonic into the decomposition (Brigham, 1988). If the data are sampled from a combination of harmonic oscillations and a noise background,
then each frequency ωp produces a Fourier-peak, broadened and lowered by the noise ξ(t) (Newland, 2005; Perotti et al., 2014).
The DPT extends the expansion (8) from the unit circle into the complex plane, ,
where z is a generic complex number. For the oscillatory component of s(t), the sum (9), extended to infinity, yields a meromorphic function,
whose poles, , and residues, , define the frequencies, the amplitudes, and the phases of the contributing harmonics (Bessis, 1996; Bessis and Perotti, 2009; Perotti et al., 2013, 2019; DeVito and Dabaghian, 2014; Perotti and Wojtylak, 2018).
The sub-diagonal Padé approximant to (10),
rapidly approaches R(z) as the degree N of the polynomials PN(z) and QN+1(z) grows, (Baker and Graves-Morris, 1996). In particular, the poles zp of R(z) are approximated by the roots ζq of the denominator in (11), QN+1(ζq) = 0 (Bessis, 1996; Bessis and Perotti, 2009; Perotti et al., 2013).
As for the z-transform of the noise component,
the Steinhaus theorem establishes that its poles appear at the unit circle with unit probability (Steinhaus, 1929). The manifestation of this effect in the Padé approximations to Ξ(z) is subtle: the “noisy” poles are the ones that not only cluster around the unit circle, but also pair with the zeroes of ΞN(z), thus forming the so-called “Froissart doublets” (Froissart, 1973; Gilewicz and Pindor, 1997; Gilewicz and Kryakin, 2003). A typical pole-zero distance in these pairs is smaller than 10−6 − 10−7 in the standard Euclidean metric in C1. Furthermore, the Froissart doublets are unstable with respect to variations of the algorithm's parameters, in contrast with the unpaired, stable poles produced by the regular part of the signal (Froissart, 1973; Bessis, 1996; Bessis and Perotti, 2009). These qualitative differences allow the separation of the regular component of the signal from its noise background, as expressed by the decomposition (1). The original study of Steinhaus (1929) presumed uniformly distributed noise series; subsequent works cited above allow generic, continuous noise distributions.
4.2. Sliding Window
Sliding window or the Short Time Padé Transform (STPT) uses a segment of the signal of length TW, centered at time ti, to extract the time-localized spectra—in full analogy with the Short Time Fourier Transform, STFT (Howell, 2001; Jacobsen and Lyons, 2003). Plotting the reconstructed frequencies along the vertical axis and arranging the times ti horizontally yields the Padé spectrogram, which we use to illustrate spectral dynamics, in direct analogy with the standard Fourier spectrograms.
4.3. Signal Processing
The mean amplitude of the input data was normalized to . The LFPs were originally recorded at the rate Sr = 8 kHz. To increase time resolution in the biologically relevant range of frequencies (f <300 Hz), we interpolated the signal to higher rates ( kHz, kHz or kHz), which did not alter the shape of the studied spectral patterns but significantly improved stability and sharpness of the results. The oversampled time series were then downsampled 2 ≤ m ≤ 4 times, which produced m interlaced subseries that were independently studied with DPT. As one would anticipate, the stable frequencies generated by each subsequence form tight clusters of m points, grouping around the frequency produced by the original sequence, while the Froissart doublets exhibit erratic behavior (Bessis, 1996; Bessis and Perotti, 2009; Perotti et al., 2013, 2019; DeVito and Dabaghian, 2014). These procedures allow using time windows as short as TW = 10 − −20 ms while keeping the order of the Padé approximants high, N = 100 or more. Shifting the time windows by a single data point ensures maximal contiguity of the reconstructed spectral waves and the oscillons' amplitudes. The Froissart distance used to identify close pole-zero pairs (Froissart doublet) is . To increase stability, the signals were filtered between 1 and 40 Hz
Funding
MZ and YD are supported by NIH grants R01NS110806 and R01AG074226. CD and DJ are supported by NIH grants R01MH106552 and R01MH112523.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The data will be shared for research purposes once the corresponding papers have been accepted for publication.
Author contributions
YD conceived of the study, developed the method, carried analyses, and wrote the manuscript. MZ developed the method, conducted analyses, and provided figures. CD and DJ provided data and helped conceptualizing the results. LP developed the mathematical foundations of the method and helped adopting it to data analyses. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
1.^Throughout the text, terminological definitions are given in italics.
2.^The MathWorks, I. (2019). Symbolic Math Toolbox. Natick, MA, USA, see https://www.mathworks.com/help/curvefit/.
References
1
AminM.ShirokoffD. (2010). Flat-top oscillons in an expanding universe. Phys. Rev. D81, 085045. 10.1103/PhysRevD.81.085045
2
ArenasA.Díaz-GuileraA.KurthsJ.MorenoY.ZhouC. (2008). Synchronization in complex networks. Phys. Rep.469, 93–153. 10.1016/j.physrep.2008.09.002
3
AruJ.AruJ.PriesemannV.WibralM.LanaL.PipaG.et al. (2015). Untangling cross-frequency coupling in neuroscience. Curr. Opin. Neurobiol.31, 51–61. 10.1016/j.conb.2014.08.002
4
BabbS. (1977). Accuracy of planetary theories, particularly for Mars. Isis68, 426–434. 10.1086/351818
5
BakerG.Graves-MorrisP. (1996). Padé Approximants, 2nd Edn.. Cambridge, UK: Cambridge University Press.
6
BattagliaF.SutherlandG.McNaughtonB. (2004). Hippocampal sharp wave bursts coincide with neocortical “up-state” transitions. Learn. Mem.11, 697–704. 10.1101/lm.73504
7
BergerH. (1933). Über das elektrenkephalogramm des menschen. Archiv für Psychiatrie und Nervenkrankheiten99, 555–574. 10.1007/BF01797193
8
BessisD. (1996). Padé approximations in noise filtering. J. Comput. Appl. Math.66, 85–88. 10.1016/0377-0427(95)00177-8
9
BessisD.PerottiL. (2009). Universal analytic properties of noise: introducing the J-matrix formalism. J. Phys. A42, 365202. 10.1088/1751-8113/42/36/365202
10
BoashashH. (2003). Time Frequency Signal Analysis and Processing: A Comprehensive Reference, 1st Edn. Boston, MA: Elsevier.
11
BogolubskyI.MakhankovV. (1976). Lifetime of pulsating solitons in certain classical models. JETP Lett.24, 12; ibid. 25, 120 (1977); ibid. 25, 107 (1977).
12
BosnyakovaD.GabovaA.KuznetsovaG.ObukhovY.MidzyanovskayaI.SaloninD.et al. (2006). Time-frequency analysis of spike-wave discharges using a modified wavelet transform. J. Neurosci. Methods154, 80–88. 10.1016/j.jneumeth.2005.12.006
13
BrighamE. (1988). The Fast Fourier Transform And Its Applications. Englewood Cliffs, NJ: Prentice Hall.
14
BrunV.OtnassM.MoldenS.SteffenachH.WitterM.MoserM. B.et al. (2002). Place cells and place recognition maintained by direct entorhinal-hippocampal circuitry. Science296, 2243–2246. 10.1126/science.1071089
15
BuzsákiG. (2011). Rhythms in the Brain. Oxford: Oxford University Press.
16
BuzsákiG.AnastassiouC.KochC. (2012). The origin of extracellular fields and currents—EEG, ECoG, LFP and spikes. Nat. Rev. Neurosci.13, 407–420. 10.1038/nrn3241
17
CannonJ.McCarthyM.LeeS.LeeJ.B'́orgersC.WhittingtonM.et al. (2014). Neurosystems: brain rhythms and cognitive processing. Eur. J. Neurosci.39, 705–719. 10.1111/ejn.12453
18
CarrM.KarlssonM.FrankL. (2012). Transient slow gamma synchrony underlies hippocampal memory replay. Neuron75, 700–713. 10.1016/j.neuron.2012.06.014
19
CerdaE.MeloF.RicaS. (1997). Model for subharmonic waves in granular materials. Phys. Rev. Lett.79, 4570–4573. 10.1103/PhysRevLett.79.4570
20
ColeS.VoytekB. (2017). Brain oscillations and the importance of waveform shape. Trends Cogn. Sci.21, 137–149. 10.1016/j.tics.2016.12.008
21
ColginL. (2015). Do slow and fast gamma rhythms correspond to distinct functional states in the hippocampal network?Brain Res.1621, 309–315. 10.1016/j.brainres.2015.01.005
22
ColginL. (2016). Rhythms of the hippocampal network. Nat. Rev. Neurosci.17, 239–249. 10.1038/nrn.2016.21
23
CopelandE.GleiserM.MullerH. (1995). Oscillons: Resonant configurations during bubble collapse. Phys. Rev. D52, 1920–1933. 10.1103/PhysRevD.52.1920
24
DeVitoJ.DabaghianY. (2014). New Signal Processing Method Reveals Discrete Structure of local Field Potential in Hippocampus. SfN Annual Meeting 2014, Abstract 751.16/UU15.
25
FollandG.SitaramA. (1997). The uncertainty principle: a mathematical survey. J. Fourier Anal. Appl.3, 207–238. 10.1007/BF02649110
26
FriesP. (2005). A mechanism for cognitive dynamics: neuronal communication through neuronal coherence. Trends Cogn. Sci.9, 474–480. 10.1016/j.tics.2005.08.011
27
FroissartM. (1973). Approximation de Padé, application á la physique des particules élémentaires. Les rencontres physiciens-mathématiciens de Strasbourg-RCP25, 16, 1–13.
28
GallavottiG. (2001). Quasi periodic motions from Hipparchus to Kolmogorov. Atti della Accademia Nazionale dei Lincei. Classe di Scienze Fisiche, Matematiche e Naturali. Rendiconti Lincei. Matematica e Applicazioni12, 125–152. 10.48550/arXiv.chao-dyn/9907004
29
GilewiczJ.KryakinY. (2003). Froissart doublets in Padé approximation in the case of polynomial noise. J. Comput. Appl. Math.153, 235–242. 10.1016/S0377-0427(02)00674-X
30
GilewiczJ.PindorM. (1997). Padé approximants and noise: a case of geometric series. J. Comput. Appl. Math.87, 199–214. 10.1016/S0377-0427(97)00185-4
31
GleiserM. (1994). Pseudostable bubbles. Phys. Rev. D49, 2978–2981. 10.1103/PhysRevD.49.2978
32
GrünbaumF. (2003). The Heisenberg inequality for the discrete Fourier transform. Appl. Comput. Harmon. Anal.15, 163–167. 10.1016/S1063-5203(03)00033-2
33
HansonN. (1960). The mathematical power of epicyclical astronomy. Isis51, 150–158. 10.1086/348869
34
HoppensteadtF.IzhikevichE. (1997). Weakly Connected Neural Networks. New York, NY: Springer.
35
HoppensteadtF.IzhikevichE. (1998). Thalamo-cortical interactions modeled by weakly connected oscillators: could the brain use FM radio principles?Biosystems48, 85–94. 10.1016/s0303-2647(98)00053-7
36
HowellK. (2001). Principles of Fourier Analysis. Portland, OR: CRC Press.
37
IzhikevichE. (1999a). Class 1 neural excitability, conventional synapses, weakly connected networks, and mathematical foundations of pulse-coupled models. IEEE Trans. Neural Netw.10, 499–507. 10.1109/72.761707
38
IzhikevichE. (1999b). Weakly pulse-coupled oscillators, FM interactions, synchronization, and oscillatory associative memory. IEEE Trans. Neural Netw.10, 508–526. 10.1109/72.761708
39
IzhikevichE. (2000). Phase equations for relaxation oscillators. Siam. J. Appl. Math.60, 1789–1805. 10.1137/s0036139999351001
40
JacobsenE.LyonsR. (2003). The sliding DFT. Signal Process. Mag. IEEE20, 74–80. 10.1109/MSP.2003.1184347
41
KasuyaS.KawasakiM.TakahashiF. (2003). I-balls. Phys. Lett. B559, 99–106. 10.1016/S0370-2693(03)00344-7
42
KesnerR. (2007). Behavioral functions of the CA3 subregion of the hippocampus. Learn. Mem.14, 771–781. 10.1101/lm.688207
43
KopellN.KramerM.MalerbaP.WhittingtonM. (2010). Are different rhythms good for different functions?Front. Hum. Neurosci.4, 187. 10.3389/fnhum.2010.00187
44
KropffE.CarmichaelJ.MoserE.MoserM. B. (2021). Frequency of theta rhythm is controlled by acceleration, but not speed, in running rats. Neuron109, 1–11. 10.1016/j.neuron.2021.01.017
45
LangstonR.StevensonC.WilsonC.SaundersI.WoodE. (2010). The role of hippocampal subregions in memory for stimulus associations. Behav. Brain Res. 215, 275–291. 10.1016/j.bbr.2010.07.006
46
Lozano-SoldevillaD.ter HuurneN.OostenveldR. (2016). Neuronal oscillations with non-sinusoidal morphology produce spurious phase-to-amplitude coupling and directionality. Front. Comput. Neurosci.10, 87. 10.3389/fncom.2016.00087
47
LuL.RenY.YuT.LiuZ.WangS.TanL.et al. (2020). Control of locomotor speed, arousal, and hippocampal theta rhythms by the nucleus incertus. Nat. Commun.11, 262. 10.1038/s41467-019-14116-y
48
NedaZ.RavaszE.BrechetY.VicsekT.BarabasiA. L. (2000b). The sound of many hands clapping. Nature403, 849–850. 10.1038/35002660
49
NedaZ.RavaszE.VicsekT.BrechetY.BarabasiA. L. (2000a). Physics of the rhythmic applause. Phys. Rev. E61, 6987–6992. 10.1103/PhysRevE.61.6987
50
NewlandD. (2005). An introduction To Random Vibrations, Spectral And Wavelet Analysis. New York, NY: Dover Publications.
51
PerottiL.DeVitoJ.BessisD.DabaghianY. (2019). Discrete spectra of brain rhythms. Sci. Rep.9, 1105. 10.1038/s41598-018-37196-0
52
PerottiL.RegimbauT.VrinceanuD.BessisD. (2014). Identification of gravitational-wave bursts in high noise using Padé filtering. Phys. Rev. D90, 124047. 10.1103/PhysRevD.90.124047
53
PerottiL.VrinceanuD.BessisD. (2013). Enhanced frequency resolution in data analysis. Am. J. Comput. Math.3, 242–251. 10.4236/ajcm.2013.33034
54
PerottiL.WojtylakM. (2018). Matrix methods for Padé approximation: Numerical calculation of poles, zeros and residues. Lin. Alg. App.548, 95–122. 10.1016/j.laa.2018.03.004
55
ProakisJ.ManolakisD. (1996). Digital Signal Processing: Principles, Algorithms and Applications, 3rd Edn. New York, NY: Prentice-Hall.
56
RichardG.TitizA.TylerA.HolmesG.ScottR.Lenck-SantiniP. (2013). Speed modulation of hippocampal theta frequency correlates with spatial memory performance. Hippocampus23, 1269–1279. 10.1002/hipo.22164
57
RoopunA.KramerM.CarracedoL.KaiserM.DaviesC.TraubR.et al. (2008). Temporal Interactions between Cortical Rhythms. Front. Neurosci.2, 145–154. 10.3389/neuro.01.034.2008
58
SeniorT.HuxterJ.AllenK.O'NeillJ.CsicsvariJ. (2008). Gamma oscillatory firing reveals distinct populations of pyramidal cells in the CA1 region of the hippocampus. J. Neurosci.28, 2274–2286. 10.1523/JNEUROSCI.4669-07.2008
59
SitnikovaE.HramovA.KoronovskyA.van LuijtelaarG. (2009). Sleep spindles and spike—wave discharges in EEG: Their generic features, similarities and distinctions disclosed with Fourier transform and continuous wavelet analysis. J. Neurosci. Methods180, 304–316. 10.1016/j.jneumeth.2009.04.006
60
SreenivasanK.D'EspositoM. (2019). The what, where and how of delay activity. Nat. Rev. Neurosci.20, 466–481. 10.1038/s41583-019-0176-7
61
SteinhausH. (1929). Über die Wahrscheinlichkeit dafür dass der Konvergenzkreis einer Potenzreihe ihre natuerliche Grenze ist. Math. Zeitschrift, 31, 408–416.
62
StrogatzS. (2000). From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D143, 1–20. 10.1016/S0167-2789(00)00094-4
63
ThutG.MiniussiC.GrossJ. (2012). The functional importance of rhythmic activity in the brain. Curr. Biol.22, R658–R663. 10.1016/j.cub.2012.06.061
64
UmbanhowarP.MeloF.SwinneyH. (1996). Localized excitations in a vertically vibrated granular layer. Lett. Nat.382, 793–796. 10.1038/382793a0
65
Van der WaerdenB. (1974). The earliest form of the epicycle theory. J. History Astron.5, 175–185. 10.1177/002182867400500303
66
Van der WaerdenB. (1982). The motion of venus, mercury and the sun in early greek astronomy. Arch. History Exact. Sci.26, 99–113. 10.1007/BF00348348
67
van LuijtelaarG.HramovA.SitnikovaE.KoronovskiiA. (2011). Spike-wave discharges in WAG/Rij rats are preceded by delta and theta precursor activity in cortex and thalamus. Clin. Neurophysiol.122, 687–695. 10.1016/j.clinph.2010.10.038
68
van VugtM.SederbergP.KahanaM. (2007). Comparison of spectral analysis methods for characterizing brain oscillations. J. Neurosci. Methods162, 49–63. 10.1016/j.jneumeth.2006.12.004
69
WelchP. (1967). 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. 10.1109/TAU.1967.1161901
70
YamamotoJ.TonegawaS. (2017). Direct medial entorhinal cortex input to hippocampal CA1 is crucial for extended quiet awake replay. Neuron96, 217.e4–227.e4. 10.1016/j.neuron.2017.09.017
Summary
Keywords
brain rhythms, oscillons, hippocampus, theta, spectral wave
Citation
Zobaer MS, Domenico CM, Perotti L, Ji D and Dabaghian Y (2022) Rapid Spectral Dynamics in Hippocampal Oscillons. Front. Comput. Neurosci. 16:880742. doi: 10.3389/fncom.2022.880742
Received
21 February 2022
Accepted
12 April 2022
Published
10 June 2022
Volume
16 - 2022
Edited by
Daya Shankar Gupta, Husson University, United States
Reviewed by
Joao Streibel, Federal University of Santa Catarina, Brazil; Marcelo Gleiser, Dartmouth College, United States
Updates
Copyright
© 2022 Zobaer, Domenico, Perotti, Ji and Dabaghian.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuri Dabaghian yuri.a.dabaghian@uth.tmc.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.