We formally introduce and study locally-balanced Markov jump processes (LBMJPs) defined on a general state space. These continuous-time stochastic processes with a user-specified limiting distribution are designed for sampling in settings involving discrete parameters and/or non-smooth distributions, addressing limitations of other processes such as the overdamped Langevin diffusion. The paper establishes the well-posedness, non-explosivity, and ergodicity of LBMJPs under mild conditions. We further explore regularity properties such as the Feller property and characterise the weak generator of the process. We then derive conditions for exponential ergodicity via spectral gaps and establish comparison theorems for different balancing functions. In particular we show an equivalence between the spectral gaps of Metropolis--Hastings algorithms and LBMJPs with bounded balancing function, but show that LBMJPs can exhibit uniform ergodicity on unbounded state spaces when the balancing function is unbounded, even when the limiting distribution is not sub-Gaussian. We also establish a diffusion limit for an LBMJP in the small jump limit, and discuss applications to Monte Carlo sampling and
This paper establishes the precise small-time asymptotic behavior of the spectral heat content for isotropic Lévy processes on bounded $C^{1,1}$ open sets of $\mathbb{R}^{d}$ with $d\ge 2$, where the underlying characteristic exponents are regularly varying at infinity with index $α\in (1,2]$, including the case $α=2$. Moreover, this asymptotic behavior is shown to be stable under an integrable perturbation of its Lévy measure. These results cover a wide class of isotropic Lévy processes, including Brownian motions, stable processes, and jump diffusions, and the proofs provide a unified approach to the asymptotic behavior of the spectral heat content for all of these processes.
Arcade processes are a class of continuous stochastic processes that interpolate in a strong sense, i.e., omega by omega, between zeros at fixed pre-specified times. Their additive randomisation allows one to match any finite sequence of target random variables, indexed by the given fixed dates, on the whole probability space. The randomised arcade processes (RAPs) can thus be interpreted as a generalisation of anticipative stochastic bridges. The filtrations generated by these processes are utilised to construct a class of martingales that interpolate between the given target random variables. These so-called filtered arcade martingales (FAMs) are almost-sure solutions to the martingale interpolation problem and reveal an underlying stochastic filtering structure. In the special case of conditionally Markov randomised arcade processes, the dynamics of FAMs are informed by Bayesian updating. The same ideas are applied to filtered arcade reverse-martingales, which are constructed in a similar fashion, using reverse-filtrations of RAPs, instead. Several explicit examples for RAPs and FAMs are provided and simulated. This paper concludes with an outlook on potential connections betwee
We introduce a class of hybrid marked point processes, which encompasses and extends continuous-time Markov chains and Hawkes processes. While this flexible class amalgamates such existing processes, it also contains novel processes with complex dynamics. These processes are defined implicitly via their intensity and are endowed with a state process that interacts with past-dependent events. The key example we entertain is an extension of a Hawkes process, a state-dependent Hawkes process interacting with its state process. We show the existence and uniqueness of hybrid marked point processes under general assumptions, extending the results of Massoulié (1998) on interacting point processes.
We consider random walks associated with conductances on Delaunay triangulations, Gabriel graphs and skeletons of Voronoi tilings which are generated by point processes in $\mathbb{R}^d$. Under suitable assumptions on point processes and conductances, we show that, for almost any realization of the point process, these random walks are recurrent if $d=2$ and transient if $d\geq 3$. These results hold for a large variety of point processes including Poisson point processes, Mat{é}rn cluster and Mat{é}rn hardcore processes which have clustering or repulsive properties. In order to prove them, we state general criteria for recurrence or almost sure transience which apply to random graphs embedded in $\mathbb{R}^d$.
Superpositions of Ornstein-Uhlenbeck type (supOU) processes provide a rich class of stationary stochastic processes for which the marginal distribution and the dependence structure may be modeled independently. We show that they can also display intermittency, a phenomenon affecting the rate of growth of moments. To do so, we investigate the limiting behavior of integrated supOU processes with finite variance. After suitable normalization four different limiting processes may arise depending on the decay of the correlation function and on the characteristic triplet of the marginal distribution. To show that supOU processes may exhibit intermittency, we establish the rate of growth of moments for each of the four limiting scenarios. The rate change indicates that there is intermittency, which is expressed here as a change-point in the asymptotic behavior of the absolute moments.
We derive a generalised Itō formula for stochastic processes which are constructed by a convolution of a deterministic kernel with a centred Lévy process. This formula has a unifying character in the sense that it contains the classical Itō formula for Lévy processes as well as recent change-of-variable formulas for Gaussian processes such as fractional Brownian motion as special cases. Our result also covers fractional Lévy processes (with Mandelbrot-Van Ness kernel) and a wide class of related processes for which such a generalised Itō formula has not yet been available in the literature.
In this paper, we study the memory properties of transformations of linear processes. Dittmann and Granger (2002) studied the polynomial transformations of Gaussian FARIMA(0,d,0) processes by applying the orthonormality of the Hermite polynomials under the measure for the standard normal distribution. Nevertheless, the orthogonality does not hold for transformations of non-Gaussian linear processes. Instead, we use the decomposition developed by Ho and Hsing (1996, 1997) to study the memory properties of nonlinear transformations of linear processes, which include the FARIMA(p,d,q) processes, and obtain consistent results as in the Gaussian case. In particular, for stationary processes, the transformations of short-memory time series still have short-memory and the transformation of long-memory time series may have different weaker memory parameters which depend on the power rank of the transformation. On the other hand, the memory properties of transformations of non-stationary time series may not depend on the power ranks of the transformations. This study has application in econometrics and financial data analysis when the time series observations have non-Gaussian heavy tails.
Chi-square processes with trend appear naturally as limiting processes in various statistical models. In this paper we are concerned with the exact tail asymptotics of the supremum taken over (0; 1) of a class of locally stationary chi-square processes with particular admissible trends. An important tool for establishing our results is a weak version of Slepian's lemma for chi-square processes. Some special cases including squared Brownian bridge and Bessel process are discussed.
For an arbitrary diffusion process $X$ with time-homogeneous drift and variance parameters $μ(x)$ and $σ^2(x)$, let $V_\varepsilon$ be $1/\varepsilon$ times the total time $X(t)$ spends in the strip $[a+bt-(1/2)\varepsilon,a+bt+(1/2)\varepsilon]$.The limit $V$ as $\varepsilon\rightarrow0$ is the full halfline version of the local time of $X(t)-a-bt$ at zero, and can be thought of as the time $X$ spends along the straight line $x=a+bt$. We prove that $V$ is either infinite with probability 1 or distributed as a mixture of an exponential and a unit point mass at zero, and we give formulae for the parameters of this distribution in terms of $μ(\cdot)$, $σ(\cdot)$, $a$, $b$, and the starting point $X(0)$. The special case ofa Brownian motion is studied in more detail, leading in particular to a full process $V(b)$ with continuous sample paths and exponentially distributed marginals. This construction leads to new families of bivariate and multivariate exponential distributions. Truncated versions of such `total relative time' variables are also studied. A relation is pointed out to a second order asymptotics problem in statistical estimation theory, recently investigated in Hjort and F
We propose a novel fused Gromov-Wasserstein alignment method to jointly learn the Hawkes processes in different event spaces, and align their event types. Given two Hawkes processes, we use fused Gromov-Wasserstein discrepancy to measure their dissimilarity, which considers both the Wasserstein discrepancy based on their base intensities and the Gromov-Wasserstein discrepancy based on their infectivity matrices. Accordingly, the learned optimal transport reflects the correspondence between the event types of these two Hawkes processes. The Hawkes processes and their optimal transport are learned jointly via maximum likelihood estimation, with a fused Gromov-Wasserstein regularizer. Experimental results show that the proposed method works well on synthetic and real-world data.
A new class of formal latent-variable stochastic processes called hidden quantum models (HQM's) is defined in order to clarify the theoretical foundations of ion channel signal processing. HQM's are based on quantum stochastic processes which formalize time-dependent observation. They allow the calculation of autocovariance functions which are essential for frequency-domain signal processing. HQM's based on a particular type of observation protocol called independent activated measurements are shown to to be distributionally equivalent to hidden Markov models yet without an underlying physical Markov process. Since the formal Markov processes are non-physical, the theory of activated measurement allows merging energy-based Eyring rate theories of ion channel behavior with the more common phenomenological Markov kinetic schemes to form energy-modulated quantum channels. Using the simplest quantum channel model consistent with neuronal membrane voltage-clamp experiments, activation eigenenergies are calculated for the Hodgkin-Huxley K+ and Na+ ion channels. It is also shown that maximizing entropy under constrained activation energy yields noise spectral densities approximating $S(f)
We characterize the finite variation property for stationary increment mixed moving averages driven by infinitely divisible random measures. Such processes include fractional and moving average processes driven by Levy processes, and also their mixtures. We establish two types of zero-one laws for the finite variation property. We also consider some examples to illustrate our results.
In this paper, we consider a long-time behavior of stable-like processes. A stable-like process is a Feller process given by the symbol $p(x,ξ)=-iβ(x)ξ+γ(x)|ξ|^{α(x)},$ where $α(x)\in(0,2)$, $β(x)\in\R$ and $γ(x)\in(0,\infty)$. More precisely, we give sufficient conditions for recurrence, transience and ergodicity of stable-like processes in terms of the stability function $α(x)$, the drift function $β(x)$ and the scaling function $γ(x)$. Further, as a special case of these results we give a new proof for the recurrence and transience property of one-dimensional symmetric stable Lévy processes with the index of stability $α eq1.$
Nonlinear phenomena are essential ingredients in many oceanic and atmospheric processes, and successful understanding of them benefits from multidisciplinary collaboration between oceanographers, meteorologists, physicists and mathematicians. The present Special Issue on ``Nonlinear Processes in Oceanic and Atmospheric Flows'' contains selected contributions from attendants to the workshop which, in the above spirit, was held in Castro Urdiales, Spain, in July 2008. Here we summarize the Special Issue contributions, which include papers on the characterization of ocean transport in the Lagrangian and in the Eulerian frameworks, generation and variability of jets and waves, interactions of fluid flow with plankton dynamics or heavy drops, scaling in meteorological fields, and statistical properties of El Niño Southern Oscillation.
The motivation of this paper is to prove verification theorems for stochastic optimal control of finite dimensional diffusion processes without control in the diffusion term, in the case that the value function is assumed to be continuous in time and once differentiable in the space variable ($C^{0,1}$) instead of once differentiable in time and twice in space ($C^{1,2}$), like in the classical results. For this purpose, the replacement tool of the Itô formula will be the Fukushima-Dirichlet decomposition for weak Dirichlet processes. Given a fixed filtration, a weak Dirichlet process is the sum of a local martingale $M$ plus an adapted process $A$ which is orthogonal, in the sense of covariation, to any continuous local martingale. The mentioned decomposition states that a $C^{0,1}$ function of a weak Dirichlet process with finite quadratic variation is again a weak Dirichlet process. That result is established in this paper and it is applied to the strong solution of a Cauchy problem with final condition. Applications to the proof of verification theorems will be addressed in a companion paper.
The detection and characterization of potentially habitable exoplanets is one of the chief goals of astrophysics for the coming decades. Imaging in reflected light is well suited for characterizing Earth-like planets, as much can be learned about these planets in this wavelength range (i.e., ~0.3-2 μm). Several studies have been conducted to determine the abilities and limitations of reflectance spectroscopy, but most previous studies assumed a homogeneous atmospheric and surface composition. Here we investigate how heterogeneities in the atmosphere and surface of an Earth-like planet impact retrieval results. We extend the ExoReL retrieval framework to include a step function for retrieving wavelength varying surface albedo. We then use it to retrieve on visible-to-near-infrared spectra of realistic 3D Earth models with different surface features in view and varying cloud types/distributions synthesized with the Planetary Spectrum Generator. Including the ability to fit for wavelength dependent albedo mitigates degeneracies that arise when using 1D models to analyze 3D planets, and we recover an Earth-like planet in all cases. We detect surface albedo steps at ~0.7 and ~1.1 μm des
We introduce the title process via a particular construction, and relate it to processes previously studied, in particular a process introduced by G. E. H. Reuter in 1969. We derive elementary properties and quantities of this processes: Markov property, transition rates, stationary distribution, and the infinitesimal generator for a case not treated by Reuter.
We obtain large deviation results for a two time-scale model of jump-diffusion processes. The processes on the two time scales are fully inter-dependent, the slow process has small perturbative noise and the fast process is ergodic. Our results extend previous large deviation results for diffusions. We provide concrete examples in their applications to finance and biology, with an explicit calculation of the large deviation rate function.
The Malliavin integration-by-parts formula is a key ingredient to develop stochastic analysis on the Wiener space. In this article we show that a suitable integration-by-parts formula also characterizes a wide class of Gaussian processes, the so-called Gaussian Fredholm processes.