Small amplitude chimeras for coupled clocks

We report the arise of small amplitude chimera states in three coupled pendulum clocks suspended on an oscillating base. Two types of chimeras are identified and described by the character of the behaviour of particular units (which can be both regular or irregular). The regions of the appearance of the dynamical patterns are determined and the scenarios of their coexistence with typical synchronization states are discussed. We investigate the chimeras’ basins of attraction, showing that the arise of complex dynamics is not straightforward and highly depends on the system’s parameters and the initial conditions. The latter is confirmed by the probability analysis, exhibiting the rare character of the observed attractors. The scenarios of bifurcations between the chimeric patterns are studied and supported using the energy balance method, which allows to describe the changes of the energy flows between particular nodes of the system. The results presented in this paper confirm the ones obtained for the previous models, extending the analysis with an additional degree of freedom.


Introduction
The spontaneous coexistence of coherence and incoherence in networks of coupled systems [1] is one of the most surprising phenomena observed in the nonlinear dynamics discipline. The so-called chimera state [2], named after a hybrid creature from the Greek mythology, combines both regular and irregular types of patterns, showing that even though the coupling scheme applied to the system is homogeneous, the oscillators can split into many groups of different behaviours, uncovering unexpected scenarios of partial synchronization and chaos.
Chimeras have been observed in various dynamical models, just to mention heterogeneous networks [3,4], chemical oscillators [5], mechanical systems [6,7] or the delayed ones [8,9]. Patterns of coexistence are reported in neuronal problems [10][11][12][13][14] (including such classical systems as the FitzHugh-Nagumo [12] or the Hindmarsh-Rose [13] ones) and can be related to the complex brain dynamics [14]. In [15], the authors present a solvable model exhibiting chimeric behaviours, while in [16], the spectral properties of chimeras are discussed. Coexisting states arise natu-rally on the road from complete synchronization to pure spatial chaos, which has been studied in [17,18] but can be also related to the transient dynamics occurring in the system (see, e.g. [19] for details). In [20], Sethia and Sen discuss the existence criteria for chimera states, while in [21], the classification scheme of the chimeric behaviour is described. A general study on the phenomenon can be found in [22].
Apart from typical large networks of coupled oscillators, as the ones referenced above, coexisting patterns can be also found in small models of connected units. This particular scenario of chimeras' appearance has been introduced as the so-called weak chimeras [23,24], which occur even for the simplest case of three coupled oscillators [25,26]. Small/weak chimeras have been observed in Kuramoto model with inertia [26], experimental pendula [25] or delay-coupled lasers [27]. Some studies on their arise in systems with broken symmetry and complex, nonlinear interactions have been described in [28] and [29], respectively.
In this paper, we report the arise of small chimeras, investigating the amplitude of the oscillations of the connected units. The amplitude chimeras [30,31], characterized by the coexistence of (at least) two groups of oscillators synchronized within different amplitudes appear universally in networks of coupled systems, e.g. in chaotic oscillators [32] or non-locally coupled maps [33]. The stability of amplitude chimeras has been studied in [34], while in [35], Banerjee et al. investigate their relation with classical, phase patterns.
The results presented in this study extend the ones discussed in the previous works [25,26], introducing an additional degree of freedom, which corresponds to the dynamics of the supporting base. When the motion of the base is included, the scenario becomes more realistic (the oscillations of the support, on which the nodes are working cannot be excluded completely in practical implementations). Moreover, the relation between the motion of the oscillators and the base enables a new way of the energy distribution (the coupling scheme within the network becomes more complex). The inclusion of the support's dynamics (i.e. the central hub node globally coupled with all of the units [36,37]) can influence the behaviour of the oscillators and uncover new, complex states, just to mention the transient ones [38]. In this study, we focus on the smallest amplitude chimeras, which have not been described previously.
Namely, we investigate the dynamics of a simple system of three coupled pendulum clocks hanged on an oscillating base, which is schematically shown in Fig. 1. The supporting base mounted at the origin of 0x y plane (the location of the mass centre) can oscillate around the axis, which is perpendicular to Figure's plane (position denoted by θ variable). The properties of the base are described by the following parameters: B 0 [kgm 2 ] (moment of inertia in regard to the origin), k θ [Nm] (the stiffness of the spring connecting the base and the unmoving support) and c θ [Nms] (the damping of the damper connecting the base and the unmoving support). The support is connected with three identical planar pendulum clocks (marked by red, blue and green), which are mounted at points S i , i = 1, 2, 3 (shown as the brown dots). The suspension of the 1st pendulum (point S 1 ) is located above the origin, while the positions of the remaining two points (S 2 and S 3 ) are chosen to preserve the angular symmetry of the model, i.e. (0S i , 0S i+1 ) = 120 • (i = 1, 2, 3; index i periodic mod 3). The distance between the origin and each of the suspension points is the same and has been denoted by parameter r = 0S i (i = 1, 2, 3). The displacement of i-th pendulum is given by ϕ i variable, while the parameters of the nodes are as follows: m (the mass), l (the length) and c ϕ (the damping, not shown in Figure). The pendula shown in Fig. 1 are coupled via the supporting base, as well as the springs with stiffness coefficient k ϕ [N/m] (shown as the pink ones). Due to the number of the oscillators (three), the topol-ogy of the latter coupling scheme is both local and global [25,26]. In the considered scenario, the motion of the nodes induces the oscillations of the base and the pendula are not parametrically excited as in classical models [39][40][41][42].
The Lagrange equations of motion determining the dynamics of the model shown in Fig. 1 can be described in the following way: where i = 1, 2, 3.
During our research, we have fixed the following parameters of model (1) [Nms] has been varied during the research to preserve fixed logarithmic decrement at = ln(2) (in reference to a simple 1 DOF system of an oscillating base with the same values of the parameters as described above).
The motion of the pendula in system (1) is induced by the build-in escapement mechanisms (not shown in Fig. 1), providing the moment of force described by functions M E i , i = 1, 2, 3 [43,44]. The i-th function depends on both pendulum's displacement (ϕ i ) and the position of the i-th mechanism's cogwheel versus the mechanism's pallet (parameter σ i ∈ {1, 2}; the cogwheel can rest on either left or the right pallet). Each of functions M E i can be expressed as follows: where M = 0.075[Nm] denotes the value of the fixed external momentum, while ε = 5.0 • is the escapement's threshold (the mechanism is working as long as the displacement of the pendulum exceeds the threshold). Components V ϕ i and V θ in eqs. (1) describe the moments of forces acting on the pendula (and indirectly on the base), which are provided by the coupling springs (see Fig. 1 for details). The values of components V θ and V ϕ i depend on both s i (the distance between the i-th and the (i + 1)-th pendulum masses for the unmoving system, i.e. in the equilibrium state) andŝ i (the distance between the masses when the nodes oscillate) and can be expressed using the following formulas: In this study, we have varied both B 0 and k ϕ (the moment of inertia of the supporting base and the stiffness of the springs, respectively), uncovering possible dynamical patterns of model (1). Our results are discussed in Sect. 2.

Results
The regions of different types of dynamics of system (1) in two-parameters plane (B 0 , k ϕ ) are presented in Fig. 2, with corresponding examples of behaviours shown in Fig. 3. The description of the possible states is as follows.
The pendula of model (1) can completely synchronize (ϕ 1 = ϕ 2 = ϕ 3 ), which is one of the most typical behaviours observed in coupled mechanical systems. In this case, the synchronization pattern is stable in the whole region shown in Fig. 2. However, we have found that when the moment of inertia B 0 is small enough (region labeled 'CS' and shown in the light grey in Fig. 2), the synchronization state becomes the general solution of the system, within which all of the pendula oscillate (the case when the escapement mechanisms of the units are partially switched off has not been investigated in this study).
(the natural period of the pendula for chosen length l = 0.24849 [m]). As one can see, in this case the lower nodes have clustered (the 2nd-ϕ 2 (blue) and the 3rd-ϕ 3 (green)) and oscillate in the anti-phase to the upper one (the 1st-ϕ 1 (red)). The amplitude of the latter pendulum is greater (different) than the amplitude of the remaining ones.
In the considered case of three coupled oscillators, the smallest amplitude chimera can be defined as the state, in which all of the nodes are frequency synchro-nized (ω 1 = ω 2 = ω 3 , where ω i denotes the frequency of the oscillations of the i-th pendulum, i = 1, 2, 3), but only two of them are oscillating with the same amplitude (there exists a solitary node with different amplitude). For the case shown in Fig. 3a, the amplitudes i of the pendula satisfy the following relation: The yellow region shown in Fig. 2 is surrounded by the orange one, where another example of chimera has been found, i.e. the Quasiperiodic Amplitude Chimera, labeled 'QAC' in Fig. 2. The latter pattern is shown in Fig. 3b In this case, the behaviour of the system becomes quasiperiodic, which can be observed in the Poincare sections in the inbox in the left panel in Fig. 3b. The motion of the base has been chosen as the reference for the creation of the maps, i.e. whenθ = 0 and θ > 0 (which means that the base has reached its maximum), both position ϕ i and velocityφ i of each pendulum are marked in the map. In this case, the dynamics of the lower pendula, i.e. the 2nd (blue) and the 3rd (green) is almost the same (see the curves in the Poincare section) but the attractors are shifted in space, forming an anti-phase pattern. On the other hand, the upper, red pendulum is moving almost periodically (the small irregularities are caused by the quasiperiodic motion of the lower nodes; see the phase-time plot in the right panel in Fig. 3b), and its amplitude is qualitatively different (smaller) that the amplitude of the remaining two oscillators. As in the previous case (the PAC in Fig. 3a), one can identify two coexisting groups of patterns (the 1st node versus the 2nd and the 3rd one), which can be distinguished by the difference in the amplitudes ( 1 < 2 = 3 ) and lead to the classification of the amplitude chimera.
Bifurcating system (1) (increasing parameter B 0 ), we have observed that the Quasiperiodic Amplitude Chimera (QAC) resides not only in the orange domain in Fig. 2 but can be also found in the yellow one, where the Periodic Amplitude Chimera (PAC) is present. The yellow region of the coexistence of both types of chimeric patterns (PAC and QAC) has been hatched in Fig. 2. However, with the decrease in the moment of inertia B 0 (or the springs' stiffness k ϕ ), the QAC state vanishes and the model bifurcates into the phase-locked synchronization profile. The latter solution appears in the dark grey domain labeled 'PLS' in Fig. 2 and is shown in Fig. 3(c) (parameters: B 0 = 3.045 [kgm 2 ], k ϕ = 11.75 [N/m]). In this scenario, each pendulum oscillates with the same frequency and different amplitudes (the motion is periodic), while the phases between the nodes are locked (see the phase-time plot in the right panel in Fig. 3c for details).
With further decrease in coefficients B 0 , k ϕ , the behaviour of the pendula loses periodicity and becomes quasiperiodic. The irregular domain is marked in Fig. 2 in the brown colour and can be observed in a wide region of the parameters values. The example of the dynamics corresponding to this scenario is presented in Fig. 3d (B 0 = 0.585 [kgm 2 ] and k ϕ = 9.35 [N/m]), where one can see that the motion of the pendula becomes more complex than in the case of the QAC state (this can be especially observed in the inbox in the left panel in Fig. 3d, where the Poincare maps for the oscillators are included). Due to the increase in the irregularity, which affects the dynamics of the 1st (red) pendulum (the oscillations of the node are more irregular, compared to an almost periodic motion discussed in Fig. 3b), we have not identified the states similar to the one presented in Fig. 3d with any kind of chimera. However, when the coupling strength k ϕ is small enough, one can observe an arise of patterns qualitatively similar to the one presented in Fig. 3b, i.e. characterized by an almost periodic motion of the upper pendulum (the 1st) and an anti-phase profile between the lower ones (the 2nd and the 3rd). The surprising appearance of the QAC states is marked in Fig. 2 as the brown, hatched region at the bottom of the map. The latter result uncovers that the additional connection of the oscillators through the springs (see Fig. 1 for details) is not necessary for the occurrence of the Quasiperiodic Amplitude Chimeras, i.e. the states are also possible when k ϕ = 0. However, the local coupling is still required for the arise of the Periodic Amplitude Chimeras (PAC), as shown in the map in Fig. 2.
Dynamical patterns described in Fig. 2 reside in different regions of parameters plane (B 0 , k ϕ ). To investigate the possible scenarios of their coexistence, we have calculated the basins of attraction for system (1), which are discussed in Fig. 4.
The phase space of system (1) is an 8-dimensional subset of the R 8 space. To present the cross sections of the basins of attraction (on the plane), we have considered the following two subsets of the initial conditions: Scenario I C L corresponds to the case, when the motion of both lower pendula (the 2nd and the 3rd) is initiated from the same conditions (ϕ 0 2 = ϕ 0 3 ), which may support the appearance of the PAC state (in this pattern, the lower pendula are synchronized). On the other hand, I C R case involves the scenario, when the escapement mechanism of the 3rd pendulum (green in Fig. 2) is switched off (ϕ 0 3 = 0) and the motion has to be induced by the remaining two oscillators (the red and the blue ones, as shown in Fig. 2).
The methodology of the initial conditions selection has a major influence on the dynamics of the system, which has been shown further. The basins of attraction corresponding to scenarios I C L and I C R (calculated in (ϕ 1 , ϕ 2 ) ∈ [−π, π) 2 cross section) are presented in the left and the right panel in Fig. 4, respectively.
In Fig. 4a, one can observe the results obtained for parameters B 0 = 5.5 [kgm 2 ] and k ϕ = 18.0 [N/m]. In this case, three dynamical patterns are possible, i.e. the Complete Synchronization (CS), the Periodic Amplitude Chimeras (PAC) and the Quasiperiodic Amplitude Chimeras (QAC), which are marked in the light grey, yellow and orange, respectively (see the legend in the bottom in Fig. 4). When the initial conditions are chosen from I C L set (left panel in Fig. 4a), one can see the basin of the PAS state (yellow), surrounded by the narrow, orange stripes corresponding to the QAS solution. However, when I C R scenario is considered (right panel in Fig. 4a), the middle of the phase-space is filled mainly by the irregular orange basin. The small area, where the yellow domain can be found (the hatched inbox) has been shown in the enlargement. In both cases, the basin of the complete synchronization (CSlight grey) is the dominant one.
A similar scenario but with the absence of the PAC state can be observed in Fig. 4b for B 0 = 4.8 [kgm 2 ] and k ϕ = 9.0 [N/m]. The orange basin (QAC) competes with the light grey one (CS), but the regime of the second state is dominating the phase space. It should be noted that the small, black region marked in the right panel in Fig. 4b corresponds to the scenario, when the escapement mechanism of the upper pendulum (the 1st one) has switched off and the dynamics is mainly induced by the oscillations of the lower nodes and the supporting base.
The dominant character of the CS attractor changes, when the general quasiperiodic dynamics occurs. The example of such a scenario is presented in Fig. 4c (parameters: where only the CS (light grey) and the Q (brown) dynamical patterns coexist. In this case, the measures of the regimes occupied by the basins are similar (in contrast to the results described in Fig. 4a and 4b) and their structure becomes more complex. The riddled nature of the domains (especially around ϕ 1 , ϕ 2 ≈ ±π ) implicates that the final behaviour of system (1) is very hard to predict and may change even for a slight perturbation of the initial conditions or parameters. The results discussed in Fig. 4a exhibit that the choice of the initial conditions (sets I C L and I C R ) has a major influence on the appearance of the PAC and the QAC states (especially the former one, when I C R scenario is applied). To investigate the probabilities of the chimeras' arise, we have performed a statistical analysis, which is discussed in Fig. 5.
The results presented in  Fig. 5a and b, respectively. Each iteration point (marked by dots in Figure) has been examined for a total number of 2000 trials (the initial conditions randomly chosen from the uniform distribution using I C L or I C R subset).
As one can see, when scenario I C L is considered (ϕ 0 2 = ϕ 0 3 ; Fig. 5a, the probability of the occurrence of the PAC state is greater than the QAC one almost in whole B 0 interval (except the smallest possible values around B 0 ≈ 1.8 [kgm 2 ]). With the increase in the moment of inertia, both probabilities also increase, reaching the stable levels around p = 0.11 (yellow) and p = 0.02 (orange) for B 0 ≈ 20 [kgm 2 ]. On the other hand, the appearance of the complete synchronization pattern (CS -light grey) decreases from p = 0.99 for The second scenario, i.e. I C R one (ϕ 0 3 = 0) is shown in Fig. 5b and uncovers that the probability of the occurrence of the QAC state (orange) exceeds the one observed for the PAC pattern (yellow) for B 0 < 10.9 [kgm 2 ]. Moreover, the Periodic Amplitude Chimeras have not been observed during the analysis for B 0 < 5.7 [kgm 2 ]. The latter phenomenon is caused by the change of the initial conditions procedure. The arise of the PAC state is more probable, when the initial conditions for the lower pendula are the same, i.e. ϕ 0 2 = ϕ 0 3 (as in scenario I C L ), since the state involves the same motion of the lower nodes (see Fig. 3a for details); when one of the pendulum is switched off at the beginning (ϕ 0 3 = 0; case I C R ), the periodic chimeric pattern may not arise from the chosen initial conditions and parameters at all. With the increase in B 0 coefficient, the probabilities of the states stabilize around p = 0.09 (yellow), p = 0.005 (orange) and p = 0.905 (light grey) for B 0 = 20 [kgm 2 ].
We have found that with further increase in the moment of inertia B 0 , the results do not change significantly, oscillating around the values obtained for B 0 = 20 [kgm 2 ]. The dominant character of the CS solution presented in Fig. 5 indicates that the basins of attraction for the parameters chosen from the orange and the yellow regions in Fig. 2 (where amplitude chimeras are possible) should be qualitatively similar to the ones presented in Fig. 4a and b (the value of B 0 parameter has a minor influence on the probability of the appearance of the complete synchronization attractor). Taking into account the obtained results, chimeras discussed above can be identified as the rare attractors [45][46][47][48] and their occurrence in the system from randomly chosen initial conditions is highly uncertain.
To study the transitions between different patterns appearing in system (1), we have performed the bifur-cation analysis, applying the energy balance method [43,44,49,50].
Multiplying the dynamical equation of the i-th pendulum in system (1) by its velocity, i.e.φ i and integrating for one full period of oscillations, one can get the energy equation, describing the flows of the energy components for the node: The left-hand side of eq. (5) describes the increase in the total energy of the i-th pendulum during one full period of oscillations and equals zero in the case of periodic motion. The components on the right-hand side of Eq. (5) can be denoted in the following way: where i = 1, 2, 3. Component (6) W i ESC [J] describes the energy supplied to the i-th pendulum by its escapement mechanism, while (7)  The results of the bifurcation and the energy analysis for model (1) are discussed in Fig. 6.
We investigate two bifurcation scenarios, i.e. (i) the transition from the PAC (Periodic Amplitude Chimera) to the QAC (Quasiperiodic Amplitude Chimera)left panel in Fig. 6 and (ii) the transition from the PLS (Phase-Locked Synchronization) to the Q state (Quasiperiodic motion) -right panel in Fig. 6. In both cases, the moment of inertia is fixed at B 0 = 4.005 [kgm 2 ] and the stiffness of the coupling springs is decreased in the following intervals: (i) k ϕ ∈ [11.0, 11.5] [N/m]; (ii) k ϕ ∈ [8.0, 9.6] [N/m] (see the horizontal axes in Fig. 6). The bifurcation diagrams are presented in Fig. 6a, where we denote the local maxima of the supporting base (θ ) and the pendula (φ i ) during the oscillations (colour-code as in Fig. 3). On the other hand, the energy diagrams for the 1st, the 2nd and the 3rd pendulum are shown in Figs. 6b, c and d, respectively. The colour code used to distinguish different energy components is as follows: The transition from the PAC to the QAC states is discussed in the left panel in Fig. 6. Beginning from k ϕ = 11.5 [N/m] and decreasing the parameter, one can observe that the periodic amplitude chimera (see Fig. 3a) is present for the coupling strength k ϕ > 11.3[N/m]. The higher amplitude of the 1st node (compared to pendula 2nd and 3rd) can be explained by the energy flows. Indeed, the energy dissipated by the 1st oscillator is higher than the energy supplied by its escapement mechanism: W 1 DAMP > W 1 ESC (Fig. 6b). In order to preserve the same frequency of the oscillations for all of the units, this lack of energy has to be compensated by the decreased amplitudes of the remaining nodes (the 2nd and the 3rd). This is possible due to the fact that the escapement mechanisms of the latter units produce more energy than is required by their dampers: W i DAMP < W i ESC , i = 2, 3 ( Fig. 6c  and d). The first pendulum consumes the energy from the system: W 1 FLOW < 0, while the other two oscillators supply it: W 2 FLOW > 0, W 3 FLOW > 0. However, this scenario begins to change around k ϕ = 11.3[N/m], when the amplitude of the 1st pendulum decreases (see the red curve in Fig. 6a). Dissipation W 1 DAMP also decreases and the consumption of the energy is smaller (W 1 FLOW increases). Simultaneously, the amplitudes of the lower pendula change, i.e. the amplitude of the 2nd node decreases, while the amplitude of the 3rd one increases (blue and green curves in Fig. 6a, respectively), which corresponds to the changes in the energy   Figs. 6c and d uncover the differences in the directions of the energy distribution for the pen-dula: the energy transferred via the base by the 2nd (the 3rd) oscillator decreases (increases), while the flow transferred via the coupling springs by the 2nd (the 3rd) node increases (decreases). The variation of the amplitudes continues until the occurrence of the sudden bifurcation to the QAC state around k ϕ = 11.18 [N/m]. Interestingly, the bifurcation arises close to the point in which W 3 DAMP = W 3 FLOW (the energy dissipated by the 3rd pendulum collides with the energy transferred to the system). With further decrease in coupling strength k ϕ , the Quasiperiodic Amplitude Chimera can be observed, until the state bifurcates to the PLS pattern (Phase-Locked Synchronization). Due to the fact that the motion is no longer periodic in this scenario, the analysis of the energy balance cannot be longer followed.
The scenario presented in the right panel in Fig. 6 corresponds to the case, when the phase-locked synchronization (the PLS state) leads to the appearance of the quasiperiodic motion (Q). The analysis begins at k ϕ = 9.36 [N/m], when the motion becomes periodic (the system leaves the region of the QAC). In the first stage (k ϕ ∈ [9.17, 9.36] [N/m]), the amplitude of the 3rd node is greater than the amplitude of the 2nd one (see the green and the blue curves in Fig. 6a). Indeed, the additional energy produced by the 2nd pendulum (W 2 DAMP < W 2 ESC ) can be transferred to the 3rd neighbour, for which W 3 DAMP > W 3 ESC (Fig. 6c and  d), respectively). Simultaneously, the amplitudes vary (blue increases, green decreases) and finally get even around k ϕ = 9.17 [N/m] (W 2 FLOW = W 3 FLOW = 0; see the hatched, horizontal lines in Fig. 6c and d).
With further decrease in parameter k ϕ < 9.17 [N/m], the energy scenario reverses. The amplitude of the 2nd node (blue) exceeds the amplitude of the 3rd one (green), which is accompanied by the following inequalities: W 2 DAMP > W 2 ESC , W 2 FLOW < 0 and W 3 DAMP < W 3 ESC , W 3 FLOW > 0. The distribution of components W i BASE , W i SPRING , i = 2, 3 shows the directions in which the energy is transferred. It should be noted that the whole process (during the PLS state) is supported by the 1st pendulum (red in Fig. 6a), which has the smallest amplitude and supplies the energy to the system: W 1 DAMP < W 1 ESC , W 1 FLOW > 0. With further decrease in coupling strength k ϕ , the periodic motion vanishes (around k ϕ = 8.3 [N/m]) and the behaviour becomes irregular. The latter bifurcation occurs around W 3 DAMP = W 3 FLOW point, similarly to the case described previously (left panel in Fig. 6; see the details above).
The energy scenarios discussed in Fig. 6 seem universal and have been observed also for other values of parameters B 0 , k ϕ during the study on system (1).

Conclusions
In this paper, we have studied the dynamics of three coupled pendulum clocks suspended on an oscillating supporting base. The inclusion of the support (an additional degree of freedom) expands the research performed before for small networks of coupled oscillators and allows to observe new types of solutions. We have distinguished two types of amplitude chimeras, characterized by the periodic and the quasiperiodic dynamics of the nodes. The regions of the occurrence of the dynamical patterns uncover that the states can coexist in different configurations, including the appearance of such classical states as the complete or the phase-locked synchronization. Typical examples of behaviours have been shown and discussed. Investigations on the basins of attraction (realized in various scenarios, depending on the methodology of the initial conditions selection) have shown that the chimeric dynamics may be hard to predict and can appear in the system spontaneously. The latter observations have been studied closely during the probability analysis, where the rare character of the patterns has been confirmed. We have shown that the investigations on the transitions between different types of attractors can be supported by the energy balance method, which follows the qualitative changes occurring in the model. Our results participate in the research on chimera states in small networks of coupled oscillators, confirming the universal appearance of complex, coexisting patterns in dynamical systems. statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/ by/4.0/.