BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ronald Dickman

Ronald Dickman contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechcond-mat.softcond-mat.dis-nnphysics.data-anPopulations and EvolutionBiological Physicsnlin.CDphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

21 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
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
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 (&#34;signals&#34;) 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&#39;s law [PNAS 100, 788 (2003)]. Direct analysis shows that the minimum cost is $\min {λ, 1-λ}$, where $λ$ determines the relative weights of speaker&#39;s and hearer&#39;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&#39;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&#39;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&#39;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&#39;)$, for successive system sizes ($L&#39; > 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&#34;) 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
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