We investigate the European call option pricing problem under the fractional stochastic volatility model. The stochastic volatility model is driven by both fractional Brownian motion and standard Brownian motion. We obtain an analytical solution of the European option price via the Itô’s formula for fractional Brownian motion, Malliavin calculus, derivative replication and the fundamental solution method. Some numerical simulations are given to illustrate the impact of parameters on option prices, and the results of comparison with other models are presented.

Mixing rates, relaxation rates, and decay of correlations for dynamics defined by potentials with summable variations are well understood, but little is known for non-summable variations. This paper exhibits upper bounds for these quantities for dynamics defined by potentials with square-summable variations. We obtain these bounds as corollaries of a new block coupling inequality between pairs of dynamics starting with different histories. As applications of our results, we prove a new weak invariance principle and a Hoeffding-type inequality.

Cogdell et al. [‘Evaluating the Mahler measure of linear forms via Kronecker limit formulas on complex projective space’, Trans. Amer. Math. Soc. (2021), to appear] developed infinite series representations for the logarithmic Mahler measure of a complex linear form with four or more variables. We establish the case of three variables by bounding an integral with integrand involving the random walk probability density $a\int _0^\infty tJ_0(at) \prod _{m=0}^2 J_0(r_m t)\,dt$, where $J_0$ is the order-zero Bessel function of the first kind and a and $r_m$ are positive real numbers. To facilitate our proof we develop an alternative description of the integral’s asymptotic behaviour at its known points of divergence. As a computational aid for numerical experiments, an algorithm to calculate these series is presented in the appendix.

In this work, we study a new model for continuum line-of-sight percolation in a random environment driven by the Poisson–Voronoi tessellation in the d-dimensional Euclidean space. The edges (one-dimensional facets, or simply 1-facets) of this tessellation are the support of a Cox point process, while the vertices (zero-dimensional facets or simply 0-facets) are the support of a Bernoulli point process. Taking the superposition Z of these two processes, two points of Z are linked by an edge if and only if they are sufficiently close and located on the same edge (1-facet) of the supporting tessellation. We study the percolation of the random graph arising from this construction and prove that a 0–1 law, a subcritical phase, and a supercritical phase exist under general assumptions. Our proofs are based on a coarse-graining argument with some notion of stabilization and asymptotic essential connectedness to investigate continuum percolation for Cox point processes. We also give numerical estimates of the critical parameters of the model in the planar case, where our model is intended to represent telecommunications networks in a random environment with obstructive conditions for signal propagation.

We derive the large-sample distribution of the number of species in a version of Kingman’s Poisson–Dirichlet model constructed from an $\alpha$-stable subordinator but with an underlying negative binomial process instead of a Poisson process. Thus it depends on parameters $\alpha\in (0,1)$ from the subordinator and $r>0$ from the negative binomial process. The large-sample distribution of the number of species is derived as sample size $n\to\infty$. An important component in the derivation is the introduction of a two-parameter version of the Dickman distribution, generalising the existing one-parameter version. Our analysis adds to the range of Poisson–Dirichlet-related distributions available for modeling purposes.

There are two types of tempered stable (TS) based Ornstein–Uhlenbeck (OU) processes: (i) the OU-TS process, the OU process driven by a TS subordinator, and (ii) the TS-OU process, the OU process with TS marginal law. They have various applications in financial engineering and econometrics. In the literature, only the second type under the stationary assumption has an exact simulation algorithm. In this paper we develop a unified approach to exactly simulate both types without the stationary assumption. It is mainly based on the distributional decomposition of stochastic processes with the aid of an acceptance–rejection scheme. As the inverse Gaussian distribution is an important special case of TS distribution, we also provide tailored algorithms for the corresponding OU processes. Numerical experiments and tests are reported to demonstrate the accuracy and effectiveness of our algorithms, and some further extensions are also discussed.

We focus on the population dynamics driven by two classes of truncated $\alpha$-stable processes with Markovian switching. Almost necessary and sufficient conditions for the ergodicity of the proposed models are provided. Also, these results illustrate the impact on ergodicity and extinct conditions as the parameter $\alpha$ tends to 2.

We investigate the impact of Knightian uncertainty on the optimal timing policy of an ambiguity-averse decision-maker in the case where the underlying factor dynamics follow a multidimensional Brownian motion and the exercise payoff depends on either a linear combination of the factors or the radial part of the driving factor dynamics. We present a general characterization of the value of the optimal timing policy and the worst-case measure in terms of a family of explicitly identified excessive functions generating an appropriate class of supermartingales. In line with previous findings based on linear diffusions, we find that ambiguity accelerates timing in comparison with the unambiguous setting. Somewhat surprisingly, we find that ambiguity may lead to stationarity in models which typically do not possess stationary behavior. In this way, our results indicate that ambiguity may act as a stabilizing mechanism.

For a determinantal point process (DPP) X with a kernel K whose spectrum is strictly less than one, André Goldman has established a coupling to its reduced Palm process $X^u$ at a point u with $K(u,u)>0$ so that, almost surely, $X^u$ is obtained by removing a finite number of points from X. We sharpen this result, assuming weaker conditions and establishing that $X^u$ can be obtained by removing at most one point from X, where we specify the distribution of the difference $\xi_u: = X\setminus X^u$. This is used to discuss the degree of repulsiveness in DPPs in terms of $\xi_u$, including Ginibre point processes and other specific parametric models for DPPs.

Rough volatility is a well-established statistical stylized fact of financial assets. This property has led to the design and analysis of various new rough stochastic volatility models. However, most of these developments have been carried out in the mono-asset case. In this work, we show that some specific multivariate rough volatility models arise naturally from microstructural properties of the joint dynamics of asset prices. To do so, we use Hawkes processes to build microscopic models that accurately reproduce high-frequency cross-asset interactions and investigate their long-term scaling limits. We emphasize the relevance of our approach by providing insights on the role of microscopic features such as momentum and mean-reversion in the multidimensional price formation process. In particular, we recover classical properties of high-dimensional stock correlation matrices.

We provide upper and lower bounds for the mean $\mathscr{M}(H)$ of $\sup_{t\geq 0} \{B_H(t) - t\}$, with $B_H(\!\cdot\!)$ a zero-mean, variance-normalized version of fractional Brownian motion with Hurst parameter $H\in(0,1)$. We find bounds in (semi-) closed form, distinguishing between $H\in(0,\frac{1}{2}]$ and $H\in[\frac{1}{2},1)$, where in the former regime a numerical procedure is presented that drastically reduces the upper bound. For $H\in(0,\frac{1}{2}]$, the ratio between the upper and lower bound is bounded, whereas for $H\in[\frac{1}{2},1)$ the derived upper and lower bound have a strongly similar shape. We also derive a new upper bound for the mean of $\sup_{t\in[0,1]} B_H(t)$, $H\in(0,\frac{1}{2}]$, which is tight around $H=\frac{1}{2}$.

Nowadays many financial derivatives, such as American or Bermudan options, are of early exercise type. Often the pricing of early exercise options gives rise to high-dimensional optimal stopping problems, since the dimension corresponds to the number of underlying assets. High-dimensional optimal stopping problems are, however, notoriously difficult to solve due to the well-known curse of dimensionality. In this work, we propose an algorithm for solving such problems, which is based on deep learning and computes, in the context of early exercise option pricing, both approximations of an optimal exercise strategy and the price of the considered option. The proposed algorithm can also be applied to optimal stopping problems that arise in other areas where the underlying stochastic process can be efficiently simulated. We present numerical results for a large number of example problems, which include the pricing of many high-dimensional American and Bermudan options, such as Bermudan max-call options in up to 5000 dimensions. Most of the obtained results are compared to reference values computed by exploiting the specific problem design or, where available, to reference values from the literature. These numerical results suggest that the proposed algorithm is highly effective in the case of many underlyings, in terms of both accuracy and speed.

We investigate various variable martingale Hardy spaces corresponding to variable Lebesgue spaces $\mathcal {L}_{p(\cdot )}$ defined by rearrangement functions. In particular, we show that the dual of martingale variable Hardy space $\mathcal {H}_{p(\cdot )}^{s}$ with $0<p_{-}\leq p_{+}\leq 1$ can be described as a BMO-type space and establish martingale inequalities among these martingale Hardy spaces. Furthermore, we give an application of martingale inequalities in stochastic integral with Brownian motion.

It is shown that for a dense $G_\delta $-subset of the subgroup of non-singular transformations (of a standard infinite $\sigma $-finite measure space) whose Poisson suspensions are non-singular, the corresponding Poisson suspensions are ergodic and of Krieger’s type III1.

In this paper we consider the pricing and hedging of financial derivatives in a model-independent setting, for a trader with additional information, or beliefs, on the evolution of asset prices. In particular, we suppose that the trader wants to act in a way which is independent of any modelling assumptions, but that she observes market information in the form of the prices of vanilla call options on the asset. We also assume that both the payoff of the derivative, and the insider’s information or beliefs, which take the form of a set of impossible paths, are time-invariant. In this way we accommodate drawdown constraints, as well as information/beliefs on quadratic variation or on the levels hit by asset prices. Our setup allows us to adapt recent work of [12] to prove duality results and a monotonicity principle. This enables us to determine geometric properties of the optimal models. Moreover, for specific types of information, we provide simple conditions for the existence of consistent models for the informed agent. Finally, we provide an example where our framework allows us to compute the impact of the information on the agent’s pricing bounds.

We consider the optimal prediction problem of stopping a spectrally negative Lévy process as close as possible to a given distance $b \geq 0$ from its ultimate supremum, under a squared-error penalty function. Under some mild conditions, the solution is fully and explicitly characterised in terms of scale functions. We find that the solution has an interesting non-trivial structure: if b is larger than a certain threshold then it is optimal to stop as soon as the difference between the running supremum and the position of the process exceeds a certain level (less than b), while if b is smaller than this threshold then it is optimal to stop immediately (independent of the running supremum and position of the process). We also present some examples.

We present closed-form solutions to some discounted optimal stopping problems for the running maximum of a geometric Brownian motion with payoffs switching according to the dynamics of a continuous-time Markov chain with two states. The proof is based on the reduction of the original problems to the equivalent free-boundary problems and the solution of the latter problems by means of the smooth-fit and normal-reflection conditions. We show that the optimal stopping boundaries are determined as the maximal solutions of the associated two-dimensional systems of first-order nonlinear ordinary differential equations. The obtained results are related to the valuation of real switching lookback options with fixed and floating sunk costs in the Black–Merton–Scholes model.

We consider a space-time random field on ${{\mathbb{R}^d} \times {\mathbb{R}}}$ given as an integral of a kernel function with respect to a Lévy basis with a convolution equivalent Lévy measure. The field obeys causality in time and is thereby not continuous along the time axis. For a large class of such random fields we study the tail behaviour of certain functionals of the field. It turns out that the tail is asymptotically equivalent to the right tail of the underlying Lévy measure. Particular examples are the asymptotic probability that there is a time point and a rotation of a spatial object with fixed radius, in which the field exceeds the level x, and that there is a time interval and a rotation of a spatial object with fixed radius, in which the average of the field exceeds the level x.

We present Lyapunov-type conditions for non-strong ergodicity of Markov processes. Some concrete models are discussed, including diffusion processes on Riemannian manifolds and Ornstein–Uhlenbeck processes driven by symmetric $\alpha$-stable processes. In particular, we show that any process of d-dimensional Ornstein–Uhlenbeck type driven by $\alpha$-stable noise is not strongly ergodic for every $\alpha\in (0,2]$.

Extended gamma processes have been seen as a flexible extension of standard gamma processes in the recent reliability literature, for the purpose of cumulative deterioration modeling. The probabilistic properties of the standard gamma process have been well explored since the 1970s, whereas those of its extension remain largely unexplored. In particular, stochastic comparisons between degradation levels modeled by standard gamma processes and ageing properties for the corresponding level-crossing times are now well understood. The aim of this paper is to explore similar properties for extended gamma processes and see which ones can be broadened to this new context. As a by-product, new stochastic comparisons for convolutions of gamma random variables are also obtained.