Diffusion disorder in the contact process

Although power-counting suggests that spatially inhomogeneous diffusion is irrelevant in the contact process, large-scale simulations reveal that such disorder destabilizes the standard directed-percolation critical point by generating effective random-mass disorder, thereby driving the transition into an infinite-randomness universality class.

The Gist

Imagine a vast, quiet city where people (let's call them "infected" or "active" citizens) can spread a rumor to their neighbors. Sometimes, a person stops spreading the rumor and goes home to sleep (they "heal"). Sometimes, they move to a different house to tell the rumor to someone new (they "diffuse").

This is the Contact Process: a simple model scientists use to understand how things like diseases, fires, or ideas spread through a population. Usually, if the rumor spreads fast enough, it takes over the whole city. If it spreads too slowly, it dies out. The exact moment where it switches from dying out to taking over is called a critical point.

For a long time, scientists thought that if the "moving" speed (diffusion) varied randomly from house to house—some people move fast, some don't move at all—it wouldn't change the rules of the game. They thought the city would still switch from "dying out" to "taking over" in the exact same predictable way, just like a clean, uniform city.

But this paper says: "Not so fast!"

Here is the story of what the researchers found, explained simply:

1. The "Wrong" Prediction (The Map vs. The Terrain)

The scientists first looked at the mathematical "map" (a theory called field theory) to predict what would happen. The map said: "If people move around randomly, it's like a gentle breeze. It shouldn't change the storm's behavior. The critical point should stay stable."

They used a method called "power counting" (a way of measuring how much a small change matters). By this math, random moving speeds were supposed to be irrelevant—meaning they wouldn't change the outcome.

2. The Reality Check (The Simulation)

So, the team built a giant digital city and ran millions of computer simulations (Monte Carlo simulations). They created a city where some people had super-fast moving speeds and others were stuck in place (zero speed).

The result? The map was wrong.
The random moving speeds completely broke the old rules. The city didn't switch from dying to thriving in a smooth, predictable way anymore. Instead, it switched in a chaotic, extremely slow, and unpredictable manner.

They discovered that the system had fallen into a new, exotic state called an Infinite-Randomness Fixed Point.

• Analogy: Imagine a clean city where a rumor spreads like a wave. Now, imagine a city with random moving speeds. It's like the rumor is trying to cross a swamp. It gets stuck in deep mud (slow spots) and rushes through dry land (fast spots). The overall progress becomes incredibly slow and depends entirely on the specific, random layout of the mud.

3. The "Magic Trick" (Why it happened)

The researchers asked: "If the math said it shouldn't matter, why did it change everything?"

They developed a clever "effective model" to explain the trick.

• The Setup: Imagine a person who cannot move (speed = 0). They are surrounded by a "pocket" of people who move infinitely fast.
• The Trick: If the stationary person gets "healed" (stops spreading the rumor), a neighbor from the super-fast pocket instantly hops over and re-infects them.
• The Result: The stationary person is effectively immune to healing as long as the fast pocket has anyone active.
• The Translation: The randomness of moving (diffusion) secretly created a randomness in healing. It made some people incredibly hard to cure, while others were easy.

The Metaphor: It's like having a game of "Tag" where the rule "you can't move" actually makes you harder to tag out because your friends are moving so fast they can instantly rescue you. The disorder in movement created a disorder in survival.

4. The Field Theory Explanation (The Deep Dive)

The researchers also went back to the math to prove this. They showed that while the "moving disorder" looked harmless at first glance, when you zoom in and look at the interactions over time (renormalization), it morphs into a "healing disorder."

Think of it like mixing paints. You start with blue paint (moving disorder). You think it will just stay blue. But under the microscope of the math, the blue paint chemically reacts and turns into red paint (healing disorder), which is known to change the game completely.

5. The Big Picture

This paper is important because it teaches us a valuable lesson about complex systems:

• Don't trust the simple math blindly: Just because a variable looks "small" or "irrelevant" in a basic equation doesn't mean it won't cause a chain reaction that changes the whole system.
• Hidden connections: In nature, things that seem unrelated (like how fast you move and how long you survive) are often secretly linked. Changing one can accidentally break the other.

In a nutshell:
The researchers found that if you make the "movement" of particles random, it doesn't just make them move weirdly; it secretly makes their "survival" random too. This turns a predictable, smooth transition into a chaotic, slow-motion struggle, proving that in the world of non-equilibrium physics, random movement is a powerful disruptor.

Technical Summary

Here is a detailed technical summary of the paper "Diffusion disorder in the contact process" by Anfray et al.

1. Problem Statement

The paper investigates the effect of quenched spatial disorder in diffusion rates on the non-equilibrium phase transition of the Contact Process (CP). The CP is a paradigmatic model for absorbing-state phase transitions belonging to the Directed Percolation (DP) universality class.

• The Paradox: Standard field-theoretic power-counting arguments suggest that uncorrelated spatial fluctuations in the diffusion coefficient are irrelevant perturbations at the DP critical point. This is because the scaling dimension of the diffusion disorder term is negative ($[\sigma_D] = -d$).
• The Contradiction: Despite this theoretical prediction, recent studies on other reaction-diffusion systems (like the Diffusive Epidemic Process) suggest that diffusion disorder can significantly alter critical behavior. The authors aim to resolve whether diffusion disorder is truly irrelevant in the CP or if it destabilizes the clean DP critical point.

2. Methodology

The authors employ a multi-pronged approach combining analytical field theory, effective modeling, and large-scale numerical simulations.

A. Field-Theoretic Analysis (Power Counting & Renormalization)

• Power Counting: They start with the dynamic functional of the CP (DP field theory) and introduce a Gaussian random variable for the diffusion coefficient $D(x) = D_0 + \chi_D(x)$. Standard scaling analysis confirms that the disorder term has a negative scaling dimension, predicting irrelevance.
• Renormalization Group (RG): To go beyond simple power counting, they analyze Feynman diagrams. They demonstrate that while one-loop diagrams involving diffusion disorder retain momentum dependence ($k^2$ factors) and do not renormalize the mass term, two-loop diagrams can integrate out these momentum factors. This generates an effective random-mass disorder term (disorder in the control parameter $r$), which is known to be a relevant perturbation.

B. Effective Model Construction

• To intuitively understand the mechanism, the authors construct an effective model where the local diffusion rate $D_i$ takes only two values: $0$(with probability$p$) or$\infty$(with probability$1-p$).
• In regions with $D=\infty$ ("pockets"), particles diffuse instantly, making the density uniform.
• They derive that an active site with $D=0$ is effectively protected from healing as long as an adjacent pocket contains an active particle. This mechanism generates spatially varying effective healing rates ($\tilde{\tau}_i$) for the $D=0$ sites, effectively mapping the diffusion disorder problem onto a CP with quenched disorder in the healing rates.

C. Monte Carlo Simulations

• System: One-dimensional Contact Process with quenched binary diffusion disorder (sites have $D=0$ or $D=10$).
• Algorithm: They use the Dickman algorithm for event-based simulation, tracking infection, healing, and diffusion events.
• Parameters: Large system sizes ($L$ up to $2 \times 10^8$) and long simulation times ($t$up to$3 \times 10^8$) with extensive disorder averaging ($10^5$to$24,600$ realizations).
• Analysis: They analyze scaling behaviors of the active site density ($\rho$), survival probability ($P_s$), and cluster size ($N_s$) to distinguish between conventional power-law scaling (clean DP) and activated scaling (infinite-randomness).

3. Key Contributions & Results

A. Destabilization of the Clean Critical Point

The simulations conclusively show that quenched diffusion disorder is relevant. It destabilizes the clean DP critical point. The transition does not follow the standard power-law scaling of the clean DP class.

B. Emergence of Infinite-Randomness Fixed Point (IRFP)

The critical behavior of the CP with diffusion disorder belongs to the same universality class as the CP with disorder in infection/healing rates (random-mass disorder).

• Scaling: The system exhibits activated scaling (logarithmic time dependence) characteristic of an Infinite-Randomness Fixed Point (IRFP).
• Exponents: The measured critical exponents match the predictions of the Strong Disorder Renormalization Group (SDRG) for random-mass disorder:
  • $\bar{\Theta}/\bar{\delta} \approx 3.24$ (vs. SDRG prediction $\approx 3.236$).
  • $\nu_\perp \psi \approx 0.96$ (vs. SDRG prediction $\approx 1.0$).
  • $\psi = 0.5$.
• Verification: The data fits the activated scaling forms:
  $$\rho(t) \sim [\ln(t/t_0)]^{-\bar{\delta}}, \quad N_s \sim [\ln(t/t_0)]^{\bar{\Theta}}$$

C. Mechanism of Relevance

The paper provides two complementary explanations for why diffusion disorder becomes relevant:

1. Effective Model: Diffusion disorder induces an effective disorder in the local healing rates. Since random healing rates are known to be relevant (leading to IRFP), the diffusion disorder inherits this relevance.
2. Field Theory: Under renormalization, the irrelevant diffusion disorder operator generates a relevant random-mass disorder operator via two-loop Feynman diagrams.

D. Consistency with Harris Criterion

The results align with the non-perturbative Harris criterion. Since tuning the disorder concentration ($p$) significantly shifts the critical infection rate ($\lambda_c$), the clean critical behavior (which violates the Harris inequality $d\nu_\perp > 2$) must change.

4. Significance

• Resolution of Theoretical Contradiction: The paper resolves the conflict between naive power-counting (predicting irrelevance) and physical intuition/numerical evidence (predicting relevance). It demonstrates that "irrelevant" operators in the bare theory can generate relevant operators under renormalization.
• Universality Class Identification: It establishes that diffusion disorder in the CP leads to the same IRFP universality class as infection/healing disorder, unifying the understanding of disorder in absorbing-state transitions.
• Broader Implications: The authors suggest that this mechanism (irrelevant operators generating relevant ones) may apply to other non-equilibrium systems. However, they note that in systems where diffusion plays a dominant role (e.g., Diffusive Epidemic Process), diffusion disorder might lead to distinct critical behaviors not captured by random-mass disorder alone, highlighting a need for further research.

Conclusion

The study demonstrates that quenched spatial disorder in diffusion rates is a relevant perturbation for the Contact Process. Contrary to simple power-counting predictions, it drives the system away from the clean Directed Percolation fixed point toward an Infinite-Randomness Fixed Point, governed by activated scaling and sharing the same universality class as the Contact Process with random-mass disorder. This is achieved through the generation of effective random healing rates and the renormalization-induced generation of random-mass terms.