# Rapid Spectral Dynamics in Hippocampal Oscillons

^{1}Department of Neurology, McGovern Medical Center at Houston, The University of Texas, Houston, TX, United States^{2}Department of Neuroscience, Baylor College of Medicine, Houston, TX, United States^{3}Department of Physics, Texas Southern University, Houston, TX, United States

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), ${\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**. 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,

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 interpolating^{2} 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*_{1}, *t*_{2}, …, *t*_{n}],

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**. 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

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**. 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, ${\stackrel{~}{\omega}}_{\theta ,0}\approx 8$ Hz (dotted line) and ${\stackrel{~}{\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 ${\stackrel{~}{\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,

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**. 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},

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, *s*_{1}, *s*_{2}, …, *s*_{N}, 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, $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,

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$,

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, ${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),

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,

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 *C*^{1}. 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 (${\stackrel{~}{S}}_{r}=30$ kHz, ${\stackrel{~}{S}}_{r}=36$ kHz or ${\stackrel{~}{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}=1{0}^{-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

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

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

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

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

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

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

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

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

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).

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

Edited by:

Daya Shankar Gupta, Husson University, United StatesReviewed by:

Joao Streibel, Federal University of Santa Catarina, BrazilMarcelo Gleiser, Dartmouth College, United States

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