Quadrature rule for Abel’s equations: Uniformly approximating fractional derivatives

An automatic quadrature method is presented for approximating fractional derivative D^qf(x)$D^{q}f\left(x\right)$ of a given function f(x)$f\left(x\right)$, which is defined by an indefinite integral involving f(x)$f\left(x\right)$. The present method interpolates f(x)$f\left(x\right)$ in terms of the Chebyshev polynomials in the range [0, 1] to approximate the fractional derivative D^qf(x)$D^{q}f\left(x\right)$ uniformly for 0≤x≤1$0\le x\le 1$, namely the error is bounded independently of x$x$. Some numerical examples demonstrate the performance of the present automatic method.     

Regularity for the fractional Gelfand problem up to dimension 7

We study the problem (−Δ)^su=λe^u${(-\mathrm{\Delta })}^{s}u=\lambda e^u$ in a bounded domain Ω⊂Rⁿ$\Omega \subset {\mathbb{R}}^n$, where λ   is a positive parameter. More precisely, we study the regularity of the extremal solution to this problem. Our main result yields the boundedness of the extremal solution in dimensions n≤7$n\le 7$ for all s∈(0,1)$s\in (0,1)$ whenever Ω   is, for every i=1,...,n$i=1,...,n$, convex in the _xi$x_i$-direction and symmetric with respect to {_xi=0}$\{x_i=0\}$. The same holds if n=8$n=8$ and s?0.28206...$s?0.28206...$, or if n=9$n=9$ and s?0.63237...$s?0.63237...$. These results are new even in the unit ball Ω=_B1$\Omega =B_1$.     

Transformation formulas for fractional Brownian motion

We derive a Molchan–Golosov-type integral transform which changes fractional Brownian motion of arbitrary Hurst index K$K$ into fractional Brownian motion of index H$H$. Integration is carried out over [0,t]$[0,t]$, t>0$t>0$. The formula is derived in the time domain. Based on this transform, we construct a prelimit which converges in L²(P)$L^2\left(\mathbb{P}\right)$-sense to an analogous, already known Mandelbrot–Van Ness-type integral transform, where integration is over (−∞,t]$(-\infty ,t]$, t>0$t>0$.     

On testing Hamiltonicity of graphs

Let us fix a function f(n)=o(nlnn)$f\left(n\right)=o(nlnn)$ and real numbers 0≤α<β≤1$0\le \alpha <\beta \le 1$. We present a polynomial time algorithm which, given a directed graph G$G$ with n$n$ vertices, decides either that one can add at most βn$\beta n$ new edges to G$G$ so that G$G$ acquires a Hamiltonian circuit or that one cannot add αn$\alpha n$ or fewer new edges to G$G$ so that G$G$ acquires at least e^−f(n)n!$e^{-f\left(n\right)}n!$ Hamiltonian circuits, or both.     

Weak convergence of the Stratonovich integral with respect to a class of Gaussian processes

For a Gaussian process X$X$ and smooth function f$f$, we consider a Stratonovich integral of f(X)$f\left(X\right)$, defined as the weak limit, if it exists, of a sequence of Riemann sums. We give covariance conditions on X$X$ such that the sequence converges in law. This gives a change-of-variable formula in law with a correction term which is an Itô integral of f^?$f^{?}$ with respect to a Gaussian martingale independent of X$X$. The proof uses Malliavin calculus and a central limit theorem from Nourdin and Nualart (2010) [8]. This formula was known for fBm with H=1/6$H=1/6$ Nourdin et al. (2010) [9]. We extend this to a larger class of Gaussian processes.     

Monotonicity of the reflected Bessel transition density on the diagonal

For α∈R$\alpha \in \mathbb{R}$, let p^R(t,x,x)$p^R(t,x,x)$ denote the diagonal of the transition density of the α$\alpha$-Bessel process in (0,1]$(0,1]$, killed at 0 and reflected at 1. As a function of x$x$, if either α≥3$\alpha \ge 3$ or α=1$\alpha =1$, then for t>0$t>0$, the diagonal is nondecreasing. This monotonicity property fails if 1≠α<3$1\ne \alpha <3$.     

Selections and their absolutely continuous invariant measures

Let I=[0,1]$I=[0,1]$ and let P be a partition of I   into a finite number of intervals. Let _τ1${\tau }_1$, _τ2${\tau }_2$; I→I$I\to I$ be two piecewise expanding maps on P  . Let G⊂I×I$G\subset I\times I$ be the region between the boundaries of the graphs of _τ1${\tau }_1$ and _τ2${\tau }_2$. Any map τ:I→I$\tau :I\to I$ that takes values in G is called a selection of the multivalued map defined by G  . There are many results devoted to the study of the existence of selections with specified topological properties. However, there are no results concerning the existence of selection with measure-theoretic properties. In this paper we prove the existence of selections which have absolutely continuous invariant measures (acim). By our assumptions we know that _τ1${\tau }_1$ and _τ2${\tau }_2$ possess acims preserving the distribution functions F⁽¹⁾$F^{(1)}$ and F⁽²⁾$F^{(2)}$. The main result shows that for any convex combination F   of F⁽¹⁾$F^{(1)}$ and F⁽²⁾$F^{(2)}$ we can find a map η   with values between the graphs of _τ1${\tau }_1$ and _τ2${\tau }_2$ (that is, a selection) such that F is the η-invariant distribution function. Examples are presented. We also study the relationship of the dynamics of our multivalued maps to random maps.     

Constructing cubature formulae of degree 5 with few points

This paper is devoted to construct a family of fifth degree cubature formulae for n$n$-cube with symmetric measure and n$n$-dimensional spherically symmetrical region. The formula forn$n$-cube contains at most n²+5n+3$n^2+5n+3$ points and for n$n$-dimensional spherically symmetrical region contains only n²+3n+3$n^2+3n+3$ points. Moreover, the numbers can be reduced to n²+3n+1$n^2+3n+1$ and n²+n+1$n^2+n+1$ if n=7$n=7$ respectively, the latter of which is minimal.     

On functional limits of short- and long-memory linear processes with GARCH(1,1) noises

This paper considers the short- and long-memory linear processes with GARCH (1,1) noises. The functional limit distributions of the partial sum and the sample autocovariances are derived when the tail index α$\alpha$ is in (0,2)$(0,2)$, equal to 2, and in (2,∞)$(2,\infty )$, respectively. The partial sum weakly converges to a functional of α$\alpha$-stable process when α<2$\alpha <2$ and converges to a functional of Brownian motion when α≥2$\alpha \ge 2$. When the process is of short-memory and α<4$\alpha <4$, the autocovariances converge to functionals of α/2$\alpha /2$-stable processes; and if α≥4$\alpha \ge 4$, they converge to functionals of Brownian motions. In contrast, when the process is of long-memory, depending on α$\alpha$ and β$\beta$ (the parameter that characterizes the long-memory), the autocovariances converge to either (i) functionals of α/2$\alpha /2$-stable processes; (ii) Rosenblatt processes (indexed by β$\beta$, 1/2<β<3/4$1/2<\beta <3/4$); or (iii) functionals of Brownian motions. The rates of convergence in these limits depend on both the tail index α$\alpha$ and whether or not the linear process is short- or long-memory. Our weak convergence is established on the space of càdlàg functions on [0,1]$[0,1]$ with either (i) the _J1$J_1$ or the _M1$M_1$ topology (Skorokhod, 1956); or (ii) the weaker form S$S$ topology (Jakubowski, 1997). Some statistical applications are also discussed.     

On the length of an external branch in the Beta-coalescent

In this paper, we consider Beta(2−α,α)$(2-\alpha ,\alpha )$ (with 1<α<2$1<\alpha <2$) and related Λ$\Lambda$-coalescents. If T⁽ⁿ⁾$T^{\left(n\right)}$ denotes the length of a randomly chosen external branch of the n$n$-coalescent, we prove the convergence of n^α−1T⁽ⁿ⁾$n^{\alpha -1}{T}^{\left(n\right)}$ when n$n$ tends to ∞$\infty$, and give the limit. To this aim, we give asymptotics for the number σ⁽ⁿ⁾${\sigma }^{\left(n\right)}$ of collisions which occur in the n$n$-coalescent until the end of the chosen external branch, and for the block counting process associated with the n$n$-coalescent.     

Boundedness and blowup for nonlinear degenerate parabolic equations

The author deals with the quasilinear parabolic equation _ut=[u^α+g(u)]Δu+bu^α+1+f(u,∇u)$u_t=[u^{\alpha }+g\left(u\right)]\Delta u+b{u}^{\alpha +1}+f(u,\nabla u)$ with Dirichlet boundary conditions in a bounded domain Ω$\Omega$, where f$f$ and g$g$ are lower-order terms. He shows that, under suitable conditions on f$f$ and g$g$, whether the solution is bounded or blows up in a finite time depends only on the first eigenvalue of −Δ$-\Delta$ in Ω$\Omega$ with Dirichlet boundary condition. For some special cases, the result is sharp.     

A note on distance-regular graphs with a small number of vertices compared to the valency

In this note we study distance-regular graphs with a small number of vertices compared to the valency. We show that for a given α>2$\alpha >2$, there are finitely many distance-regular graphs Γ$\Gamma$ with valency k$k$, diameter D≥3$D\ge 3$ and v$v$ vertices satisfying v≤αk$v\le \alpha k$ unless (D=3$D=3$ and Γ$\Gamma$ is imprimitive) or (D=4$D=4$ and Γ$\Gamma$ is antipodal and bipartite). We also show, as a consequence of this result, that there are finitely many distance-regular graphs with valency k≥3$k\ge 3$, diameter D≥3$D\ge 3$ and _c2≥εk$c_2\ge \epsilon k$ for a given 0<ε<1$0<\epsilon <1$ unless (D=3$D=3$ and Γ$\Gamma$ is imprimitive) or (D=4$D=4$ and Γ$\Gamma$ is antipodal and bipartite).     

A boundary blow-up elliptic problem with an inhomogeneous term

By a perturbation method and constructing comparison functions, we reveal how the inhomogeneous term h$h$ affects the exact asymptotic behaviour of solutions near the boundary to the problem △u=b(x)g(u)+λh(x)$△u=b\left(x\right)g\left(u\right)+\lambda h\left(x\right)$, u>0$u>0$ in Ω$\Omega$, u_|∂Ω=∞$u{|}_{\partial \Omega }=\infty$, where Ω$\Omega$ is a bounded domain with smooth boundary in R^N${\mathbb{R}}^N$, λ>0$\lambda >0$, g∈C¹[0,∞)$g\in C^1[0,\infty )$ is increasing on [0,∞)$[0,\infty )$, g(0)=0$g\left(0\right)=0$, g^′$g^{\prime }$ is regularly varying at infinity with positive index ρ$\rho$, the weight b$b$, which is non-trivial and non-negative in Ω$\Omega$, may be vanishing on the boundary, and the inhomogeneous term h$h$ is non-negative in Ω$\Omega$ and may be singular on the boundary.     

Estimating the ground state energy of the Schrödinger equation for convex potentials

In 2011, the fundamental gap conjecture for Schrödinger operators was proven. This can be used to estimate the ground state energy of the time-independent Schrödinger equation with a convex potential and relative error ε$\epsilon$. Classical deterministic algorithms solving this problem have cost exponential in the number of its degrees of freedom d$d$. We show a quantum algorithm, that is based on a perturbation method, for estimating the ground state energy with relative error ε$\epsilon$. The cost of the algorithm is polynomial in d$d$ and ε⁻¹${\epsilon }^{-1}$, while the number of qubits is polynomial in d$d$ and logε⁻¹$log{\epsilon }^{-1}$. In addition, we present an algorithm for preparing a quantum state that overlaps within 1−δ,δ∈(0,1)$1-\delta ,\delta \in (0,1)$, with the ground state eigenvector of the discretized Hamiltonian. This algorithm also approximates the ground state with relative error ε$\epsilon$. The cost of the algorithm is polynomial in d$d$, ε⁻¹${\epsilon }^{-1}$ and δ⁻¹${\delta }^{-1}$, while the number of qubits is polynomial in d$d$, logε⁻¹$log{\epsilon }^{-1}$ and logδ⁻¹$log{\delta }^{-1}$.     

Local asymptotic powers of nonparametric and semiparametric tests for fractional integration

The paper concerns testing long memory for fractionally integrated nonlinear processes. We show that the exact local asymptotic power is of order O[(logn)⁻¹]$O\left[{(logn)}^{-1}\right]$ for four popular nonparametric tests and is O(m^−1/2)$O\left(m^{-1/2}\right)$, where m$m$ is the bandwidth which is allowed to grow as fast as n^κ$n^{\kappa }$, κ∈(0,2/3)$\kappa \in (0,2/3)$, for the semiparametric Lagrange multiplier (LM) test proposed by Lobato and Robinson [I. Lobato, P.M. Robinson, A nonparametric test for I(0)$I\left(0\right)$, Rev. Econom. Stud. 68 (1998) 475–495]. Our theory provides a theoretical justification for the empirical findings in finite sample simulations by Lobato and Robinson [I. Lobato, P.M. Robinson, A nonparametric test for I(0)$I\left(0\right)$, Rev. Econom. Stud. 68 (1998) 475–495] and Giraitis et al. [L. Giraitis, P. Kokoszka, R. Leipus, G. Teyssiére, Rescaled variance and related tests for long memory in volatility and levels, J. Econometrics 112 (2003) 265–294] that nonparametric tests have lower power than LM tests in detecting long memory.     

Parametric inference for discretely observed multidimensional diffusions with small diffusion coefficient

We consider a multidimensional diffusion X$X$ with drift coefficient b(α,_Xt)$b(\alpha ,X_t)$ and diffusion coefficient ?σ(β,_Xt)$?\sigma (\beta ,X_t)$. The diffusion sample path is discretely observed at times _tk=kΔ$t_k=k\Delta$ for k=1…n$k=1\dots n$ on a fixed interval [0,T]$[0,T]$. We study minimum contrast estimators derived from the Gaussian process approximating X$X$ for small ?$?$. We obtain consistent and asymptotically normal estimators of α$\alpha$ for fixed Δ$\Delta$ and ?→0$?\to 0$ and of (α,β)$(\alpha ,\beta )$ for Δ→0$\Delta \to 0$ and ?→0$?\to 0$ without any condition linking ?$?$ and Δ$\Delta$. We compare the estimators obtained with various methods and for various magnitudes of Δ$\Delta$ and ?$?$ based on simulation studies. Finally, we investigate the interest of using such methods in an epidemiological framework.     

Blow-up and propagation of disturbances for fast diffusion equations

This paper is concerned with the Cauchy problem for the fast diffusion equation _ut−Δu^m=αu^_p1$u_t-\Delta u^m=\alpha u^{p_1}$ in R^N${\mathbb{R}}^N$ (N≥1$N\ge 1$), where m∈(0,1)$m\in (0,1)$, _p1>1$p_1>1$ and α>0$\alpha >0$. The initial condition _u0$u_0$ is assumed to be continuous, nonnegative and bounded. Using a technique of subsolutions, we set up sufficient conditions on the initial value _u0$u_0$ so that u(t,x)$u(t,x)$ blows up in finite time, and we show how to get estimates on the profile of u(t,x)$u(t,x)$ for small enough values of t>0$t>0$.