arxiv: 2605.10644 · v1 · submitted 2026-05-11 · ❄️ cond-mat.stat-mech · q-bio.PE

Recognition: 2 theorem links

· Lean Theorem

Susceptible-Infected-Susceptible Model with Mitigation on Scale-Free Networks

Jo\~ao Gabriel Sim\~oes Delboni , M. O. Hase

Authors on Pith no claims yet

Pith reviewed 2026-05-12 04:38 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech q-bio.PE

keywords SIS modelscale-free networksmitigation factorepidemic spreadingdegree exponentheterogeneous mean-fieldendemic prevalencenonlinear saturation

0 comments

The pith

Incorporating mitigation into the SIS model on scale-free networks inverts the dependence of prevalence on the degree exponent at high transmission rates.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies the spread of infection using the susceptible-infected-susceptible model on scale-free networks with an added mitigation factor that reduces transmission due to behavioral responses or external interventions. This factor introduces a nonlinear saturation effect. The analysis reveals that the probability a random link points to an infected node reaches a maximum at some finite transmission rate, while overall disease prevalence keeps rising with the rate. Most notably, the mitigation causes an inversion: at high transmission rates, networks with larger degree exponents have higher prevalence, opposite to the standard model where more heterogeneous networks with smaller exponents spread disease more. This matters because it shows how mitigation can change which network structures are most vulnerable to epidemics.

Core claim

The authors claim that the mitigation mechanism leads to an inversion in the dependence of epidemic observables on the degree exponent at sufficiently high transmission rates. In the modified model, larger exponents produce higher prevalence and increased infection probability along network links, whereas the standard model shows the opposite trend with smaller exponents yielding higher prevalence.

What carries the argument

A nonlinear mitigation factor inspired by Malthus-Verhulst constraints that reduces the effective transmission rate from infected individuals, incorporated into the heterogeneous mean-field equations for the SIS model on scale-free networks.

If this is right

The probability that a link points to an infected node develops a maximum at finite infection rates.
The overall prevalence remains a monotonically increasing function of the transmission rate.
At high transmission rates, larger degree exponents lead to higher endemic prevalence.
In the modified model, infection probability along links increases with the degree exponent at high rates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Control measures might be more impactful on highly heterogeneous networks if this inversion holds.
The model suggests testing behavioral mitigation effects in real contact networks with varying heterogeneity.
Extensions could explore how different forms of mitigation affect the inversion point.

Load-bearing premise

The heterogeneous mean-field approximation remains valid when the nonlinear mitigation term is added, and the specific form of the mitigation factor accurately represents the combined behavioral and external effects.

What would settle it

Performing stochastic simulations on scale-free networks with varying degree exponents and observing whether the endemic prevalence becomes higher for larger exponents above a certain transmission rate would test the predicted inversion.

Figures

Figures reproduced from arXiv: 2605.10644 by Jo\~ao Gabriel Sim\~oes Delboni, M. O. Hase.

**Figure 1.** Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗

read the original abstract

We investigate infectious disease spreading on scale-free networks using a heterogeneous mean-field approach applied to the susceptible-infected-susceptible model, incorporating a mitigation factor. Individual heterogeneity is incorporated through a power-law distribution, while a mitigation factor accounts for behavioral responses and external effects that effectively reduce transmission from infected individuals. This mechanism, inspired by Malthus-Verhulst-type constraints, introduces a nonlinear saturation effect that encodes self-limiting dynamics in a tractable way. Analytical results are supported by stochastic simulations. We find that the mitigation factor induces a nontrivial behavior in the probability that a link points to an infected node, which develops a maximum at finite infection rates. In contrast, the overall prevalence remains a monotonically increasing function of the transmission rate. Additionally, the mitigation mechanism leads to an inversion in the dependence of epidemic observables on the degree exponent at sufficiently high transmission rates. While in the standard model smaller exponents yield higher endemic prevalence, in the modified model this trend reverses, with larger exponents producing higher prevalence and increased infection probability along network links.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper reports that a mitigation term in SIS on scale-free networks inverts the usual degree-exponent dependence of prevalence at high transmission rates, but the heterogeneous mean-field closure with the added nonlinearity is the part that needs the most scrutiny.

read the letter

The central new observation is that the mitigation factor, which adds a nonlinear reduction to transmission, produces a peak in the probability that a random link points to an infected node at intermediate rates, while overall prevalence still rises monotonically. At high enough transmission the dependence on the degree exponent gamma flips: larger gamma now gives higher prevalence and higher link-infection probability, opposite to the standard SIS case where more heterogeneous networks (smaller gamma) spread more readily. They derive this from heterogeneous mean-field equations and say stochastic simulations on the networks confirm the trends. That is a clean, specific extension of the usual model, and the fact that they tie the new behaviors directly to the mitigation term is useful for anyone who wants a simple way to add self-limiting effects without full individual-based rules. The equations stay tractable, which is the main practical gain here. The simulations are mentioned as support, which is better than pure theory alone. The mitigation form itself is an external modeling choice rather than something derived from the network, but that is standard for this kind of phenomenological addition and does not make the math circular. The soft spot is exactly the one the stress test flags. Once the mitigation couples the degree classes through a global infected density, the mean-field closure can deviate more on scale-free graphs where hubs already strain the approximation, especially near gamma = 3. If the reported inversion only appears cleanly in the closed equations and the simulations are not run at large enough sizes or with enough independent realizations to check the reversal across gamma, then the inversion could be partly an artifact of the closure rather than a robust network effect. The abstract claims agreement, but without the actual simulation protocols or error bars it is hard to judge how strong that agreement is. This is aimed at people who model epidemics on heterogeneous networks and want to include behavioral saturation in a mean-field setting. A reader working on extensions of SIS or on mean-field validity for nonlinear terms would get concrete new behaviors to test. It is not broad enough or rigorous enough to be a landmark, but the specific claims are falsifiable and worth checking, so it deserves a serious referee who can ask for tighter simulation-theory comparisons and perhaps an exact check on small networks. I would send it to review rather than desk reject, with the main request being clearer evidence that the gamma inversion survives beyond the mean-field level.

Referee Report

2 major / 1 minor

Summary. The manuscript applies heterogeneous mean-field theory to the SIS epidemic model on scale-free networks with a power-law degree distribution. It introduces an ad-hoc mitigation factor, inspired by Malthus-Verhulst saturation, that nonlinearly reduces transmission from infected nodes to account for behavioral and external effects. Analytical expressions are derived for endemic prevalence and the probability that a random link points to an infected node; these are compared to stochastic simulations. The key claims are that mitigation produces a maximum in link-infection probability at finite transmission rates (while prevalence remains monotonic), and that at high transmission rates the dependence of both observables on the degree exponent gamma inverts relative to the standard SIS model, with larger gamma yielding higher prevalence and link-infection probability.

Significance. If the inversion result holds under the stated approximations, the work would be significant for showing how nonlinear self-limiting mechanisms can reverse the usual role of degree heterogeneity in epidemic spreading. The combination of tractable HMF analysis with simulation support is a strength, as is the explicit contrast to the standard model. However, the ad-hoc character of the mitigation term and the untested validity of the mean-field closure after its inclusion limit the immediate impact.

major comments (2)

[model and analytical results (as described in abstract)] The heterogeneous mean-field closure is applied after inserting the nonlinear mitigation term that couples all degree classes through a global or average infected density. Standard SIS HMF is already known to have limitations for gamma near 3; the added nonlinearity can shift the regime of validity differently for small versus large gamma, making the reported reversal in prevalence and link-infection probability versus gamma potentially an artifact of the closure rather than a network-level effect. This is load-bearing for the central claim in the abstract.
[mitigation mechanism and simulation support] The functional form of the mitigation factor is introduced as an external modeling choice rather than derived from the network equations. The abstract states that analytical results are supported by stochastic simulations, but without explicit details on parameter regimes, network sizes, or direct comparisons that isolate the effect of the nonlinear term, it is difficult to confirm that the gamma-inversion is robust and not due to post-hoc adjustments.

minor comments (1)

[Abstract] The abstract would benefit from specifying the range of gamma and transmission rates at which the inversion occurs, as well as the precise functional form of the mitigation factor.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us improve the clarity and robustness of our presentation. We address each major comment below, indicating the revisions made to the manuscript.

read point-by-point responses

Referee: The heterogeneous mean-field closure is applied after inserting the nonlinear mitigation term that couples all degree classes through a global or average infected density. Standard SIS HMF is already known to have limitations for gamma near 3; the added nonlinearity can shift the regime of validity differently for small versus large gamma, making the reported reversal in prevalence and link-infection probability versus gamma potentially an artifact of the closure rather than a network-level effect. This is load-bearing for the central claim in the abstract.

Authors: We agree that the standard HMF approximation has known limitations near γ ≈ 3 and that the global nature of the mitigation term warrants careful scrutiny. The mitigation depends on the average infected density, which is already a mean-field quantity, so the closure structure remains formally the same as in the standard SIS-HMF equations. To strengthen the claim, we have added a dedicated paragraph in the revised discussion section analyzing the regime of validity and performed additional stochastic simulations on configuration-model networks (N = 10^4) for γ = 2.5, 3.0, 3.5. These simulations reproduce the γ-inversion at high transmission rates, indicating that the effect is not an artifact of the closure alone. We have also included a brief comparison of HMF predictions versus direct simulations to quantify the approximation error. revision: partial
Referee: The functional form of the mitigation factor is introduced as an external modeling choice rather than derived from the network equations. The abstract states that analytical results are supported by stochastic simulations, but without explicit details on parameter regimes, network sizes, or direct comparisons that isolate the effect of the nonlinear term, it is difficult to confirm that the gamma-inversion is robust and not due to post-hoc adjustments.

Authors: The mitigation factor is introduced phenomenologically, as stated in the model section, to capture behavioral and external saturation effects in a tractable manner inspired by Malthus-Verhulst dynamics; we do not claim a microscopic derivation. We have revised the methods and results sections to supply the requested details: networks of size N = 10^4 generated via the configuration model, γ values 2.5/3.0/3.5, transmission rates β ∈ [0.01, 20] with μ = 1, and averages over 50 independent realizations. New supplementary figures directly compare the full model against the standard SIS case (mitigation coefficient set to zero) to isolate the nonlinear term, confirming both the maximum in link-infection probability and the high-β γ-inversion. These additions demonstrate robustness within the simulated regimes. revision: yes

Circularity Check

0 steps flagged

No significant circularity; mitigation term is explicit external ansatz yielding derived network behaviors

full rationale

The paper sets up the heterogeneous mean-field equations for the SIS model on scale-free networks and adds a mitigation factor (nonlinear saturation term) as an explicit modeling choice inspired by Malthus-Verhulst constraints. The inversion of prevalence dependence on the degree exponent, the maximum in link-infection probability, and monotonicity of overall prevalence are obtained by solving these modified equations. No self-definitional reductions, no fitted parameters renamed as predictions, no load-bearing self-citations, and no ansatz smuggled via prior work appear in the provided derivation outline. Analytical results are cross-checked against stochastic simulations, confirming the chain remains independent of its inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The central results rest on the validity of the heterogeneous mean-field closure with the added nonlinear term and on the chosen functional form of mitigation; these are modeling choices rather than derived quantities.

free parameters (2)

mitigation factor
Parameter controlling the strength of transmission reduction, introduced to encode behavioral and external effects; specific functional dependence not detailed in abstract.
degree exponent gamma
Power-law tail exponent of the degree distribution, varied to study heterogeneity effects.

axioms (2)

domain assumption Heterogeneous mean-field theory provides a closed set of equations for the SIS dynamics on uncorrelated scale-free networks even after adding the mitigation nonlinearity.
Invoked to obtain analytical expressions for link infection probability and prevalence.
ad hoc to paper The mitigation factor takes a form inspired by Malthus-Verhulst saturation that effectively reduces transmission from infected nodes.
Chosen to introduce self-limiting dynamics in a tractable manner.

invented entities (1)

mitigation factor no independent evidence
purpose: To account for behavioral responses and external effects that reduce transmission from infected individuals.
New modeling construct added to the standard SIS equations; no independent empirical calibration or falsifiable prediction outside the model is mentioned.

pith-pipeline@v0.9.0 · 5489 in / 1644 out tokens · 44880 ms · 2026-05-12T04:38:33.788040+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We replace (3) by Θ_mHMF = 1/⟨k⟩ ∑ k P(k) ρ_k (1-ρ_k) ... inspired by Malthus-Verhulst-type constraints, introduces a nonlinear saturation effect
IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the mitigation mechanism leads to an inversion in the dependence of epidemic observables on the degree exponent

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

[1]

C. J. Murray and A. D. Lopez, Mortality by cause for eight regions of the world: Global burden of disease study, The Lancet349, 1269 (1997)

work page 1997
[2]

W. H. O. WHO, Global health estimates (2019), https://www.who.int/data/gho/data/themes/mortality- and-global-health-estimates. Accessed: May 4, 2025

work page 2019
[3]

W. O. Kermack and A. G. Mckendrick, A contribution to the mathematical theory of epidemics, Proceedings of The Royal Society A: Mathematical, Physical and Engi- neering Sciences115, 700 (1927)

work page 1927
[4]

M. J. Keeling and P. Rohani, Modeling infectious diseases in humans and animals (Princeton University Press, 2008)

work page 2008
[5]

M. E. J. Newman, Networks (OUP Oxford, 2018)

work page 2018
[6]

A. L. Barab´ asi,Network Science (Cambridge University Press, 2016)

work page 2016
[7]

Joo and J

J. Joo and J. L. Lebowitz, Behavior of susceptible- infected-susceptible epidemics on heterogeneous net- works with saturation, Physical Review E69, 066105 (2004)

work page 2004
[8]

Dias and M

C. Dias and M. O. Hase, Weakening connections in het- erogeneous mean-field models, Journal of Statistical Me- chanics: Theory and Experiment , 013404 (2021)

work page 2021
[9]

K. M. Kim, C. Dias, and M. O. Hase, Susceptible- infected-susceptible dynamics with mitigation in connec- tion of infected population, Brazilian Journal of Physics 53, 94 (2023)

work page 2023
[10]

Malthus, An Essay on the Principle of Population:, as It Aspects the Future Improvement of Society, with Remarks on the Speculations of Mr

T. Malthus, An Essay on the Principle of Population:, as It Aspects the Future Improvement of Society, with Remarks on the Speculations of Mr. Goodwin, M. Con- dorcet, and Other Writers (J. Johnson, London, 1798)

work page
[11]

P. F. Verhulst, Notice sur la loi que la population suit dans son accroissement, Correspondence Mathematique et Physique10, 113 (1838)

work page
[12]

M. O. Hase and H. L. Casa Grande, Carrying capacity in growing networks, Journal of Statistical Mechanics: Theory and Experiment , 043304 (2016)

work page 2016
[13]

A. L. Barab´ asi and R. Albert, Emergence of scaling in random networks, Science286, 509–512 (1999)

work page 1999
[14]

Pastor-Satorras and A

R. Pastor-Satorras and A. Vespignani, Epidemic dynam- ics and endemic states in complex networks, Physical Re- view Letters86, 3200 (2001)

work page 2001
[15]

D. T. Gillespie, A general method for numerically sim- ulating the stochastic time evolution of coupled chemi- cal reactions, Journal of Computational Physics22, 403 (1976)

work page 1976
[16]

Boguna, R

M. Boguna, R. Pastor-Satorras, and A. Vespignani, Epi- demic spreading in complex networks with degree corre- lations, Lecture Notes in Physics625, 127–147 (2003)

work page 2003
[17]

Liljeros, C

F. Liljeros, C. R. Edling, L. A. N. Amaral, H. E. Stanley, and Y. ˚Aberg, The web of human sexual contacts, Nature 411, 907 (2001)

work page 2001
[18]

Barrat, M

A. Barrat, M. Barthelemy, and A. Vespignani, Dynamical processes on complex networks (Cambridge University Press, 2008)

work page 2008
[19]

I. Z. Kiss, J. C. Miller, and P. L. Simon, Mathematics of epidemics on networks (Springer, 2017)

work page 2017
[20]

Nold, Heterogeneity in disease-transmission modeling, Mathematical Biosciences52, 227 (1980)

A. Nold, Heterogeneity in disease-transmission modeling, Mathematical Biosciences52, 227 (1980)

work page 1980
[21]

H. W. Hethcote, J. A. Yorke, and A. Nold, Gonorrhea modeling: a comparison of control methods, Mathemat- ical Biosciences58, 93 (1982)

work page 1982
[22]

Hethcote, J

H. Hethcote, J. Yorke, and A. Nold, Gonorrhea transmission dynamics and control, Vol. 56 (1984)

work page 1984
[23]

R. M. Anderson and R. M. May, Infectious diseases of humans: dynamics and control (Oxford University Press, 1991)

work page 1991
[24]

Pastor-Satorras and A

R. Pastor-Satorras and A. Vespignani, Epidemic dynam- ics and endemic states in complex networks, Physical Re- view E63, 066117 (2001)

work page 2001
[25]

Eames and M

K. Eames and M. Keeling, Modeling dynamic and net- work heterogeneities in the spread of sexually transmit- ted diseases, Proceedings of the National Academy of Sci- ences of the United States of America99, 13330 (2002)

work page 2002
[26]

Ashcroft and N

N. Ashcroft and N. Mermin, Solid State Physics (Holt, Rinehart and Winston, 1976)

work page 1976
[27]

Castellano and R

C. Castellano and R. Pastor-Satorras, Thresholds for epi- demic spreading in networks, Physical Review Letters 105, 218701 (2010)

work page 2010
[28]

Y. Wang, D. Chakrabarti, C. Wang, and C. Falout- sos, Epidemic spreading in real networks: An eigenvalue viewpoint, in 22nd International SRDS, Vol. 10 (2003) pp. 25 – 34

work page 2003
[29]

W. Cota, A. S. Mata, and S. C. Ferreira, Robustness and fragility of the susceptible-infected-susceptible epi- demic models on complex networks, Physical Review E 98, 012310 (2018)

work page 2018
[30]

Ferguson, D

N. Ferguson, D. Laydon, G. Nedjati-Gilani, et al., Im- pact of non-pharmaceutical interventions (npis) to reduce covid19 mortality and healthcare demand, Imperial Col- lege London, Report 9 (2020)

work page 2020
[31]

Bollob´ as, A probabilistic proof of an asymptotic for- mula for the number of labelled regular graphs, European Journal of Combinatorics1, 311 (1980)

B. Bollob´ as, A probabilistic proof of an asymptotic for- mula for the number of labelled regular graphs, European Journal of Combinatorics1, 311 (1980)

work page 1980
[32]

Catanzaro, M

M. Catanzaro, M. Boguna, and R. Pastor-Satorras, Gen- eration of uncorrelated random scale-free networks, Phys- ical Review E71, 027103 (2005)

work page 2005