Sungchul Kwon

Condensation phenomena of a conserved-mass aggregation model on weighted complex networks

We investigate the condensation phase transitions of the conserved-mass aggregation (CA) model on weighted scale-free networks (WSFNs). In WSFNs, the weight w_{ij} is assigned to the link between the nodes i and j . We consider the symmetric weight given by w_{ij}=(k_{i}k_{j});{alpha} . On WSFNs, we numerically show that a certain critical alpha_{c} exists below which the CA model undergoes the same type of condensation transitions as those of the CA model on regular lattices. However, for alpha > or = alpha_{c} , the condensation always occurs for any density rho and omega . We analytically find alpha_{c}=(gamma-3)/2 on the WSFN with the degree exponent gamma . To obtain alpha_{c} , we analytically derive the scaling behavior of the stationary probability distribution P_{k};{infinity} of finding a walker at nodes with degree k , and the probability D(k) of finding two walkers simultaneously at the same node with degree k . We find P_{k};{infinity} approximately k;{alpha+1-gamma} and D(k) approximately k;{2(alpha+1)-gamma} , respectively. With P_{k};{infinity} , we also show analytically and numerically that the average mass m(k) on a node with degree k scales as k;{alpha+1} without any jumps at the maximal degree of the network for any rho as in SFNs with alpha=0 .

Condensation phase transitions of symmetric conserved-mass aggregation model on complex networks

We investigate condensation phase transitions of the symmetric conserved-mass aggregation (SCA) model on random networks (RNs) and scale-free networks (SFNs) with degree distribution P(k) approximately k(-gamma). In the SCA model, masses diffuse with unit rate, and unit mass chips off from mass with rate omega. The dynamics conserves total mass density rho. In the steady state, on RNs and SFNs with gamma > 3 for omega is not equal to infinity, we numerically show that the SCA model undergoes the same type of condensation transitions as those on regular lattices. However, the critical line rho(c)(omega) depends on network structures. On SFNs with gamma < or = 3, the fluid phase of exponential mass distribution completely disappears and no phase transitions occurs. Instead, the condensation with exponentially decaying background mass distribution always takes place for any nonzero density. For the existence of the condensed phase for gamma < or = 3 at the zero density limit, we investigate one lamb-lion problem on RNs and SFNs. We numerically show that a lamb survives indefinitely with finite survival probability on RNs and SFNs with gamma > 3, and dies out exponentially on SFNs with gamma< or = 3. The finite lifetime of a lamb on SFNs with gamma < or = 3 ensures the existence of the condensation at the zero density limit on SFNs with gamma < or = 3, at which direct numerical simulations are practically impossible. At omega = infinity, we numerically confirm that complete condensation takes place for any rho > 0 on RNs. Together with the recent study on SFNs, the complete condensation always occurs on both RNs and SFNs in zero range process with constant hopping rate.

Reentrant phase diagram of branching annihilating random walks with one and two offspring

We investigate the phase diagram of branching annihilating random walks with one and two offspring in one dimension. A walker can hop to a nearest neighbor site or branch with one or two offspring with relative ratio. Two walkers annihilate immediately when they meet. In general, this model exhibits a continuous phase transition from an active state into the absorbing state (vacuum) at a finite hopping probability. We map out the phase diagram by Monte Carlo simulations that shows a reentrant phase transition from vacuum to an active state and finally into vacuum again as the relative rate of the two-offspring branching process increases. This reentrant property apparently contradicts the conventional wisdom that increasing the number of offspring will tend to make the system more active. We show that the reentrant property is due to the static reflection symmetry of two-offspring branching processes and the conventional wisdom is recovered when the dynamic reflection symmetry is introduced instead of the static one. (c) 1995 The American Physical Society

Does hard core interaction change absorbing-type critical phenomena?

It has been generally believed that hard core interaction is irrelevant to absorbing-type critical phenomena because the particle density is so low near an absorbing phase transition. We study the effect of hard core interaction on the N-species branching annihilating random walks with two offspring and report that hard core interaction drastically changes the absorbing-type critical phenomena in a nontrivial way. Through a Langevin equation-type approach, we predict analytically the values of the scaling exponents, nu( perpendicular) = 2, z = 2, alpha = 1/2, and beta = 2 in one dimension for all N>1. Direct numerical simulations confirm our prediction. When the diffusion coefficients for different species are not identical, nu( perpendicular) and beta vary continuously with the ratios between the coefficients.

Critical phenomena of nonequilibrium dynamical systems with two absorbing states

We study nonequilibrium dynamical models with two absorbing states: interacting monomer-dimer models, probabilistic cellular automata models, nonequilibrium kinetic Ising models. These models exhibit a continuous phase transition from an active phase into an absorbing phase which belongs to the universality class of the models with the parity conservation. However, when we break the symmetry between the absorbing states by introducing a symmetry-breaking field, Monte Carlo simulations show that the system goes back to the conventional directed percolation universality class. In terms of domain wall language, the parity conservation is not affected by the presence of the symmetry-breaking field. So the symmetry between the absorbing states rather than the conservation laws plays an essential role in determining the universality class. We also perform Monte Carlo simulations for the various interface dynamics between different absorbing states, which yield new universal dynamic exponents. With the symmetry-breaking field, the interface moves, in average, with a constant velocity in the direction of the unpreferred absorbing state and the dynamic scaling exponents apparently assume trivial values. However, we find that the hyperscaling relation for the directed percolation universality class is restored if one focuses on the dynamics of the interface on the side of the preferred absorbing state only.

Continuously varying exponents in A+B-->0 reaction with long-ranged attractive interaction.

We investigate kinetics of A+B-->0 reaction with long-range attractive interaction V(r) approximately -r(-2sigma) between A and B or with drift velocity v approximately r(-sigma) in one dimension, where r is the closest distance between A and B . It is analytically shown that dynamical exponents for density of particles (rho) and size of domains (l) continuously vary with sigma when sigma < sigma(c) = 1/2 , while that for the distance between adjacent opposite species (l(AB)) varies when sigma < sigma(c)AB = 7/6 . For sigma > sigma(c)AB diffusive motions dominate the kinetics. These anomalous behaviors with the two crossover values of sigma are supported by numerical simulations.

Hyperuniformity of initial conditions and critical decay of a diffusive epidemic process belonging to the Manna class

For a fixed-energy Manna sandpile model belonging to a Manna class in one dimension (d=1), we recently showed that the critical decay is different for random and regular initial conditions (ICs). Compared with previous results of natural IC for several models, we suggested for the Manna class that the critical decay depends on the characteristics of the three ICs. But the dependence on the random and regular ICs was shown only for a single model. In this work, we study the critical decay for the random and regular ICs for another model of the Manna class in d=1, a diffusive epidemic process. It is shown that the critical decay exponent agrees with the previous result for each IC, which verifies that IC dependence is a common feature of the Manna class. In addition, for the random and regular ICs, we measure the variance σ^{2}(r) of total particle density in a region of size r by increasing r up to system size and investigate its temporal evolution toward the value σ_{q}^{2}(r) of the quasisteady state at criticality. In d=1,σ^{2}(r) scales as σ^{2}(r)∼r^{-ψ} with ψ=1 for random distributions and 1<ψ≤2 for hyperuniform ones. The temporal evolution shows that σ^{2}(r) of the two ICs differently relax toward σ_{q}^{2}(r) and the regular IC becomes a hyperuniform distribution of ψ=2 in the beginning of the evolution. We estimate ψ=1.45(3) for both the quasisteady state and absorbing states, so the quasisteady state is also as hyperuniform as absorbing states. The hyperuniformity of the quasisteady state shows that the natural IC also should be hyperuniform as much as the quasisteady state, because the natural IC is obtained from particle configurations close to the quasisteady state. Consequently, the different ψ of the three ICs suggest that σ^{2}(r) can classify the characteristics of the three ICs in a unified way and the different degree of hyperuniformity of the ICs provides another explanation for the observed IC-dependent critical decay in a point of view of initial fluctuations and correlations.

Anomalous kinetics of attractive A +B→0 reactions

We investigate the kinetics of the A + B-->0 reaction with the attractive interaction between opposite species in one spatial dimension. The attractive interaction leads to isotropic diffusions inside segregated single species domains, and accelerates the reactions of opposite species at the domain boundaries. At equal initial densities of and , we analytically and numerically show that the density of particles (rho), the size of domains (l), the distance between the closest neighbor of same species (lAA), and the distance between adjacent opposite species (lAB) scale in time as rho approximately t(-1/3), lAA approximately t(1/3), and l approximately lAB approximately lAB(2/3), respectively. These dynamical exponents define critical behavior distinguished from the class of uniformly driven systems of hard-core particles.

Absence of absorbing phase transitions in a conserved lattice-gas model in one dimension.

A one-dimensional conserved lattice-gas model is known to undergo continuous absorbing phase transitions where some of the critical exponents are exactly known. In one dimension, we recently showed that the model is mapped onto a two species reaction A+B→0 with diffusion rate of D_{A}>0 and D_{B}=0. In this work, it is explicitly shown from the scaling theory for A+B→0 that the observed scaling behavior of the conserved lattice-gas model is not associated with the absorbing phase transitions. Instead, the model indeed undergoes a crossover between two different scaling behaviors of A+B→0, the scaling behaviors of equal and unequal initial densities of two species. The crossover is similar to the absorbing transitions in many respects but some important features of continuous transitions such as the diverging fluctuations of an order parameter are absent.

Critical behavior for random initial conditions in the one-dimensional fixed-energy Manna sandpile model.

A fixed-energy Manna sandpile model undergoes an absorbing phase transition at a critical ρ_{c}, where an order parameter ϕ(t) decays as t^{-α} in time t. As the prototype of the Manna class, the model has been extensively studied in one dimension. However, the previous estimates of ρ_{c} and some critical exponents are different, depending on the types of initial conditions; random, natural, and regular conditions. The estimates of ρ_{c} for the random and the regular conditions are the lower and the upper bound among currently known estimates, respectively. In this work, for the random conditions, ρ_{c} and α are measured by taking into account finite-size (FS) effects. At the previous estimate of ρ_{c}, simulation results show that the temporal decay of ϕ(t) is strongly affected by the FS effects up to much larger system size (∼10^{6}). For the sizes for which ϕ(t) is independent up to t=2×10^{7}, we estimate ρ_{c}=0.8925(1) and α=0.110(5), which clearly differ from the previous results for the random conditions, ρ_{c}=0.89199(5) and α=0.141(24). Instead, the present ρ_{c} agrees with ρ_{c}=0.89255(2) of the regular conditions. In addition, the present α is substantially distinguishable from the results of the other types of initial conditions, α=0.159(3) and 0.146(2) for the natural and the regular conditions, respectively, which supports the claim of the initial condition dependence of dynamical exponents in the Manna class.