SciPost Submission Page
Smoothly vanishing density in the contact process by an interplay of disorder and long-distance dispersal
by Róbert Juhász, István A. Kovács
This is not the latest submitted version.
Submission summary
| Authors (as registered SciPost users): | Róbert Juhász |
| Submission information | |
|---|---|
| Preprint Link: | scipost_202405_00018v1 (pdf) |
| Date submitted: | May 14, 2024, 10:28 a.m. |
| Submitted by: | Róbert Juhász |
| Submitted to: | SciPost Physics Core |
| Ontological classification | |
|---|---|
| Academic field: | Physics |
| Specialties: |
|
| Approaches: | Theoretical, Computational |
Abstract
Realistic modeling of ecological population dynamics requires spatially explicit descriptions that can take into account spatial heterogeneity as well as long-distance dispersal. Here, we present Monte Carlo simulations and numerical renormalization group results for the paradigmatic model, the contact process, in the combined presence of these factors in both one and two-dimensional systems. Our results confirm our analytic arguments stating that the density vanishes smoothly at the extinction threshold, in a way characteristic of infinite-order transitions. This extremely smooth vanishing of the global density entails an enhanced exposure of the population to extinction events. At the same time, a reverse order parameter, the local persistence displays a discontinuity characteristic of mixed-order transitions, as it approaches a non-universal critical value algebraically with an exponent $\beta_p'<1$.
Current status:
Reports on this Submission
Report #2 by Anonymous (Referee 1) on 2024-6-7 (Invited Report)
- Cite as: Anonymous, Report on arXiv:scipost_202405_00018v1, delivered 2024-06-07, doi: 10.21468/SciPost.Report.9196
Strengths
1) The work tackles an interesting problem: understanding the combined effect of quenched disorder and long-range dispersal on the dynamical transition in a prototypical lattice model of population dynamics, specifically the Contact Process. 2) The numerical analysis is performed with appropriate methods, which are compared with each others, and the results are presented rather clearly.
Weaknesses
Report
The authors conduct a numerical analysis of this model using both Monte Carlo simulations and Strong Disorder Renormalization Group (SDRG) methods, in one and two dimensions. They compute the stationary density (order parameter of the transition) and the local persistence (the probability that a given site is never activated). In the 1D case, their analysis corroborates existing predictions based on SDRG arguments, confirming the infinite-order vanishing of the stationary density at the extinction transition, (Eq. 2), and the mixed-order transition exhibited by local persistence, (Eq. 4). The numerical results further suggest that these behaviors extend to the 2D case.
ASSESSMENT. This work provides a consistent numerical verification of theoretical predictions in 1D and extends these insights to 2D, contributing to the understanding of the effects of quenched disorder and long-range dispersal on the Contact Process. The results are quite clearly presented. Although most of the results verified in this paper have been theoretically derived by the authors in previous works, this paper provides a meaningful and due numerical validation. I therefore recommend publication in SciPost Physics Core.
Requested changes
I think that the presentation of certain results could be slightly improved, to make the manuscript more self-consistent. My suggestions in this direction are reported below.
(i) The Authors stress that the infinite-order behavior in Eq. (2) is a consequence of the combined presence of quenched randomness in the local rates and long-range dispersal, while the two ingredients taken separately lead to an algebraic vanishing of the density at the extinction transition. Do the Authors have some insight on the mechanism behind this?
(ii) The simplified, analytically tractable SDRG analysis of Ref. [31] is mentioned in several points of the manuscript: I would suggest to briefly describe what are the assumptions behind that simplified treatment, in order to make the manuscript self-contained. I would do the same regarding the numerical algorithm to perform SDRG mentioned on p.3 and introduced in Ref. [39].
(iii) What justifies the choice of c, the dilution parameter, used in the Monte Carlo numerics?
(iv) Could the Authors comment on what are the expectations in the case of truly long-range dispersal, i.e. when the exponent alpha is smaller than the dimensionality d?
(v) In the discussion, the Authors mention the model with an environmental gradient as a perspective, and refer to comparison with satellite image data; I would add some reference or clarify this context: which type of data should one think of, which one hopes to describe in terms of a contact process? How should the environmental gradient be implemented?
Recommendation
Ask for minor revision