2024-10-22 | | Total: 31
We consider a degenerate infinite dimensional stochastic Hamiltonian system with multiplicative noise and establish the essential m-dissipativity on $L^2(\mu^{\Phi})$ of the corresponding Kolmogorov (backwards) operator. Here, $\Phi$ is the potential and $\mu^{\Phi}$ the invariant measure with density $e^{-\Phi}$ with respect to an infinite dimensional non-degenerate Gaussian measure. The main difficulty, besides the non-sectorality of the Kolmogorov operator, is the coverage of a large class of potentials. We include potentials that have neither a bounded nor a Lipschitz continuous gradient. The essential m-dissipativity is the starting point to establish the hypocoercivity of the strongly continuous contraction semigroup $(T_t)_{t\geq 0}$ generated by the Kolmogorov operator. By using the refined abstract Hilbert space hypocoercivity method of Grothaus and Stilgenbauer, originally introduced by Dolbeault, Mouhot and Schmeiser, we construct a $\mu^{\Phi}$-invariant Hunt process with weakly continuous paths and infinite lifetime, whose transition semigroup is associated with $(T_t)_{t\geq 0}$. This process provides a stochastically and analytically weak solution to the degenerate infinite dimensional stochastic Hamiltonian system with multiplicative noise. The hypocoercivity of $(T_t)_{t\geq 0}$ and the identification of $(T_t)_{t\geq 0}$ with the transition semigroup of the process leads to the exponential ergodicity. Finally, we apply our results to degenerate second order in time stochastic reaction-diffusion equations with multiplicative noise. A discussion of the class of applicable potentials and coefficients governing these equations completes our analysis.
In this note, we connect two seemingly unrelated objects: On the one hand is a two-dimensional drift-diffusion process $X$ with divergence-free and time-independent drift $b$. The drift is given by a stationary Gaussian ensemble, and we focus on the critical case where a small-scale cut-off is necessary for well-posedness and the large-scale cancellations lead to a borderline super-diffusive behavior. On the other hand is the natural diffusion $F$ on the Lie group $\textbf{SL}(2)$ of matrices of determinant one. As a consequence of this connection, the strongly non-Gaussian character of $F$ transmits to how $X$ depends on its starting point.
The object of this paper is to study and develop a Poisson distribution in generalized Wright function form.
We study a system of stochastic differential equations with singular drift which describes the dynamics of signed particles in two dimensions interacting by the Coulomb potential. In contrast to the well-studied cases of identical particles that either all repel each other or all attract each other, this system contains both `positive' and `negative' particles. Equal signs repel and opposite signs attract each other; apart from the sign, the potential is the same. We derive results on well-posedness of the system, on the type of collisions that can occur, and on the mean-field limit as the number of particles tends to infinity. Our results demonstrate that the signed system shares features of both the fully repulsive and the fully attractive cases. Our proof method is inspired by the work of Fournier and Jourdain (The Annals of Applied Probability, 27, pp. 2807-2861, 2017) on the fully attractive case; we construct an approximate system of equations, establish uniform estimates, and use tightness to pass to limits.
Examples of stochastic processes whose state space representations involve functions of an integral type structure $$I_{t}^{(a,b)}:=\int_{0}^{t}b(Y_{s})e^{-\int_{s}^{t}a(Y_{r})dr}ds, \quad t\ge 0$$ are studied under an ergodic semi-Markovian environment described by an $S$ valued jump type process $Y:=(Y_{s}:s\in\mathbb{R}^{+})$ that is ergodic with a limiting distribution $\pi\in\mathcal{P}(S)$. Under different assumptions on signs of $E_{\pi}a(\cdot):=\sum_{j\in S}\pi_{j}a(j)$ and tail properties of the sojourn times of $Y$ we obtain different long time limit results for $I^{(a,b)}_{}:=(I^{(a,b)}_{t}:t\ge 0).$ In all cases mixture type of laws emerge which are naturally represented through an affine stochastic recurrence equation (SRE) $X\stackrel{d}{=}AX+B,\,\, X\perp\!\!\!\perp (A, B)$. Examples include explicit long-time representations of pitchfork bifurcation, and regime-switching diffusions under semi-Markov modulated environments, etc.
We study two continuous and isotropic analogues of the model of greedy lattice animals introduced by Cox, Gandolfi, Griffin and Kesten in 1993. In our framework, animals collect masses scattered on a Poisson point process on $\mathbb R^d$, and are allowed to have vertices outside the process or not, depending on the model. The author recently proved in a more general setting that for all $u$ in the Euclidean open unit ball, the maximal mass of animals with length $\ell$, containing $0$ and $\ell u$ satisfies a law of large numbers. We prove some additional properties in the Poissonian case, including an extension of the functional law of large numbers to the closed unit ball, and study strict monotonicity of the limit function along a radius. Moreover, we prove that a third, penalized model is a suitable interpolation between the former two.
In this paper, we study a model of long-range site percolation on graphs of bounded degree, namely the Boolean percolation model. In this model, each vertex of an infinite connected graph is the center of a ball of random radius, and vertices are said to be active independently with probability $p \in [0, 1]$. We consider $W$ to be the reunion of random balls with an active center. In certain circumstances, the model does not exhibit a phase transition, in the sense that $W$ almost surely contains an infinite component for all $p > 0$, or even $W$ covers the entire graph. In this paper, we give a sufficient condition on the radius distribution for the existence of a subcritical phase, namely a regime such that almost surely all the connected components of $W$ are finite. Additionally, we provide a sufficient condition for the exponential decay of the size of a typical component.
We obtain Hanson-Wright inequalities for the quadratic form of a random vector with independent sparse random variables. Specifically, we consider cases where the components of the random vector are sparse $\alpha$-sub-exponential random variables with $\alpha>0$. Our proof relies on a novel combinatorial approach to estimate the moments of the random quadratic form.
Given suitable functions $V, \Psi:[0,\infty) \to [0,\infty)$, we obtain necessary and sufficient conditions on $V,\Psi$ for the existence of a metric measure space and a symmetric diffusion process that satisfies sub-Gaussian heat kernel estimates with volume growth profile $V$ and escape time profile $\Psi$. We prove sufficiency by constructing a new family of diffusions. Special cases of this construction also leads to a new family of infinite graphs whose simple random walks satisfy sub-Gaussian heat kernel estimates with prescribed volume growth and escape time profiles. In particular, these random walks on graphs generalizes earlier results of Barlow who considered the case $V(r)=r^\alpha$ and $\Psi(r)=r^\beta$ (Rev Mat Iberoam 2004). The family of diffusions we construct have martingale dimension one but can have arbitrarily high spectral dimension. Therefore our construction shows the impossibility of obtaining non-trivial \emph{lower} bounds on martingale dimension in terms of spectral dimension which is in contrast with upper bounds on martingale dimension using spectral dimension obtained by Hino (Probab Theory Relat Fields 2013).
We present a unifying approach to studying the replica symmetric solution in general diluted spin glass models on random $p$-uniform hypergraphs with sparsity parameter $\alpha$. Our result shows that there exist two key regimes in which the model exhibits replica symmetry and the free energy can be explicitly represented as the evaluation of an energy functional at the unique fixed point of a recursive distributional equation. One is called the high temperature regime, where the temperature and the sparsity parameter are essentially inversely proportional to each other; the other is the subcritical regime defined as $\alpha p (p-1)\leq 1$. In particular, the fact that the second regime is independent of the temperature parameter further allows us to deduce an analogous representation of the ground state energy in the subcritical regime. Along the way, we revisit several well-known formulas and also derive new ones for the free and ground state energies in the constraint satisfaction problem, Potts model, XY model, and continuous hardcore model.
We study standard processes with no negative jumps under the entrance boundary condition. Similarly to one-dimensional diffusions, we show that the process can be made into a Feller process by attaching the boundary point to the state space and extending the transition semigroup. We investigate the spectrum of the generator in detail via the scale function, characterizing it as the zeros of an entire function. As an application, we prove that under the strong Feller property, the convergence to the Yaglom limit of the process killed at the boundary is exponentially fast.
We present a computer assisted proof for a result concerning a three player betting game, introduced by Angel and Holmes. The three players start with initial capital $x, y, z > 0$ respectively. At each step of this game two players are selected at random to bet on the outcome of a fair coin toss, with the size of the bet being the largest possible, namely the total capital held by the poorer of the two players at that time. The main quantity of interest is the probability of player 1 being eliminated (reaching 0 capital) first. Angel and Holmes have shown that this probability is not monotone decreasing as a function of the initial capital $x$ of player 1. They conjecture that if $x < y < z$ then player 1 would be better off (less likely to be eliminated first) by swapping their capital with another player. In this paper we present a computer-assisted proof of this conjecture. To achieve this, we introduce the theoretical framework MeshItUp, and then perform a two-stage reduction to make MeshItUp computationally feasible, through the use of mixed-integer programming.
This paper contains a study on stochastic Volterra integral equations with fuzzy sets-values and involving on a constant retardation. Moreover, the form of the equation is symmetric in the sense that fuzzy stochastic integrals are placed on both sides of the equation. We show that the considered initial value problem formulated in terms of symmetric fuzzy stochastic Volterra integral equation is well-posed. In particular, we show that there exists a unique solution and this solution depends continuously on the parameters of the equation. The results are achieved with the conditions of Lipschitz continuity of drift and diffusion coefficients, and continuity of kernels
Consider a high-dimensional Wishart matrix $\bd{W}=\bd{X}^T\bd{X}$ where the entries of $\bd{X}$ are i.i.d. random variables with mean zero, variance one, and a finite fourth moment $\eta$. Motivated by problems in signal processing and high-dimensional statistics, we study the maximum of the largest eigenvalues of any two-by-two principal minors of $\bd{W}$. Under certain restrictions on the sample size and the population dimension of $\bd{W}$, we obtain the limiting distribution of the maximum, which follows the Gumbel distribution when $\eta$ is between 0 and 3, and a new distribution when $\eta$ exceeds 3. To derive this result, we first address a simpler problem on a new object named a deformed Gaussian orthogonal ensemble (GOE). The Wishart case is then resolved using results from the deformed GOE and a high-dimensional central limit theorem. Our proof strategy combines the Stein-Poisson approximation method, conditioning, U-statistics, and the Hájek projection. This method may also be applicable to other extreme-value problems. Some open questions are posed.
This paper investigates the structure of product systems of Hilbert spaces derived from Banach space-valued Lévy processes. We establish conditions under which these product systems are completely spatial and show that Gaussian Lévy processes with non-degenerate covariance always give rise to product systems of type I. Furthermore, we construct a continuum of non-isomorphic product systems of type \(\rm{II}\sb\infty\) from pure jump Lévy processes.
We prove equivalence between nonnegative distributional solutions of the fractional heat equation and caloric functions, i.e., functions satisfying the mean value property with respect to the space-time isotropic $\alpha$-stable process. We also provide sufficient conditions for the boundary and exterior data under which the solutions are classical and we give off-diagonal estimates for the derivatives of the Dirichlet heat kernel and the lateral Poisson kernel, which might be of their own interest.
We show that for a given DAG $G$, among all observational distributions of Bayesian networks over $G$ with arbitrary outcome spaces, the faithful distributions are `typical': they constitute a dense, open set with respect to the total variation metric. As a consequence, the set of faithful distributions is non-empty, and the unfaithful distributions are nowhere dense. We extend this result to the space of Bayesian networks, where the properties hold for Bayesian networks instead of distributions of Bayesian networks. As special cases, we show that these results also hold for the faithful parameters of the subclasses of linear Gaussian -- and discrete Bayesian networks, giving a topological analogue of the measure-zero results of Spirtes et al. (1993) and Meek (1995). Finally, we extend our topological results and the measure-zero results of Spirtes et al. and Meek to Bayesian networks with latent variables.
The Robbins-Siegmund theorem is one of the most important results in stochastic optimization, where it is widely used to prove the convergence of stochastic algorithms. We provide a quantitative version of the theorem, establishing a bound on how far one needs to look in order to locate a region of metastability in the sense of Tao. Our proof involves a metastable analogue of Doob's theorem for $L_1$-supermartingales along with a series of technical lemmas that make precise how quantitative information propagates through sums and products of stochastic processes. In this way, our paper establishes a general methodology for finding metastable bounds for stochastic processes that can be reduced to supermartingales, and therefore for obtaining quantitative convergence information across a broad class of stochastic algorithms whose convergence proof relies on some variation of the Robbins-Siegmund theorem. We conclude by discussing how our general quantitative result might be used in practice.
The Wright-Fisher diffusion is a fundamentally important model of evolution encompassing genetic drift, mutation, and natural selection. Suppose you want to infer the parameters associated with these processes from an observed sample path. Then to write down the likelihood one first needs to know the mutual arrangement of two path measures under different parametrizations; that is, whether they are absolutely continuous, equivalent, singular, and so on. In this paper we give a complete answer to this question by finding the separating times for the diffusion - the stopping time before which one measure is absolutely continuous with respect to the other and after which the pair is mutually singular. In one dimension this extends a classical result of Dawson on the local equivalence between neutral and non-neutral Wright-Fisher diffusion measures. Along the way we also develop new zero-one type laws for the diffusion on its approach to, and emergence from, the boundary. As an application we derive an explicit expression for the joint maximum likelihood estimator of the mutation and selection parameters and show that its convergence properties are closely related to the separating time.
In this work, we characterise the statistics of Markov chains by constructing an associated sequence of periodic differential operators. Studying the density of states of these operators reveals the absolutely continuous invariant measure of the Markov chain. This approach also leads to a direct proof of convergence to the invariant measure, along with explicit convergence rates. We show how our method can be applied to a class of related Markov chains including the logistic map, the tent map and Chebyshev maps of arbitrary order.
In this paper, we address three Principal--Agent problems in a moral hazard context and show that they are connected. We start by studying the problem of Principal with multiple Agents in cooperation. The term cooperation is manifested here by the fact that the agents optimize their criteria through Pareto equilibria. We show that as the number of agents tends to infinity, the principal's value function converges to the value function of a McKean--Vlasov control problem. Using the solution to this McKean--Vlasov control problem, we derive a constructive method for obtaining approximately optimal contracts for the principal's problem with multiple agents in cooperation. In a second step, we show that the problem of Principal with multiple Agents turns out to also converge, when the number of agents goes to infinity, towards a new Principal--Agent problem which is the Principal--Agent problem with Mckean--Vlasov dynamics. This is a Principal--Agent problem where the agent--controlled production follows a Mckean-Vlasov dynamics and the contract can depend of the distribution of the production. The value function of the principal in this setting is equivalent to that of the same McKean--Vlasov control problem from the multi--agent scenario. Furthermore, we show that an optimal contract can be constructed from the solution to this McKean--Vlasov control problem. We conclude by discussing, in a simple example, the connection of these problems with the multitask Principal--Agent problem which is a situation when a principal delegates multiple tasks that can be correlated to a single agent.
We study the mixing properties of the white-forced Navier--Stokes system in an infinite periodic pipe $\mathbb{R}\times\mathbb{T}$. Assuming that the noise is sufficiently non-degenerate, we prove the uniqueness of stationary measure and polynomial mixing in the dual-Lipschitz metric. The proof combines the coupling method with a Foiaş--Prodi type estimate and weighted growth estimates.
We study large deviations and mean-squared error (MSE) guarantees of a general framework of nonlinear stochastic gradient methods in the online setting, in the presence of heavy-tailed noise. Unlike existing works that rely on the closed form of a nonlinearity (typically clipping), our framework treats the nonlinearity in a black-box manner, allowing us to provide unified guarantees for a broad class of bounded nonlinearities, including many popular ones, like sign, quantization, normalization, as well as component-wise and joint clipping. We provide several strong results for a broad range of step-sizes in the presence of heavy-tailed noise with symmetric probability density function, positive in a neighbourhood of zero and potentially unbounded moments. In particular, for non-convex costs we provide a large deviation upper bound for the minimum norm-squared of gradients, showing an asymptotic tail decay on an exponential scale, at a rate $\sqrt{t} / \log(t)$. We establish the accompanying rate function, showing an explicit dependence on the choice of step-size, nonlinearity, noise and problem parameters. Next, for non-convex costs and the minimum norm-squared of gradients, we derive the optimal MSE rate $\widetilde{\mathcal{O}}(t^{-1/2})$. Moreover, for strongly convex costs and the last iterate, we provide an MSE rate that can be made arbitrarily close to the optimal rate $\mathcal{O}(t^{-1})$, improving on the state-of-the-art results in the presence of heavy-tailed noise. Finally, we establish almost sure convergence of the minimum norm-squared of gradients, providing an explicit rate, which can be made arbitrarily close to $o(t^{-1/4})$.
We prove the global well-posedness of the dynamical sine-Gordon model up to the third threshold, i.e., for parameters $\beta^2 < 6\pi$. The key novelty in our approach is the introduction of the so-called resonant equation, whose solution is entirely deterministic and completely captures the size of the solution to the dynamical sine-Gordon model. The probabilistic fluctuations in the dynamical sine-Gordon model are then controlled using uniform estimates for modified stochastic objects.
The Stefan problem with surface tension is well known to exhibit discontinuities in the associated moving aggregate (i.e., in the domain occupied by the solid), whose structure has only been understood under translational or radial symmetry so far. In this paper, we derive an auxiliary partial differential equation of second-order hyperbolic type, referred to as the cascade equation, that captures said discontinuities in the absence of any symmetry assumptions. Specializing to the one-phase setting, we introduce a novel (global) notion of weak solution to the cascade equation, which is defined as a limit of mean-field game equilibria. For the spatial dimension two, we show the existence of such a weak solution and prove a natural perimeter estimate on the associated moving aggregate.