-
Long-term behavior of stochastic SIQRS epidemic models
Authors:
Alexandru Hening,
Dang H. Nguyen,
Tran Ta,
Sergiu C. Ungureanu
Abstract:
In this paper we analyze and classify the dynamics of SIQRS epidemiological models with susceptible, infected, quarantined, and recovered classes, where the recovered individuals can become reinfected. We are able to treat general incidence functional responses. Our models are more realistic than what has been studied in the literature since they include two important types of random fluctuations.…
▽ More
In this paper we analyze and classify the dynamics of SIQRS epidemiological models with susceptible, infected, quarantined, and recovered classes, where the recovered individuals can become reinfected. We are able to treat general incidence functional responses. Our models are more realistic than what has been studied in the literature since they include two important types of random fluctuations. The first type is due to small fluctuations of the various model parameters and leads to white noise terms. The second type of noise is due to significant environment regime shifts in the that can happen at random. The environment switches randomly between a finite number of environmental states, each with a possibly different disease dynamic. We prove that the long-term fate of the disease is fully determined by a real-valued threshold $λ$. When $λ< 0$ the disease goes extinct asymptotically at an exponential rate. On the other hand, if $λ> 0$ the disease will persist indefinitely. We end our analysis by looking at some important examples where $λ$ can be computed explicitly, and by showcasing some simulation results that shed light on real-world situations.
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Stability and Fluctuations in Complex Ecological Systems
Authors:
Eric Forgoston,
Sarah Day,
Peter C. de Ruiter,
Arjen Doelman,
Nienke Hartemink,
Alan Hastings,
Lia Hemerik,
Alexandru Hening,
Josef Hofbauer,
Sonia Kefi,
David A. Kessler,
Toni Klauschies,
Christian Kuehn,
Xiaoxiao Li,
John C. Moore,
Elly Morrien,
Anje-Margriet Neutel,
Jelena Pantel,
Sebastian J. Schreiber,
Leah B. Shaw,
Nadav Shnerb,
Eric Siero,
Laura S. Storch,
Michael A. S. Thorne,
Ingrid van de Leemput
, et al. (2 additional authors not shown)
Abstract:
From 08-12 August, 2022, 32 individuals participated in a workshop, Stability and Fluctuations in Complex Ecological Systems, at the Lorentz Center, located in Leiden, The Netherlands. An interdisciplinary dialogue between ecologists, mathematicians, and physicists provided a foundation of important problems to consider over the next 5-10 years. This paper outlines eight areas including (1) improv…
▽ More
From 08-12 August, 2022, 32 individuals participated in a workshop, Stability and Fluctuations in Complex Ecological Systems, at the Lorentz Center, located in Leiden, The Netherlands. An interdisciplinary dialogue between ecologists, mathematicians, and physicists provided a foundation of important problems to consider over the next 5-10 years. This paper outlines eight areas including (1) improving our understanding of the effect of scale, both temporal and spatial, for both deterministic and stochastic problems; (2) clarifying the different terminologies and definitions used in different scientific fields; (3) developing a comprehensive set of data analysis techniques arising from different fields but which can be used together to improve our understanding of existing data sets; (4) having theoreticians/computational scientists collaborate closely with empirical ecologists to determine what new data should be collected; (5) improving our knowledge of how to protect and/or restore ecosystems; (6) incorporating socio-economic effects into models of ecosystems; (7) improving our understanding of the role of deterministic and stochastic fluctuations; (8) studying the current state of biodiversity at the functional level, taxa level and genome level.
△ Less
Submitted 12 June, 2023;
originally announced June 2023.
-
Stochastic nutrient-plankton models
Authors:
Alexandru Hening,
Nguyen Trong Hieu,
Dang Hai Nguyen,
Nhu Ngoc Nguyen
Abstract:
We analyze plankton-nutrient food chain models composed of phytoplankton, herbivorous zooplankton and a limiting nutrient. These models have played a key role in understanding the dynamics of plankton in the oceanic layer. Given the strong environmental and seasonal fluctuations that are present in the oceanic layer, we propose a stochastic model for which we are able to fully classify the longter…
▽ More
We analyze plankton-nutrient food chain models composed of phytoplankton, herbivorous zooplankton and a limiting nutrient. These models have played a key role in understanding the dynamics of plankton in the oceanic layer. Given the strong environmental and seasonal fluctuations that are present in the oceanic layer, we propose a stochastic model for which we are able to fully classify the longterm behavior of the dynamics. In order to achieve this we had to develop new analytical techniques, as the system does not satisfy the regular dissipativity conditions and the analysis is more subtle than in other population dynamics models.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
Population dynamics under demographic and environmental stochasticity
Authors:
Alexandru Hening,
Weiwei Qi,
Zhongwei Shen,
Yingfei Yi
Abstract:
The present paper is devoted to the study of the long term dynamics of diffusion processes modelling a single species that experiences both demographic and environmental stochasticity. In our setting, the long term dynamics of the diffusion process in the absence of demographic stochasticity is determined by the sign of $Λ_0$, the external Lyapunov exponent, as follows: $Λ_0<0$ implies (asymptotic…
▽ More
The present paper is devoted to the study of the long term dynamics of diffusion processes modelling a single species that experiences both demographic and environmental stochasticity. In our setting, the long term dynamics of the diffusion process in the absence of demographic stochasticity is determined by the sign of $Λ_0$, the external Lyapunov exponent, as follows: $Λ_0<0$ implies (asymptotic) extinction and $Λ_0>0$ implies convergence to a unique positive stationary distribution $μ_0$. If the system is of size $\frac{1}{ε^{2}}$ for small $ε>0$ (the intensity of demographic stochasticity), demographic effects will make the extinction time finite almost surely. This suggests that to understand the dynamics one should analyze the quasi-stationary distribution (QSD) $μ_ε$ of the system. The existence and uniqueness of the QSD is well-known under mild assumptions.
We look at what happens when the population size is sent to infinity, i.e., when $ε\to 0$. We show that the external Lyapunov exponent still plays a key role: 1) If $Λ_0<0$, then $μ_ε\to δ_0$, the mean extinction time is of order $|\ln ε|$ and the extinction rate associated with the QSD $μ_ε$ has a lower bound of order $\frac{1}{|\lnε|}$; 2) If $Λ_0>0$, then $μ_ε\to μ_0$, the mean extinction time is polynomial in $\frac{1}{ε^{2}}$ and the extinction rate is polynomial in $ε^{2}$. Furthermore, when $Λ_0>0$ we are able to show that the system exhibits multiscale dynamics: at first the process quickly approaches the QSD $μ_ε$ and then, after spending a polynomially long time there, it relaxes to the extinction state. We give sharp asymptotics in $ε$ for the time spent close to $μ_ε$.
△ Less
Submitted 8 July, 2024; v1 submitted 18 July, 2022;
originally announced July 2022.
-
Stochastic epidemic SIR models with hidden states
Authors:
Nguyen Du,
Alexandru Hening,
Nhu Nguyen,
George Yin
Abstract:
This paper focuses on and analyzes realistic SIR models that take stochasticity into account. The proposed systems are applicable to most incidence rates that are used in the literature including the bilinear incidence rate, the Beddington-DeAngelis incidence rate, and a Holling type II functional response. Given that many diseases can lead to asymptomatic infections, we look at a system of stocha…
▽ More
This paper focuses on and analyzes realistic SIR models that take stochasticity into account. The proposed systems are applicable to most incidence rates that are used in the literature including the bilinear incidence rate, the Beddington-DeAngelis incidence rate, and a Holling type II functional response. Given that many diseases can lead to asymptomatic infections, we look at a system of stochastic differential equations that also includes a class of hidden state individuals, for which the infection status is unknown. We assume that the direct observation of the percentage of hidden state individuals that are infected, $α(t)$, is not given and only a noise-corrupted observation process is available. Using the nonlinear filtering techniques in conjunction with an invasion type analysis (or analysis using Lyapunov exponents from the dynamical system point of view), this paper proves that the long-term behavior of the disease is governed by a threshold $λ\in \mathbb{R}$ that depends on the model parameters. It turns out that if $λ<0$ the number $I(t)$ of infected individuals converges to zero exponentially fast, or the extinction happens. In contrast, if $λ>0$, the infection is endemic and the system is permanent. We showcase our results by applying them in specific illuminating examples. Numerical simulations are also given to illustrate our results.
△ Less
Submitted 17 January, 2022;
originally announced January 2022.
-
Random switching in an ecosystem with two prey and one predator
Authors:
Alexandru Hening,
Dang Nguyen,
Nhu Nguyen,
Harrison Watts
Abstract:
In this paper we study the long term dynamics of two prey species and one predator species. In the deterministic setting, if we assume the interactions are of Lotka-Volterra type (competition or predation), the long term behavior of this system is well known. However, nature is usually not deterministic. All ecosystems experience some type of random environmental fluctuations. We incorporate these…
▽ More
In this paper we study the long term dynamics of two prey species and one predator species. In the deterministic setting, if we assume the interactions are of Lotka-Volterra type (competition or predation), the long term behavior of this system is well known. However, nature is usually not deterministic. All ecosystems experience some type of random environmental fluctuations. We incorporate these into a natural framework as follows. Suppose the environment has two possible states. In each of the two environmental states the dynamics is governed by a system of Lotka-Volterra ODE. The randomness comes from spending an exponential amount of time in each environmental state and then switching to the other one. We show how this random switching can create very interesting phenomena. In some cases the randomness can facilitate the coexistence of the three species even though coexistence is impossible in each of the two environmental states. In other cases, even though there is coexistence in each of the two environmental states, switching can lead to the loss of one or more species. We look into how predators and environmental fluctuations can mediate coexistence among competing species.
△ Less
Submitted 24 November, 2021;
originally announced November 2021.
-
The effects of random and seasonal environmental fluctuations on optimal harvesting and stocking
Authors:
Alexandru Hening,
Ky Q. Tran,
Sergiu C. Ungureanu
Abstract:
We analyze the harvesting and stocking of a population that is affected by random and seasonal environmental fluctuations. The main novelty comes from having three layers of environmental fluctuations. The first layer is due to the environment switching at random times between different environmental states. This is similar to having sudden environmental changes or catastrophes. The second layer i…
▽ More
We analyze the harvesting and stocking of a population that is affected by random and seasonal environmental fluctuations. The main novelty comes from having three layers of environmental fluctuations. The first layer is due to the environment switching at random times between different environmental states. This is similar to having sudden environmental changes or catastrophes. The second layer is due to seasonal variation, where there is a significant change in the dynamics between seasons. Finally, the third layer is due to the constant presence of environmental stochasticity -- between the seasonal or random regime switches, the species is affected by fluctuations which can be modelled by white noise. This framework is more realistic because it can capture both significant random and deterministic environmental shifts as well as small and frequent fluctuations in abiotic factors. Our framework also allows for the price or cost of harvesting to change deterministically and stochastically, something that is more realistic from an economic point of view.
The combined effects of seasonal and random fluctuations make it impossible to find the optimal harvesting-stocking strategy analytically. We get around this roadblock by developing rigorous numerical approximations and proving that they converge to the optimal harvesting-stocking strategy. We apply our methods to multiple population models and explore how prices, or costs, and environmental fluctuations influence the optimal harvesting-stocking strategy. We show that in many situations the optimal way of harvesting and stocking is not of threshold type.
△ Less
Submitted 27 September, 2021; v1 submitted 27 September, 2021;
originally announced September 2021.
-
Optimal ergodic harvesting under ambiguity
Authors:
Asaf Cohen,
Alexandru Hening,
Chuhao Sun
Abstract:
We consider an ergodic harvesting problem with model ambiguity that arises from biology. To account for the ambiguity, the problem is constructed as a stochastic game with two players: the decision-maker (DM) chooses the `best' harvesting policy and an adverse player chooses the `worst' probability measure. The main result is establishing an optimal strategy (also referred to as a control) of the…
▽ More
We consider an ergodic harvesting problem with model ambiguity that arises from biology. To account for the ambiguity, the problem is constructed as a stochastic game with two players: the decision-maker (DM) chooses the `best' harvesting policy and an adverse player chooses the `worst' probability measure. The main result is establishing an optimal strategy (also referred to as a control) of the DM and showing that it is a threshold policy. The optimal threshold and the optimal payoff are obtained by solving a free-boundary problem emerging from the Hamilton--Jacobi--Bellman (HJB) equation. As part of the proof, we fix a gap that appeared in the HJB analysis of [Alvarez and Hening, {\em Stochastic Process. Appl.}, 2019], a paper that analyzed the risk-neutral version of the ergodic harvesting problem. Finally, we study the dependence of the optimal threshold and the optimal payoff on the ambiguity parameter and show that if the ambiguity goes to 0, the problem converges, to the risk-neutral problem.
△ Less
Submitted 21 April, 2021;
originally announced April 2021.
-
Quasi-stationary distributions of multi-dimensional diffusion processes
Authors:
Alexandru Hening,
Weiwei Qi,
Zhongwei Shen,
Yingfei Yi
Abstract:
The present paper is devoted to the investigation of the long term behavior of a class of singular multi-dimensional diffusion processes that get absorbed in finite time with probability one. Our focus is on the analysis of quasi-stationary distributions (QSDs), which describe the long term behavior of the system conditioned on not being absorbed. Under natural Lyapunov conditions, we construct a…
▽ More
The present paper is devoted to the investigation of the long term behavior of a class of singular multi-dimensional diffusion processes that get absorbed in finite time with probability one. Our focus is on the analysis of quasi-stationary distributions (QSDs), which describe the long term behavior of the system conditioned on not being absorbed. Under natural Lyapunov conditions, we construct a QSD and prove the sharp exponential convergence to this QSD for compactly supported initial distributions. Under stronger Lyapunov conditions ensuring that the diffusion process comes down from infinity, we show the uniqueness of a QSD and the exponential convergence to the QSD for all initial distributions. Our results can be seen as the multi-dimensional generalization of Cattiaux et al (Ann. Prob. 2009) as well as the complement to Hening and Nguyen (Ann. Appl. Prob. 2018) which looks at the long term behavior of multi-dimensional diffusions that can only become extinct asymptotically.
The centerpiece of our approach concerns a uniformly elliptic operator that we relate to the generator, or the Fokker-Planck operator, associated to the diffusion process. This operator only has singular coefficients in its zeroth-order terms and can be handled more easily than the generator. For this operator, we establish the discreteness of its spectrum, its principal spectral theory, the stochastic representation of the semigroup generated by it, and the global regularity for the associated parabolic equation. We show how our results can be applied to most ecological models, among which cooperative, competitive, and predator-prey Lotka-Volterra systems.
△ Less
Submitted 10 February, 2021;
originally announced February 2021.
-
A general theory of coexistence and extinction for stochastic ecological communities
Authors:
Alexandru Hening,
Dang H. Nguyen,
Peter Chesson
Abstract:
We analyze a general theory for coexistence and extinction of ecological communities that are influenced by stochastic temporal environmental fluctuations. The results apply to discrete time (stochastic difference equations), continuous time (stochastic differential equations), compact and non-compact state spaces and degenerate or non-degenerate noise. In addition, we can also include in the dyna…
▽ More
We analyze a general theory for coexistence and extinction of ecological communities that are influenced by stochastic temporal environmental fluctuations. The results apply to discrete time (stochastic difference equations), continuous time (stochastic differential equations), compact and non-compact state spaces and degenerate or non-degenerate noise. In addition, we can also include in the dynamics auxiliary variables that model environmental fluctuations, population structure, eco-environmental feedbacks or other internal or external factors.
We are able to significantly generalize the recent discrete time results by Benaim and Schreiber (Journal of Mathematical Biology '19) to non-compact state spaces, and we provide stronger persistence and extinction results. The continuous time results by Hening and Nguyen (Annals of Applied Probability '18) are strengthened to include degenerate noise and auxiliary variables.
Using the general theory, we work out several examples. In discrete time, we classify the dynamics when there are one or two species, and look at the Ricker model, Log-normally distributed offspring models, lottery models, discrete Lotka-Volterra models as well as models of perennial and annual organisms. For the continuous time setting we explore models with a resource variable, stochastic replicator models, and three dimensional Lotka-Volterra models.
△ Less
Submitted 17 February, 2021; v1 submitted 17 July, 2020;
originally announced July 2020.
-
A classification of the dynamics of three-dimensional stochastic ecological systems
Authors:
Alexandru Hening,
Dang H. Nguyen,
Sebastian J. Schreiber
Abstract:
The classification of the long-term behavior of dynamical systems is a fundamental problem in mathematics. For both deterministic and stochastic dynamics specific classes of models verify Palis' conjecture: the long-term behavior is determined by a finite number of stationary distributions. In this paper we consider the classification problem for stochastic models of interacting species. For a lar…
▽ More
The classification of the long-term behavior of dynamical systems is a fundamental problem in mathematics. For both deterministic and stochastic dynamics specific classes of models verify Palis' conjecture: the long-term behavior is determined by a finite number of stationary distributions. In this paper we consider the classification problem for stochastic models of interacting species. For a large class of three-species, stochastic differential equation models, we prove a variant of Palis' conjecture: the long-term statistical behavior is determined by a finite number of stationary distributions and, generically, three general types of behavior are possible: 1) convergence to a unique stationary distribution that supports all species, 2) convergence to one of a finite number of stationary distributions supporting two or fewer species, 3) convergence to convex combinations of single species, stationary distributions due to a rock-paper-scissors type of dynamic. Moreover, we prove that the classification reduces to computing Lyapunov exponents (external Lyapunov exponents) that correspond to the average per-capita growth rate of species when rare. Our results stand in contrast to the deterministic setting where the classification is incomplete even for three-dimensional, competitive Lotka--Volterra systems. For these SDE models, our results also provide a rigorous foundation for ecology's modern coexistence theory (MCT) which assumes the external Lyapunov exponents determine long-term ecological outcomes.
△ Less
Submitted 18 May, 2021; v1 submitted 1 April, 2020;
originally announced April 2020.
-
Stationary distributions of persistent ecological systems
Authors:
Alexandru Hening,
Yao Li
Abstract:
We analyze ecological systems that are influenced by random environmental fluctuations. We first provide general conditions which ensure that the species coexist and the system converges to a unique invariant probability measure (stationary distribution). Since it is usually impossible to characterize this invariant probability measure analytically, we develop a powerful method for numerically app…
▽ More
We analyze ecological systems that are influenced by random environmental fluctuations. We first provide general conditions which ensure that the species coexist and the system converges to a unique invariant probability measure (stationary distribution). Since it is usually impossible to characterize this invariant probability measure analytically, we develop a powerful method for numerically approximating invariant probability measures. This allows us to shed light upon how the various parameters of the ecosystem impact the stationary distribution.
We analyze different types of environmental fluctuations. At first we study ecosystems modeled by stochastic differential equations. In the second setting we look at piecewise deterministic Markov processes. These are processes where one follows a system of differential equations for a random time, after which the environmental state changes, and one follows a different set of differential equations -- this procedure then gets repeated indefinitely. Finally, we look at stochastic differential equations with switching, which take into account both the white noise fluctuations and the random environmental switches.
As applications of our theoretical and numerical analysis, we look at competitive Lotka--Volterra, Beddington-DeAngelis predator-prey, and rock-paper-scissors dynamics. We highlight new biological insights by analyzing the stationary distributions of the ecosystems and by seeing how various types of environmental fluctuations influence the long term fate of populations.
△ Less
Submitted 18 May, 2021; v1 submitted 9 March, 2020;
originally announced March 2020.
-
Coexistence, extinction, and optimal harvesting in discrete-time stochastic population models
Authors:
Alexandru Hening
Abstract:
We analyze the long term behavior of interacting populations which can be controlled through harvesting. The dynamics is assumed to be discrete in time and stochastic due to the effect of environmental fluctuations. We present extinction and coexistence criteria when there are one or two interacting species. We then use these tools in order to see when harvesting leads to extinction or persistence…
▽ More
We analyze the long term behavior of interacting populations which can be controlled through harvesting. The dynamics is assumed to be discrete in time and stochastic due to the effect of environmental fluctuations. We present extinction and coexistence criteria when there are one or two interacting species. We then use these tools in order to see when harvesting leads to extinction or persistence of species, as well as what the optimal harvesting strategies, which maximize the expected long term yield, look like. For single species systems, we show under certain conditions that the optimal harvesting strategy is of bang-bang type: there is a threshold under which there is no harvesting, while everything above this threshold gets harvested.
The second part of the paper is concerned with the analysis of ecosystems that have two interacting species which can be harvested. In particular, we carefully study predator-prey and competitive Ricker models when there are two species. In this setting we show how to find the optimal proportional harvesting strategy. If the system is of predator-prey type the optimal proportional harvesting strategy is, depending on the interaction parameters and the price of predators relative to prey, either to harvest the predator to extinction and maximize the asymptotic yield of the prey or to not harvest the prey and to maximize the asymptotic harvesting yield of the predators. If the system is competitive, in certain instances it is optimal to drive one species extinct and to harvest the other one. In other cases it is best to let the two species coexist and harvest both species while maintaining coexistence.
In the setting of the competitive Ricker model we show that if one competitor is dominant and pushes the other species to extinction, the harvesting of the dominant species can lead to coexistence.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Harvesting and seeding of stochastic populations: analysis and numerical approximation
Authors:
Alexandru Hening,
Ky Tran
Abstract:
It is well known that excessive harvesting or hunting has driven species to extinction both on local and global scales. This leads to one of the fundamental problems of conservation ecology: how should we harvest a population so that economic gain is maximized, while also ensuring that the species is safe from extinction? We study an ecosystem of interacting species that are influenced by random e…
▽ More
It is well known that excessive harvesting or hunting has driven species to extinction both on local and global scales. This leads to one of the fundamental problems of conservation ecology: how should we harvest a population so that economic gain is maximized, while also ensuring that the species is safe from extinction? We study an ecosystem of interacting species that are influenced by random environmental fluctuations. At any point in time, we can either harvest or seed (repopulate) species. Harvesting brings an economic gain while seeding incurs a cost. The problem is to find the optimal harvesting-seeding strategy that maximizes the expected total income from harvesting minus the cost one has to pay for the seeding of various species. We consider what happens when one, or both, of the seeding and harvesting rates are bounded. The focus of this paper is the analysis of these three novel settings: bounded seeding and infinite harvesting, bounded seeding and bounded harvesting, and infinite seeding and bounded harvesting.
We prove analytical results and develop numerical approximation methods. By implementing these approximations, we are able to gain qualitative information about how to best harvest and seed species. We are able to show that in the single species setting there are thresholds $0<L_1<L_2<\infty$ such that: 1) if the population size is `low', so that it lies in $(0, L_1]$, there is seeding using the maximal seeding rate; 2) if the population size `moderate', so that it lies in $(L_1,L_2)$, there is no harvesting or seeding; 3) if the population size is `high', so that it lies in the interval $[L_2, \infty)$, there is harvesting using the maximal harvesting rate. Once we have a system with at least two species, numerical experiments show that constant threshold strategies are not optimal anymore.
△ Less
Submitted 2 November, 2019;
originally announced November 2019.
-
Two Long-Period Cataclysmic Variable Stars: ASASSN-14ho and V1062 Cyg
Authors:
L. Claire Gasque,
Callum A. Hening,
Raphael E. Hviding,
John R. Thorstensen,
Kerry Paterson,
Hannes Breytenbach,
Mokhine Motsoaledi,
Patrick A. Woudt
Abstract:
We report spectroscopy and photometry of the cataclysmic variable stars ASASSN-14ho and V1062 Cyg. Both are dwarf novae with spectra dominated by their secondary stars, which we classify as approxomately K4 and M0.5, respectively. Their orbital periods, determined mostly from the secondary stars' radial velociites, proved to be nearly identical, respectively 350.14 +- 0.15 and 348.25 +- 0.60 min.…
▽ More
We report spectroscopy and photometry of the cataclysmic variable stars ASASSN-14ho and V1062 Cyg. Both are dwarf novae with spectra dominated by their secondary stars, which we classify as approxomately K4 and M0.5, respectively. Their orbital periods, determined mostly from the secondary stars' radial velociites, proved to be nearly identical, respectively 350.14 +- 0.15 and 348.25 +- 0.60 min. The H-alpha emission line in V1062 Cyg displays a relatively sharp emission component that tracks the secondary's motion, which may arise on the irradiated face of the secondary; tihs is not often seen and may indicate an unusually strong flux of ionizing radiation. Both systems exhibit double-peaked orbital modulation consistent with ellipsoidal variation from the changing aspect of the secondary. We model these variations to constrain the orbital inclination i, and estimate approximate component masses based oni and the secondary velocity amplitude K2.
△ Less
Submitted 28 September, 2019;
originally announced September 2019.
-
Dynamical systems with fast switching and slow diffusion: Hyperbolic equilibria and stable limit cycles
Authors:
Nguyen H. Du,
Alexandru Hening,
Dang H. Nguyen,
George Yin
Abstract:
We study the long-term qualitative behavior of randomly perturbed dynamical systems. More specifically, we look at limit cycles of stochastic differential equations (SDE) with Markovian switching, in which the process switches at random times among different systems of SDEs, when the switching is fast and the diffusion (white noise) term is small. The system is modeled by…
▽ More
We study the long-term qualitative behavior of randomly perturbed dynamical systems. More specifically, we look at limit cycles of stochastic differential equations (SDE) with Markovian switching, in which the process switches at random times among different systems of SDEs, when the switching is fast and the diffusion (white noise) term is small. The system is modeled by $$ dX^{ε,δ}(t)=f(X^{ε,δ}(t), α^ε(t))dt+\sqrtδσ(X^{ε,δ}(t), α^ε(t))dW(t) , \ X^ε(0)=x, $$ where $α^ε(t)$ is a finite state space Markov chain with irreducible generator $Q=(q_{ij})$. The relative changing rates of the switching and the diffusion are highlighted by the two small parameters $ε$ and $δ$. We associate to the system the averaged ODE \[ d\bar X(t)=\bar f(\bar X(t))dt, \ X(0)=x, \] where $\bar f(\cdot)=\sum_{i=1}^{m_0}f(\cdot, i)ν_i$ and $(ν_1,\dots,ν_{m_0})$ is the unique invariant probability measure of the Markov chain with generator $Q$. Suppose that for each pair $(ε,δ)$ of parameters, the process has an invariant probability measure $μ^{ε,δ}$, and that the averaged ODE has a limit cycle in which there is an averaged occupation measure $μ^0$ for the averaged equation. We are able to prove that if $\bar f$ has finitely many unstable or hyperbolic fixed points, then $μ^{ε,δ}$ converges weakly to $μ^0$ as $ε\to 0$ and $δ\to 0$. Our results generalize to the setting of state-dependent switching \[ \mathbb{P}\{α^ε(t+Δ)=j~|~α^ε=i, X^{ε,δ}(s),α^ε(s), s\leq t\}=q_{ij}(X^{ε,δ}(t))Δ+o(Δ),~~ i\neq j \] as long as the generator $Q(\cdot)=(q_{ij}(\cdot))$ is bounded, Lipschitz, and irreducible for all $x\in\mathbb{R}^d$. We conclude our analysis by studying a predator-prey model.
△ Less
Submitted 18 May, 2021; v1 submitted 5 February, 2019;
originally announced February 2019.
-
Harvesting of interacting stochastic populations
Authors:
Alexandru Hening,
Ky Quan Tran,
Tien Trong Phan,
George Yin
Abstract:
We analyze the optimal harvesting problem for an ecosystem of species that experience environmental stochasticity. Our work generalizes the current literature significantly by taking into account non-linear interactions between species, state-dependent prices, and species injections. The key generalization is making it possible to not only harvest, but also `seed' individuals into the ecosystem. T…
▽ More
We analyze the optimal harvesting problem for an ecosystem of species that experience environmental stochasticity. Our work generalizes the current literature significantly by taking into account non-linear interactions between species, state-dependent prices, and species injections. The key generalization is making it possible to not only harvest, but also `seed' individuals into the ecosystem. This is motivated by how fisheries and certain endangered species are controlled. The harvesting problem becomes finding the optimal harvesting-seeding strategy that maximizes the expected total income from the harvest minus the lost income from the species injections. Our analysis shows that new phenomena emerge due to the possibility of species injections.
It is well-known that multidimensional harvesting problems are very hard to tackle. We are able to make progress, by characterizing the value function as a viscosity solution of the associated Hamilton-Jacobi-Bellman (HJB) equations. Moreover, we provide a verification theorem, which tells us that if a function has certain properties, then it will be the value function. This allows us to show heuristically, as was shown in Lungu and $Ø$ksendal (Bernoulli '01), that it is almost surely never optimal to harvest or seed from more than one population at a time.
We approximate the continuous-time systems by Markov chains and show that the optimal harvesting-seeding strategies of the Markov chain approximations converge to the correct optimal harvesting strategy. This is used to provide numerical approximations to the optimal harvesting-seeding strategies and is a first step towards a full understanding of the intricacies of how one should harvest and seed interacting species. In particular, we look at three examples: one species modeled by a Verhulst-Pearl diffusion, two competing species and a two-species predator-prey system.
△ Less
Submitted 12 October, 2018;
originally announced October 2018.
-
The competitive exclusion principle in stochastic environments
Authors:
Alexandru Hening,
Dang H. Nguyen
Abstract:
In its simplest form, the competitive exclusion principle states that a number of species competing for a smaller number of resources cannot coexist. However, it has been observed empirically that in some settings it is possible to have coexistence. One example is Hutchinson's `paradox of the plankton'. This is an instance where a large number of phytoplankton species coexist while competing for a…
▽ More
In its simplest form, the competitive exclusion principle states that a number of species competing for a smaller number of resources cannot coexist. However, it has been observed empirically that in some settings it is possible to have coexistence. One example is Hutchinson's `paradox of the plankton'. This is an instance where a large number of phytoplankton species coexist while competing for a very limited number of resources. Both experimental and theoretical studies have shown that temporal fluctuations of the environment can facilitate coexistence for competing species. Hutchinson conjectured that one can get coexistence because nonequilibrium conditions would make it possible for different species to be favored by the environment at different times.
In this paper we show in various settings how a variable (stochastic) environment enables a set of competing species limited by a smaller number of resources or other density dependent factors to coexist. If the environmental fluctuations are modeled by white noise, and the per-capita growth rates of the competitors depend linearly on the resources, we prove that there is competitive exclusion. However, if either the dependence between the growth rates and the resources is not linear or the white noise term is nonlinear we show that coexistence on fewer resources than species is possible. Even more surprisingly, if the temporal environmental variation comes from switching the environment at random times between a finite number of possible states, it is possible for all species to coexist even if the growth rates depend linearly on the resources. We show in an example (a variant of which first appeared in Benaim and Lobry '16) that, contrary to Hutchinson's explanation, one can switch between two environments in which the same species is favored and still get coexistence.
△ Less
Submitted 19 September, 2019; v1 submitted 1 October, 2018;
originally announced October 2018.
-
The Inverse First Passage Time Problem for killed Brownian motion
Authors:
Boris Ettinger,
Alexandru Hening,
Tak Kwong Wong
Abstract:
The classical inverse first passage time problem asks whether, for a Brownian motion $(B_t)_{t\geq 0}$ and a positive random variable $ξ$, there exists a barrier $b:\mathbb{R}_+\to\mathbb{R}$ such that $\mathbb{P}\{B_s>b(s), 0\leq s \leq t\}=\mathbb{P}\{ξ>t\}$, for all $t\geq 0$. We study a variant of the inverse first passage time problem for killed Brownian motion. We show that if $λ>0$ is a kil…
▽ More
The classical inverse first passage time problem asks whether, for a Brownian motion $(B_t)_{t\geq 0}$ and a positive random variable $ξ$, there exists a barrier $b:\mathbb{R}_+\to\mathbb{R}$ such that $\mathbb{P}\{B_s>b(s), 0\leq s \leq t\}=\mathbb{P}\{ξ>t\}$, for all $t\geq 0$. We study a variant of the inverse first passage time problem for killed Brownian motion. We show that if $λ>0$ is a killing rate parameter and $\mathbb{1}_{(-\infty,0]}$ is the indicator of the set $(-\infty,0]$ then, under certain compatibility assumptions, there exists a unique continuous function $b:\mathbb{R}_+\to\mathbb{R}$ such that $\mathbb{E}\left[-λ\int_0^t \mathbb{1}_{(-\infty,0]}(B_s-b(s))\,ds\right] = \mathbb{P}\{ζ>t\}$ holds for all $t\geq 0$. This is a significant improvement of a result of the first two authors (Annals of Applied Probability 24(1):1--33, 2014).
The main difficulty arises because $\mathbb{1}_{(-\infty,0]}$ is discontinuous. We associate a semi-linear parabolic partial differential equation (PDE) coupled with an integral constraint to this version of the inverse first passage time problem. We prove the existence and uniqueness of weak solutions to this constrained PDE system. In addition, we use the recent Feynman-Kac representation results of Glau (Finance and Stochastics 20(4):1021--1059, 2016) to prove that the weak solutions give the correct probabilistic interpretation.
△ Less
Submitted 26 August, 2019; v1 submitted 14 July, 2018;
originally announced July 2018.
-
Optimal sustainable harvesting of populations in random environments
Authors:
Luis H. R. Alvarez E.,
Alexandru Hening
Abstract:
We study the optimal sustainable harvesting of a population that lives in a random environment. The novelty of our setting is that we maximize the asymptotic harvesting yield, both in an expected value and almost sure sense, for a large class of harvesting strategies and unstructured population models. We prove under relatively weak assumptions that there exists a unique optimal harvesting strateg…
▽ More
We study the optimal sustainable harvesting of a population that lives in a random environment. The novelty of our setting is that we maximize the asymptotic harvesting yield, both in an expected value and almost sure sense, for a large class of harvesting strategies and unstructured population models. We prove under relatively weak assumptions that there exists a unique optimal harvesting strategy characterized by an optimal threshold below which the population is maintained at all times by utilizing a local time push-type policy. We also discuss, through Abelian limits, how our results are related to the optimal harvesting strategies when one maximizes the expected cumulative present value of the harvesting yield and establish a simple connection and ordering between the values and optimal boundaries. Finally, we explicitly characterize the optimal harvesting strategies in two different cases, one of which is the celebrated stochastic Verhulst Pearl logistic model of population growth.
△ Less
Submitted 6 July, 2018;
originally announced July 2018.
-
Asymptotic harvesting of populations in random environments
Authors:
Alexandru Hening,
Dang H. Nguyen,
Sergiu C. Ungureanu,
Tak Kwong Wong
Abstract:
We consider the harvesting of a population in a stochastic environment whose dynamics in the absence of harvesting is described by a one dimensional diffusion. Using ergodic optimal control, we find the optimal harvesting strategy which maximizes the asymptotic yield of harvested individuals. To our knowledge, ergodic optimal control has not been used before to study harvesting strategies. However…
▽ More
We consider the harvesting of a population in a stochastic environment whose dynamics in the absence of harvesting is described by a one dimensional diffusion. Using ergodic optimal control, we find the optimal harvesting strategy which maximizes the asymptotic yield of harvested individuals. To our knowledge, ergodic optimal control has not been used before to study harvesting strategies. However, it is a natural framework because the optimal harvesting strategy will never be such that the population is harvested to extinction -- instead the harvested population converges to a unique invariant probability measure.
When the yield function is the identity, we show that the optimal strategy has a bang-bang property: there exists a threshold $x^*>0$ such that whenever the population is under the threshold the harvesting rate must be zero, whereas when the population is above the threshold the harvesting rate must be at the upper limit. We provide upper and lower bounds on the maximal asymptotic yield, and explore via numerical simulations how the harvesting threshold and the maximal asymptotic yield change with the growth rate, maximal harvesting rate, or the competition rate.
We also show that, if the yield function is $C^2$ and strictly concave, then the optimal harvesting strategy is continuous, whereas when the yield function is convex the optimal strategy is of bang-bang type. This shows that one cannot always expect bang-bang type optimal controls.
△ Less
Submitted 30 March, 2018; v1 submitted 3 October, 2017;
originally announced October 2017.
-
On a predator-prey system with random switching that never converges to its equilibrium
Authors:
Alexandru Hening,
Edouard Strickler
Abstract:
We study the dynamics of a predator-prey system in a random environment. The dynamics evolves according to a deterministic Lotka-Volterra system for an exponential random time after which it switches to a different deterministic Lotka-Volterra system. This switching procedure is then repeated. The resulting process is a Piecewise Deterministic Markov Process (PDMP). In the case when the equilibriu…
▽ More
We study the dynamics of a predator-prey system in a random environment. The dynamics evolves according to a deterministic Lotka-Volterra system for an exponential random time after which it switches to a different deterministic Lotka-Volterra system. This switching procedure is then repeated. The resulting process is a Piecewise Deterministic Markov Process (PDMP). In the case when the equilibrium points of the two deterministic Lotka--Volterra systems coincide we show that almost surely the trajectory does not converge to the common deterministic equilibrium. Instead, with probability one, the densities of the prey and the predator oscillate between $0$ and $\infty$. This proves a conjecture of Takeuchi et al (J. Math. Anal. Appl 2006).
The proof of the conjecture is a corollary of a result we prove about linear switched systems. Assume $(Y_t, I_t)$ is a PDMP that evolves according to $\frac{dY_t}{dt}=A_{I_t} Y_t$ where $A_0,A_1$ are $2\times2$ matrices and $I_t$ is a Markov chain on $\{0,1\}$ with transition rates $k_0,k_1>0$. If the matrices $A_0$ and $A_1$ are not proportional and have purely imaginary eigenvalues, then there exists $λ>0$ such that $\lim_{t \to \infty} \frac{\log \| Y_t \|}{t} = λ.$
△ Less
Submitted 26 August, 2019; v1 submitted 3 October, 2017;
originally announced October 2017.
-
Persistence in Stochastic Lotka--Volterra food chains with intraspecific competition
Authors:
Alexandru Hening,
Dang H. Nguyen
Abstract:
This paper is devoted to the analysis of a simple Lotka-Volterra food chain evolving in a stochastic environment. It can be seen as the companion paper of Hening and Nguyen (J. of Math. Biol. `18) where we have characterized the persistence and extinction of such a food chain under the assumption that there is no intraspecific competition among predators. In the current paper we focus on the case…
▽ More
This paper is devoted to the analysis of a simple Lotka-Volterra food chain evolving in a stochastic environment. It can be seen as the companion paper of Hening and Nguyen (J. of Math. Biol. `18) where we have characterized the persistence and extinction of such a food chain under the assumption that there is no intraspecific competition among predators. In the current paper we focus on the case when all the species experience intracompetition. The food chain we analyze consists of one prey and $n-1$ predators. The $j$th predator eats the $j-1$st species and is eaten by the $j+1$st predator; this way each species only interacts with at most two other species - the ones that are immediately above or below it in the trophic chain. We show that one can classify, based on the invasion rates of the predators (which we can determine from the interaction coefficients of the system via an algorithm), which species go extinct and which converge to their unique invariant probability measure. We obtain stronger results than in the case with no intraspecific competition because in this setting we can make use of the general results of Hening and Nguyen (Ann. of Appl. Probab.). Unlike most of the results available in the literature, we provide an in depth analysis for both non-degenerate and degenerate noise.
We exhibit our general results by analysing trophic cascades in a plant--herbivore--predator system and providing persistence/extinction criteria for food chains of length $n\leq 3$.
△ Less
Submitted 2 April, 2018; v1 submitted 24 April, 2017;
originally announced April 2017.
-
Coexistence and extinction for stochastic Kolmogorov systems
Authors:
Alexandru Hening,
Dang H. Nguyen
Abstract:
In recent years there has been a growing interest in the study of the dynamics of stochastic populations. A key question in population biology is to understand the conditions under which populations coexist or go extinct. Theoretical and empirical studies have shown that coexistence can be facilitated or negated by both biotic interactions and environmental fluctuations. We study the dynamics of…
▽ More
In recent years there has been a growing interest in the study of the dynamics of stochastic populations. A key question in population biology is to understand the conditions under which populations coexist or go extinct. Theoretical and empirical studies have shown that coexistence can be facilitated or negated by both biotic interactions and environmental fluctuations. We study the dynamics of $n$ populations that live in a stochastic environment and which can interact nonlinearly (through competition for resources, predator-prey behavior, etc.). Our models are described by $n$-dimensional Kolmogorov systems with white noise (stochastic differential equations - SDE). We give sharp conditions under which the populations converge exponentially fast to their unique stationary distribution as well as conditions under which some populations go extinct exponentially fast.
The analysis is done by a careful study of the properties of the invariant measures of the process that are supported on the boundary of the domain. To our knowledge this is one of the first general results describing the asymptotic behavior of stochastic Kolmogorov systems in non-compact domains.
We are able to fully describe the properties of many of the SDE that appear in the literature. In particular, we extend results on two dimensional Lotka-Volterra models, two dimensional predator-prey models, $n$ dimensional simple food chains, and two predator and one prey models. We also show how one can use our methods to classify the dynamics of any two-dimensional stochastic Kolmogorov system satisfying some mild assumptions.
△ Less
Submitted 23 April, 2017;
originally announced April 2017.
-
Stochastic Lotka-Volterra food chains
Authors:
Alexandru Hening,
Dang H. Nguyen
Abstract:
We study the persistence and extinction of species in a simple food chain that is modelled by a Lotka-Volterra system with environmental stochasticity. There exist sharp results for deterministic Lotka-Volterra systems in the literature but few for their stochastic counterparts. The food chain we analyze consists of one prey and $n-1$ predators. The $j$th predator eats the $j-1$th species and is e…
▽ More
We study the persistence and extinction of species in a simple food chain that is modelled by a Lotka-Volterra system with environmental stochasticity. There exist sharp results for deterministic Lotka-Volterra systems in the literature but few for their stochastic counterparts. The food chain we analyze consists of one prey and $n-1$ predators. The $j$th predator eats the $j-1$th species and is eaten by the $j+1$th predator; this way each species only interacts with at most two other species - the ones that are immediately above or below it in the trophic chain. We show that one can classify, based on an explicit quantity depending on the interaction coefficients of the system, which species go extinct and which converge to their unique invariant probability measure. Our work can be seen as a natural extension of the deterministic results of Gard and Hallam '79 to a stochastic setting.
As one consequence we show that environmental stochasticity makes species more likely to go extinct. However, if the environmental fluctuations are small, persistence in the deterministic setting is preserved in the stochastic system. Our analysis also shows that the addition of a new apex predator makes, as expected, the different species more prone to extinction.
Another novelty of our analysis is the fact that we can describe the behavior the system when the noise is degenerate. This is relevant because of the possibility of strong correlations between the effects of the environment on the different species.
△ Less
Submitted 14 March, 2017;
originally announced March 2017.
-
Markov processes conditioned on their location at large exponential times
Authors:
Steven N. Evans,
Alexandru Hening
Abstract:
Suppose that $(X_t)_{t \ge 0}$ is a one-dimensional Brownian motion with negative drift $-μ$. It is possible to make sense of conditioning this process to be in the state $0$ at an independent exponential random time and if we kill the conditioned process at the exponential time the resulting process is Markov. If we let the rate parameter of the random time go to $0$, then the limit of the killed…
▽ More
Suppose that $(X_t)_{t \ge 0}$ is a one-dimensional Brownian motion with negative drift $-μ$. It is possible to make sense of conditioning this process to be in the state $0$ at an independent exponential random time and if we kill the conditioned process at the exponential time the resulting process is Markov. If we let the rate parameter of the random time go to $0$, then the limit of the killed Markov process evolves like $X$ conditioned to hit $0$, after which time it behaves as $X$ killed at the last time $X$ visits $0$. Equivalently, the limit process has the dynamics of the killed "bang--bang" Brownian motion that evolves like Brownian motion with positive drift $+μ$ when it is negative, like Brownian motion with negative drift $-μ$ when it is positive, and is killed according to the local time spent at $0$.
An extension of this result holds in great generality for Borel right processes conditioned to be in some state $a$ at an exponential random time, at which time they are killed. Our proofs involve understanding the Campbell measures associated with local times, the use of excursion theory, and the development of a suitable analogue of the "bang--bang" construction for general Markov processes.
As examples, we consider the special case when the transient Borel right process is a one-dimensional diffusion. Characterizing the limiting conditioned and killed process via its infinitesimal generator leads to an investigation of the $h$-transforms of transient one-dimensional diffusion processes that goes beyond what is known and is of independent interest.
△ Less
Submitted 12 July, 2016;
originally announced July 2016.
-
Stochastic population growth in spatially heterogeneous environments: The density-dependent case
Authors:
Alexandru Hening,
Dang H. Nguyen,
George Yin
Abstract:
This work is devoted to studying the dynamics of a structured population that is subject to the combined effects of environmental stochasticity, competition for resources, spatio-temporal heterogeneity and dispersal. The population is spread throughout $n$ patches whose population abundances are modelled as the solutions of a system of nonlinear stochastic differential equations living on…
▽ More
This work is devoted to studying the dynamics of a structured population that is subject to the combined effects of environmental stochasticity, competition for resources, spatio-temporal heterogeneity and dispersal. The population is spread throughout $n$ patches whose population abundances are modelled as the solutions of a system of nonlinear stochastic differential equations living on $[0,\infty)^n$.
We prove that $r$, the stochastic growth rate of the total population in the absence of competition, determines the long-term behaviour of the population. The parameter $r$ can be expressed as the Lyapunov exponent of an associated linearized system of stochastic differential equations. Detailed analysis shows that if $r>0$, the population abundances converge polynomially fast to a unique invariant probability measure on $(0,\infty)^n$, while when $r<0$, the population abundances of the patches converge almost surely to $0$ exponentially fast. This generalizes and extends the results of Evans et al (2014 J. Math. Biol.) and proves one of their conjectures.
Compared to recent developments, our model incorporates very general density-dependent growth rates and competition terms. Furthermore, we prove that persistence is robust to small, possibly density dependent, perturbations of the growth rates, dispersal matrix and covariance matrix of the environmental noise. Our work allows the environmental noise driving our system to be degenerate. This is relevant from a biological point of view since, for example, the environments of the different patches can be perfectly correlated. As an example we fully analyze the two-patch case, $n=2$, and show that the stochastic growth rate is a decreasing function of the dispersion rate. In particular, coupling two sink patches can never yield persistence, in contrast to the results from the non-degenerate setting treated by Evans et al.
△ Less
Submitted 17 June, 2017; v1 submitted 6 May, 2016;
originally announced May 2016.
-
Transient one-dimensional diffusions conditioned to converge to a different limit point
Authors:
Alexandru Hening
Abstract:
Let $(X_t)_{t\geq 0}$ be a regular one-dimensional diffusion that models a biological population. If one assumes that the population goes extinct in finite time it is natural to study the $Q$-process associated to $(X_t)_{t\geq 0}$. This is the process one gets by conditioning $(X_t)_{t\geq 0}$ to survive into the indefinite future. The motivation for this paper comes from looking at populations t…
▽ More
Let $(X_t)_{t\geq 0}$ be a regular one-dimensional diffusion that models a biological population. If one assumes that the population goes extinct in finite time it is natural to study the $Q$-process associated to $(X_t)_{t\geq 0}$. This is the process one gets by conditioning $(X_t)_{t\geq 0}$ to survive into the indefinite future. The motivation for this paper comes from looking at populations that are modeled by diffusions which do not go extinct in finite time but which go `extinct asymptotically' as $t\rightarrow \infty$. We look at transient one-dimensional diffusions $(X_t)_{t \geq 0}$ with state space $I=(\ell, \infty)$ such that $X_t\rightarrow \ell$ as $t\rightarrow \infty$, $\mathbb{P}^x$-almost surely for all $x\in I$. We `condition' $(X_t)_{t \geq 0}$ to go to $\infty$ as $t\rightarrow \infty$ and show that the resulting diffusion is the Doob $h$-transform of $(X_t)_{t\geq 0}$ with $h=s$ where $s$ is the scale function of $(X_t)_{t\geq 0}$. Finally, we explore what this conditioning does in two examples.
△ Less
Submitted 1 October, 2015;
originally announced October 2015.
-
The free path in a high velocity random flight process associated to a Lorentz gas in an external field
Authors:
Alexandru Hening,
Douglas Rizzolo,
Eric S. Wayman
Abstract:
We investigate the asymptotic behavior of the free path of a variable density random flight model in an external field as the initial velocity of the particle goes to infinity. The random flight models we study arise naturally as the Boltzmann-Grad limit of a random Lorentz gas in the presence of an external field. By analyzing the time duration of the free path, we obtain exact forms for the asym…
▽ More
We investigate the asymptotic behavior of the free path of a variable density random flight model in an external field as the initial velocity of the particle goes to infinity. The random flight models we study arise naturally as the Boltzmann-Grad limit of a random Lorentz gas in the presence of an external field. By analyzing the time duration of the free path, we obtain exact forms for the asymptotic mean and variance of the free path in terms of the external field and the density of scatterers. As a consequence, we obtain a diffusion approximation for the joint process of the particle observed at reflection times and the amount of time spent in free flight.
△ Less
Submitted 8 March, 2015;
originally announced March 2015.
-
Quasistationary distributions for one-dimensional diffusions with singular boundary points
Authors:
Alexandru Hening,
Martin Kolb
Abstract:
In the present work we characterize the existence of quasistationary distributions for diffusions on $(0,\infty)$ allowing singular behavior at $0$ and $\infty$. If absorption at 0 is certain, we show that there exists a quasistationary distribution as soon as the spectrum of the generator is strictly positive. This complements results of Collet et al. (Ann. Probab. 2009) and Kolb and Steinsaltz (…
▽ More
In the present work we characterize the existence of quasistationary distributions for diffusions on $(0,\infty)$ allowing singular behavior at $0$ and $\infty$. If absorption at 0 is certain, we show that there exists a quasistationary distribution as soon as the spectrum of the generator is strictly positive. This complements results of Collet et al. (Ann. Probab. 2009) and Kolb and Steinsaltz (Ann. Probab. 2012) for $0$ being a regular boundary point and extends results by Collet et al. (Ann. Probab. 2009) on singular diffusions.
△ Less
Submitted 30 March, 2018; v1 submitted 8 September, 2014;
originally announced September 2014.
-
When do skew-products exist?
Authors:
Steven N. Evans,
Alexandru Hening,
Eric S. Wayman
Abstract:
The classical skew-product decomposition of planar Brownian motion represents the process in polar coordinates as an autonomously Markovian radial part and an angular part that is an independent Brownian motion on the unit circle time-changed according to the radial part. Theorem 4 of Liao (2009) gives a broad generalization of this fact to a setting where there is a diffusion on a manifold $X$ wi…
▽ More
The classical skew-product decomposition of planar Brownian motion represents the process in polar coordinates as an autonomously Markovian radial part and an angular part that is an independent Brownian motion on the unit circle time-changed according to the radial part. Theorem 4 of Liao (2009) gives a broad generalization of this fact to a setting where there is a diffusion on a manifold $X$ with a distribution that is equivariant under the smooth action of a Lie group $K$. Under appropriate conditions, there is a decomposition into an autonomously Markovian "radial" part that lives on the space of orbits of $K$ and an "angular" part that is an independent Brownian motion on the homogeneous space $K/M$, where $M$ is the isotropy subgroup of a point of $x$, that is time-changed with a time-change that is adapted to the filtration of the radial part. We present two apparent counterexamples to Theorem 4 of Liao (2009). In the first counterexample the angular part is not a time-change of any Brownian motion on $K/M$, whereas in the second counterexample the angular part is the time-change of a Brownian motion on $K/M$ but this Brownian motion is not independent of the radial part. In both of these examples $K/M$ has dimension $1$. The statement and proof of Theorem 4 from Liao (2009) remain valid when $K/M$ has dimension greater than $1$. Our examples raise the question of what conditions lead to the usual sort of skew-product decomposition when $K/M$ has dimension $1$ and what conditions lead to there being no decomposition at all or one in which the angular part is a time-changed Brownian motion but this Brownian motion is not independent of the radial part.
△ Less
Submitted 22 July, 2015; v1 submitted 27 August, 2014;
originally announced August 2014.
-
Protected polymorphisms and evolutionary stability of patch-selection strategies in stochastic environments
Authors:
Steven N. Evans,
Alexandru Hening,
Sebastian J. Schreiber
Abstract:
We consider a population living in a patchy environment that varies stochastically in space and time. The population is composed of two morphs (that is, individuals of the same species with different genotypes). In terms of survival and reproductive success, the associated phenotypes differ only in their habitat selection strategies. We compute invasion rates corresponding to the rates at which th…
▽ More
We consider a population living in a patchy environment that varies stochastically in space and time. The population is composed of two morphs (that is, individuals of the same species with different genotypes). In terms of survival and reproductive success, the associated phenotypes differ only in their habitat selection strategies. We compute invasion rates corresponding to the rates at which the abundance of an initially rare morph increases in the presence of the other morph established at equilibrium. If both morphs have positive invasion rates when rare, then there is an equilibrium distribution such that the two morphs coexist; that is, there is a protected polymorphism for habitat selection. Alternatively, if one morph has a negative invasion rate when rare, then it is asymptotically displaced by the other morph under all initial conditions where both morphs are present. We refine the characterization of an evolutionary stable strategy for habitat selection from [Schreiber, 2012] in a mathematically rigorous manner. We provide a necessary and sufficient condition for the existence of an ESS that uses all patches and determine when using a single patch is an ESS. We also provide an explicit formula for the ESS when there are two habitat types. We show that adding environmental stochasticity results in an ESS that, when compared to the ESS for the corresponding model without stochasticity, spends less time in patches with larger carrying capacities and possibly makes use of sink patches, thereby practicing a spatial form of bet hedging.
△ Less
Submitted 20 September, 2014; v1 submitted 27 April, 2014;
originally announced April 2014.
-
Killed Brownian motion with a prescribed lifetime distribution and models of default
Authors:
Boris Ettinger,
Steven N. Evans,
Alexandru Hening
Abstract:
The inverse first passage time problem asks whether, for a Brownian motion $B$ and a nonnegative random variable $ζ$, there exists a time-varying barrier $b$ such that $\mathbb{P}\{B_s>b(s),0\leq s\leq t\}=\mathbb{P}\{ζ>t\}$. We study a "smoothed" version of this problem and ask whether there is a "barrier" $b$ such that $ \mathbb{E}[\exp(-λ\int_0^tψ(B_s-b(s))\,ds)]=\mathbb{P}\{ζ>t\}$, where $λ$ i…
▽ More
The inverse first passage time problem asks whether, for a Brownian motion $B$ and a nonnegative random variable $ζ$, there exists a time-varying barrier $b$ such that $\mathbb{P}\{B_s>b(s),0\leq s\leq t\}=\mathbb{P}\{ζ>t\}$. We study a "smoothed" version of this problem and ask whether there is a "barrier" $b$ such that $ \mathbb{E}[\exp(-λ\int_0^tψ(B_s-b(s))\,ds)]=\mathbb{P}\{ζ>t\}$, where $λ$ is a killing rate parameter, and $ψ:\mathbb{R}\to[0,1]$ is a nonincreasing function. We prove that if $ψ$ is suitably smooth, the function $t\mapsto \mathbb{P}\{ζ>t\}$ is twice continuously differentiable, and the condition $0<-\frac{d\log\mathbb{P}\{ζ>t\}}{dt}<λ$ holds for the hazard rate of $ζ$, then there exists a unique continuously differentiable function $b$ solving the smoothed problem. We show how this result leads to flexible models of default for which it is possible to compute expected values of contingent claims.
△ Less
Submitted 15 January, 2014; v1 submitted 12 November, 2011;
originally announced November 2011.
-
Non-existence of Markovian time dynamics for graphical models of correlated default
Authors:
Steven N. Evans,
Alexandru Hening
Abstract:
Filiz et al. (2008) proposed a model for the pattern of defaults seen among a group of firms at the end of a given time period. The ingredients in the model are a graph, where the vertices correspond to the firms and the edges describe the network of interdependencies between the firms, a parameter for each vertex that captures the individual propensity of that firm to default, and a parameter for…
▽ More
Filiz et al. (2008) proposed a model for the pattern of defaults seen among a group of firms at the end of a given time period. The ingredients in the model are a graph, where the vertices correspond to the firms and the edges describe the network of interdependencies between the firms, a parameter for each vertex that captures the individual propensity of that firm to default, and a parameter for each edge that captures the joint propensity of the two connected firms to default. The correlated default model can be re-rewritten as a standard Ising model on the graph by identifying the set of defaulting firms in the default model with the set of sites in the Ising model for which the spin is +1. We ask whether there is a suitable continuous time Markov chain taking values in the subsets of the vertex set such that the initial state of the chain is the empty set, each jump of the chain involves the inclusion of a single extra vertex, the distribution of the chain at some fixed time horizon time is the one given by the default model, and the distribution of the chain for other times is described by a probability distribution in the same family as the default model. We show for three simple but financially natural special cases that this is not possible outside of the trivial case where there is complete independence between the firms.
△ Less
Submitted 12 August, 2010;
originally announced August 2010.
