Abstract

Spatially uniform cross-correlated sine-Wiener (CCSW) noises are imposed on each Hodgkin-Huxley (HH) neuron in a two-dimensional (2D) regular network. Noise-induced spiral wave destruction and order-disorder spatial pattern transition can be observed by adjusting cross-correlation time or cross-correlation intensity of CCSW noises. The sudden change of the curve of synchronization factor and the time series of average membrane voltage F can be used to semiquantitatively assess this spatial pattern transition and spiral wave destruction induced by CCSW noises, respectively. In addition, moderate cross-correlation time and strong cross-correlation intensity of CCSW noises are both detrimental to survival of spiral wave. Comparing with non-cross-correlated sine-Wiener noises, CCSW noises can destroy organized spiral wave with lesser noise amplitude.

1. Introduction

A hot issue in neurosciences is to investigate synchronized collective behaviors of a group of interacting neurons due to their potential applications in information processing, information transmission, and so on. Synchronized collective behaviors of neurons in a network often manifest spatiotemporal pattern. Spiral wave is one of the most important spatiotemporal patterns and widely observed in various real systems. For example, spiral wave is observed in the cardiac tissue and it relates to a kind of heart disease [13]. Therefore, many schemes have been presented to remove and suppress spiral wave in the cardiac tissue [46]. Additionally, spiral wave in the neocortex may relate to seizure by extending the duration of evoked activity in pathological conditions [7, 8]. However, spiral wave in the neocortex can be used as a spatial framework to coordinate cortical oscillations over a group of neurons so that it may contribute to signal communication [7, 8]. In consideration of importance of spiral wave, it is especially interesting and important to investigate the dynamics of spiral waves in various different systems.

Random fluctuations (also known as noises) are omnipresent in nature [9]. Choosing reasonable model to mimic these random fluctuations is very important for the study of the dynamics of nonlinear systems [10]. Gaussian noise is widely used to model random fluctuations because of its convenience of analysis. However, considering that some real physical quantities, such as reaction rates, must be strictly positive particularly attention must be paid in order to avoid the fact that some undesired and unrealistic model-based results may occur due to unboundedness of Gaussian noise [11]. Unlike unbounded Gaussian noise, bounded noise takes values in bounded intervals and bounded noise is naturally consistent with the boundedness of real physical quantities in real systems [10, 11]. Furthermore, bounded noise is a convenient and natural mathematical tool for modeling either narrow band or broad band random process by appropriately adjusting corresponding parameters [12]. Hence, bounded stochastic processes are increasingly employed in various research fields involving electrical engineering, mechanical and structural engineering, and biological systems [13]. For example, impacts of bounded noise on the formation and instability of spiral waves in regular and small-world networks of Hodgkin-Huxley (HH) neurons were investigated in our previous works [14, 15] and the constructive role of bounded noise in facilitating formation and stability of spiral waves was discovered [14, 15].

It is noteworthy that the majority of studies about the impact of bounded noises including our previous works assume that a unique noise or multiple uncorrelated noises are perturbing an otherwise deterministic system. However, there may be correlation between multiple noises if they share a common origin. A few studies on cross-correlated bounded noises have been reported recently. For example, cross-correlated sine-Wiener (CCSW) noises can induce transitions in tumor-immune system interplay [16], genotype selection model [17], and the modified FitzHugh–Nagumo (FHN) neuron model [18]. More interestingly, CCSW noises-induced transition vanishes if both noises are not correlated [18]. Moreover, CCSW noises-induced coherence resonance (CR) phenomenon is discovered in the FHN neuron. Furthermore, CR is sensitive to the amplitudes and correlation times of CCSW noises [19]. CCSW noises [20] can also greatly enhance weak signal detection and transmission through resonance mechanism like periodic force [2123], time delay [24, 25], electromagnetic induction [2631], and sine-Wiener noise [32]. Although some studies on CCSW noises have been achieved, to our best knowledge, the effects of CCSW noises on the spiral wave and pattern dynamics have not been so far reported. To do this, we firstly construct a two-dimensional (2D) regular network of HH neurons subjected to spatially uniform CCSW noises and then investigate spiral wave destruction and spatial pattern transition induced by spatially uniform CCSW noises.

2. Model and Simulation

As reported previously, CCSW noises of and are described by following formulas [16, 18, 33]:Here and are amplitudes of and , while and are self-correlation times of and , respectively. Unless other specified, and are set to 15 throughout this work. and are two cross-correlated standard Wiener processes with cross-correlation functions defined asHere denotes correlation intensity. The symbol of represents an ensemble average and means smaller value between and . For the same correlation times (), the cross-correlation function of and can be expressed as follows [16]:withHere    and denote cross-correlation intensity and cross-correlation time of and , respectively. It is very difficult to directly treat CCSW noises due to the presence of coupling between noises. As reported previously, an existing decoupling procedure is used to equivalently transform (1) into the following formulas [16]:Herein and denote two independent standard Wiener processes [34]. They can be numerically simulated by the Euler-Maruyama algorithm in conjunction with the Box-Muller algorithm to generate normally distributed variables [16, 35]:Here , , , and are four independent random numbers uniformly distributed between 0 and 1.

The mathematical model for a 2D regular network of HH neurons subjected to CCSW noises is given by the following coupled differential equations [36, 37]:Here denotes the voltage of the cellular membrane of the neuron at the node , while , , and stand for gate channel parameters of the neuron at the node . All parameter values used throughout this paper are listed in Table 1.

As reported previously, to measure spatial synchronization, a synchronization factor is introduced as follows [14, 15, 38]:Here represents averaging over time. The value of indicates the degree of synchronization. Larger denotes better synchronization. In particular, that is close to 0 means no synchronization, while that is close to 1 indicates perfect synchronization. As reported previously, the sudden change of can be used to indicate synchronization transition [14, 39]. In the presence of CCSW noises, a sufficiently large time window of 500 time units is used in the calculation of synchronization factor .

According to Euler forward difference algorithm, equations of (7) and (8) are solved with the no-flux boundary condition and time step in our numerical simulations. A wedge-shaped initial value method of (41 : 43, 1 : 50) = −40.2,  (44 : 46, 1 : 50) = 0,  (47 : 49, 1 : 50) = 40.0,  (41 : 43, 1 : 50) = 0.1203,  (44 : 46, 1 : 50) = 0.5203,  (47 : 49, 1 : 50) = 0.98203,  (41 : 43, 1 : 50) = 0.9, (44 : 46, 1 : 50) = 0.7,  (47 : 49, 1 : 50) = 0.5,  (41 : 43, 1 : 50) = 0.9,  (44 : 46, 1 : 50) = 0.7,  and (47 : 49, 1 : 50) = 0.5 is used to induce a spiral seed, while the initial values for other neurons are selected with ,  ,  , and  . In the absence of CCSW noises, a sufficiently large time window of 500 time units is used in order to ensure obtaining a perfect spiral wave. CCSW noises are introduced into deterministic system at a given time instant of time units, and then CCSW noises are persistently imposed on each HH neuron in the system.

3. Main Results

This study focuses primarily on spiral wave destruction and spatial pattern transition induced by spatially uniform CCSW noises. Therefore, the effect of cross-correlation time of CCSW noises is firstly investigated. As shown in Figure 1, for small cross-correlation time (such as , 0, and 0.5), spiral wave can survive under the disturbance of CCSW noises (Figures 1(a) and 1(b)). However, with the increase of cross-correlation time , spiral wave starts to be destroyed (Figure 1(c)). If cross-correlation time is large enough, CCSW noises can completely destroy spiral wave and result in the occurrence of irregular pattern (Figures 1(d)-1(e)). Interestingly, if cross-correlation time of CCSW noises is further increased, spiral wave actually keeps its alive (Figures 1(h)-1(i)). Therefore, moderate cross-correlation time of CCSW noises is detrimental to survival of spiral wave. Moreover, CCSW noises-induced order-disorder spatial pattern transitions from spiral wave to irregular pattern can be obtained by adjusting cross-correlation time of CCSW noises. Besides qualitative observation from snapshots of spatial pattern, this transition process can be semiquantitatively assessed by the sudden change of the curve of . As displayed in Figure 2, that sharply increases at first and then decreases rapidly as cross-correlation time increases results in an obvious peak in the curve of -, which suggests spatial pattern transition from spiral wave to irregular pattern (Figure 2).

It is also interesting and important to observe dynamic process of CCSW noises-induced spiral wave destruction. The snapshots in Figure 3 show the breakup of spiral wave under the disturbance of CCSW noises. Without CCSW noises (Figure 3(a)), a perfect spiral wave can be observed in the network (Figure 3(a)). In the presence of CCSW noises, spiral wave starts to be destroyed by CCSW noises (Figure 3). Finally, spiral wave is almost completely destroyed and irregular spatial pattern appears in the network under the persistent disturbance of CCSW noises (Figure 3). Similarly, in addition to qualitative observation from snapshots of spatial pattern, spiral wave destruction can be identified by analyzing the time series of average membrane voltage . Before CCSW noises are applied to deterministic FHN model, stable spiral wave (indicated by stable periodical oscillation) can be observed in the network over a period of about 200 time units from running (Figure 4). After CCSW noises are imposed on each neuron, spatial pattern is under disturbance of CCSW noises and spiral wave may be destroyed to some extent. For small cross-correlation time of CCSW noises, time series of presents quasi-periodical oscillation, which indicates that spiral wave can survive (Figure 4). For large enough cross-correlation of CCSW noises, time series of become arbitrary, which denotes CCSW noises-induced spiral wave destruction (Figure 4).

Besides cross-correlation time, cross-correlation intensity of CCSW noises is a key factor to significantly influence spatial pattern. Therefore, it is also interesting and important to investigate the effect of cross-correlation intensity of CCSW noises. For weak cross-correlation intensity of CCSW noises, spiral wave can keep it alive under noise disturbance (Figures 5(a) and 5(b)). With the increase of cross-correlation intensity of CCSW noises, CCSW noises may destroy regular spiral wave and cause formation of irregular pattern (Figure 5(c)). Interestingly, further increasing cross-correlation intensity of CCSW noises does not always destroy spiral wave. For some certain moderate , spiral wave can survive under the noise disturbance (Figures 5(d) and 5(e)). Of course, spiral wave can be completely destroyed for a strong enough cross-correlation intensity (Figures 5(f)5(i)). In other words, strong cross-correlation intensity of CCSW noises is detrimental to survival of spiral wave. In brief, the results illustrated in Figure 5 show noise-induced spiral wave destruction and pattern transition from regular spiral wave to irregular pattern due to the increase of cross-correlation intensity of CCSW noises. As described above, synchronization factor can be used to semiquantificationally assess noise-induced spatial pattern transition. As shown in Figure 6, small denotes the survival of spiral wave in the condition of weak cross-correlation intensity , while large indicates the breakup of spiral wave and the formation of irregular pattern (Figure 6).

Similarly, we also consider evolution process of spiral wave under the disturbance of CCSW noises. Spiral wave can be clearly observed in the initial transient process after CCSW noises are applied to network (Figures 7(a) and 7(b)). As time goes on, spiral wave pattern is gradually destroyed and irregular pattern is gradually generated (Figures 7(c) and 7(d)). Finally, spiral wave disappears (Figures 7(e) and 7(f)). The corresponding plot of is drawn in Figure 8. Time series of presents stable and regular oscillation if without CCSW noises in system, while time series of starts to become arbitrary if with CCSW noises in system, which suggests spiral wave is being destroyed by CCSW noises (Figure 8).

To get a global view, the contour plot of in the - plane is depicted in Figure 9. As indicated by the blue area with smaller value, too short or long cross-correlation time of CCSW noises is not fatal for spiral wave. On the whole, for moderate cross-correlation time , increases gradually along with the increase of cross-correlation intensity of CCSW noises (Figure 9). However, when or is fixed, the change of is not always monotonous, which indicates the occurrence of CCSW noises-induced spatial transition (Figure 9). As shown in Figure 9, strong cross-correlation intensity between noises results in large , which suggests strong cross-correlation intensity is detrimental to the survival of spiral wave. This conclusion can be further confirmed by comparing CCSW noises with non-cross-correlated sine-Wiener (UNCCSW) noises. Dependence of on increasing noise amplitudes and is plotted in Figure 10 for CCSW and UNCCSW noises. is assumed to be equal to in Figure 10 for the sake of the convenient comparison. The value of for CCSW noises is larger than that for UNCCSW noises (Figure 10). The dividing line between regular spiral wave and irregular pattern is near 0.00085. For CCSW noises, the threshold of noise amplitude for destroying the organized spiral wave is distinctly less than that for UNCCSW noises (Figure 10). That is to say, cross-correlation between CCSW noises may induce instability of spiral wave and formation of irregular pattern.

In what follows, we investigate how noise amplitudes of CCSW noises affect spatial pattern. Synchronization factor is calculated for different noise amplitude when is set to 15. When noise amplitude is small, the value of is smaller, which denotes that spiral wave can survive under the disturbance of CCSW noises (Figure 11(a)). When noise amplitude is large enough, the value of is larger, which indicates that spiral wave can be destroyed by CCSW noises (Figure 11(a)). However, the change of with increasing noise amplitude is nonmonotonic (Figure 11(a)). Similarly, a contour plot of in the - plane is drawn to provide a global view. The blue area with smaller appears in the lower left corner of the contour plot of , which indicates the survival of spiral wave. The area with larger in the upper right corner of the contour plot of indicates disappearance of spiral wave and appearance of irregular pattern due to too large noise disturbance (Figure 11(b)). As a whole, increases gradually along with the increase of noise amplitudes and of CCSW noises (Figure 11(b)), but the increase of is not monotonous with the increase of or (Figure 11(b)). Although the two cross-correlated bounded noises are in the symmetric form as shown in formulas of (1), the contour plot is not strictly symmetrical along the diagonal, which can be explained by the asymmetry of amplitudes and in the formulas of (5). The asymmetry is introduced when we indirectly and equivalently deal with CCSW noises by an existing decoupling method.

4. Discussions and Conclusions

In summary, spatially uniform CCSW noises are persistently imposed on each HH neuron in 2D regular network after spiral wave is developed in the network by a wedge-shaped initial value scheme. Then the effects of CCSW noises on spiral wave destruction and spatial pattern transition are investigated. The obtained main results are concluded below. Spatial pattern transitions from spiral wave to irregular pattern can be observed by adjusting cross-correlation time or cross-correlation intensity of CCSW noises and this CCSW noises-induced spatial pattern transition can be semiquantitatively assessed by the sudden change of the curve of . The regular and irregular oscillation of time series of average membrane voltage correspond, respectively, to regular spiral wave and irregular pattern. Therefore, CCSW noises-induced spiral wave destruction can be indicated by time series of statistical parameter . For small or large cross-correlation time of CCSW noises, spiral wave may survive. For moderate cross-correlation time of CCSW noises, spiral wave may be destroyed. Strong cross-correlation intensity between CCSW noises is detrimental to the survival of spiral wave. In particular, comparing with non-cross-correlated sine-Wiener noises, CCSW noises can destroy spiral wave with lesser noise amplitude.

The network model in this work is a most widely adopted one. However, some limitations should be pointed out. The first one is that the coupling between neurons is bidirectional. In fact, neuronal networks are mainly coupled by chemical synapses, which is one directional coupling from presynaptic to postsynaptic neuron. The second one is that temporal delay is negligible in this study. In reality, time delay is inherent in realistic neuronal network due to finite propagation velocities of information. There is a time delay of ~1 ms between the presynaptic action potential and postsynaptic depolarization. However, in cardiac tissue the coupling is mainly gap junctions, the coupling is bidirectional, and there is less time delay. The third one is that topology structure of network is not considered in our model. Therefore, the network model we used here is a generalized and idealized mathematical model for neuronal network. We believe that idealization procedure of model does not affect the main results we report here. In the future work we will extend our current research to a more complex network model that includes more details such as time delay and complex network topology.

We only considered the temporal correlation between bounded noises. CCSW noises are assumed to be spatially uniform in 2D regular network, which implies a strong unrealistic assumption: the correlation length of the noises is infinite. It is more realistic and reasonable to introduce spatiotemporal Gaussian or bounded noises when investigating dynamics of spatially extended systems subjected to random disturbance. There have been several landmark studies about spatiotemporal noises [11, 4046]. In particular, spatiotemporal sine-Wiener bounded noise was proposed and its effect on Ginzburg–Landau model was investigated [43]. The emergence of noise-induced reentrant transitions was revealed through numerical simulations [43]. Furthermore, there have been some representative works that apply spatiotemporal bounded noises to a spatial biological problem. For instance, cell polarization may be destroyed when external spatiotemporal bounded noise is introduced into model [44].

In this work, we used the sudden change of the curve of to semiquantitatively assess pattern transition. More detailed quantitative comparison could be obtained by using the bivariate Moran index [42, 47]. Note that since and are two statistics, more information could be inferred by plotting, at a given time instant and for various values of the parameters of interest, the histograms of the whole distribution of the fields as in [42, 43]. The procedure allows better understanding of how the parameters of a spatiotemporal noise impact on the whole patterns of the output of the perturbed system. In future studies we will extend our current research from CCSW noises to spatiotemporal bounded noises and adopt more detailed quantitative methods for better presentation of research findings.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 31601071, 11605014, and 11675060) and the Huazhong Agricultural University Scientific & Technological Self-Innovation Foundation (Grant no. 2015RC021).