## Brazilian Journal of Physics

##
*Print version* ISSN 0103-9733*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.30 no.4 São Paulo Dec. 2000

#### http://dx.doi.org/10.1590/S0103-97332000000400007

**Damage spreading in the one-dimensional Hinrichsen-Domany model **

Everaldo Arashiro and J. R. Drugowich de Felício

*Departamento de Física e Matemática*

*Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto*

*Universidade de São Paulo*

*Av. dos Bandeirantes, 3900 - 14040-901 - Ribeirão Preto, SP, Brasil*

**Received on 1 November, 2000**

In this paper we study damage spreading in a one-dimensional model under two dynamics introduced by Hinrichsen and Domany. In particular, we study the effects of synchronous and asynchronous updating on the spreading properties. We show that the damage does not spread when the second dynamic is implemented in a synchronous way. We find that the rules for updating the damage produced by this dynamics, as the temperature goes to infinity and a certain parameter l is zero, are equivalent to those of Grassberger's well-known model A cellular automaton.

**I Introduction **

The damage spreading technique was proposed by Kauffman thirty years ago in the context of biologically motivated cellular automata [1]. More recently, this concept was extended to analyse kinetic Ising and Potts models which are described by well-known Hamiltonians and exhibit critical properties [2]. The technique requires following the behavior of two similar samples (slightly different initial conditions) under same thermal noise. Usually, the evolution is done by means of Glauber, Metropolis or heat-bath updating. However, in order to understand the role of the dynamics, other kinds of updating have been studied. In fact, this was the goal of a recent paper by Hinrichsen and Domany which describes the spreading properties in a one-dimensional model submitted to two kinds of updating (henceforth named I and II). The first case (dynamics I) is very similar to Glauber dynamics but depends on the sign of the two neighbours instead of the site itself. The new value of the variable at the site *i* is:

where s_{i}( *t*) = ±1, *z* is a random number (between zero and 1) and the probability *q*_{i}(*t*) is given by

with

and *J* the exchange coupling constant (*H* = -*J * s_{i}s_{j}). Dynamics II combines the Glauber updating with information about three sites (two first neighbours and the site itself). In this case, the new value for s is obtained by

if a random number *r* is greater than a given parameter l. If *r* < l, the update follows Glauber dynamics

The two rules above can be combined in one updating, by introducing another parameter, *y*. Thus,

where

z is a random number between 0 and 1. In both cases there is a damage spreading transition, but the critical properties are totally different. Whereas the transition in case I belongs to the directed percolation universality class, dynamics II leads to a "parity-conserving" phase transition. In the first case, the evolution of the damage (difference between two samples) can be described by Domany-Kinzel cellular automata rules, at least when *T* = ¥ [3]. Such a mapping is not known for the second case. In this paper we address two questions: first, What is the effect of synchronization on the spreading of the damage; second, Is there a known model that dynamics II can be mapped to?

**II The role of updating **

To implement the dynamics presented above, we have at least two possibilities: updating one site at a time (asynchronous or "continuous dynamics") as is usual in simulation of Ising-like models, or changing all the sites at the same time (parallel or synchronous updating), usually done in cellular automata. Our simulations for dynamics I show that damage spreading occurs in both cases. We checked this conclusion by simulating the Domany-Kinzel cellular automaton (related to dynamics I, at *T* = ¥, via a mapping found by Hinrichsen and Domany) with continuous dynamics. Evolution of the damage is the same in both cases (see Fig. 1). However, for dynamics II, the situation is completely different. Depending on the kind of the updating (synchronous or asynchronous) damage spreading may or may not occur (see Fig. 2). In their study Hinrichsen and Domany [3] used only parallel updating [4]. If they had used continuous dynamics they would not have observed damage spreading, as we have shown in our simulations. In the next section, we will obtain a mapping of dynamics II to the model A cellular automaton of Grassberger et al.[5]. This will allow a better understanding of our results on damage spreading.

**III Mapping to Grassberger's model A **

We begin by setting l = 0 (which implies working only with the equation (3.a) for updating the spins) and taking high temperature limit (*T *® ¥). In this case, equation (3.a) becomes

which means that the evolution of the damage ( D ) follows deterministic rules which depend on the site itself and the two nearest neighbours:

Now we observe that these rules correspond to the one-dimensional cellular automaton obeying rule 22 (in Wolfram's classification), which is known to exhibit chaotic behavior. On the other hand, we note that the rule 22 is exactly the limit of the Grassberger's model A below with *p* = 1. Thus the damage spreads.

In this sense, we can conclude that Hinrichsen-Domany's dynamics II exhibits damage spreading at least when l = 0 and *T* = ¥. On the other hand, for l = 1, and *T* = ¥, equation 4 can be simplified to

which implies another kind of updating for the damage. We find that the local damage in this case evolves exactly as the one-dimensional cellular automaton governed by the rule 204 (again in Wolfram's notation).

We remark that damage does not spread in this limit because the value of D is conserved in the evolution (0 ® 0 and 1 ® 1 ). The absence of the damage spreading in this case is not surprising, since for l = 1 we recover Glauber dynamics.

Figure 1. Patterns created with Domany-Kinzel automaton rules (with p_{1} = 0.83 e p_{2} = 0). (a) Synchronous updating and (b) Asynchronous updating. |

Figure 2. Patterns created with model A cellular automaton rules (with p = 0.3) from initial states containing a single kink. (a) Synchronous updating and (b) Asynchronous updating. |

**IV Conclusions**

We studied damage spreading in a one-dimensional model evolving under two different dynamics, proposed by Hinrichsen and Domany [3]. We showed that synchronous updating can destroy the damage spreading when the dynamics II is used. In addition, we obtained a mapping of dynamics II to the model A cellular automaton introduced by Grassberger et al.[5]. In this light, it is not strange that critical exponents obtained by HD were in the parity conserving universality class [6].

**References **

[1] S. Kauffman, J. Theor. Biol. **22**, 437 (1969). [ Links ]

[2] H.E. Stanley et al., Phys. Rev. Lett. **59**, 2326 (1987) [ Links ]

[3] H. Hinrichsen and E. Domany, Phys. Rev. E**56**, 94 (1997). [ Links ]

[4] H. Hinrichsen, private communication.

[5] P. Grassberger, F. Krause and T. von der Twer, J. Phys. A: Math. Gen. **17**, (1984). [ Links ]

[6] G. Ódor and N. Menyhárd, Phys. Rev. E **57**, 5168 (1998). [ Links ]