# Probability (math.PR)

• We propose a deep generative Markov State Model (DeepGenMSM) learning framework for inference of metastable dynamical systems and prediction of trajectories. After unsupervised training on time series data, the model contains (i) a probabilistic encoder that maps from high-dimensional configuration space to a small-sized vector indicating the membership to metastable (long-lived) states, (ii) a Markov chain that governs the transitions between metastable states and facilitates analysis of the long-time dynamics, and (iii) a generative part that samples the conditional distribution of configurations in the next time step. The model can be operated in a recursive fashion to generate trajectories to predict the system evolution from a defined starting state and propose new configurations. The DeepGenMSM is demonstrated to provide accurate estimates of the long-time kinetics and generate valid distributions for molecular dynamics (MD) benchmark systems. Remarkably, we show that DeepGenMSMs are able to make long time-steps in molecular configuration space and generate physically realistic structures in regions that were not seen in training data.
• We construct a statistical indicator for the detection of short-term asset price bubbles based on the information content of bid and ask market quotes for plain vanilla put and call options. Our construction makes use of the martingale theory of asset price bubbles and the fact that such scenarios where the price for an asset exceeds its fundamental value can in principle be detected by analysis of the asymptotic behavior of the implied volatility surface. For extrapolating this implied volatility, we choose the SABR model, mainly because of its decent fit to real option market quotes for a broad range of maturities and its ease of calibration. As main theoretical result, we show that under lognormal SABR dynamics, we can compute a simple yet powerful closed-form martingale defect indicator by solving an ill-posed inverse calibration problem. In order to cope with the ill-posedness and to quantify the uncertainty which is inherent to such an indicator, we adopt a Bayesian statistical parameter estimation perspective. This allows us to incorporate any available prior knowledge in a fully explicit way. We probe the resulting posterior densities with a combination of optimization and adaptive Markov chain Monte Carlo methods, thus providing a full-blown uncertainty estimation of all the underlying parameters and the martingale defect indicator. Finally, we provide a real-market test of the proposed option-based indicator.
• Firstly, we shall introduce the so-called snapping out Walsh's Brownian motion and present its relation with Walsh's Brownian motion. Then the stiff problem related to Walsh's Brownian motion will be described and we shall build a phase transition for it. The snapping out Walsh's Brownian motion corresponds to the so-called semi-permeable pattern of this stiff problem.
• Following recent developments in the application of Stein's method in queueing theory, this paper is intended to be a short treatment showing how Stein's method can be developed and applied to the single server queue in heavy traffic. Here we provide two approaches to this approximation: one based on equilibrium couplings and another involving comparison of generators.
• We propose a new explicit numerical scheme for stochastic differential equation with super-linearly growing drift and linearly growing diffusion coefficients which are also twice continuously differentiable. The rate of strong convergence in $\mathcal{L}^p$-norm is shown to be equal to $1.5$. Moreover, the scheme is computationally more efficient that the corresponding scheme available in the literature.
• Consider an intersection measure $\ell_t ^{\mathrm{IS}}$ of $p$ independent (possibly different) $m$-symmetric Hunt processes up to time $t$ in a metric measure space $E$ with a Radon measure $m$. We derive a Donsker-Varadhan type large deviation principle for the normalized intersection measure $t^{-p}\ell_t ^{\mathrm{IS}}$ on the set of finite measures on $E$ as $t \rightarrow \infty$, under the condition that $t$ is smaller than life times of all processes. This extends earlier work by W. König and C. Mukherjee (2013), in which the large deviation principle was established for the intersection measure of $p$ independent $N$-dimensional Brownian motions before exiting some bounded open set $D \subset \mathbb{R}^N$. We also obtain the asymptotic behaviour of logarithmic moment generating function, which is related to the results of X. Chen and J. Rosen (2005) on the intersection measure of independent Brownian motions or stable processes. Our results rely on assumptions about the heat kernels and the 1-order resolvents of the processes, hence include rich examples. For example, the assumptions hold for $p\in \mathbb{Z}$ with $2\leq p < p_*$ when the processes enjoy (sub-)Gaussian type or jump type heat kernel estimates, where $p_*$ is determined by the Hausdorff dimension of $E$ and the so-called walk dimensions of the processes.
• We study the mixing time of the $(n,k)$ Bernoulli--Laplace urn model, where $k\in\{0,1,\ldots,n\}$. Consider two urns, each containing $n$ balls, so that when combined they have precisely $n$ red balls and $n$ white balls. At each step of the process choose uniformly at random $k$ balls from the left urn and $k$ balls from the right urn and switch them simultaneously. We show that if $k=o(n)$, this Markov chain exhibits mixing time cutoff at $\frac{n}{4k}\log n$ and window of the order $\frac{n}{k}\log\log n$. This is an extension of a classical theorem of Diaconis and Shahshahani who treated the case $k=1$.
• Dunkl processes are multidimensional Markov processes defined through the use of Dunkl operators. These processes have discontinuities, and they can be separated into their continuous (radial) part, and their discontinuous (jump) part. While radial Dunkl processes have been studied thoroughly due to their relationship to families of stochastic particle systems such as the Dyson model and Wishart-Laguerre processes, Dunkl jump processes have gone largely unnoticed after the initial work of Gallardo, Yor and Chybiryakov. We study the dynamical properties of these processes, and we derive their master equation. By calculating the asymptotic behavior of their total jump rate, we find that the jump processes of types $A_{N-1}$ and $B_N$ undergo a phase transition when the parameter $\beta$ decreases toward one in the bulk scaling limit. In addition, we show that the relaxation behavior of these processes is given by a power law, and formulate a conjecture for the jump rate asymptotics based on numerical simulations.
• We show how an improvement on Doob's inequality leads to new inequalities for several stochastic processes: Lévy processes, processes with independent increments in general, branching processes, and time homogeneous Markov processes.
• We consider a class of multitype Galton-Watson branching processes with a countably infinite type set $\mathcal{X}_d$ whose mean progeny matrices have a block lower Hessenberg form. For these processes, the probability $\boldsymbol{q}(A)$ of extinction in subsets of types $A\subseteq \mathcal{X}_d$ may differ from the global extinction probability $\boldsymbol{q}$ and the partial extinction probability $\tilde{\boldsymbol{q}}$. After deriving partial and global extinction criteria, we develop conditions for $\boldsymbol{q}<\boldsymbol{q}(A)<\tilde{\boldsymbol{q}}$. We then present an iterative method to compute the vector $\boldsymbol{q}(A)$ for any set $A$. Finally, we investigate the location of the vectors $\boldsymbol{q}(A)$ in the set of fixed points of the progeny generating vector.
• We prove existence and uniqueness of strong solutions for a class of semilinear stochastic evolution equations driven by general Hilbert space-valued semimartingales, with drift equal to the sum of a linear maximal monotone operator in variational form and of the superposition operator associated to a random time-dependent monotone function defined on the whole real line. Such a function is only assumed to satisfy a very mild symmetry-like condition, but its rate of growth towards infinity can be arbitrary. Moreover, the noise is of multiplicative type and can be path-dependent. The solution is obtained via a priori estimates on solutions to regularized equations, interpreted both as stochastic equations as well as deterministic equations with random coefficients, and ensuing compactness properties. A key role is played by an infinite-dimensional Doob-type inequality due to Métivier and Pellaumail.
• This paper presents an algorithm to simulate Gaussian random vectors whose precision matrix can be expressed as a polynomial of a sparse matrix. This situation arises in particular when simulating Gaussian Markov random fields obtained by the discretization by finite elements of the solutions of some stochastic partial derivative equations. The proposed algorithm uses a Chebyshev polynomial approximation to compute simulated vectors with a linear complexity. This method is asymptotically exact as the approximation order grows. Criteria based on tests of the statistical properties of the produced vectors are derived to determine minimal orders of approximation.

Danial Dervovic Mar 01 2018 12:08 UTC

Hello again Māris, many thanks for your patience. Your comments and questions have given me much food for thought, and scope for an amended version of the paper -- please see my responses below.

Please if any of the authors of [AST17 [arXiv:1712.01609](https://arxiv.org/abs/1712.01609)] have any fu

...(continued)
Danial Dervovic Feb 05 2018 15:03 UTC

Thank you Māris for the extremely well thought-out and articulated points here.

I think this very clearly highlights the need to think explicitly about the precompute time if using the lifting to directly simulate the quantum walk, amongst other things.

I wish to give a well-considered respons

...(continued)
Māris Ozols Feb 01 2018 17:53 UTC

This paper considers the problem of using "lifted" Markov chains to simulate the mixing of coined quantum walks. The Markov chain has to approximately (in the total variational distance) sample from the distribution obtained by running the quantum walk for a randomly chosen time $t \in [0,T]$ follow

...(continued)
Danial Dervovic Dec 10 2017 15:25 UTC

Thank you for the insightful observations, Simon.

In response to the first point, there is a very short comment in the Discussion section to this effect. I felt an explicit dependence on $T$ as opposed to the diameter would make the implications of the result more clear. Namely, lifting can mix

...(continued)
Simon Apers Dec 09 2017 07:54 UTC

Thanks for the comment, Simone. A couple of observations:

- We noticed that Danial's result can in fact be proved more directly using the theorem that is used from ([arXiv:1705.08253][1]): by choosing the quantum walk Cesaro average as the goal distribution, it can be attained with a lifted Markov

...(continued)
Simone Severini Dec 07 2017 02:51 UTC

Closely related to

Simon Apers, Alain Sarlette, Francesco Ticozzi, Simulation of Quantum Walks and Fast Mixing with Classical Processes, https://scirate.com/arxiv/1712.01609

In my opinion, lifting is a good opportunity to put on a rigorous footing the relationship between classical and quantu

...(continued)
Mile Gu Nov 20 2015 05:04 UTC

Good question! There shouldn't be any contradiction with the correspondence principle. The reason here is that the quantum models are built to simulate the output behaviour of macroscopic, classical systems, and are not necessarily macroscopic themselves. When we compare quantum and classical comple

...(continued)
hong Nov 20 2015 00:40 UTC

Interesting results. But, just wondering, does it contradict to the correspondence principle?

Richard Kueng Jul 28 2015 07:01 UTC

fyi: our quantum implications are presented in Subsection 2.2 (pp 7-9).

Zoltán Zimborás May 28 2014 04:42 UTC

It's a bit funny to look at a formally verified proof of the CLT :), here it is online: