Fractal Parameter Boundaries

• Fractal Parameter Boundaries are transition zones in LV competition models that exhibit complex, sensitive dependence on parameters and intricate geometric structures.
• Analytical techniques such as Lyapunov functionals and bifurcation theory elucidate unique equilibrium, coexistence regimes, and spatial segregation.
• Extensions to stochastic and nonlocal diffusion models demonstrate how fractal boundaries influence invasion speeds, extinction thresholds, and pattern formation.

The Lotka-Volterra competition equations describe the population dynamics of two or more species competing for common resources and are foundational in mathematical biology, nonlinear dynamics, and ecological theory. The classically studied deterministic ODE takes the form for populations $n_i(t) \ge 0$, $i=1,\dots,N$,

$$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$

where $r_i$ is the intrinsic growth rate of species $i$ and $b_{ij}$ quantifies the competitive impact of species $j$ on $i$. These equations have been extensively generalized to include spatial structure, stochasticity, time scales, interfacial mixing, anomalous diffusion, and higher-order interaction structures.

1. General Mathematical Formulation and Model Derivation

The canonical formulation for $N$ interacting populations in well-mixed habitats is as above; the system generalizes to resource-based models, spatial PDEs, systems on graphs, and more exotic kinetic and stochastic models. In the mean-field limit, stochastic microscopic models of pairwise competition, where individuals of species $m$ and $i=1,\dots,N$0 interact at rates $i=1,\dots,N$1, yield deterministic equations for the fractions $i=1,\dots,N$2 (with $i=1,\dots,N$3) in the form (Zia, 2010):

$i=1,\dots,N$4

with $i=1,\dots,N$5 antisymmetric. Inclusion of linear birth and death recovers the standard Lotka-Volterra (LV) competition form:

$i=1,\dots,N$6

where $i=1,\dots,N$7 (Zia, 2010).

2. Key Dynamical Regimes and Global Behavior

Existence and Uniqueness of Equilibria

In the deterministic $i=1,\dots,N$8-species LV system, the existence of a unique globally attracting positive equilibrium is guaranteed under strict competition ($i=1,\dots,N$9), symmetry ($$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$0), irreducibility, and suitable growth-rate inequalities. The key tool is a strict Lyapunov functional whose convexity implies global convergence to the unique equilibrium in the positive orthant and excludes nontrivial periodic orbits or chaotic dynamics (Champagnat et al., 2010).

Hypothesis	Requirement	Consequence
Strict competition	$$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$1 strictly increasing, $$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$2	Orbits bounded, no blow-up
Symmetry	$$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$3 for constants $$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$4	Lyapunov function convex, ESS unique
Irreducibility	$$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$5 for all $$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$6	Unique strongly positive equilibrium
Non-extinction	Sufficient condition on $$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$7 vs. resource	Population persists

Phase-Plane Classification (Two Species)

In the planar (two species) case:

• If self-limitation dominates ($$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$8 large relative to $$\dot n_i = n_i \left( r_i - \sum_{j=1}^N b_{ij} n_j \right)$$9), both species may coexist.
• If interspecific competition is too strong, one species excludes the other depending on parameter inequalities.
• For symmetric, irreducible matrices, the unique equilibrium is a sink and all positive solutions converge toward it (Champagnat et al., 2010, Streipert et al., 27 May 2026).

In time-scale models (arbitrary discrete, continuous, or hybrid) the classification persists, but local Jacobians and the dynamic phase-plane can depend on time-dependent "graininess" $r_i$0 (Streipert et al., 27 May 2026). The separatrix between attraction basins is explicitly tracked via root-operator analysis.

3. Spatiotemporal Extensions and Propagation Phenomena

Diffusive and Graph-based PDE Models

On $r_i$1 or finite graphs, competition is modeled by systems such as:

$r_i$2

On graphs, Laplacian operators are replaced by weighted discrete operators with Neumann, Dirichlet, or no boundary conditions (Hu et al., 2022). The equilibrium structure is preserved, with the graph-Laplacian eigenvalues providing threshold quantification for persistence or extinction. The long-time fates (one wins, coexistence, bistability) are determined by explicit inequalities comparing the various competition and growth parameters.

Asymptotic Spreading and Directional Invasion

Recent work in high dimensions characterizes the precise, direction-dependent spreading sets and invasion speeds when initial supports are general measurable sets rather than compact balls (Guo, 25 Feb 2026, Bao et al., 2024). The key is the construction of variational formulas for spreading speeds in each direction based on "bounded” and “unbounded” directions of the support and geometric projection paths. In strong-competition regimes (i.e., $r_i$3 in

$r_i$4

) bistable traveling-wave solutions determine the "glass ceiling" for invasion and coexistence zones in space.

4. Non-classical Competition Structures and Bifurcation Phenomena

Segregation and Pattern Formation Under Strong Competition

In PDE models with strong competition rates ($r_i$5), the equilibrium configurations become spatially segregated: limiting profiles are Lipschitz, and each density occupies disjoint regions (Verzini et al., 2013, Li et al., 2024). When $r_i$6 is used as a bifurcation parameter, one observes the destabilization of spatially homogeneous coexistence and the creation of branches of spatially inhomogeneous solutions, governed by the underlying Laplacian eigenvalues and non-linear free-boundary equations.

Connected and Interface-mediated Populations

In spatial domains segmented by interfaces with partial permeability, as in $r_i$7-species LV-competition on two subdomains sharing a membrane, equilibrium structure can include coexistence regimes and rescue effects even when one region is a "sink" (Álvarez-Caudevilla et al., 2024). The segmentations induce bifurcations that cannot occur in classic homogeneous models.

5. Stochasticity, Jumps, and Anomalous Mobility

Stochastic Competition with Diffusion and Jumps

Generalizations to stochastic differential equations with diffusion and Poissonian jumps are analytically tractable for the LV system (Bao et al., 2011). Explicit SDEs

$r_i$8

admit explicit solutions in 1-dimension. Extinction and persistence criteria are given by sign conditions on compensated growth rates inclusive of noise and jump terms. Sample Lyapunov exponents provide almost-sure asymptotics and exhibit stabilization or extinction mechanisms absent in deterministic settings.

Fractional Diffusion and Nonlocal Interactions

Models with anomalous diffusion (fractional Laplacian) and strong LV competition yield quasi-optimal Hölder regularity and spatial segregation in the singular competition limit (Verzini et al., 2013). When $r_i$9 densities compete, limiting profiles are Lipschitz and satisfy complementary slackness and generalized free-boundary conditions, reflecting the sharp interface dynamics induced by strong nonlocal competition.

6. Oscillatory, Cooperative-Competitive, and Bifurcation Dynamics

In planar LV systems with nonlinear intraspecific effects, e.g. predator-prey models with $i$0 switching sign (competition at high density, cooperation at low), a Hopf bifurcation can generate periodic cycles. The limit cycle amplitude can undergo "blow-up" as a system parameter traverses a finite interval, connecting equilibrium to arbitrarily large oscillatory outbreaks (Bouse et al., 2010). This phenomenon is robust to various functional forms of $i$1 and is confirmed by both Lyapunov function constructions and numerical computations.

Dynamical Phenomenon	Key Mathematical Feature	Biological Interpretation
Hopf bifurcation & blow-up	$i$2, branch of cycles persists and grows unbounded	Extreme cyclical population outbreaks or extinction
Segregation under strong competition	$i$3 limit, piecewise disjoint supports	Spatial exclusion and pattern formation
Multistability with barriers	Multiple segregated equilibria, interface-induced bifurcation	Coexistence enabled by migration/partial permeability

7. Methodologies: Analytical and Constructive Techniques

• Construction of strict Lyapunov and energy functionals for existence and global attractivity (Champagnat et al., 2010, Hu et al., 2022).
• Bifurcation theory, including Crandall–Rabinowitz and Rabinowitz global alternative, for inhomogeneous branches and stability exchange (Li et al., 2024, Álvarez-Caudevilla et al., 2024).
• Augmented, dynamic phase-plane analysis on arbitrary time scales for unified treatment of continuous, discrete, and hybrid systems (Streipert et al., 27 May 2026).
• Variational and monotonicity methods, sub/supersolution barrier constructions for spreading phenomena and threshold classification (Guo, 25 Feb 2026, Hu et al., 2022).
• Rescaling and limiting profile analysis for segregation and free-boundary formation in the strong competition limit (Verzini et al., 2013).

These methods provide a coherent analytical framework bridging deterministic, stochastic, spatial, and hybrid Lotka–Volterra competition theory.