Ordering dynamics in the voter model with aging
Antonio F. Peralta, Nagi Khalil, Raul Toral

TL;DR
This paper investigates the voter model with aging, incorporating memory-dependent activation probabilities, and derives analytical conditions for different dynamical regimes including consensus, coexistence, and frozen states.
Contribution
It introduces a mean-field framework with age-dependent activation probabilities, providing analytical results for the system's behavior near absorbing states.
Findings
Coexistence of opinions occurs with increasing activation probability.
Consensus or frozen states depend on aging and activation probability limits.
Time to consensus can be exponential or power-law depending on parameters.
Abstract
The voter model with memory-dependent dynamics is theoretically and numerically studied at the mean-field level. The `internal age', or time an individual spends holding the same state, is added to the set of binary states of the population, such that the probability of changing state (or activation probability ) depends on this age. A closed set of integro-differential equations describing the time evolution of the fraction of individuals with a given state and age is derived, and from it analytical results are obtained characterizing the behavior of the system close to the absorbing states. In general, different age-dependent activation probabilities have different effects on the dynamics. When the activation probability is an increasing function of the age , the system reaches a steady state with coexistence of opinions. In the case of aging, with being a…
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11
Figure 12
Figure 13
Figure 14
Figure 15
Figure 16Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Ordering dynamics in the voter model with aging
Antonio F. Peralta
Nagi Khalil
Raúl Toral
IFISC (CSIC-UIB), Instituto de Física Interdisciplinar y Sistemas Complejos, Campus Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain
Abstract
The voter model with memory-dependent dynamics is theoretically and numerically studied at the mean-field level. The “internal age”, or time an individual spends holding the same state, is added to the set of binary states of the population, such that the probability of changing state (or activation probability ) depends on this age. A closed set of integro-differential equations describing the time evolution of the fraction of individuals with a given state and age is derived, and from it analytical results are obtained characterizing the behavior of the system close to the absorbing states. In general, different age-dependent activation probabilities have different effects on the dynamics. When the activation probability is an increasing function of the age , the system reaches a steady state with coexistence of opinions. In the case of aging, with being a decreasing function, either the system reaches consensus or it gets trapped in a frozen state, depending on the value of (zero or not) and the velocity of approaching . Moreover, when the system reaches consensus, the time ordering of the system can be exponential () or power-law like (). Exact conditions for having one or another behavior, together with the equations and explicit expressions for the exponents, are provided.
I Introduction
Agent-based binary-state models are commonly used to study the effect that simple individual dynamical rules may have on the collective behavior of a system. Typical examples include modeling the spreading of diseases Pastor-Satorras and Vespignani (2001); Pastor-Satorras et al. (2015), the time evolution of the number of speakers of a determined language Abrams and Strogatz (2003); Castelló et al. (2006); Stauffer et al. (2007); Vazquez et al. (2010), the evolution of prices in financial markets Kirman (1993); Alfarano et al. (2005, 2008); Alfarano and Milaković (2009); Gontis and Kononovicius (2017); Carro et al. (2015); Kononovicius and Ruseckas (2019), or the dynamics of opinion formation in a population of individuals Barrat et al. (2008); Castellano et al. (2009); Fernández-Gracia et al. (2014). As a prominent case of the latter, the voter model Clifford and Sudbury (1973); Holley and Liggett (1975); Vazquez and Eguíluz (2008) and its variants thereof Redner (2018), by which an agent adopts with some probability the state of a randomly chosen neighbor, are widely used when the mechanism at hand is that of imitation. Beyond its potential applications, the voter model has become a paradigm in the field of non-equilibrium statistical physics. It is defined by very simple rules, and explicit solutions can be given in some particular cases, yet it is capable of showing very rich behavior, including phase transitions, characterized by critical exponents, power-law correlations, finite-size scaling, tricritical behavior, etc. Krapivsky et al. (2010); Liggett (2012). Recent studies of the voter model include the effect of a network structure Peralta et al. (2018a); Carro et al. (2016), non-linear or group interactions Nyczka et al. (2012); Nyczka and Sznajd-Weron (2013); J\polhkedrzejewski (2017); Peralta et al. (2018b); Vieira and Anteneodo (2018); Raducha et al. (2018); Min and Miguel (2019), more than two states Herrerías-Azcué and Galla (2019); Vazquez et al. (2019), the effect of zealots Khalil et al. (2018) and contrarians Khalil and Toral (2019), etc.
A fundamental aspect of the agent-based binary-state models is the memory the constituents of the system have J\polhkedrzejewski and Sznajd-Weron (2018), and how this affects the dynamical rules that define the model. A clear example is the framework of non-Poissonian infection processes in epidemic spreading Min et al. (2013); Boguñá et al. (2014); Starnini et al. (2017), where the infection of an individual depends on the times each one of its neighbors became infected Blythe and Anderson (1988). This is because the rate at which one transmits a disease is not independent of the time one has been infected. In most cases of binary-state models a Markovian assumption is postulated which implies that the probability of changing state is independent of the past history of events of the individuals, that is to say neglecting memory effects by including Poissonian processes. This is of course a simplification that can be justified in many phenomena and it makes a lot easier the mathematical treatment van Kampen (1998); Łucza (2005), but it may be lacking of realistic features in other cases. In this context, one of the objectives of the present work is to provide a consistent mathematical description of the voter model including memory effects, namely by making the probability of changing opinion of agents to depend on their age or time elapsed until their last update Fernández-Gracia et al. (2011); Stark et al. (2008).
One of the main questions addressed by the voter model is whether the imitation rule leads to a situation of consensus, with all agents (or at least, a significant majority) adopting the same value of the binary state Ben-Naim et al. (1996); Sood and Redner (2005); Suchecki et al. (2005). The answer to this question is sometimes counter-intuitive. Although it would appear, for example, that an increase in the connections amongst the agents should favor consensus, a detailed analysis shows that for an effective dimensionality of the connectivity network greater than two, the system is not able to sustain a consensus state in the thermodynamic limit. The problem has been also addressed when including memory effects. Stark et. al Stark et al. (2008) introduced increasing inertia, otherwise known as “aging”, as a mechanism that reduces the probability of an agent to copy the state of another one. Aging can roughly be related to an increase of the inflexibility in the copying mechanism Martins and Galam (2013). Again, it would seem naively that reducing the number of interactions amongst agents should impede the reaching of consensus, but it was shown that exactly the opposite happens, namely the consensus is reached more easily in such an scenario of aging. In a recent paper by some of us Artime et al. (2018, 2018), it was shown that the inclusion of aging in the noisy version of the voter model Considine et al. (1989); Kirman (1993); Granovsky and Madras (1995) induces a state of imperfect consensus that remains in the thermodynamic limit. Another study of the noiseless voter model Fernández-Gracia et al. (2011) focused on the distribution of the time between individual changes of states , which in turn is related with the activation probability , or the age-dependent probability of interaction. The conclusion of Fernández-Gracia et al. (2011), sustained from the analysis of extensive numerical simulations, was that a particular functional form of induces a power-law dependence , which has been observed in real datasets, for example in human communications Wu et al. (2010); Candia et al. (2008). This power law is similar to the ones observed in the approach to consensus in low-dimensional lattices, but it shows a strong dependence on the details of the activation probability. The main objective of this paper is to analyze this problem and to derive the relations in Fernández-Gracia et al. (2011); Stark et al. (2008), among others, from an analytical point of view. We also aim to give a well constructed theoretical framework for aging, which has been studied together with network structure in several recent works Pérez et al. (2016); Artime et al. (2017a, b) but mainly computationally, due to difficulties encountered in the mathematical description.
The outline of the paper is as follows. In section II we introduce the voter model including memory effects. Two specific forms of the activation probability, rational and exponential functions of the age, are considered which depend on several parameters, so as to cover phenomenology observed in the literature. Section III includes a summary of the main results present in the literature about dynamical behavior of the system close to the absorbing states. In section IV the dynamical equations of the mean fraction of agents with a given state and age are obtained under a mean field description and its steady-state analysed. Section V contains the main results of the work. First, a closed set of integro-differential equations describing the system close to the absorbing states is derived. Second, an approximate analysis of the previous description is carried out, providing explanation to the numerical findings. Complementary results are given in Appendix A and Appendix B. Finally, section VI includes a summary of the main results.
II Model
The voter model implements in arguably the simplest way the herding mechanism for the evolution of binary state systems. In the original version one considers a system formed by individuals or agents. Agents are located in the nodes of a network and are connected by undirected, bidirectional, links. Two connected individuals are said to be “neighbors”. The network is single connected, meaning that every node can be reached from any other node by a sequence of links. Each individual holds a binary-state (spin) variable . Several interpretations can be given to this binary variable, but its exact meaning does not concern us in this paper. For example, the voter model or variants of it have been used to represent the optimistic/pessimistic state of a stock market broker Kirman (1993), the language A/B used by a speaker Abrams and Strogatz (2003); Vazquez et al. (2010), or the direction of the velocity right/left in a one-dimensional model of active particles Escaff et al. (2018). Those variables evolve over time by the following (stochastic) rules:
(i) An individual is selected at random amongst the possibilities.
(ii) The selected individual copies the state of another individual chosen also at random between the set of neighbors of .
In this sequential updating scheme, every time an individual is chosen for updating, time increases by one unit, while updates constitute one Monte Carlo step (MCS).
In the aging version of the voter model the above rules are modified such that the herding mechanism occurs only with an activation probability that depends on the internal age of the selected individual. The internal age of individual stands for the number of update attempts elapsed since its last change of state.
More explicitly, the model is reproduced in the simulations by modifying the last step as follows:
(ii) The selected individual copies with probability the state of another individual chosen at random between the set of neighbors of . If the selected individual changes state , then its internal age resets to zero ; otherwise it increases in one unit . Initially, all internal ages are set to zero.
The standard (no aging) voter model is recovered taking . Here, for the sake of concreteness and mathematical treatment, we have considered the following two functional forms for the activation probability:
- 1.-
The first functional form is a rational function of the age,
[TABLE]
where and are constants. If , the probability of interaction decreases with age, a typical aging situation. The opposite occurs if , as the activation probability increases with age. This can be interpreted as “anti-aging” or “rejuvenating” where as nodes get older they become more prone to change state. This aging form, with , is basically the form considered in Fernández-Gracia et al. (2011) and that has been shown to induce features observed in several real-word systems, such as power-law inter-event time distributions. If and , this reproduces qualitatively the form used in Ref. Stark et al. (2008). In that paper, the activation probability , equivalent to with being the “inertia”, decreases linearly from up to a final constant, non-null value , but at a finite value of .
- 2.-
The second form we adopt for mathematical simplicity is an exponential dependence of the age,
[TABLE]
with the same meaning as before for and , and is a parameter.
III Dynamical behavior. Absorbing states
When one individual changes its state after copying the state of one of its neighbors, it might occur that it decreases the total number of neighbors with which it shares the same state. Therefore, a question of interest is whether by iteration of the dynamical rules a “consensus”, or fully ordered, state in which all agents share the same state, either or , is finally achieved. One can argue that the stochastic rules permit eventually to reach any possible configuration of states. Therefore, when any of the particular consensus states or is reached, then there can not be further evolution and the dynamics stops. It seems, though, that the ultimate fate of the system is to reach consensus, with all nodes sharing the same value. While this is certainly true for any system with a finite number of agents, it is interesting to analyze the way the consensus state is reached and the dependence with system size of the average time to reach that state. Two different measures have been introduced to study the approach to order: the magnetization that takes values in the ordered states, and the density of active links, defined as the fraction of links that join nodes in different states. In the ordered states, the magnetization can be or , while in both absorbing states. The existence of absorbing states then ensures that and for a finite system.
There is a significant difference between the aging and non-aging versions of the model. In the non-aging version the approach to the absorbing state depends crucially on the spatial dimension . If it is observed that the density of active links 111A similar qualitative behavior is observed for the magnetization. In fact, for the all-to-all connectivity considered later in this paper, one has simply , although in other networks must be considered as an independent variable. has a fast transient to a plateau value (for , i.e. all to all connectivity, ) and it oscillates around it until a large fluctuation brings the system to the absorbing state at a time . This time is a stochastic variable whose mean value scales as for , see also Artime et al. (2018) where the probability distribution of is obtained for . Therefore, larger systems stay longer in an active state until a large fluctuation brings them to the absorbing state. As we take the thermodynamic limit , it is and the system does not order at all in a finite time. If we average over initial conditions and realizations of the dynamics, we find an exponential decay , with . As mentioned, in the thermodynamic limit when , we obtain , or absence of order. If we take first the limit and then the limit we obtain otherwise . If the spatial dimension is the decay to the absorbing state occurs very differently. One does not observe that fluctuates around a plateau value, but rather decreases on average for any realization even for large system sizes. As a consequence, in the thermodynamic limit, it is predicted that . The exact time dependence depends on the dimension as for and for Ben-Naim et al. (1996); Sood and Redner (2005); Suchecki et al. (2005).
The phenomenology in the aging version was described in Fernández-Gracia et al. (2011) for particular forms of the activation probability. It was observed that if follows a similar dependence as given by Eq. (1) with , then the approach to the absorbing states is described by a power law , for large spatial dimensions. Furthermore, in a fully-connected network it was observed numerically that the exponent of the power-law decay agreed within the numerical accuracy with the parameter of Eq. (1) 222Strictly speaking, the aforementioned power-law dependence was reported in Fernández-Gracia et al. (2011) for the cumulative inter-event time distribution . One of the authors of that reference (J. F. Gracia, private communication) has confirmed to us that the aforementioned asymptotic time dependence of was also observed, but not displayed.. This indicates that in the presence of aging, the system always orders, contrarily to the non-aging version. In this paper we offer an explanation of this behavior.
IV The dynamical equations
We now derive the mean-field dynamical equations of our stochastic process, as defined in section II. In the derivation we restrict ourselves to the all-to-all (or fully connected) network in which all nodes are neighbors. The mathematical description in the case of a more complex network structure in the interactions between nodes is a further complication Peralta et al. (2018a) and it is left for future studies Artime et al. (2017a, b). In the all-to-all setup, all the information needed to implement the stochastic update rules is contained in the set of the numbers of individuals with internal age in states , Stark et al. (2008). The global variables for the total number of up and down spins are and . Note, however, that not all variables of the state are independent, since . Hence, it is useful to consider an alternative representation of the system in terms of independent variables, for instance by using the variable and obviating , as . In this representation are given by and .
In order to derive the dynamical equations of the global state of the system we must include with their respective probabilities all events that induce changes in . Every time an individual changes its state or age, there is a change in the set of variables . The change can occur in different ways and with different probabilities. We now spell in detail the possibilities and their effects on the variables.
(1) Consider that at time the chosen individual is in state and has age and that as a result of the interaction it switches to and, hence, its internal age is reset to . The probability of this event is equal to the probability, , of choosing and individual with age and state , multiplied by the age-dependent probability, , that it activates the copying (herding) mechanism, and multiplied by the probability that the randomly selected neighbor is in the opposite state. Altogether, the probability is . When the switching of state occurs, if we have and , while if the only change is .
(2) This is similar to the previous case but now the chosen individual is initially in state . The probability of switching is . When the switching of state occurs, if we have and , while if the only change is .
(3) Consider that at time the chosen individual has and age , but that now it keeps its current state . This event happens with a probability equal to the probability of choosing an individual in state with age , , multiplied by the probability that it does not switch, which can arise either because the copying mechanism was not activated, with probability , or it was activated but the selected neighbor was in the same state, with probability . Altogether, the probability is . In this case the variables change as and if and if .
(4) Finally, we consider a similar case to the previous one but the chosen individual keeps its state. The switching probability is now . The changes in the state are , if and if .
If time is measured in MCS, the rates (probability per unit time) of the four possible processes are
[TABLE]
where and , thus , .
This derivation has considered a discrete time process in the limit of a large number of individuals. It is also possible to consider from the beginning a continuous time process in which an agent of age and state has a rate of switching state and set its internal age to [math], and a rate of increasing its internal age keeping the state; similarly an agent of age and state has a rate of switching state and set its internal age to [math], and a rate of increasing its internal age and keep the state.
Following standard techniques van Kampen (2007); Toral and Colet (2014); Peralta and Toral (2018), we could derive a master equation for the probability of having state at time , with the set of rates and changes in the variables explained before. As we restrict ourself to average values in this paper, we will only obtain the evolution equations for the ensemble average of the fraction of nodes with a given state and age , as well as for . The time evolution of the average values can be computed as a weighted sum of all the rates that contribute to its variation, the weights being the variation of the regarded variable in that process Peralta and Toral (2018). Using the standard mean-field approximation Toral and Colet (2014) which neglects correlations as we end up with a closed but infinite system of equations:
[TABLE]
Similar equations were obtained in Ref. Stark et al. (2008); Artime et al. (2018). In these equations, variables , whenever they appear, should be expressed in terms of the independent variables,
[TABLE]
The explicit time evolution of these variables is
[TABLE]
By equating all time derivatives to zero, we can identify the steady-state solutions for the mean-field description. From Eqs. (5,6,8) we find
[TABLE]
where
[TABLE]
and the convention . Obviously, the validity of these relations require that the series Eq. (12) defining is convergent. This is certainly not the case for , as always. As we will see, another case in which the series might diverge occurs when , when d’Alembert’s criterion does not ensure convergence as . Using Eqs. (7,9,10) in the steady state, one obtains easily , or
[TABLE]
The solutions to this equation provide the possible steady-state values of and those values of determine the other quantities, through Eqs. (11). Note that is always a trivial solution that corresponds to a symmetric steady state, with the same mean number of nodes with a given age having opposite states for . In any case, the steady-state solutions describe situations where are decreasing functions of the age. Note also that are always steady state solutions (absorbing states) of the dynamics, and indeed satisfy Eq. (13) as and .
For the particular case , we obtain the function . Eq. (13) is satisfied for any and, in fact, it can be shown that Eq. (7) becomes and, hence, . This is the behavior of the non-aging voter model in the thermodynamic limit that was described before. If , the distribution of ages in this steady state follows a geometric distribution: and . For , see section V.
For the activation probabilities Eq. (1), one finds 333Simpler expressions are obtained if is an integer number, e.g. , if .
[TABLE]
where is the Gauss hypergeometric function. This function, and hence the series Eq. (12), is always convergent for . If , we already know that the series is divergent. If , the series diverges whenever . Although these divergences may be problematic in the case , one can still study Eq. (13) as a limit or, alternatively regularize the sum Eq. (12) with a cut-off and study the dependence with 444In some cases this limit is non-trivial and must be performed carefully. Specially when decays very fast as in the exponential example. A possibility is to consider as time dependent and a phenomenological reasoning based on imposing convergence of with suggests . This is similar to imposing only for lower than the simulation time .. We checked that the predicted behavior of the fixed points of Eq. (13) in all cases is correct. In Figure (1a) we plot in a case of aging and anti-aging . We find in both cases that provides the three mentioned fixed points, . We also hypothesize that the sign of determines the stability or the direction of movement of , , thus for aging we expect the symmetric solution to be unstable and the absorbing states stable and the other way around for anti-aging Ozaita (2018). The time-dependence and ordering dynamics depending on the form of is studied in detail in the next section.
V Ordering dynamics
It is clear from the rules of the process that the absorbing states and must be solutions of the dynamical equations in all cases, whatever the particular form of . We now discuss these solutions. For the sake of concreteness we consider only the absorbing state, but an equivalent analysis could be performed at as well. Consider, then, that the system is at at time with all internal ages set to [math]. As there are no nodes in the state, the herding mechanism implies that a selected node can never switch state to and will remain in the state, hence increasing its internal age in one unit according to the rules of the model. Therefore we have , , and thus Eq. (10) reduces to , or , where the superscript indicates that those values are valid at the absorbing state . The solution of Eqs. (6) with (implying ) and the above mentioned initial conditions can be readily found as
[TABLE]
which indicates that the distribution of internal ages in the population follows a Poisson distribution with . We now study the stability of this solution.
We linearize Eqs. (5,9) around the solution at the absorbing state:
[TABLE]
The solution of Eq. (16) with the initial condition , , is
[TABLE]
Imposing , and replacing in Eq. (17) we obtain
[TABLE]
with
[TABLE]
The integro-differential Eq. (20) has to be solved with the initial condition .
A simple case in which the linearized equations (19,20) can be solved explicitly is , constant. In this case we have and , which after a tedious but straightforward calculation 555A possibility is to use the Laplace transform in Eqs. (19,20). leads to and , the well known solution of the voter model in the thermodynamic limit , 666Note here that in apparent contradiction to the result of section IV, . But this is a product of the linearization in Eqs. (19,20)..
As it follows from Eqs. (19,20) that , the stability of requires to study the sign of which after a simple calculation, turns out to be . Therefore, at this level, we obtain that in the anti-aging case, , the absorbing solution is unstable, while in the aging case, , the absorbing solution is stable, as suggested at the end of the previous section. Since for both functional forms Eqs. (1,2) in the anti-aging case there are no other solutions than , we predict that the final state of the dynamics is that of coexistence of states. This is indeed observed in the numerical simulations, as shown in Fig. (1b).
To study in full detail the return to the absorbing state of the variable we would need to solve the linearized equations (19,20). Although it does not seem to be possible to obtain their full time dependence analytically 777It is possible to try a power-series expansion and and obtain recurrence relations for the coefficients , but we have not been able to sum analytically the resulting series., simple arguments allow us to obtain the asymptotic behaviors for large . All we need for this purpose is the asymptotic expressions of and . We split the discussion in the cases and , as the difference is of crucial importance.
V.1 Aging case, Eq. (1) with .
This case includes the one studied in Fernández-Gracia et al. (2011) where the activation probability decayed as . With our explicit expression, Eq. (1), the calculation of Appendix A leads to an asymptotic behavior valid for ,
[TABLE]
Assuming power-law dependences as , and that the integral in Eq. (20) can be approximated as , with , which is justified if the time evolution of is slower than that of , we obtain consistently . Finally, using this value with Eq. (19) we get . This proves that the density of nodes in the state and, consequently, the density of active links go to zero as , explaining the numerical results of Ref. Fernández-Gracia et al. (2011).
To check the previous predictions, we plot in Fig. (2) the results coming from numerical simulations for and for different cases of and . As it is apparent, the figures confirm the validity of the asymptotic expansion and , as derived theoretically. Besides the asymptotic slope, we also plot in the figures the values of and obtained from the numerical solution of Eqs. (19,20) using a suitable integration algorithm described in Day (1967).
This discussion assumes that . If , it is and a direct integration of Eqs. (19,20) leads to an initial decay , followed by an intermediate regime . In the numerical simulations we observe that this is followed by a crossover to the same asymptotic power-law behavior . The reason why the integro-differential Eqs. (19,20) do not capture this crossover is because they are linearized. The solution of the full set of equations Eqs. (5-7) can be written with the different orders of the initial condition . Linearizing, we were able to find the first term . For this first term is a power law and it prevails over the rest for all times. For it decays exponentially and, apparently, the second order term is a power law. Thus for sufficiently large times it can happen that .
We now consider for . The kernel inside the integral of Eq. (18) is a Poisson distribution with , which localizes the integral of in a window of position and width determined by . This tells us that in order to find the qualitative shape of the age distribution for , we can replace the early dynamics approximate solution inside the integral at Eq. (18) by and by its asymptotic expansion form for :
[TABLE]
This is essentially a function localized around whose amplitude decreases as . This feature of the distribution of ages is analyzed in Fig. (3).
For , however, the major contribution to the integral at Eq. (18) comes from the region and we can use the asymptotic expression as a constant inside the integral
[TABLE]
whose validity is checked in Fig. (3).
V.2 Aging case, Eq. (1) with
In this case the activation probability starts at and decays to a non-null value . This is qualitatively similar to the case studied in Stark et al. (2008), although in that reference the final value was reached after a finite number of aging steps, but we do not expect this to be a significant difference.
As shown in Appendix A, tends to a non-null asymptotic value while decays exponentially:
[TABLE]
If we replace by this constant value in Eq. (20) we can use the Laplace transform of to find the explicit solution:
[TABLE]
and from Eq. (19) we get the Laplace transform of :
[TABLE]
Note that, due to the asymptotic behavior of , Eq. (27), the Laplace transform does not converge for . It turns out that has a zero located at . Hence, inverting the Laplace transform we find the asymptotic behavior . The exact value of has to be found solving numerically and it is found to depend non-trivially on the parameters of the model, see Figure (4a). As it is shown in this figure, we find for (aging) and for (anti-aging), thus the asymptotic behavior is in accordance to the linear stability analysis performed before. There are two cases where : (i) for which corresponds to the non-aging voter model with , and (ii) for where the time dependence changes from being exponential to a power law as explained in the previous section. The case that concern us in this section is , where we have for small and as approaches , increases until . We have checked the proposed asymptotic behavior in the particular case , , in Figure 5. The theory predicts an exponential decay with which matches perfectly compared to numerical simulation and also compared to the numerical solution of Eqs. (19,20). In the same figure we also show the special problematic case , with the same other parameters. In this case it is thus according to Eqs. (19,20) the solution is initially an exponential function as followed by the regime . This behavior is correct in the early states of the dynamics but after it there is a crossover to a slower exponential decay . We expect this effect to be also a result of linearization, where in , is slower than which is the result that we obtained. Thus for a large enough time it can happen that .
V.3 Aging case, Eq. (2) with
From expression Eq. (2) in Appendix A we find the asymptotic behavior for large
[TABLE]
where is the q-Pochhammer symbol qpo (2019). Note the remarkable property that which is something that we did not observe in the other cases. As explained in Appendix B, we can solve Eqs. (19,20) in this case replacing Eqs. (31,32). Approximately, we have the exponential decay and decays exponentially until it reaches a constant value that depends on the initial condition . Here is a constant that depends on and in an intricate way, see Appendix B. In summary, the system does not order in this case. We can deduce that this will be the case as long as or, equivalently, if converges. Obviously for this is never the case, while for if decays faster than the system does no order, but it does order if decays slower than . The case is somehow a critical value, where the system orders very slowly as a power law as explored in section V.1. We checked this prediction against the results of numerical simulations in Figure 6 for , . For there is again the failure of the linearization, where the theory predicts that the system orders exponentially, while the numerical simulation shows the same behavior initially until a crossover is reached to a constant value , which is smaller than in the case . According to this argument , and . The fixed point value is expected to scale as as the linear theory predicts . This and the linear prediction is checked with numerical simulation in Figure 7.
V.4 Aging case, Eq. (2) with
In this case, from expression Eq. (2) in Appendix A we find the asymptotic expressions
[TABLE]
Thus, this case is similar to the one studied in section V.2 where decays exponentially and goes to a constant. Following the procedure of that section we use the Laplace transform in Eq. (19,20) after replacing by its limit , see Appendix A for an explicit expression for these functions. In Figure (4b) we explored the root of the denominator of Eq. (29) as a function of . We find the same phenomenology, with for and for , for corresponding to the non-aging voter model. For we also obtain , in accordance to the result of the previous section where goes to a constant. In Figure 8 we compare the results of numerical simulations of with the numerical integration of Eqs. (19,20) and the asymptotic , for , , and . The theory predicts which is in good agreement to what it is observed in the simulation. The case is again singular as the linear approximation fails to reproduce the long term dynamics. In the linear theory we have which leads to . This exponential decay happens to be faster than the second order contribution, which is also exponential , and thus for long enough time the second order prevails.
VI Summary and conclusions
In this work we have studied memory effects (aging and anti-aging) in an opinion dynamics model, by means of analytical and numerical approaches. More specifically, we have focused on a mean-field (all-to-all) description of the voter model, where the influence of the aging of the individuals has been modelled through the activation probability, , or the copying probability of the individuals as a function of their internal age . Overall, we have shown that the functional form of has a crucial impact on how and whether the system reaches consensus or not. For the sake of concreteness, we have employed in our analysis monotonic functions with three parameters: (or for zero age), (or for age), and a third one, or , tuning the transition .
When the activation probability is constant , we recover the voter model. In this case, the system keeps its initial magnetization, provided the number of agents is large enough. If the activation probability is not constant and tends to a nonzero value for large ages, , the aging and anti-aging cases produce opposite effects in the dynamics. On the one hand, for the case of aging, with , the system orders with an exponential decay of the global opinion towards one of the consensus states . On the other hand, for the anti-aging case, with , the system reaches coexistence with .
As explained in Ref. Stark et al. (2008), the results for can be understood heuristically: the aging mechanism amplifies any small asymmetry in the initial conditions. Take an initial condition close to consensus, say , and all the agents with zero age. Then, the systems evolves in a first, very short stage as if no aging was present since all agents have the same herding rate. After the first transitions, most of the agents become older in the state , while a small fraction are younger with states . From a dynamical point of view, the only difference between the latter situation and the voter model is in the presence of the small fraction of the young agents. In the case of aging, since the rate of changing state is bigger for the young agents, the most probable event is that a young agent copies the state of an old one. Altogether, the magnetization decreases. On the contrary, for the case of anti-aging with the young agents having smaller rates, the most probable event is that of an old copying a young. Now, the magnetization increases.
For the case of aging with zero , , the system may or may not order, depending on how “fast” is the decay of to zero as the age increases. From the theoretical analysis of the integro-differential equations close to the consensus state, together with the results of numerical simulations, we infer the following criterion: if decays faster than then the system gets trapped in a frozen state given by the initial conditions, while the system orders in other cases. The decay can be then regarded as a critical functional form where we proved that the system orders with a power-law time dependence.
The results of can also be understood using qualitative arguments. The evolution of the magnetization can be explained as for the case , the only difference now being the time needed for the system to reach its final value. For decaying slow enough, the system reaches consensus but on a time much larger than that for , since it takes more time to convince all agents, specially the older ones. For the other case of decaying faster, almost all agents become zealots Khalil et al. (2018) on a short time scale, and the system gets frozen on an active state, , in general.
Appendix A Asymptotic expansions
A.1
If we use the expression Eq. (1) for the activation probabilities, we obtain the explicit expressions for the functions Eqs. (21,22)
[TABLE]
where is Kummer’s confluent hypergeometric function. Using the expansion for large : we obtain in the asymptotic limit the expansions:
[TABLE]
The expressions for assume . For we find
[TABLE]
It is also possible to obtain the Laplace transform of as:
[TABLE]
A.2
In the case of activation probabilities given by Eq. (2) it is possible to obtain analytically:
[TABLE]
although we have not been able to find a closed expression for ,
[TABLE]
with the q-Pochhammer symbol. It is possible, however, to obtain the asymptotic expansion for as
[TABLE]
In the particular case it yields .
Finally, we mention the Laplace transform of as:
[TABLE]
For the numerical calculation of and in this case we use Eqs.(45,47) using a sufficiently large number of terms in the sums.
Appendix B Solution of Eqs. (19,20) using Eqs. (31,32)
If is a constant, we can differentiate Eq. (19) and obtain a system of differential equations:
[TABLE]
From the last equation we get which, replaced in Eq.(48) leads to a closed second order linear differential equation for . The change of variables makes the coefficients of the equation to be simple polynomials and a standard technique based of Frobenius power-series expansion leads to the general solution , , where and are constants imposed by satisfying the initial condition , and
[TABLE]
In the limit , the asymptotic behavior is
[TABLE]
leading to
[TABLE]
We do not reproduce the explicit expressions for the constants and as they are too long and not very illuminating. It suffices to say that in the case , used in subsection V.3 it is .
Acknowledgements
Partial financial support has been received from the Agencia Estatal de Investigacion (AEI, Spain) and Fondo Europeo de Desarrollo Regional under Project PACSS RTI2018-093732-B-C21 (AEI/FEDER,UE) and the Spanish State Research Agency, through the Maria de Maeztu Program for units of Excellence in R&D (MDM-2017-0711). A. F. P. acknowledges support by the Formacion de Profesorado Universitario (FPU14/00554) program of Ministerio de Educacion, Cultura y Deportes (MECD) (Spain). We thank J. F. Gracia for valuable discussions concerning the results of Ref. Fernández-Gracia et al. (2011).
References
- Pastor-Satorras and Vespignani (2001)
R. Pastor-Satorras, A. Vespignani, Epidemic spreading in scale-free networks, Physical review letters 86 (14) (2001) 3200.
- Pastor-Satorras et al. (2015)
R. Pastor-Satorras, C. Castellano, P. Van Mieghem, A. Vespignani, Epidemic processes in complex networks, Reviews of modern physics 87 (3) (2015) 925.
- Abrams and Strogatz (2003)
D. M. Abrams, S. H. Strogatz, Linguistics: Modelling the dynamics of language death, Nature 424 (6951) (2003) 900.
- Castelló et al. (2006)
X. Castelló, V. M. Eguíluz, M. San Miguel, Ordering dynamics with two non-excluding options: bilingualism in language competition, New Journal of Physics 8 (12) (2006) 308.
- Stauffer et al. (2007)
D. Stauffer, X. Castelló, V. M. Eguiluz, M. San Miguel, Microscopic Abrams–Strogatz model of language competition, Physica A: Statistical Mechanics and its Applications 374 (2) (2007) 835–842.
- Vazquez et al. (2010)
F. Vazquez, X. Castelló, M. S. Miguel, Agent based models of language competition: Macroscopic descriptions and order-disorder transition, J. Stat. Mech. Theory Exp 2010.
- Kirman (1993)
A. Kirman, Ants, rationality, and recruitment, Quart. J. Econ 108 (1993) 137.
- Alfarano et al. (2005)
S. Alfarano, T. Lux, F. Wagner, Time Estimation of agent-based models: The case of an asymmetric Herding model, Comput. Econ 26 (2005) 19.
- Alfarano et al. (2008)
S. Alfarano, T. Lux, F. Wagner, Time variation of higher moments in a financial market with heterogeneous agents: An analytical approach, Econ. Dyn. Control 32 (2008) 101.
- Alfarano and Milaković (2009)
S. Alfarano, M. Milaković, Network structure and N-dependence in agent-based herding models, Journal of Economic Dynamics and Control 33 (2009) 78.
- Gontis and Kononovicius (2017)
V. Gontis, A. Kononovicius, Spurious Memory in Non-Equilibrium Stochastic Models of Imitative Behavior, Entropy 19 (2017) 387.
- Carro et al. (2015)
A. Carro, R. Toral, M. San Miguel, Markets, Herding and Response to External Information, PloS one 10 (7) (2015) e0133287.
- Kononovicius and Ruseckas (2019)
A. Kononovicius, J. Ruseckas, Order book model with herd behavior exhibiting long-range memory, Physica A: Statistical Mechanics and its Applications 525 (2019) 171 – 191.
- Barrat et al. (2008)
A. Barrat, M. Barthelemy, A. Vespignani, Dynamical processes on complex networks, Cambridge university press, 2008.
- Castellano et al. (2009)
C. Castellano, S. Fortunato, V. Loreto, Statistical physics of social dynamics, Reviews of modern physics 81 (2) (2009) 591.
- Fernández-Gracia et al. (2014)
J. Fernández-Gracia, K. Suchecki, J. J. Ramasco, M. San Miguel, V. M. Eguíluz, Is the Voter Model a Model for Voters?, Phys. Rev. Lett. 112 (2014) 158701.
- Clifford and Sudbury (1973)
P. Clifford, A. Sudbury, A model for spatial conflict, Biometrika 60 (3) (1973) 581–588.
- Holley and Liggett (1975)
R. A. Holley, T. M. Liggett, Ergodic theorems for weakly interacting infinite systems and the voter model, The annals of probability 3 (4) (1975) 643–663.
- Vazquez and Eguíluz (2008)
F. Vazquez, V. M. Eguíluz, Analytical solution of the voter model on uncorrelated networks, New Journal of Physics 10 (6) (2008) 063011.
- Redner (2018)
S. Redner, Reality Inspired Voter Models: A Mini-Review, arXiv e-prints (2018) arXiv:1811.11888.
- Krapivsky et al. (2010)
P. L. Krapivsky, S. Redner, E. Ben-Naim, A kinetic view of statistical physics, Cambridge University Press, 2010.
- Liggett (2012)
T. M. Liggett, Interacting particle systems, vol. 276, Springer Science & Business Media, 2012.
- Peralta et al. (2018a)
A. F. Peralta, A. Carro, M. S. Miguel, R. Toral, Stochastic pair approximation treatment of the noisy voter model, New Journal of Physics 20 (2018a) 103045.
- Carro et al. (2016)
A. Carro, R. Toral, M. S. Miguel, The noisy voter model on complex networks, Scientific Reports 6 (2016) 24775.
- Nyczka et al. (2012)
P. Nyczka, K. Sznajd-Weron, J. Cisło, Phase transitions in the -voter model with two types of stochastic driving, Phys. Rev. E 86 (2012) 011105.
- Nyczka and Sznajd-Weron (2013)
P. Nyczka, K. Sznajd-Weron, Anticonformity or Independence?—Insights from Statistical Physics, Journal of Statistical Physics 151 (1) (2013) 174–202.
- J\polhkedrzejewski (2017)
A. J\polhkedrzejewski, Pair approximation for the -voter model with independence on complex networks, Phys. Rev. E 95 (2017) 012307.
- Peralta et al. (2018b)
A. F. Peralta, A. Carro, M. S. Miguel, R. Toral, Analytical and numerical study of the non-linear noisy voter on complex networks, Chaos 28 (2018b) 075516.
- Vieira and Anteneodo (2018)
A. R. Vieira, C. Anteneodo, Threshold -voter model, Phys. Rev. E 97 (2018) 052106.
- Raducha et al. (2018)
T. Raducha, B. Min, M. S. Miguel, Coevolving nonlinear voter model with triadic closure, EPL (Europhysics Letters) 124 (3) (2018) 30001.
- Min and Miguel (2019)
B. Min, M. S. Miguel, Multilayer coevolution dynamics of the nonlinear voter model, New Journal of Physics 21 (3) (2019) 035004.
- Herrerías-Azcué and Galla (2019)
F. Herrerías-Azcué, T. Galla, Consensus and diversity in multi-state noisy voter models, arXiv e-prints (2019) arXiv:1903.09198.
- Vazquez et al. (2019)
F. Vazquez, E. S. Loscar, G. Baglietto, A multi-state voter model with imperfect copying, arXiv e-prints (2019) arXiv:1902.07253.
- Khalil et al. (2018)
N. Khalil, M. S. Miguel, R. Toral, Zealots in the mean-field noisy voter model, Physical Review 97.
- Khalil and Toral (2019)
N. Khalil, R. Toral, The noisy voter model under the influence of contrarians, Physica A 515 (2019) 81–92.
- J\polhkedrzejewski and Sznajd-Weron (2018)
A. J\polhkedrzejewski, K. Sznajd-Weron, Impact of memory on opinion dynamics, Physica A: Statistical Mechanics and its Applications 505 (2018) 306 – 315.
- Min et al. (2013)
B. Min, K. i. Goh, I. m. Kim, Suppression of epidemic outbreaks with heavy-tailed contact dynamics, Europhys. Lett 103 (2013) 50002.
- Boguñá et al. (2014)
M. Boguñá, L. Lafuerza, R. Toral, M. Á. Serrano, Simulating non-Markovian stochastic processes, Physical Review E 90 (2014) 042108.
- Starnini et al. (2017)
M. Starnini, J. P. Gleeson, M. Boguñá, Equivalence between Non-Markovian and Markovian Dynamics in Epidemic Spreading Processes, Phys. Rev. Lett 118 (2017) 128301.
- Blythe and Anderson (1988)
S. P. Blythe, R. M. Anderson, Variable Infectiousness in HIV Transmission Models, Math. Med. Biol 5 (1988) 181.
- van Kampen (1998)
N. G. van Kampen, Remarks on Non-Markov Processes, Braz. J. Phys 28 (1998) 90.
- Łucza (2005)
J. Łucza, Non-markovian stochastic processes: Colored noise, Chaos 15 (2005) 026107.
- Fernández-Gracia et al. (2011)
J. Fernández-Gracia, V. M. Eguíluz, M. San Miguel, Update rules and interevent time distributions: Slow ordering versus no ordering in the voter model, Physical Review E 84 (1) (2011) 015103.
- Stark et al. (2008)
H.-U. Stark, C. J. Tessone, F. Schweitzer, Decelerating Microdynamics Can Accelerate Macrodynamics in the Voter Model, Phys. Rev. Lett 101 (2008) 018701.
- Ben-Naim et al. (1996)
E. Ben-Naim, L. Frachebourg, P. L. Krapivsky, Coarsening and persistence in the voter model, Physical Review E 53 (4) (1996) 3078.
- Sood and Redner (2005)
V. Sood, S. Redner, Voter model on heterogeneous graphs, Physical review letters 94 (17) (2005) 178701.
- Suchecki et al. (2005)
K. Suchecki, V. M. Eguíluz, M. San Miguel, Voter model dynamics in complex networks: Role of dimensionality, disorder, and degree distribution, Physical Review E 72 (3) (2005) 036132.
- Martins and Galam (2013)
A. C. R. Martins, S. Galam, Building up of individual inflexibility in opinion dynamics, Physical Review E 87 (2013) 042807.
- Artime et al. (2018)
O. Artime, A. F. Peralta, R. Toral, J. J. Ramasco, M. S. Miguel, Aging-induced phase transition, Phys. Rev. E 98 (2018) 032104.
- Artime et al. (2018)
O. Artime, A. Carro, A. F. Peralta, J. J. Ramasco, M. San Miguel, R. Toral, Herding and idiosyncratic choices: Nonlinearity and aging-induced transitions in the noisy voter model, arXiv e-prints arXiv:1812.05378.
- Considine et al. (1989)
D. Considine, S. Redner, H. Takayasu, Comment on “Noise-induced bistability in a Monte Carlo surface-reaction model”, Phys. Rev. Lett. 63 (1989) 2857–2857.
- Granovsky and Madras (1995)
B. L. Granovsky, N. Madras, The noisy voter model, Stoch. Proc. Appl. 55 (1995) 23.
- Wu et al. (2010)
Y. Wu, C. Zhou, J. Xiao, J. Kurths, H. J. Schellnhuber, Evidence for a bimodal distribution in human communication, PNAS 107 (2010) 18803.
- Candia et al. (2008)
J. Candia, M. C. González, P. Wang, T. Schoenharl, G. Madey, A. Barabási, Uncovering individual and collective human dynamics from mobile phone records, J. Phys. A: Math. Theor 41 (2008) 224015.
- Pérez et al. (2016)
T. Pérez, K. Klemm, V. M. Eguíluz, Competition in the presence of aging: dominance, coexistence, and alternation between states, Sci. Rep 6 (2016) 21128.
- Artime et al. (2017a)
O. Artime, J. J. Ramasco, M. S. Miguel, Dynamics on networks: competition of temporal and topological correlations, Sci. Rep 7 (2017a) 41627.
- Artime et al. (2017b)
O. Artime, J. F. Gracia, J. J. Ramasco, M. S. Miguel, Joint effect of ageing and multilayer structure prevents ordering in the voter model, Sci. Rep 7 (2017b) 7166.
- Escaff et al. (2018)
D. Escaff, R. Toral, C. Van Den Broeck, K. Lindenberg, A continuous-time persistent random walk model for flocking, Chaos 28 (7) (2018) 075507.
- Artime et al. (2018)
O. Artime, N. Khalil, R. Toral, M. S. Miguel, First-passage distributions for the one-dimensional Fokker-Planck equation, Physical Review 98.
- van Kampen (2007)
N. van Kampen, Stochastic Processes in Physics and Chemistry, North-Holland, Amsterdam, third edn., 2007.
- Toral and Colet (2014)
R. Toral, P. Colet, Stochastic Numerical Methods: An Introduction for Students and Scientists, Wiley, URL http://eu.wiley.com/WileyCDA/WileyTitle/productCd-3527411496.html, 2014.
- Peralta and Toral (2018)
A. F. Peralta, R. Toral, System-size expansion of the moments of a master equation, Chaos 28 (2018) 106303.
- Ozaita (2018)
J. Ozaita, Noisy voter model with partial aging and anti-aging, Master’s thesis, University of the Balearic Islands, Palma, Spain, 2018.
- Day (1967)
J. T. Day, Note on the Numerical Solution of Integro-Differential Equations, The Computer Journal 9 (4) (1967) 394–395.
- qpo (2019)
URL https://en.wikipedia.org/wiki/Q-Pochhammer_symbol, 2019.
- Note1 (????)
A similar qualitative behavior is observed for the magnetization. In fact, for the all-to-all connectivity considered later in this paper, one has simply , although in other networks must be considered as an independent variable.
- Note2 (????)
Strictly speaking, the aforementioned power-law dependence was reported in Fernández-Gracia et al. (2011) for the cumulative inter-event time distribution . One of the authors of that reference (J. F. Gracia, private communication) has confirmed to us that the aforementioned asymptotic time dependence of was also observed, but not displayed.
- Note3 (????)
Simpler expressions are obtained if is an integer number, e.g. , if .
- Note4 (????)
In some cases this limit is non-trivial and must be performed carefully. Specially when decays very fast as in the exponential example. A possibility is to consider as time dependent and a phenomenological reasoning based on imposing convergence of with suggests . This is similar to imposing only for lower than the simulation time .
- Note5 (????)
A possibility is to use the Laplace transform in Eqs. (19,20).
- Note6 (????)
Note here that in apparent contradiction to the result of section IV, . But this is a product of the linearization in Eqs. (19,20).
- Note7 (????)
It is possible to try a power-series expansion and and obtain recurrence relations for the coefficients , but we have not been able to sum analytically the resulting series.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Pastor-Satorras and Vespignani (2001) R. Pastor-Satorras, A. Vespignani, Epidemic spreading in scale-free networks, Physical review letters 86 (14) (2001) 3200.
- 2Pastor-Satorras et al. (2015) R. Pastor-Satorras, C. Castellano, P. Van Mieghem, A. Vespignani, Epidemic processes in complex networks, Reviews of modern physics 87 (3) (2015) 925.
- 3Abrams and Strogatz (2003) D. M. Abrams, S. H. Strogatz, Linguistics: Modelling the dynamics of language death, Nature 424 (6951) (2003) 900.
- 4Castelló et al. (2006) X. Castelló, V. M. Eguíluz, M. San Miguel, Ordering dynamics with two non-excluding options: bilingualism in language competition, New Journal of Physics 8 (12) (2006) 308.
- 5Stauffer et al. (2007) D. Stauffer, X. Castelló, V. M. Eguiluz, M. San Miguel, Microscopic Abrams–Strogatz model of language competition, Physica A: Statistical Mechanics and its Applications 374 (2) (2007) 835–842.
- 6Vazquez et al. (2010) F. Vazquez, X. Castelló, M. S. Miguel, Agent based models of language competition: Macroscopic descriptions and order-disorder transition, J. Stat. Mech. Theory Exp 2010.
- 7Kirman (1993) A. Kirman, Ants, rationality, and recruitment, Quart. J. Econ 108 (1993) 137.
- 8Alfarano et al. (2005) S. Alfarano, T. Lux, F. Wagner, Time Estimation of agent-based models: The case of an asymmetric Herding model, Comput. Econ 26 (2005) 19.
