We consider a problem introduced by Mossel and Ross [Shotgun assembly of labeled graphs, arXiv:1504.07682]. Suppose a random $n\times n$ jigsaw puzzle is constructed by independently and uniformly choosing the shape of each "jig" from $q$ possibilities. We are given the shuffled pieces. Then, depending on $q$, what is the probability that we can reassemble the puzzle uniquely? We say that two solutions of a puzzle are similar if they only differ by permutation of duplicate pieces, and rotation of rotationally symmetric pieces. In this paper, we show that, with high probability, such a puzzle has at least two non-similar solutions when $2\leq q \leq \frac{2}{\sqrt{e}}n$, all solutions are similar when $q\geq (2+\varepsilon)n$, and the solution is unique when $q=\omega(n)$.

We prove a result on non-clustering of particles in a two-dimensional Coulomb plasma, which holds provided that the inverse temperature $\beta$ satisfies $\beta>1$. As a consequence we obtain a result on crystallization as $\beta\to\infty$: the particles will, on a microscopic scale, appear at a certain distance from each other. The estimation of this distance is connected to Abrikosov's conjecture that the particles should freeze up according to a honeycomb lattice when $\beta\to\infty$.

In this paper we consider an interacting particle system in $\mathbb{R}^d$ modelled as a system of $N$ stochastic differential equations driven by Lévy processes. The limiting behaviour as the size $N$ grows to infinity is achieved as a law of large numbers for the empirical density process associated with the interacting particle system. We prove that the empirical process converges, uniformly in the space variable, to the solution of the $d$-dimensional fractal conservation law.

Jan 18 2017

math.PR arXiv:1701.04660v1

Consider a reaction-diffusion equation of the form \[\dotu(t\,,x)=\tfrac12 u"(t\,,x) + b(u(t\,,x)) + \sigma(u(t,x)) \xi(t,x),\]on $\R_+\times[0\,,1]$, with the Dirichlet boundary condition and a nice initial condition, where $\xi(t,x)$ is a space-time white noise, in the case that there exists $\varepsilon>0$ such that $| b(z)| \ge|z|(\log|z|)^{1+\varepsilon}$ for all sufficiently-large values of $|z|$. When $\sigma\equiv 0$, it is well known that such PDEs frequently have non-trivial stationary solutions. By contrast, Bonder and Groisman have recently shown that when $\sigma$ is constant and $\sigma \neq 0$, there is %the addition of any amount of space-time white noise to the reaction-diffusion equation often results in finite-time blowup. In this paper, we prove that the Bonder--Groisman condition is unimproveable by showing that the reaction-diffusion equation with noise is "typically" well posed when $| b(z)| =O(|z|\log_+|z|)$ as $|z|\to\infty$. We interpret the word "typically" in two essentially-different ways without altering the conclusions of our assertions

Jan 18 2017

math.PR arXiv:1701.04618v1

We introduce the class of continuous-time autoregressive moving-average (CARMA) processes in Hilbert spaces. As driving noises of these processes we consider Levy processes in Hilbert space. We provide the basic definitions, show relevant properties of these processes and establish the equivalents of CARMA processes on the real line. Finally, CARMA processes in Hilbert space are linked to the stochastic wave equation and functional autoregressive processes.

Jan 18 2017

math.PR arXiv:1701.04607v1

We consider noise perturbations of delay differential equations (DDE) experiencing Hopf bifurcation. The noise is assumed to be exponentially ergodic, i.e. transition density converges to stationary density exponentially fast uniformly in the initial condition. We show that, under an appropriate change of time scale, as the strength of the perturbations decreases to zero, the law of the critical eigenmodes converges to the law of a diffusion process (without delay). We prove the result only for scalar DDE. For vector-valued DDE without proofs see Phys.Rev.E v93, 062104.

Jan 18 2017

math.PR arXiv:1701.04598v1

Motivated by truncated EM method introduced by Mao (2015), a new explicit numerical method named modified truncated Euler-Maruyama method is developed in this paper. Strong convergence rates of the given numerical scheme to the exact solutions to stochastic differential equations are investigated under given conditions in this paper. Compared with truncated EM method, the given numerical simulation strongly converges to the exact solution at fixed time $T$ and over a time interval $[0,T]$ under weaker sufficient conditions. Meanwhile, the convergence rates are also obtained for both cases. Two examples are provided to support our conclusions.

In the work [Bull, Austr. Math. Soc. 85 (2012), 315-234], S.R. Moghadasi has shown how the decomposition of the $N$-fold product of Lebesgue measure on $\mathbb R^n$ implied by matrix polar decomposition can be used to derive the Blaschke-Petkantschin decomposition of measure formula from integral geometry. We use known formulas from random matrix theory to give a simplified derivation of the decomposition of Lebesgue product measure implied by matrix polar decomposition, applying too to the cases of complex and real quaternion entries, and we give corresponding generalisations of the Blaschke--Petkantschin formula. A number of applications to random matrix theory and integral geometry are given, including to the calculation of the moments of the volume content of the convex hull of $k \le N+1$ points in $\mathbb R^N$, $\mathbb C^N$ or $\mathbb H^N$ with a Gaussian or uniform distribution.

We consider randomly distributed mixtures of bonds of ferromagnetic and antiferromagnetic type in a two-dimensional square lattice with probability $1-p$ and $p$, respectively, according to an i.i.d. random variable. We study minimizers of the corresponding nearest-neighbour spin energy on large domains in ${\mathbb Z}^2$. We prove that there exists $p_0$ such that for $p\le p_0$ such minimizers are characterized by a majority phase; i.e., they take identically the value $1$ or $-1$ except for small disconnected sets. A deterministic analogue is also proved.

A recent paper \citeKMMO introduced the stochastic U_q(A_n^(1)) vertex model. The stochastic S-matrix is related to the R-matrix of the quantum group U_q(A_n^(1)) by a gauge transformation. We will show that a certain function D^+_\mu intertwines with the transfer matrix and its space reversal. When interpreting the transfer matrix as the transition matrix of a discrete-time totally asymmetric particle system on the one-dimensional lattice Z, the function D^+_\mu becomes a Markov duality function D_\mu which only depends on q and the vertical spin parameters \mu_x. By considering degenerations in the spectral parameter, the duality results also hold on a finite lattice with closed boundary conditions, and for a continuous-time degeneration. This duality function had previously appeared in a multi-species ASEP(q,j) process. The proof here uses that the R-matrix intertwines with the co-product, but does not explicitly use the Yang-Baxter equation. It will also be shown that the stochastic U_q(A_n^(1)) is a multi-species version of a stochastic vertex model studied in \citeBP,CP. This will be done by generalizing the fusion process of \citeCP and showing that it matches the fusion of \citeKRL up to the gauge transformation. We also show, by direct computation, that the multi-species q-Hahn Boson process (which arises at a special value of the spectral parameter) also satisfies duality with respect to D_0, generalizing the single-species result of \citeC.

We consider a sparse linear regression model Y=X\beta^*+W where X has a Gaussian entries, W is the noise vector with mean zero Gaussian entries, and \beta^* is a binary vector with support size (sparsity) k. Using a novel conditional second moment method we obtain a tight up to a multiplicative constant approximation of the optimal squared error \min_\beta\|Y-X\beta\|_2, where the minimization is over all k-sparse binary vectors \beta. The approximation reveals interesting structural properties of the underlying regression problem. In particular, a) We establish that n^*=2k\log p/\log (2k/\sigma^2+1) is a phase transition point with the following "all-or-nothing" property. When n exceeds n^*, (2k)^-1\|\beta_2-\beta^*\|_0≈0, and when n is below n^*, (2k)^-1\|\beta_2-\beta^*\|_0≈1, where \beta_2 is the optimal solution achieving the smallest squared error. With this we prove that n^* is the asymptotic threshold for recovering \beta^* information theoretically. b) We compute the squared error for an intermediate problem \min_\beta\|Y-X\beta\|_2 where minimization is restricted to vectors \beta with \|\beta-\beta^*\|_0=2k \zeta, for \zeta∈[0,1]. We show that a lower bound part \Gamma(\zeta) of the estimate, which corresponds to the estimate based on the first moment method, undergoes a phase transition at three different thresholds, namely n_\textinf,1=\sigma^2\log p, which is information theoretic bound for recovering \beta^* when k=1 and \sigma is large, then at n^* and finally at n_\textLASSO/CS. c) We establish a certain Overlap Gap Property (OGP) on the space of all binary vectors \beta when n\le ck\log p for sufficiently small constant c. We conjecture that OGP is the source of algorithmic hardness of solving the minimization problem \min_\beta\|Y-X\beta\|_2 in the regime n<n_\textLASSO/CS.

Based on numerical simulation and local stability analysis we describe the structure of the phase space of the edge/triangle model of random graphs. We support simulation evidence with mathematical proof of continuity and discontinuity for many of the phase transitions. All but one of themany phase transitions in this model break some form of symmetry, and we use this model to explore how changes in symmetry are related to discontinuities at these transitions.