Apéndice A: Demostración del teorema 4.1
Si el número de especies es estrictamente mayor que el número de recursos, \(n>m\), el sistema
admite una solución no trivial \((c_1, \dots , c_n)\N).
Teorema A.1
Supongamos que \(\lim _{Vert \mathbf {x}Vert \rightarrow \infty } R_j(\mathbf {x})=-\infty , j=1,\dots ,m\). Supongamos además que n especies interactúan de acuerdo con (4.3), el número de especies es mayor que el número de recursos \(n>m\) y los recursos dependen de las densidades de las especies de acuerdo con (2.2) de modo que finalmente se agotan. Supongamos además que \i(g_i(\mathbf {x})=1\i) y
para \(j=1,\dots ,m\). Sea \((c_1,\dots ,c_n)\Nuna solución no trivial de (A.1) y supongamos que \N(\sum _{i=1}^n c_i\left( \alpha _i+\frac{sigma _{ii}}{2}right) \ne 0\). Entonces, para cualquier densidad de partida \N(\mathbf {x}(0)\Nen (0,\infty )^n\N) con probabilidad 1
Demostración del Teorema 4.1
Supongamos que \(g_i(\mathbf {x})=1\N), \(i=1,\dots , n\) y \(\sum _jb_{ij}>b_m>0\) para cualquier i y algún \(b_m>0\). Nótese que si \(\sum _j b_{ij}=0\) entonces podemos eliminar \(R_j\) de la ecuación. Supongamos que
Entonces, ya que \(\lim _{Vert \mathbf {x}\Vert \rightarrow \infty } R_j(\mathbf {x})=-\infty \), tenemos cuando \ {\mathbf {x}||\\\\) grande que:
que junto con la linealidad de la parte de difusión implica que la suposición 1.1 del trabajo de HN16 se cumple con \mathbf {c}=(1,\dots , 1)\c). Como resultado, para cualquier punto de partida \(\mathbf {x}(0)\in (0,\infty )^n\) la EDE (4.3) tiene una única solución positiva y por Hening y Nguyen (2018a) (ecuación (5.22)) con probabilidad 1
Sustituyendo eventualmente todos los \(c_i\) por \(-c_i\) podemos suponer que \(\sum _{i=1}^n c_i\left( \alpha _i+\frac{\sigma _{ii}{2}right) >0\). Usando esto en conjunción con (4.3), (A.1) y el Lemma de Itô vemos que
Dejando \N(t\\Nderecha \Ninfty \N) y utilizando que \Nlim _{t\\Nderecha \Ninfty }dfrac{suma _{i=1}^n c_iE_i(t)}{t}=0\} con probabilidad 1
En vista de (A.2) esto implica que con probabilidad 1
(\square \)
Apéndice B: Prueba del Teorema 4.2
Supongamos que dos especies interactúan de acuerdo con
El recurso R depende linealmente de las densidades de las especies
y \N(b_i{overline{R}}>\Nalpha _i, i=1,2\N). Entonces existen \N(\beta _1, \beta _2>0\) tales que las dos especies coexisten.
Demostración
Considera
Si la especie \(x_2\) está ausente la especie \(x_1\) tiene ladinámica unidimensional
Dado que \(b_1{overline{R}}>{alpha _1\}), podemos utilizar a Hening y Nguyen (2018a) para mostrar que el proceso x(t) tiene una medida invariante única en \((0,\infty )\), digamos \(\mu _1\). Además, (Hening y Nguyen 2018a, Lemma 2.1 ) muestra que
o
La tasa de invasión de \(x_2\) con respecto a \(x_1\) puede calcularse mediante (4.2) como
Similarmente, se puede calcular la tasa de invasión de \_1 con respecto a \_2 como
Dado que \N(b_i{overline{R}-\Nalpha _i>0), se puede ver fácilmente que \( \NLambda _{x_}>0), y \(\Lambda _{x_1}>0\) si \( \beta _1>dfrac{b_2a_1(b_1{sobrelínea{R}-alfa _1)}{b_2{sobrelínea{R}-alfa _2}-b_1a_1. \N y \N (\Nbeta _2>dfrac{b_1a_2(b_2\\Nsobrelínea{R}-alfa _2)}{b_1{sobrelínea{R}}-alfa _1}-b_2a_2. \) Si ambas tasas de invasión son positivas obtenemos por Hening y Nguyen (2018a) que las especies coexisten. \(\square \)
Apéndice C: Prueba del Teorema 4.3
Construimos un ejemplo de SDE de dos especies que compiten por un recurso abiótico y coexisten. Observamos que esto ocurre únicamente por el término de variación temporal aleatoria del medio ambiente.
Teorema C.1
Supongamos que la dinámica de las dos especies viene dada por
donde f es una función de Lipschitz continuamente diferenciable que satisface \(\lim _{x\rightarrow -\infty }f(x)=-\infty \), \(\dfrac{df(x)}{dx}>0, \0) para todo \N(x en R) y \N(dfrac{d^2f(x)}{dx^2}<0) para x en algún subintervalo de \N(\Nleft( -\infty ,\frac{{\overline{R}}}{a_1}\right) \). Sea \(a_1, a_2, \sigma _1, \alpha _1,\sigma _1, {\overline{R}} cualesquiera constantes positivas fijas que satisfacen \(f({\overline{R}})>\alpha _1+\dfrac{\sigma _1^2}{2}). Entonces existe un intervalo \((c_0,c_1)\Nsubconjunto (0,\infty )\Ntal que las dos especies coexisten para todo \N(\Nalfa _2\Nen (c_0,c_1)\N).
Prueba
La dinámica de la especie \(x_1\) en ausencia de la especie \(x_2\) viene dada por la EDE unidimensional
Dado que (\lim _{x\rightarrow \infty }f({\overline{R}-a_1x)=-\infty \}), y \f(f(\overline{R})> {alfa _1+dfrac{{sigma _1^2}{2}), esta difusión tiene una única medida de probabilidad invariante \(\mu \) en \((0,\infty )\) cuya densidad es estrictamente positiva en \((0,\infty )\) (ver Borodin y Salminen (2016) o Mao (1997)). Además, observando que \(\lim _{t\arrow \infty }\frac{\ln x(t)}{t}=0\) con probabilidad 1 (usando el Lemma 5.1 de Hening y Nguyen (2018a)) y usando la fórmula de Itô se ve que
Dado que f es una función cóncava y que \(\dfrac{d^2f(x)}{dx^2}>0\ para todo x en algún subintervalo de \(\left( -\infty ,\frac{overline{R}}{a_1}{right) \) debemos tener por la desigualdad de Jensen que
El hecho de que la función f sea estrictamente creciente junto con (C.2) y (C.3) obliga
donde \(f^{-1}\) es la inversa de f- existe porque f es estrictamente creciente. Como resultado, la tasa de invasión de las especies \(x_2\) con respecto a \(x_1\), de la medida de probabilidad invariante \(\mu \), se puede calcular utilizando (4.2) como
Esto implica que \N(\NLambda _{x_2}>0\Nsi y sólo si
La dinámica de la especie \_2 en ausencia de la especie \_1 es
Las soluciones positivas de esta ecuación convergen al punto \(y^*=\dfrac{overline{R}}-\alpha _2}{a_2}) si y sólo si
La tasa de invasión de \(x_1\) con respecto a \(x_2\) será
Nota que como la función f es creciente obtenemos \(\Lambda _{x_1}>0\) si y sólo si
Nótese que \(f^{-1}(\alpha _1+\dfrac{{sigma _1^2}{2}\a la derecha) <{sobrelínea{R}}) ya que por suposición \(f({sobrelínea{R}})>{alfa _1+\dfrac{sigma _1^2}{2}). Como resultado, haciendo uso de las desigualdades (C.5), (C.6) y (C.7) obtenemos que \N(\NLambda _{x_2}>0, \NLambda _{x_1}>0) si y sólo si
Esto implica por el Teorema 3.1 o por Benaim (2018) que las dos especies coexisten. \N-(\N-cuadrado)
Apéndice D: Demostración del teorema 5.1
Supongamos que
Supongamos además que existe un vector \((c_1,\dots ,c_n)\n) que es simultáneamente solución de los sistemas (5.5) para todo \n(u en {1,\dots ,N\}}. Entonces, con probabilidad 1,
Excepto posiblemente para el caso crítico cuando
donde \((\nu _k)_{k en \mathcal {N}}) es la medida de probabilidad invariante de la cadena de Markov (r(t)).
Prueba
Bajo la condición (D.1), existe un \N(M>0\N) tal que el conjunto \N(K_M:={mathbf {x}\Nen \Nmathbb {R}^n: \NVert \Nmathbf {x}\NVert \le M\}) es un atractor global de (5.4). Como resultado, la solución de (5.4) eventualmente entra y nunca sale del conjunto compacto \le (K_M\). En particular, esto demuestra que el proceso \(\mathbf {x}(t)\) está acotado. A continuación, observemos que podemos suponer que
En caso contrario, si \(\sum _{i=1}^n c_i\sum _{k=1}^N\alpha _i(k)\pi _k<0\), podemos sustituir \(c_i\) por \(-c_i\), \(i=1,\dots , n\) y entonces obtener (D.3). Usando (5.1) y el hecho de que \(c_i)’s resuelve (5.5) simultáneamente obtenemos
Poniendo \(t\franja derecha \infty \) y usando la ergodicidad de la cadena de Markov (r(t)) obtenemos que con probabilidad 1
{ {{&
}\sum _{i=1}^n c_i\sum _{k=1}^N\alpha _i(k)\pi _k< 0 \text { a.s.}
Dado que \N(\mathbf {x}(t)\Nestá acotado, esto implica que con probabilidad 1
(\cuadrado{\i})
Apéndice E: Demostración del Teorema 5.2
Según Benaïm y Lobry (2016), Malrieu y Zitt (2017), y Malrieu y Phu (2016) basta con encontrar un ejemplo para el que las tasas de invasión \(\Lambda _{x_1}, \Lambda _{x_}2) sean positivas. Seguiremos a Benaïm y Lobry (2016) para calcular las tasas de invasión de las dos especies. Para \(u=1,2\)\(\mu _u=-\alpha _1(u)+b_1(u){\overline{R}}, \nu _u=-\alpha _2(u)+b_2(u){\overline{R}})\({\overline{a}_u=\dfrac{b_1(u)a_1(u)}{\mu _u}, \overline{b}_u=\dfrac{b_1(u)a_2(u)}{\mu _u}, \overline{c}_u=\dfrac{b_2(u)a_1(u)}{\nu _u}, \overline{d}(u)=\dfrac{b_2(u)a_2(u)}{\nu _u}\), \(p_u=\dfrac{1}{{\overline{a}}_u}, q_u=dfrac{1}{sobrelínea{d}_u}, \gamma _1=dfrac{q_{12}{mu _u}, \gamma _2=dfrac{q_{21}{nu _u}). Si \(p_1\ne p_2), supongamos sin pérdida de generalidad que \(p_1<p_2\). Definir las funciones
y
Por Benaïm y Lobry (2016) tenemos
La expresión para \(\Lambda _{x_1}} se puede obtener intercambiando \(\mu _i\) y \(\nu _i\), \({\overline{a}_i,{\overline{c}_i)\Npor \N({\overline{d}_i,\overline{b}_i)\Ny \Npor \N(p_i) por \N(q_i).
Para el ejemplo de las figuras 3 y 4 hemos utilizado la ecuación integral de (E.1) junto con el paquete de integración numérica de Mathematica para encontrar \(\Lambda _{x_1}>0\) y \(\Lambda _{x_2}>0\). Esto implica por el Teorema 3.1, Benaïm y Lobry (2016) o por Benaim (2018) que las dos especies coexisten.