Source author record

Alexandru Hening

Alexandru Hening appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2021arXiv

Quasi-stationary distributions of multi-dimensional diffusion processes

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.

preprint2019arXiv

Coexistence, extinction, and optimal harvesting in discrete-time stochastic population models

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.

preprint2019arXiv

Harvesting and seeding of stochastic populations: analysis and numerical approximation

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.

preprint2019arXiv

The competitive exclusion principle in stochastic environments

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.

preprint2019arXiv

The Inverse First Passage Time Problem for killed Brownian motion

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.

preprint2015arXiv

The free path in a high velocity random flight process associated to a Lorentz gas in an external field

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.

preprint2015arXiv

Transient one-dimensional diffusions conditioned to converge to a different limit point

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.

preprint2015arXiv

When do skew-products exist?

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.

preprint2014arXiv

Killed Brownian motion with a prescribed lifetime distribution and models of default

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.

preprint2014arXiv

Protected polymorphisms and evolutionary stability of patch-selection strategies in stochastic environments

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.

preprint2010arXiv

Non-existence of Markovian time dynamics for graphical models of correlated default

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.