We prove some estimates for the variations of transition probabilities of the (1+1)-affine process. From these estimates we deduce the strong Feller and the ergodic properties of the total variation distance of the process. The key strategy is the pathwise construction and analysis of several Markov couplings using strong solutions of stochastic equations.

We study a stochastic differential equation with an unbounded drift and general Hölder continuous noise of order $\lambda \in (0,1)$ . The corresponding equation turns out to have a unique solution that, depending on a particular shape of the drift, either stays above some continuous function or has continuous upper and lower bounds. Under some mild assumptions on the noise, we prove that the solution has moments of all orders. In addition, we provide its connection to the solution of some Skorokhod reflection problem. As an illustration of our results and motivation for applications, we also suggest two stochastic volatility models which we regard as generalizations of the CIR and CEV processes. We complete the study by providing a numerical scheme for the solution.

We consider solutions of Lévy-driven stochastic differential equations of the form $\textrm{d} X_t=\sigma(X_{t-})\textrm{d} L_t$ , $X_0=x$ , where the function $\sigma$ is twice continuously differentiable and the driving Lévy process $L=(L_t)_{t\geq0}$ is either vector or matrix valued. While the almost sure short-time behavior of Lévy processes is well known and can be characterized in terms of the characteristic triplet, there is no complete characterization of the behavior of the solution X. Using methods from stochastic calculus, we derive limiting results for stochastic integrals of the form $t^{-p}\int_{0+}^t\sigma(X_{t-})\,\textrm{d} L_t$ to show that the behavior of the quantity $t^{-p}(X_t-X_0)$ for $t\downarrow0$ almost surely reflects the behavior of $t^{-p}L_t$ . Generalizing $t^{{\kern1pt}p}$ to a suitable function $f\colon[0,\infty)\rightarrow\mathbb{R}$ then yields a tool to derive explicit law of the iterated logarithm type results for the solution from the behavior of the driving Lévy process.

We prove existence and uniqueness for the solution of a class of mixed fractional stochastic differential equations with discontinuous drift driven by both standard and fractional Brownian motion. Additionally, we establish a generalized Itô rule valid for functions with an absolutely continuous derivative and applicable to solutions of mixed fractional stochastic differential equations with Lipschitz coefficients, which plays a key role in our proof of existence and uniqueness. The proof of such a formula is new and relies on showing the existence of a density of the law under mild assumptions on the diffusion coefficient.

We study approximations for the Lévy area of Brownian motion which are based on the Fourier series expansion and a polynomial expansion of the associated Brownian bridge. Comparing the asymptotic convergence rates of the Lévy area approximations, we see that the approximation resulting from the polynomial expansion of the Brownian bridge is more accurate than the Kloeden–Platen–Wright approximation, whilst still only using independent normal random vectors. We then link the asymptotic convergence rates of these approximations to the limiting fluctuations for the corresponding series expansions of the Brownian bridge. Moreover, and of interest in its own right, the analysis we use to identify the fluctuation processes for the Karhunen–Loève and Fourier series expansions of the Brownian bridge is extended to give a stand-alone derivation of the values of the Riemann zeta function at even positive integers.

In this paper we study a class of optimal stopping problems under g-expectation, that is, the cost function is described by the solution of backward stochastic differential equations (BSDEs). Primarily, we assume that the reward process is $L\exp\bigl(\mu\sqrt{2\log\!(1+L)}\bigr)$ -integrable with $\mu>\mu_0$ for some critical value $\mu_0$ . This integrability is weaker than $L^p$ -integrability for any $p>1$ , so it covers a comparatively wide class of optimal stopping problems. To reach our goal, we introduce a class of reflected backward stochastic differential equations (RBSDEs) with $L\exp\bigl(\mu\sqrt{2\log\!(1+L)}\bigr)$ -integrable parameters. We prove the existence, uniqueness, and comparison theorem for these RBSDEs under Lipschitz-type assumptions on the coefficients. This allows us to characterize the value function of our optimal stopping problem as the unique solution of such RBSDEs.

In this paper, we study the optimal multiple stopping problem under the filtration-consistent nonlinear expectations. The reward is given by a set of random variables satisfying some appropriate assumptions, rather than a process that is right-continuous with left limits. We first construct the optimal stopping time for the single stopping problem, which is no longer given by the first hitting time of processes. We then prove by induction that the value function of the multiple stopping problem can be interpreted as the one for the single stopping problem associated with a new reward family, which allows us to construct the optimal multiple stopping times. If the reward family satisfies some strong regularity conditions, we show that the reward family and the value functions can be aggregated by some progressive processes. Hence, the optimal stopping times can be represented as hitting times.

In the current work, we study a stochastic parabolic problem. The presented problem is motivated by the study of an idealised electrically actuated MEMS (Micro-Electro-Mechanical System) device in the case of random fluctuations of the potential difference, a parameter that actually controls the operation of MEMS device. We first present the construction of the mathematical model, and then, we deduce some local existence results. Next for some particular versions of the model, relevant to various boundary conditions, we derive quenching results as well as estimations of the probability for such singularity to occur. Additional numerical study of the problem in one dimension follows, which also allows the further investigation the problem with respect to its quenching behaviour.

We prove a rate of convergence for the N-particle approximation of a second-order partial differential equation in the space of probability measures, such as the master equation or Bellman equation of the mean-field control problem under common noise. The rate is of order $1/N$ for the pathwise error on the solution v and of order $1/\sqrt{N}$ for the $L^2$ -error on its L-derivative $\partial_\mu v$ . The proof relies on backward stochastic differential equation techniques.

A system of mutually interacting superprocesses with migration is constructed as the limit of a sequence of branching particle systems arising from population models. The uniqueness in law of the superprocesses is established using the pathwise uniqueness of a system of stochastic partial differential equations, which is satisfied by the corresponding system of distribution function-valued processes.

Motivated by applications to wireless networks, cloud computing, data centers, etc., stochastic processing networks have been studied in the literature under various asymptotic regimes. In the heavy traffic regime, the steady-state mean queue length is proved to be $\Theta({1}/{\epsilon})$ , where $\epsilon$ is the heavy traffic parameter (which goes to zero in the limit). The focus of this paper is on obtaining queue length bounds on pre-limit systems, thus establishing the rate of convergence to heavy traffic. For the generalized switch, operating under the MaxWeight algorithm, we show that the mean queue length is within $\textrm{O}({\log}({1}/{\epsilon}))$ of its heavy traffic limit. This result holds regardless of the complete resource pooling (CRP) condition being satisfied. Furthermore, when the CRP condition is satisfied, we show that the mean queue length under the MaxWeight algorithm is within $\textrm{O}({\log}({1}/{\epsilon}))$ of the optimal scheduling policy. Finally, we obtain similar results for the rate of convergence to heavy traffic of the total queue length in load balancing systems operating under the ‘join the shortest queue’ routeing algorithm.

In this paper we propose a general framework for modeling an insurance liability cash flow in continuous time, by generalizing the reduced-form framework for credit risk and life insurance. In particular, we assume a nontrivial dependence structure between the reference filtration and the insurance internal filtration. We apply these results for pricing and hedging non-life insurance liabilities in hybrid financial and insurance markets, while taking into account the role of inflation under the benchmarked risk-minimization approach. This framework offers at the same time a general and flexible structure, and an explicit and treatable pricing-hedging formula.

We study an intertemporal consumption and portfolio choice problem under Knightian uncertainty in which agent’s preferences exhibit local intertemporal substitution. We also allow for market frictions in the sense that the pricing functional is nonlinear. We prove existence and uniqueness of the optimal consumption plan, and we derive a set of sufficient first-order conditions for optimality. With the help of a backward equation, we are able to determine the structure of optimal consumption plans. We obtain explicit solutions in a stationary setting in which the financial market has different risk premia for short and long positions.

In this paper, we construct operator fractional Lévy motion (ofLm), a broad class of infinitely divisible stochastic processes that are covariance operator self-similar and have wide-sense stationary increments. The ofLm class generalizes the univariate fractional Lévy motion as well as the multivariate operator fractional Brownian motion (ofBm). OfLm can be divided into two types, namely, moving average (maofLm) and real harmonizable (rhofLm), both of which share the covariance structure of ofBm under assumptions. We show that maofLm and rhofLm admit stochastic integral representations in the time and Fourier domains, and establish their distinct small- and large-scale limiting behavior. We also characterize time-reversibility for ofLm through parametric conditions related to its Lévy measure. In particular, we show that, under non-Gaussianity, the parametric conditions for time-reversibility are generally more restrictive than those for the Gaussian case (ofBm).

We provide constructions of age-structured branching processes without or with immigration as pathwise-unique solutions to stochastic integral equations. A necessary and sufficient condition for the ergodicity of the model with immigration is also given.

This paper discusses a general class of replicator–mutator equations on a multidimensional fitness space. We establish a novel probabilistic representation of weak solutions of the equation by using the theory of Fokker–Planck–Kolmogorov (FPK) equations and a martingale extraction approach. We provide examples with closed-form probabilistic solutions for different fitness functions considered in the existing literature. We also construct a particle system and prove a general convergence result to the unique solution of the FPK equation associated with the extended replicator–mutator equation with respect to a Wasserstein-like metric adapted to our probabilistic framework.

In this paper we employ a Gaussian-type heat kernel estimate to establish Krylov’s estimate and Khasminskii’s estimate for the Euler–Maruyama (EM) algorithm. For applications, by taking Zvonkin’s transformation into account, we investigate the convergence rate of the EM algorithm for a class of multidimensional stochastic differential equations (SDEs) with low regular drifts, which need not be piecewise Lipschitz.

We obtain a generalisation of the Stroock–Varadhan support theorem for a large class of systems of subcritical singular stochastic partial differential equations driven by a noise that is either white or approximately self-similar. The main problem that we face is the presence of renormalisation. In particular, it may happen in general that different renormalisation procedures yield solutions with different supports. One of the main steps in our construction is the identification of a subgroup $\mathcal {H}$ of the renormalisation group such that any renormalisation procedure determines a unique coset $g\circ \mathcal {H}$ . The support of the solution then depends only on this coset and is obtained by taking the closure of all solutions obtained by replacing the driving noises by smooth functions in the equation that is renormalised by some element of $g\circ \mathcal {H}$ .

One immediate corollary of our results is that the $\Phi ^4_3$ measure in finite volume has full support, and the associated Langevin dynamic is exponentially ergodic.

In this paper an exact rejection algorithm for simulating paths of the coupled Wright–Fisher diffusion is introduced. The coupled Wright–Fisher diffusion is a family of multivariate Wright–Fisher diffusions that have drifts depending on each other through a coupling term and that find applications in the study of networks of interacting genes. The proposed rejection algorithm uses independent neutral Wright–Fisher diffusions as candidate proposals, which are only needed at a finite number of points. Once a candidate is accepted, the remainder of the path can be recovered by sampling from neutral multivariate Wright–Fisher bridges, for which an exact sampling strategy is also provided. Finally, the algorithm’s complexity is derived and its performance demonstrated in a simulation study.

Taking the consideration of two-dimensional stochastic Navier–Stokes equations with multiplicative Lévy noises, where the noises intensities are related to the viscosity, a large deviation principle is established by using the weak convergence method skillfully, when the viscosity converges to 0. Due to the appearance of the jumps, it is difficult to close the energy estimates and obtain the desired convergence. Hence, one cannot simply use the weak convergence approach. To overcome the difficulty, one introduces special norms for new arguments and more careful analysis.