Memory-Delay Stability Switching and Ecological Thresholds in a Dimensionally Balanced Fractional-Order Predator-Prey Model

Sintesi Tecnica: Commutazione della Stabilità tra Memoria-Ritardo e Soglie Ecologiche in un Modello Predatore-Preda a Ordine Frazionario Bilanciato Dimensionalmente

1. Problematica e Motivazione

Il documento affronta un difetto metodologico critico in molti modelli ecologici a ordine frazionario esistenti: l'incoerenza dimensionale. Nei modelli classici a ordine intero, la variazione $dx/dt$ ha dimensioni di $[Massa][Tempo]^{-1}$. Tuttavia, la derivata frazionaria di Caputo di ordine $0 < \alpha < 1$, indicata con $D_t^\alpha$, possiede intrinsecamente dimensioni di $[Massa][Tempo]^{-\alpha}$. Sostituire direttamente la derivata a ordine intero con una frazionaria senza regolare i parametri del sistema provoca un disallineamento dimensionale, rendendo il modello privo di significato fisico.

Inoltre, gli autori sostengono che i sistemi ecologici siano governati da molteplici meccanismi interagenti spesso studiati isolatamente:

• Effetti di memoria: Le popolazioni biologiche rispondono alla storia passata, non solo agli stati attuali.
• Ritardi temporali: La crescita dei predatori spesso ritarda rispetto al consumo delle prede a causa della gestazione o della digestione.
• Fattori comportamentali: La paura della predazione e il rifugio delle prede alterano significativamente la dinamica delle popolazioni.
• Gestione: Il prelievo non lineare (ad esempio, la pesca) è un meccanio di controllo realistico.

Lo studio mira a formulare un modello predatore-preda a ordine frazionario bilanciato dimensionalmente che integri questi fattori (paura, rifugio, risposta funzionale di tipo Holling III, prelievo non lineare e ritardo temporale) garantendo al contempo la coerenza fisica.

2. Metodologia

Formulazione del Modello:
Gli autori costruiscono un sistema bilanciato dimensionalmente ricostruendo i parametri di tasso dipendenti dal tempo. Invece di sostituire semplicemente $r$ con $r^\alpha$, definiscono coefficienti di tasso frazionari (ad esempio, $r_\alpha = r^\alpha$, $d_\alpha = d^\alpha$, $q_\alpha = q^\alpha$) per garantire che ogni termine nelle equazioni differenziali abbia la dimensione $[Massa][Tempo]^{-\alpha}$. Particolare attenzione è rivolta ai parametri a dimensione mista (come il coefficiente di auto-limitazione $\sigma$), ridefinendo essi come $\sigma_\alpha = \sigma d^{\alpha-1}$ per mantenere l'omogeneità.

Il sistema risultante (Equazione 8) è:
$$\begin{aligned} D_t^\alpha X(t) &= \frac{r_\alpha X(t)}{1 + \kappa Y(t)} \left(1 - \frac{X(t)}{K}\right) - H_\alpha(X(t))Y(t) - \frac{q_\alpha X^2(t)}{c + X(t)} \\ D_t^\alpha Y(t) &= e H_\alpha(X(t-\tau))Y(t) - d_\alpha Y(t) - \sigma_\alpha Y^2(t) \end{aligned}$$
dove $H_\alpha$ è la risposta funzionale di tipo Holling III bilanciata dimensionalmente, $\tau$ è il ritardo temporale e $\alpha$ è l'ordine frazionario.

Framework Analitico:

• Ben determinazione (Well-posedness): Gli autori dimostrano l'esistenza locale, l'unicità, la positività e l'ultimate boundednes (limitatezza ultima) delle soluzioni utilizzando la formulazione integrale di Volterra e i principi di confronto frazionario.
• Analisi dell'Equilibrio: Derivano l'equilibrio della sola preda e l'equilibrio di coesistenza. Viene definito un numero di invasione del predatore ($R_0$) per determinare la soglia di persistenza rispetto all'estinzione del predatore.
• Analisi della Stabilità:
  • Senza ritardo ($\tau=0$): La stabilità è analizzata utilizzando il criterio di Matignon, che richiede che gli autovalori $\mu$ dello Jacobiano soddisfino $|\arg(\mu)| > \alpha\pi/2$.
  • Con ritardo ($\tau>0$): Viene derivata l'equazione caratteristica. Gli autori risolvono per le frequenze critiche e i valori di ritardo critico ($\tau_j$) in cui le radici attraversano l'asse immaginario. Sono stabiliti i termini di trasversalità per determinare la direzione della commutazione della stabilità.

Simulazione Numerica:
Le simulazioni sono eseguite utilizzando un metodo predittore-correttore per equazioni differenziali differenziali a ritardo di Caputo. Gli autori ricostruiscono esplicitamente i coefficienti bilanciati dimensionalmente per ogni simulazione anziché utilizzare una ingenua sostituzione di parametri.

3. Contributi Chiave

1. Coerenza Dimensionale: L'articolo fornisce un framework rigoroso per costruire modelli ecologici frazionari che siano fisicamente coerenti dal punto di vista dimensionale, correggendo la pratica comune della sostituzione diretta della derivata.
2. Classificazione Completa della Stabilità: Lo studio classifica e verifica numericamente sei distinti scenari di commutazione della stabilità indotti dal ritardo per l'equilibrio di coesistenza:
   • Rimane stabile
   • Da stabile a instabile
   • Stabile-instabile-stabile (Bolla di stabilità)
   • Rimane instabile
   • Da instabile a stabile
   • Instabile-stabile-instabile (Bolla di instabilità)
3. Interazione Memoria-Ritardo: L'analisi dimostra come l'ordine frazionario $\alpha$ (memoria) interagisca con il ritardo temporale $\tau$. Il ritardo critico $\tau_0(\alpha)$ si dimostra essere una funzione di $\alpha$, indicando che gli effetti di memoria possono stabilizzare o destabilizzare il sistema a seconda dell'entità del ritardo.
4. Soglie di Gestione Ecologica: Il documento collega la stabilità matematica alla gestione ecologica analizzando l'impatto del prelievo non lineare ($q$) sul numero di invasione del predatore ($R_0$).

4. Risultati

• Esistenza e Persistenza: Il modello garantisce soluzioni positive e limitate. Viene stabilita una condizione di soglia $R_0(q, \alpha) > 1$ per l'invasione del predatore. Gli autori dimostrano che l'aumento dell'intensità del prelievo $q$ riduce strettamente $R_0$, portando potenzialmente all'estinzione del predatore.
• Commutazione della Stabilità (Stability Switching):
  • Per il sistema senza ritardo, la stabilità dipende dalla traccia e dal determinante dello Jacobiano e dall'ordine frazionario $\alpha$.
  • Per il sistema con ritardo, l'introduzione di $\tau$ può indurre biforcazioni di Hopf. Lo studio conferma che la stabilità non è monotona rispetto al ritardo; il sistema può transitare da stabile a instabile e nuovamente a stabile (o viceversa) all'aumentare di $\tau$.
• Stabilizzazione Indotta dalla Memoria: I risultati numerici mostrano che per un ritardo fisso che causa oscillazioni in un modello a ordine intero, la riduzione dell'ordine frazionario $\alpha$ (aumento della memoria) può sopprimere le oscillazioni e stabilizzare l'equilibrio di coesistenza.
• Sensibilità: L'analisi di sensibilità indica che la dinamica del sistema è altamente reattiva alle variazioni nei tassi di prelievo e nelle densità iniziali di popolazione.

5. Significato e Rivendicazioni

L'articolo sostiene che la sua importanza primaria risieda nel colmare il divario tra coerenza fisica e realismo ecologico nella modellazione frazionaria. Garantendo il bilanciamento dimensionale, il modello evita gli artefatti matematici associati alla frazionazione ingenua.

Gli autori sottolineano che il loro lavoro fornisce un quadro unificato per studiare come memoria, ritardo, paura, rifugio e prelievo influenzino congiuntamente la persistenza della popolazione. Nello specifico, affermano di aver:

• Dimostrato che la memoria frazionaria può agire come un meccanismo di stabilizzazione contro le oscillazioni indotte dal ritardo.
• Fornito una classificazione completa degli scenari di commutazione della stabilità, mostrando che il ritardo può avere un doppio ruolo (destabilizzare gli stati stabili e stabilizzare quelli instabili).
• Offerto spunti gestionali, mostrando che l'intensità del prelievo è un parametro critico che può essere calibrato per controllare l'invasione del predatore e la coesistenza.

Lo studio conclude che i lavori futuri potrebbero estendere questo framework per includere la diffusione spaziale, perturbazioni stocastiche o strategie di prelievo ottimali, ma il contributo attuale è focalizzato sulla rigorosa formulazione e l'analisi della stabilità del sistema ritardato bilanciato dimensionalmente.