Deterministic Multicultural Dynamic Networks : Seeking a Balance between Attractive and Repulsive Forces

An important issue in society is the attempt to balance communities working cooperatively 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.


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.

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 ( ) where a, b and c are positive constants and n y ∈  .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 x ∈  is a cultural position at time t > 0. The vector ( ) i x t 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 ik i k x x x = − , and a center of cultural state of the network Consider the network whose dynamic is given by ( ) .
The constant coefficient parameters, a, b, c and q + ∈  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 0 t T t ≤ ≤ 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 1, i I m ∈ .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 ( ) ( ) ( ) x t x t t x x t t x = − , the distance between the cultural vector states of members i x and k x .

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, x as defined in Equation (2).We note that 0 and x is a stationary center of the network.We define the transformation of the net- 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 i am z and the magnitude of the long range attractive force is described by 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 spectively.Using Lyapunov's Second Method and the comparison method [16], we will then use these estimates to consider the behavior of the agents over time t.For instance, we consider the stability of the network, and establish invariant sets for the network.

i LV z
In this subsection, we seek constraints on a, b, c, and q such that for i z outside of a given ball, we can establish an upper estimate of ( ) i LV z .We will then use these as- sumptions in conjunction with the Lyapunov method and comparison theorem [16] to establish the case for which where ( ) ( ) is the maximal solution of a comparison differential equation through ( ) By considering the derivative of the function, ( ) We note that has a global maximum when with a maximum value of From Equation ( 12), Equation ( 13), and the fact that sin 1 i z ≤ , for ( ) reduces to: Assumption H_1: Suppose 1 0 qm − > .Let us define ( ) and let ( ) { } . Further, let us denote the compliment of the ( ) , Equation ( 14) yields the following inequality: 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), i z be a solution of the initial value problem defined in Equation (6).For each ( ) satisfying the differential inequality Equation ( 16), it follows that the network is cohesive and where ( ) r t is the maximal solution of the initial value problem ( ) 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 H 1 is an alternative sufficient condition as: From Equation ( 14), we have where ( ) { } ( ) for any 1 0 r ≥ .

i LV z
Next, we look to establish a lower estimate of where ( ) is the minimal solution to a comparison equation through ( ) Imitating the argument used to arrive at Equation ( 14) and noting that, for 0, x α α .
Lemma 2. Let V be the energy function defined in Equation (7) and z i a solution of the initial value problem defined in Equation (6).For each ( ) where ( ) is the minimal solution of the initial value problem ( ) 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 ( 0, : , where ( ) For any 2 0 r > .

Long and Short Term Behavior of Members and Invariant Sets
After First let us note from ( ) t ρ , the minimal solution to the initial value problem in Equation (23) in Lemma 2, we find ( ) ( ) 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 2 1 β β ≤ .That is, let us suppose that ( ) Let us further suppose that it is the case that ( ) For 1 β and 2 β , let us define the following sets (Figure 1): 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.
 , a neighborhood of the center x , then Equation (9) reduces to the inequality Further, if ( ) 0 2 a q u > , there exists 0 Considering the non-homogeneous comparison equation, ( ) ( ) where ( ) u t is the minimal solution of Equation (39) when ( ) û t be the solution of the homogeneous differential equation ( ) Then, by using the method of variation of parameters, the solution to the non-homogenous differential equation given in Equation ( 39) is given by where ( ) ( ) Using separation of variables, the solution of the homogeneous differential equation is given by ( ) ( ) Therefore, from Equation ( 42), where ( ) Let ( ) g t be the function defined as We note that ( ) As the limit as t → ∞ of ( ) As ( ) ≤ , it must be the case that 0 M > and hence Equation (42) reduces to ( ) ( ) Suppose that it is the case that ( ) 0 2 a q u > and so the solution ( ) it follows that Equations (51) has the lower bound ( ) Let the hypotheses of Lemmas 1 and 2 be satisfied.Then 1) the set ( ) ( ) the hypotheses of Lemmas 1 and 2 are satisfied.Thus by the application of these Lemmas, we ) and ( ) 0 0 , , r t t r are the minimal and maximal solutions of the comparison differential equations given in Equations ( 27) and ( 18) respectively.Moreover, for ( ) , for 0 t t ≥ .From Equation (54), 0 i z E ∈ and the definitions of self-invariant and con- ditionally invariant [18], it follows that statement 1) is valid.The proofs of 2), 3) and 4) follow by imitating the argument used in the proof of 1).For proving statement 2).For decreasing and the proof of 3) and 4) follows from 1) and 2).Similarly, the proof for statement 5) also follows by imitating the argument used in 1).For 0 i z B ∈ , the solu- tions to the comparison equation given by Equation ( 27), ( ) a qm as t → ∞ .Therefore, by Equation ( 54), ( ) □ 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 in- dividual versus community level interaction domains of the cultural state.Suppose Therefore, if the cultural state of the ith member i x of the network is such that the relative cultural affinity between i x and the center, x of the network is sufficiently close to zero, then the agent's cultural state is repulsed from the center.That is to say, the membership of the social network will obtain and then maintain a relative cultural affinity between members and the center that is bounded below by a value strictly greater than zero.Once the state of the ith member i z has moved away from the cen- ter, it may be the case that i z remains in A or the case that the state i z moves to the cultural set B , at which time the agent's cultural state behavior will follow that of another category of membership described by the cultural state set B discussed below.
Suppose the initial value, that is the function of the magnitude of the cultural state, , then the solution decays and asymptotically approaches to the threshold limit ( ) 2 a qm from above.Therefore, if i z is a member of the transformed social network such that i z B ∈ , then by Theorem 2, over time, i z moves to the cultural bound of the set C .It may also cross the cultural bound or it may be the case that i z approaches asymptotically to the cultural bound of C .Similarly, if i z C ∈ , i z may stay in C , approaching the cultural bounds of sets B and/or D or it may be the case that i z crosses the cultural bound of D at which point the member will behave as other members of D .However, if i z C ∈ , even though it may approach the cultural bound of B , it will never cross the bound.In terms of a given social network, this implies that members with a distinct enough cultural state from the weighted average of cultural states will retain that distinctiveness of culture.Thus, if the relative cultural affinity between a member i x and the center of the network is at least ( ) a qm initially, then the relative cultural affinity will always be at least that value., , r t t r grows and approaches asymptotically to the value ( ) , the solution decays and approaches asymptotically to the limit from above.Therefore, if i z D ∈ , i z may approach and cross the cultural boundary of C (but will remain in C D ∪ ) or i z may approach but not cross the cultural boundary of E .For i z E ∈ , i z may either cross the cultural boundary of D or the member's cultural state will approach asymptotically to the cultural boundary of D .Thus, for agents i x within the network whose initial relative cultural affinity with respect to the center is sufficiently large, as t → ∞ , the relative cultural affinity will remain large and the although the agent is attracted back towards the center of the network, the relative affinity is bounded below by ( ) 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 multicultural networks.In the modeling for members whose cultural state is in  (so one aspect of culture/interest being considered), we see the network dividing into two subgroups with agents converging to states that are symmetric with respect to the time axis.
One can think of situations like a large urban environment in which there exist communities each with a strong cultural identity.In such a case, agents within the community seek to retain their cultural diversity.Thus, it is expected that a large relative cultural affinity between agents of different communities, but a small relative cultural affinity between agents within the same community, is expected to exist.

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 con-sisting of 50 members with parameters a = 0.5, q = 0.04, b = 0.41, and c = 2. Further, we note that in this case, 1 In this example, the conditions for the invariant sets given in Theorems 1 and 2 in Section (4) are satisfied.Hence, for i z such that 0.5 i z ≤ , it is the case that after some time, 0.5 i z ≥ ; that is, the member does not move towards the center of the network.
Further, for i z such that 2.3 i z ≤ , after some time, 0.5 2 is a plot of the approximate solutions for the full membership of the In order to make the dynamics of the network clearer, Figure 3 is a plot of the approximate solution of Equation ( 6) for six of the members of the network.
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 to above, we have plotted the approximate solution for the full network in Figure 4 and the approximate solution for the same six members as Figure 3 in Figure 5.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 i z such that 0.22 i z ≥ , the member does not move towards the center of the network and for i z such that 0.84 i z ≥ , after some time, 0.22 5 i z ≤ ≤ .Similar to above we have plotted the approximate solution for the full network in Figure 6 and the approximate solution for the same six members in Figure 7.

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).  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 ij x and not just the cultural affinity between an agent and the center of the network.Our hopes are to better understand what factors may lead to preserving a lower bound on the relative cultural affinity ij x that is strictly greater than zero as t → ∞ .In modeling such a network, we are looking to better understand how diversity between all members may be maintained over the long term within a culturally diverse network.Further, we are exploring the effect of noise on the network by considering similar dynamics and stochastic differential equations.The goal for both such questions is to better understand the impact of perturbations/impulses, both internal and external, on the behavior and diversity of multi-cultural networks.
sin i z is the sine-cyclical influence due to the 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" 23) reduces to the following differential inequality:

Figure 1 .
Figure 1.An example in 2 of the sets defined in Equation (36).Under the assumptions in Equation (34), the sets form concentric annuli.
are both monotonically decreasing and approaching to ( ) 2 am and ( ) 2 a qm respectively.This yields

∈
. It can be shown that there exists a neighborhood, ( ) individual member cultural state is pushed out/re- pulsed from the cultural state center x at some time T depending on 0 >  .

.
Then the solution to the lower comparison equation grows as t grows and approaches asymptotically to the threshold limit ( ) 2 a qm from below resulting in stronger ties with the community center state, x .If the initial value of the lower comparison equation is such that Turning to the upper comparison equation, we can consider the behavior of the transformed network members whose initial positions are in the sets D and E .Let 0 r be the initial position of the solution ( )

Figure 2 .
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 .
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 .
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 .
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.

Figure 6 .
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 .
Figure 7. Euler approximation of the solution to the differential equation given by Equation(6) we note that