International Journal of Communications, Network and System Sciences
Vol.09 No.12(2016), Article ID:72966,21 pages
10.4236/ijcns.2016.912045
Deterministic Multicultural Dynamic Networks: Seeking a Balance between Attractive and Repulsive Forces
Kristina B. Hilton, G. S. Ladde
Department of Mathematics and Statistics, University of South Florida, Tampa, FL, USA

Copyright © 2016 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/



Received: October 20, 2016; Accepted: December 20, 2016; Published: December 23, 2016
ABSTRACT
An important issue in society is the attempt to balance communities working co- operatively and cohesively with one another while allowing members the ability to retain individuality and fostering an environment of cultural diversity. We seek to study the cohesive properties of a culturally diverse dynamic social network. By considering a multi-agent dynamic network, we seek to model a social structure and find conditions under which cohesion and coexistence are maintained. We present a specific illustration that serves to establish the framework in which explicit sufficient conditions in terms of system parameters are found for which the network is cohesive. By utilizing Lyapunov’s Second Method and comparison equations, we are able to find such conditions for the given illustration. Further, for the illustration, we decompose the cultural state domain into invariant sets and consider the behavior of members within each set. Moreover, we analyze the relative cultural affinity between individual members relative to the center of the social network. We also demonstrate how conservative the estimates are using Euler type numerical approximation schemes based on the given illustration. We are then able to consider how changes in the various parameters affect the dynamics of the illustrated network. By gaining such insight into the behavior of the illustrated network, we are able to better understand the impact of both attractive and repulsive influences on the network. This leads to establishing a schema for helping when creating policies and practices catered to promoting both diversity and cohesion within a cultural network.
Keywords:
Multi-Agent Network, Cohesiveness, Lyapunov Second Method, Invariant Sets, Social Network

1. Introduction
The goal of this work is to explore the cohesive properties and behavior of individual members of a dynamic social network. Dynamic network systems are often used to model the behavior patterns of animals, autonomous vehicles, the spread of a contagion, traffic flow, and many other types of situations. In this work, we are interested in modeling the behavior of members within a social network. In particular, we are interested in the cultural shifts of members within culturally diverse groups. We seek to better understand the internal and environmental factors that may foster a sense of cooperation between members of the network while allowing individuality and diversity to be maintained and enhanced.
One of the concepts studied using network dynamics is that of consensus [1] [2] [3] [4] . In such models, the conditions under which a group collectively comes to an agreement are studied. Another question of interest for such a network is when the group might subdivide into smaller subgroups each converging to a consensus but never reaching a consensus as an overall group. Such dynamic network models are useful in many areas beyond just social networks. For example, work in both biological networks and control theory considers such large-scale dynamic models in the context of connectivity, stability and convergence [5] [6] [7] . Using these ideas, much of the work done in these areas look to develop consensus seeking algorithms and consider long term co-existence, cohesion, and stability of the network under consideration [8] [9] [10] [11] .
Cohesion within a social network is a current topic of great interest and many authors have done research within this area [12] [13] . The concepts of cohesion and cooperation within a group are often multi-faceted, dynamic, and complex but are important concepts when trying to better understand how nations or human groups interact and function [14] . As Knoke and Yang note [15] , it is social cohesion that enables information to spread and allows a group to act as a unit rather than individuals.
We often seek to create situations for which people of different backgrounds and beliefs are able to coexists and create a thriving sense of community. We seek to better understand the group dynamics of such a society in order to create policies and practices that encourage a sense of community among individuals from a variety of cultural backgrounds. We use the term multicultural social network to describe a social network in which the agents have a diverse cultural and/or ethical background and are actively seeking to enhance and to maintain diversity with harmony and prosperity. In such a network, the goal of agents is not in approaching a consensus but rather the ability to live and work cooperatively with one another for a common good and goal. For example, consider a population in an area for which there exists a sub-populace of immigrants. In such a situation, the subgroups or sub-communities of immigrants desire to be an integral part of the community and seek to be respectable productive members of the community and the society in general while retaining their cultural diversity. In exploring the dynamics of a multi-cultural network, we are looking to better understand the delicate balance between a culturally diverse cohesive social structure and a social structure for which cohesion does not exists. For when cohesion is lacking in the social network, cooperation may not be as prevalent and we begin to see features such as segregation, violence, economic destabilization and crime within the network.
We seek to model such a situation and better understand the social dynamics of a group seeking to find such a balance. In particular, we are looking to model a dynamic social network for which there is a balance between consensus and cohesion. We wish to model a network that is cohesive but for which there is not a consensus of culture, that is to say the network does not develop a singular cultural identity. In doing so, we are interested in better understanding the cohesive properties of a multicultural social network. We present a prototype of a dynamic model for which we explore the features of such a network. The presented example is used to exhibit the quantitative and qualitative properties of the network. Further, the techniques used are computationally attractive, easy to verify and algebraically simple. In addition, the presented results are in terms of network parameters that characterize the attributes of the network. The byproduct of this provides tools for planning and making policies regarding a dynamic network.
In Section 2, we present an example of such a network, as well as assumptions and notations used throughout this work. In Section 3, using Lyapunov’s Second Method and the comparison method [16] , we consider the dynamics of individual members within the network. In Section 4, long and short term behaviors of group members and invariant cultural state sets are investigated. In Section 5, we consider numerical simulations of the network to better understand the extent that conservative estimates in Sections 4 and 5 are in a given example. Finally, in Section 5, we consider parametric variations within the model affecting the dynamics of the network. Further, we will consider how the model relates to a multicultural network.
2. Problem Formulation
We wish to model a multi-cultural social network and therefore desire to capture the behavior of individual agents who are seeking to belong to the group but also wanting to retain individuality and diversity from other agents. In order to do so, we therefore consider dynamic equations subjected to both attractive and repulsive forces. In [17] , one such function considered when modeling biological dynamic networks is given by
(1)
where a, b and c are positive constants and
. The function g has long-range attraction and short-range repulsion.
In the following, we formulate a modified version of a network dynamic model in which individuals seek to retain a balance between individual member identity and a group/community membership. We consider a network whose dynamics are described by incorporating a long-range attraction and short-range repulsion similar to that in Equation (1).
Let us consider a network of m members. For each member
,
is a cultural position at time t > 0. The vector
can be representative of many various aspects of culture such as beliefs, behaviors, ways of life, etc. depending on the network being considered. Further, let us define a relative cultural state of ith member with a kth member of the community as
, and a center of cultural state of the network
(2)
Consider the network whose dynamic is given by
(3)
The constant coefficient parameters, a, b, c and
represent the weight of the social moderation attractiveness (q), the repulsive forces (a), the rate of decay of the long range attractiveness (c), and the long-range attractiveness (b) between individual members and social groups. We say that the network is cohesive if there exist constants T and M such that
and
that is to say that the members of the network after some point in time remain within a certain distance of the network center. Further, we say that the network reaches a consensus if
as
for all
. In this case each member of the network draws closer to each other and the network center. Moreover, we define the term relative cultural affinity to be
, the distance between the cultural vector states of members
and
.
3. Characteristics of the Network
In this section, we wish to explore the dynamics of the agents with the network dynamic described by Equation (3). We will be considering the cohesion, qualitative and quantitative properties such as the overall stability of the network center, and various types of invariant sets. While exploring these ideas, we will also consider what happens as the size of the network increases and what roles the parameters a, b and c play within the model. Moreover, the presented example is utilized to exhibit the quantitative and qualitative properties of the network. In order to accomplish such a task, we utilize Lyapunov’s Second Method [16] . This method is algebraically simple, easy to verify and computationally attractive. Furthermore, the results depend on the system parameters a, b, c and q.
Let us first consider the dynamic of the network center,
as defined in Equation (2). We note that

and 


Therefore, the dynamics of the transformed network are given by

Dynamic Equation (6) can be useful in modeling a variety of multicultural social networks. In Equation (6), the magnitude of the repulsive force is represented by 


magnitude of the deviation of the ith agent’s cultural state from the center of the network. Attractive influences can be thought of as attributes that bring people to active membership within the group. Social acceptance, gaining social status, economic opportunity, career growth, common purpose and membership, personal development, and a sense of mutual respect, trust and understanding are examples of attractive influences within a social cultural network. Repulsive forces are attributes that create some desire for individuals to leave or be less involved in the group or to preserve some personal identity from one other with their individual magnitude of inner repulsive force. A desire to retain a sense of individuality, economic or emotional cost, interpersonal conflict within the group, or disagreement with parts of the overall philosophies of the group are forces that may be considered as repulsive forces. In short, economic, educational, and social inequalities coupled with the race, gender, ethical and religious bias are sources of repulsive forces. A balance between the total attraction and repulsive forces attributes to a general sense of individual agents maintaining a “live and let live” philosophy for the greater benefit of the community and the common good of society.
In order to better understand the dynamics of Equation (6), we will use Lyapunov’s Second Method in conjunction with the comparison method [16] . These methods will provide a computationally attractive means to better understand the movement of members within the network. To that end, let us consider a candidate for energy function defined by

Then the differential of V along the vector field generated by Equation (7) is given by

where

In Subsections 3.1 and 3.2, we will find upper and lower estimates for 
3.1. Upper Estimate of
In this subsection, we seek constraints on a, b, c, and q such that for 


where 

By considering the derivative of the function,

We note that

has a global maximum when 

From Equation (12), Equation (13), and the fact that


Assumption H_1: Suppose

and let




Using Equation (16) along with the comparison theorems [16] , we establish the following result.
Lemma 1. Let V be the energy function defined in Equation (7), 


where 

provided that

Proof: Under the assumptions of the lemma and using the standard argument [16] combined with the above discussion, the proof of the lemma follows from Equation (16). The cohesiveness of the network follows by definition as the solution to Equation (18) is bounded. □
Remark 1: We remark that the assumption H1 is an alternative sufficient condition as: From Equation (14), we have

where

for any
3.2. Lower Estimate of
Next, we look to establish a lower estimate of 

where 

Imitating the argument used to arrive at Equation (14) and noting that, for 



Assumption H2: Let us define

and



Using Equation (25) along with the comparison theorems [16] , we establish the following result.
Lemma 2. Let V be the energy function defined in Equation (7) and zi a solution of the initial value problem defined in Equation (6). For each 

where 

Provided

Proof: Under the assumptions of the lemma and using the standard argument [16] combined with the above discussion, the proof of the lemma follows from Equation (25). □
Remark 2: A remark similar to remark 1 is as follows: From Equation (23), we have

where

For any
4. Long and Short Term Behavior of Members and Invariant Sets
After First let us note from

Similarly, from the solution of the comparison differential Equation (17) and Lemma 1, we note that

Therefore, by Lemmas 1 and 2, when

From Equations (31), (32) and (33), we consider one case and the associated invariant sets. First, let us consider the case for which

Let us further suppose that it is the case that

For 


In the following, we state and prove a few qualitative properties of the solution process of the center of the multi-agent determinist dynamic network described by Equation (3). The following result exhibits the major influence of long range attractive forces.
Figure 1. An example in 
Theorem 1. For



Further, if




Proof. Let 



Considering the non-homogeneous comparison equation,

it follows that

where 



Then, by using the method of variation of parameters, the solution to the non-ho- mogenous differential equation given in Equation (39) is given by

where

Using separation of variables, the solution of the homogeneous differential equation is given by

and

Therefore, from Equation (42),

where

Let 

We note that 


As the limit as 




and so 










Suppose that it is the case that 






Thus, 



Theorem 2. Let the hypotheses of Lemmas 1 and 2 be satisfied. Then
1) the set 

2) the set 


3) set C is either self-invariant or 
4) the set 
5) the set 

Proof. For


for












for















Let us expand upon the results of Theorems 1 and 2. First, let us note that these two theorems provide the qualitative and quantitative requirements on the cultural state parameters to insure that the model is cohesive (Theorem 2) and simultaneously does not reach a cultural consensus (Theorem 1). We introduce the definition of cultural bound to describe the boundary between two cultural sets, dividing the degree of individual versus community level interaction domains of the cultural state. Suppose














Suppose the initial value, that is the function of the magnitude of the cultural state, 























Turning to the upper comparison equation, we can consider the behavior of the transformed network members whose initial positions are in the sets 




















Further, from Lemmas 1 and 2, if all parameters other than the size of the network are held constant, then as the size of the network increases, so also the difference between the upper and lower bounds on the relative cultural affinity between agents and the center of the network increases. Naturally, increasing the size of the network leads to the concept of the crowding effect. Competition over ideology or cultural traits creates a stronger desire for agents to retain more of their individuality within the society or group. Cultural subgroups that have a high degree of separation in terms of their relative cultural affinities are an emergent characteristic of such large scale multi- cultural networks. In the modeling for members whose cultural state is in 
5. Numerical Simulation
In this section, using Euler’s type numerical to approximation scheme applied to Equation (6), we consider the numerical simulations for the network dynamics governed by Equation (6). The goal is to compare the long term behavior of the simulated solution with the theoretical long term behavior given in Section 4. We consider a network consisting of 50 members with parameters a = 0.5, q = 0.04, b = 0.41, and c = 2. Further, we note that in this case, 

In this example, the conditions for the invariant sets given in Theorems 1 and 2 in Section (4) are satisfied. Hence, for 





Next, we consider the network with the same initial values with the parameters a = 0.25, b = 0.14 and q = 0.04. In this case 

For 




Figure 2. Euler approximation of the solution to the differential equation given by Equation (6) with parameters a = 0.5, b = 0.41, and q = 0.04 yielding the cultural positions over time t for the full 50 members of the network.
Figure 3. Euler approximation of the solution to the differential equation given by Equation (6) with parameters a = 0.5, b = 0.41, and q = 0.04 yielding the cultural positions over time t for six of the network members.
Figure 4. Euler approximation of the solution to the differential equation given by Equation (6) with parameters a = 0.25, b = 0.14, and q = 0.04 yielding the cultural positions over time t for the full 50 network members.
Figure 5. Euler approximation of the solution to Equation (6) with parameters a = 0.25, b = 0.14, and q = 0.04 yielding the cultural positions over time t for six of the network members.
The last case we considered is the network with the same initial positions with the parameters a = 0.5, b = 0.18 and q = 0.2. Thus, with the given parameters, 

For 




6. Conclusions
We have considered requirements on network parameters for long term qualitative properties of the network. We develop a model and establish conditions on the parameters that ensure a balance between cohesion and consensus. Further, we have considered how the initial cultural state of a network member affects the behavior of that member over time. The presented conditions of the system are algebraically simple, easily verifiable and computationally attractive. The developed results provide a tool for planning, decision making, and performance. Furthermore, the presented sufficient conditions are conservative but robust, verifiable, and reliable. From the above conditions, we are able to consider certain dynamic properties of the social networks governed by Equation (3).
Figure 6. Euler approximation of the solution to the differential equation given by Equation (6) with parameters a = 0.5, b = 0.18, and q = 0.2 yielding the cultural positions over time t for full 50 network members.
Figure 7. Euler approximation of the solution to the differential equation given by Equation (6) with parameters a = 0.5, b = 0.18, and q = 0.2 yielding the cultural positions over time t for six network members.
In this paper, we explored the features of a multi-cultural network with dynamics described by a specific differential equation and the long term stability and behaviors of individual members within such a network. We are interested in further exploring social networks in the context of better understanding the relative cultural affinity between agents 


Acknowledgements
The authors would like to acknowledge the research support by the Mathematical Sciences Division, US Army Research Office, Grants No. W911NF-12-1-0090 and W911NF-15-1- 0182.
Cite this paper
Hilton, K.B. and Ladde, G.S. (2016) Deterministic Multicultural Dynamic Networks: Seeking a Balance between Attractive and Repulsive Forces. Int. J. Communications, Network and System Sciences, 9, 582-602. http://dx.doi.org/10.4236/ijcns.2016.912045
References
- 1. Degroot, M.H. (1974) Reaching a Consensus. Journal of the American Statistical Association, 69, 291-293.
https://doi.org/10.1080/01621459.1974.10480137 - 2. Friedkin, N.E. (2014) Complex Objects in the Polytopes of the Linear State-Space. arXiv:1401.5339.
- 3. Acemoglu, D., Como, G., Fagnani, F. and Ozdaglar, A. (2013) Opinion Fluctuations and Disagreement in Social Networks. Mathematics of Operations Research, 38, 1-27.
https://doi.org/10.1287/moor.1120.0570 - 4. Ma, H. (2013) Literature Survey of Stability of Dynamical Multi-Agent Systems with Applications in Rural-Urban Migration. American Journal of Engineering and Technology Research, 13, 131-140.
- 5. Ladde, Gangaram S. and Siljak, Dragoslav D. (1975) Connective Stability of Large-scale Stochastic Systems. International Journal of Systems Science, 6, 713-721.
https://doi.org/10.1080/00207727508941855 - 6. Ladde, G.S. and Lawrence, B.A. (1995) Stability and Convergence of Large-Scale Stochastic Approximation Procedures. International Journal of Systems Science, 26, 595-618.
https://doi.org/10.1080/00207729508929055 - 7. Anabtawi, M., Sathananthan, S. and Ladde, G.S. (2000) Convergence and Stability Analysis of Large-Scale Parabolic Systems under Markovian Structural Perturbations-I. International Journal of Applied Mathematics, 2, 57-85.
- 8. Cao, Y., Yu, W., Ren, W. and Chen, G. (2013) An Overview of Recent Progress in the Study of Distributed Mult-Agent Coordination. IEEE Transactions on Industrial Informatics, 9, 427-438.
https://doi.org/10.1109/TII.2012.2219061 - 9. Zhu, Y.-K., Guan, X.-P. and Luo, X.-Y. (2013) Finite-Time Consensus for Multi-Agent Systems via Nonlinear Control Protocols. International Journal of Automation and Computing, 10, 455-462.
https://doi.org/10.1007/s11633-013-0742-9 - 10. Hu, H.-X., Yu, L., Zhang, W.-A. and Song, H. (2013) Group Consensus in Multi-Agent Systems with Hybrid Protocol. Journal of the Franklin Institute, 350, 575-597.
https://doi.org/10.1016/j.jfranklin.2012.12.020 - 11. Huang, M. and Manton, J.H. (2009) Coordination and Consensus of Networked Agents with Noisy Measurements: Stochastic Algorithms and Asymptotic Behavior. SIAM Journal on Control and Optimization, 48, 134-161.
https://doi.org/10.1137/06067359X - 12. Friedkin, N.E. (2004) Social Cohesion. Annual Review of Sociology, 30, 409-425.
https://doi.org/10.1146/annurev.soc.30.012703.110625 - 13. Bruhn, J. (2009) The Group Effect, Social Cohesion and Health Outcomes. Springer, New York.
- 14. Axelrod, R.M. (1997) The Complexity of Cooperation: Agent-Based Models of Competition and Collaboration. Princeton University Press, Princeton.
- 15. Knoke, D. and Yang, S. (2008) Social Network Analysis. Sage, Thousand Oaks.
- 16. Lakshmikantham, V. and Leela, S. (1969) Differential and Integral Inequalities: Ordinary Differential Equations. Academic Press, Pittsburgh.
- 17. Gazi, V. and Passino, K.M. (2003) Stability Analysis of Swarms. IEEE Transactions on Automatic Control, 48, 692-697.
https://doi.org/10.1109/TAC.2003.809765 - 18. Ladde, G.S. and Leela, S. (1972) Analysis of Invariant Sets. Annali di Matematica Pura ed Applicata, 94, 283-289.
https://doi.org/10.1007/BF02413615










