2025-06-19 | | Total: 21
We study the discrete Bak-Sneppen model introduced by Barbay and Kenyon (2001) "On the discrete Bak-Sneppen model of self-organized criticality". We extend their results as well as the non-triviality result of Meester and Znamenskiy (2002) for a finite segment of $\mathbb{Z}^1$ with the periodic boundary condition to a large class of graphs, by using coupling between the Bak-Sneppen model and the oriented percolation in a quadrant. This allows us to avoid dealing with the so-called avalanches, thus simplifying many arguments.
This work concerns about stochastic Burgers type equations with reflection. First of all, by means of the equicontinuous uniform Laplace principle, we prove the Freidlin-Wentzell uniform large deviation principle for these equations uniformly on bounded sets. Then based on this result, we establish the Dembo-Zeitouni uniform large deviation principle for these equations uniformly on compact sets. Finally, an averaging principle result for these equations is obtained through the time discretization approach.
In this paper, we study the connection between a general class of Conditional Mckean-Vlasov Stochastic Differential Equations (CMVSDEs) and its corresponding (infinite dimensional) Conditional Fokker-Planck Equation. The CMVSDE under consideration is similar to the one studied in [4], which is a non-trivial generalization of the McKean-Vlasov SDE with common noise and is closely related to a new type of non-linear Zakai equation that has not been studied in the literature. The main purpose of this paper is to establish the superposition principles among the three subjects so that their well-posedness can imply each other. More precisely, we shall first prove the superposition principle between the non-linear Zakai equation, a non-linear measure-valued stochastic PDE, and a CMVSDE; and then prove the superposition principle between an infinite dimensional conditional Fokker-Planck equation and the nonlinear Zakai equations. It is worth noting that none of the (weak) well-posedness of these SDEs/SPDEs in such generality have been investigated in the literature.
In this paper, we study the Exponential Random Graph Models (ERGMs) conditioning on the number of edges. In subcritical region of model parameters, we prove a conditional Central Limit Theorem (CLT) with explicit mean and variance for the number of two stars. This generalizes the corresponding result in the literature for the Erdős--Rényi random graph. To prove our main result, we develop a new conditional CLT via exchangeable pairs based on the ideas of Dey and Terlov. Our key technical contributions in the application to ERGMs include establishing a linearity condition for an exchangeable pair involving two star counts, a local CLT for edge counts, as well as new higher-order concentration inequalities. Our approach also works for general subgraph counts, and we give a conjectured form of their conditional CLT.
We investigate the large deviation probabilities of first passage times (FPT) of discrete-time supercritical non-lattice branching random walks (BRWs) in $\mathbb{R}^d$ where $d\geq 1$. Here, the FPT refers to the first time the BRW enters a ball of radius one with a distance $x$ from the origin, conditioned upon survival. We then apply the spine decomposition technique to construct an asymptotically optimal polynomial-time algorithm for computing the lower large deviation probabilities of the FPT. Our algorithm is also verified numerically. Our analysis provides a further understanding of the microstructural feature that can help characterize the strength of polymers.
In the existing works, stochastic sets $\mathbb{B}$ of interval type, along with $\mathbb{B}$-stochastic processes, were introduced within the framework of stochastic analysis. In this paper, we undertake the construction of $\mathbb{B}$-stochastic integration by exploring three novel types of $\mathbb{B}$-stochastic integrals: Stieltjes integrals of $\mathbb{B}$-predictable processes with respect to $\mathbb{B}$-adapted processes with finite variation, stochastic integrals of $\mathbb{B}$-predictable processes with respect to $\mathbb{B}$-inner local martingales, and stochastic integrals of $\mathbb{B}$-predictable processes with respect to $\mathbb{B}$-inner semimartingales. These $\mathbb{B}$-stochastic integrals are exclusively defined on subsets $\mathbb{B}$, with values outside the scope of $\mathbb{B}$ being deemed irrelevant. Additionally, we present several notable consequences, including the relationship between $\mathbb{B}$-stochastic integrals and existing stochastic integrals, as well as Itô's formula for $\mathbb{B}$-inner semimartingales. In the context of models pertaining to uncertain time-horizons in mathematical finance, we establish essentials of mathematical finance for general markets characterized by sudden-stop horizons. This is achieved by defining self-financing strategies, admissible strategies, and no-arbitrary conditions. In such financial markets, the exclusivity characteristic inherent in $\mathbb{B}$-stochastic integrals offers investors a viable alternative approach. This approach enables them to effectively filter out unnecessary information pertaining to asset price dynamics and portfolio strategies that extend beyond the predefined time-horizons.
We introduce a novel percolation model that generalizes the classical Random Connection Model (RCM) to a random simplicial complex, allowing for a more refined understanding of connectivity and emergence of large-scale structures in random topological spaces. Regarding percolation with respect to the notion of up-connectivity, we establish the existence of a sharp phase transition for the appearance of a giant component, akin to the well-known threshold behavior in random graphs. This sharp phase transition is, in its generality, new even for the classical RCM as a random graph. As special cases, we obtain sharp phase transitions for the Vietoris-Rips complex, the Cech complex, and the Boolean model, allowing us to identify which properties of these well-known percolation models are actually required.
We investigate the existence and uniqueness of non-Markovian second-order backward stochastic differential equations with an uncertain terminal horizon and establish comparison principles under the assumption that the driver is Lipschitz continuous. The terminal time is both random and exogenous, and it may not be adapted to the Brownian filtration, leading to a singular jump in the 2BSDE decomposition. We also provide a connection between this new class of 2BSDE and a fully nonlinear PDE in a Markovian setting. Our theoretical results are applied to non-Markovian stochastic control problems in two settings: (1) when an agent seeks to maximize utility from a payoff received at an uncertain terminal time by controlling both the drift and volatility of a diffusion process; and (2) when the agent contends with volatility uncertainty stemming from an external source, referred to as Nature, and optimizes the drift in a worst-case scenario for the ambiguous volatility. We term this class of problems erratic stochastic control, reflecting the dual uncertainty in both model parameters and the timing of the terminal horizon.
We study coarse-graining methods for stochastic differential equations. In particular we consider averaging and a type of projection operator method, sometimes referred to as effective dynamic via conditional expectations. The projection method (PM) we consider is related to the ``mimicking marginals'' coarse graining approach proposed by Gyöngy. The first contribution of this paper is to provide further theoretical background for the PM and a rigorous link to the Gyöngy method. Moreover, we compare PM and averaging. While averaging applies to systems with time scale separation, the PM can in principle be applied irrespective of this. However it is often assumed that the two methods coincide in presence of scale separation. The second contribution of this paper is to make this statement precise, provide sufficient conditions under which these two methods coincide and then show -- via examples and counterexamples -- that this needs not be the case.
Consider a Pólya urn where a drawn ball of colour $i$ is replaced together with a fixed number $m_i$ of balls of the same colour. We give a simple proof that if, for example, there are two colours and the urn starts with more balls of colour 1 than 2, and $m_1\ge m_2$, then there is a positive probability that there always will be more balls of colour 1 than colour 2.
This paper delves into the well-posedness and the numerical approximation of non-autonomous stochastic differential algebraic equations (SDAEs) with nonlinear local Lipschitz coefficients that satisfy the more general monotonicity condition called Khasminskii condition. The key challenge is the presence of a singular matrix which makes the numerical integration hard and heavy. To address this challenge, we propose a novel numerical scheme based on semi-implicit method for the drift component of the SDAEs. More precisely we split the drift term as the sum of a linear term and a nonlinear term. The linear part is approximated implicitly, while the nonlinear part is approximated explicitly. The linear component's role is to handle the singularity issues during the numerical integration without the resolution of nonlinear algebraic equations in the constraint equations. This novel scheme is therefore very efficient for SDAEs in high dimension that come after the spatial discretisation of stochastic partial differential algebraic equations (SPDAEs). To prove the pathwise convergence of our novel scheme, we first derive a equivalent scheme called dual scheme, suitable for mathematical analysis and linked to the inherent stochastic differential equation resulting from the elimination of constraints in the initial SDAEs. We prove that our novel scheme converges to the exact solution with rate $\frac{1}{2}-\epsilon$, for arbitrary $\epsilon>0$ in the pathwise sense. Numerical simulations are performed to demonstrate the efficiency of the scheme in high dimension and to show that our theoretical results are in agreement with numerical experiments.
In the paper average reward per unit time and average risk sensitive reward functionals are considered for controlled nonhomogeneous Markov processes. Existence of solutions to suitable Bellman equations is shown. Continuity of the value functions with respect to risk parameter is also proved. Finally stability of functionals with respect to pointwise convergence of Markov controls is studied.
We establish conditioned stochastic stability of equilibrium states for Hölder potentials on uniformly hyperbolic sets.
In this paper we introduce and study the Maximum-Average Subtensor ($p$-MAS) problem, in which one wants to find a subtensor of size $k$ of a given random tensor of size $N$, both of order $p$, with maximum sum of entries. We are motivated by recent work on the matrix case of the problem in which several equilibrium and non-equilibrium properties have been characterized analytically in the asymptotic regime $1 \ll k \ll N$, and a puzzling phenomenon was observed involving the coexistence of a clustered equilibrium phase and an efficient algorithm which produces submatrices in this phase. Here we extend previous results on equilibrium and algorithmic properties for the matrix case to the tensor case. We show that the tensor case has a similar equilibrium phase diagram as the matrix case, and an overall similar phenomenology for the considered algorithms. Additionally, we consider out-of-equilibrium landscape properties using Overlap Gap Properties and Franz-Parisi analysis, and discuss the implications or lack-thereof for average-case algorithmic hardness.
We consider real tensors of order $D$, that is $D$-dimensional arrays of real numbers $T_{a^1a^2 \dots a^D}$, where each index $a^c$ can take $N$ values. The tensor entries $T_{a^1a^2 \dots a^D}$ have no symmetry properties under permutations of the indices. The invariant polynomials built out of the tensor entries are called trace invariants. We prove that for a Gaussian random tensor with $D\ge 3$ indices (that is such that the entries $T_{a^1a^2 \dots a^D}$ are independent identically distributed Gaussian random variables) the cumulant, or connected expectation, of a product of trace invariants is not always suppressed in scaling in $N$ with respect to the product of the expectations of the individual invariants. Said otherwise, not all the multi-trace expectations factor at large $N$ in terms of the single-trace ones and the Gaussian scaling is not subadditive on the connected components. This is in stark contrast to the $D=2$ case of random matrices in which the multi-trace expectations always factor at large $N$. The best one can do for $D\ge 3$ is to identify restricted families of invariants for which the large $N$ factorization holds and we check that this indeed happens when restricting to the family of melonic observables, the dominant family in the large $N$ limit.
We study the energy landscape of the Random Energy model (REM) integrated along trajectories of the simple random walk on the hypercube. We show that the quenched cumulant generating function of the time integral of the REM energy undergoes phase transitions in the large $N$ limit for trajectories of any time extent, and identify phases distinguished by the activity and value of the time integral. This is achieved by relating the dynamical behavior to the spectral properties of Hamiltonians associated with the Quantum Random Energy Model (QREM). Of independent interest are deterministic $ \ell^p $-properties of the resolvents of such Hamiltonians, which we establish.
Representation results for absolutely continuous curves $\mu:[0,T]\to \mathcal{P}_p(\mathbb{R}^d)$, $p>1$, with values in the Wasserstein space $(\mathcal{P}_p(\mathbb{R}^d),W_p)$ of Borel probability measures in $\mathbb{R}^d$ with finite $p$-moment, provide a crucial tool to study evolutionary PDEs in a measure-theoretic setting. They are strictly related to the superposition principle for measure-valued solutions to the continuity equation. This paper addresses the extension of these results to the case $p=1$, and to curves $\mu:[0,+\infty)\to\mathcal{P}_1(\mathbb{R}^d)$ that are only of bounded variation in time: in the corresponding continuity equation, the flux measure $\nu\in\mathcal{M}_{loc}([0,+\infty)\times\mathbb{R}^{d};\mathbb{R}^{d})$ thus possesses a non-trivial singular part w.r.t. $\mu$ in addition to the absolutely continuous part featuring the velocity field. Firstly, we carefully address the relation between curves in ${\rm BV}_{loc}([0,+\infty);\mathcal{P}_1(\mathbb{R}^d))$ and solutions to the associated continuity equation, among which we select those with minimal singular (contribution to the) flux $\nu$. We show that, with those distinguished solutions it is possible to associate an `auxiliary' continuity equation, in an augmented phase space, solely driven by its velocity field. For that continuity equation, a standard version of the superposition principle can be thus obtained. In this way, we derive a first probabilistic representation of the pair $(\mu,\nu)$ solutions by projection over the time and space marginals. This representation involves Lipschitz trajectories in the augmented phase space, reparametrized in time and solving the characteristic system of ODEs. Finally, for the same pair $(\mu,\nu)$ we also prove a superposition principle in terms of BV curves on the actual time interval, providing a fine description of their behaviour at jump points.
Algorithms for jointly obtaining projection estimates of the density and distribution function of a random variable using the Legendre polynomials are proposed. For these algorithms, a problem of the conditional optimization is solved. Such an optimization allows one increasing the approximation accuracy with a minimum computational costs. The proposed algorithms are tested on examples with different degree of smoothness of the density.
In classical probability theory, the convergence of empirical frequencies to theoretical probabilities: as captured by the Law of Large Numbers (LLN): is treated as axiomatic and emergent from statistical assumptions such as independence and identical distribution. In this work, a novel dynamical framework is constructed in which convergence arises as a consequence of structured evolution in outcome space, rather than a statistical postulate. Through this formalism, statistical convergence is derived dynamically, revealing an internal structure to randomness and exposing entanglement between successive trials. The system recovers classical LLN behaviour in the large-number limit, while predicting deviations, transient fluctuations, and geometric asymmetries in the early regime. This work inaugurates a new paradigm: dynamical probability mechanics: in which randomness is modelled not as a sequence of disconnected stochastic events, but as a physically structured, feedback-driven process. The theory provides a novel explanatory layer beneath statistical laws and opens pathways toward a mechanistic foundation of probability itself.
In this paper, we propose an accelerated version for the Sinkhorn algorithm, which is the reference method for computing the solution to Entropic Optimal Transport. Its main draw-back is the exponential slow-down of convergence as the regularization weakens $\varepsilon \rightarrow 0$. Thanks to spectral insights on the behavior of the Hessian, we propose to mitigate the problem via an original spectral warm-start strategy. This leads to faster convergence compared to the reference method, as also demonstrated in our numerical experiments.
In this paper we consider the problem of computing the stationary distribution of nearly completely decomposable Markov processes, a well-established area in the classical theory of Markov processes with broad applications in the design, modeling, analysis and optimization of computer systems. We design general classes of algorithmic solution approaches that exploit forms of mixed-precision computation to significantly reduce computation times and that exploit forms of iterative approximate methods to mitigate the impact of inaccurate computations, further reduce computation times, and ensure convergence. Then we derive a mathematical analysis of our general algorithmic approaches that establishes theoretical results on approximation errors, convergence behaviors, and other algorithmic properties. Numerical experiments demonstrate that our general algorithmic approaches provide significant improvements in computation times over the most efficient existing numerical methods.