1. Introduction
Cooperation is a common phenomenon in nature, which is the cornerstone of the stable development of economy and society. However, according to the theory of natural selection [
1], cooperation that benefits others at the cost of one’s own interests is bound to disappear, which inevitably leads to social dilemmas. Therefore, it has become a hot topic to reveal the deep reasons for the spontaneous emergence and maintenance of cooperation in nature and society [
2,
3,
4,
5,
6,
7,
8,
9,
10].
Over the past few decades, scholars in various fields have made great efforts to alleviate social dilemmas from both theoretical and experimental aspects (please see the references [
11,
12,
13,
14,
15] for more information). Among them, evolutionary game theory assumes that players have bounded retionality, providing a general and effective theoretical framework for modeling the conflict between individual interests and interests of collectives [
16,
17]. Furthermore, Nowak summarized these studies into five mechanisms to facilitate the evolution of cooperation, such as kin selection, direct reciprocity, group selection, indirect reciprocity, and spatial or network reciprocity [
18]. In the seminal contribution, Nowak and May found that cooperators in the spatial game model can resist the attacks from defectors by forming tight and compact clusters, thus improving the evolution of cooperation [
19]. Henceforth, spatial or network reciprocity is favored as a novel channel to enrich emergence and maintenance of cooperation, and along this line, a lot of related works have been triggered, for instance, reputation [
20,
21], reward and punishment [
22,
23], incomplete interaction [
24,
25,
26,
27], interactive diversity [
28,
29], and reinforcement learning [
30,
31], to name but a few.
Although these efforts have made great progress, most of them simply suppose that players can only play games with others in a single network, which is obviously inconsistent with actual observation. In fact, players may belong to different subsystems in most cases [
32,
33], although not all of them. In addition, even slight changes in one network may lead to disastrous consequences on another one [
34], which is called cascade failure of interdependence. Therefore, it is instructive to propose a new framework beyond the simple isolated network theory. In recent years, interdependent networks, which are defined as nontrivial combination classes of interconnected networks, have become a primary instrument to quantitatively portray the interactions between a network of network, as well as their constituent elements. In particular, interdependent networks can model various interactions in nature, which makes more and more scholars realize that it is more practical to explore the evolution of cooperation on this framework [
35,
36,
37,
38]. Yet, how to build the interdependence between networks seems to be crucial. In general, players’ fitness directly affects the diffusion and dissemination of their strategies. Thus, rescaling the fitness of a player may be a potential method to achieve the interdependence between a point on one network and the corresponding partner on another network [
39,
40]. As a classical case, Wang et al. investigated the public goods game on interdependent networks and found the evolution of cooperation is positively correlated with bias [
41]. Meanwhile, strategy persistence represents the ability of each player to obtain the highest fitness in the competitive environment. So, it is also a way to establish the interdependent relationship between networks by affecting the imitation of strategies. To this end, Szolnoki and Perc proposed the mechanism of information sharing, that is, players had the intention to slow down the change in strategies when the strategies of the corresponding players in both networks were consistent, which led to the observation of prosocial behavior beyond simple network reciprocity in both populations [
42]. Moreover, heterogeneous coupled relationships also have an important effect on promoting the evolution of cooperation on interdependent networks [
43,
44].
However, the study of the co-evolution of interdependent network structure and strategy imitation can provide more beneficial enlightenment for exploring effective mechanisms to facilitate the evolution of cooperation. From the point of utility, Wang et al. studied the impact of the co-evolution of interdependent networks and strategy on the evolution of cooperation [
45]. In this work, if a player’s utility met or exceeded the given threshold
E, they were allowed to establish connection with the corresponding partner on the other network and obtained an additional reward. Furthermore, based on the perspective of the focal player, Jia et al. introduced the ability to investigate the impact of self-organization on the evolution of cooperation on interdependent networks [
46]. In their model, once the focal player changed their strategy, their ability would be reset to 1; otherwise, their ability would be added by one unit. Then, a player located on one network could build coupled links with the corresponding partner on another network, and then gain additional benefit as reward if their ability is beyond the given threshold
. In particular, we find that although the models of [
45,
46] are similar in design, the resulting dynamics of the evolution of cooperation are quite different. In fact, there are many factors that can prompt two networks to establish effective interdependence in a self-organizing manner. In reality, once a person’s thoughts, behaviors, and strategies are learned or imitated by others, it means that the former has a substantial influence on the latter. Thus, on the one hand, we seek to uncover whether influence-based self-organization can inspire the sense of cooperation on interdependent networks; on the other hand, we are also curious about what kind of dynamics of the evolution of cooperation can be generated on an interdependent network by influence-based self-organization. Herein, we provide a quantitative description of the player’s influence. In concrete terms, if the neighbor
strategy is learned by the focal player
x, the influence of neighbor
y will increase by one unit; otherwise, their influence will decrease by one unit. It should be noted that the co-evolution of influence and strategy in our model is fundamentally different from the co-evolution in [
45,
46], which leads to the fact that the dynamics of the evolution of cooperation are not completely consistent. After that, the player on one network whose influence exceeds the threshold will connect with the corresponding partner on another network. In this way, the interdependence between networks is dynamically built in a self-organized manner. Moreover, due to the fact that the interdependence between networks caused by the co-evolution of influence and strategy is time-varying, it is difficult to accurately predict the effectiveness of the setup to alleviate social dilemmas in advance, which makes this research quite interesting.
In the rest of this thesis, the models and methods are first described in
Section 2. After that,
Section 3 presents the results obtained through a large number of numerical simulations and analyzes them in detail, which is sufficient to show that the simple setup of the model can promote cooperation going beyond isolated network reciprocity. At last, the summative reviews are concluded in
Section 4.
2. Models and Methods
To highlight the influence-based self-organization and eliminate the impact of network heterogeneity on the evolution of cooperation, we assume that players play the prisoner’s dilemma game (PDG, symmetric two-player game) on two interdependent lattices (i.e., network A (bottom) and network B (top)) with periodic boundaries, where the size of each lattice is . In both networks, each point can only denote a player and has four direct neighbors, that is, a Von Neumann neighborhood.
Initially, each player x on the bottom network and on the top network will be set to take a cooperative (C) or defective (D) strategy at random. In the PDG, if two players adopt the same strategies, cooperation results in reward R for each other, while defection makes each obtain punishment P; however, if two players take different strategies, it will enable the cooperator to gain the sucker’s payoff S, while the defector will get lucky and reap the temptation to defect T. In the original PDG, the payoff elements meet the order of and > . The former inequality shows that the point of Nash equilibrium in the PDG is mutual defection (), while the latter inequality indicates that mutual cooperation () can make the system reach Pareto Optimality, which results in the choice conflict between individual interests and interests of collective, i.e., social dilemmas. Within this paper, we just consider the weak version for simplicity, namely , and . Although the ranking of payoff elements in the weak version does not strictly meet the requirements of the original PDG, it can still depict the social dilemmas and the pattern conflicts in the system. After the corresponding players x and sequentially interact with their neighbors in their own network, they can obtain the cumulative payoffs and , respectively.
In order to build the interdependence between networks, we introduce the concept of influence, which is defined as the external expression of a person’s effect on others. The influence of players changes dynamically with the spread of strategies, which is described in more detail later. The utilities of players on both networks are calculated according to the following protocol. Taking the bottom network as an example, if player
influence
exceeds the given threshold
, they will connect with the corresponding partner
on the top network and gain additional benefit, which means that player
utility
is not only determined by the interaction in their own network but is also related to that of the corresponding partner
on another network. However, for player
x who does not build links with the corresponding point
on another network, their utility
only depends on the interaction in their own network. Thus, the utility calculation for player
x can be indicated as
where
indicates the coupled strength of the interdependent networks, while the parameter
determines the influence threshold. If a player’s influence goes beyond
, then they can build a relationship with the the corresponding partner
on another network, and vice versa for player
. Note, unlike the cases of only bidirectional coupled and decoupled relationships in [
44], there are also unidirectional coupled relationships in our model. That is, only players in the bottom network build a coupled relationship with the corresponding partners in the top network, or vice versa.
After player
x finishes the utility
calculation, they will enter the strategy update phase. It is important to point out that strategy imitation can only take place between immediate neighbors in any given network, while prohibition occurs between players on diverse networks, despite the existence of additional connections between them. In other words, player
x on network A cannot imitate the strategy of the corresponding player
on network B at any time. During the procedure of strategy imitation, player
x on network A (and similarly player
on network B) can randomly select one of their nearest neighbor, named
y, to implement strategy imitation with the probability [
47,
48]:
where
K indicates the noises or decision error related to the strategy imitation process. For the sake of simplicity, the value of
K in this paper is set to be 0.1, in accordance with previous works [
49]. The same process applies to
on another network. Moreover, it should be noted that, during the process of strategy update, whether the bottom network A or the top network B is first updated is randomly determined rather than explicitly specified.
In addition, the players
strategy and their influence co-evolve. To be specific, once the nearest neighbor
y is picked out by the focal player
x, they obtain the opportunity to expand their influence. To develop a step further, if the strategy of player
y is successfully imitated, it means their influence will be added by one unit; otherwise, they give up the opportunity to enhance their influence, which will lead their existing influence to be reduced by one unit. Therefore, the dynamic change in player
influence can be expressed as
Here, it should be pointed out that all players have the homogeneous influence value of 1 at the initial stage, which can only be dynamically changed in the interval [0, 100] during the process of evolution; otherwise, it will make no sense to set the threshold . Moreover, if the value of is given at the beginning, it will remain unchanged throughout the process of evolution. However, even so, it still makes the dynamics of the evolution of cooperation difficult to predict, because the co-evolution of player influence and strategy imitation makes it possible for the coupling relationships between networks to change dynamically.
Computer numerical simulations are carried out in random order, which means that each player on both networks has an average opportunity (but not all) to update their strategies within a full Monte Carlo (MC) generation. The creation or termination of additional links between interdependent networks through the propagation of strategies are simultaneously performed, however, independently of each other. Moreover, the results presented later are mainly obtained on a lattice with the size of . To eliminate the impact of finite size, some much smaller () or much larger () lattices are also taken into account. For most data points, the system is able to reach a steady state after MC generations; while for some phase transition data points, the procedure will need to execute longer MC generations (such as or much longer) to guarantee that the system enters these states. Moreover, the proportion of cooperation, , the key indicator to measure the effectiveness of the mechanism in promoting the evolution of cooperation, is the average of the final MC generations. It is known that when cooperators are completely cleared up on the isolated lattice, the threshold of temptation to defect , which is used as a benchmark to test the effectiveness of facilitating the evolution of cooperation on interdependent networks.
3. Results and Analysis
Within this part, we try our best to comprehensively present the relevant simulation results from different angles and deeply analyze the impact of influence-based self-organization on cooperation on interdependent networks in depth, mainly including macroscopic inquiry, time series and pattern analysis, the coupled relationships, and verification of robustness.
3.1. Macroscopic Inquiry
To begin with, under the condition of influence-based self-organization,
Figure 1 shows the
phase diagrams encoding the proportion of cooperation,
, where the values of the coupled strength
in subgraphs from (a) to (c) are 0.2, 0.5, and 0.8, respectively. It should be noted that only the results of the bottom network (i.e., network A) are presented here, and the ones of the top network (i.e., network B) are qualitatively consistent. From
Figure 1, it can be found that no matter what the value of the coupled strength
is, the level of cooperation undeniably decreases with the increase in the temptation to defect
b. However, some interesting phenomena can also be clearly observed. In particular, on the one hand, the level of cooperation is not gradually improved with the enhancement of coupled strength
. It can be found that when the coupling strength
is small (
) or large (
), the system can hardly achieve full cooperation. On the contrary, the state of complete cooperation emerges when the coupled strength
is medium (
). These results imply that the higher coupled strength
is not always better for promoting the evolution of cooperation under the mechanism of influence-based self-organization. On the other hand, the effect of influence threshold
on the evolution of cooperation is nonmonotonic for some values of
b, which is similar to those reported in [
46]. Even if the proportion of cooperation,
, does not significantly change with
when the value of
is quite lower (
), the trend towards nonmonotonic change is more pronounced for larger values of
. In particular, with the increase in
, the proportion of cooperation,
, has gone through a process of first rapid increase, then keeping constant, and quickly decreasing when
and
. That is to say, there is an appropriate influence threshold
making the system achieve the optimal evolution of cooperation.
Interestingly, when the value of
is about within the interval [0, 35], the transformation of the proportion of cooperation,
, with
is sensitive. When
goes beyond the above interval, the proportion of cooperation,
, hardly changes with
. The reason for the these phenomena is that when a player’s influence exceeds the given threshold
, they can build a link with the corresponding partner on another network and gain additional payoff, which enhances their ability to spread the strategy. According to [
44], it is neither that more coupled connections between networks is more conducive to evolutionary cooperation, nor that certainly less is more conducive to evolutionary cooperation. In fact, there is an optimal interdependency ratio leading the best evolution of cooperation. Without doubt, it is easy to build coupled links between networks when the influence threshold
is small. As the influence threshold
increases, the implementation of the coupled relationship becomes more and more difficult, especially for larger values of
b. Thus, it is not difficult to comprehend the trend of the proportion of cooperation,
, changing with
. In particular, apart from the bidirectional coupling and decoupled relationships mentioned in [
44], there are also unidirectional coupled relationships in this research, which might be the reason that the evolution of cooperation does not increase completely with the enhancement of coupled strength
.
It is important to note that the dynamics of the evolution of cooperation in our model are slightly different from those reported in [
45]. In the latter, the proportion of cooperation,
, with the increase in the threshold
E, first rose slowly in the form of a ladder and then fell rapidly for any given values of
b. In addition, within the model of [
45], if the threshold
E is large enough, no matter what the values of
b, a player located on one network cannot build an interdependent relationship with the corresponding partner on another network, thus making two networks become isolated networks. However, from
Figure 1, we can observe the level of cooperation beyond what can be achieved through simple spatial reciprocity on an isolated lattice for quite small values of
b. In particular, these phenomena are even more remarkable when the value of
is medium or large (i.e.,
and
). Thus, in our model, there is always interdependency between networks, regardless of the threshold
.
Anyway, the results presented in
Figure 1 suggest that the mechanism of influence-based self-organization, i.e., the co-evolution of players’ influence and strategy imitation, does promote the evolution of cooperation on interdependent networks.
3.2. Time Series and Pattern Analysis
To reveal the potential reasons for the evolution of cooperation on interdependent networks promoted by influence-based self-organization, we need to conduct a discussion in depth from the microperspective. First of all, we depict the time series of the proportion of cooperation,
, under the temptation to defect
in
Figure 2, where the blue, cyan, and purple curves correspond to the results when the values of the influence threshold
are 0, 8, and 40, respectively.
By examining the time series in
Figure 2, we can further divide these lines into END phase and EXP phase [
50,
51,
52]. It is well known that cooperators cannot survive on simple network reciprocity on an isolated lattice when
. Therefore, due to the fact that defectors have incomparable advantages over the cooperators, this leads to a significant downward trend among all curves at the initial stage, which can be called the enduring (END) stage [
50,
51,
52]. After that, the fate of these curves clearly shows differentiation. To be specific, when
, the proportion of cooperation,
, slowly decreases over time until it completely disappears. When
and
, however, the curves start to rebound after falling to the bottom, which can be called the expanding (EXP) stage [
50,
51,
52]. When
, the blue curve climbs to the highest point, known as the metastable state, and then slightly drops before leveling off, while for
, although the proportion of cooperation,
, will decrease even lower than the case of
, it can eventually control the whole system at last. It should be pointed out that this time dependence of first down and later up is a typical trademark of network reciprocity. In fact, this is a robust phenomenon that applies not only to paired interactions such as the prisoner’s dilemma game [
53] but also to multipoint interactions such as the public goods game [
54]. Furthermore, it could be observed many times on a simple graph [
54,
55] and multilayer structures as well [
45,
46,
50,
51,
52].
Furthermore, the reasons behind these phenomena in
Figure 2 are not difficult to explain. If the influence threshold
is too large, it will lead the influence of most players to hardly reach the goal, so that almost none of them can obtain additional payoffs from peers on another network. Thus, the proportion of cooperation,
, continues to decline until it disappears altogether. If the influence threshold
is too small, it will make it easy for players to achieve their goals, so that most of them can obtain extra benefits from their partners on another network. In neither case can there be significant heterogeneity, which is also not conducive to the evolution of cooperation. Nevertheless, an appropriate influence threshold
can only enable some players to meet the conditions for establishing coupled relationships with partners on another network, and the resulting heterogeneity is a necessary requirement for achieving a higher level of cooperation or even full cooperation.
Next, to explore the self-organization pattern of strategy evolution, we further draw some representative snapshots of the strategies evolution when
in
Figure 3. Initially, cooperators are carefully arranged in the center of the network with the radius of 30 rather than randomly distributed, which will help us to observe the evolution of the strategies more clearly. Each two rows from top to bottom in
Figure 3 are grouped. In each group, the upper line shows the snapshots of the top network and the lower one displays the results of the bottom network. For the groups from top to bottom, the values of the influence threshold
are 0, 8, and 40, respectively. For the columns from left to right in
Figure 3, these typical snapshots of the strategies evolution are obtained at the MC generations of 0, 20, 100, 1000, and 50,000, respectively. Note that players can be further divided into upper class and lower class according to whether their influence exceeds the given threshold
. Thus, cooperators can be distinguished as upper class (green) and lower class (red). Similarly, defectors can also be defined as upper class (yellow) and lower class (blue). Whether cooperators or defectors, only upper class players located on one network can establish additional connections with corresponding partners on another network.
It can be found that no matter what value of the given threshold is, even if the initial setting of the strategies in the bottom and top networks are completely consistent, there are still asymmetries in the process of evolution. Moreover, since the value of influence is initially set to be 1, only upper cooperators and defectors exist in the system at the beginning when . However, it can be found that with the cooperation strategy spreading out from the middle of both networks, lower classes will appear in the process of evolution. In other words, although the value the influence threshold is set quite low, there are still some players who are unable to achieve the spread of strategies, resulting in the loss of influence. Meanwhile, we can also find that full cooperation has not been achieved, even if most players can establish coupled relationships with their partners on another network, suggesting that nearly homogeneous coupled connections between networks are insufficient to support a higher level of cooperation. Nevertheless, when the value of the influence threshold is set to be 8, only lower class players appear in the system at the initial stage, i.e., the circular red lower class cooperators are located at the center of the quadrate of blue lower class defectors. Over time, green spots and yellow flecks begin to appear, indicating that some players’ influence has gone beyond the given threshold , and they have built connections with their partners on another network. It is particularly exciting that full cooperation has finally been achieved. In the last two lines, even the value of the influence threshold is too high (). Although a few players achieve the dream of transcending class, cooperation is not maintained. To be more precise, the system ends up with only defectors, indicating that too high an influence threshold is more unfavorable to the evolution of cooperation. After all, once the influence threshold is set too high, most players can hardly achieve class transcendence, and ultimately, can not establish coupled relationships with their partners on another network. Note, if there is only one strategy (cooperation or defection) in the system, then the propagation of the strategy will not be involved, and the player’s influence will be frozen. At this point, the case of is similar to the one of .
The time series of the proportion of cooperation, , and the spatial pattern of the strategy evolution are not only in line but further reveal the reason that influence-based self-organization promotes cooperation on interdependent networks.
3.3. The Coupled Relationships
It is necessary to study coupled relationships between networks, which can help us to further clarify the details of how influence-based self-organization inspires the sense of cooperation on interdependent networks. To this end, the time series of different types of relationships between networks are plotted in
Figure 4, where from subgraphs (a) to (c) the values of influence threshold
are 0, 8, and 40, respectively. It should be pointed out that, in
Figure 4, the coupled relationships also include the special case of decoupled links between corresponding points on interdependent networks. From
Figure 4, it can be clearly found that when the influence threshold
is low (
) or moderate (
), the relationship of DD type has a temporary rise at the initial stage, while the others have different degrees of decline. In addition, with time evolution, the relationships of DD type or CC type almost reach the peak or bottom at the same time. After then, the relationship of DD type begins to decline and that of CC type starts to rise. The difference between subgraphs (a) and (b) is that, on the one hand, the crest of the relationship of DD type or the trough of the relationship of CC type in subgraph (b) are higher or deeper than those in subgraph (a) in the initial stage; on the other hand, the relationships of all types (i.e., CC, CD, DC, and DD) in subgraph (a) are balanced at last, while only the relationship of CC type in subgraph (b) exists and occupies the whole system in the end. However, in subgraph (c), the value of influence threshold
is set too high (
), and the coupled relationship of DD type continues to increase, while the others inevitably decline, and finally, only the coupled relationship of DD type remains in the system. By comparing these results, it can be seen that the relationship of CC type plays a incomparable role in improving the level of the evolution of cooperation during the process of evolution.
Since the relationships between networks includes bidirectional, unidirectional, and decoupled types, we need to further refine the relationships in
Figure 4, which will help us to more clearly uncover the role of influence-based self-organization in promoting the evolution of cooperation on interdependent networks. For this point, we plot the time series of each specific relationship in
Figure 5, where from subgraphs (a) to (c) the values of influence threshold
are 0, 8, and 40, respectively. When
, although the curves representing coupled or decoupled relationships of different types show varying degrees of alternating oscillations over time, only the relationships of bidirectional type (i.e., bi-CC, bi-CD, bi-DC, and bi-DD) are preserved in the end. When
, the curves of different types also exhibit similar alternating oscillations over time evolution. However, it is interesting to note that, even if the system finally enters the full cooperation state, only the coupled bidirectional relationships of CC type (i.e., bi-CC) and the decoupled relationships of CC type (i.e., de-CC) remain at last, and the bi-CC is much higher than de-CC. Different from the results in subgraphs (a) and (b), the curves representing relationships of various types do not exhibit alternating oscillations when
. Specifically, except that the decoupled relationship of DD type (i.e., de-DD) continues to climb and the coupled unidirectional relationships of DD type (i.e., un-DD) slightly rises, not only do the other curves of decoupled types continue to decline until they completely disappear, but also the curves of other types never emerge. These phenomena indicate that the bidirectional coupled relationships, especially the bidirectional coupled relationship of CC type, play an incomparable role in motivating the sense of cooperation under the mechanism of influence-based self-organization.
Through the above analysis, we learned that, on the one hand, for the prepared initial settings (in the middle two columns of
Figure 3), the system can achieve full cooperation, even if the cooperators are far lower than defectors; on the other hand, for random initial settings, the relationship of CC type ( in
Figure 4b), especially the bidirectional coupled relationship of CC type (i.e., bi-CC in
Figure 5b), plays a key role in promoting the evolution of cooperation. In order to further confirm the above assertion, we investigate the characteristics snapshots under the conditions of the special strategies distribution settings, as shown in
Figure 6. According to Hauert [
56], in the isolated lattice of a Von Neumann neighborhood, the
cooperative cluster is the minimum setting to ensure that cooperation can flourish under the PDG. For this purpose, in subfigure (a) of
Figure 6, we randomly set three
minimum cooperative clusters with different relationships at the initial stage, such as unidirectional coupled relationships (position A), bidirectional coupled relationships (position B), and decoupled relationships (position C). In addition, the generations of MC from subgraphs (a) to (f) are 0, 20, 300, 1000, 5000, and 50,000, respectively.
It can be found that, over time, only the cooperative cluster with the decoupled relationships (position C) disappears at last, but the ones with bidirectional (position B) and unidirectional (position A) coupled relationships can survive and thrive, which is markedly different from the phenomenon reported in [
46]. Moreover, bidirectional coupled relationships of CC type can emerge in unidirectional coupled relationships of CC type clusters, which further strengthens the competitive advantage of the cooperative strategy. From
Figure 6, especially in subfigures (c) and (d), it can be found that defectors with decoupled relationships (light yellow) first turn into decoupled relationships (light cyan or blue) with only one cooperator, then into weak (unidirectional) coupled relationships (medium cyan or blue) with only one cooperator, then into cooperators with weak (unidirectional) coupled relationships (medium red), and finally into cooperators with strong (bidirectional) coupled relationships (dark red). This is because once players establish coupled relationships with their corresponding partners on another network, they can obtain additional payoffs, thus enhancing the spread of the cooperative strategy. So, defectors at the boundaries of cooperative clusters (positions A and B) can be transformed into cooperators. Thus, clusters in these snapshots have very typical characteristics, namely, bidirectional coupled clusters of CC type are primarily in the center, while the ones of DD type are mainly at the boundaries. Over time, the cooperation gradually diffuses from the positions A and B until it occupies the whole system, as shown in subfigures (e) and (f).
To date, it seems that we can come to the conclusion that, under the mechanism of influence-based self-organization, the bidirectional or unidirectional coupled relationship of CC type plays a pivotal role in alleviating social dilemmas and improving the evolution of cooperation on interdependent networks.
3.4. Verification of Robustness
All the above results assume that, if a player’s influence is beyond the given threshold
, they can build the coupled relationship with the corresponding partner on another network. A question that naturally arises is how cooperation dynamics evolve on interdependent networks, if the players who meet the conditions are randomly rather than correspondingly coupled to partners on another network. To this end, we relax the assumption that players in the two-layer network establish the coupled relationships. As an example,
Figure 7 presents the results of these two cases at the coupled strength
, where subgraph (a) is the result of the corresponding coupled relationship, and subgraph (b) is the case of a random coupled relationship.
It is not difficult to find that, in both cases, under the given thresholds, , except for the proportion of cooperation, , nonlinearly decreases with the increase in b (the temptation to defect), and there are moderate values of influence threshold to achieve the optimal evolution of cooperation in both cases. If the condition for players to build a coupled relationship with partners in another network is relaxed, however, the differences in the trend of the evolution of cooperation are still obvious. On the one hand, for some given thresholds , the range of parameters b that can maintain full cooperation is significantly reduced; on the other hand, for all given thresholds , the proportion of cooperation, , and the threshold for the disappearance of cooperators are dramatically lower. Nevertheless, under the mechanism of influence-based self-organization, the level of cooperation on interdependent networks is still better than that on an isolated lattice.
4. Conclusions and Discussion
To summarize, we creatively introduce influence-based self-organization into interdependent networks and investigate the effect of this simple setup on the evolution of cooperation in a spatial PDG. Within this work, influence represents the ability of players to act on others. To be specific, if the player’s strategy is successfully learned by others in the process of strategy imitation, their influence will increase by one unit; otherwise, their influence will be damaged and lose one unit. Then, once the player’s influence exceeds the given threshold , they can build a coupled relationship with the corresponding partners on another network. The co-evolution of influence and strategy imitation may lead to unstable connection between networks, which makes it difficult and unattainable to predict the dynamics of the evolution of cooperation on interdependent networks. Thus, our research is quite interesting and meaningful.
Through sufficient simulations, it can be found that there is a moderate value of influence threshold
for cooperation to evolve the best, which is similar to the results reported in [
46]. According to [
44], too many or too few coupled links between networks are not conducive to heighten the evolution of cooperation. In this model, if the value of
is set too small, it is easy for most players to build coupled relationships with the corresponding partners on another network, while if the setting of
is too large, a lot of players can hardly build connections with the corresponding partners on another network. However, we still found a phenomenon that seems to go against cognition, that is, the level of cooperation is not positively related to the coupled strength
, which is different from the results in [
44]. In fact, this phenomenon is not difficult to explain. In [
44], there are are only bidirectional and decoupled relationships in the system, while unidirectional coupled relationships can also appear in our model. The diversification of coupled relationships may be the root cause of the above observation.
Then, through rich microscopic statistical analysis, we are surprised to find that the coupled relationships of CC type, especially bidirectional or unidirectional links of CC type, play an irreplaceable role in alleviating social dilemmas. This is because the existence of bidirectional or unidirectional coupled relationships of CC type enables cooperators to obtain additional benefits from their partners on another network, which significantly enhances their competitiveness in the process of strategy imitation. Thus, it directly causes defectors around the cooperative clusters to be assimilated. Furthermore, if players randomly rather than correspondingly construct coupled relationships with partners on another network, even if the evolution of cooperation is compromised, it still exceeds the level of cooperation that the isolated network can guarantee.
The current results are convincing enough that the mechanism of influence-based self-organization can indeed improve the evolution of cooperation on interdependent networks, which enriches and develops the theoretical research on the co-evolution of network structure to alleviate social dilemmas. In view of the fact that coupling among different subsystems is widespread in nature and society, and the dynamics of the evolution of cooperation obtained in our model are more or less strikingly different from those reported in [
45,
46], these models are still similar in design. In reality, there are many more factors affecting the interdependence of different subsystems. Thus, we hope that our study will be valuable and can arouse the interest of researchers in related areas.