BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Zbigniew Palmowski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider a discrete time parallel queue, which is two-queue network, where at each time-slot there is a the same batch arrival to both queues and at each queue there is a random service available. The service law at each time-slot for each queue is different. Let $(Q_n^1, Q_n^2)$ be the queue length at $n$th time-slot. We present several open questions related to the steady-state of this queue.
In this paper we develop the Gerber-Shiu theory for the classic and dual discrete risk processes in a Markovian (regime switching) environment. In particular, by expressing the Gerber-Shiu function in terms of potential measures of an upward (downward) skip-free discrete-time and discrete-space Markov Additive Process (MAP), we derive closed form expressions for the Gerber-Shiu function in terms of the so-called (discrete) $\boldsymbol{W}_v$ and $\boldsymbol{Z}_v$ scale matrices, which were introduced in arXiv:2008.06697. We show that the discrete scale matrices allow for a unified approach for identifying the Gerber-Shiu function as well as the value function of the associated constant dividend barrier problems.
Stochastic fluid-fluid models (SFFMs) offer powerful modeling ability for a wide range of real-life systems of significance. The existing theoretical framework for this class of models is in terms of operator-analytic methods. For the first time, we establish matrix-analytic methods for the efficient analysis of SFFMs. We illustrate the theory with numerical examples.
In this paper we give a solution to the quickest drift change detection problem for a multivariate Lévy process consisting of both continuous (Gaussian) and jump components in the Bayesian approach. We do it for a general 0-modified continuous prior distribution of the change point. Classically, our criterion of optimality is based on a probability of false alarm and an expected delay of the detection, which is then reformulated in terms of a posterior probability of the change point. We find a generator of the posterior probability, which in case of general prior distribution is inhomogeneous in time. The main solving technique uses the optimal stopping theory and is based on solving a certain free-boundary problem. We also construct a Generelized Shiryaev-Roberts statistic, which can be used for applications. The paper is supplemented by two examples, one of which is further used to analyze Polish life tables (after proper calibration) and detect the drift change in the correlated force of mortality of men and women jointly.
We study critical GI/G/1 queues under finite second moment assumptions. We show that the busy period distribution is regularly varying with index half. We also review previously known M/G/1/ and M/M/1 derivations, yielding exact asymptotics as well as a similar derivation for GI/M/1. The busy period asymptotics determine the growth rate of moments of the renewal process counting busy cycles. We further use this to demonstrate a BRAVO phenomenon (Balancing Reduces Asymptotic Variance of Outputs) for the work-output process (namely the busy-time). This yields new insight on the BRAVO effect. A second contribution of the paper is in settling previous conjectured results about GI/G/1 and GI/G/s BRAVO. Previously, infinite buffer BRAVO was generally only settled under fourth-moment assumptions together with an assumption about the tail of the busy-period. In the current paper we strengthen the previous results by reducing to assumptions to existence of $2+ε$ moments.
We consider the sample average of a centered random walk in $\mathbb{R}^d$ with regularly varying step size distribution. For the first exit time from a compact convex set $A$ not containing the origin, we show that its tail is of lognormal type. Moreover, we show that the typical way for a large exit time to occur is by having a number of jumps growing logarithmically in the scaling parameter.
In this paper, we build on the techniques developed in Albrecher et al. (2013), to generate initial-boundary value problems for ruin probabilities of surplus-dependent premium risk processes, under a renewal case scenario, Erlang (2) claim arrivals, and an exponential claims scenario, Erlang (2) claim sizes. Applying the approximation theory of solutions of linear ordinary differential equations developed in Fedoryuk (1993), we derive the asymptotics of the ruin probabilities when the initial reserve tends to infinity. When considering premiums that are {\it linearly} dependent on reserves, representing for instance returns on risk-free investments of the insurance capital, we firstly derive explicit formulas for the ruin probabilities, from which we can easily determine their asymptotics, only to match the ones obtained for general premiums dependent on reserves. We compare them with the asymptotics of the equivalent ruin probabilities when the premium rate is fixed over time, to measure the gain generated by this additional mechanism of binding the premium rates with the amount of reserve own by the insurance company.
In this paper we consider the following optimal stopping problem $$V^ω_{\rm A}(s) = \sup_{τ\in\mathcal{T}} \mathbb{E}_{s}[e^{-\int_0^τω(S_w) dw} g(S_τ)],$$ where the process $S_t$ is a jump-diffusion process, $\mathcal{T}$ is a family of stopping times while $g$ and $ω$ are fixed payoff function and discount function, respectively. In a financial market context, if $g(s)=(K-s)^+$ or $g(s)=(s-K)^+$ and $\mathbb{E}$ is the expectation taken with respect to a martingale measure, $V^ω_{\rm A}(s)$ describes the price of a perpetual American option with a discount rate depending on the value of the asset process $S_t$. If $ω$ is a constant, the above problem produces the standard case of pricing perpetual American options. In the first part of this paper we find sufficient conditions for the convexity of the value function $V^ω_{\rm A}(s)$. This allows us to determine the stopping region as a certain interval and hence we are able to identify the form of $V^ω_{\rm A}(s)$. We also prove a put-call symmetry for American options with asset-dependent discounting. In the case when $S_t$ is a geometric Lévy process we give exact expressions using the so-called omega scale functions introduced in Li and Palmowski (2018). We prove that the analysed value function satisfies the HJB equation and we give sufficient conditions for the smooth fit property as well. Finally, we present a few examples for which we obtain the analytical form of the value function $V^ω_{\rm A}(s)$.
The main objective of this paper is to present an algorithm of pricing perpetual American put options with asset-dependent discounting. The value function of such an instrument can be described as \begin{equation*} V^ω_{\text{A}^{\text{Put}}}(s) = \sup_{τ\in\mathcal{T}} \mathbb{E}_{s}[e^{-\int_0^τω(S_w) dw} (K-S_τ)^{+}], \end{equation*} where $\mathcal{T}$ is a family of stopping times, $ω$ is a discount function and $\mathbb{E}$ is an expectation taken with respect to a martingale measure. Moreover, we assume that the asset price process $S_t$ is a geometric Lévy process with negative exponential jumps, i.e. $S_t = s e^{ζt + σB_t - \sum_{i=1}^{N_t} Y_i}$. The asset-dependent discounting is reflected in the $ω$ function, so this approach is a generalisation of the classic case when $ω$ is constant. It turns out that under certain conditions on the $ω$ function, the value function $V^ω_{\text{A}^{\text{Put}}}(s)$ is convex and can be represented in a closed form; see Al-Hadad and Palmowski (2021). We provide an option pricing algorithm in this scenario and we present exact calculations for the particular choices of $ω$ such that $V^ω_{\text{A}^{\text{Put}}}(s)$ takes a simplified form.
We analyse an additive-increase and multiplicative-decrease (aka growth-collapse) process that grows linearly in time and that experiences downward jumps at Poisson epochs that are (deterministically) proportional to its present position. This process is used for example in modelling of Transmission Control Protocol (TCP) and can be viewed as a particular example of the so-called shot noise model, a basic tool in modeling earthquakes, avalanches and neuron firings. For this process, and also for its reflected versions, we consider one- and two-sided exit problems that concern the identification of the laws of exit times from fixed intervals and half-lines. All proofs are based on a unified first-step analysis approach at the first jump epoch, which allows us to give explicit, yet involved, formulas for their Laplace transforms. All the eight Laplace transforms can be described in terms of two so-called scale functions $Z_{\uparrow}$ and $L_{\uparrow}$. Here $Z_{\uparrow}$ is described in terms of multiple explicit sums, and $L_{\uparrow}$ in terms of an explicit recursion formula. All other Laplace transforms can be obtained from $Z_{\uparrow}$ and $L_{\uparrow}$ by taking limits, derivatives, integrals and combinations of these.
This paper presents an analysis of the stochastic recursion $W_{i+1} = [V_iW_i+Y_i]^+$ that can be interpreted as an autoregressive process of order 1, reflected at 0. We start our exposition by a discussion of the model's stability condition. Writing $Y_i=B_i-A_i$, for independent sequences of non-negative i.i.d.\ random variables $\{A_i\}_{i\in N_0}$ and $\{B_i\}_{i\in N_0}$, and assuming $\{V_i\}_{i\in N_0}$ is an i.i.d. sequence as well (independent of $\{A_i\}_{i\in N_0}$ and $\{B_i\}_{i\in N_0}$), we then consider three special cases: (i) $V_i$ attains negative values only and $B_i$ has a rational LST, (ii) $V_i$ equals a positive value $a$ with certain probability $p\in (0,1)$ and is negative otherwise, and both $A_i$ and $B_i$ have a rational LST, (iii) $V_i$ is uniformly distributed on $[0,1]$, and $A_i$ is exponentially distributed. In all three cases we derive transient and stationary results, where the transient results are in terms of the transform at a geometrically distributed epoch.
This paper describes an application of dynamic programming to determine the optimal strategy for assigning grapes to pressing tanks in one of the largest Portuguese wineries. To date, linear programming has been employed to generate proposed solutions to analogous problems, but this approach lacks robustness and may, in fact, result in severe losses in cases of sudden changes, which frequently occur in weather-dependent wine factories. Hence, we endowed our model with stochasticity, thereby rendering it less vulnerable to such changes. Our analysis, which is based on real-world data, demonstrates that the proposed algorithm is highly efficient and, after calibration, can be used to support winery's decision-making. The solution proposed herein could also be applied in numerous other contexts where production processes rely on outside supplies.
We consider the Lévy model of the perpetual American call and put options with a negative discount rate under Poisson observations. Similar to the continuous observation case as in De Donno et al. [24], the stopping region that characterizes the optimal stopping time is either a half-line or an interval. The objective of this paper is to obtain explicit expressions of the stopping and continuation regions and the value function, focusing on spectrally positive and negative cases. To this end, we compute the identities related to the first Poisson arrival time to an interval via the scale function and then apply those identities to the computation of the optimal strategies. We also discuss the convergence of the optimal solutions to those in the continuous observation case as the rate of observation increases to infinity. Numerical experiments are also provided.
We consider the distributional fixed-point equation: $$R \stackrel{\mathcal{D}}{=} Q \vee \left( \bigvee_{i=1}^N C_i R_i \right),$$ where the $\{R_i\}$ are i.i.d.~copies of $R$, independent of the vector $(Q, N, \{C_i\})$, where $N \in \mathbb{N}$, $Q, \{C_i\} \geq 0$ and $P(Q > 0) > 0$. By setting $W = \log R$, $X_i = \log C_i$, $Y = \log Q$ it is equivalent to the high-order Lindley equation $$W \stackrel{\mathcal{D}}{=} \max\left\{ Y, \, \max_{1 \leq i \leq N} (X_i + W_i) \right\}.$$ It is known that under Kesten assumptions, $$P(W > t) \sim H e^{-αt}, \qquad t \to \infty,$$ where $α>0$ solves the Cramér-Lundberg equation $E \left[ \sum_{j=1}^N C_i ^α\right] = E\left[ \sum_{i=1}^N e^{αX_i} \right] = 1$. The main goal of this paper is to provide an explicit representation for $P(W > t)$, which can be directly connected to the underlying weighted branching process where $W$ is constructed and that can be used to construct unbiased and strongly efficient estimators for all $t$. Furthermore, we show how this new representation can be directly analyzed using Alsmeyer's Markov renewal theorem, yielding an alternative representation for the constant $H$. We provide numerical examples illustrating the use of this new algorithm.
In this work, we adapt a Monte Carlo algorithm introduced by Broadie and Glasserman (1997) to price a $π$-option. This method is based on the simulated price tree that comes from discretization and replication of possible trajectories of the underlying asset's price. As a result this algorithm produces the lower and the upper bounds that converge to the true price with the increasing depth of the tree. Under specific parametrization, this $π$-option is related to relative maximum drawdown and can be used in the real-market environment to protect a portfolio against volatile and unexpected price drops. We also provide some numerical analysis.
This paper considers an optimal dividend distribution problem for an insurance company where the dividends are paid in a foreign currency. In the absence of dividend payments, our risk process follows a spectrally negative Lévy process. We assume that the exchange rate is described by a an exponentially Lévy process, possibly containing the same risk sources like the surplus of the insurance company under consideration. The control mechanism chooses the amount of dividend payments. The objective is to maximise the expected dividend payments received until the time of ruin and a penalty payment at the time of ruin, which is an increasing function of the size of the shortfall at ruin. A complete solution is presented to the corresponding stochastic control problem. Via the corresponding Hamilton--Jacobi--Bellman equation we find the necessary and sufficient conditions for optimality of a single dividend barrier strategy. A number of numerical examples illustrate the theoretical analysis.
This paper discusses the valuation of credit default swaps, where default is announced when the reference asset price has gone below certain level from the last record maximum, also known as the high-water mark or drawdown. We assume that the protection buyer pays premium at fixed rate when the asset price is above a pre-specified level and continuously pays whenever the price increases. This payment scheme is in favour of the buyer as she only pays the premium when the market is in good condition for the protection against financial downturn. Under this framework, we look at an embedded option which gives the issuer an opportunity to call back the contract to a new one with reduced premium payment rate and slightly lower default coverage subject to paying a certain cost. We assume that the buyer is risk neutral investor trying to maximize the expected monetary value of the option over a class of stopping time. We discuss optimal solution to the stopping problem when the source of uncertainty of the asset price is modelled by Lévy process with only downward jumps. Using recent development in excursion theory of Lévy process, the results are given explicitly in terms of scale function of the Lévy process. Furthermore, the value function of the stopping problem is shown to satisfy continuous and smooth pasting conditions regardless of regularity of the sample paths of the Lévy process. Optimality and uniqueness of the solution are established using martingale approach for drawdown process and convexity of the scale function under Esscher transform of measure. Some numerical examples are discussed to illustrate the main results.
We consider in this paper a risk reserve process where the claims and gains arrive according to two independent Poisson processes. While the gain sizes are phase-type distributed, we assume instead that the claim sizes are phase-type perturbed by a heavy-tailed component; that is, the claim size distribution is formally chosen to be phase-type with large probability $1-ε$ and heavy-tailed with small probability $ε$. We analyze the seminal Gerber-Shiu function coding the joint distribution of the time to ruin, the surplus immediately before ruin, and the deficit at ruin. We derive its value as an expansion with respect to powers of $ε$ with known coefficients and we construct approximations from the first two terms of the aforementioned series. The main idea is based on the so-called fluid embedding that allows to put the considered risk process into the framework of spectrally negative Markov-additive processes and use its fluctuation theory developed in Ivanovs and Palmowski (2012).
We revisit the optimal capital structure model with endogenous bankruptcy first studied by Leland \cite{Leland94} and Leland and Toft \cite{Leland96}. Differently from the standard case, where shareholders observe continuously the asset value and bankruptcy is executed instantaneously without delay, we assume that the information of the asset value is updated only at intervals, modeled by the jump times of an independent Poisson process. Under the spectrally negative Lévy model, we obtain the optimal bankruptcy strategy and the corresponding capital structure. A series of numerical studies are given to analyze the sensitivity of observation frequency on the optimal solutions, the optimal leverage and the credit spreads.
In this paper we provide the analysis of the limiting conditional distribution (Yaglom limit) for stochastic fluid models (SFMs), a key class of models in the theory of matrix-analytic methods. So far, transient and stationary analyses of the SFMs have been only considered in the literature. The limiting conditional distribution gives useful insights into what happens when the process has been evolving for a long time, given its busy period has not ended yet. We derive expressions for the Yaglom limit in terms of the singularity $s^*$ such that the key matrix of the SFM, ${\bfΨ}(s)$, is finite (exists) for all $s\geq s^*$ and infinite for $s<s^*$. We show the uniqueness of the Yaglom limit and illustrate the application of the theory with simple examples.
In this paper we identify the asymptotic tail of the distribution of the exit time $τ_C$ from a cone $C$ of an isotropic $α$-self-similar Markov process $X_t$ with a skew-product structure, that is $X_t$ is a product of its radial process and independent time changed angular component $Θ_t$. Under some additional regularity assumptions, the angular process $Θ_t$ killed on exiting from the cone $C$ has the transition density that could be expressed in terms of a complete set of orthogonal eigenfunctions with corresponding eigenvalues of an appropriate generator. Using this fact and some asymptotic properties of the exponential functional of a killed Lévy process related with Lamperti representation of the radial process, we prove that $$\mathbb{P}_x(τ_C>t)\sim h(x)t^{-κ_1}$$ as $t\rightarrow\infty$ for $h$ and $κ_1$ identified explicitly. The result extends the work of DeBlassie (1988) and Bañuelos and Smits (1997) concerning the Brownian motion.
This paper concerns an optimal dividend distribution problem for an insurance company with surplus-dependent premium. In the absence of dividend payments, such a risk process is a particular case of so-called piecewise deterministic Markov processes. The control mechanism chooses the size of dividend payments. The objective consists in maximazing the sum of the expected cumulative discounted dividend payments received until the time of ruin and a penalty payment at the time of ruin, which is an increasing function of the size of the shortfall at ruin. A complete solution is presented to the corresponding stochastic control problem. We identify the associated Hamilton-Jacobi-Bellman equation and find necessary and sufficient conditions for optimality of a single dividend-band strategy, in terms of particular Gerber-Shiu functions. A number of concrete examples are analyzed.
This paper concerns the dual risk model, dual to the risk model for insurance applications, where premiums are surplus-dependent. In such a model premiums are regarded as costs, while claims refer to profits. We calculate the mean of the cumulative discounted dividends paid until ruin, if the barrier strategy is applied. We formulate associated Hamilton-Jacobi-Bellman equation and identify sufficient conditions for a barrier strategy to be optimal. Some numerical examples are provided when profits have exponential law.
In this paper we consider a classical risk process perturbed by a Brownian motion. We analyze the value function describing the mean of the cumulative discounted dividend payments paid up to Parisian ruin time and further discounted by the number of claims appeared up to this ruin time. We identify this value function for the barrier strategy and find the sufficient conditions for this strategy to be optimal. We also consider few particular examples.
In recent years there has been some focus on quasi-stationary behaviour of an one-dimensional Lévy process $X$, where we ask for the law $P(X_t\in dy | τ^-_0>t)$ for $t\to\infty$ and $τ_0^-=\inf\{t\geq 0: X_t<0\}$. In this paper we address the same question for so-called Parisian ruin time $τ^θ$, that happens when process stays below zero longer than independent exponential random variable with intensity $θ$.
We give a review of the state of the art with regard to the dividend problem.
We analyse the ruin probabilities for a renewal insurance risk process with inter-arrival time distributions depending on the claims that arrived within a fixed (past) time window. This dependence could be explained through a regenerative structure. The main inspiration of the model comes from the Bonus-Malus feature. We discuss first asymptotic results of ruin probabilities for different regimes of claim distributions. For numerical results, we recognise an embedded Markov additive process. Via an appropriate change of measure, ruin probabilities could be computed to a closed form formulae. Additionally, we present simulated results via the importance sampling method, which further permit an in-depth analysis of a few concrete cases.
In this paper we propose new iterative algorithm of calculating the joint distribution of the Parisian ruin time and the number of claims until Parisian ruin for the classical risk model. Examples are provided when the generic claim size is exponentially distributed.
We give the asymptotics of the tail of the distribution of the first exit time of the isotropic $α$-stable Lévy process from the Lipschitz cone in $\mathbb{R}^d$. We obtain the Yaglom limit for the killed stable process for the cone. We construct and estimate entrance laws for the process from the vertex into the cone. For the symmetric Cauchy process and the positive half-line we give a spectral representation of the Yaglom limit. Our approach relies on the scalings of the stable process and the cone, which allow to express the temporal asymptotics of the distribution of the process at infinity by means of the spatial asymptotics of harmonic functions of the process at the vertex; on the representation of the probability of survival of the process in the cone as a Green potential; and on the approximate factorization of the heat kernel of the cone, which secures compactness and yields a limiting (Yaglom) measure by means of Prokhorov's theorem.
We investigate the tail behaviour of the steady state distribution of a stochastic recursion that generalises Lindley's recursion. This recursion arises in queuing systems with dependent interarrival and service times, and includes alternating service systems and carousel storage systems as special cases. We obtain precise tail asymptotics in three qualitatively different cases, and compare these with existing results for Lindley's recursion and for alternating service systems.
In this article we derive formulas for the probability $P(\sup_{t\leq T} X(t)>u)$ $T>0$ and $P(\sup_{t<\infty} X(t)>u)$ where $X$ is a spectrally positive Lévy process with infinite variation. The formulas are generalizations of the well-known Takács formulas for stochastic processes with non-negative and interchangeable increments. Moreover, we find the joint distribution of $\inf_{t\leq T} Y(t)$ and $Y(T)$ where $Y$ is a spectrally negative Lévy process.
We consider a two-dimensional stochastic fluid model with $N$ ON-OFF inputs and temporary assistance, which is an extension of the same model with $N = 1$ in Mahabhashyam et al. (2008). The rates of change of both buffers are piecewise constant and dependent on the underlying Markovian phase of the model, and the rates of change for Buffer 2 are also dependent on the specific level of Buffer 1. This is because both buffers share a fixed output capacity, the precise proportion of which depends on Buffer 1. The generalization of the number of ON-OFF inputs necessitates modifications in the original rules of output-capacity sharing from Mahabhashyam et al. (2008) and considerably complicates both the theoretical analysis and the numerical computation of various performance measures.
In this note we give, for a spectrally negative Levy process, a compact formula for the Parisian ruin probability, which is defined by the probability that the process exhibits an excursion below zero, with a length that exceeds a certain fixed period r. The formula involves only the scale function of the spectrally negative Levy process and the distribution of the process at time r.
In this paper we analyze a branching process with immigration defined recursively by $X_t=θ_t\circ X_{t-1}+B_t$ for a sequence $(B_t)$ of i.i.d. random variables and random mappings $ θ_t\circ x:=θ_t(x)=\sum_{i=1}^xA_i^{(t)}, $ with $(A_i^{(t)})_{i\in \mathbb{N}_0}$ being a sequence of $\mathbb{N}_0$-valued i.i.d. random variables independent of $B_t$. We assume that one of generic variables $A$ and $B$ has a regularly varying tail distribution. We identify the tail behaviour of the distribution of the stationary solution $X_t$. We also prove CLT for the partial sums that could be further generalized to FCLT. Finally, we also show that partial maxima have a Fréchet limiting distribution.
We consider a queuing model with the workload evolving between consecutive i.i.d.\ exponential timers $\{e_q^{(i)}\}_{i=1,2,...}$ according to a spectrally positive Lévy process $Y_i(t)$ that is reflected at zero, and where the environment $i$ equals 0 or 1. When the exponential clock $e_q^{(i)}$ ends, the workload, as well as the Lévy input process, are modified; this modification may depend on the current value of the workload, the maximum and the minimum workload observed during the previous cycle, and the environment $i$ of the Lévy input process itself during the previous cycle. We analyse the steady-state workload distribution for this model. The main theme of the analysis is the systematic application of non-trivial functionals, derived within the framework of fluctuation theory of Lévy processes, to workload and queuing models.
We prove the Wiener-Hopf factorization for Markov Additive processes. We derive also Spitzer-Rogozin theorem for this class of processes which serves for obtaining Kendall's formula and Fristedt representation of the cumulant matrix of the ladder epoch process. Finally, we also obtain the so-called ballot theorem.
In this paper we consider dividend problem for an insurance company whose risk evolves as a spectrally negative Lévy process (in the absence of dividend payments) when Parisian delay is applied. The objective function is given by the cumulative discounted dividends received until the moment of ruin when so-called barrier strategy is applied. Additionally we will consider two possibilities of delay. In the first scenario ruin happens when the surplus process stays below zero longer than fixed amount of time $ζ>0$. In the second case there is a time lag $d$ between decision of paying dividends and its implementation.
In this paper we develop a symbolic technique to obtain asymptotic expressions for ruin probabilities and discounted penalty functions in renewal insurance risk models when the premium income depends on the present surplus of the insurance portfolio. The analysis is based on boundary problems for linear ordinary differential equations with variable coefficients. The algebraic structure of the Green's operators allows us to develop an intuitive way of tackling the asymptotic behavior of the solutions, leading to exponential-type expansions and Cramér-type asymptotics. Furthermore, we obtain closed-form solutions for more specific cases of premium functions in the compound Poisson risk model.
In this paper we consider a general Lévy process $X$ reflected at downward periodic barrier $A_t$ and constant upper barrier $K$ giving a process $V^K_t=X_t+L^A_t-L^K_t$. We find the expression for a loss rate defined by $l^K=\mathbb{E} L^K_1$ and identify its asymptotics as $K\to\infty$ when $X$ has light-tailed jumps and $\mathbb{E}X_1<0$.
This paper solves exit problems for spectrally negative Markov additive processes and their reflections. A so-called scale matrix, which is a generalization of the scale function of a spectrally negative \levy process, plays a central role in the study of exit problems. Existence of the scale matrix was shown in Thm. 3 of Kyprianou and Palmowski (2008). We provide a probabilistic construction of the scale matrix, and identify the transform. In addition, we generalize to the MAP setting the relation between the scale function and the excursion (height) measure. The main technique is based on the occupation density formula and even in the context of fluctuations of spectrally negative Lévy processes this idea seems to be new. Our representation of the scale matrix $W(x)=e^{-Λx}\eL(x)$ in terms of nice probabilistic objects opens up possibilities for further investigation of its properties.
We study the time reversal of a general PDMP. The time reversed process is defined as $X_{(T-t)-}$, where $T$ is some given time and $X_t$ is a stationary PDMP. We obtain the parameters of the reversed process, like the jump intensity and the jump measure.
In this paper we analyze the quasi-stationary workload of a Lévy-driven storage system. More precisely, assuming the system is in stationarity, we study its behavior conditional on the event that the busy period $T$ in which time 0 is contained has not ended before time $t$, as $t\to\infty$. We do so by first identifying the double Laplace transform associated with the workloads at time 0 and time $t$, on the event $\{T>t\}.$ This transform can be explicitly computed for the case of spectrally one-sided jumps. Then asymptotic techniques for Laplace inversion are relied upon to find the corresponding behavior in the limiting regime that $t\to\infty.$ Several examples are treated; for instance in the case of Brownian input, we conclude that the workload distribution at time 0 and $t$ are both Erlang(2).
The areas under workload process and under queuing process in a single server queue over the busy period have many applications not only in queuing theory but also in risk theory or percolation theory. We focus here on the tail behaviour of distribution of these two integrals. We present various open problems and conjectures, which are supported by partial results for some special cases.
In this paper we analyze so-called Parisian ruin probability that happens when surplus process stays below zero longer than fixed amount of time $ζ>0$. We focus on general spectrally negative Lévy insurance risk process. For this class of processes we identify expression for ruin probability in terms of some other quantities that could be possibly calculated explicitly in many models. We find its Cramér-type and convolution-equivalent asymptotics when reserves tends to infinity. Finally, we analyze few explicit examples.
We consider a queuing model with the workload evolving between consecutive i.i.d. exponential timers $\{e_q^{(i)}\}_{i=1,2,...}$ according to a spectrally positive Lévy process $Y(t)$ which is reflected at 0. When the exponential clock $e_q^{(i)}$ ends, the additional state-dependent service requirement modifies the workload so that the latter is equal to $F_i(Y(e_q^{(i)}))$ at epoch $e^{(1)}_q+...+e^{(i)}_q$ for some random nonnegative i.i.d. functionals $F_i$. In particular, we focus on the case when $F_i(y)=(B_i-y)^+$, where $\{B_i\}_{i=1,2,...}$ are i.i.d. nonnegative random variables. We analyse the steady-state workload distribution for this model.