A soluble model for synchronized rhythmic activity in ant colonies.

Synchronization is one of the most striking instances of collective behavior, occurring in many natural phenomena. For example, in some ant species, ants are inactive within the nest most of the time, but their bursts of activity are highly synchronized and involve the entire nest population. Here we revisit a simulation model that generates this synchronized rhythmic activity through autocatalytic behavior, i.e., active ants can activate inactive ants, followed by a period of rest. We derive a set of delay differential equations that provide an accurate description of the simulations for large ant colonies. Analysis of the fixed-point solutions, complemented by numerical integration of the equations, indicates the existence of stable limit-cycle solutions when the rest period is greater than a threshold and the event of spontaneous activation of inactive ants is very unlikely, so that most of the arousal of ants is done by active ants. Furthermore, we argue that the persistent oscillations observed in the simulations for colonies of finite size are due to resonant amplification of demographic noise.

The dynamics of casual groups can keep free-riders at bay.

Understanding the conditions for maintaining cooperation in groups of unrelated individuals despite the presence of non-cooperative members is a major research topic in contemporary biological, sociological, and economic theory. The N-person snowdrift game models the type of social dilemma where cooperative actions are costly, but there is a reward for performing them. We study this game in a scenario where players move between play groups following the casual group dynamics, where groups grow by recruiting isolates and shrink by losing individuals who then become isolates. This describes the size distribution of spontaneous human groups and also the formation of sleeping groups in monkeys. We consider three scenarios according to the probability of isolates joining a group. We find that for appropriate choices of the cost-benefit ratio of cooperation and the aggregation-disaggregation ratio in the formation of casual groups, free-riders can be completely eliminated from the population. If individuals are more attracted to large groups, we find that cooperators persist in the population even when the mean group size diverges. We also point out the remarkable similarity between the replicator equation approach to public goods games and the trait group formulation of structured demes.

A stochastic model for the influence of social distancing on loneliness.

The short-term economic consequences of the critical measures employed to curb the transmission of Covid-19 are all too familiar, but the consequences of isolation and loneliness resulting from those measures on the mental well-being of the population and their ensuing long-term economic effects are largely unknown. Here we offer a stochastic agent-based model to investigate social restriction measures in a community where the feelings of loneliness of the agents dwindle when they are socializing and grow when they are alone. In addition, the intensity of those feelings, which are measured by a real variable that we term degree of loneliness, determines whether the agent will seek social contact or not. We find that decrease of the number, quality or duration of social contacts lead the community to enter a regime of burnout in which the degree of loneliness diverges, although the number of lonely agents at a given moment amounts to only a fraction of the total population. This regime of mental breakdown is separated from the healthy regime, where the degree of loneliness is finite, by a continuous phase transition. We show that the community dynamics is described extremely well by a simple mean-field theory so our conclusions can be easily verified for different scenarios and parameter settings. The appearance of the burnout regime illustrates neatly the side effects of social distancing, which give to many of us the choice between physical infection and mental breakdown.

The paradox of productivity during quarantine: an agent-based simulation.

ABSTRACT: Economies across the globe were brought to their knees due to lockdowns and social restriction measures to contain the spread of the SARS-CoV-2, despite the quick switch to remote working. This downfall may be partially explained by the "water cooler effect", which holds that higher levels of social interaction lead to higher productivity due to a boost in people's mood. Somewhat paradoxically, however, there are reports of increased productivity in the remote working scenario. Here we address quantitatively this issue using a variety of experimental findings of social psychology that address the interplay between mood, social interaction and productivity to set forth an agent-based model for a workplace composed of extrovert and introvert agent stereotypes that differ solely on their propensities to initiate a social interaction. We find that the effects of curtailing social interactions depend on the proportion of the stereotypes in the working group: while the social restriction measures always have a negative impact on the productivity of groups composed predominantly of introverts, they may actually improve the productivity of groups composed predominantly of extroverts. Our results offer a proof of concept that the paradox of productivity during quarantine can be explained by taking into account the distinct effects of the social distancing measures on extroverts and introverts.

Influence of technological progress and renewability on the sustainability of ecosystem engineers populations.

Overpopulation and environmental degradation due to inadequate resource-use are outcomes of human's ecosystem engineering that has profoundly modified the world's landscape. Despite the age-old concern that unchecked population and economic growth may be unsustainable, the prospect of societal collapse remains contentious today. Contrasting with the usual approach to modeling human-nature interactions, which are based on the Lotka-Volterra predator-prey model with humans as the predators and nature as the prey, here we address this issue using a discrete-time population dynamics model of ecosystem engineers. The growth of the population of engineers is modeled by the Beverton-Holt equation with a density-dependent carrying capacity that is proportional to the number of usable habitats. These habitats (e.g., farms) are the products of the work of the individuals on the virgin habitats (e.g., native forests), hence the denomination engineers of ecosystems to those agents. The human-made habitats decay into degraded habitats, which eventually regenerate into virgin habitats. For slow regeneration resources, we find that the dynamics is dominated by rounds of prosperity and collapse, in which the population reaches vanishing small densities. However, increase of the efficiency of the engineers to explore the resources eliminates the dangerous oscillatory patterns of feast and famine and leads to a stable equilibrium that balances population growth and resource availability. This finding supports the viewpoint of growth optimists that technological progress may avoid collapse.

Mobility helps problem-solving systems to avoid groupthink.

Groupthink occurs when everyone in a group starts thinking alike, as when people put unlimited faith in a leader. Avoiding this phenomenon is a ubiquitous challenge to problem-solving enterprises and typical countermeasures involve the mobility of group members. Here we use an agent-based model of imitative learning to study the influence of the mobility of the agents on the time they require to find the global maxima of NK-fitness landscapes. The agents cooperate by exchanging information on their fitness and use this information to copy the fittest agent in their influence neighborhoods, which are determined by face-to-face interaction networks. The influence neighborhoods are variable since the agents perform random walks in a two-dimensional space. We find that mobility is slightly harmful for solving easy problems, i.e., problems that do not exhibit suboptimal solutions or local maxima. For difficult problems, however, mobility can prevent the imitative search being trapped in suboptimal solutions and guarantees a better performance than the independent search for any system size.

Agent-based models of collective intelligence.

Collective or group intelligence is manifested in the fact that a team of cooperating agents can solve problems more efficiently than when those agents work in isolation. Although cooperation is, in general, a successful problem solving strategy, it is not clear whether it merely speeds up the time to find the solution, or whether it alters qualitatively the statistical signature of the search for the solution. Here we review and offer insights on two agent-based models of distributed cooperative problem-solving systems, whose task is to solve a cryptarithmetic puzzle. The first model is the imitative learning search in which the agents exchange information on the quality of their partial solutions to the puzzle and imitate the most successful agent in the group. This scenario predicts a very poor performance in the case imitation is too frequent or the group is too large, a phenomenon akin to Groupthink of social psychology. The second model is the blackboard organization in which agents read and post hints on a public blackboard. This brainstorming scenario performs the best when there is a stringent limit to the amount of information that is exhibited on the board. Both cooperative scenarios produce a substantial speed up of the time to solve the puzzle as compared with the situation where the agents work in isolation. The statistical signature of the search, however, is the same as that of the independent search.

The spatial dynamics of ecosystem engineers.

The changes on abiotic features of ecosystems have rarely been taken into account by population dynamics models, which typically focus on trophic and competitive interactions between species. However, understanding the population dynamics of organisms that must modify their habitats in order to survive, the so-called ecosystem engineers, requires the explicit incorporation of abiotic interactions in the models. Here we study a model of ecosystem engineers that is discrete both in space and time, and where the engineers and their habitats are arranged in patches fixed to the sites of regular lattices. The growth of the engineer population is modeled by Ricker equation with a density-dependent carrying capacity that is given by the number of modified habitats. A diffusive dispersal stage ensures that a fraction of the engineers move from their birth patches to neighboring patches. We find that dispersal influences the metapopulation dynamics only in the case that the local or single-patch dynamics exhibit chaotic behavior. In that case, it can suppress the chaotic behavior and avoid extinctions in the regime of large intrinsic growth rate of the population.

Impact of centrality on cooperative processes.

The solution of today's complex problems requires the grouping of task forces whose members are usually connected remotely over long physical distances and different time zones. Hence, understanding the effects of imposed communication patterns (i.e., who can communicate with whom) on group performance is important. Here we use an agent-based model to explore the influence of the betweenness centrality of the nodes on the time the group requires to find the global maxima of NK-fitness landscapes. The agents cooperate by broadcasting messages, informing on their fitness to their neighbors, and use this information to copy the more successful agents in their neighborhood. We find that for easy tasks (smooth landscapes), the topology of the communication network has no effect on the performance of the group, and that the more central nodes are the most likely to find the global maximum first. For difficult tasks (rugged landscapes), however, we find a positive correlation between the variance of the betweenness among the network nodes and the group performance. For these tasks, the performances of individual nodes are strongly influenced by the agents' dispositions to cooperate and by the particular realizations of the rugged landscapes.

Effect of long-range interactions on the phase transition of Axelrod's model.

Axelrod's model with F=2 cultural features, where each feature can assume k states drawn from a Poisson distribution of parameter q, exhibits a continuous nonequilibrium phase transition in the square lattice. Here we use extensive Monte Carlo simulations and finite-size scaling to study the critical behavior of the order parameter ρ, which is the fraction of sites that belong to the largest domain of an absorbing configuration averaged over many runs. We find that it vanishes as ρ∼(q_{c}^{0}-q)^{ß} with ß≈0.25 at the critical point q_{c}^{0}≈3.10 and that the exponent that measures the width of the critical region is ν^{0}≈2.1. In addition, we find that introduction of long-range links by rewiring the nearest-neighbors links of the square lattice with probability p turns the transition discontinuous, with the critical point q_{c}^{p} increasing from 3.1 to 27.17, approximately, as p increases from 0 to 1. The sharpness of the threshold, as measured by the exponent ν^{p}≈1 for p>0, increases with the square root of the number of nodes of the resulting small-world network.

Influence of network topology on cooperative problem-solving systems.

The idea of a collective intelligence behind the complex natural structures built by organisms suggests that the organization of social networks is selected so as to optimize problem-solving competence at the group level. Here we study the influence of the social network topology on the performance of a group of agents whose task is to locate the global maxima of NK fitness landscapes. Agents cooperate by broadcasting messages informing on their fitness and use this information to imitate the fittest agent in their influence networks. In the case those messages convey accurate information on the proximity of the solution (i.e., for smooth fitness landscapes), we find that high connectivity as well as centralization boosts the group performance. For rugged landscapes, however, these characteristics are beneficial for small groups only. For large groups, it is advantageous to slow down the information transmission through the network to avoid local maximum traps. Long-range links and modularity have marginal effects on the performance of the group, except for a very narrow region of the model parameters.

Phenotypic plasticity, the Baldwin effect, and the speeding up of evolution: the computational roots of an illusion.

An increasing number of dissident voices claim that the standard neo-Darwinian view of genes as 'leaders' and phenotypes as 'followers' during the process of adaptive evolution should be turned on its head. This idea is older than the rediscovery of Mendel's laws of inheritance, with the turn-of-the-twentieth-century notion eventually labeled as the 'Baldwin effect' as one of the many ways in which the standard neo-Darwinian view can be turned around. A condition for this effect is that environmentally induced variation such as phenotypic plasticity or learning is crucial for the initial establishment of a trait. This gives the additional time for natural selection to act on genetic variation and the adaptive trait can be eventually encoded in the genotype. An influential paper published in the late 1980s claimed the Baldwin effect to happen in computer simulations, and avowed that it was crucial to solve a difficult adaptive task. This generated much excitement among scholars in various disciplines that regard neo-Darwinian accounts to explain the evolutionary emergence of high-order phenotypic traits such as consciousness or language almost hopeless. Here, we use analytical and computational approaches to show that a standard population genetics treatment can easily crack what the scientific community has granted as an unsolvable adaptive problem without learning. Evolutionary psychologists and linguists have invoked the (claimed) Baldwin effect to make wild assertions that should not be taken seriously. What the Baldwin effect needs are plausible case-histories.

Imitative learning as a connector of collective brains.

The notion that cooperation can aid a group of agents to solve problems more efficiently than if those agents worked in isolation is prevalent in computer science and business circles. Here we consider a primordial form of cooperation - imitative learning - that allows an effective exchange of information between agents, which are viewed as the processing units of a social intelligence system or collective brain. In particular, we use agent-based simulations to study the performance of a group of agents in solving a cryptarithmetic problem. An agent can either perform local random moves to explore the solution space of the problem or imitate a model agent - the best performing agent in its influence network. There is a trade-off between the number of agents N and the imitation probability p, and for the optimal balance between these parameters we observe a thirtyfold diminution in the computational cost to find the solution of the cryptarithmetic problem as compared with the independent search. If those parameters are chosen far from the optimal setting, however, then imitative learning can impair greatly the performance of the group.

Effect of migration in a diffusion model for template coexistence in protocells.

The compartmentalization of distinct templates in protocells and the exchange of templates between them (migration) are key elements of a modern scenario for prebiotic evolution. Here we use the diffusion approximation of population genetics to study analytically the steady-state properties of such a prebiotic scenario. The coexistence of distinct template types inside a protocell is achieved by a selective pressure at the protocell level (group selection) favoring protocells with a mixed template composition. In the degenerate case, where the templates have the same replication rate, we find that a vanishingly small migration rate suffices to eliminate the segregation effect of random drift and so to promote coexistence. In the nondegenerate case, a small migration rate greatly boosts coexistence as compared with the situation where there is no migration. However, increase of the migration rate beyond a critical value leads to the complete dominance of the more efficient template type (homogeneous regime). In this case, we find a continuous phase transition separating the homogeneous and the coexistence regimes, with the order parameter vanishing linearly with the distance to the transition point.

Nonlinear group survival in Kimura's model for the evolution of altruism.

Establishing the conditions that guarantee the spreading or the sustenance of altruistic traits in a population is the main goal of intergroup selection models. Of particular interest is the balance of the parameters associated to group size, migration and group survival against the selective advantage of the non-altruistic individuals. Here we use Kimura's diffusion model of intergroup selection to determine those conditions in the case the group survival rate is a nonlinear non-decreasing function of the proportion of altruists in a group. In the case this function is linear, there are two possible steady states which correspond to the non-altruistic and the altruistic phases. At the discontinuous transition line separating these phases there is a non-ergodic coexistence phase. For a continuous concave survival function, we find an ergodic coexistence phase that occupies a finite region of the parameter space in between the altruistic and the non-altruistic phases, and is separated from these phases by continuous transition lines. For a convex survival function, the coexistence phase disappears altogether but a bistable phase appears for which the choice of the initial condition determines whether the evolutionary dynamics leads to the altruistic or the non-altruistic steady state.

Reinforcement and inference in cross-situational word learning.

Cross-situational word learning is based on the notion that a learner can determine the referent of a word by finding something in common across many observed uses of that word. Here we propose an adaptive learning algorithm that contains a parameter that controls the strength of the reinforcement applied to associations between concurrent words and referents, and a parameter that regulates inference, which includes built-in biases, such as mutual exclusivity, and information of past learning events. By adjusting these parameters so that the model predictions agree with data from representative experiments on cross-situational word learning, we were able to explain the learning strategies adopted by the participants of those experiments in terms of a trade-off between reinforcement and inference. These strategies can vary wildly depending on the conditions of the experiments. For instance, for fast mapping experiments (i.e., the correct referent could, in principle, be inferred in a single observation) inference is prevalent, whereas for segregated contextual diversity experiments (i.e., the referents are separated in groups and are exhibited with members of their groups only) reinforcement is predominant. Other experiments are explained with more balanced doses of reinforcement and inference.

Effect of external fields in Axelrod's model of social dynamics.

The study of the effects of spatially uniform fields on the steady-state properties of Axelrod's model has yielded plenty of counterintuitive results. Here, we reexamine the impact of this type of field for a selection of parameters such that the field-free steady state of the model is heterogeneous or multicultural. Analyses of both one- and two-dimensional versions of Axelrod's model indicate that the steady state remains heterogeneous regardless of the value of the field strength. Turning on the field leads to a discontinuous decrease on the number of cultural domains, which we argue is due to the instability of zero-field heterogeneous absorbing configurations. We find, however, that spatially nonuniform fields that implement a consensus rule among the neighborhood of the agents enforce homogenization. Although the overall effects of the fields are essentially the same irrespective of the dimensionality of the model, we argue that the dimensionality has a significant impact on the stability of the field-free homogeneous steady state.

A structural model of emotions of cognitive dissonances.

Cognitive dissonance is the stress that comes from holding two conflicting thoughts simultaneously in the mind, usually arising when people are asked to choose between two detrimental or two beneficial options. In view of the well-established role of emotions in decision making, here we investigate whether the conventional structural models used to represent the relationships among basic emotions, such as the Circumplex model of affect, can describe the emotions of cognitive dissonance as well. We presented a questionnaire to 34 anonymous participants, where each question described a decision to be made among two conflicting motivations and asked the participants to rate analogically the pleasantness and the intensity of the experienced emotion. We found that the results were compatible with the predictions of the Circumplex model for basic emotions.

Statistics of opinion domains of the majority-vote model on a square lattice.

The existence of juxtaposed regions of distinct cultures in spite of the fact that people's beliefs have a tendency to become more similar to each other's as the individuals interact repeatedly is a puzzling phenomenon in the social sciences. Here we study an extreme version of the frequency-dependent bias model of social influence in which an individual adopts the opinion shared by the majority of the members of its extended neighborhood, which includes the individual itself. This is a variant of the majority-vote model in which the individual retains its opinion in case there is a tie among the neighbors' opinions. We assume that the individuals are fixed in the sites of a square lattice of linear size L and that they interact with their nearest neighbors only. Within a mean-field framework, we derive the equations of motion for the density of individuals adopting a particular opinion in the single-site and pair approximations. Although the single-site approximation predicts a single opinion domain that takes over the entire lattice, the pair approximation yields a qualitatively correct picture with the coexistence of different opinion domains and a strong dependence on the initial conditions. Extensive Monte Carlo simulations indicate the existence of a rich distribution of opinion domains or clusters, the number of which grows with L2 whereas the size of the largest cluster grows with ln L2. The analysis of the sizes of the opinion domains shows that they obey a power-law distribution for not too large sizes but that they are exponentially distributed in the limit of very large clusters. In addition, similarly to other well-known social influence model--Axelrod's model--we found that these opinion domains are unstable to the effect of a thermal-like noise.

Social interaction as a heuristic for combinatorial optimization problems.

We investigate the performance of a variant of Axelrod's model for dissemination of culture--the Adaptive Culture Heuristic (ACH)--on solving an NP-Complete optimization problem, namely, the classification of binary input patterns of size F by a Boolean Binary Perceptron. In this heuristic, N agents, characterized by binary strings of length F which represent possible solutions to the optimization problem, are fixed at the sites of a square lattice and interact with their nearest neighbors only. The interactions are such that the agents' strings (or cultures) become more similar to the low-cost strings of their neighbors resulting in the dissemination of these strings across the lattice. Eventually the dynamics freezes into a homogeneous absorbing configuration in which all agents exhibit identical solutions to the optimization problem. We find through extensive simulations that the probability of finding the optimal solution is a function of the reduced variable F/N(») so that the number of agents must increase with the fourth power of the problem size, N∝F4, to guarantee a fixed probability of success. In this case, we find that the relaxation time to reach an absorbing configuration scales with F6 which can be interpreted as the overall computational cost of the ACH to find an optimal set of weights for a Boolean binary perceptron, given a fixed probability of success.