Next Article in Journal
Fractional Fourier Transform: Main Properties and Inequalities
Next Article in Special Issue
Bursting Dynamics of Spiking Neural Network Induced by Active Extracellular Medium
Previous Article in Journal
Dynamic Network Resource Autonomy Management and Task Scheduling Method
Previous Article in Special Issue
A Moving Target Detection Model Inspired by Spatio-Temporal Information Accumulation of Avian Tectal Neurons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Model of a Main Rhythm in Limbic Seizures

by
Maksim V. Kornilov
1,2,*,† and
Ilya V. Sysoev
1,2,†
1
Saratov Branch of Kotel’nikov Institute of Radio Engineering and Electronics of Russian Academy of Sciences, 38 Zelenaya Street, 410019 Saratov, Russia
2
Saratov State University, 83 Astrakhanskaya Street, 410012 Saratov, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(5), 1233; https://doi.org/10.3390/math11051233
Submission received: 3 January 2023 / Revised: 12 February 2023 / Accepted: 23 February 2023 / Published: 3 March 2023

Abstract

:
While synchronization in the brain neural networks has been studied, the emergency of the main oscillation frequency and its evolution at different normal and pathological states remains poorly investigated. We propose a new concept of the formation of a main frequency in limbic epilepsy. The idea is that the main frequency is not a result of the activity of a single cell, but is formed due to collective dynamics in a ring of model neurons connected with delay. The individual cells are in an excitable mode providing no self-oscillations without coupling. We considered the ring of a different number of Hodgkin–Huxley neurons connected with synapses with time delay. We have shown that the proposed circuit can generate oscillatory activity with frequencies close to those experimentally observed. The frequency can be varied by changing the number of model neurons, time delay in synapses and coupling strength. The linear dependence of the oscillation period on both coupling delay and the number of neurons in the ring was hypothesized and checked by means of fitting the values obtained from the numerical experiments to the empirical formula for a constant value of coupling coefficient.

1. Introduction

Synchronization in neural networks traditionally attracts great interest when studying the brain. In particular, for epilepsy research and modelling, a mathematical description of synchronization is very important [1,2], since epilepsy is characterized by the patterns of highly synchronous activity [3]. However, studying only the synchronization is not enough to understand the nature of epilepsy. Limbic epilepsy (also known as temporal lobe epilepsy [4]) is the most common type of epilepsy in humans [5,6]. In secondary generalized forms, such as limbic epilepsy, epileptic rhythms including the main frequency of epileptiform activity in the brain occurs in very small groups of neurons [7]. Only then does the pathological activity spread to large brain networks.
For primary generalized epilepsies, there is a long history of modelling. In some early papers [8,9] epileptic regime was considered to be a result of bifurcation. Then, such ideas became less popular due to them not explaining the observed phenomena well. In particular, bifurcation mechanisms should mostly occur due to changes in ion conductivity [10] or due to significant changes in the connectivity architecture of cells [11]. But for many forms of epilepsy, including limbic epilepsy in rodents, there could be tens of seizures per hour [12]. If we accept a “bifurcation theory of epilepsy” this means that there are reversible chemical reactions or reversible anatomical changes with the time scale of some seconds to tens of seconds which are repeated 100 or more times during 8 h of recording. This seems to be very unlikely. Instead, the ideas of coexisting attractors have became more popular [13,14]: low-amplitude (possibly noise-induced) oscillations with a wide basin (this regime corresponds to a normal behaviour) and high-amplitude oscillations with a narrow basin (this is the pathological regime). Transition between these regimes can occur by means of either noise or external input (e.g., sensory stimulus from trigeminal nerve). Some recent papers also provide a third idea that the epileptic state may be a long-lasting transient which can occur due to a specific organization of the brain network near the limit cycle saddle-node bifurcation [15,16].
Little has been conducted till now for modelling the generator of the principal rhythm of limbic epileptic [17,18]. Considering the larger problem, i.e., the generation of main frequencies of different brain rhythms, the recent paper [19] is of a great interest. In general, there could be two approaches to the construction of mathematical models of brain rhythms. The first approach is to construct a network based on some physiological ideas, and then to investigate it in order to find parameter values at which the network exhibits the desired behaviour. Previously, such an approach was used for modelling thalamocortical system [15,20,21]. The other approach is to construct the desirable behaviour manually in the phase space using well-known subsystems as building blocks as it was performed in [22,23]. This second approach is used here.
In this work, we construct a dynamical system for which the oscillatory regime of interest is directly specified by its properties (assumed by construction). The system is a ring of pyramidal neurons of the hippocampus connected unidirectionally. For individual neurons, the Hodgkin–Huxley equations [24] are written with parameters taken from physiological considerations [25]. When an pulse is generated by a single neuron, it is transmitted to the next in the ring. If the ring is long enough for the pulse to return after the first neuron has left the refractory state, the first neuron fires again and thus periodic oscillations occur. The possibility of the existence of such an oscillatory regime should be determined by the coupling strength and the time of signal propagation along the ring. Proving the existence and investigation of this regime is the primary task of this work.

2. Model

The main model in this work is a ring of identical Hodgkin–Huxley neurons [24] unidirectionally connected with a delay. An example of the coupling architecture for a system of 11 elements is shown in Figure 1. The neurons No. 1–10 make up the main oscillatory ring. Neuron No. 11 was used as an external driving force to begin oscillations; with its help, the system was transferred from the stable point to the oscillatory mode. We also considered the option when, instead of a single neuron, a chain of elements was used to switch to the oscillatory mode; but such a technique did not provide any significant advantages. The equation for each neuron can be written as follows.
d V i d t = [ I s t i m g N a m i 3 h i ( V i V N a ) g K n 4 ( V i V K ) g L ( V i V L ) + κ i ] / C , d m i d t = α m ( V i ) ( 1 m i ) β m ( V i ) m i , d n i d t = α n ( V i ) ( 1 n i ) β n ( V i ) n i , d h i d t = α h ( V i ) ( 1 h i ) β h ( V i ) h i ,
where V is the membrane potential measured across the membrane of the neuron; C is the membrane capacitance; g N a , g K and g L are the maximum conductivities of sodium and potassium channels and leakage, respectively; V N a and V K are reciprocal potentials for sodium and potassium ions; and m, n and h are the gate characteristics of various types of ion channels: m and n describe the activation of sodium and potassium channels, respectively, and h denotes the deactivation of the sodium channel.
In the numerical experiment, we mainly used the parameters of neurons from [25]: C = 1 mF/cm 2 , g N a = 40 mS/cm 2 , g K = 35 mS/cm 2 , g L = 0.3mS/cm 2 , V N a = 55 mV, V K = −77 mV, and V L = −65 mV. These parameters make it possible to obtain dynamics corresponding to self-oscillations in rodent pyramidal neurons. Since we were not interested in self-oscillations in a single neuron, except for one neuron used for network excitation (in Figure 1 this is neuron No. 11), we changed the value of V L = −66.8 mV for all other neurons in the ring to obtain a sub-threshold excitable regime. A special neuron was in oscillatory mode for a short time, then it was disabled (it’s oscillations is shown in Figure 2a)
The functions α i and β p ( p { n , m , h } ) were taken from the same work [25] to reproduce the dynamics of an axon of a pyramid cell.
α n ( V ) = 0.8 ( V 25 ) 1 e ( V 25 ) / 9 , β n ( V ) = 0.002 ( V 25 ) 1 e ( V 25 ) / 9 , α m ( V ) = 0.182 ( V + 35 ) 1 e ( V + 35 ) / 9 , β m ( V ) = 0.124 ( V + 35 ) 1 e ( V + 35 ) / 9 , α h ( V ) = 0.25 e ( V + 90 ) / 12 , β h ( V ) = 0.25 e ( V + 62 ) / 6 e ( V + 90 ) / 12 .
The total effect of the system on each i-th neuron was set using the coupling function κ i as follows:
κ i ( t ) = j = 1 D k i , j ( 1 + tanh ( V j ( t τ ) ) ) ,
where D is the total number of neurons in the system, k i , j is the coefficient of impact of the j neuron on i, and τ is the time delay parameter. The connection was introduced in such a way that, at high-voltage values, the neuron ceased to respond to the driving. All coupling coefficients were set to either k i , j = k for enabled connections or k i , j = 0 for disabled ones.
The network was studied in a bistable regime. The first attractor was a stable fixed point. This regime was inherited from the phase space of individual neurons. The second attractor was a limit cycle, it was of particular interest for us. This cycle appeared due to coupling. It was theoretically hypothesized based on properties of the Hodgkin–Huxley model as described in the Introduction. Since the oscillatory regime emerges with the finite amplitude at some value of k and coexists with the fixed point, the only way it can be born is via a saddle-node bifurcation of a limit cycle [26,27]. In the parameter space, both V L and k can be bifurcation parameters.

3. Simulation

We performed simulations for four different values of k: 30, 40, 50 and 60. For each k value we examined the system for 0.2 ms τ 1 ms changes with a step of 0.05 ms (let us denote them τ i , i = 1 , , i max , i max = 17 ) and 10 D 100 changes by a step of 5 (let us denote them D j , j = 1 , , j max , j max = 19 ). Since the considered systems has a time delay, it was numerically integrated using second-order Runge–Kutta method with time step Δ t = 0.01 ms.
The network demonstrated two regimes. The first one was a fixed point which was inherited from uncoupled oscillators due to the synaptic coupling (3) tending towards zero for V j < < 0 . The second regime was periodic oscillations. This regime was born due to saddle-node bifurcation of a limit cycle (not due to Andronov–Hopf bifurcation); therefore, it coexists with a fixed point. The bifurcation parameters are D, τ and k and an increase in any of them for some reasonable values of others leads to the appearance of the oscillatory regime. Due to bistability in the system, to reach the oscillatory regime we used limited (and short) time external driving from the external neuron No. ( D + 1 ) . The end time moment of this driving is shown by a blue line in Figure 2 (the starting time of the driving is t = 0 ).
The studied network was constructed to produce periodical oscillations, modelling the main rhythm of epileptic activity. Therefore, the simulation mostly aimed to prove the existence of the oscillatory regime and to obtain the dependence of the main oscillation frequency f on the number of neurons D and the coupling delay τ . So, we studied the dynamics at those values of k, D and τ at which the oscillatory regime appears, see Figure 3 where f is plotted by colour.
Considering Figure 3, one can suggest that the dependency of f on D and τ is hyperbolic. Let us formulate this hypothesis and test it.

4. Approximation

It is easier to consider linear dependency rather than hyperbolic dependency. Therefore, let us formulate our hypothesis about the period of oscillations T rather than about the frequency f, assuming f = 1 / T .
Hypothesis 1. 
For the fixed coupling coefficient k the oscillation period in the ring of unidirectionally coupled Hodgkin–Huxley neurons (1) depends only on the number D of oscillators in the ring and the delay time τ, with this dependence being the linear superposition of D and τ D with constant parameters.
To prove or disprove the formulated hypothesis analytically, it is necessary to obtain an exact solution for the Equation (1). For much simpler systems of phase oscillators similar problem was studied analytically in [28]. But for a ring of Hodgkin–Huxley neurons this does not seem achievable in the current state of science. Therefore, we chose another way. We solved the system (1) numerically for different D and τ values and calculated the oscillation period T in each case. Then, we fitted the dependency T ( D , τ ) reported above to the following Formula (4):
T ˜ ( D , τ , k ) = T 0 ( k ) + γ ( k ) τ D + ε ( k ) D ,
where T ˜ is the approximation of the period T measured in the simulation, and T 0 , ε and γ are constants depending only on k, with T 0 and ε possessing the dimension of time and γ being dimensionless. We should note that ε is a parameter of inertia of a single neuron with an axon (as in most time delay systems, see [29,30]), i.e., the time shift which the neuron provides itself without coupling delay τ . The formula was fitted for each k separately using all Q = i max j max values of T measured for different D and τ values. To estimate the goodness of fit, we calculated the normalized error of approximation σ appr 2 as follows:
σ appr 2 ( k ) = 1 Q σ T 2 ( k ) i = 1 i max j = 1 j max T ( k ) T ˜ ( k ) 2 ,
where σ T ( k ) is the empirical standard deviation for T at the fixed k.
Values of both reconstructed T ˜ and measured T periods of oscillations are represented in Figure 4. These plots show that the reconstructed values are much closer to the measured ones with a greater number of neurons in the ring and with higher values of the coupling coefficient. The results of approximation are briefly summarized in the Table 1. One can see that the approximation error is quite moderate and in all considered cases it is less than 10 3 . This error can be expressed by both model imperfection and calculation errors. First, the period T was measured in simulation with a precision of up to Δ t = 0.01 ms, i.e., the approximated value was measured with a relative error ∼ 10 4 . Second, the method of numerical integration was imperfect providing similar error O ( Δ t ) . The approximation procedure also has some impact on the error, though the use of the LAPACK procedure from Python code (see [31]) should provide much less error than ∼ 10 4 .

5. Conclusions and Discussion

In this study, a system of Hodgkin–Huxley neuron oscillators connected in a ring was considered. The coupling was unidirectional. The oscillations were excited using the short-term activity of a starting neuron. The non-zero coupling coefficients between neurons were equal. It was shown to be possible to control the signal frequency f by changing the number of neurons in the ring D and the coupling delay time τ . Changing D affects the frequency change more strongly than τ . In addition, increasing the values of the coupling coefficient between neurons in the ring made it possible to obtain modes with higher oscillation frequencies. The whole range of frequencies characteristic to limbic seizures [12] can be reproduced by the proposed approach.
Then, the dependence of T (the oscillation period) on D and τ was analysed. From this we formulated the hypothesis that this dependence is linear. To check this hypothesis, a mathematical model was constructed, the coefficients of which were fitted to the numerically obtained solutions. This fitting showed that the model has sufficiently high accuracy (relative error σ appr 2 < 10 3 ).
We did not formulate our hypothesis as a theorem since we were not completely sure whether the Formula (4) was exact and an analytical proof or refutation is difficulty to perform. Increasing the calculation precision could provide some additional information if the approximation error reduces with a decrease in Δ t . However, for modelling purposes the obtained precision is enough.
When fitting the Formula (4) to the simulated data, the three parameters behave very differently. The parameter γ is mostly constant for all values of k, which seems to be natural. This parameter shows that the impact of the delay in the ring as it should not depend on the coupling from its physical sense. The parameter ε has a physical sense of an inertia parameter of an individual neuron. One can see that it is reduced with an increase in k. This means that unlike simple time delay systems [29], the Hodgkin–Huxley neuron the inertia is non-linear and the inertia coefficient depends on the signal amplitude and shape. In particular, a decrease in ε with an increase in k should be the result of more narrow impulses in the group for larger k values. It is also interesting that we obtained a non-zero value of T 0 in the model (4), which decreases with an increase in the coupling coefficient k. This constant cannot be described physically, since it gives a non-zero period for D = 0 . However, for small values of D the oscillatory regime, which was of particular interest to us, does not appear at all since the neurons are deeply depolarized and do not respond to the external impulse of the previous neuron.
In our simulations we observed the only oscillatory mode for each given set of parameters, i.e., there were only two stable regimes in the system: the stable point and the observed regime. However, one could theoretically expect that with a sufficiently large number of oscillators in the system, not one, but several oscillatory modes with multiple frequencies could exist. In addition to the main regime observed in the framework of this work, there could be regimes with two, three times, etc. lower frequencies. These modes corresponded to two, three, etc. simultaneously moving pulses along the ring.

Author Contributions

Conceptualization, I.V.S.; methodology, I.V.S.; software, M.V.K.; validation, I.V.S.; formal analysis, I.V.S.; writing—original draft preparation, I.V.S. and M.V.K.; visualization, M.V.K.; supervision, I.V.S.; project administration, I.V.S.; funding acquisition, I.V.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant number 19-72-10030-P, https://rscf.ru/project/19-72-10030/.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jefferys, J.G.; Jiruska, P.; de Curtis, M.; Avoli, M. Limbic network synchronization and temporal lobe epilepsy. In Jasper’s Basic Mechanisms of the Epilepsies [Internet], 4th ed.; National Center for Biotechnology Information (US): Bethesda, MD, USA, 2012. [Google Scholar]
  2. Jiruska, P.; De Curtis, M.; Jefferys, J.G.; Schevon, C.A.; Schiff, S.J.; Schindler, K. Synchronization and desynchronization in epilepsy: Controversies and hypotheses. J. Physiol. 2013, 591, 787–797. [Google Scholar] [CrossRef] [Green Version]
  3. Avoli, M.; D’Antuono, M.; Louvel, J.; Köhling, R.; Biagini, G.; Pumain, R.; D’Arcangelo, G.; Tancredi, V. Network and pharmacological mechanisms leading to epileptiform synchronization in the limbic system in vitro. Prog. Neurobiol. 2002, 68, 167–207. [Google Scholar] [CrossRef]
  4. Berg, A.T.; Berkovic, S.F.; Brodie, M.J.; Buchhalter, J.; Cross, J.H.; Van Emde Boas, W.; Engel, J.; French, J.; Glauser, T.A.; Mathern, G.W.; et al. Revised terminology and concepts for organization of seizures and epilepsies: Report of the ILAE Commission on Classification and Terminology, 2005–2009. Epilepsia 2010, 51, 676–685. [Google Scholar] [CrossRef]
  5. Wiebe, S. Epidemiology of temporal lobe epilepsy. Can. J. Neurol. Sci. 2000, 27, S6–S10. [Google Scholar] [CrossRef]
  6. Téllez-Zenteno, J.F.; Hernández-Ronquillo, L. A review of the epidemiology of temporal lobe epilepsy. Epilepsy Res. Treat. 2012, 2012, 630853. [Google Scholar] [CrossRef] [Green Version]
  7. Paz, J.T.; Huguenard, J.R. Microcircuits and their interactions in epilepsy: Is the focus out of focus? Nat. Neurosci. 2015, 18, 351–359. [Google Scholar] [CrossRef] [Green Version]
  8. Babloyantz, A.; Destexhe, A. Low-dimensional chaos in an instance of epilepsy. Proc. Natl. Acad. Sci. USA 1986, 83, 3513–3517. [Google Scholar] [CrossRef] [Green Version]
  9. Lopes Da Silva, F.; Blanes, W.; Kalitzin, S.N.; Parra, J.; Suffczynski, P.; Velis, D.N. Epilepsies as dynamical diseases of brain systems: Basic models of the transition between normal and epileptic activity. Epilepsia 2003, 44, 72–83. [Google Scholar] [CrossRef]
  10. Destexhe, A.; Babloyantz, A.; Sejnowski, T. Ionic mechanisms for intrinsic slow oscillations in thalamic relay neurons. Biophys J. 1993, 65, 1538–1552. [Google Scholar] [CrossRef]
  11. Spencer, S. Neural Networks in Human Epilepsy: Evidence of and Implications for Treatment. Epilepsia 2002, 43, 219–227. [Google Scholar] [CrossRef]
  12. Bertram, E.H. Functional Anatomy of Spontaneous Seizures in a Rat Model of Limbic Epilepsy. Epilepsia 1997, 38, 95–105. [Google Scholar] [CrossRef]
  13. Suffczynski, P.; Kalitzin, S.; Lopes Da Silva, F. Dynamics of non-convulsive epileptic phenomena modeled by a bistable neuronal network. Neuroscience 2004, 126, 467–484. [Google Scholar] [CrossRef]
  14. Taylor, P.N.; Wang, Y.; Goodfellow, M.; Dauwels, J.; Moeller, F.; Stephani, U.; Baier, G. A Computational Study of Stimulus Driven Epileptic Seizure Abatement. PLoS ONE 2014, 9, e114316. [Google Scholar] [CrossRef] [Green Version]
  15. Medvedeva, T.M.; Sysoeva, M.V.; Lüttjohann, A.; van Luijtelaar, G.; Sysoev, I.V. Dynamical mesoscale model of absence seizures in genetic models. PLoS ONE 2020, 15, e239125. [Google Scholar] [CrossRef]
  16. Kapustnikov, A.A.; Sysoeva, M.V.; Sysoev, I.V. Transient dynamics in a class of mathematical models of epileptic seizures. Commun. Nonlinear Sci. Numer. Simul. 2022, 109, 106284. [Google Scholar] [CrossRef]
  17. Klimesch, W. Memory processes, brain oscillations and EEG synchronization. Int. J. Psychophysiol. 1996, 24, 61–100. [Google Scholar] [CrossRef]
  18. Bennett, M.V.L.; Zukin, R.S. Electrical Coupling and Neuronal Synchronization in the Mammalian Brain. Neuron 2004, 41, 495–511. [Google Scholar] [CrossRef] [Green Version]
  19. Kurkin, S.A.; Kulminskiy, D.D.; Ponomarenko, V.I.; Prokhorov, M.D.; Astakhov, S.V.; Hramov, A.E. Central pattern generator based on self-sustained oscillator coupled to a chain of oscillatory circuits. Chaos Interdiscip. J. Nonlinear Sci. 2022, 32, 033117. [Google Scholar] [CrossRef]
  20. Bazhenov, M.; Timofeev, I.; Steriade, M.; Sejnowski, T.J. Model of thalamocortical slow-wave sleep oscillations and transitions to activated states. J. Neurosci. 2002, 22, 8691–8704. [Google Scholar] [CrossRef] [Green Version]
  21. Hill, S.; Tononi, G. Modeling sleep and wakefulness in the thalamocortical system. J. Neurophysiol. 2005, 93, 1671–1698. [Google Scholar] [CrossRef]
  22. Afraimovich, V.S.; Rabinovich, M.I.; Varona, P. Heteroclinic Contours In Neural Ensembles And The Winnerless Competition Principle. Int. J. Bifurc. Chaos 2004, 14, 1195–1208. [Google Scholar] [CrossRef] [Green Version]
  23. Rabinovich, M.I.; Varona, P.; Tristan, I.; Afraimovich, V.S. Chunking dynamics: Heteroclinics in mind. Front. Comput. Neurosci. 2014, 8, 22. [Google Scholar] [CrossRef] [Green Version]
  24. Hodgkin, A.; Huxley, A. A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 1952, 117, 500–544. [Google Scholar] [CrossRef]
  25. Talukdar, S.; Shrivastava, R.; Ghosh, S. Modeling activity-dependent reduction in after hyper-polarization with Hodgkin-Huxley equation of action potential. Biomed. Phys. Eng. Express 2019, 5, 047001. [Google Scholar] [CrossRef]
  26. Kuznetsov, Y.A. Elements of Applied Bifurcation Theory; Springer: New York, NY, USA, 2004; p. 632. [Google Scholar] [CrossRef]
  27. Gonchenko, A.S.; Gonchenko, S.V.; Kazakov, A.O.; Kozlov, A.D.; Bakhanova, Y.V. Mathematical theory of dynamical chaos and its applications: Review Part 2. Spiral chaos of three-dimensional flows. Izv. Vuz. Appl. Nonlinear Dyn. 2019, 27, 7–52. [Google Scholar] [CrossRef]
  28. Omelchenko, O.E.; Hauptmann, C.; Maistrenko, Y.L.; Tass, P.A. Collective dynamics of globally coupled phase oscillators under multisite delayed feedback stimulation. Phys. D 2008, 237, 365–384. [Google Scholar] [CrossRef]
  29. Kuznetsov, S.P. Complex dynamics of oscillators with delayed feedback (review). Radiophys. Quantum Electron. 1982, 25, 996–1009. [Google Scholar] [CrossRef]
  30. Wernecke, H.; Sándor, B.; Gros, C. Chaos in time delay systems, an educational review. Phys. Rep. 2019, 824, 1–40. [Google Scholar] [CrossRef] [Green Version]
  31. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
Figure 1. An example of a connection diagram for a system of 11 Hodgkin–Huxley [24] neurons, 10 of which are connected by a unidirectional couplings with a delay; the neuron No. 11 is used to begin oscillations.
Figure 1. An example of a connection diagram for a system of 11 Hodgkin–Huxley [24] neurons, 10 of which are connected by a unidirectional couplings with a delay; the neuron No. 11 is used to begin oscillations.
Mathematics 11 01233 g001
Figure 2. Time series of the transmembrane voltage (V variable of the model) by D = 10 , τ = 0.2 ms and k = 40 for the starter neuron No. 11 (a) and for the ring neurons No. 1 (b) and No. 6 (c). Vertical blue line represents the oscillating time of neuron No. 11.
Figure 2. Time series of the transmembrane voltage (V variable of the model) by D = 10 , τ = 0.2 ms and k = 40 for the starter neuron No. 11 (a) and for the ring neurons No. 1 (b) and No. 6 (c). Vertical blue line represents the oscillating time of neuron No. 11.
Mathematics 11 01233 g002
Figure 3. Dependence of the frequency f on the number of neuron-like oscillators and coupling delay time for different values of the coupling coefficient: k = 30 (a), k = 40 (b), k = 50 (c) and k = 60 (d) correspondingly.
Figure 3. Dependence of the frequency f on the number of neuron-like oscillators and coupling delay time for different values of the coupling coefficient: k = 30 (a), k = 40 (b), k = 50 (c) and k = 60 (d) correspondingly.
Mathematics 11 01233 g003
Figure 4. Comparison of the time series period of Hodgkin–Huxley neurons obtained as a result of the numerical experiments T (black lines) and reconstructed (red dashed lines) in accordance with the model (4) plotted at fixed D for the coupling coefficients k = 30 (ac), k = 40 (df), k = 50 (gi), and k = 60 (jl), respectively.
Figure 4. Comparison of the time series period of Hodgkin–Huxley neurons obtained as a result of the numerical experiments T (black lines) and reconstructed (red dashed lines) in accordance with the model (4) plotted at fixed D for the coupling coefficients k = 30 (ac), k = 40 (df), k = 50 (gi), and k = 60 (jl), respectively.
Mathematics 11 01233 g004
Table 1. Values of the parameters T 0 , ε and γ in the Formula (4) fitted for different k with the mean-squared-normalized error of approximation σ appr 2 .
Table 1. Values of the parameters T 0 , ε and γ in the Formula (4) fitted for different k with the mean-squared-normalized error of approximation σ appr 2 .
k T 0 , ms γ ε , ms σ appr 2
3014.270.9853.08 9.8 × 10 4
4013.490.9602.00 6.4 × 10 4
5010.810.9601.49 4.0 × 10 4
608.790.9501.21 2.1 × 10 4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kornilov, M.V.; Sysoev, I.V. Mathematical Model of a Main Rhythm in Limbic Seizures. Mathematics 2023, 11, 1233. https://doi.org/10.3390/math11051233

AMA Style

Kornilov MV, Sysoev IV. Mathematical Model of a Main Rhythm in Limbic Seizures. Mathematics. 2023; 11(5):1233. https://doi.org/10.3390/math11051233

Chicago/Turabian Style

Kornilov, Maksim V., and Ilya V. Sysoev. 2023. "Mathematical Model of a Main Rhythm in Limbic Seizures" Mathematics 11, no. 5: 1233. https://doi.org/10.3390/math11051233

APA Style

Kornilov, M. V., & Sysoev, I. V. (2023). Mathematical Model of a Main Rhythm in Limbic Seizures. Mathematics, 11(5), 1233. https://doi.org/10.3390/math11051233

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop