Q factor: A measure of competition between the topper and the average in percolation and in self-organized criticality.

We define the Q factor in the percolation problem as the quotient of the size of the largest cluster and the average size of all clusters. As the occupation probability p is increased, the Q factor for the system size L grows systematically to its maximum value Q_{max}(L) at a specific value p_{max}(L) and then gradually decays. Our numerical study of site percolation problems on the square, triangular, and simple cubic lattices exhibits that the asymptotic values of p_{max}, though close, are distinct from the corresponding percolation thresholds of these lattices. We also show, using scaling analysis, that at p_{max} the value of Q_{max}(L) diverges as L^{d} (d denoting the dimension of the lattice) as the system size approaches its asymptotic limit. We further extend this idea to nonequilibrium systems such as the sandpile model of self-organized criticality. Here the Q(ρ,L) factor is the quotient of the size of the largest avalanche and the cumulative average of the sizes of all the avalanches, with ρ the drop density of the driving mechanism. This study was prompted by some observations in sociophysics.

Non-Abelian sandpile automata with height restrictions.

We have studied the properties of a sandpile automata under the constraint of height restriction of sand columns. In this sandpile, an active site transfers a grain to a neighboring site if and only if the height of the sand column at the destination site is less than a preassigned value n_{c}. This sandpile was studied by Dickman et al. [Phys. Rev. E 66, 016111 (2002)1063-651X10.1103/PhysRevE.66.016111] in a conserved system with a fixed number of sand grains. In contrast, we have studied the avalanche dynamics of the driven sandpile under the open boundary conditions. The deterministic dynamics of the Bak, Tang, and Wiesenfeld (BTW) sandpile under the height restriction is found to be non-Abelian. Using numerical results, we argue that the steady states of the sandpile are exactly the recurrent states of the BTW sandpile, but occur with nonuniform probabilities. A detailed analysis of the cluster size distributions indicates that the associated exponent values are likely to be different from those of the BTW sandpile. The other differences include that the drop number distribution decays as a power law, and the largest avalanche size grows as the fourth power of the system size.

Nonstationary but quasisteady states in self-organized criticality.

The notion of self-organized criticality (SOC) was conceived to interpret the spontaneous emergence of long-range correlations in nature. Since then many different models have been introduced to study SOC. All of them have a few common features: externally driven dynamical systems self-organize themselves to nonequilibrium stationary states exhibiting fluctuations of all length scales as the signatures of criticality. In contrast, we have studied here in the framework of the sandpile model a system that has mass inflow but no outflow. There is no boundary, and particles cannot escape from the system by any means. Therefore, there is no current balance, and consequently it is not expected that the system would arrive at a stationary state. In spite of that, it is observed that the bulk of the system self-organizes to a quasisteady state where the grain density is maintained at a nearly constant value. Power law distributed fluctuations of all lengths and time scales have been observed, which are the signatures of criticality. Our detailed computer simulation study gives the set of critical exponents whose values are very close to their counterparts in the original sandpile model. This study indicates that (i) a physical boundary and (ii) the stationary state, though sufficient, may not be the necessary criteria for achieving SOC.

Bond percolation between k separated points on a square lattice.

We consider a percolation process in which k points separated by a distance proportional the system size L simultaneously connect together (k>1), or a single point at the center of a system connects to the boundary (k=1), through adjacent connected points of a single cluster. These processes yield new thresholds p[over ¯]_{ck} defined as the average value of p at which the desired connections first occur. These thresholds not sharp, as the distribution of values of p_{ck} for individual samples remains broad in the limit of L→∞. We study p[over ¯]_{ck} for bond percolation on the square lattice and find that p[over ¯]_{ck} are above the normal percolation threshold p_{c}=1/2 and represent specific supercritical states. The p[over ¯]_{ck} can be related to integrals over powers of the function P_{∞}(p) equal to the probability a point is connected to the infinite cluster; we find numerically from both direct simulations and from measurements of P_{∞}(p) on L×L systems that for L→∞, p[over ¯]_{c1}=0.51755(5), p[over ¯]_{c2}=0.53219(5), p[over ¯]_{c3}=0.54456(5), and p[over ¯]_{c4}=0.55527(5). The percolation thresholds p[over ¯]_{ck} remain the same, even when the k points are randomly selected within the lattice. We show that the finite-size corrections scale as L^{-1/ν_{k}} where ν_{k}=ν/(kß+1), with ß=5/36 and ν=4/3 being the ordinary percolation critical exponents, so that ν_{1}=48/41, ν_{2}=24/23, ν_{3}=16/17, ν_{4}=6/7, etc. We also study three-point correlations in the system and show how for p>p_{c}, the correlation ratio goes to 1 (no net correlation) as L→∞, while at p_{c} it reaches the known value of 1.022.

Brittle to quasibrittle transition in a compound fiber bundle.

The brittle to quasibrittle transition has been studied for a compound of two different kinds of fibrous materials, having distinct difference in their breaking strengths under the framework of the fiber bundle model. A random fiber bundle model has been devised with a bimodal distribution of the breaking strengths of the individual fibers. The bimodal distribution is assumed to consist of two symmetrically placed rectangular probability distributions of strengths p and 1-p, each of width d, and separated by a gap 2s. Different properties of the transition have been studied varying these three parameters and using the well-known equal load-sharing dynamics. Our study exhibits a brittle to quasibrittle transition at the critical width d_{c}(s,p)=p(1/2-s)/(1+p) confirmed by our numerical results.

Colored percolation.

A model called "colored percolation" has been introduced with its infinite number of versions in two dimensions. The sites of a regular lattice are randomly occupied with probability p and are then colored by one of the n distinct colors using uniform probability q=1/n. Denoting different colors by the letters of the Roman alphabet, we have studied different versions of the model like AB,ABC,ABCD,ABCDE,... etc. Here, only those lattice bonds having two different colored atoms at the ends are defined as connected. The percolation threshold p_{c}(n) asymptotically converges to its limiting value of p_{c} as 1/n. The model has been generalized by introducing a preference towards a subset of colors when m out of n colors are selected with probability q/m each and the rest of the colors are selected with probability (1-q)/(n-m). It has been observed that p_{c}(q,m) depends nontrivially on q and has a minimum at q_{min}=m/n. In another generalization the fractions of bonds between similarly and dissimilarly colored atoms have been treated as independent parameters. Phase diagrams in this parameter space have been drawn exhibiting percolating and nonpercolating phases.

Double transition in a model of oscillating percolation.

Two distinct transition points have been observed in a problem of lattice percolation studied using a system of pulsating disks. Sites on a regular lattice are occupied by circular disks whose radii vary sinusoidally within [0,R_{0}] starting from a random distribution of phase angles. A lattice bond is said to be connected when its two end disks overlap with each other. Depending on the difference of the phase angles of these disks, a bond may be termed as dead or live. While a dead bond can never be connected, a live bond is connected at least once in a complete time period. Two different time scales can be associated with such a system, leading to two transition points. Namely, a percolation transition occurs at R_{0c}=0.908(5) when a spanning cluster of connected bonds emerges in the system. Here, information propagates across the system instantly, i.e., with infinite speed. Secondly, there exists another transition point R_{0}^{*}=0.5907(3) where the giant cluster of live bonds spans the lattice. In this case the information takes finite time to propagate across the system through the dynamical evolution of finite-size clusters. This passage time diverges as R_{0}→R_{0}^{*} from above. Both the transitions exhibit the critical behavior of ordinary percolation transition. The entire scenario is robust with respect to the distribution of frequencies of the individual disks. This study may be relevant in the context of wireless sensor networks.

Brittle-to-quasibrittle transition in bundles of nonlinear elastic fibers.

Properties of the fiber bundle model have been studied using equal load-sharing dynamics where each fiber obeys a nonlinear stress (s)-strain (x) characteristic function s=G(x) till its breaking threshold. In particular, four different functional forms have been studied: G(x)=e^{αx}, 1+x^{α}, x^{α}, and xe^{αx} where α is a continuously tunable parameter of the model in all cases. Analytical studies, supported by extensive numerical calculations of this model, exhibit a brittle to quasibrittle phase transition at a critical value of α_{c} only in the first two cases. This transition is characterized by the weak power law modulated logarithmic (brittle) and logarithmic (quasibrittle) dependence of the relaxation time on the two sides of the critical point. Moreover, the critical load σ_{c}(α) for the global failure of the bundle depends explicitly on α in all cases. In addition, four more cases have also been studied, where either the nonlinear functional form or the probability distribution of breaking thresholds has been suitably modified. In all these cases similar brittle to quasibrittle transitions have been observed.

Percolation model with an additional source of disorder.

The ranges of transmission of the mobiles in a mobile ad hoc network are not uniform in reality. They are affected by the temperature fluctuation in air, obstruction due to the solid objects, even the humidity difference in the environment, etc. How the varying range of transmission of the individual active elements affects the global connectivity in the network may be an important practical question to ask. Here a model of percolation phenomena, with an additional source of disorder, is introduced for a theoretical understanding of this problem. As in ordinary percolation, sites of a square lattice are occupied randomly with probability p. Each occupied site is then assigned a circular disk of random value R for its radius. A bond is defined to be occupied if and only if the radii R_{1} and R_{2} of the disks centered at the ends satisfy a certain predefined condition. In a very general formulation, one divides the R_{1}-R_{2} plane into two regions by an arbitrary closed curve. One defines a point within one region as representing an occupied bond; otherwise it is a vacant bond. The study of three different rules under this general formulation indicates that the percolation threshold always varies continuously. This threshold has two limiting values, one is p_{c}(sq), the percolation threshold for the ordinary site percolation on the square lattice, and the other is unity. The approach of the percolation threshold to its limiting values are characterized by two exponents. In a special case, all lattice sites are occupied by disks of random radii R∈{0,R_{0}} and a percolation transition is observed with R_{0} as the control variable, similar to the site occupation probability.

Fiber bundle model with highly disordered breaking thresholds.

We present a study of the fiber bundle model using equal load-sharing dynamics where the breaking thresholds of the fibers are drawn randomly from a power-law distribution of the form p(b)∼b-1 in the range 10-ß to 10ß. Tuning the value of ß continuously over a wide range, the critical behavior of the fiber bundle has been studied both analytically as well as numerically. Our results are: (i) The critical load σc(ß,N) for the bundle of size N approaches its asymptotic value σc(ß) as σc(ß,N)=σc(ß)+AN-1/ν(ß), where σc(ß) has been obtained analytically as σc(ß)=10ß/(2ßeln10) for ß≥ßu=1/(2ln10), and for ß<ßu the weakest fiber failure leads to the catastrophic breakdown of the entire fiber bundle, similar to brittle materials, leading to σ_{c}(ß)=10-ß; (ii) the fraction of broken fibers right before the complete breakdown of the bundle has the form 1-1/(2ßln10); (iii) the distribution D(Δ) of the avalanches of size Δ follows a power-law D(Δ)∼Δ-ξ with ξ=5/2 for Δ≫Δc(ß) and ξ=3/2 for Δ≪Δc(ß), where the crossover avalanche size Δc(ß)=2/(1-e10-2ß)2.

Topological-distance-dependent transition in flocks with binary interactions.

We have studied a flocking model with binary interactions (binary flock), where the velocity of an agent depends on the velocity of only another agent and its own velocity, topped by the angular noise. The other agent is selected as the nth topological neighbor; the specific value of n being a fixed parameter of the problem. On the basis of extensive numerical simulation results, we argue that for n = 1, the phase transition from the ordered to the disordered phase of the flock is a special kind of discontinuous transition. Here, the order parameter does not flip-flop between multiple metastable states. It continues its initial disordered state for a period t(c), then switches over to the ordered state and remains in this state ever after. For n = 2, it is the usual discontinuous transition between two metastable states. Beyond this range, the continuous transitions are observed for n≥3. Such a system of binary flocks has been further studied using the hydrodynamic equations of motion. Linear stability analysis of the homogeneous polarized state shows that such a state is unstable close to the critical point and above some critical speed, which increases as we increase n. The critical noise strengths, which depend on the average correlation between a pair of topological neighbors, are estimated for five different values of n, which match well with their simulated values.

Interacting particles in a periodically moving potential: traveling wave and transport.

We study a system of interacting particles in a periodically moving external potential, within the simplest possible description of paradigmatic symmetric exclusion process on a ring. The model describes diffusion of hardcore particles where the diffusion dynamics is locally modified at a uniformly moving defect site, mimicking the effect of the periodically moving external potential. The model, though simple, exhibits remarkably rich features in particle transport, such as polarity reversal and double peaks in particle current upon variation of defect velocity and particle density. By tuning these variables, the most efficient transport can be achieved in either direction along the ring. These features can be understood in terms of a traveling density wave propagating in the system. Our results could be experimentally tested, e.g., in a system of colloidal particles driven by a moving optical tweezer.

Space-filling percolation.

A region of two-dimensional space has been filled randomly with a large number of growing circular disks allowing only a "slight" overlapping among them just before their growth stops. More specifically, each disk grows from a nucleation center that is selected at a random location within the uncovered region. The growth rate δ is a continuously tunable parameter of the problem which assumes a specific value while a particular pattern of disks is generated. When a growing disk overlaps for the first time with at least one other disk, its growth is stopped and is said to be frozen. In this paper we study the percolation properties of the set of frozen disks. Using numerical simulations we present evidence for the following: (i) The order parameter appears to jump discontinuously at a certain critical value of the area coverage; (ii) the width of the window of the area coverage needed to observe a macroscopic jump in the order parameter tends to vanish as δ→0; and on the contrary (iii) the cluster size distribution has a power-law-decaying functional form. While the first two results are the signatures of a discontinuous transition, the third result is indicative of a continuous transition. Therefore we refer to this transition as a sharp but continuous transition similar to what has been observed in the recently introduced Achlioptas process of explosive percolation. It is also observed that in the limit of δ→0, the critical area coverage at the transition point tends to unity, implying that the limiting pattern is space filling. In this limit, the fractal dimension of the pore space at the percolation point has been estimated to be 1.42(10) and the contact network of the disk assembly is found to be a scale-free network.

Scaling forms for relaxation times of the fiber bundle model.

Using extensive numerical analysis of the fiber bundle model with equal load sharing dynamics we studied the finite-size scaling forms of the relaxation times against the deviations of applied load per fiber from the critical point. Our most crucial result is we have not found any ln(N) dependence of the average relaxation time in the precritical state. The other results are as follows: (i) The critical load σ(c)(N) for the bundle of size N approaches its asymptotic value σ(c)(∞) as σ(c)(N)=σ(c)(∞)+AN(-1/ν). (ii) Right at the critical point the average relaxation time scales with the bundle size N as ~N(η) and this behavior remains valid within a small window of size |Δσ|~N(-ζ) around the critical point. (iii) When 1/N<|Δσ|<100N(-ζ) the finite-size scaling takes the form /N(η)~G[{σ(c)(N)-σ}N(ζ)] so in the limit of N→∞ one has ~(σ-σ(c))(-τ). The high precision of our numerical estimates led us to verify that ν=3/2, conjecture that η=1/3, ζ=2/3, and, therefore, τ=1/2.

Information sharing and sorting in a community.

We present the results of a detailed numerical study of a model for the sharing and sorting of information in a community consisting of a large number of agents. The information gathering takes place in a sequence of mutual bipartite interactions where randomly selected pairs of agents communicate with each other to enhance their knowledge and sort out the common information. Although our model is less restricted compared to the well-established naming game, the numerical results strongly indicate that the whole set of exponents characterizing this model are different from those of the naming game and they assume nontrivial values. Finally, it appears that in analogy to the emergence of clusters in the phenomenon of percolation, one can define clusters of agents here having the same information. We have studied in detail the growth of the largest cluster in this article and performed its finite-size scaling analysis.

Weighted trade network in a model of preferential bipartite transactions.

Using a model of wealth distribution where traders are characterized by quenched random saving propensities and trade among themselves by bipartite transactions, we mimic the enhanced rates of trading of the rich by introducing the preferential selection rule using a pair of continuously tunable parameters. The bipartite trading defines a growing trade network of traders linked by their mutual trade relationships. With the preferential selection rule this network appears to be highly heterogeneous characterized by the scale-free nodal degree and the link weight distributions and presents signatures of nontrivial strength-degree correlations. With detailed numerical simulations and using finite-size scaling analysis we present evidence that the associated critical exponents are continuous functions of the tuning parameters. However the wealth distribution has been observed to follow the well-known Pareto law robustly for all positive values of the tuning parameters.

Branching process in a stochastic extremal model.

We considered a stochastic version of the Bak-Sneppen model (SBSM) of ecological evolution where the number M of sites mutated in a mutation event is restricted to only two. Here the mutation zone consists of only one site and this site is randomly selected from the neighboring sites at every mutation event in an annealed fashion. The critical behavior of the SBSM is found to be the same as the BS model in dimensions d=1 and 2. However on the scale-free graphs the critical fitness value is nonzero even in the thermodynamic limit but the critical behavior is mean-field like. Finally M has been made even smaller than two by probabilistically updating the mutation zone, which also shows the original BS model behavior. We conjecture that a SBSM on any arbitrary graph with any small branching factor greater than unity will lead to a self-organized critical state.

A possible participation of gonadotropin-releasing hormone in the neuroleptic and cataleptic effect of haloperidol.

Haloperidol, an antipsychotic agent, stimulates the release of gonadotropin-releasing hormone (GnRH), and this hormone is known to mimic some of the behavioral effects of haloperidol. Hence, the present study was carried out to find out the contribution of GnRH in the behavioral effects of haloperidol. The studies revealed that haloperidol (0.15, 0.25 and 0.5 mg/kg, i.p.) and leuprolide (GnRH agonist; 50, 100, 200 and 400 microg/kg, s.c.) dose-dependently inhibited conditioned avoidance response (CAR) in male Sprague-Dawley rats. In higher doses, haloperidol (0.5, 1 mg/kg, i.p.) and leuprolide (200, 400 microg/kg, s.c.) produced catalepsy in rats. Co-administration of sub-effective dose of leuprolide (50 or 100 microg/kg, s.c.) and haloperidol (0.15 or 0.5 mg/kg, i.p.) similarly inhibited CAR and induced catalepsy. Pre-treatment of rats with antide (GnRH antagonist; 10 microg/rat, s.c.), attenuated the inhibitory effect of both the agents on CAR; blocked leuprolide-induced catalepsy; and attenuated the intensity and reduced the duration of haloperidol-induced catalepsy. In conclusion, the studies suggest a possible role of GnRH in the neuroleptic and cataleptic effect of haloperidol.

Weighted scale-free networks in Euclidean space using local selection rule.

A spatial scale-free network is introduced and studied, whose motivation originated in the growing Internet as well as airport networks. We argue that in these real-world networks a new node necessarily selects one of its neighboring local nodes for connection and is not controlled by preferential attachment as in the Barabási-Albert (BA) model. This observation is mimicked in our model where the nodes pop up at randomly located positions in the Euclidean space and are connected to one end of the nearest link. In spite of this crucial difference it is observed that the leading behavior of our network is like that of the BA model. Defining the link weight as an algebraic power of its Euclidean length, the weight distribution and the nonlinear dependence of the nodal strength on the degree are analytically calculated. It is claimed that a power law decay of the link weights with time ensures such nonlinear behavior. Switching off the Euclidean space from the same model yields a much simpler definition of the BA model where numerical effort grows linearly with N.

Phase transition in a directed traffic flow network.

The generic feature of traffic in a network of flowing electronic data packets is a phase transition from a stationary free-flow phase to a continuously growing congested nonstationary phase. In the most simple network of directed oriented square lattice we have been able to observe all crucial features of such flow systems having nontrivial critical behavior near the critical point of transition. The network here is in the shape of a square lattice and data packets are randomly posted with a rate rho at one side of the lattice. Each packet executes a directed diffusive motion toward the opposite boundary where it is delivered. Packets accumulated at a particular node form a queue and a maximum of m such packets randomly jump out of this node at every time step to its neighbors on a first-in-first-out basis. The phase transition occurs at rho(c) = m. The distribution of travel times through the system is found to have a log-normal behavior and the power spectrum of the load time series shows 1/f-like noise similar to the scenario of Internet traffic.