BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alexey Kuznetsov

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

ResearcherUnclaimed source record
math.PRmath.NAastro-ph.SRmath.CAApplicationsmath.CVNumerical Analysismath.APmath.DSmath.NTq-fin.PR

Catalog footprint

What is connected

32works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

32 published item(s)

preprint2022arXiv

Computing the Barnes $G$-function and the gamma function in the entire complex plane

We present an algorithm for generating approximations for the logarithm of Barnes $G$-function in the half-plane $Re(z)\ge 3/2$. These approximations involve only elementary functions and are easy to implement. The algorithm is based on a two-point Padé approximation and we use it to provide two approximations to $\ln(G(z))$, accurate to $3 \times 10^{-16}$ and $3 \times 10^{-31}$ in the half-plane $Re(z)\ge 3/2$; a reflection formula is then used to compute Barnes $G$-function in the entire complex plane. A by-product of our algorithm is that it also produces accurate approximations to the gamma function.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Prediction of Personal Protective Equipment Use in Hospitals During COVID-19

Demand for Personal Protective Equipment (PPE) such as surgical masks, gloves, and gowns has increased significantly since the onset of the COVID-19 pandemic. In hospital settings, both medical staff and patients are required to wear PPE. As these facilities resume regular operations, staff will be required to wear PPE at all times while additional PPE will be mandated during medical procedures. This will put increased pressure on hospitals which have had problems predicting PPE usage and sourcing its supply. To meet this challenge, we propose an approach to predict demand for PPE. Specifically, we model the admission of patients to a medical department using multiple independent queues. Each queue represents a class of patients with similar treatment plans and hospital length-of-stay. By estimating the total workload of each class, we derive closed-form estimates for the expected amount of PPE required over a specified time horizon using current PPE guidelines. We apply our approach to a data set of 22,039 patients admitted to the general internal medicine department at St. Michael's hospital in Toronto, Canada from April 2010 to November 2019. We find that gloves and surgical masks represent approximately 90% of predicted PPE usage. We also find that while demand for gloves is driven entirely by patient-practitioner interactions, 86% of the predicted demand for surgical masks can be attributed to the requirement that medical practitioners will need to wear them when not interacting with patients.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the rate of convergence of the Gaver-Stehfest algorithm

The Gaver-Stehfest algorithm is widely used for numerical inversion of Laplace transform. In this paper we provide the first rigorous study of the rate of convergence of the Gaver-Stehfest algorithm. We prove that Gaver-Stehfest approximations converge exponentially fast if the target function is analytic in a neighbourhood of a point and they converge at a rate $o(n^{-k})$ if the target function is $(2k+3)$-times differentiable at a point.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A short proof of duality relations for hypergeometric functions

Identities involving finite sums of products of hypergeometric functions and their duals have been studied since 1930s. Recently Beukers and Jouhet have used an algebraic approach to derive a very general family of duality relations. In this paper we provide an alternative way of obtaining such results. Our method is very simple and it is based on the non-local derangement identity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Approximating Lévy processes with completely monotone jumps

Lévy processes with completely monotone jumps appear frequently in various applications of probability. For example, all popular stock price models based on Lévy processes (such as the Variance Gamma, CGMY/KoBoL and Normal Inverse Gaussian) belong to this class. In this paper we continue the work started in [Int. J. Theor. Appl. Finance 13 (2010) 63-91, Quant. Finance 10 (2010) 629-644] and develop a simple yet very efficient method for approximating processes with completely monotone jumps by processes with hyperexponential jumps, the latter being the most convenient class for performing numerical computations. Our approach is based on connecting Lévy processes with completely monotone jumps with several areas of classical analysis, including Padé approximations, Gaussian quadrature and orthogonal polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constructing measures with identical moments

The Nevanlinna parametrization establishes a bijection between the class of all measures having a prescribed set of moments and the class of Pick functions. The fact that all measures constructed through the Nevanlinna parametrization have identical moments follows from the theory of orthogonal polynomials and continued fractions. In this paper we explore the opposite direction: we construct a set of measures and we show that they all have identical moments, and then we establish a Nevanlinna-type parametrization for this set of measures. Our construction does not require the theory of orthogonal polynomials and it exposes the analytic structure behind the Nevanlinna parametrization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exponential functionals of Levy processes and variable annuity guaranteed benefits

Exponential functionals of Brownian motion have been extensively studied in financial and insurance mathematics due to their broad applications, for example, in the pricing of Asian options. The Black-Scholes model is appealing because of mathematical tractability, yet empirical evidence shows that geometric Brownian motion does not adequately capture features of market equity returns. One popular alternative for modeling equity returns consists in replacing the geometric Brownian motion by an exponential of a Levy process. In this paper we use this latter model to study variable annuity guaranteed benefits and to compute explicitly the distribution of certain exponential functionals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Synthetic radio views on simulated solar flux ropes

In this paper, we produce synthetic radio views on simulated flux ropes in the solar corona, where finite-beta magnetohydrodynamic (MHD) simulations serve to mimic the flux rope formation stages, as well as their stable endstates. These endstates represent twisted flux ropes where balancing Lorentz forces, gravity and pressure gradients determine the full thermodynamic variation throughout the flux rope. The obtained models are needed to quantify radiative transfer in radio bands, and allow us to contrast weak to strong magnetic field conditions. Field strengths of up to 100 G in the flux rope yield the radio views dominated by optically thin free-free emission. The forming flux rope shows clear morphological changes in its emission structure as it deforms from an arcade to a flux rope, both on disk and at the limb. For an active region filament channel with a field strength of up to 680 G in the flux rope, gyroresonance emission (from the third-fourth gyrolayers) can be detected and even dominates over free-free emission at the frequencies of up to 7 GHz. Finally, we also show synthetic views on a simulated filament embedded within a (weak-field) flux rope, resulting from an energetically consistent MHD simulation. For this filament, synthetic views at the limb show clear similarities with actual observations, and the transition from optically thick (below 10 GHz) to optically thin emission can be reproduced. On the disk, its dimension and temperature conditions are as yet not realistic enough to yield the observed radio brightness depressions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tail dependence of the Gaussian copula revisited

Tail dependence refers to clustering of extreme events. In the context of financial risk management, the clustering of high-severity risks has a devastating effect on the well-being of firms and is thus of pivotal importance in risk analysis.When it comes to quantifying the extent of tail dependence, it is generally agreed that measures of tail dependence must be independent of the marginal distributions of the risks but rather solely copula-dependent. Indeed, all classical measures of tail dependence are such, but they investigate the amount of tail dependence along the main diagonal of copulas, which has often little in common with the concentration of extremes in the copulas' domain of definition.In this paper we urge that the classical measures of tail dependence may underestimate the level of tail dependence in copulas. For the Gaussian copula, however, we prove that the classical measures are maximal. The implication of the result is two-fold: On the one hand, it means that in the Gaussian case, the (weak) measures of tail dependence that have been reported and used are of utmost prudence, which must be a reassuring news for practitioners. On the other hand, it further encourages substitution of the Gaussian copula with other copulas that are more tail dependent.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Explicit Hermite-type eigenvectors of the discrete Fourier transform

The search for a canonical set of eigenvectors of the discrete Fourier transform has been ongoing for more than three decades. The goal is to find an orthogonal basis of eigenvectors which would approximate Hermite functions -- the eigenfunctions of the continuous Fourier transform. This eigenbasis should also have some degree of analytical tractability and should allow for efficient numerical computations. In this paper we provide a partial solution to these problems. First, we construct an explicit basis of (non-orthogonal) eigenvectors of the discrete Fourier transform, thus extending the results of [7]. Applying the Gramm-Schmidt orthogonalization procedure we obtain an orthogonal eigenbasis of the discrete Fourier transform. We prove that the first eight eigenvectors converge to the corresponding Hermite functions, and we conjecture that this convergence result remains true for all eigenvectors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fractional Laplace operator and Meijer G-function

We significantly expand the number of functions whose image under the fractional Laplace operator can be computed explicitly. In particular, we show that the fractional Laplace operator maps Meijer G-functions of |x|^2, or generalized hypergeometric functions of -|x|^2, multiplied by a solid harmonic polynomial, into the same class of functions. As one important application of this result, we produce a complete system of eigenfunctions of the operator (1-|x|^2)_+^{alpha/2} (-Delta)^{alpha/2} with the Dirichlet boundary conditions outside of the unit ball. The latter result will be used to estimate the eigenvalues of the fractional Laplace operator in the unit ball in a companion paper "Eigenvalues of the fractional Laplace operator in the unit ball".

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulations of Gyrosynchrotron Microwave Emission from an Oscillating 3D Magnetic Loop

Radio observations of solar flares often reveal various periodic or quasi-periodic oscillations. Most likely, these oscillations are caused by magnetohydrodynamic (MHD) oscillations of flaring loops which modulate the emission. Interpretation of the observations requires comparing them with simulations. We simulate the gyrosynchrotron radio emission from a semi-circular (toroidal-shaped) magnetic loop containing sausage-mode MHD oscillations. The aim is to detect the observable signatures specific to the considered MHD mode and to study their dependence on the various source parameters. The MHD waves are simulated using a linear three-dimensional model of a magnetized plasma cylinder; both standing and propagating waves are considered. The curved loop is formed by replicating the MHD solutions along the plasma cylinder and bending the cylinder; this model allows us to study the effect of varying the viewing angle along the loop. The radio emission is simulated using a three-dimensional model and its spatial and temporal variations are analyzed. We consider several loop orientations and different parameters of the magnetic field, plasma, and energetic electrons in the loop. In the model with low plasma density, the intensity oscillations at all frequencies are synchronous (with the exception of a narrow spectral region below the spectral peak). In the model with high plasma density, the emission at low frequencies (where the Razin effect is important) oscillates in anti-phase with the emissions at higher frequencies. The oscillations at high and low frequencies are more pronounced in different parts of the loop (depending on the loop orientation). The layers where the line-of-sight component of the magnetic field changes sign can produce additional peculiarities in the oscillation patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computing the truncated theta function via Mordell integral

Hiary [3] has presented an algorithm which allows to evaluate the truncated theta function $\sum_{k=0}^n \exp(2πı(zk+τk^2))$ to within $\pm ε$ in $O(\ln(\tfrac{n}ε)^κ)$ arithmetic operations for any real $z$ and $τ$. This remarkable result has many applications in Number Theory, in particular it is the crucial element in Hiary's algorithm for computing $ζ(\tfrac{1}{2}+ıt)$ to within $\pm t^{-λ}$ in $O_λ(t^{\frac{1}{3}}\ln(t)^κ)$ arithmetic operations, see [2]. We present a significant simplification of Hiary's algorithm for evaluating the truncated theta function. Our method avoids the use of the Poisson summation formula, and substitutes it with an explicit identity involving the Mordell integral. This results in an algorithm which is efficient, conceptually simple and easy to implement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

New families of subordinators with explicit transition probability semigroup

There exist only a few known examples of subordinators for which the transition probability density can be computed explicitly along side an expression for its Lévy measure and Laplace exponent. Such examples are useful in several areas of applied probability, for example, they are used in mathematical finance for modeling stochastic time change, they appear in combinatorial probability to construct sampling formulae, which in turn is related to a variety of issues in the theory of coalescence models, moreover, they have also been extensively used in the potential analysis of subordinated Brownian motion in dimension greater than or equal to 2. In this paper, we show that Kendall's classic identity for spectrally negative Lévy processes can be used to construct new families of subordinators with explicit transition probability semigroups. We describe the properties of these new subordinators and emphasise some interesting connections with explicit and previously unknown Laplace transform identities and with complete monotonicity properties of certain special functions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatially-resolved Energetic Electron Properties for the 21 May 2004 Flare from Radio Observations and 3D Simulations

We investigate in detail the 21 May 2004 flare using simultaneous observations of the Nobeyama Radioheliograph, Nobeyama Radiopolarimeters, Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and Solar and Heliospheric Observatory (SOHO). The flare images in different spectral ranges reveal the presence of a well-defined single flaring loop in this event. We have simulated the gyrosynchrotron microwave emission using the recently developed interactive IDL tool GX Simulator. By comparing the simulation results with the observations, we have deduced the spatial and spectral properties of the non-thermal electron distribution. The microwave emission has been found to be produced by the high-energy electrons ($>100$ keV) with a relatively hard spectrum ($δ\simeq 2$); the electrons were strongly concentrated near the loop top. At the same time, the number of high-energy electrons near the footpoints was too low to be detected in the RHESSI images and spatially unresolved data. The SOHO Extreme-ultraviolet Imaging Telescope images and the low-frequency microwave spectra suggest the presence of an extended "envelope" of the loop with lower magnetic field. Most likely, the energetic electron distribution in the considered flare reflects the localized (near the loop top) particle acceleration (injection) process accompanied by trapping and scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strongly reinforced Pólya urns with graph-based competition

We introduce a class of reinforcement models where, at each time step $t$, one first chooses a random subset $A_t$ of colours (independent of the past) from $n$ colours of balls, and then chooses a colour $i$ from this subset with probability proportional to the number of balls of colour $i$ in the urn raised to the power $α>1$. We consider stability of equilibria for such models and establish the existence of phase transitions in a number of examples, including when the colours are the edges of a graph, a context which is a toy model for the formation and reinforcement of neural connections.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The hitting time of zero for a stable process

For any two-sided jumping $α$-stable process, where $1 < α< 2$, we find an explicit identity for the law of the first hitting time of the origin. This complements existing work in the symmetric case and the spectrally one-sided case; cf. Yano-Yano-Yor (2009) and Cordero (2010), and Peskir (2008) respectively. We appeal to the Lamperti-Kiu representation of Chaumont-Pantí-Rivero (2011) for real-valued self-similar Markov processes. Our main result follows by considering a vector-valued functional equation for the Mellin transform of the integrated exponential Markov additive process in the Lamperti-Kiu representation. We conclude our presentation with some applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A note on the series representation for the density of the supremum of a stable process

An absolutely convergent double series representation for the density of the supremum of $α$-stable Levy process is given in [3, Theorem 2] for almost all irrational $α$. This result cannot be made stronger in the following sense: the series does not converge absolutely when $α$ belongs to a certain subset of irrational numbers of Lebesgue measure zero (see [6, Theorem 2]). Our main result in this note shows that for every irrational $α$ there is a way to rearrange the terms of the double series, so that it converges to the density of the supremum. We show how one can establish this stronger result by introducing a simple yet non-trivial modification in the original proof of [3,Theorem 2].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Asian options and meromorphic Levy processes

One method to compute the price of an arithmetic Asian option in a Levy driven model is based on the exponential functional of the underlying Levy process: If we know the distribution of the exponential functional, we can calculate the price of the Asian option via the inverse Laplace transform. In this paper we consider pricing Asian options in a model driven by a general meromorphic Levy process. We prove that the exponential functional is equal in distribution to an infinite product of indepedent beta random variables, and its Mellin transform can be expressed as an infinite product of gamma functions. We show that these results lead to an efficient algorithm for computing the price of the Asian option via the inverse Mellin-Laplace transform, and we compare this method with some other techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the convergence of the Gaver-Stehfest algorithm

The Gaver-Stehfest algorithm for numerical inversion of Laplace transform was developed in the late 1960s. Due to its simplicity and good performance it is becoming increasingly more popular in such diverse areas as Geophysics, Operations Research and Economics, Financial and Actuarial Mathematics, Computational Physics and Chemistry. Despite the large number of applications and numerical studies, this method has never been rigorously investigated. In particular, it is not known whether the Gaver-Stehfest approximations converge and what is the rate of convergence. In this paper we answer the first of these two questions: We prove that the Gaver-Stehfest approximations converge for functions of bounded variation and functions satisfying an analogue of Dini criterion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A convergent series representation for the density of the supremum of a stable process

We study the density of the supremum of a strictly stable Lévy process. We prove that for almost all values of the index $α$ -- except for a dense set of Lebesgue measure zero -- the asymptotic series which were obtained in A. Kuznetsov (2010) "On extrema of stable processes" are in fact absolutely convergent series representations for the density of the supremum.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic approximations to the Hardy-Littlewood function

The function $Q(x):=\sum_{n\ge 1} (1/n) \sin(x/n)$ was introduced by Hardy and Littlewood [10] in their study of Lambert summability, and since then it has attracted attention of many researchers. In particular, this function has made a surprising appearance in the recent disproof by Alzer, Berg and Koumandos [1] of a conjecture by Clark and Ismail [2]. More precisely, Alzer et. al. [1] have shown that the Clark and Ismail conjecture is true if and only if $Q(x)\ge -π/2$ for all $x>0$. It is known that $Q(x)$ is unbounded in the domain $x \in (0,\infty)$ from above and below, which disproves the Clark and Ismail conjecture, and at the same time raises a natural question of whether we can exhibit at least one point $x$ for which $Q(x) < -π/2$. This turns out to be a surprisingly hard problem, which leads to an interesting and non-trivial question of how to approximate $Q(x)$ for very large values of $x$. In this paper we continue the work started by Gautschi in [7] and develop several approximations to $Q(x)$ for large values of $x$. We use these approximations to find an explicit value of $x$ for which $Q(x)<-π/2$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Wiener-Hopf factorization for a family of Levy processes related to theta functions

In this paper we study the Wiener-Hopf factorization for a class of Lévy processes with double-sided jumps, characterized by the fact that the density of the Lévy measure is given by an infinite series of exponential functions with positive coefficients. We express the Wiener-Hopf factors as infinite products over roots of a certain transcendental equation, and provide a series representation for the distribution of the supremum/infimum process evaluated at an independent exponential time. We also introduce five eight-parameter families of Lévy processes, defined by the fact that the density of the Lévy measure is a (fractional) derivative of the theta-function, and we show that these processes can have a wide range of behavior of small jumps. These families of processes are of particular interest for applications, since the characteristic exponent has a simple expression, which allows efficient numerical computation of the Wiener-Hopf factors and distributions of various functionals of the process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Meromorphic Levy processes and their fluctuation identities

The last couple of years has seen a remarkable number of new, explicit examples of the Wiener-Hopf factorization for Levy processes where previously there had been very few. We mention in particular the many cases of spectrally negative Levy processes, hyper-exponential and generalized hyper-exponential Levy processes, Lamperti-stable processes, Hypergeometric processes, Beta-processes and Theta-processes. In this paper we introduce a new family of Levy processes, which we call Meromorphic Levy processes, or just M-processes for short, which overlaps with many of the aforementioned classes. A key feature of the M-class is the identification of their Wiener-Hopf factors as rational functions of infinite degree written in terms of poles and roots of the Levy-Khintchin exponent, all of which appear on the imaginary axis of the complex plane. The specific structure of the M-class Wiener-Hopf factorization enables us to explicitly handle a comprehensive suite of fluctuation identities that concern first passage problems for finite and infinite intervals for both the process itself as well as the resulting process when it is reflected in its infimum. Such identities are of fundamental interest given their repeated occurrence in various fields of applied probability such as mathematical finance, insurance risk theory and queuing theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On extrema of stable processes

We study the Wiener--Hopf factorization and the distribution of extrema for general stable processes. By connecting the Wiener--Hopf factors with a certain elliptic-like function we are able to obtain many explicit and general results, such as infinite series representations and asymptotic expansions for the density of supremum, explicit expressions for the Wiener--Hopf factors and the Mellin transform of the supremum, quasi-periodicity and functional identities for these functions, finite product representations in some special cases and identities in distribution satisfied by the supremum functional.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the density of the supremum of a stable process

We study the density of the supremum of a strictly stable Lévy process. As was proved recently in F. Hubalek and A. Kuznetsov "A convergent series representation for the density of the supremum of a stable process" (Elect. Comm. in Probab., 16, 84-95, 2011), for almost all irrational values of the stability parameter $α$ this density can be represented by an absolutely convergent series. We show that this result is not valid for all irrational values of $α$: we construct a dense uncountable set of irrational numbers $α$ for which the series does not converge absolutely. Our second goal is to investigate in more detail the important case when $α$ is rational. We derive an explicit formula for the Mellin transform of the supremum, which is given in terms of Gamma function and dilogarithm. In order to illustrate the usefulness of these results we perform several numerical experiments and discuss their implications. Finally, we state some interesting connections that this problem has to other areas of Mathematics and Mathematical Physics, such as q-series, Diophantine approximations and quantum dilogarithms, and we also suggest several open problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the Wiener-Hopf factorization for Levy processes with bounded positive jumps

We study the Wiener-Hopf factorization for Levy processes with bounded positive jumps and arbitrary negative jumps. Using the results from the theory of entire functions of Cartwright class we prove that the positive Wiener-Hopf factor can be expressed as an infinite product in terms of the solutions to the equation $ψ(z)=q$, where $ψ$ is the Laplace exponent of the process. Under some additional regularity assumptions on the Levy measure we obtain an asymptotic expression for these solutions, which is important for numerical computations. In the case when the process is spectrally negative with bounded jumps, we derive a series representation for the scale function in terms of the solutions to the equation $ψ(z)=q$. To illustrate possible applications we discuss the implementation of numerical algorithms and present the results of several numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The theory of scale functions for spectrally negative Le vy processes

The purpose of this review article is to give an up to date account of the theory and application of scale functions for spectrally negative Levy processes. Our review also includes the first extensive overview of how to work numerically with scale functions. Aside from being well acquainted with the general theory of probability, the reader is assumed to have some elementary knowledge of Levy processes, in particular a reasonable understanding of the Levy-Khintchine formula and its relationship to the Levy-Ito decomposition. We shall also touch on more general topics such as excursion theory and semi-martingale calculus. However, wherever possible, we shall try to focus on key ideas taking a selective stance on the technical details. For the reader who is less familiar with some of the mathematical theories and techniques which are used at various points in this review, we note that all the necessary technical background can be found in the following texts on Levy processes; Bertoin (1996), Sato (1999), Applebaum (2004), Kyprianou (2006) and Doney (2007).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fluctuations of stable processes and exponential functionals of hypergeometric Levy processes

We study the distribution and various properties of exponential functionals of hypergeometric Levy processes. We derive an explicit formula for the Mellin transform of the exponential functional and give both convergent and asymptotic series expansions of its probability density function. As applications we present a new proof of some of the results on the density of the supremum of a stable process, which were recently obtained by Kuznetsov "On extrema of stable processes" (2010) and Hubalek and Kuznetsov "A convergent series representation for the density of the supremum of a stable process" (2010). We also derive some new results related to (i) the entrance law of the stable process conditioned to stay positive, (ii) the entrance law of the excursion measure of the stable process reflected at its past infimum and (iii) the entrance law and the last passage time of the radial part of n-dimensional symmetric stable process.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Series representations and asymptotic expansions for the density of the supremum of a stable process

We derive explicit asymptotic expansions of the density of the supremum of a strictly stable process when the index $α$ is not rational. In the case when parameters $α$ and $ρ=\p(X_1>0)$ satisfy $ρ+k=l/α$ for some integers $k,l \ge 1$ we prove that these asymptotic expansions are in fact convergent series representations of the density of supremum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Wiener-Hopf factorization and distribution of extrema for a family of Lévy processes

In this paper we introduce a ten-parameter family of Lévy processes for which we obtain Wiener-Hopf factors and distribution of the supremum process in semi-explicit form. This family allows an arbitrary behavior of small jumps and includes processes similar to the generalized tempered stable, KoBoL and CGMY processes. Analytically it is characterized by the property that the characteristic exponent is a meromorphic function, expressed in terms of beta and digamma functions. We prove that the Wiener-Hopf factors can be expressed as infinite products over roots of a certain transcendental equation, and the density of the supremum process can be computed as an exponentially converging infinite series. In several special cases when the roots can be found analytically, we are able to identify the Wiener-Hopf factors and distribution of the supremum in closed form. In the general case we prove that all the roots are real and simple, and we provide localization results and asymptotic formulas which allow an efficient numerical evaluation. We also derive a convergence acceleration algorithm for infinite products and a simple and efficient procedure to compute the Wiener-Hopf factors for complex values of parameters. As a numerical example we discuss computation of the density of the supremum process.

0 citations0 reviewsNo code linkedOpen access