Paradox of integration—Cellular automata approach

Abstract

We discuss the self-deprecating strategy introduced by Peter Blau as one of stages of the process of social integration. Recently we have introduced a two-dimensional space of status, real and surface one ($A$ and $B$), and we have demonstrated that with this setup, the self-deprecating strategy efficiently prevents the rejection K. Malarz and K. Kułakowski, International Journal of Modern Physics C 30, 1950040 (2019). There, the process of reducing the conflict was described by master equations, i.e. a set of differential equations describing evolution of density $v(A,B)$ of actors of status $(A,B)$. Here we reformulate the problem in terms of probabilistic cellular automata. The obtained results for number $n(A,B)$ of actors of status $(A,B)$ are qualitatively the same as in the previous approach, both for synchronous and asynchronous version of the automaton. Namely, an enhancement of the surface status compensates a deficiency of the real one.

I Introduction

According to Peter Blau, social integration of a group of adults contains two stages Blau 1964. In first stage, actors demonstrate their strongest points, to achieve social status as high as possible. This stage can be painful for some, who have no virtues to present; they respond with fear and hostility. As Blau puts it: “The more successful A is in impressing B and earning B’s high regard, the more displeasure he causes to C, whose relative standing in the eyes of B has suffered. All group members simultaneously play the role of A, B, and C in this schema, which greatly complicates the competitive process” (Blau 1964, p. 44). This leads to a paradoxical phenomenon: persons most skilful, intelligent and physically attractive meet with rejection and hostility. To neutralize this, in the second stage intelligent persons know and apply a clever strategy: they seemingly reduce their advantage, demonstrating their weak points in less important aspects of status. To cite Blau again: “Having first impressed us with his Harvard accent and Beacon Hill friends, he may later tell a story that reveals his immigrant background” (Blau 1964, p. 48). It is worth to mention that this strategy is not efficient if applied by an actor too weak or too strong. In first case, it is read as a fake suggestion of non-existent strong points (Blau 1964, pp. 48-49), in the second—as an arrogant demonstration of lack of understanding of difficulties of the others’ life Benen June 16, 2011. As a fairly complex social process, the social integration is a promising playground for an interdisciplinary research, between social sciences and computational modelling.

The self-deprecating strategy (SDS) depicted by Blau has been the subject of a series of recent papers Krawczyk and Kułakowski 2017; Kułakowski et al. 2017; Malarz and Kułakowski 2019. Yet, it is only in the last Malarz and Kułakowski 2019 where the model results successfully reproduced the efficiency of SDS as a tool to reduce the fear-driven rejection. This was achieved by an introduction of a two-dimensional space of status, with real and surface axes. Along these two axes, two processes were competing: fear-induced rejection equivalent to a shift down along the real axis, and preventive praising, which drives the status of an opponent up along the surface axis. The probabilities of these two strategies ($\alpha$ and $1-\alpha$, respectively) have been used as parameters; below we keep the same notation. The formalism applied was a set of differential equations, with both probability distribution of agents in the space of status and the related cumulative distribution involved as variables. In this sense, the description presented in Ref. Malarz and Kułakowski 2019 was non-local; the time evolution of the status of actors depended not only on their direct neighbours, but also on those fairly distant on the status plane.

To which extent social interactions are active between individuals of clearly different social status depends on the context. Excerpts from the Blau book Blau 1964 given above indicate a living room or another informal gathering where new acquaintances are made. Even there, the answer relies on local culture (Weber 1978, p. 1401), Goody 1971; Ridgeway and Smith-Lovin 1999. On the other hand, both in Ref. Malarz and Kułakowski 2019 and here we have in mind SDS applied in a scale of groups and not only individuals. One of examples quoted in Ref. Malarz and Kułakowski 2019 is the glorification of working class in communist countries, attributing the role of dictators to proletariat Marx and Engels 1969 [1848]; the scale of this manipulation was transnational. Yet, it is obvious that in this and similar cases SDS is expected to be more efficient when performed by groups of status at least nominally close to the one of the target group. In particular, the Soviet Politburo composed of aristocrats would be much less credible for working classes.

Here we are interested on the efficiency of SDS if applied only to nearest neighbours in the space of status. According to the note in the first paragraph, if the difference of statuses of two individual is large the mechanism is less effective. Our aim is to check, how the results of the simulation depend on the assumption on the local character of interactions. To achieve this purpose, the problem is reformulated in the frames of probabilistic cellular automata, the formalism local by definition Adamatzky 2018; *Wolfram-2002; *Ilachinski-2001; *Hegselmann-2000; *Hegselmann-1998. Here, the time evolution of the positions of actors in the status plane depends only on their direct neighbours in this plane. Besides this, we keep the time-dependent distribution of actors in the status space as the variable, as was done in Ref. Malarz and Kułakowski 2019.

Table 1: Videos showing the status–temporal evolution of the number $n(A,B)$ of actors in status $(A,B)$ for synchronous and asynchronous sites update and various values of $\alpha$, $L_{A}=40$, $L_{B}=60$. The final system states are presented in Figs. 1 and 2.

$\alpha$	URL
synchronous update
0.00	http://www.zis.agh.edu.pl/files/synchro000.gif
0.25	http://www.zis.agh.edu.pl/files/synchro025.gif
0.50	http://www.zis.agh.edu.pl/files/synchro050.gif
0.75	http://www.zis.agh.edu.pl/files/synchro075.gif
1.00	http://www.zis.agh.edu.pl/files/synchro100.gif
asynchronous update
0.00	http://www.zis.agh.edu.pl/files/asynchro000.gif
0.25	http://www.zis.agh.edu.pl/files/asynchro025.gif
0.50	http://www.zis.agh.edu.pl/files/asynchro050.gif
0.75	http://www.zis.agh.edu.pl/files/asynchro075.gif
1.00	http://www.zis.agh.edu.pl/files/asynchro100.gif

The automaton rule is checked both for synchronous (parallel) and asynchronous (sequential) version. This is done because we can expect that for intermediate values of the coefficient $\alpha$ the results of the simulations within synchronous and asynchronous scheme are different. This expectation is in agreement with literature. Since publication of Blok and Bergersen 1999 paper, we are aware that the updating scheme influence the results of simulation based on cellular automata technique. For famous Conway’s “Game of life” Berlekamp et al. 1982 automaton changing synchronous to asynchronous updating scheme results in final picture of the lattice similar to the maze instead of well-known lattice of structures with density of life not exceeding 3% Bagnoli et al. 1991; Malarz et al. 1998. Also for Ising model the application parallel or sequential spins updating leads to different results Skorupa et al. 2012. As Skorupa et al. 2012 put it: “the problem of updating methods is widely discussed in a recent work on cellular automata, Boolean networks, neural networks, and the so-called agent-based modelling in ecology and sociology Bandini et al. 2012; Caron-Lormier et al. 2008; Castellano et al. 2009. It has been shown that the updating scheme can have an enormous influence on the model output Cornforth et al. 2005.” The list of examples of papers devoted to difference between synchronous and asynchronous scheme of sites updating may be extended further, for instance for paper published in the journal devoted to biosystems Schönfisch and de Roos 1999. Also in Ref. Skorupa et al. 2012 authors speculate, that differences in system behaviours in sequential or synchronous updating may depend on existence (or not) of the equilibrium, i.e. on satisfying (or not) the detailed balance condition by the updating rules. As automaton rules are asymmetric—the flow of actors may be from bottom to top (Eq. (1)) and from right to left (Eq. (2))—the detailed balance conditions are violated and thus the order of site updates may influence the results. These arguments incline us to check how the results of our simulations depend on the details of updating the cell states.

In the next section, the automaton is described in detail. Third section is devoted to our numerical results, presented in form of computer animations. The same way of presentation was used in Ref. Malarz and Kułakowski 2019. A short discussion is closing the text.

II Model

Figure 1: The distribution of the final number $n(A,B)$ of actors with status $(A,B)$ after $t$ time steps and various values of $\alpha$, $L_{A}=40$, $L_{B}=60$ and synchronous sites update.

Two processes are competing in the time evolution. First is the fear-driven rejection, which takes place by a lowering of status $A$ (real status) of an actor by his direct neighbour with lower status $A$ and the same status $B$ (surface status). This process takes place with probability $\alpha$. The second process is an enhancement of status $B$ (surface status) of an actor by his direct neighbour with higher status $A$ and the same status $B$. This process takes place with probability $1-\alpha$. Actors with different status $B$ do not interact. The latter rule assures, that the action to increase the status $B$ of a neighbour neutralises his rejection and allows to preserve own status $A$. This rule is not possible if we have only one-dimensional status, as was assumed earlier Krawczyk and Kułakowski 2017; Kułakowski et al. 2017.

Every site of rectangle lattice $\mathcal{G}=\{(A,B):1\leq A\leq L_{A},1\leq B\leq L_{B}\}$ represents the number $n(A,B)$ of actors with real status $A$ and surface status $B$. Initially, all sites for $B\leq 15$ are occupied by twenty actors. Every time step $t$, every pair of actors at position $\{(A,B)\cup(A+1,B)\}$ for which $n(A,B;t)n(A+1,B;t)>0$ apply either SDS (with probability of $1-\alpha$):

$$n(A,B+1;t+1)=n(A,B+1;t)+1,$$ (1a)

$$n(A,B;t+1)=n(A,B;t)-1$$ (1b)

or the fear-driven rejection process (with probability $\alpha$):

$$n(A+1,B;t+1)=n(A+1,B;t)-1,$$ (2a)

$$n(A,B;t+1)=n(A,B;t)+1.$$ (2b)

In the asynchronous version, one time step is equivalent to an update of all pairs in random order. In both versions, the random number (to apply SDS or not) is selected for each pair separately.

III Results

In the upper part of Tab. 1 the links to videos showing status-time system evolution obtained with parallel update scheme are provided.

The final states of system evolution for various values of probabilities $\alpha$ and synchronous sites update are shown in Figs. 1(b)–(f). In Fig. 1(a) the common initial state of the system is presented. For $\alpha=0$ (Fig. 1(b)) all actors apply SDS, which yields generation of plenty actors with high surface status $B$. On contrary, assuming $\alpha=1$ pushes system to the final state presented in Fig. 1(f) with all actors with minimal real status $A=1$ and the surface status $B$ the same as initial. We note that the system evolution for $\alpha=0$ and $\alpha=1$ presented in Fig. 1(b) and 1(f) are common for both schemes discussed in this paper.

Figure 2: The distribution of the final number $n(A,B)$ of actors with status $(A,B)$ after $t$ time steps and various values of $\alpha$, $L_{A}=40$, $L_{B}=60$ and asynchronous sites update.

In Figs. 1(b)–(f) we see that the final outcome depends on the probability $\alpha$; indeed, the asymptotic (stationary) state depends only on one parameter $\alpha/(1-\alpha)$. For $\alpha$ close to one, SDS is not active, and the rejection reduces the real status A of all actors to the minimal value. In the opposite limit $\alpha$ close to zero, a straight line appears in the plane $(A,B)$, with a slope $\pi/4$ to both axes. This result is the same as the one obtained in Ref. Malarz and Kułakowski 2019. As we argued there, the angle has no real mining, because it is the consequence of the assumed scales of both statuses, both arbitrary and unverifiable.

In Figs. 2(b)–(f) the results of simulations for the asynchronous version of sites update are presented. As we see, the results for both schemes are the same up to details due to the randomness of the system. In both cases we observe characteristic chess-board-like pattern of neighboring cells. These configurations remain unchanged during the evolution; this is due to the rule that cells with different status $B$ do not interact anymore.

IV Discussion

The results on the number $n(A,B)$ of actors of status $(A,B)$ obtained within the synchronous (parallel) and asynchronous schemes are qualitatively the same as those produced by the differential equations Malarz and Kułakowski 2019. Namely, an enhancement of the surface status $B$ compensates a deficiency of the real one $A$. It appears that the localness is not crucial. This result is consistent with our previous tests Malarz and Kułakowski 2019, where the interaction strength was assumed to decrease with the difference of statuses. Both the present calculations for local interactions and the tests mentioned above take into account the hint by Blau Blau 1964, that the interaction strength decreases with the distance in the status space. Yet, the final results remain not influenced by the localness of the automaton rules.

As we mentioned in the Introduction, due to an asymmetry of automaton rules, the flow of actors is directional: actors move only downwards along $A$ axis and upwards along $B$ axis. This assumption leads to vanishing the fittest part of the population, i.e. those with the highest real status $A=L_{x}$. This phenomenon is described in the literature of genetics as an error catastrophe. The error catastrophe is the first step in the so-called Muller’s ratchet Muller 1964 observed also in the Eigen’s quasispecies Eigen et al. 1989; *Tiggemann1998. Next, also the secondary most fitted (i.e. with real status $A=L_{x}-1$) may disappear, what is the second step in the Muller’s ratchet, etc. Due to the asymmetry of automaton rules the fraction with given $A$, once vanished, will never appear again. As we can see in Fig. 1(b) only taking by actors the pure SDS strategy (for $\alpha=0$) may prevent the error catastrophe. On the other hand, for $\alpha=1$ all but the last available step of Muller’s ratchet take place when SDS is avoided. On the contrary to genetics, here the effect is due to a conscious action of actors with low status A. Notwithstanding, for $\alpha>0$ the catastrophe is unavoidable (see Figs. 1(c)–(f), 2(c)–(f)).

Concluding, the outcome of the simulations indicates that a quick application of SDS by a smart actor blocks the fear-driven rejection. We can expect that in a consequence, SDS is simultaneously putting the subject actor in dependence of surface praising.