Rapid Spectral Dynamics in Hippocampal Oscillons

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 ω₀ (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)¹, 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.

$\begin{array}{ll}s(t)={\displaystyle \sum _{q=1}^M}{A}_q{e}^{i{\Phi }_q(t)}+\xi (t).& (1)\end{array}$

The individual oscillatory terms in (1), ${\vartheta }_q(t)=A_q{e}^{i{\Phi }_q(t)}$, are referred below as brain wave oscillons (DeVito and Dabaghian, 2014; Perotti et al., 2019).

FIGURE 1

Figure 1. Spectral waves. (A) A second-long segment of the Discrete Pade' Transform (DPT) spectrogram computed for Local Field Potentials (LFPs) recorded in the hippocampal CA1 area of an actively moving rat exhibits a series of traces—the spectral waves, which can be viewed as timelines of time-dependent frequencies ω_q(t). The dot colors designate the instantaneous amplitudes A_q(t) of the corresponding oscillons (colorbar). Here, the sliding window width is T_W ≈ 50 ms, and the full number of DPT harmonics is N = 200, of which 1.7% are stable and produce spectral waves (pie diagram). (B) Oscillatory part of the LFP signal reconstructed from the stable frequencies (red) differs from the original signal (blue) by ≈ 9% of the signal's power (pie diagram). The mismatch is due to the discarded unstable frequencies, i.e., to the removed noise component ξ(t) (black curve). (C) At higher time resolutions (T_W ≈ 20 ms), spectral waves exhibit a quasiperiodic pattern. Shown are the θ (below) and the slow-γ (above) spectral waves. (D) The spectral power profiles constructed using Fourier (black line) and Welch's (red line) techniques show a set of peaks indicating the individual embedded frequencies.

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,

$\begin{array}{l}{\dot{\Phi }}_q={\omega }_q(t)={\omega }_{q,0}+{\omega }_{q,1}sin({\Omega }_{q,1}t\\ +{\phi }_{q,1})+{\omega }_{q,2}sin({\Omega }_{q,2}t+{\phi }_{q,2})+\dots ,& (2)\end{array}$

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

$\begin{array}{ll}s(t)=A_{\theta }{e}^{i{\Phi }_{\theta }(t)}+A_{{\gamma }_1}{e}^{i{\Phi }_{{\gamma }_1}(t)}+A_{{\gamma }_2}{e}^{i{\Phi }_{{\gamma }_2}(t)}+\dots +\xi (t).& (3)\end{array}$

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 interpolating² 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, ${\omega }_{\theta }^1,{\omega }_{\theta }^2,\dots ,{\omega }_{\theta }^n$, centered at times 𝔗 ≡ [t₁, t₂, …, t_n],

$\begin{array}{l}{\omega }_{\theta }^i={\omega }_{\theta }(t),\text{for }t\in \left[t_i-\Delta t/2,t_i+\Delta t/2\right],\end{array}$

with t_i+1 − t_i ≈ 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

Figure 2. Hippocampal theta wave. (A) Stable frequencies falling under 40 Hz form intermittent traces occupying the θ-domain (lower trace) and the slow-γ domain (upper trace), with the means ω_θ,0 ≈ 8 Hz and ω_γ1,0 ≈ 32 Hz (doted and dashed lines, respectively). Color of the dots represents the instantaneous amplitude of the corresponding oscillons (colorbar). (B) Interpolating the raw θ-trace (dimmed pattern in the background) over uniformly spaced time points yields the reconstructed spectral wave ω_θ(t) (solid colored line). (C) Welch's spectrogram of ω_θ(t) exhibits domains of “peak ranges” (within the domains ΔΩ_θ1 ≈ 4 − 7 Hz and ΔΩ_θ2 ≈ 10 − 14 Hz), separated by a “valley” extending over ${\overline{\Delta \Omega }}_{{\theta }_1}=7-10$ Hz.

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 ${\overline{\Delta \Omega }}_{{\theta }_1}=7--10$ Hz, second over ${\overline{\Delta \Omega }}_{{\theta }_2}=14-16$ 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

$\begin{array}{ll}{\omega }_{\theta }(t)={\omega }_{\theta ,0}+{\omega }_{\theta ,1}sin({\Omega }_{\theta ,1}t)+{\omega }_{\theta ,2}sin({\Omega }_{\theta ,2}t)+\dots ,& (\mathrm{4\theta })\end{array}$

$\begin{array}{ll}{\omega }_{\gamma }(t)={\omega }_{\gamma ,0}+{\omega }_{\gamma ,1}sin({\Omega }_{\gamma ,1}t)+{\omega }_{\gamma ,2}sin({\Omega }_{\gamma ,2}t)+\dots ,& (\mathrm{4\gamma })\end{array}$

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

Figure 3. Simulated signal with constant embedded frequencies. (A) The power spectra of the simulated LFP wave—a combination of two synthetic oscillons produced at the same sampling rate as the recorded data. The timelines of the reconstructed stable frequencies undulate around the imputed mean values, ${\tilde{\omega }}_{\theta ,0}\approx 8$ Hz (dotted line) and ${\tilde{\omega }}_{{\gamma }_1,0}\approx 32$ Hz (dashed line). (B) Interpolating the “raw” θ-trace (both frequencies and amplitudes) over the full set of timepoints yields the reconstructed spectral wave (solid colored line), which also matches the inputted ${\tilde{\omega }}_{\theta }$ (dashed line), with the embedded frequencies Ω_θ,1 ≈ 1.9 Hz, Ω_θ,2 ≈ 4.4 Hz, Ω_θ,3 ≈ 6.2 Hz, and Ω_θ,3 ≈ 8.2. (C) The “evolvent” Welch spectrogram reveals the embedded frequencies Ω_θ,1, Ω_θ,2 Hz, and Ω_θ,3 that are sharply defined and approximately constant.

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,

$\begin{array}{ll}{\omega }_{\theta }(t)={\omega }_{\theta ,0}+{\widehat{\omega }}_{\theta ,1}(t)sin({\Omega }_{\theta ,1}t)+{\widehat{\omega }}_{\theta ,2}(t)sin({\Omega }_{\theta ,2}t)+\dots ,& (\mathrm{5\theta })\end{array}$

$\begin{array}{ll}{\omega }_{\gamma }(t)={\omega }_{\gamma ,0}+{\widehat{\omega }}_{\gamma ,1}(t)sin({\Omega }_{\gamma ,1}t)+{\widehat{\omega }}_{\gamma ,2}(t)sin({\Omega }_{\gamma ,2}t)+\dots ,& (\mathrm{5\gamma })\end{array}$

where ${\widehat{\omega }}_{\theta ,i}(t)$ are narrow (${\sigma }_t^2\approx 70-80$ 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

Figure 4. Perturbed spectral waves. (A) Series of Gaussian pulses applied to a simulated spectral harmonic. (B) Spectral pulses on Welch's spectrogram. (C) The corresponding spectral θ-wave (continuous line) and the underlying raw traces of the DPT-reconstructed stable θ and γ frequencies. (D) Welch's spectrogram of the spectral θ-waves reproduce the locations of the spectral pulses.

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,

$\begin{array}{ll}{\dot{\varphi }}_i={\omega }_i+{\displaystyle \sum _j}{\lambda }_{ij}cos({\varphi }_i-{\varphi }_j).& (6)\end{array}$

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

$\begin{array}{ll}\Phi ={\displaystyle \sum _i}{\varphi }_i.& (7)\end{array}$

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, s₁, s₂, …, s_N, with a discrete set of harmonics with fixed frequencies,

$\begin{array}{ll}{A}_l={\displaystyle \sum _n}{s}_n{e}^{i\frac{2\pi l}{N}n},& (8)\end{array}$

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, $r(t)={\sum }_p{A}_p{e}^{i{\omega }_{p}t}$, 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,

$\begin{array}{l}s(t)={\displaystyle \sum _p}{A}_p{e}^{i{\omega }_{p}t}+\xi (t)\end{array}$

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, $e^{i\frac{2\pi l}{N}}\to z$,

$\begin{array}{ll}S(z)={\displaystyle \sum _n}{s}_n{z}^n& (9)\end{array}$

where z is a generic complex number. For the oscillatory component of s(t), the sum (9), extended to infinity, yields a meromorphic function,

$\begin{array}{ll}R(z)={\displaystyle \sum _n}{\displaystyle \sum _p}{A}_p{z}^n{e}^{i{\omega }_{p}n\tau }={\displaystyle \sum _p}\frac{a_p{e}^{i{\phi }_p}}{1-z/z_p},& (10)\end{array}$

whose poles, $z_p=e^{-i{\omega }_p\tau }$, and residues, $a_p{e}^{i{\phi }_p}$, 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),

$\begin{array}{ll}{R}_N(z)=\frac{P_{N-1}(z)}{Q_N(z)},& (11)\end{array}$

rapidly approaches R(z) as the degree N of the polynomials P_N(z) and Q_N+1(z) grows, $R(z)=R_N(z)+O(z^{2N})$ (Baker and Graves-Morris, 1996). In particular, the poles z_p of R(z) are approximated by the roots ζ_q of the denominator in (11), Q_N+1(ζ_q) = 0 (Bessis, 1996; Bessis and Perotti, 2009; Perotti et al., 2013).

As for the z-transform of the noise component,

$\Xi (z)={\displaystyle \sum _n}{\xi }_n{z}^n,$

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⁻⁶ − 10⁻⁷ in the standard Euclidean metric in C¹. 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 T_W, centered at time t_i, 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 t_i 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 $\stackrel{-}{s}(t)=2$. The LFPs were originally recorded at the rate S_r = 8 kHz. To increase time resolution in the biologically relevant range of frequencies (f <300 Hz), we interpolated the signal to higher rates (${\tilde{S}}_r=30$ kHz, ${\tilde{S}}_r=36$ kHz or ${\tilde{S}}_r=44$ 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 T_W = 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 $d_F=10^{-6}$. To increase stability, the signals were filtered between 1 and 40 Hz

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.

Funding

MZ and YD are supported by NIH grants R01NS110806 and R01AG074226. CD and DJ are supported by NIH grants R01MH106552 and R01MH112523.

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.

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.

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

Amin, M., and Shirokoff, D. (2010). Flat-top oscillons in an expanding universe. Phys. Rev. D 81, 085045. doi: 10.1103/PhysRevD.81.085045

CrossRef Full Text | Google Scholar

Arenas, A., Díaz-Guilera, A., Kurths, J., Moreno, Y., and Zhou, C. (2008). Synchronization in complex networks. Phys. Rep. 469, 93–153. doi: 10.1016/j.physrep.2008.09.002

CrossRef Full Text | Google Scholar

Aru, J., Aru, J., Priesemann, V., Wibral, M., Lana, L., Pipa, G., et al. (2015). Untangling cross-frequency coupling in neuroscience. Curr. Opin. Neurobiol. 31, 51–61. doi: 10.1016/j.conb.2014.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Babb, S. (1977). Accuracy of planetary theories, particularly for Mars. Isis 68, 426–434. doi: 10.1086/351818

CrossRef Full Text | Google Scholar

Baker, G., and Graves-Morris, P. (1996). Padé Approximants, 2nd Edn.. Cambridge, UK: Cambridge University Press.

Battaglia, F., Sutherland, G., and McNaughton, B. (2004). Hippocampal sharp wave bursts coincide with neocortical “up-state” transitions. Learn. Mem. 11, 697–704. doi: 10.1101/lm.73504

PubMed Abstract | CrossRef Full Text | Google Scholar

Berger, H. (1933). Über das elektrenkephalogramm des menschen. Archiv für Psychiatrie und Nervenkrankheiten 99, 555–574. doi: 10.1007/BF01797193

CrossRef Full Text | Google Scholar

Bessis, D. (1996). Padé approximations in noise filtering. J. Comput. Appl. Math. 66, 85–88. doi: 10.1016/0377-0427(95)00177-8

CrossRef Full Text | Google Scholar

Bessis, D., and Perotti, L. (2009). Universal analytic properties of noise: introducing the J-matrix formalism. J. Phys. A 42, 365202. doi: 10.1088/1751-8113/42/36/365202

CrossRef Full Text | Google Scholar

Boashash, H. (2003). Time Frequency Signal Analysis and Processing: A Comprehensive Reference, 1st Edn. Boston, MA: Elsevier.

Google Scholar

Bogolubsky, I., and Makhankov, V. (1976). Lifetime of pulsating solitons in certain classical models. JETP Lett. 24, 12; ibid. 25, 120 (1977); ibid. 25, 107 (1977).

Google Scholar

Bosnyakova, D., Gabova, A., Kuznetsova, G., Obukhov, Y., Midzyanovskaya, I., Salonin, D., et al. (2006). Time-frequency analysis of spike-wave discharges using a modified wavelet transform. J. Neurosci. Methods 154, 80–88. doi: 10.1016/j.jneumeth.2005.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Brigham, E. (1988). The Fast Fourier Transform And Its Applications. Englewood Cliffs, NJ: Prentice Hall.

Google Scholar

Brun, V., Otnass, M., Molden, S., Steffenach, H., Witter, M., Moser, M. B., et al. (2002). Place cells and place recognition maintained by direct entorhinal-hippocampal circuitry. Science 296, 2243–2246. doi: 10.1126/science.1071089

PubMed Abstract | CrossRef Full Text | Google Scholar

Buzsáki, G. (2011). Rhythms in the Brain. Oxford: Oxford University Press.

Google Scholar

Buzsáki, G., Anastassiou, C., and Koch, C. (2012). The origin of extracellular fields and currents—EEG, ECoG, LFP and spikes. Nat. Rev. Neurosci. 13, 407–420. doi: 10.1038/nrn3241

PubMed Abstract | CrossRef Full Text | Google Scholar

Cannon, J., McCarthy, M., Lee, S., Lee, J., B'́orgers, C., Whittington, M., et al. (2014). Neurosystems: brain rhythms and cognitive processing. Eur. J. Neurosci. 39, 705–719. doi: 10.1111/ejn.12453

PubMed Abstract | CrossRef Full Text | Google Scholar

Carr, M., Karlsson, M., and Frank, L. (2012). Transient slow gamma synchrony underlies hippocampal memory replay. Neuron 75, 700–713. doi: 10.1016/j.neuron.2012.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerda, E., Melo, F., and Rica, S. (1997). Model for subharmonic waves in granular materials. Phys. Rev. Lett. 79, 4570–4573. doi: 10.1103/PhysRevLett.79.4570

CrossRef Full Text | Google Scholar

Cole, S., and Voytek, B. (2017). Brain oscillations and the importance of waveform shape. Trends Cogn. Sci. 21, 137–149. doi: 10.1016/j.tics.2016.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Colgin, L. (2015). Do slow and fast gamma rhythms correspond to distinct functional states in the hippocampal network? Brain Res. 1621, 309–315. doi: 10.1016/j.brainres.2015.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Colgin, L. (2016). Rhythms of the hippocampal network. Nat. Rev. Neurosci. 17, 239–249. doi: 10.1038/nrn.2016.21

PubMed Abstract | CrossRef Full Text | Google Scholar

Copeland, E., Gleiser, M., and Muller, H. (1995). Oscillons: Resonant configurations during bubble collapse. Phys. Rev. D 52, 1920–1933. doi: 10.1103/PhysRevD.52.1920

PubMed Abstract | CrossRef Full Text | Google Scholar

DeVito, J., and Dabaghian, Y. (2014). New Signal Processing Method Reveals Discrete Structure of local Field Potential in Hippocampus. SfN Annual Meeting 2014, Abstract 751.16/UU15.

Folland, G., and Sitaram, A. (1997). The uncertainty principle: a mathematical survey. J. Fourier Anal. Appl. 3, 207–238. doi: 10.1007/BF02649110

CrossRef Full Text | Google Scholar

Fries, P. (2005). A mechanism for cognitive dynamics: neuronal communication through neuronal coherence. Trends Cogn. Sci. 9, 474–480. doi: 10.1016/j.tics.2005.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Froissart, M. (1973). Approximation de Padé, application á la physique des particules élémentaires. Les rencontres physiciens-mathématiciens de Strasbourg-RCP25, 16, 1–13.

Google Scholar

Gallavotti, G. (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 Applicazioni 12, 125–152. doi: 10.48550/arXiv.chao-dyn/9907004

CrossRef Full Text | Google Scholar

Gilewicz, J., and Kryakin, Y. (2003). Froissart doublets in Padé approximation in the case of polynomial noise. J. Comput. Appl. Math. 153, 235–242. doi: 10.1016/S0377-0427(02)00674-X

CrossRef Full Text | Google Scholar

Gilewicz, J., and Pindor, M. (1997). Padé approximants and noise: a case of geometric series. J. Comput. Appl. Math. 87, 199–214. doi: 10.1016/S0377-0427(97)00185-4

CrossRef Full Text | Google Scholar

Gleiser, M. (1994). Pseudostable bubbles. Phys. Rev. D 49, 2978–2981. doi: 10.1103/PhysRevD.49.2978

PubMed Abstract | CrossRef Full Text | Google Scholar

Grünbaum, F. (2003). The Heisenberg inequality for the discrete Fourier transform. Appl. Comput. Harmon. Anal. 15, 163–167. doi: 10.1016/S1063-5203(03)00033-2

CrossRef Full Text | Google Scholar

Hanson, N. (1960). The mathematical power of epicyclical astronomy. Isis 51, 150–158. doi: 10.1086/348869

CrossRef Full Text | Google Scholar

Hoppensteadt, F., and Izhikevich, E. (1997). Weakly Connected Neural Networks. New York, NY: Springer.

Google Scholar

Hoppensteadt, F., and Izhikevich, E. (1998). Thalamo-cortical interactions modeled by weakly connected oscillators: could the brain use FM radio principles? Biosystems 48, 85–94. doi: 10.1016/s0303-2647(98)00053-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Howell, K. (2001). Principles of Fourier Analysis. Portland, OR: CRC Press.

Google Scholar

Izhikevich, E. (1999a). Class 1 neural excitability, conventional synapses, weakly connected networks, and mathematical foundations of pulse-coupled models. IEEE Trans. Neural Netw. 10, 499–507. doi: 10.1109/72.761707

PubMed Abstract | CrossRef Full Text | Google Scholar

Izhikevich, E. (1999b). Weakly pulse-coupled oscillators, FM interactions, synchronization, and oscillatory associative memory. IEEE Trans. Neural Netw. 10, 508–526. doi: 10.1109/72.761708

PubMed Abstract | CrossRef Full Text | Google Scholar

Izhikevich, E. (2000). Phase equations for relaxation oscillators. Siam. J. Appl. Math. 60, 1789–1805. doi: 10.1137/s0036139999351001

CrossRef Full Text | Google Scholar

Jacobsen, E., and Lyons, R. (2003). The sliding DFT. Signal Process. Mag. IEEE 20, 74–80. doi: 10.1109/MSP.2003.1184347

CrossRef Full Text | Google Scholar

Kasuya, S., Kawasaki, M., and Takahashi, F. (2003). I-balls. Phys. Lett. B 559, 99–106. doi: 10.1016/S0370-2693(03)00344-7

CrossRef Full Text | Google Scholar

Kesner, R. (2007). Behavioral functions of the CA3 subregion of the hippocampus. Learn. Mem. 14, 771–781. doi: 10.1101/lm.688207

PubMed Abstract | CrossRef Full Text | Google Scholar

Kopell, N., Kramer, M., Malerba, P., and Whittington, M. (2010). Are different rhythms good for different functions? Front. Hum. Neurosci. 4, 187. doi: 10.3389/fnhum.2010.00187

PubMed Abstract | CrossRef Full Text | Google Scholar

Kropff, E., Carmichael, J., Moser, E., and Moser, M. B. (2021). Frequency of theta rhythm is controlled by acceleration, but not speed, in running rats. Neuron 109, 1–11. doi: 10.1016/j.neuron.2021.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Langston, R., Stevenson, C., Wilson, C., Saunders, I., and Wood, E. (2010). The role of hippocampal subregions in memory for stimulus associations. Behav. Brain Res. 215, 275–291. doi: 10.1016/j.bbr.2010.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozano-Soldevilla, D., ter Huurne, N., and Oostenveld, R. (2016). Neuronal oscillations with non-sinusoidal morphology produce spurious phase-to-amplitude coupling and directionality. Front. Comput. Neurosci. 10, 87. doi: 10.3389/fncom.2016.00087

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, L., Ren, Y., Yu, T., Liu, Z., Wang, S., Tan, L., et al. (2020). Control of locomotor speed, arousal, and hippocampal theta rhythms by the nucleus incertus. Nat. Commun. 11, 262. doi: 10.1038/s41467-019-14116-y

PubMed Abstract | CrossRef Full Text

Neda, Z., Ravasz, E., Brechet, Y., Vicsek, T., and Barabasi, A. L. (2000b). The sound of many hands clapping. Nature 403, 849–850. doi: 10.1038/35002660

PubMed Abstract | CrossRef Full Text | Google Scholar

Neda, Z., Ravasz, E., Vicsek, T., Brechet, Y., and Barabasi, A. L. (2000a). Physics of the rhythmic applause. Phys. Rev. E 61, 6987–6992. doi: 10.1103/PhysRevE.61.6987

PubMed Abstract | CrossRef Full Text | Google Scholar

Newland, D. (2005). An introduction To Random Vibrations, Spectral And Wavelet Analysis. New York, NY: Dover Publications.

Google Scholar

Perotti, L., DeVito, J., Bessis, D., and Dabaghian, Y. (2019). Discrete spectra of brain rhythms. Sci. Rep. 9, 1105. doi: 10.1038/s41598-018-37196-0

PubMed Abstract | CrossRef Full Text

Perotti, L., Regimbau, T., Vrinceanu, D., and Bessis, D. (2014). Identification of gravitational-wave bursts in high noise using Padé filtering. Phys. Rev. D 90, 124047. doi: 10.1103/PhysRevD.90.124047

CrossRef Full Text | Google Scholar

Perotti, L., Vrinceanu, D., and Bessis, D. (2013). Enhanced frequency resolution in data analysis. Am. J. Comput. Math. 3, 242–251. doi: 10.4236/ajcm.2013.33034

CrossRef Full Text | Google Scholar

Perotti, L., and Wojtylak, M. (2018). Matrix methods for Padé approximation: Numerical calculation of poles, zeros and residues. Lin. Alg. App. 548, 95–122. doi: 10.1016/j.laa.2018.03.004

CrossRef Full Text | Google Scholar

Proakis, J., and Manolakis, D. (1996). Digital Signal Processing: Principles, Algorithms and Applications, 3rd Edn. New York, NY: Prentice-Hall.

Richard, G., Titiz, A., Tyler, A., Holmes, G., Scott, R., and Lenck-Santini, P. (2013). Speed modulation of hippocampal theta frequency correlates with spatial memory performance. Hippocampus 23, 1269–1279. doi: 10.1002/hipo.22164

PubMed Abstract | CrossRef Full Text | Google Scholar

Roopun, A., Kramer, M., Carracedo, L., Kaiser, M., Davies, C., Traub, R., et al. (2008). Temporal Interactions between Cortical Rhythms. Front. Neurosci. 2, 145–154. doi: 10.3389/neuro.01.034.2008

PubMed Abstract | CrossRef Full Text

Senior, T., Huxter, J., Allen, K., O'Neill, J., and Csicsvari, J. (2008). Gamma oscillatory firing reveals distinct populations of pyramidal cells in the CA1 region of the hippocampus. J. Neurosci. 28, 2274–2286. doi: 10.1523/JNEUROSCI.4669-07.2008

PubMed Abstract | CrossRef Full Text | Google Scholar

Sitnikova, E., Hramov, A., Koronovsky, A., and van Luijtelaar, G. (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. Methods 180, 304–316. doi: 10.1016/j.jneumeth.2009.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Sreenivasan, K., and D'Esposito, M. (2019). The what, where and how of delay activity. Nat. Rev. Neurosci. 20, 466–481. doi: 10.1038/s41583-019-0176-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Steinhaus, H. (1929). Über die Wahrscheinlichkeit dafür dass der Konvergenzkreis einer Potenzreihe ihre natuerliche Grenze ist. Math. Zeitschrift, 31, 408–416.

Google Scholar

Strogatz, S. (2000). From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D 143, 1–20. doi: 10.1016/S0167-2789(00)00094-4

CrossRef Full Text | Google Scholar

Thut, G., Miniussi, C., and Gross, J. (2012). The functional importance of rhythmic activity in the brain. Curr. Biol. 22, R658–R663. doi: 10.1016/j.cub.2012.06.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Umbanhowar, P., Melo, F., and Swinney, H. (1996). Localized excitations in a vertically vibrated granular layer. Lett. Nat. 382, 793–796. doi: 10.1038/382793a0

CrossRef Full Text | Google Scholar

Van der Waerden, B. (1974). The earliest form of the epicycle theory. J. History Astron. 5, 175–185. doi: 10.1177/002182867400500303

CrossRef Full Text | Google Scholar

Van der Waerden, B. (1982). The motion of venus, mercury and the sun in early greek astronomy. Arch. History Exact. Sci. 26, 99–113. doi: 10.1007/BF00348348

CrossRef Full Text | Google Scholar

van Luijtelaar, G., Hramov, A., Sitnikova, E., and Koronovskii, A. (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. doi: 10.1016/j.clinph.2010.10.038

PubMed Abstract | CrossRef Full Text | Google Scholar

van Vugt, M., Sederberg, P., and Kahana, M. (2007). Comparison of spectral analysis methods for characterizing brain oscillations. J. Neurosci. Methods 162, 49–63. doi: 10.1016/j.jneumeth.2006.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Welch, P. (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. doi: 10.1109/TAU.1967.1161901

CrossRef Full Text | Google Scholar

Yamamoto, J., and Tonegawa, S. (2017). Direct medial entorhinal cortex input to hippocampal CA1 is crucial for extended quiet awake replay. Neuron 96, 217.e4–227.e4. doi: 10.1016/j.neuron.2017.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar