Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

The mechanism of Parkinson oscillation in the cortex: Possible evidence in a feedback model projecting from the globus pallidus to the cortex


  • The origin, location and cause of Parkinson's oscillation are not clear at present. In this paper, we establish a new cortex-basal ganglia model to study the origin mechanism of Parkinson beta oscillation. Unlike many previous models, this model includes two direct inhibitory projections from the globus pallidus external (GPe) segment to the cortex. We first obtain the critical calculation formula of Parkinson's oscillation by using the method of Quasilinear analysis. Different from previous studies, the formula obtained in this paper can include the self-feedback connection of GPe. Then, we use the bifurcation analysis method to systematically explain the influence of some key parameters on the oscillation. We find that the bifurcation principle of different cortical nuclei is different. In general, the increase of the discharge capacity of the nuclei will cause oscillation. In some special cases, the sharp reduction of the discharge rate of the nuclei will also cause oscillation. The direction of bifurcation simulation is consistent with the critical condition curve. Finally, we discuss the characteristics of oscillation amplitude. At the beginning of the oscillation, the amplitude is relatively small; with the evolution of oscillation, the amplitude will gradually strengthen. This is consistent with the experimental phenomenon. In most cases, the amplitude of cortical inhibitory nuclei (CIN) is greater than that of cortical excitatory nuclei (CEX), and the two direct inhibitory projections feedback from GPe can significantly reduce the amplitude gap between them. We calculate the main frequency of the oscillation generated in this model, which basically falls between 13 and 30 Hz, belonging to the typical beta frequency band oscillation. Some new results obtained in this paper can help to better understand the origin mechanism of Parkinson's disease and have guiding significance for the development of experiments.

    Citation: Minbo Xu, Bing Hu, Weiting Zhou, Zhizhi Wang, Luyao Zhu, Jiahui Lin, Dingjiang Wang. The mechanism of Parkinson oscillation in the cortex: Possible evidence in a feedback model projecting from the globus pallidus to the cortex[J]. Mathematical Biosciences and Engineering, 2023, 20(4): 6517-6550. doi: 10.3934/mbe.2023281

    Related Papers:

    [1] Omar Hassan Mahmood, Mustafa Sh. Aljanabi, Farouk M. Mahdi . Effect of Cu nanoparticles on microhardness and physical properties of aluminum matrix composite prepared by PM. AIMS Materials Science, 2025, 12(2): 245-257. doi: 10.3934/matersci.2025013
    [2] Supriya Rattan, Derek Fawcett, Gerrard Eddy Jai Poinern . Williamson-Hall based X-ray peak profile evaluation and nano-structural characterization of rod-shaped hydroxyapatite powder for potential dental restorative procedures. AIMS Materials Science, 2021, 8(3): 359-372. doi: 10.3934/matersci.2021023
    [3] Xuan Luc Le, Nguyen Dang Phu, Nguyen Xuan Duong . Enhancement of ferroelectricity in perovskite BaTiO3 epitaxial thin films by sulfurization. AIMS Materials Science, 2024, 11(4): 802-814. doi: 10.3934/matersci.2024039
    [4] Saif S. Irhayyim, Hashim Sh. Hammood, Hassan A. Abdulhadi . Effect of nano-TiO2 particles on mechanical performance of Al–CNT matrix composite. AIMS Materials Science, 2019, 6(6): 1124-1134. doi: 10.3934/matersci.2019.6.1124
    [5] Saif S. Irhayyim, Hashim Sh. Hammood, Anmar D. Mahdi . Mechanical and wear properties of hybrid aluminum matrix composite reinforced with graphite and nano MgO particles prepared by powder metallurgy technique. AIMS Materials Science, 2020, 7(1): 103-115. doi: 10.3934/matersci.2020.1.103
    [6] Denise Arrozarena Portilla, Arturo A. Velázquez López, Rosalva Mora Escobedo, Hernani Yee Madeira . Citrate coated iron oxide nanoparticles: Synthesis, characterization, and performance in protein adsorption. AIMS Materials Science, 2024, 11(5): 991-1012. doi: 10.3934/matersci.2024047
    [7] Christoph Janiak . Inorganic materials synthesis in ionic liquids. AIMS Materials Science, 2014, 1(1): 41-44. doi: 10.3934/matersci.2014.1.41
    [8] Abdulkader A. Annaz, Saif S. Irhayyim, Mohanad L. Hamada, Hashim Sh. Hammood . Comparative study of mechanical performance between Al–Graphite and Cu–Graphite self-lubricating composites reinforced by nano-Ag particles. AIMS Materials Science, 2020, 7(5): 534-551. doi: 10.3934/matersci.2020.5.534
    [9] Purabi R. Ghosh, Derek Fawcett, Michael Platten, Shashi B. Sharma, John Fosu-Nyarko, Gerrard E. J. Poinern . Sustainable green chemical synthesis of discrete, well-dispersed silver nanoparticles with bacteriostatic properties from carrot extracts aided by polyvinylpyrrolidone. AIMS Materials Science, 2020, 7(3): 269-287. doi: 10.3934/matersci.2020.3.269
    [10] Mohsen Safaei, Mohammad Salmani Mobarakeh, Bahram Azizi, Ehsan Shoohanizad, Ling Shing Wong, Nafiseh Nikkerdar . Optimization of synthesis of cellulose/gum Arabic/Ag bionanocomposite for antibacterial applications. AIMS Materials Science, 2025, 12(2): 278-300. doi: 10.3934/matersci.2025015
  • The origin, location and cause of Parkinson's oscillation are not clear at present. In this paper, we establish a new cortex-basal ganglia model to study the origin mechanism of Parkinson beta oscillation. Unlike many previous models, this model includes two direct inhibitory projections from the globus pallidus external (GPe) segment to the cortex. We first obtain the critical calculation formula of Parkinson's oscillation by using the method of Quasilinear analysis. Different from previous studies, the formula obtained in this paper can include the self-feedback connection of GPe. Then, we use the bifurcation analysis method to systematically explain the influence of some key parameters on the oscillation. We find that the bifurcation principle of different cortical nuclei is different. In general, the increase of the discharge capacity of the nuclei will cause oscillation. In some special cases, the sharp reduction of the discharge rate of the nuclei will also cause oscillation. The direction of bifurcation simulation is consistent with the critical condition curve. Finally, we discuss the characteristics of oscillation amplitude. At the beginning of the oscillation, the amplitude is relatively small; with the evolution of oscillation, the amplitude will gradually strengthen. This is consistent with the experimental phenomenon. In most cases, the amplitude of cortical inhibitory nuclei (CIN) is greater than that of cortical excitatory nuclei (CEX), and the two direct inhibitory projections feedback from GPe can significantly reduce the amplitude gap between them. We calculate the main frequency of the oscillation generated in this model, which basically falls between 13 and 30 Hz, belonging to the typical beta frequency band oscillation. Some new results obtained in this paper can help to better understand the origin mechanism of Parkinson's disease and have guiding significance for the development of experiments.



    Consider the semilinear elliptic equation

    Δu=f(u),inRn. (1.1)

    A solution u is (linearized) stable if

    Rn[|φ|2f(u)φ2]0,φC0(Rn). (1.2)

    The Morse index of a solution is defined to be the maximal dimension of the negative space for this quadratic form. A solution with finite Morse index is therefore not too unstable.

    In 2000s, Dancer wrote a series of papers on stable and finite Morse index solutions, [13,14,15,16,17,18,19,21] (see also his survey [20] at 2010 ICM and the summary in Du [35,Section 8]). He obtained various classification results about these solutions and applied them to the study of equations with small parameters in bounded domains and global bifurcation problems. Many results on stable and finite Morse index solutions have appeared since then. We refer the reader to the monograph of Dupaigne [40] for an extensive list of results and references up to 2010.

    In this paper we review some recent results about stable and finite Morse index solutions, mostly between 2010-2020. We will restrict our attention to the Liouville (-Bernstein-De Giorgi) type results on stable and finite Morse index solutions defined on the entire Rn, although there are many studies devoted to problems on bounded or unbounded domains in Rn and on complete Riemannian manifolds. As in Dancer's work, by a blow up procedure, these Liouville (-Bernstein-De Giorgi) type results should be useful in many other problems. We also refer the reader to a recent survey of Takahashi [83], which is focused on the a priori estimate aspect for stable solutions.

    In [20], Dancer proposed the following conjecture:

    Assume n8 and u is a bounded, linearized stable solutions of (1.1).

    Then either u is constant on Rn or, after a rotation of axis, u=u(xn).

    If u=u(xn), the solution is called one dimensional. This is equivalent to saying that all level sets of u are hyperplanes. In the following, we say the dimension N is critical if for n<N there are only trivial (constant or one dimensional) stable solutions on Rn, while for nN, there exist nontrivial stable solutions. This critical dimension always exists and it is usually finite.

    For nonegative nonlinearities f0, in a recent breakthrough work [7], Cabre et. al. obtained several a priori estimates on stable solutions, and used them to settle Brezis's conjecture on the L regularity of "extremal solutions" [5]. In [41], Dupaigne and Farina used these estimates to establish the following Liouville theorem:

    Theorem 2.1. Assume that f0 and it is locally Lipschitz. Suppose u is a stable solution of (1.1) and it is bounded below. If n10, then u is constant.

    Hence in this case, the critical dimension is 11. This theorem can be trivially extended to nonpositive nonlinearities f0 (subharmonic functions) if the solution is bounded from above. But in the subharmonic case, in view of results on the singular elliptic equation

    Δu=up (2.1)

    (see Esposito-Ghoussoub-Guo [42], Meadows [72], Ma-Wei [71] and Du-Guo [36]), to consider Liouville property for stable solutions, a more natural class is those solutions with a suitable polynomial growth at infinity (or even without any growth condition). In this case, it seems that the critical dimension should be 7. (In fact, (2.1) with p=1 is almost the worst case.) More precisely, we want to know if the following Liouville property holds.

    Problem 1. Assume f0 is locally Lipschitz on (0,+), and it satisfies

    10f(u)du=. (2.2)

    Then there is no positive, stable solution on Rn, provided that n<7.

    Here the condition (2.2) implies that there does not exist one dimensional stable solution. Note that by the above remark, no assumption on the growth of solutions is added.

    If f is sign-changing but there is no one dimensional stable solution of (1.1), currently there is still not too much progress on this conjecture in dimensions n3.

    If there exists a one dimensional stable solution of (1.1), then a double well structure is associated to this equation, that is, there exist two constants a<a+ such that a primitive function of f, F satisfies

    F(a±)=0,F(t)<0in(a,a+).

    Here a±=limxn±g(xn), where g is a one dimensional stable (equivalently, monotone) solution.

    Under this assumption, Dancer's conjecture is closely related to De Giorgi conjecture (De Giorgi [28]) and stable De Giorgi conjecture about Allen-Cahn equation

    Δu=uu3. (2.3)

    By the way, based on the result of Pacard-Wei [76], in this case, the critical dimension should be 8.

    Recall that De Giorgi conjecture states that

    Suppose u is a solution of the Allen-Cahn equation (2.3) in Rn. If xnu>0 and n8, then after a rotation of axis, u=u(xn).

    While the stable De Giorgi conjecture states that

    Suppose u is a stable solution of the Allen-Cahn equation (2.3) in Rn. If n7, then after a rotation of axis, u=u(xn).

    In the Allen-Cahn equation, the standard double well potential (1u2)2/4 contains only one critical point in the interval (1,1). But in Dancer's conjecture, besides stating that if there is no one dimensional stable profile, then stable solutions must be constant, it also includes the claim when F contains more than one critical points in the interval (a,a+). Here it is possible that a sub-double-well structure exists in (a,a+), which makes the situation more complicate.

    The De Giorgi conjecture has been solved by Ghoussoub-Gui [55] in dimension 2, by Ambrosio-Cabre [4] in dimension 3 and under an additional assumption in dimension 4n8 by Savin [78]. It was also shown to be not true in n9 by Del Pino-Kowalczyk-Wei [32].

    By an observation of Dancer, it is known that the method introduced in Ghoussoub-Gui [55] and Ambrosio-Cabre [4] can be used to prove the stable De Giorgi conjecture in dimension 2. By the time of this writing, this is still the only proven case for stable De Giorgi conjecture. The main difficulty seems to be that, by now the only tool to tackle such problems is the Liouville property for the degenerate equation

    div(σ2φ)=0, (2.4)

    but by Gazzola [54], there is no hope to use this in dimensions n3.

    Even the following weaker version of the stable De Giorgi conjecture is still open.

    Problem 2. Suppose u is a stable solution of the Allen-Cahn equation (2.3) in Rn, satisfying the natural energy growth bound

    BR(0)[12|u|2+14(1u2)2]CRn1,R>0. (2.5)

    If 4n7, then u is one dimensional.

    We say the energy growth condition (2.5) is natural because it is satisfied by minimizing solutions. (Whether this condition holds for stable solutions is another unknown point in the stable De Giorgi conjecture.) If n=3, the energy growth bound is quadratic and we can repeat the above proof of stable De Giorgi conjecture in dimension 2. This fact has been used by Chodosh-Mantoulidis in [12] to establish curvature estimates for Allen-Cahn equation on three dimensional manifolds. But as explained above, this approach does not work without such a quadratic bound. A possible way to solve Problem 2 is to establish a sheeting theorem similar to the ones for stable minimal hypersurfaces (see Schoen-Simon [80] and Wickramasekera [93]), but a geometric estimate (i.e. the oscillation estimate in [80]) is not known in this semilinear setting.

    Another missing geometric estimate for semilinear elliptic equations is the one correspondent to the famous Simons inequality (Simons [82]) for minimal hypersurfaces, which is a fundamental tool in the study of stable minimal hypersurfaces (see e.g. [81], [79]). This difficulty is also encountered in the study of Bernstein property for the Alt-Caffarelli one phase free boundary problem (see Alt-Caffarelli [2]) and nonlocal minimal surfaces (see Caffarelli-Roquejoffre-Savin [8]).

    Recall that the Alt-Caffarelli one phase free boundary problem is

    {Δu=0in{u>0},|u|=1on{u>0}. (2.6)

    There are some variants, such as the problem studied by Phillips [77] and Alt-Phillips [3]

    Δu=upχ{u>0},0<p<1, (2.7)

    as well as various approximations to these problems, e.g.

    Δu=fε(u),

    where fε is an approximate Dirac at 0 (see Fernández-Real and Ros-Oton [51]). These problems are similar to the ones considered in Problem 1, but with an integrable condition rather than the divergence condition (2.2). This integrable condition implies that there is a one dimensional monotone solution, where the solution has a non-empty dead core {u=0}, leading to a free boundary condition on {u=0}.

    Although most studies on these free boundary problems are focused on minimizing solutions or viscosity solutions, recently there arises some interest in higher energy critical points, see Jerison-Perera [64]. To understand these solutions, the stability condition should play an important role. Indeed, even for minimizing solutions, to prove the optimal partial regularity of free boundaries, one needs the classification of stable, homogeneous solutions just as in the Bernstein problem for minimal hypersurfaces, see Weiss [92]. For the Alt-Caffarelli one phase free boundary problem, it is conjectured that the critical dimension is 7. Currently, the Bernstein property is known to be true if n4 (by Caffarelli-Jerison-Kenig [10] and Jerison-Savin [65]), and not true if n7 (by De Silva-Jerison [29]).

    In recent years, we see also much progress on De Giorgi conjecture for fractional Allen-Cahn equation

    (Δ)su=uu3. (2.8)

    In particular, Figalli-Serra [52] solved the stable De Giorgi conjecture for the s=1/2 (half Laplacian) case in dimension 3 and the corresponding De Giorgi conjecture in dimension 4. Later their energy estimates were generalized to other fractional Laplacians by Gui-Li [56]. These estimates are optimal for s<1/2. In this case, the fractional Allen-Cahn equation is related to the nonlocal minimal surfaces. But for nonlocal minimal surfaces, it is still not known in general what is the critical dimension for the stable Bernstein property. Only some perturbative results have been obtained for those s sufficiently close to 1/2, see Caffarelli-Valdinoci [9] and Cabre-Cinti-Serra [6].

    For equations enjoying a scaling invariance, much progress has been obtained in the last decade. This is because in this case, usually there exists a monotonicity formula. As in the Bernstein problem for minimal hypersurfaces (see Fleming [53]), we can use the blowing down analysis and then the classification of homogeneous solutions to prove Liouville type results.

    This approach was first undertaken by the author in [86] to study the partial regularity of stable solutions to the Lane-Emden equation (see also Davila-Dupaigne-Farina [25] for related results)

    Δu=|u|p1u,p>1. (3.1)

    The scaling invariance for this equation says, if u is a solution of (3.1), then for any λ>0,

    uλ(x):=λ2p1u(λx)

    is also a solution of (3.1).

    The optimal Liouville theorem for stable and finite Morse index solutions to this equation was established in Farina [43] by a Moser type iteration argument. The blowing down analysis can be used to give another proof. By employing Pacard's monotonicity formula ([74,75]) and Federer's dimension reduction principle ([50]), a sharp dimension estimate on the singular set of stable solutions is given in [86]. In this approach, usually we need only an estimate up to the energy level.

    This approach was further developed in Davila et. al. [26] and Du-Guo-Wang [39]. In [26], a monotonicity formula is derived for the fourth order Lane-Emden equation

    Δ2u=|u|p1u. (3.2)

    Then by the blowing down analysis and the classification of stable, homogeneous solutions, an optimal Liouville theorem for stable solutions of (3.2) is established. In [39], a similar result is obtained for the weighted equation

    div(|x|θu)=|x||u|p1u. (3.3)

    For this equation, an ε-regularity theorem is needed to analyse the convergence of blowing down sequences. See also Dancer-Du-Guo [22], Du-Guo [37,38] and Wang-Ye [84] for related results on this equation.

    By now this approach has been applied to many other problems, for example, fourth order weighted equations or weighted systems [60,61], polyharmonic equations [70], nonlinear elliptic system [94], Toda system [89] and various elliptic equations involving fractional Laplacians [27,48,45,46,47,49,63].

    One may be tempted to believe that this approach works well once the equation enjoys a scaling invariance. However, there are several important exceptions.

    Problem 3. What is the optimal dimension for the Liouville theorem for stable solutions to the equation

    Δ2u=up,u0. (3.4)

    This equation arises from the MEMS problem, see Esposito-Ghoussoub-Guo [42]. For some p, a monotonicity formula was given in [26], but it is not known whether there exists a monotonicity formula for all p. See Guo-Wei [59], Huang-Ye-Zhou [62], Lai [67], Lai-Ye [68] for some recent results on this problem.

    We also encounter the same difficulty with the possible lack of a monotonicity formula in some other problems, which include the equation with p-Laplacians

    Δpu=|u|m1u (3.5)

    and its fourth order version

    Δ(|Δu|m1Δu)=|u|p1u. (3.6)

    When both u and Δu are positive, (3.6) is equivalent to the Lane-Emden system

    {Δu=vq,Δv=up, (3.7)

    see Mtiri-Ye [73]. For these problems, a mysterious problem is the role of homogeneous (or radial) solutions in the classification of stable solutions. In particular, is the radial, homogenous solution mostly unstable (in a suitable sense) among all solutions?

    Finally, if the blowing down analysis approach works, usually we could obtain a radial symmetry result about stable solutions when the space dimension is critical. By the moving plane method, this claim can be reduced to the classification of stable, homogeneous solutions in the critical dimension. This is similar to the classification of stable minimal hypercones in R8. For example, the following problem is still not completely solved.

    Problem 4. Suppose u is a stable, homogeneous solution of (3.1) in Rn, where n is the critical dimension. Is u radially symmetric (after a translation)?

    By the radial symmetry criterion of Guo [58], this is reduced to the classification of solutions wC2(Sn1) to the equation

    ΔSn1w+2p1(n22p1)w=|w|p1w, (3.8)

    satisfying the stability condition

    Sn1[|φ|2+(n2)24φ2]pSn1|w|p1φ2,φC1(Sn1). (3.9)

    If n is critical, we need to show that w must be constant.

    In [86], this is wrongly claimed to be true. But the proof therein works only for a small range p[pJL(n),pJL(n)+ε(n)). (Here pJL is the Joseph-Lundgren exponent, see [66].) The remaining case is still not known.

    In Dancer-Guo-Wei [24], infinitely many solutions to (3.8) are constructed. However, it seems difficult to verify (3.9) because it involves a spectral bound condition. To the best knowledge of the author, there is still no known smooth stable (in the sense of (3.9)) solutions of (3.8) other than the constant solutions.

    Problem 5. Take k,1 so that k+=n. Does there exist a nontrivial, stable solution of (3.1) invariant under the group SO(k)×SO()?

    Concerning finite Morse index solutions on Rn, for nonegative nonlinearities f0, Dupaigne and Farina in [41], under the same assumptions of Theorem 2.1 (plus some technical conditions), proved that finite Morse index solutions to (1.1) are radially symmetric.

    Next, for those scaling invariant equations discussed in Section 3, the blowing down analysis still works. So if the Liouville theorem holds for stable solutions, then it also holds for finite Morse index solutions, except in an exceptional dimension. (This is the "Sobolev" critical dimension. For example, for (3.1), it is well known that when p=(n+2)/(n2), there exist infinitely many solutions with finite Morse index.)

    In the rest of this section we review some recent results on finite Morse index solutions of Allen-Cahn equation (2.3). The author in [88], and later jointly with Wei in [90], studied the structure of finite Morse index solutions in R2. For example, the main result in [90] states that

    Theorem 4.1. A finite Morse index solutions of the Allen-Cahn equation (2.3) in R2 has finitely many ends.

    Here an end is a connected component of the nodal set {u=0} at infinity. But in fact more information such as the refined asymptotic behavior of u at infinity follows from the proof.

    By applying a reverse version of the infinite dimensional Lyapnunov-Schimdt reduction method (see Del Pino et. al. [30], [31], [32], [33], [34]), when the interfaces in the Allen-Cahn equation are clustering, we were able to reduce the stability condition in the Allen-Cahn equation to a corresponding one for Toda system

    Δfk=e2(fkfk1)e2(fk+1fk). (4.1)

    With such a connection between Allen-Cahn equation and Toda system, various results about stable solutions of Toda system can be transferred to the Allen-Cahn equation. For example, in Wang-Wei [90,91], Farina's integral estimate in [44] and the ε-regularity theorem in Wang [85,87], both for the Liouville equation

    Δu=eu, (4.2)

    were used to establish a curvature estimate for level sets of solutions to singularly perturbed Allen-Cahn equation. (More precisely, we need the corresponding results for Toda system (4.1), but which are direct generalizations of the results about Liouville equation, see [89].) In Gui-Wang-Wei [57], the Liouville theorem in Dancer-Farina [23] about finite Morse index solutions to (4.2) was used to establish the following result.

    Theorem 4.2. Suppose u is an entire solution of the Allen-Cahn equation (2.3) in Rn+1, and it is axially symmetric (i.e. u depends only on x21+x2n and xn+1) and stable outside a cylinder {x21++x2n<R2}. If 3n9, then u=u(xn+1).

    Here we do not state the n=2 case (i.e. for Allen-Cahn equation in R3), which can be viewed as the "Sobolev" critical case for Allen-Cahn equation. It was proved in [57] that in this case we have "finite Morse index implies finite ends" as in Theorem 4.1.

    The "finite Morse index implies finite ends" property is expected to be true in any dimension, but by now only this axially symmetric case has been proven. Note that, although Theorem 4.2 is a result in high dimensions, the axial symmetry makes the problem essentially two dimensional. This allows us to prove the quadratic decay for curvatures of {u=0} as in [90].

    In Cao-Shen-Zhu [11], it was shown that stable minimal hypersurfaces in Rn have only one end, that is, they are connected at infinity. Later in [69], Li-Wang also showed that minimal hypersurfaces with finite Morse index have finitely many ends (still in the topological sense). We want to know whether the corresponding results hold for stable or finite Morse index solutions of Allen-Cahn equation (2.3).

    Problem 6. Can we prove a topological finiteness result for ends of finite Morse index solutions to Allen-Cahn equation (2.3)?

    Since we only want a topological finiteness, this should hold for any dimension n. However, the one end result for stable solutions is not true for Allen-Cahn equation, at least in dimensions n11. In fact, the examples constructed in [1] are stable solutions of Allen-Cahn equation (2.3) in Rn, where n11, but they have two ends. Whether this one end result holds for n10 is open.

    In general, our understanding of finite Morse index solutions in higher dimensions is very lacking. Anyway, we do not know too much about stable solutions in dimensions n3. But even assuming the stable De Giorgi conjecture, we can not deduce that "finite Morse index implies finite ends", that is,

    Problem 7. Assuming the stable De Giorgi conjecture in dimensions 4n7, can we prove that finite Morse index solutions of (2.3) in these dimensions have finitely many ends?

    The three dimensional case can be proved as in [90] and [57], but this approach does not work in dimensions 4n7. This is because the techniques used in [90] relies on a finiteness result for nodal domains of directional derivatives ξu, which in turn needs the Liouville property for the degenerate elliptic equation (2.4), and as explained before, it cannot be applied in high dimensions. (In dimension 3, as in [57], a localization around each end allows us to use this technique, so we could expect the above result.) Note that the proof of Theorem 4.2 in [57] relies crucially on the axially symmetric assumption and the idea in Dancer-Farina [23]. Here the stability condition is used only in a reductive way as in Wang-Wei [90,91]. Therefore the power of the stability condition is not fully utilized. But now it is still not clear how to use the stability condition more efficiently in high dimensions.

    The author's interest in stable and finite Morse index solutions was largely intrigued by N. Dancer and Y. Du about ten years ago. They taught me a lot about stable and finite Morse index solutions when I was a postdoc at Sydney University. Several problems collected in this paper were communicated to the author during various occasions from Juan Davila, Louis Dupaigne, Zongming Guo, Xia Huang, Yong Liu, Yoshihiro Tonegawa, Juncheng Wei and Dong Ye over a long period. I am grateful to them for sharing their insights on these problems. My research has been supported by the National Natural Science Foundation of China (No. 11871381).



    [1] E. Wressle, C. Engstrand, A. K. Granérus, Living with Parkinson's disease: elderly patients' and relatives' perspective on daily living, Aust. Occup. Ther. J., 54 (2007), 131–139. https://doi.org/10.1111/j.1440-1630.2006.00610.x doi: 10.1111/j.1440-1630.2006.00610.x
    [2] P. Mahlknecht, A. Gasperi, P. Willeit, S. Kiechl, H. Stockner, J. Willeit, et al., Prodromal Parkinson's disease as defined per MDS research criteria in the general elderly community, Mov. Disord., 31 (2016), 1405–1408. https://doi.org/10.1002/mds.26674 doi: 10.1002/mds.26674
    [3] M. Politis, K. Wu, S. Molloy, P. G. Bain, K. R. Chaudhuri, P. Piccini, Parkinson's disease symptoms: the patient's perspective, Mov. Disord., 25 (2010), 1646–1651. https://doi.org/10.1002/mds.23135 doi: 10.1002/mds.23135
    [4] C. Váradi, Clinical features of Parkinson's disease: the evolution of critical symptoms, Biology, 9 (2020), 103. https://doi.org/10.3390/biology9050103 doi: 10.3390/biology9050103
    [5] D. J. Surmeier, Determinants of dopaminergic neuron loss in Parkinson's disease, FEBS J., 285 (2018), 3657–3668. https://doi.org/10.1111/febs.14607 doi: 10.1111/febs.14607
    [6] C. Raza, R. Anjum, Parkinson's disease: Mechanisms, translational models and management strategies, Life Sci., 226 (2019), 77–90. https://doi.org/10.1016/j.lfs.2019.03.057 doi: 10.1016/j.lfs.2019.03.057
    [7] A. B. Holt, E. Kormann, A. Gulberti, M. Pötter-Nerger, C. G. McNamara, H. Cagnan, et al., Phase-dependent suppression of beta oscillations in Parkinson's disease patients, J. Neurosci., 39 (2019), 1119–1134. https://doi.org/10.1523/JNEUROSCI.1913-18.2018 doi: 10.1523/JNEUROSCI.1913-18.2018
    [8] A. Singh, R. C. Cole, A. I. Espinoza, D. Brown, J. F. Cavanagh, N. S. Narayanana, Frontal theta and beta oscillations during lower-limb movement in Parkinson's disease, Clin. Neurophysiol., 131 (2020), 694-702. https://doi.org/10.1016/j.clinph.2019.12.399 doi: 10.1016/j.clinph.2019.12.399
    [9] M. H. Trager, M. M. Koop, A. Velisar, Z. Blumenfeld, J. S. Nikolau, E. J. Quinn, et al., Subthalamic beta oscillations are attenuated after withdrawal of chronic high frequency neurostimulation in Parkinson's disease, Neurobiol. Dis., 96 (2016), 22–30. https://doi.org/10.1016/j.nbd.2016.08.003 doi: 10.1016/j.nbd.2016.08.003
    [10] C. Hammond, H. Bergman, P. Brown, Pathological synchronization in Parkinson's disease: networks, models and treatments, Trends Neurosci., 30 (2007), 357–364. https://doi.org/10.1016/j.tins.2007.05.004 doi: 10.1016/j.tins.2007.05.004
    [11] Z. Wang, B. Hu, W. Zhou, M. Xu, D. Wang, Hopf bifurcation mechanism analysis in an improved cortex-basal ganglia network with distributed delays: An application to Parkinson's disease, Chaos, Solitons Fractals, 166 (2023), 113022. https://doi.org/10.1016/j.chaos.2022.113022 doi: 10.1016/j.chaos.2022.113022
    [12] B. Hu, X. Diao, H. Guo, et al., The beta oscillation conditions in a simplified basal ganglia network, Cogn Neurodynamics, 13(2019), 201-217. https://doi.org/10.1007/s11571-018-9514-0 doi: 10.1007/s11571-018-9514-0
    [13] A. B. Holt, T. I. Netoff, Origins and suppression of oscillations in a computational model of Parkinson's disease, J. Comput. Neurosci., 37 (2014), 505–521. https://doi.org/10.1007/s10827-014-0523-7 doi: 10.1007/s10827-014-0523-7
    [14] A. B. Holt, E. Kormann, A. Gulberti, M. Pötter-Nerger, C. G. McNamara, H. Cagnan, et al., Phase-dependent suppression of beta oscillations in Parkinson's disease patients, J. Neurosci., 39 (2019), 1119–1134. https://doi.org/10.1523/JNEUROSCI.1913-18.2018 doi: 10.1523/JNEUROSCI.1913-18.2018
    [15] L. L. Grado, M. D. Johnson, T. I. Netoff, Bayesian adaptive dual control of deep brain stimulation in a computational model of Parkinson's disease, PLoS Comput. Biol., 14 (2018), e1006606. https://doi.org/10.1371/journal.pcbi.1006606 doi: 10.1371/journal.pcbi.1006606
    [16] J. E. Fleming, E. Dunn, M. M. Lowery, Simulation of closed-loop deep brain stimulation control schemes for suppression of pathological beta oscillations in Parkinson's disease, Front. Neurosci., 14 (2020), 166. https://doi.org/10.3389/fnins.2020.00166 doi: 10.3389/fnins.2020.00166
    [17] A. B. Holt, D. Wilson, M. Shinn, J. Moehlis, T. I. Netoff, Phasic burst stimulation: a closed-loop approach to tuning deep brain stimulation parameters for Parkinson's disease, PLoS Comput. Biol., 12 (2016), e1005011. https://doi.org/10.1371/journal.pcbi.1005011 doi: 10.1371/journal.pcbi.1005011
    [18] K. Kumaravelu, D. T. Brocker, W. M. Grill, A biophysical model of the cortex-basal ganglia-thalamus network in the 6-OHDA lesioned rat model of Parkinson's disease, J. Comput. Neurosci., 40 (2016), 207–229. https://doi.org/10.1007/s10827-016-0593-9 doi: 10.1007/s10827-016-0593-9
    [19] A. J. N. Holgado, J. R. Terry, R. Bogacz, Conditions for the generation of beta oscillations in the subthalamic nucleus-globus pallidus network, J. Neurosci., 30 (2010), 12340–12352. https://doi.org/10.1523/JNEUROSCI.0817-10.2010 doi: 10.1523/JNEUROSCI.0817-10.2010
    [20] A. Pavlides, S. J. Hogan, R. Bogacz, Improved conditions for the generation of beta oscillations in the subthalamic nucleus-globus pallidus network, Eur. J. Neurosci., 36 (2012), 2229–2239. https://doi.org/10.1111/j.1460-9568.2012.08105.x doi: 10.1111/j.1460-9568.2012.08105.x
    [21] Gillies, D. Willshaw, Z. Li, Subthalamic-pallidal interactions are critical in determining normal and abnormal functioning of the basal ganglia, Proc. R. Soc. London, Ser. B, 269 (2002), 545–551. https://doi.org/10.1098/rspb.2001.1817 doi: 10.1098/rspb.2001.1817
    [22] J. E. Rubin, D. Terman, High frequency stimulation of the subthalamic nucleus eliminates pathological thalamic rhythmicity in a computational model, J. Comput. Neurosci., 16 (2004), 211–235. https://doi.org/10.1023/B:JCNS.0000025686.47117.67 doi: 10.1023/B:JCNS.0000025686.47117.67
    [23] D. Terman, J. E. Rubin, A. C. Yew, C. J. Wilson, Activity patterns in a model for the subthalamopallidal network of the basal ganglia, J. Neurosci., 22 (2002), 2963–2976. https://doi.org/10.1523/JNEUROSCI.22-07-02963.2002 doi: 10.1523/JNEUROSCI.22-07-02963.2002
    [24] B. Hu, M. Xu, L. Zhu, J. Lin, Z. Wang, D. Wang, et al., A bidirectional Hopf bifurcation analysis of Parkinson's oscillation in a simplified basal ganglia model, J. Theor. Biol., 536 (2022), 110979. https://doi.org/10.1016/j.jtbi.2021.110979 doi: 10.1016/j.jtbi.2021.110979
    [25] S. R. Cole, R. van der Meij, E. J. Peterson, C. de Hemptinne, P. A. Starr, B. Voytek, Nonsinusoidal beta oscillations reflect cortical pathophysiology in Parkinson's disease, J. Neurosci., 37 (2017), 4830–4840. https://doi.org/10.1523/JNEUROSCI.2208-16.2017 doi: 10.1523/JNEUROSCI.2208-16.2017
    [26] B. Pollok, V. Krause, W. Martsch, C. Wach, A. Schnitzler, M. Südmeyer, Motor‐cortical oscillations in early stages of Parkinson's disease, J.Physiol., 590 (2012), 3203–3212. https://doi.org/10.1113/jphysiol.2012.231316 doi: 10.1113/jphysiol.2012.231316
    [27] S. J. van Albada, P. A. Robinson, Mean-field modeling of the basal ganglia-thalamocortical system. I: Firing rates in healthy and parkinsonian states, J. Theor. Biol., 257 (2009), 642–663. https://doi.org/10.1016/j.jtbi.2008.12.018 doi: 10.1016/j.jtbi.2008.12.018
    [28] S. J. van Albada, R. T. Gray, P. M. Drysdale, P. A. Robinson, Mean-field modeling of the basal ganglia-thalamocortical system. Ⅱ: dynamics of parkinsonian oscillations, J. Theor. Biol., 257 (2009), 664–688. https://doi.org/10.1016/j.jtbi.2008.12.013 doi: 10.1016/j.jtbi.2008.12.013
    [29] G. W. Arbuthnott, M. Garcia-Munoz, Are the symptoms of parkinsonism cortical in origin?, Comput. Struct. Biotechnol. J., 15 (2017), 21–25. https://doi.org/10.1016/j.csbj.2016.10.006 doi: 10.1016/j.csbj.2016.10.006
    [30] C. F. Underwood, L. C. Parr-Brownlie, Primary motor cortex in Parkinson's disease: Functional changes and opportunities for neurostimulation, Neurobiol. Dis., 147 (2021), 105159. https://doi.org/10.1016/j.nbd.2020.105159 doi: 10.1016/j.nbd.2020.105159
    [31] M. D. Bevan, P. J. Magill, D. Terman, J. P. Bolam, C. J Wilson, Move to the rhythm: oscillations in the subthalamic nucleus-external globus pallidus network, Trends Neurosci., 25 (2002), 525–531. https://doi.org/10.1016/S0166-2236(02)02235-X doi: 10.1016/S0166-2236(02)02235-X
    [32] A. Pavlides, S. J. Hogan, R. Bogacz, Computational models describing possible mechanisms for generation of excessive beta oscillations in Parkinson's disease, PLoS Comput. Biol., 11 (2015), e1004609. https://doi.org/10.1371/journal.pcbi.1004609 doi: 10.1371/journal.pcbi.1004609
    [33] Y. Chen, J. Wang, Y. Kang, M. B. Ghori, Emergence of beta oscillations of a resonance model for Parkinson's disease, Neural Plast., 2020 (2020), 1–15. https://doi.org/10.1155/2020/8824760 doi: 10.1155/2020/8824760
    [34] M. M. McGregor, A. B. Nelson, Circuit mechanisms of Parkinson's disease, Neuron, 101 (2019), 1042–1056. https://doi.org/10.1016/j.neuron.2019.03.004 doi: 10.1016/j.neuron.2019.03.004
    [35] M. D. Humphries, J. A. Obeso, J. K. Dreyer, Insights into Parkinson's disease from computational models of the basal ganglia, J. Neurol., Neurosurg. Psychiatry, 89 (2018), 1181–1188. https://doi.org/10.1136/jnnp-2017-315922 doi: 10.1136/jnnp-2017-315922
    [36] B. C. M. van Wijk, H. Cagnan, V. Litvak, V. Litvak, A. A. Kühn, K. J. Friston, et al., Generic dynamic causal modelling: An illustrative application to Parkinson's disease, NeuroImage, 181 (2018), 818–830. https://doi.org/10.1016/j.neuroimage.2018.08.039 doi: 10.1016/j.neuroimage.2018.08.039
    [37] M. C. Chen, L. Ferrari, M. D. Sacchet, L. C. Foland-Ross, M. Qiu, I. H. Gotlib, et al., Identification of a direct GABA ergic pallidocortical pathway in rodents, Eur. J. Neurosci., 41 (2015), 748–759. https://doi.org/10.1111/ejn.12822 doi: 10.1111/ejn.12822
    [38] A. Saunders, I. A. Oldenburg, V. K. Berezovskii, C. A. Johnson, N. D. Kingery, H. L. Elliott, et al., A direct GABAergic output from the basal ganglia to frontal cortex, Nature, 521 (2015), 85–89. https://doi.org/10.1038/nature14179 doi: 10.1038/nature14179
    [39] P. R. Castillo, E. H. Middlebrooks, S. S. Grewal, L. Okromelidze, J. F. Meschia, A. Quinones-Hinojosa, et al., Globus pallidus externus deep brain stimulation treats insomnia in a patient with Parkinson disease, in Mayo Clinic Proceedings, Elsevier, 95 (2020), 419–422. https://doi.org/10.1016/j.mayocp.2019.11.020
    [40] J. Dong, S. Hawes, J. Wu, W. Le, H. Cai, Connectivity and functionality of the globus pallidus externa under normal conditions and Parkinson's disease, Front. Neural Circuits, 15 (2021), 8. https://doi.org/10.3389/fncir.2021.645287 doi: 10.3389/fncir.2021.645287
    [41] T. Tsuboi, M. Charbel, D. T. Peterside, M. Rana, A. Elkouzi, W. Deeb, et al., Pallidal connectivity profiling of stimulation-induced dyskinesia in Parkinson's disease, Mov. Disord., 36 (2021), 380–388. https://doi.org/10.1002/mds.28324 doi: 10.1002/mds.28324
    [42] R. G. Burciu, D. E. Vaillancourt, Imaging of motor cortex physiology in Parkinson's disease, Mov. Disord., 33 (2018), 1688–1699. https://doi.org/10.1002/mds.102 doi: 10.1002/mds.102
    [43] G. Foffani, J. A. Obeso, A cortical pathogenic theory of Parkinson's disease, Neuron, 99 (2018), 1116–1128. https://doi.org/10.1016/j.neuron.2018.07.028 doi: 10.1016/j.neuron.2018.07.028
    [44] A. Guerra, D. Colella, M. Giangrosso, A. Cannavacciuolo, G. Paparella, G. Fabbrini, et al., Driving motor cortex oscillations modulates bradykinesia in Parkinson's disease, Brain, 145 (2022), 224–236. https://doi.org/10.1093/brain/awab257 doi: 10.1093/brain/awab257
    [45] Z. Wang, B. Hu, L. Zhu, J, Lin, M. Xu, D. Wang, Hopf bifurcation analysis for Parkinson oscillation with heterogeneous delays: A theoretical derivation and simulation analysis, Commun. Nonlinear. Sci., 114 (2022), 106614. https://doi.org/10.1016/j.cnsns.2022.106614 doi: 10.1016/j.cnsns.2022.106614
    [46] T. P. Vogels, K. Rajan, L. F. Abbott, Neural network dynamics, Annu. Rev. Neurosci., 28 (2005), 357–376. https://doi.org/10.1146/annurev.neuro.28.061604.135637 doi: 10.1146/annurev.neuro.28.061604.135637
    [47] H. Kita, Y. Tachibana, A. Nambu, S. Chiken, Balance of monosynaptic excitatory and disynaptic inhibitory responses of the globus pallidus induced after stimulation of the subthalamic nucleus in the monkey, J. Neurosci., 25 (2005), 8611–8619. https://doi.org/10.1523/JNEUROSCI.1719-05.2005 doi: 10.1523/JNEUROSCI.1719-05.2005
    [48] J. T. Paz, J. M. Deniau, S. Charpier, Rhythmic bursting in the cortico-subthalamo-pallidal network during spontaneous genetically determined spike and wave discharges, J. Neurosci., 25 (2005), 2092–2101. https://doi.org/10.1523/JNEUROSCI.4689-04.2005 doi: 10.1523/JNEUROSCI.4689-04.2005
    [49] H. Kita, S. T. Kitai, Intracellular study of rat globus pallidus neurons: membrane properties and responses to neostriatal, subthalamic and nigral stimulation, Brain Res., 564 (1991), 296–305. https://doi.org/10.1016/0006-8993(91)91466-E doi: 10.1016/0006-8993(91)91466-E
    [50] M. A. Lebedev, S. P. Wise, Oscillations in the premotor cortex: single-unit activity from awake, behaving monkeys, Exp. Brain Res., 130 (2000), 195–215. https://doi.org/10.1007/s002210050022 doi: 10.1007/s002210050022
    [51] W. Schultz, R. Romo, Neuronal activity in the monkey striatum during the initiation of movements, Exp. Brain Res., 71 (1988), 431–436. https://doi.org/10.1007/BF00247503 doi: 10.1007/BF00247503
    [52] N. E. Hallworth, C. J. Wilson, M. D. Bevan, Apamin-sensitive small conductance calcium-activated potassium channels, through their selective coupling to voltage-gated calcium channels, are critical determinants of the precision, pace, and pattern of action potential generation in rat subthalamic nucleus neurons in vitro, J. Neurosci., 23 (2003), 7525–7542. https://doi.org/10.1523/JNEUROSCI.23-20-07525.2003 doi: 10.1523/JNEUROSCI.23-20-07525.2003
    [53] H. Kita, Globus pallidus external segment, Prog. Brain Res., 160 (2007), 111–133. https://doi.org/10.1016/S0079-6123(06)60007-1 doi: 10.1016/S0079-6123(06)60007-1
    [54] H. Kita, A. Nambu, K. Kaneda, Y. Tachibana, M. Takada, Role of ionotropic glutamatergic and GABAergic inputs on the firing activity of neurons in the external pallidum in awake monkeys, J. Neurophysiol., 92 (2004), 3069–3084. https://doi.org/10.1152/jn.00346.2004 doi: 10.1152/jn.00346.2004
    [55] H. Nakanishi, H. Kita, S. T. Kitai, Intracellular study of rat substantia nigra pars reticulata neurons in an in vitro slice preparation: electrical membrane properties and response characteristics to subthalamic stimulation, Brain Res., 437 (1987), 45–55. https://doi.org/10.1016/0006-8993(87)91525-3 doi: 10.1016/0006-8993(87)91525-3
    [56] K. Fujimoto, H. Kita, Response characteristics of subthalamic neurons to the stimulation of the sensorimotor cortex in the rat, Brain Res., 609 (1993), 185–192. https://doi.org/10.1016/0006-8993(93)90872-K doi: 10.1016/0006-8993(93)90872-K
    [57] A. Gillies, D. Willshaw, Membrane channel interactions underlying rat subthalamic projection neuron rhythmic and bursting activity, J. Neurophysiol., 95 (2006), 2352–2365. https://doi.org/10.1152/jn.00525.2005 doi: 10.1152/jn.00525.2005
    [58] H. Kita, S. T. Kitai, Intracellular study of rat globus pallidus neurons: membrane properties and responses to neostriatal, subthalamic and nigral stimulation, Brain Res., 564 (1991), 296–305. https://doi.org/10.1016/0006-8993(91)91466-E doi: 10.1016/0006-8993(91)91466-E
    [59] Y. Hirai, M. Morishima, F. Karube, Y. Kawaguchi, Specialized cortical subnetworks differentially connect frontal cortex to parahippocampal areas, J. Neurosci., 32 (2012), 1898–1913. https://doi.org/10.1523/JNEUROSCI.2810-11.2012 doi: 10.1523/JNEUROSCI.2810-11.2012
    [60] Y. H. Tanaka, Y. Tanaka, F. Fujiyama, T. Furuta, Y. Yanagawa, T. Kaneko, Local connections of layer 5 GABAergic interneurons to corticospinal neurons, Front. Neural Circuits, 5 (2011), 12. https://doi.org/10.3389/fncir.2011.00012 doi: 10.3389/fncir.2011.00012
    [61] M. H. Higgs, S. J. Slee, W. J. Spain, Diversity of gain modulation by noise in neocortical neurons: regulation by the slow afterhyperpolarization conductance, J. Neurosci., 26 (2006), 8787–8799. https://doi.org/10.1523/JNEUROSCI.1792-06.2006 doi: 10.1523/JNEUROSCI.1792-06.2006
    [62] A. L. Barth, J. F. A. Poulet, Experimental evidence for sparse firing in the neocortex, Trends Neurosci., 35 (2012), 345–355. https://doi.org/10.1016/j.tins.2012.03.008 doi: 10.1016/j.tins.2012.03.008
    [63] A. P. Prudnikov, Y. A. Brychkov, O. I. Marichev, Integrals and series: direct laplace transforms, Routledge, 2018. https://doi.org/10.1201/9780203750643 doi: 10.1201/9780203750643
    [64] K. Udupa, N. Bahl, Z. Ni, C. Gunraj, F. Mazzella, E. Moro, et al., Cortical plasticity induction by pairing subthalamic nucleus deep-brain stimulation and primary motor cortical transcranial magnetic stimulation in Parkinson's disease, J. Neurosci., 36 (2016), 396–404. https://doi.org/10.1523/JNEUROSCI.2499-15.2016 doi: 10.1523/JNEUROSCI.2499-15.2016
    [65] M. Dagan, T. Herman, R. Harrison, J. Zhou, N. Giladi, G. Ruffini, et al., Multitarget transcranial direct current stimulation for freezing of gait in Parkinson's disease, Mov. Disord., 33 (2018), 642–646. https://doi.org/10.1002/mds.27300 doi: 10.1002/mds.27300
    [66] E. Lattari, S. S. Costa, C. Campos, A. J. Oliveira, S. Machado, G. A. M. Neto, Can transcranial direct current stimulation on the dorsolateral prefrontal cortex improves balance and functional mobility in Parkinson's disease?, Neurosci. Lett., 636 (2017), 165–169. https://doi.org/10.1016/j.neulet.2016.11.019 doi: 10.1016/j.neulet.2016.11.019
  • This article has been cited by:

    1. Sk Irsad Ali, Anjan Das, Apoorva Agrawal, Shubharaj Mukherjee, Maudud Ahmed, P M G Nambissan, Samiran Mandal, Atis Chandra Mandal, Characterization, spectroscopic investigation of defects by positron annihilation, and possible application of synthesized PbO nanoparticles* , 2021, 30, 1674-1056, 026103, 10.1088/1674-1056/abd2a9
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2241) PDF downloads(110) Cited by(0)

Figures and Tables

Figures(22)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog