---
title: Multiplex Bi-Virus Reaction-Diffusion (MBRD)
url: https://www.emergentmind.com/topics/multiplex-bi-virus-reaction-diffusion-framework-mbrd
type: topic
---

# Multiplex Bi-Virus Reaction-Diffusion (MBRD)

The Multiplex Bi-Virus Reaction-Diffusion framework (MBRD) is a class of multiplex metapopulation network models for the spatio-temporal evolution of two interacting pathogens, with diffusion and cross-diffusion acting on layer-specific graphs and nonlinear reaction terms governing transmission, recovery, death, super-infection, or co-infection. In the 2025 formulations, MBRD appears in two concrete closures: the super-infection model (MBRD-SI) and the co-infection model (MBRD-CI), both of which were introduced to capture epidemic pattern formation, Turing and Turing-Hopf instabilities, hotspot growth, coexistence, and strain dominance on multiplex networks [2508.15740].

## 1. Conceptual definition and modeling scope

MBRD treats a population as a collection of size-normalized patches or regions indexed by \(i=1,\dots,N\). Each patch carries susceptible hosts and two pathogen-associated compartments. In the super-infection closure, the state is \((S_i,I_i,J_i)\), where \(I_i\) and \(J_i\) denote mono-infection by virus 1 and virus 2. In the co-infection closure, a fourth compartment \(C_i\) tracks hosts carrying both pathogens simultaneously. The framework combines compartmental epidemic kinetics with multiplex transport: susceptibles and infected subpopulations move on distinct network layers, and the motion of one compartment can depend on gradients of another through cross-diffusion.

This construction is intended for settings in which pathogen interaction is not reducible to a single effective strain. In MBRD-SI, virus 2 can displace virus 1 through a super-infection coefficient \(\sigma\). In MBRD-CI, mono-infected and co-infected classes interact through the parameters \(\beta_{10}\), \(\beta_{02}\), \(\beta_{12}\), and \(\alpha_{12}\), allowing direct co-transmission and co-infection mortality. The reported applications extend beyond epidemiology to information propagation, malware diffusion, and urban transportation networks, where the same reaction-diffusion structure is reinterpreted in non-biological terms [2508.15740].

A distinct but related line of work in molecular communication developed a two-species nonlinear reaction-diffusion PDE system and explicitly outlined how such a system could be generalized to an \(M\)-species multiplexed reaction-diffusion setting. This suggests a methodological antecedent for multi-species reaction-diffusion computation, although the epidemic MBRD framework itself is formulated in the 2025 multiplex bi-virus network papers rather than in the molecular communication setting [1902.11152].

## 2. General multiplex reaction-diffusion structure

In its network-diffusion form, the general MBRD framework is written for each compartment \(X\in\{S,I,J,(C)\}\) as
\[
\frac{dX_i}{dt}
=
\mathrm{Reaction}_i(X_i,\mathrm{others}_i)
+
\sum_{Y\in \mathrm{layers}} d_{X,Y}\sum_{j=1}^N L_{ij}^{(Y)}Y_j.
\]
Here \(L^{(Y)}\) is the graph Laplacian for the layer governing movement of compartment \(Y\), and \(d_{X,Y}\) is the cross-diffusion coefficient by which compartment \(X\) moves up or down the gradient of \(Y\). The same framework is expressed in continuum form by replacing graph Laplacians with diffusion and cross-diffusion operators such as \(\nabla^2\) and \(\nabla\cdot(S\nabla I_j)\) [2509.03374].

The susceptible dynamics in both closures include an Allee-logistic growth term
\[
r\,S_i\Bigl(1-\frac{S_i}{K}\Bigr)\Bigl(\frac{S_i}{A}-1\Bigr),
\]
with per-capita growth rate \(r\), carrying capacity \(K\), and Allee threshold \(A\). Infection terms are normalized by local total population expressions such as \(S_i+I_i+J_i\) or \(S_i+I_i+J_i-C_i\), while recovery, disease-induced death, and background removal enter through \(\gamma\)-, \(\alpha\)-, and \(\mu\)-terms. Boundary conditions are implemented by graph Laplacians with zero-row sums, which enforce the no-flux condition in the discrete network formulation. Initial data are typically chosen as the disease-endemic steady state plus small noise.

The multiplex structure is not merely a notational refinement. The framework assigns separate layers to \(S\), \(I\), and \(J\), so spatial transport is intrinsically heterogeneous across epidemiological classes. Cross-diffusion coefficients \(d_{12}\) and \(d_{13}\) in the super-infection notation, or \(d_{SI}\) and \(d_{SJ}\) in the co-infection notation, allow susceptible motion to be biased by infected gradients. Negative cross-diffusion is especially important in the reported instability mechanisms because it can draw susceptibles toward infected regions and thereby amplify spatial inhomogeneities.

## 3. MBRD-SI: super-infection closure

The MBRD-SI model tracks susceptible density \(S_i(t)\), hosts infected by strain 1 \(I_i(t)\), and hosts infected by strain 2 \(J_i(t)\). Its network ODE system is
\[
\begin{aligned}
\frac{dS_i}{dt}
&=
r\,S_i\Bigl(1-\frac{S_i}{K}\Bigr)\Bigl(\frac{S_i}{A}-1\Bigr)
-\frac{(\beta_1 I_i+\beta_2 J_i)\,S_i}{S_i+I_i+J_i}
+\gamma_1 I_i+\gamma_2 J_i
-\mu\,S_i \\
&\quad
+d_{11}\sum_{j}L_{ij}^{(S)}\,S_j
+d_{12}\sum_{j}L_{ij}^{(I)}\,I_j
+d_{13}\sum_{j}L_{ij}^{(J)}\,J_j,\\[4pt]
\frac{dI_i}{dt}
&=
I_i\biggl(\frac{\beta_1\,S_i}{S_i+I_i+J_i}
-\mu-\alpha_1-\gamma_1
-\frac{\sigma\,\beta_2\,J_i}{S_i+I_i+J_i}\biggr)
+d_{22}\sum_{j}L_{ij}^{(I)}\,I_j,\\[4pt]
\frac{dJ_i}{dt}
&=
J_i\biggl(\frac{\beta_2\,S_i}{S_i+I_i+J_i}
-\mu-\alpha_2-\gamma_2
+\frac{\sigma\,\beta_2\,I_i}{S_i+I_i+J_i}\biggr)
+d_{33}\sum_{j}L_{ij}^{(J)}\,J_j.
\end{aligned}
\]

The biological interpretation is asymmetric. Virus 2 gains an advantage through the super-infection term proportional to \(\sigma\beta_2 I_iJ_i/(S_i+I_i+J_i)\), while virus 1 is correspondingly penalized. In the parameterization adopted in the papers, \(\sigma>1\) represents how much more easily virus 2 infects hosts already infected by virus 1. This asymmetry produces a coexistence-versus-exclusion structure that is central to the model’s long-time behavior. Specifically, large \(\sigma\) leads to strain-2 dominance: \(I_i\to 0\) and \(J_i\to\) endemic in each hotspot [2509.03374].

The continuum version replaces the graph Laplacian terms by \(D_S\nabla^2S\), \(D_{I_1}\nabla^2I_1\), \(D_{I_2}\nabla^2I_2\), and cross-diffusion operators \(D_{S,I_1}\nabla\cdot(S\nabla I_1)\), \(D_{S,I_2}\nabla\cdot(S\nabla I_2)\). This allows the same model class to be studied either as a PDE system or as a multiplex network dynamical system. In the published analysis, both views are used: the continuum description supports dispersion-relation calculations, while the network version supports numerical experiments on lattice, Watts-Strogatz, and Barabási-Albert layers.

## 4. MBRD-CI: co-infection closure

The MBRD-CI model augments the system with a co-infected class \(C_i(t)\). Its reaction terms are organized around the transmission channels \(\beta_1 I\), \(\beta_2 J\), \(\beta_{10}C\), \(\beta_{02}C\), and \(\beta_{12}C\), where \(\beta_{10}\) and \(\beta_{02}\) are baseline co-transmission modifiers, \(\beta_{12}\) is the direct co-transmission rate, and \(\alpha_{12}\) is the disease-induced death rate of the co-infected class. The network system is
\[
\begin{aligned}
\frac{dS_i}{dt}
=&\;r\,S_i\Bigl(1-\frac{S_i}{K}\Bigr)\Bigl(\frac{S_i}{A}-1\Bigr)
-\frac{\beta_1I_i+\beta_2J_i}{S_i+I_i+J_i-C_i}\,S_i \\
&\;-\frac{(\beta_{10}+\beta_{02}+\beta_{12}-\beta_1-\beta_2)\,C_i}{S_i+I_i+J_i-C_i}\,S_i
+\gamma_1I_i+\gamma_2J_i
-(\gamma_1+\gamma_2)\,C_i
-\mu\,S_i
+\sum_{Y=S,I,J}d_{1,Y}\sum_jL_{ij}^{(Y)}Y_j,\\[4pt]
\frac{dI_i}{dt}
=&\;\big[\beta_1I_i+(\beta_{10}+\beta_{12}-\beta_1)C_i\big]
\frac{S_i+J_i-C_i}{S_i+I_i+J_i-C_i}
-\gamma_1I_i-\alpha_1(I_i-C_i)-\alpha_{12}C_i-\mu\,I_i
+d_{22}\sum_jL_{ij}^{(I)}I_j,\\[4pt]
\frac{dJ_i}{dt}
=&\;\big[\beta_2J_i+(\beta_{02}+\beta_{12}-\beta_2)C_i\big]
\frac{S_i+I_i-C_i}{S_i+I_i+J_i-C_i}
-\gamma_2J_i-\alpha_2(J_i-C_i)-\alpha_{12}C_i-\mu\,J_i
+d_{33}\sum_jL_{ij}^{(J)}J_j,\\[4pt]
\frac{dC_i}{dt}
=&\;\frac{\beta_{12}C_i\,S_i}{S_i+I_i+J_i-C_i}
+\frac{[\beta_2J_i+(\beta_{02}+\beta_{12}-\beta_2)C_i](I_i-C_i)}{S_i+I_i+J_i-C_i} \\
&\;+\frac{[\beta_1I_i+(\beta_{10}+\beta_{12}-\beta_1)C_i](J_i-C_i)}{S_i+I_i+J_i-C_i}
-(\gamma_1+\gamma_2+\alpha_{12})\,C_i
-\mu\,C_i .
\end{aligned}
\]

The co-infection closure changes both the state space and the interpretation of persistence. Moderate \(\beta_{12}\) can stabilize coexistence of \(I\), \(J\), and \(C\), whereas very large \(\beta_{12}\) or \(\alpha_{12}\) drive collapse of mono-infection clusters in favor of homogenized co-infection or extinction. In the linearized analysis, MBRD-CI yields a four-morphogen problem and a quartic characteristic polynomial, in contrast to the cubic characteristic equation of MBRD-SI [2509.03374].

A potentially misleading intuition is that the addition of a co-infected class necessarily implies an independently diffusing fourth layer. In the derivation presented for the 2025 model, the network-multiplex form adds layer-specific Laplacians for \(S\), \(I\), and \(J\) diffusion, while the co-infection layer \(C\) carries no independent diffusion because \(C\) is algebraically determined once \(S\), \(I\), and \(J\) are known [2508.15740]. The general MBRD notation can accommodate \(X\in\{S,I,J,(C)\}\), but the concrete MBRD-CI closure introduced in the model papers omits a separate diffusion operator for \(C\).

## 5. Instability analysis, hotspot formation, and multiplex effects

The analytical core of MBRD is a Turing and Turing-Hopf instability analysis about a spatially homogeneous endemic steady state. For MBRD-SI, small perturbations \(\delta x=(\delta S,\delta I,\delta J)^T\) satisfy the eigenvalue problem
\[
(J-k^2D)v=\lambda v,
\qquad
\det[J-k^2D-\lambda I]=0,
\]
which yields a cubic
\[
\lambda^3-b(k^2)\lambda^2+c(k^2)\lambda-d(k^2)=0.
\]
For MBRD-CI, the corresponding characteristic equation is quartic,
\[
\lambda^4-a'\lambda^3+b'\lambda^2-c'\lambda+d'=0.
\]
The reported zero-diffusion stability conditions are \(p_1\equiv \mathrm{tr}(J)<0\), \(p_2\equiv\) sum of minors \(>0\), \(p_3\equiv \det(J)<0\), and \(p_1p_2<p_3\). The Routh-Hurwitz and discriminant tests are then applied as functions of \(k^2\) to separate diffusion-driven static patterning from oscillatory Turing-Hopf onset [2508.15740].

In qualitative terms, hotspots emerge when the usual inhibitor-fast/activator-slow separation is realized in the multiplex setting. Large differences between \(d_{33}\) and \(d_{11}\), together with negative \(d_{12}\) and \(d_{13}\), destabilize the endemic steady state and produce spatial localization. For MBRD-SI, the instability threshold is described roughly by
\[
\frac{d_{33}}{d_{11}}>\Phi(\beta_1,\beta_2,\alpha_1,\alpha_2,\dots),
\]
with \(\Phi\) increasing as \(\sigma\) grows. Analogous threshold expressions hold in the co-infection model in terms of \(\beta_{12}\) and \(\alpha_{12}\) [2509.03374].

The multiplex architecture materially alters these thresholds. If the average degrees of the \(S\), \(I\), and \(J\) layers differ substantially, pattern-forming modes are damped. Balanced degree-regular layers produce the strongest hotspots, whereas mismatched or heterogeneous layers inhibit them. This point addresses a common simplification in multistrain network epidemics: the 2025 MBRD results do not treat layer choice as incidental. Instead, layer regularity, mean degree, and topology enter directly into the onset and persistence of spatial structure.

## 6. Numerical regimes, spread metrics, and extensions

The published numerical experiments make the model classes concrete. For MBRD-SI on lattice, Watts-Strogatz, and Barabási-Albert networks, one example uses
\[
\beta_1=0.3,\;\beta_2=0.15,\;\sigma=3,\;\gamma_1=0.02,\;\gamma_2=0.05,\;\alpha_1=0.02,\;\alpha_2=0.15,
\]
with
\[
d_{11}=0.1,\;d_{12}=d_{13}=-0.2,\;d_{22}=0.01,\;d_{33}=4.8.
\]
Starting near the endemic steady state plus white noise, spatial “dots & stripes” in \(I\) and \(J\) emerge and amplify from \(t\approx 700\to 1900\), and the amplitude \(A(t)\) grows unbounded until local collapse \(S_i\to 0\), \(I_i\to 0\) in hotspots. In the companion instability paper, a closely related super-infection example on a \(100\times 100\) LA12 lattice yields stable spot-like clusters in the \(I_1\)-layer at \(t\approx 1800\) [2509.03374].

For MBRD-CI, a representative experiment uses
\[
\beta_1=0.3,\;\beta_2=0.15,\;\beta_{10}=\beta_{02}=0.1,\;\beta_{12}=0.05,\;
\gamma_1=0.02,\;\gamma_2=0.05,\;\alpha_1=0.02,\;\alpha_2=0.15,\;\alpha_{12}=0.1,
\]
with
\[
d_{11}=0.4,\;d_{12}=d_{13}=-0.2,\;d_{22}=0.01,\;d_{33}=4.8.
\]
The reported outcome is maze-like Turing spots in \(C_i\), together with persistent mono-infection structure in \(I\) and \(J\), amplifying until collapse by \(t\approx 700\). In the earlier model paper, a mixed LA12/LA4 lattice under comparable parameters produces stripe and labyrinth patterns in the \(I_1\)-layer at \(t\approx 550\) [2508.15740].

Parameter sensitivity is explicitly nontrivial. In MBRD-SI, \(A(50)\sim 1848\cdot \sigma^{-11.43}\), so larger \(\sigma\) suppresses pattern reversion. In MBRD-CI, \(A(40)\) versus \(\beta_{12}\) is non-monotonic and peaks at \(\beta_{12}\approx 0.2\), indicating that intermediate co-transmission is the most destabilizing regime. Layer-degree effects are similarly selective: patterns form only when
\[
(\deg_S,\deg_I,\deg_J)\in\{(4,4,4),(12,12,4),(12,12,12),(24,24,24)\},
\]
and higher overall degree gives larger but slower hotspots. Barabási-Albert layers yield faster saturation and weaker dependence on layer degree, suggesting hub-driven super-spread during “holiday” mobility [2509.03374].

The point-source experiments introduce a spread index
\[
\delta_X(t)=\#\{i:X_i(t)>\text{threshold}\}/N,\qquad X\in\{I,J,C\}.
\]
Within MBRD-SI, the \(I\)-spread peak is weakly sensitive to source distance, whereas \(J\)-saturation slows for closer sources. With varying time lag \(t_d\), the reported fits are \(I_{\text{peak}}\sim a_0+a_1\cos(wt_d)+\dots\) and \(J\)-saturation time \(\sim 1.03\,t_d+183\). In MBRD-CI, all three indices saturate later as \(t_d\) grows, and increasing \(\sigma\) lowers \(I_{\text{peak}}\) while advancing its peak time; increasing \(\beta_{12}\) speeds \(C\)-saturation, with a mutual-enhancement regime beyond a co-transmission threshold [2509.03374].

The broader significance claimed for MBRD is not limited to two-pathogen epidemiology. The same formalism is proposed for rumor-versus-rumor dynamics, malware-versus-malware propagation, and multiplex transportation congestion. Negative cross-diffusion is interpreted there as the tendency of naïve nodes to seek influential clusters or, in other application domains, as movement shaped by competing signals or congested routes. Proposed extensions include optimal-control formulations with vaccination or quarantine rates, environmental forcing through time-periodic reaction terms, higher-order networks such as hyperedges and simplicial complexes, and data-driven parameter estimation via physics-informed neural nets from spatial incidence maps. The simulations also suggest that limiting mobility of infected subpopulations through high \(d_{22}\) and \(d_{33}\) is more effective in slowing multi-pathogen spread than blanket reductions in \(d_{11}\), although this remains a numerical implication of the reported experiments rather than a general theorem [2509.03374].

Source: https://www.emergentmind.com/topics/multiplex-bi-virus-reaction-diffusion-framework-mbrd