Reply to "Comment on `Performance of different synchronization measures in real data: A case study on electroencephalographic signals'"

We agree with the Comment by Nicolaou and Nasuto about the utility of mutual information (MI) when properly estimated and we also concur with their view that the estimation based on k nearest neighbors gives optimal results. However, we claim that embedding parameters can indeed change MI results, as we show for the electroencephalogram data sets of our original study and for coupled chaotic systems. Furthermore, we show that proper embedding can actually improve the estimation of MI with the k nearest neighbors algorithm

Epidemic analysis of the second-order transition in the Ziff-Gulari-Barshad surface-reaction model

We study the dynamic behavior of the Ziff-Gulari-Barshad (ZGB) irreversible surface-reaction model around its kinetic second-order phase transition, using both epidemic and poisoning-time analyses. We find that the critical point is given by p_1 = 0.3873682 \pm 0.0000015, which is lower than the previous value. We also obtain precise values of the dynamical critical exponents z, \delta, and \eta which provide further numerical evidence that this transition is in the same universality class as directed percolation.Comment: REVTEX, 4 pages, 5 figures, Submitted to Physical Review

Percolation in Media with Columnar Disorder

We study a generalization of site percolation on a simple cubic lattice, where not only single sites are removed randomly, but also entire parallel columns of sites. We show that typical clusters near the percolation transition are very anisotropic, with different scaling exponents for the sizes parallel and perpendicular to the columns. Below the critical point there is a Griffiths phase where cluster size distributions and spanning probabilities in the direction parallel to the columns have power law tails with continuously varying non-universal powers. This region is very similar to the Griffiths phase in subcritical directed percolation with frozen disorder in the preferred direction, and the proof follows essentially the same arguments as in that case. But in contrast to directed percolation in disordered media, the number of active ("growth") sites in a growing cluster at criticality shows a power law, while the probability of a cluster to continue to grow shows logarithmic behavior.Comment: 9 pages, 9 figure

Event synchronization: a simple and fast method to measure synchronicity and time delay patterns

We propose a simple method to measure synchronization and time delay patterns between signals. It is based on the relative timings of events in the time series, defined e.g. as local maxima. The degree of synchronization is obtained from the number of quasi-simultaneous appearances of events, and the delay is calculated from the precedence of events in one signal with respect to the other. Moreover, we can easily visualize the time evolution of the delay and synchronization level with an excellent resolution. We apply the algorithm to short rat EEG signals, some of them containing spikes. We also apply it to an intracranial human EEG recording containing an epileptic seizure, and we propose that the method might be useful for the detection of foci and for seizure prediction. It can be easily extended to other types of data and it is very simple and fast, thus being suitable for on-line implementations.Comment: 6 pages, including 6 figures, RevTe

Compact parity conserving percolation in one-dimension

Compact directed percolation is known to appear at the endpoint of the directed percolation critical line of the Domany-Kinzel cellular automaton in 1+1 dimension. Equivalently, such transition occurs at zero temperature in a magnetic field H, upon changing the sign of H, in the one-dimensional Glauber-Ising model with well known exponents characterising spin-cluster growth. We have investigated here numerically these exponents in the non-equilibrium generalization (NEKIM) of the Glauber model in the vicinity of the parity-conserving phase transition point of the kinks. Critical fluctuations on the level of kinks are found to affect drastically the characteristic exponents of spreading of spins while the hyperscaling relation holds in its form appropriate for compact clusters.Comment: 7 pages, 7 figures embedded in the latex, final form before J.Phys.A publicatio

Sandpile model on a quenched substrate generated by kinetic self-avoiding trails

Kinetic self-avoiding trails are introduced and used to generate a substrate of randomly quenched flow vectors. Sandpile model is studied on such a substrate with asymmetric toppling matrices where the precise balance between the net outflow of grains from a toppling site and the total inflow of grains to the same site when all its neighbors topple once is maintained at all sites. Within numerical accuracy this model behaves in the same way as the multiscaling BTW model.Comment: Four pages, five figure

Phase transitions and critical behaviour in one-dimensional non-equilibrium kinetic Ising models with branching annihilating random walk of kinks

One-dimensional non-equilibrium kinetic Ising models evolving under the competing effect of spin flips at zero temperature and nearest-neighbour spin exchanges exhibiting directed percolation-like parity conserving(PC) phase transition on the level of kinks are now further investigated, numerically, from the point of view of the underlying spin system. Critical exponents characterising its statics and dynamics are reported. It is found that the influence of the PC transition on the critical exponents of the spins is strong and the origin of drastic changes as compared to the Glauber-Ising case can be traced back to the hyperscaling law stemming from directed percolation(DP). Effect of an external magnetic field, leading to DP-type critical behaviour on the level of kinks, is also studied, mainly through the generalised mean field approximation.Comment: 15 pages, using RevTeX, 13 Postscript figures included, submitted to J.Phys.A, figures 12 and 13 fixe

The three species monomer-monomer model in the reaction-controlled limit

We study the one dimensional three species monomer-monomer reaction model in the reaction controlled limit using mean-field theory and dynamic Monte Carlo simulations. The phase diagram consists of a reactive steady state bordered by three equivalent adsorbing phases where the surface is saturated with one monomer species. The transitions from the reactive phase are all continuous, while the transitions between adsorbing phases are first-order. Bicritical points occur where the reactive phase simultaneously meets two adsorbing phases. The transitions from the reactive to an adsorbing phase show directed percolation critical behaviour, while the universal behaviour at the bicritical points is in the even branching annihilating random walk class. The results are contrasted and compared to previous results for the adsorption-controlled limit of the same model.Comment: 12 pages using RevTeX, plus 4 postscript figures. Uses psfig.sty. accepted to Journal of Physics

One Dimensional Nonequilibrium Kinetic Ising Models with Branching Annihilating Random Walk

Nonequilibrium kinetic Ising models evolving under the competing effect of spin flips at zero temperature and nearest neighbour spin exchanges at $T=\infty$ are investigated numerically from the point of view of a phase transition. Branching annihilating random walk of the ferromagnetic domain boundaries determines the steady state of the system for a range of parameters of the model. Critical exponents obtained by simulation are found to agree, within error, with those in Grassberger's cellular automata.Comment: 10 pages, Latex, figures upon request, SZFKI 05/9