Dynamics of oscillators globally coupled via two mean fields

Many studies of synchronization properties of coupled oscillators, based on the classical Kuramoto approach, focus on ensembles coupled via a mean field. Here we introduce a setup of Kuramoto-type phase oscillators coupled via two mean fields. We derive stability properties of the incoherent state and find traveling wave solutions with different locking patterns; stability properties of these waves are found numerically. Mostly nontrivial states appear when the two fields compete, i.e. one tends to synchronize oscillators while the other one desynchronizes them. Here we identify normal branches which bifurcate from the incoherent state in a usual way, and anomalous branches, appearance of which cannot be described as a bifurcation. Furthermore, hybrid branches combining properties of both are described. In the situations where no stable traveling wave exists, modulated quasiperiodic in time dynamics is observed. Our results indicate that a competition between two coupling channels can lead to a complex system behavior, providing a potential generalized framework for understanding of complex phenomena in natural oscillatory systems.

the neurons are actually coupled by two mean fields. Interaction through different channels is also characteristic for physiological problems, where, e.g., cardiac and respiratory systems show coexisting couplings [16][17][18][19] . Another situation with two mean-field couplings, is related to the attempts to control synchrony in a population of coupled oscillators 20,21 . Such an approach has been e.g. discussed in the context of suppression of collective brain oscillations at Parkinson's disease [22][23][24] . Here the two channels of coupling are the internal (uncontrollable) one, and  the external coupling due to an imposed feedback. Both couplings can be considered to a good approximation as global ones, acting on the whole population of the involved neurons. Potentially, the most simple experimental realization of the two mean field coupling would be an extension of a recent experimental setup where Kuramoto-Sakaguchi coupling scheme with linear and nonlinear couplings have been experimentally verified 25 . The scheme for the two-field coupling is presented in Fig. 1 of the supplementary material. Instead of one common resistive load, like in the experiment 25 , one could implement two such loads, and additionally modify the phase shifts of actions of the second mean field on the oscillators. This would exactly correspond to the particular model we study theoretically below.
We focus below on the simplest possible setup, where all oscillators that are described in the phase approximation have the same natural frequency; furthermore, the coupling is assumed to involve the first harmonics of the oscillations only, like in the standard Kuramoto-Sakaguchi model. We show that mostly nontrivial dynamics of this system is observed if the two mean fields act differently, i.e. one is mainly attractive while another one is repulsive. We find uniformly rotating (traveling wave) solutions and study their stability in the thermodynamic limit, where also the Ott-Antonsen approach can be applied 26 . In the case traveling waves are unstable, modulated traveling waves are observed. The findings are supported by direct numerical simulations of finite ensembles.

Results
Model formulation. In this paper we consider a simple model of a population of phase oscillators subject to a coupling through two mean fields. We start by formulating the problem in a rather general context, but will make several simplifications to achieve a tractable but still nontrivial model. We assume that all the phase oscillators have the same natural frequency, and differ only by the way how they contribute to the mean fields and how they are forced by them. Furthermore, we assume that these differences are only in the phases of the coupling, not in the amplitudes. Physically, this can be modeled by an ensemble of acoustic oscillators, sounds of which are collected by two microphones, and which are subject to the forcing emitted by two loudspeakers. If the positions of the oscillators are different, they contribute to the mean fields with different phase shifts and get signals that are also differently phase shifted. On the other hand, attenuation of signals can be neglected, thus only the phase relations are important.
We denote the phases of the oscillators φ k , k = 1, …, N, and define two complex mean fields Y (1,2) according to (1,2) where N is the number of oscillators in the system, and γ j (1,2) are the phase shifts with which the oscillators contribute to the mean fields. The dynamics of the phases, driven by these fields, is given by equations (written in the reference frame rotating with the common frequency) where ε 1,2 are the coupling constants of the two fields, and δ k (1,2) are the phase shifts with which the fields act on oscillators.
Generally, the model above would be fully defined if the joint distribution density of the phase shifts γ δ W( , ) (1,2) (1,2) is given. To simplify, we assume that the parameters γ δ , (1,2) (1,2) are independent on each other in the populations. Due to this, as one can see from (2), the dynamics of φ k does not depend on the phase shifts γ (1,2) and one can simplify the expressions for the mean fields (1) as 0 5, as functions of Θ 2 for several values of λ. Large markers: states with all P l = 1 (all oscillators locked); medium markers: states that are partially locked < < P 0 1 l ; small markers: all oscillators are unlocked P l = 0. Arrows at the x-axis show the linear stability boundaries for λ = 0.1 and for λ = 0.2, according to Eq. (16).
Scientific RepoRts | 7: 2104 | DOI:10.1038/s41598-017-02283-1 i i (1,2) ( 1,2) (1,2) (1,2) where we introduced the usual Kuramoto complex mean field and two complex constants characterizing the distributions of the phase shifts γ (1,2) : The dynamics of the phases (2) remains depending on the distributions of δ k (1,2) . Below we consider a minimal nontrivial case, when only one set of the phase shifts δ k (1,2) has a nontrivial distribution (say, the second one), while another set consists of equal shifts δ (1) . Then we can rewrite Eq. (2) as We renormalize time ε → t w t 2 (2) and obtain finally the basic model that we will study in this paper:

2
(2) describing the relative strengths of couplings of the two mean fields; ( 1) is the effective constant phase shift for the first mean field; Θ 2 is the phase at which the distribution of the phase shift of the second field µ δ − k ( 2) has a maximum; and finally the parameter 2 describes the deviation of the phase shift from this most probable one. Below we use the von Mises distribution of phase shifts α 0 characterized by the parameter Δ: Δ = 0 corresponds to a uniform distribution of the phase shifts, while in the limit Δ → ∞ one gets a delta-function.
It is instructive to rewrite Eq. (6) in the real form: The physical interpretation of this two-mean-field-coupling model (8) is as follows: each of the couplings is of Sakaguchi-Kuramoto type, i.e. it contains the sin term only. The two couplings have different phase shifts, for the first coupling it is fixed to Θ 1 , for the second coupling the phase shift is different for different oscillators, but is concentrated around Θ 2 . Parameter λ defines the relative weight of the two couplings. The main parameters of the problem are λ Θ Θ , , 1 2 , and Δ. One can expect mostly nontrivial effects if the two couplings act in opposite directions: one tries to synchronize the oscillators, while another one is repulsive and tends to desynchronize them. Below we will keep all the parameters in the theoretical considerations, but in numerical examples we will mainly set either Θ = 0 1 , that corresponds to a purely attractive first coupling, or Θ = 0 2 , that corresponds to a purely attractive second coupling.

Stability of incoherent state. First, we consider stability properties of the completely incoherent state
where the phases of all the oscillators are uniformly distributed and the mean field vanishes. Analytical expression for the growth rate of potentially unstable perturbations, derived in the section Methods, reads , this boundary is depicted in Fig. 1(a). One can see that stability of the incoherent state is only possible if the second coupling is repulsive (i.e. Θ 2 is close to π) and relatively strong (λ is small). In panel (b) we show the case when Θ = 0 2 is fixed, here the stability boundary of the incoherent state is given by . In this figure there is another nontrivial line to be discussed below in the next section.
Traveling wave solutions and their stability. Here we discuss nontrivial regimes of partial synchrony in the model. We use the self-consistent approach to find the solutions; reformulate the system in terms of Ott-Antonsen (OA) equations; and then determine stability of the solutions by calculating their stability spectra. The details of these methods are described in the Methods Section, here we summarize the results.
The traveling wave synchronous solution have the mean field with a constant amplitude R and a uniformly rotating phase Ψ = Ω + Ψ t 0 . The distribution of the phases is stationary in the reference frame rotating with Ω. Driven by the two mean fields, some oscillators are locked (i.e. they rotate with the same frequency Ω), while others are not locked and rotate (although non-uniformly) with some other frequencies. The main parameters we use to characterize the state of the system, are the amplitude R and the locking parameter P l (cf. Eq. (25)), determining which portion of oscillators in the population is locked (P l = 1 means all are locked; P l = 0 means no one is locked). For determination of stability, we use a reformulation of the dynamics in terms of the Ott-Antonsen integro-differential system of equations (see description in Methods Section, in particular Eq. (26)); the latter equation is linearized and discretized to find the spectrum of eigenvalues.
In Fig. 2 we present some typical stationary synchronous states together with their stability spectra. In these figures we visualize the stationary density via a dispersion of a finite set of points, for better visibility. One can easily distinguish domains of locked (all points representing oscillators collapse to a line) and rotating (scattered) oscillators. We must emphasize that the partial phase locked state here is not a chimera state, as the oscillators react to the fields with different phase shifts, and therefore there is no structural symmetry in the system which has to be broken in the chimera phenomenon. Panels (a), (c), (e) and (g) show the profile of normal branches which can also be found in one mean field coupling cases, while the anomalous branches shown in panel (b), (d), (f) and (h), do not appear in the traditional one mean field case. One remarkable thing is that, the ratio of phase locked oscillators in (b) is larger than in (a), but the order parameter in (b) is smaller than in (a). This is also a specific effect of the two mean field coupling, compared to the usual one mean field coupling.
One can see from Fig. 2 that the spectrum generally consists of a continuous part and several discrete eigenvalues. The continuous part is related to existence of the branch of the locked oscillators (a set of purely real eigenvalues) and to the unlocked oscillators (a set of purely imaginary eigenvalues). In the case both branches are present (panels (a,d) of Fig. 2), one has a characteristic T-shaped continuous (essential) spectrum, as argued in refs 27 and 28 for the chimera states, where also locked and unlocked oscillators are present. The state in panel (b) of Fig. 2 is fully locked, here only the real continuous spectrum is present. In the case (c) of Fig. 2, where all the oscillators are unlocked, the purely imaginary continuous spectrum is observed.
As argued in refs 27 and 28 stability is determined by the discrete spectrum, which is clearly seen for all cases: panels (a,d) depict stable solutions, while panels (b,c) depict unstable ones. In some cases (not illustrated here) it is difficult to distinguish the discrete spectrum from the continuous one, as the discrete eigenvalues have very small real part and are "smeared" (i.e. depend significantly on the offset parameter α 0 ). Here refining the resolution (i.e. increasing the number of discretization points L) helps, but the computation time increases rapidly.
Normal, anomalous, and hybrid states. Basing on the approach described above, we have determined the uniformly rotating states for different values of parameters, and characterized their stability. Here below we describe the main types of the solutions.
Attractive first coupling. In this subsection we present the results for the case Θ = 0 1 , i.e. the first (non-distributed) coupling is purely attractive. We present the dependencies of the order parameter R on Θ 2 , obtained using the self-consistent method above, in Fig. 3. To fix Θ = 0 1 , in relations (21) we varied parameter β in the range β π ≤ < 0 2 and found all the points with = Θ = Q 0 1 . These points yield a parametric representation of the order parameter R and of Θ 2 as functions of the remaining auxiliary parameter a. Additionally, we calculated the portion of locked oscillators according to (25) and coded it in Fig. 3 with the size of the markers into three types: all oscillators locked (P l = 1), all oscillators unlocked (P l = 0), and partial locking ( < < P 0 1 l ). First, one can notice the symmetry of the diagram Fig. 3 Θ → −Θ 2 2 , this is the consequence of the selected value Θ = 0 1 . Therefore, for the analysis of the dependence of the order parameter R on λ, presented in Fig. 4, we choose only negative values of the phase shift Θ 2 . In Fig. 4 we also use the style of the lines (solid, dashed, dotted) to distinguish fully locked, partially locked, and fully unlocked states.
The situation is quite simple for Θ ≈ 0 2 : here both mean fields act attractively and the dependence on the parameter λ is monotonous (see the curve for π Θ = − . 0 4 2 in Fig. 4(a)). Remarkably, here solutions with small value of λ, where the second mean field dominates, are fully locked, while at medium values of λ, where both fields compete, only partial locking is observed. At large λ the attractive first coupling dominates again, and here P l = 1.
For phase shifts Θ 2 that are closer to π − , instead of a monotonous dependence of the order parameter R on λ, we obtain a dependence characteristic for a subcritical transition of the first-order type, see the curve for π Θ = − . 0 6 2 in Fig. 4(b). Stability consideration reveals that in the region where two nontrivial solutions exist, solutions with lower value of R are unstable and those with higher value of R are stable. Thus, we have a typical diagram for a hysteresic transition of the first order: there is a range of parameter λ where the incoherent and coherent solutions are both stable, while "in between" there is an unstable solution. When this unstable solution meets the incoherent one (see arrows below the Θ 2 -axis in Fig. 3) the latter loses its stability. It is worth mentioning that the large part of the unstable coherent solution is fully unlocked, although a part of it can be partially locked (e.g., this happens for π Θ = − . 0 9 2 , see Fig. 4(c)). The situation for Θ 2 even closer to −π, illustrated by the results for Θ 2 = −0.9π in Fig. 4(c), is rather unusual and deserves a more detailed description. First, one can see here the branch which is very much similar to the solution for Θ 2 = −0.6π, only it is shifted to larger values of parameter λ. Here the stability properties are similar to that at Θ 2 = −0.6π: solutions with λ > dR d / 0 are stable while the state with λ < dR d / 0 is unstable. The latter unstable part of this branch appears exactly at the linear stability border point. We call this branch "normal branch", as it is very similar to the situation observed at a usual subcritical transition.
One can see that additionally there exists another branch of solutions with relatively small values of R. Noteworthy, the values of λ at which this branch emerges from R = 0, have nothing to do with the change of linear stability of the incoherent state. Therefore we call this branch "anomalous". The anomalous solutions can also be seen in Fig. 3, they do not end at the points (marked by arrows) where linear stability of the incoherent state changes.
To clarify the behavior of the solutions in the limit R → 0, we need to take this limit for the solutions constructed by Self-consistent approach presented in the Methods Section. There are two possibilities to have vanishing order parameter R: one is related to vanishing of r in the equation for the phases (19), and another is related to the case Ω = 0. The former, "normal" case corresponds to a uniform distribution of the phases for all α in the limit R → 0, i.e. to the incoherent state. The values of parameters at which this happens are exactly those at which the linear stability of the incoherent state changes. Thus the "normal" branch bifurcates from the non-coherent state.
In contradistinction, the states with vanishing R that correspond to Ω = 0, are not close to the incoherent state with a uniform distribution of phases. Moreover, as Fig. 4(c) shows, one part of the "anomalous" branch consists of states with P l = 1, i.e. all the oscillators are locked. Therefore we cannot describe the appearance of the "anomalous" branch as a bifurcation from the incoherent solution.
Stability analysis of the anomalous branch shows that the solutions with P l = 1 are unstable, while another part with < < P 0 1 l consists of unstable and stable solutions. The change of stability happens without appearance of new traveling wave solutions, although we cannot exclude that some other solutions (e.g., modulated traveling waves) may appear at these points.
Attractive second coupling. Here we discuss a situation where the phase shift of the second coupling is fixed Θ = 0 2 , while the phase shift of the first coupling Θ 1 varies. Figures 5 and 6 are similar to the Figs 3 and 4 in the previous subsection. Figure 5 shows nontrivial states when the first and the second couplings act in different directions, i.e. when Θ 1 is close to π. Let us focus on the dependence of the order parameter R on the coupling constant λ in Fig. 6. For small values of Θ 1 , the incoherent state is unstable for all λ, and a traveling wave solution exists -however it is unstable for intermediate values of λ, see Fig. 6(a). For larger values of π Θ > /2 1 , where the first coupling is repulsive, the incoherent state becomes stable at large λ (see Fig. 6(b), here π Θ = . 0 6 1 ). When λ decreases, at the instability border the normal branch of solutions softly appears with small values of R. With decreasing of λ this branch makes a loop, and only for small λ a part of this branch becomes stable. For larger values of Θ 1 , the loop becomes larger and eventually crosses the axis R = 0 at two points (which correspond to the situations where anomalous solutions have vanishing order parameter, these points are marked red in Fig. 1(b)); see also panels (c,d) in Fig. 6. Remarkably, as a result the branch which starts as a normal branch, now ends as an anomalous one, thus it can be termed "hybrid branch". Another anomalous branch that consists of solutions with P l = 1, ends at λ = 0. The anomalous and the hybrid branches exist also for π Θ = . Simulation of the Ott-Antonsen system. In this and the next section we report on the numerical tests of the found solutions. The first test deals with the system in the thermodynamic limit. Here, according to the Scientific RepoRts | 7: 2104 | DOI:10.1038/s41598-017-02283-1 Ott-Antonsen ansatz 26 , one can write for each α a closed equation for the local order parameter ψ α e i and thus to represent the whole system as an integro-differential equation (eq. (26) in the Method Section). The simulation is accomplished via discretization of the integral and by solving the resulting finite-dimensional system with the Runge-Kutta method. In the cases, when at least one stable traveling wave solution is present, direct numerical simulations confirm the stability results above. Thus we present only the nontrivial cases where no stable traveling wave solution exists (see Fig. 6 above). In Fig. 7 we show the patterns appearing in these situations. The dynamics of the order parameter α z t ( , ) is periodic in time. Together with rotation of the phase this gives quasiperiodic dynamics (see the snapshots of the modulated synchronous states below where the same patterns are checked in finite-size ensemble simulations), i.e. a modulated traveling wave. Remarkably, patterns in panels (a,b) of Fig. 7 show two locked and two unlocked regions. Modulated wave in panel (c) of Fig. 7 can be characterized as a "blinking locked state", as here for large time intervals the local order parameter is close to one everywhere, and only for relatively short time intervals, weakly unlocked regions, where this order parameter is reduced to ≈ . z 0 9, are observed.

Direct numerical simulations of oscillator population.
In this section, we present results for numerical simulations of a finite size oscillator system, Eq. (8). First, we present results for a large system: N = 10000.   Fig. 4, and also confirm our previous results. Figure 9 shows the simulation results for traveling wave states depicted in Fig. 6 above. The parameters are Θ = 0 2 , ∆ = . π . 0 85 . The simulation results shows the normal, the anomalous, and the hybrid branches in Fig. 6, and also confirm our stability analysis. In the cases where no stable  Fig. 10 the details of the system behavior for cases presented in Fig. 7. Panels (a), (d) and (g) show periodic in time evolution of the global order parameter R. Panels (b), (e) and (h) demonstrate quasiperiodic rotations of the complex order parameter, confirming that the observed regimes are modulated traveling waves. In panel (c), (f) and (i), one can see that the oscillators are partial locked. In fact, in panel (i) one cannot recognize the unlocked pattern, because at the moment of time when the snapshot has been performed, the local order parameter is everywhere close to one. Only during small time epochs over the period the unlocked region is visible (see discussion of "blinking unlocking" above.) Next, we study systems with a small number oscillators. Here we focus mainly on situations of multistability, to check whether all the states stable in the thermodynamic limit can be also observed for small populations. For each λ, we plot the final order parameters generated from 1500 random initial conditions. Comparing with the results for a large ensemble N = 10000 (Fig. 8), one can see that for the cases of π Θ = − . 2 , a small system shows the same behavior as a large one. However, for the case of π Θ = − . 0 9 2 , we miss the anomalous solution in the small system, only the normal branch is observed. We conclude that the basin of this anomalous solution is rather small, so possibly finite-size fluctuations in a small system lead to a transition to the normal branch. Figure 4 in the Supplementary Information shows the probability of a finite-time system to end up in the incoherent state, calculated for 50 realizations with randomly chosen initial conditions. Parameters are Θ = 0 1 , ∆ = .
0 5 and λ = . 0 2 (see the green curve in Fig. 3 for the solutions in the thermodynamic limit).  (16). One can see that with the increasing of system size, the probability of getting incoherent state are closer to 1 and the stability boundaries are also closer to the theoretical results. For system size ensembles, the probability of getting incoherent state becomes lower. This can be explained by the fluctuations of the order parameter, which increase with the decreasing system size, so the system can jump out of the basin of the incoherent state due to the fluctuations. The nontrivial traveling wave state is relatively stable for these parameters and survives finite-size fluctuations.

Discussion
In this paper, we studied the dynamics of a population of oscillators driven by two mean fields via two different channels. We have focused on the simplest setup, where all the oscillators have identical frequencies, and the actions of the two mean fields differ by the coupling strengths and by the distributions of coupling phase shifts. The most close situation, previously discussed in the literature, is that of the standard globally coupled ensemble of oscillators (e.g., in the Kuramoto-Sakaguchi formulation) with an additional external feedback [20][21][22][23][24] . In these studies, however, it has been assumed that both channels act with constant, non-distributed coupling phases, what allowed for an essential reduction to a single-channel coupling. Another situation close to the considered in this paper, is that of star-coupled oscillators (i.e. all oscillators are coupled to one central, which then represents one channel of global coupling) with an additional global mean field 29 .
We have presented a general theory of traveling wave solutions in the thermodynamic limit, including study of the stability in the Ott-Antonsen approximation. Mostly nontrivial states are those where the two mean fields act oppositely, i.e. one tends to synchronize the oscillators and another desynchronizes them. Here we have found different types of solutions: states where all oscillators are locked, states where all oscillators are unlocked, and combined states where a part of oscillators are locked by driving fields, and another part is unlocked (we stress that these states cannot be called chimeras, as here there is no symmetry in the population, because the fields act on different oscillators with different phase shifts).
Stability of the found states has been studied using the Ott-Antonsen approach. The problem can be reduced to a linear integral equation, which after discretization can be formulated as a matrix eigenvalue problem. We have demonstrated existence of the essential continuous spectrum, two parts of which lie on the imaginary axis (for unlocked oscillators) and on the real axis in the stable domain (for locked oscillators). Stability is determined by the discrete spectrum which in most cases can be found rather reliably.
By comparing amplitude behavior and stability of the found traveling wave solutions with the stability properties of the incoherent state, we identified normal and anomalous branches of nontrivial solutions. Normal branches bifurcate from the incoherent state when the latter becomes unstable, typically in a subcritical way. The corresponding transitions are hysteretic and discontinuous, which is similar to recent findings of the discontinuous transition in other synchronization setups [30][31][32][33][34][35] . Additionally, we have found that some branches of solutions appear with vanishing order parameter at points where no change of stability of the incoherent state occurs. These anomalous branches thus do not bifurcate from the incoherent state. Some parts of the anomalous branches may be stable, and they are observed in the simulations of the Ott-Antonsen equations describing the thermodynamic limit. However, in simulations of small finite ensembles the population "prefers" more robust normal stable branches. Furthermore, we have found hybrid branches, which combine properties of the normal and anomalous ones: such a branch starts from a vanishing order parameter at a point where the incoherent state becomes unstable, but ends at a point where stability of the incoherent state does not change. Additionally, we have found situations where no stable traveling waves exist. Here numerical simulations, both in the thermodynamic limit and in finite ensembles, demonstrated modulated traveling waves with periodically in time pulsating order parameter. Remarkably, in some cases such a wave looks like blinking unlocking: there are epochs when all oscillators look like being locked, while in other epochs one can clearly see a wide distribution of the phases.
Actually, in the situation which the oscillators are coupled by three or more mean fields, and all the channels act with constant. The system can be reduce to our framework by summing the attractive/repulsive couplings into resultant attractive/repulsive mean field, and then it can be studied base on our results. But if the coupling channels have more complicated forms, and can not be summed into resultant mean fields, then one needs to develop a specific approach for the situation.
In conclusion, the competition between two coupling mean fields can make the system behavior much more complicated than the case with only one channel of global coupling. Further study in this direction would help to understand nontrivial synchronization regimes in complex systems in nature, such as coupled cellular systems and neuronal systems.

Methods
Linear stability of the incoherent state. In this section we perform a linear stability analysis of the fully incoherent (i.e. with a uniform distribution of phases) state of system (8) in the thermodynamic limit of infinite number of oscillators. In this limit one starts with a continuity equation for the density function ρ(θ, t, α) with the velocity υ 1 2 and the mean field   Self-consistent approach. We seek for partially synchronous solutions in the model described by Eq. (6) in the form of uniformly rotating states, with some frequency Ω, to be defined in the procedure. In terms of the distribution density ρ ϕ α t ( , , ) these solutions are traveling waves. It is convenient to transform all the variables to the rotating frame, where we then will look for stationary solutions of the distribution of the phases. Furthermore, because we treat the problem in the thermodynamic limit, it is suitable to parametrize the phases by the value of parameter α (we will write it as a lower index for the phase variables). We thus introduce new phases ψ θ = − Ψ − Θ 1 and assume that Ψ = Ω  . In these variables the ensemble driven by two mean fields (8)

∫ ∫
One can see that the essential parameters for the dynamics of oscillators in (19) are r and s, which depend explicitly on a, λ, β, α. This suggests the following strategy to find the solutions of the model:  Because the found solutions are steady states of system (26), the linearization is straightforward. Furthermore, we discretize the resulting linear integral equations using a discrete (size L) representation of the parameter α on the circle: α α π = + k L 2 / k 0 , = … − k L 0, , 1 . Here the offset α π ≤ < L 0 2 / 0 is a free parameter. It can be used to distinguish the continuous and the discrete parts of the spectrum (cf. ref. 37).
Solutions with a vanishing order parameter. Here we analyze the limit R → 0 of the constructed above traveling wave solutions. Because, according to (21), R = Ωa, there are two possibilities for order parameter R to vanish: a = 0 and Ω = 0.
For a → 0 we have also r → 0. Eq. (23) in the limit r → 0 reduces to