BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ronald Dickman

Ronald Dickman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.dis-nncond-mat.softPopulations and EvolutionBiological Physicscond-matphysics.data-anchem-phnlin.CDphysics.chem-phphysics.comp-phq-bio

Catalog footprint

What is connected

42works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2021arXiv

Contact processes on fragmented domains

Motivated by recent findings of enhanced species survival when fragmented habitats are reconnected through narrow strips of land [S. Pimm, and C. N. Jenkins, Am. Sci. {\bf 107}(3), 162 (2019).], we study the effect of a corridor connecting two square regions on the survival time of three models exhibiting extinction/survival phase transitions: the basic contact process (CP), the diffusive contact process, and the two-species symbiotic contact process (2SCP). We find that connecting the regions generally increases the survival time for $λ\ge λ_c$, where $λ$ is the reproduction rate and $λ_c$ its critical value. The enhancement of the survival time increases with $λ$, and is largest in the 2SCP.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multirange Ising model on the square lattice

We study the Ising model on $\mathbb{Z}^{2}$ and show, via numerical simulation, that allowing interactions between spins separated by distances $1$ and $m$ (two ranges), the critical temperature, $ T_c (m) $, converges monotonically to the critical temperature of the Ising model on $\mathbb{Z}^4$ as $ m \to \infty $. Only interactions between spins located in directions parallel to each coordinate axis are considered. We also simulated the model with interactions between spins at distances of $ 1 $, $ m $ and $ u $ (three ranges), with $ u $ a multiple of $ m $; in this case our results indicate that $ T_c(m, u) $ converges to the critical temperature of the model on $ \mathbb{Z}^6$. For percolation, analogous results were proven for the critical probability $p_c$ [B. N. B. de Lima, R. P. Sanchis and R. W. C. Silva, Stochastic Process. Appl. {\bf 121}, 2043 (2011)].

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase coexistence far from equilibrium

Investigation of simple far-from-equilibrium systems exhibiting phase separation leads to the conclusion that phase coexistence is not well defined in this context. This is because the properties of the coexisting nonequilibrium systems depend on how they are placed in contact, as verified in the driven lattice gas with attractive interactions, and in the two-temperature lattice gas, under (a) weak global exchange between uniform systems, and (b) phase-separated (nonuniform) systems. Thus, far from equilibrium, the notions of universality of phase coexistence (i.e., independence of how systems exchange particles and/or energy), and of phases with intrinsic properties (independent of their environment) are lost.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spatio-temporal generalization of the Harris criterion and its application to diffusive disorder

We investigate how a clean continuous phase transition is affected by spatio-temporal disorder, i.e., by an external perturbation that fluctuates in both space and time. We derive a generalization of the Harris criterion for the stability of the clean critical behavior in terms of the space-time correlation function of the external perturbation. As an application, we consider diffusive disorder, i.e, an external perturbation governed by diffusive dynamics, and its effects on a variety of equilibrium and nonequilibrium critical points. We also discuss the relation between diffusive disorder and diffusive dynamical degrees of freedom on the example of model C of the Hohenberg-Halperin classification, and we comment on Griffiths singularities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Griffiths phases and localization in hierarchical modular networks

We study variants of hierarchical modular network models suggested by Kaiser and Hilgetag [Frontiers in Neuroinformatics, 4 (2010) 8] to model functional brain connectivity, using extensive simulations and quenched mean-field theory (QMF), focusing on structures with a connection probability that decays exponentially with the level index. Such networks can be embedded in two-dimensional Euclidean space. We explore the dynamic behavior of the contact process (CP) and threshold models on networks of this kind, including hierarchical trees. While in the small-world networks originally proposed to model brain connectivity, the topological heterogeneities are not strong enough to induce deviations from mean-field behavior, we show that a Griffiths phase can emerge under reduced connection probabilities, approaching the percolation threshold. In this case the topological dimension of the networks is finite, and extended regions of bursty, power-law dynamics are observed. Localization in the steady state is also shown via QMF. We investigate the effects of link asymmetry and coupling disorder, and show that localization can occur even in small-world networks with high connectivity in case of link disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Activity, diffusion, and correlations in a two-dimensional conserved stochastic sandpile

We perform large-scale simulations of a two-dimensional restricted-height conserved stochastic sandpile, focusing on particle diffusion and mobility, and spatial correlations. Quasistationary (QS) simulations yield the critical particle density to high precision [$p_c = 0.7112687(2)$], and show that the diffusion constant scales in the same manner as the activity density, as found previously in the one-dimensional case. Short-time scaling is characterized by subdiffusive behavior (mean-square displacement $\sim t^γ$ with $γ< 1$), which is easily understood as a consequence of the initial decay of activity, $ρ(t) \sim t^{-δ}$, with $γ= 1- δ$. We verify that at criticality, the activity correlation function $C(r) \sim r^{-β/ν_\perp}$, as expected at an absorbing-state phase transition. Our results for critical exponents are consistent with, and somewhat more precise than, predictions derived from the Langevin equation for stochastic sandpiles in two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Failure of steady state thermodynamics in lattice gases under nonuniform drive

To be useful, steady state thermodynamics (SST) must be self-consistent and have predictive value. Although consistency of SST was recently verified for driven lattice gases under global weak exchange, I show here that it does not predict the coexisting densities in athermal stochastic lattice gases under a nonuniform drive. I consider the lattice gas with nearest-neighbor exclusion on the square lattice, with nearest-neighbor hopping (NNE dynamics), and with hopping to both nearest and next-nearest neighbors (NNE2 dynamics). Part of the system is subject to a drive $D$ that favors hopping along one direction, while the other part is free of the drive. Thus the steady state represents coexistence between two subsystems, one far from equilibrium and the other in equilibrium, which exchange particles along the interfaces. The dimensionless chemical potential $μ^*(ρ,D) = μ/k_B T$ for the lattice gas with density $ρ$ is readily determined in studies of a uniform system. Under the nonuniform drive, however, equating the chemical potentials of the coexisting subsystems does not yield the coexisting densities. The steady state chemical potential is, moreover, different in the coexisting bulk regions, contrary to the basic principles of thermodynamics. These results cast serious doubt on the predictive value of SST.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inconsistencies in steady state thermodynamics

We address the issue of extending thermodynamics to nonequilibrium steady states. Using driven stochastic lattice gases, we ask whether consistent definitions of an effective chemical potential mu, and an effective temperature T_e, are possible. mu and T_e are determined via coexistence, i.e., zero flux of particles and energy between the driven system and a reservoir. In the lattice gas with nearest-neighbor exclusion, temperature is not relevant, and we find that the effective chemical potential, a function of density and drive strength, satisfies the zeroth law, and correctly predicts the densities of coexisting systems. In the Katz-Lebowitz-Spohn driven lattice gas both mu and T_e need to be defined. We show analytically that in this case the zeroth law is violated for Metropolis exchange rates, and determine the size of the violations numerically. The zeroth law appears to be violated for generic exchange rates. Remarkably, the system-reservoir coupling proposed by Sasa and Tasaki [J. Stat. Phys. 125, 125 (2006)], is free of inconsistencies, and the zeroth law holds. This is because the rate depends only on the state of the donor system, and is independent of that of the acceptor.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Intrinsic convergence properties of entropic sampling algorithms

We study the convergence of the density of states and thermodynamic properties in three flat-histogram simulation methods, the Wang-Landau (WL) algorithm, the 1/t algorithm, and tomographic sampling (TS). In the first case the refinement parameter f is rescaled (f -> f/2) each time the flat-histogram condition is satisfied, in the second f ~ 1/t after a suitable initial phase, while in the third f is constant (t corresponds to Monte Carlo time). To examine the intrinsic convergence properties of these methods, free of any complications associated with a specific model, we study a featureless entropy landscape, such that for each allowed energy E = 1,...,L, there is exactly one state, that is, g(E) = 1 for all E. Convergence of sampling corresponds to g(E,t) -> const. as t -> infinity, so that the standard deviation sigma_g of g over energy values is a measure of the overall sampling error. Neither the WL algorithm nor TS converge: in both cases sigma_g saturates at long times. In the 1/t algorithm, by contrast, sigma_g decays \propto 1/\sqrt{t}. Modified TS and 1/t procedures, in which f \propto 1/t^alpha, converge for alpha values between 0 and 1. There are two essential facets to convergence of flat-histogram methods: elimination of initial errors in g(E), and correction of the sampling noise accumulated during the process. For a simple example, we demonstrate analytically, using a Langevin equation, that both kinds of errors can be eliminated, asymptotically, if f ~ 1/t^alpha with 0 < α\leq 1. Convergence is optimal for alpha = 1. For alpha \leq 0 the sampling noise never decays, while for alpha > 1 the initial error is never completely eliminated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase diagram of the symbiotic two-species contact process

We study the two-species symbiotic contact process (2SCP), recently proposed in [de Oliveira, Santos and Dickman, Phys. Rev. E {\bf 86}, 011121 (2012)] . In this model, each site of a lattice may be vacant or host single individuals of species A and/or B. Individuals at sites with both species present interact in a symbiotic manner, having a reduced death rate, $μ< 1$. Otherwise, the dynamics follows the rules of the basic CP, with individuals reproducing to vacant neighbor sites at rate $λ$ and dying at a rate of unity. We determine the full phase diagram in the $λ-μ$ plane in one and two dimensions by means of exact numerical quasistationary distributions, cluster approximations, and Monte Carlo simulations. We also study the effects of asymmetric creation rates and diffusion of individuals. In two dimensions, for sufficiently strong symbiosis (i.e., small $μ$), the absorbing-state phase transition becomes discontinuous for diffusion rates $D$ within a certain range. We report preliminary results on the critical surface and tricritical line in the $λ-μ-D$ space. Our results raise the possibility that strongly symbiotic associations of mobile species may be vulnerable to sudden extinction under increasingly adverse conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conserved sandpile with a variable height restriction

We study a restricted-height version of the one-dimensional Oslo sandpile with conserved density, using periodic boundary conditions. Each site has a limiting height which can be either two or three. When a site reaches its limiting height it becomes active and may topple, loosing two particles, which move randomly to nearest-neighbor sites. After a site topples it is randomly assigned a new limiting height. We study the model using mean-field theory and Monte Carlo simulation, focusing on the quasi-stationary state, in which the number of active sites fluctuates about a stationary value. Using finite-size scaling analysis, we determine the critical particle density and associated critical exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Survival of the Scarcer in Space

The dynamics leading to extinction or coexistence of competing species is of great interest in ecology and related fields. Recently a model of intra- and interspecific competition between two species was proposed by Gabel et al. [Phys. Rev. E 87 (2013) 010101], in which the scarcer species (i.e., with smaller stationary population size) can be more resistant to extinction when it holds a competitive advantage; the latter study considered populations without spatial variation. Here we verify this phenomenon in populations distributed in space. We extend the model of Gabel et al. to a d-dimensional lattice, and study its population dynamics both analytically and numerically. Survival of the scarcer in space is verified for situations in which the more competitive species is closer to the threshold for extinction than is the less competitive species, when considered in isolation. The conditions for survival of the scarcer species, as obtained applying renormalization group analysis and Monte Carlo simulation, differ in detail from those found in the spatially homogeneous case. Simulations highlight the speed of invasion waves in determining the survival times of the competing species.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A symbiotic two-species contact process

We study a contact process (CP) with two species that interact in a symbiotic manner. In our model, each site of a lattice may be vacant or host individuals of species A and/or B; multiple occupancy by the same species is prohibited. Symbiosis is represented by a reduced death rate, mu < 1, for individuals at sites with both species present. Otherwise, the dynamics is that of the basic CP, with creation (at vacant neighbor sites) at rate lambda and death of (isolated) individuals at a rate of unity. Mean-field theory and Monte Carlo simulation show that the critical creation rate, lambda_c (mu), is a decreasing function of mu, even though a single-species population must go extinct for lambda < lambda_c(1), the critical point of the basic CP. Extensive simulations yield results for critical behavior that are compatible with the directed percolation (DP) universality class, but with unusually strong corrections to scaling. A field-theoretical argument supports the conclusion of DP critical behavior. We obtain similar results for a CP with creation at second-neighbor sites and enhanced survival at first neighbors, in the form of an annihilation rate that decreases with the number of occupied first neighbors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An epidemic process mediated by a decaying diffusing signal

We study a stochastic epidemic model consisting of elements (organisms in a community or cells in tissue) with fixed positions, in which damage or disease is transmitted by diffusing agents ("signals") emitted by infected individuals. The signals decay as well as diffuse; since they are assumed to be produced in large numbers, the signal concentration is treated deterministically. The model, which includes four cellular states (susceptible, transformed, depleted, and removed), admits various interpretations: spread of an infection or infectious disease, or of damage in a tissue in which injured cells may themselves provoke further damage, and as a description of the so-called radiation-induced bystander effect, in which the signals are molecules capable of inducing cell damage and/or death in unirradiated cells. The model exhibits a continuous phase transition between spreading and nonspreading phases. We formulate two mean-field theory (MFT) descriptions of the model, one of which ignores correlations between the cellular state and the signal concentration, and another that treats such correlations in an approximate manner. Monte Carlo simulations of the spread of infection on the square lattice yield values for the critical exponents and the fractal dimension consistent with the dynamic percolation universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analysis of an information-theoretic model for communication

We study the cost-minimization problem posed by Ferrer i Cancho and Solé in their model of communication that aimed at explaining the origin of Zipf's law [PNAS 100, 788 (2003)]. Direct analysis shows that the minimum cost is $\min {λ, 1-λ}$, where $λ$ determines the relative weights of speaker's and hearer's costs in the total, as shown in several previous works using different approaches. The nature and multiplicity of the minimizing solution changes discontinuously at $λ=1/2$, being qualitatively different for $λ< 1/2$, $λ> 1/2$, and $λ=1/2$. Zipf's law is found only in a vanishing fraction of the minimum-cost solutions at $λ= 1/2$ and therefore is not explained by this model. Imposing the further condition of equal costs yields distributions substantially closer to Zipf's law, but significant differences persist. We also investigate the solutions reached by the previously used minimization algorithm and find that they correctly recover global minimum states at the transition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Discontinuous phase transition in a dimer lattice gas

I study a dimer model on the square lattice with nearest-neighbor exclusion as the only interaction. Detailed simulations using tomographic entropic sampling show that as the chemical potential is varied, there is a strongly discontinuous phase transition, at which the particle density jumps by about 18% of its maximum value, 1/4. The transition is accompanied by the onset of orientational order, to an arrangement corresponding to the {1/2,0,1/2} structure identified by Phares et al. [Physica B {\bf 409}, 1096 (2011)] in a dimer model with {\it finite} repulsion at {\it fixed} density. Using finite-size scaling and Binder's cumulant, the expected scaling behavior at a discontinuous transition is verified in detail. The discontinuous transition can be understood qualitatively given that the model possesses eight equivalent maximum-density configurations, so that its coarse-grained description corresponds to that of the q=8 Potts model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An infinite-period phase transition versus nucleation in a stochastic model of collective oscillations

A lattice model of three-state stochastic phase-coupled oscillators has been shown by Wood et al (2006 Phys. Rev. Lett. 96 145701) to exhibit a phase transition at a critical value of the coupling parameter, leading to stable global oscillations. We show that, in the complete graph version of the model, upon further increase in the coupling, the average frequency of collective oscillations decreases until an infinite-period (IP) phase transition occurs, at which point collective oscillations cease. Above this second critical point, a macroscopic fraction of the oscillators spend most of the time in one of the three states, yielding a prototypical nonequilibrium example (without an equilibrium counterpart) in which discrete rotational (C_3) symmetry is spontaneously broken, in the absence of any absorbing state. Simulation results and nucleation arguments strongly suggest that the IP phase transition does not occur on finite-dimensional lattices with short-range interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Complete high-precision entropic sampling

Monte Carlo simulations using entropic sampling to estimate the number of configurations of a given energy are a valuable alternative to traditional methods. We introduce {\it tomographic} entropic sampling, a scheme which uses multiple studies, starting from different regions of configuration space, to yield precise estimates of the number of configurations over the {\it full range} of energies, {\it without} dividing the latter into subsets or windows. Applied to the Ising model on the square lattice, the method yields the critical temperature to an accuracy of about 0.01%, and critical exponents to 1% or better. Predictions for systems sizes L=10 - 160, for the temperature of the specific heat maximum, and of the specific heat at the critical temperature, are in very close agreement with exact results. For the Ising model on the simple cubic lattice the critical temperature is given to within 0.003% of the best available estimate; the exponent ratios $β/ν$ and $γ/ν$ are given to within about 0.4% and 1%, respectively, of the literature values. In both two and three dimensions, results for the {\it antiferromagnetic} critical point are fully consistent with those of the ferromagnetic transition. Application to the lattice gas with nearest-neighbor exclusion on the square lattice again yields the critical chemical potential and exponent ratios $β/ν$ and $γ/ν$ to good precision.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Conserved directed percolation: exact quasistationary distribution of small systems and Monte Carlo simulations

We study symmetric sleepy random walkers, a model exhibiting an absorbing-state phase transition in the conserved directed percolation (CDP) universality class. Unlike most examples of this class studied previously, this model possesses a continuously variable control parameter, facilitating analysis of critical properties. We study the model using two complementary approaches: analysis of the numerically exact quasistationary (QS) probability distribution on rings of up to 22 sites, and Monte Carlo simulation of systems of up to 32000 sites. The resulting estimates for critical exponents beta, beta/nu_perp, and z, and the moment ratio m_{211} = < rho^2 >/< rho >^2 (rho is the activity density), based on finite-size scaling at the critical point, are in agreement with previous results for the CDP universality class. We find, however, that the approach to the QS regime is characterized by a different value of the dynamic exponent z than found in the QS regime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Contact process with sublattice symmetry breaking

We study a contact process with creation at first- and second-neighbor sites and inhibition at first neighbors, in the form of an annihilation rate that increases with the number of occupied first neighbors. Mean-field theory predicts three phases: inactive (absorbing), active symmetric, and active asymmetric, the latter exhibiting distinct sublattice densities on a bipartite lattice. These phases are separated by continuous transitions; the phase diagram is reentrant. Monte Carlo simulations in two dimensions verify these predictions qualitatively, except for a first-neighbor creation rate of zero. (In the latter case one of the phase transitions is discontinuous.) Our numerical results confirm that the symmetric-asymmetric transition belongs to the Ising universality class, and that the active-absorbing transition belongs to the directed percolation class, as expected from symmetry considerations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Glassy dynamics and hysteresis in a linear system of orientable hard rods

We study the dynamics of a one-dimensional fluid of orientable hard rectangles with a non-coarse-grained microscopic mechanism of facilitation. The length occupied by a rectangle depends on its orientation, which is coupled to an external field. The equilibrium properties of our model are essentially those of the Tonks gas, but at high densities, the orientational degrees of freedom become effectively frozen due to jamming. This is a simple analytically tractable model of glassy phase. Under a cyclic variation of the pressure, hysteresis is observed. Following a pressure quench, the orientational persistence exhibits a two-stage decay characteristic of glassy systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pair contact process with diffusion of pairs

The pair contact process (PCP) is a nonequilibrium stochastic model which, like the basic contact process (CP), exhibits a phase transition to an absorbing state. The two models belong to the directed percolation (DP) universality class, despite the fact that the PCP possesses infinitely many absorbing configurations whereas the CP has but one. The critical behavior of the PCP with hopping by particles (PCPD) is as yet unclear. Here we study a version of the PCP in which nearest-neighbor particle {\it pairs} can hop but individual particles cannot. Using quasistationary simulations for three values of the diffusion probability ($D=0.1$, 0.5 and 0.9), we find convincing evidence of DP-like critical behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Systematic enumeration of configuration classes for entropic sampling of Ising models

We describe a systematic method for complete enumeration of configuration classes (CCs) of the spin-1/2 Ising model in the energy-magnetization plane. This technique is applied to the antiferromagnetic Ising model in an external magnetic field on the square lattice, which is simulated using the tomographic entropic sampling algorithm. We estimate the number of configurations, $Ω(E,m,L)$, and related microcanonical averages, for all allowed energies $E$ and magnetizations $m$ for $L = 10$ to 30, with $ΔL = 2$. With prior knowledge of the CCs, we can be sure that all allowed classes are sampled during the simulation. Complete enumeration of CCs also enables us to use the final estimate of $Ω(E,m,L)$ to obtain good initial estimates, $Ω_0(E,m,L')$, for successive system sizes ($L' > L$) through a two-dimensional interpolation. Using these results we calculate canonical averages of the thermodynamic quantities of interest as continuous functions of temperature $T$ and external field $h$. In addition, we determine the critical line in the $h$-$T$ plane using finite-size scaling analysis, and compare these results with several approximate theoretical expressions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymmetric Exclusion Process in a System of Interacting Brownian Particles

We study a continuous-space version of the totally asymmetric simple exclusion process (TASEP), consisting of interacting Brownian particles subject to a driving force in a periodic external potential. Particles are inserted at the leftmost site at rate $α$, hop to the right at unit rate, and are removed at the rightmost site at rate $β$. Our study is motivated by recent experiments on colloidal particles in optical tweezer arrays. The external potential is of the form generated by such an array. Particles spend most of the time near potential minima, approximating the situation in the lattice gas; a short-range repulsive interaction prevents two particles from occupying the same potential well. A constant driving force, representing Stokes drag on particles suspended in a moving fluid, leads to biased motion. Our results for the density profile and current, obtained via numerical integration of the Langevin equation and dynamic Monte Carlo simulations, indicate that the continuous-space model exhibits phase transitions analogous to those observed in the lattice model. The correspondence is not exact, however, due to the lack of particle-hole symmetry in our model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical behavior of hard-core lattice gases: Wang-Landau sampling with adaptive windows

Critical properties of lattice gases with nearest-neighbor exclusion are investigated via the adaptive-window Wang-Landau algorithm on the square and simple cubic lattices, for which the model is known to exhibit an Ising-like phase transition. We study the particle density, order parameter, compressibility, Binder cumulant and susceptibility. Our results show that it is possible to estimate critical exponents using Wang-Landau sampling with adaptive windows. Finite-size-scaling analysis leads to results in fair agreement with exact values (in two dimensions) and numerical estimates (in three dimensions).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagram and critical behavior of the pair annihilation model

We study the critical behavior of the pair annihilation model (PAM) with diffusion in one, two and three dimensions, using the pair approximation (PA) and Monte Carlo simulation. Of principal interest is the dependence of the critical creation rate, lambda_c, on the diffusion probability D, in particular, whether survival is possible at arbitrarily small creation rates, for sufficiently rapid diffusion. Whilst the PA predicts that in any spatial dimension d \geq 1, lambda_c -> 0 at some diffusion probability D* < 1, Katori and Konno [Physica A {\bf 186}, 578 (1992)] showed rigorously that for d \leq 2, one has lambda_c > 0 for any D<1. Our simulation results are consistent with this theorem. In two dimensions, the extinction region becomes narrow as D approaches unity, following lambda_c \propto \exp[- \rm{const.}/(1-D)^γ], with gamma = 1.41(2). In three dimensions we find D* = 0.333(3). The simulation results confirm that the PAM belongs to the directed percolation universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Contact process with mobile disorder

I study the absorbing-state phase transition in the one-dimensional contact process with mobile disorder. In this model the dilution sites, though permanently inactive, diffuse freely, exchanging positions with the other sites, which host a basic contact process. Even though the disorder variables are not quenched, the critical behavior is affected: the critical exponents delta and z, the ratio beta/nu_perp and the moment ratio m= <rho^2>/rho^2 take values different from those of directed percolation, and appear to vary with the vacancy diffusion rate. While the survival probability starting from a single active seed follows the usual scaling, P(t) ~ t^{-delta}, at the critical point, the mean number of active sites and mean-square spread grow more slowly than power laws. The critical creation rate increases with the vacancy density v and diverges at a value v_c < 1. The scaling behavior at this point appears to be simpler than for smaller vacancy densities.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Diffusion in stochastic sandpiles

We study diffusion of particles in large-scale simulations of one-dimensional stochastic sandpiles, in both the restricted and unrestricted versions. The results indicate that the diffusion constant scales in the same manner as the activity density, so that it represents an alternative definition of an order parameter. The critical behavior of the unrestricted sandpile is very similar to that of its restricted counterpart, including the fact that a data collapse of the order parameter as a function of the particle density is only possible over a very narrow interval near the critical point. We also develop a series expansion, in inverse powers of the density. for the collective diffusion coefficient in a variant of the stochastic sandpile in which the toppling rate at a site with $n$ particles is $n(n-1)$, and compare the theoretical prediction with simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the absorbing-state phase transition in the one-dimensional triplet creation model

We study the lattice reaction diffusion model 3A -> 4A, A -> 0 (``triplet creation") using numerical simulations and n-site approximations. The simulation results provide evidence of a discontinuous phase transition at high diffusion rates. In this regime the order parameter appears to be a discontinuous function of the creation rate; no evidence of a stable interface between active and absorbing phases is found. Based on an effective mapping to a modified compact directed percolation process, shall nevertheless argue that the transition is continuous, despite the seemingly discontinuous phase transition suggested by studies of finite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Functional-integral based perturbation theory for the Malthus-Verhulst process

We apply a functional-integral formalism for Markovian birth and death processes to determine asymptotic corrections to mean-field theory in the Malthus-Verhulst process (MVP). Expanding about the stationary mean-field solution, we identify an expansion parameter that is small in the limit of large mean population, and derive a diagrammatic expansion in powers of this parameter. The series is evaluated to fifth order using computational enumeration of diagrams. Although the MVP has no stationary state, we obtain good agreement with the associated {\it quasi-stationary} values for the moments of the population size, provided the mean population size is not small. We compare our results with those of van Kampen's $Ω$-expansion, and apply our method to the MVP with input, for which a stationary state does exist.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Asymmetric dynamics and critical behavior in the Bak-Sneppen model

We investigate, using mean-field theory and simulation, the effect of asymmetry on the critical behavior and probability density of Bak-Sneppen models. Two kinds of anisotropy are investigated: (i) different numbers of sites to the left and right of the central (minimum) site are updated and (ii) sites to the left and right of the central site are renewed in different ways. Of particular interest is the crossover from symmetric to asymmetric scaling for weakly asymmetric dynamics, and the collapse of data with different numbers of updated sites but the same degree of asymmetry. All non-symmetric rules studied fall, independent of the degree of asymmetry, in the same universality class. Conversely, symmetric variants reproduce the exponents of the original model. Our results confirm the existence of two symmetry-based universality classes for extremal dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

On singular probability densities generated by extremal dynamics

Extremal dynamics is the mechanism that drives the Bak-Sneppen model into a (self-organized) critical state, marked by a singular stationary probability density $p(x)$. With the aim of understanding this phenomenon, we study the BS model and several variants via mean-field theory and simulation. In all cases, we find that $p(x)$ is singular at one or more points, as a consequence of extremal dynamics. Furthermore we show that the extremal barrier $x_i$ always belongs to the `prohibited' interval, in which $p(x)=0$. Our simulations indicate that the Bak-Sneppen universality class is robust with regard to changes in the updating rule: we find the same value for the exponent $π$ for all variants. Mean-field theory, which furnishes an exact description for the model on a complete graph, reproduces the character of the probability distribution found in simulations. For the modified processes mean-field theory takes the form of a functional equation for $p(x)$.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Quasi-Stationary Distributions for Models of Heterogeneous Catalysis

We construct the quasi-stationary (QS) distribution for two models of heterogeneous catalysis having two absorbing states: the ZGB model for the oxidation of CO, and a version with noninstantaneous reactions. Using a mean-field-like approximation, we study the quasi-stationary surface coverages, moment ratios and the lifetime of the QS state. We also derive an improved, consistent one-site mean-field theory for the ZGB model.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Recoiling DNA Molecule: Simulation & Experiment

Single molecule DNA experiments often generate data from force versus extension measurements involving the tethering of a microsphere to one end of a single DNA molecule while the other is attached to a substrate. We show that the persistence length of single DNA molecules can also be measured based on the recoil dynamics of these DNA-microsphere complexes if appropriate corrections are made to the friction coefficient of the microsphere in the vicinity of the substrate. Comparison between computer simulated recoil curves, generated from the corresponding Langevin equation, and experimental recoils is used to assure the validity of data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Variable survival exponents in history-dependent random walks: hard movable reflector

We review recent studies demonstrating a nonuniversal (continuously variable) survival exponent for history-dependent random walks, and analyze a new example, the hard movable partial reflector. These processes serve as a simplified models of infection in a medium with a history-dependent susceptibility, and for spreading in systems with an infinite number of absorbing configurations. The memory may take the form of a history-dependent step length, or be the result of a partial reflector whose position marks the maximum distance the walker has ventured from the origin. In each case, a process with memory is rendered Markovian by a suitable expansion of the state space. Asymptotic analysis of the probability generating function shows that, for large t, the survival probability decays as S(t) \sim t^{-delta}, where δvaries with the parameters of the model. We report new results for a hard partial reflector, i.e., one that moves forward only when the walker does. When the walker tries to jump to the site R occupied by the reflector, it is reflected back with probability r, and stays at R with probability 1-r; only in the latter case does the reflector move (R \to R+1). For this model, delta = 1/2(1-r), and becomes arbitrarily large as r approaches 1. This prediction is confirmed via iteration of the transfer matrix, which also reveals slowly-decaying corrections to scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Asymptotic analysis of a random walk with a history-dependent step length

We study an unbiased, discrete time random walk on the nonnegative integers, with the origin absorbing. The process has a history-dependent step length: the walker takes steps of length v while in a region which has been visited before, and steps of length n when entering a region that has never been visited. The process provides a simplified model of spreading in systems with an infinite number of absorbing configurations. Asymptotic analysis of the probability generating function shows that, for large t, the survival probability decays as S(t) \sim t^{-delta}, with delta = v/2n. Our expression for the decay exponent is in agreement with results obtained via numerical iteration of the transition matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Path Integrals and Perturbation Theory for Stochastic Processes

We review and extend the formalism introduced by Peliti, that maps a Markov process to a path-integral representation. After developing the mapping, we apply it to some illustrative examples: the simple decay process, the birth-and-death process, and the Malthus-Verhulst process. In the first two cases we show how to obtain the exact probability generating function using the path integral. We show how to implement a diagrammatic perturbation theory for processes that do not admit an exact solution. Analysis of a set of coupled Malthus-Verhulst processes on a lattice leads, in the continuum limit, to a field theory for directed percolation and allied models.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Scaling Exponents of Rough Surfaces Generated by the Domany-Kinzel Cellular Automaton

The critical behavior at the frozen/active transition in the Domany-Kinzel stochastic cellular automaton (DKCA) is studied {\it via} a surface growth process in (1+1) dimensions. At criticality, this process presents a kinetic roughening transition; we measure the critical exponents in simulations. Two update schemes are considered: in the symmetric scheme, the growth surfaces belong to the Directed Percolation (DP) universality class, except at one terminal point. At this point, the phase transition is discontinuous and the surfaces belong to the Compact Directed Percolation universality class. The relabeling of space-time points in the nonsymmetric scheme alters the surface growth dramatically. The critical behavior of rough surfaces at the nonchaotic/chaotic transition is also studied using the damage spreading technique; the exponents confirm DP values for the symmetric scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Absorbing-state phase transitions in fixed-energy sandpiles

We study sandpile models as closed systems, with conserved energy density $ζ$ playing the role of an external parameter. The critical energy density, $ζ_c$, marks a nonequilibrium phase transition between active and absorbing states. Several fixed-energy sandpiles are studied in extensive simulations of stationary and transient properties, as well as the dynamics of roughening in an interface-height representation. Our primary goal is to identify the universality classes of such models, in hopes of assessing the validity of two recently proposed approaches to sandpiles: a phenomenological continuum Langevin description with absorbing states, and a mapping to driven interface dynamics in random media. Our results strongly suggest that there are at least three distinct universality classes for sandpiles.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Paths to Self-Organized Criticality

We present a pedagogical introduction to self-organized criticality (SOC), unraveling its connections with nonequilibrium phase transitions. There are several paths from a conventional critical point to SOC. They begin with an absorbing-state phase transition (directed percolation is a familiar example), and impose supervision or driving on the system; two commonly used methods are extremal dynamics, and driving at a rate approaching zero. We illustrate this in sandpiles, where SOC is a consequence of slow driving in a system exhibiting an absorbing-state phase transition with a conserved density. Other paths to SOC, in driven interfaces, the Bak-Sneppen model, and self-organized directed percolation, are also examined. We review the status of experimental realizations of SOC in light of these observations.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Sandpiles and absorbing-state phase transitions: recent results and open problems

We review some recent results on the relations between sandpiles and a class of absorbing state phase transitions. We use the concept of fixed energy sandpiles (FES), in which external driving and dissipation are absent. FES are shown to exhibit an absorbing state transition with critical properties coinciding with those of the corresponding sandpile model. We propose a set of Langevin equations capturing the relevant features of this transition. These equations characterize the universality class of systems with an infinite number of absorbing states and a static conserved field coupled to the order parameter. Different models in this class are identified, and strong evidence is presented showing that the Manna sandpile, as well as some other stochastic sandpiles, belong in this universality class. Finally some open problems and questions are discussed.

0 citations0 reviewsNo code linkedOpen access