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}$.     

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 probabilistic results for the discrepancy of a hybrid-Monte Carlo sequence

In many applications it has been observed that hybrid-Monte Carlo sequences perform better than Monte Carlo and quasi-Monte Carlo sequences, especially in difficult problems. For a mixed s$s$-dimensional sequence m$m$, whose elements are vectors obtained by concatenating d$d$-dimensional vectors from a low-discrepancy sequence q$q$ with (s−d)$(s-d)$-dimensional random vectors, probabilistic upper bounds for its star discrepancy have been provided. In a paper of G. Ökten, B. Tuffin and V. Burago [G. Ökten, B. Tuffin, V. Burago, J. Complexity 22 (2006), 435–458] it was shown that for arbitrary ε>0$\epsilon >0$ the difference of the star discrepancies of the first N$N$ points of m$m$ and q$q$ is bounded by ε$\epsilon$ with probability at least 1−2exp(−ε²N/2)$1-2exp(-{\epsilon }^{2}N/2)$ for N$N$ sufficiently large. The authors did not study how large N$N$ actually has to be and if and how this actually depends on the parameters s$s$ and ε$\epsilon$. In this note we derive a lower bound for N$N$, which significantly depends on s$s$ and ε$\epsilon$. Furthermore, we provide a probabilistic bound for the difference of the star discrepancies of the first N$N$ points of m$m$ and q$q$, which holds without any restrictions on N$N$. In this sense it improves on the bound of Ökten, Tuffin and Burago and is more helpful in practice, especially for small sample sizes N$N$. We compare this bound to other known bounds.     

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$.     

Zero separation results for solutions of second order linear differential equations

The oscillation of solutions of f^″+Af=0$f^{″}+Af=0$ is discussed by focusing on four separate situations. In the complex case A$A$ is assumed to be either analytic in the unit disc D$\mathbb{D}$ or entire, while in the real case A$A$ is continuous either on (−1,1)$(-1,1)$ or on (0,∞)$(0,\infty )$. In all situations A$A$ is expected to grow beyond bounds that ensure finite oscillation for all (non-trivial) solutions, and the separation between distinct zeros of solutions is considered.     

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$.     

Large deviations of infinite intersections of events in Gaussian processes

Consider events of the form {_Zs≥ζ(s),s∈S}$\{Z_s\ge \zeta \left(s\right),s\in S\}$, where Z$Z$ is a continuous Gaussian process with stationary increments, ζ$\zeta$ is a function that belongs to the reproducing kernel Hilbert space R$R$ of process Z$Z$, and S⊂R$S\subset \mathbb{R}$ is compact. The main problem considered in this paper is identifying the function β^∗∈R${\beta }^{\ast }\in R$ satisfying β^∗(s)≥ζ(s)${\beta }^{\ast }\left(s\right)\ge \zeta \left(s\right)$ on S$S$ and having minimal R$R$-norm. The smoothness (mean square differentiability) of Z$Z$ turns out to have a crucial impact on the structure of the solution. As examples, we obtain the explicit solutions when ζ(s)=s$\zeta \left(s\right)=s$ for s∈[0,1]$s\in [0,1]$ and Z$Z$ is either a fractional Brownian motion or an integrated Ornstein–Uhlenbeck process.     

Approximation of solutions of Hammerstein equations with bounded strongly accretive nonlinear operators

Suppose X$X$ is a real q$q$-uniformly smooth Banach space and F,K:X→X$F,K:X\to X$ are bounded strongly accretive maps with D(K)=F(X)=X$D\left(K\right)=F\left(X\right)=X$. Let u^∗$u^{\ast }$ denote the unique solution of the Hammerstein equation u+KFu=0$u+KFu=0$. A new explicit coupled iteration process is shown to converge strongly to u^∗$u^{\ast }$. No invertibility assumption is imposed on K$K$ and the operators K$K$ and F$F$ need not be defined on compact subsets of X$X$. Furthermore, our new technique of proof is of independent interest. Finally, some interesting open questions are included.     

Asymptotic average shadowing property on compact metric spaces

We prove that if for a continuous map f$f$ on a compact metric space X$X$, the chain recurrent set, R(f)$R\left(f\right)$ has more than one chain component, then f$f$ does not satisfy the asymptotic average shadowing property. We also show that if a continuous map f$f$ on a compact metric space X$X$ has the asymptotic average shadowing property and if A$A$ is an attractor for f$f$, then A$A$ is the single attractor for f$f$ and we have A=R(f)$A=R\left(f\right)$. We also study diffeomorphisms with asymptotic average shadowing property and prove that if M$M$ is a compact manifold which is not finite with dimM=2$dimM=2$, then the C¹$C^1$ interior of the set of all C¹$C^1$ diffeomorphisms with the asymptotic average shadowing property is characterized by the set of Ω$\Omega$-stable diffeomorphisms.     

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.     

An upper bound for a valence of a face in a parallelohedral tiling

Consider a face-to-face parallelohedral tiling of R^d${\mathbb{R}}^d$ and a (d−k)$(d-k)$-dimensional face F$F$ of the tiling. We prove that the valence of F$F$ (i.e. the number of tiles containing F$F$ as a face) is not greater than 2^k$2^k$. If the tiling is affinely equivalent to a Voronoi tiling for some lattice (the so called Voronoi case), this gives a well-known upper bound for the number of vertices of a Delaunay k$k$-cell. Yet we emphasize that such an affine equivalence is not assumed in the proof.     

Convergence of invariant measures for singular stochastic diffusion equations

It is proved that the solutions to the singular stochastic p$p$-Laplace equation, p∈(1,2)$p\in (1,2)$ and the solutions to the stochastic fast diffusion equation with nonlinearity parameter r∈(0,1)$r\in (0,1)$ on a bounded open domain Λ⊂R^d$\Lambda \subset {\mathbb{R}}^d$ with Dirichlet boundary conditions are continuous in mean, uniformly in time, with respect to the parameters p$p$ and r$r$ respectively (in the Hilbert spaces L²(Λ)$L^2\left(\Lambda \right)$, H⁻¹(Λ)$H^{-1}\left(\Lambda \right)$ respectively). The highly singular limit case p=1$p=1$ is treated with the help of stochastic evolution variational inequalities, where P$\mathbb{P}$-a.s. convergence, uniformly in time, is established.     

Modules satisfying the weak Nakayama property

Let R$R$ be a commutative ring with identity. We will say that an R$R$-module M$M$ satisfies the weak Nakayama property, if IM=M$IM=M$, where I$I$ is an ideal of R$R$, implies that for any x∈M$x\in M$ there exists a∈I$a\in I$ such that (a−1)x=0$(a-1)x=0$. In this paper, we will study modules satisfying the weak Nakayama property. It is proved that if R$R$ is a local ring, then R$R$ is a Max ring if and only if J(R)$J\left(R\right)$, the Jacobson radical of R$R$, is T$T$-nilpotent if and only if every R$R$-module satisfies the weak Nakayama property.     

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.     

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).     

Ordinary and symbolic powers are Golod

Let S$S$ be a positively graded polynomial ring over a field of characteristic 0$0$, and I⊂S$I\subset S$ a proper graded ideal. In this note it is shown that S/I$S/I$ is Golod if ∂(I)²⊂I$\partial {\left(I\right)}^2\subset I$. Here ∂(I)$\partial \left(I\right)$ denotes the ideal generated by all the partial derivatives of elements of I$I$. We apply this result to find large classes of Golod ideals, including powers, symbolic powers, and saturations of ideals.     

On the effectiveness of a generalization of Miller’s primality theorem

Berrizbeitia and Olivieri showed in a recent paper that, for any integer r$r$, the notion of ω$\omega$-prime to base a$a$ leads to a primality test for numbers n≡1$n\equiv 1$ mod r$r$, that under the Extended Riemann Hypothesis (ERH) runs in polynomial time. They showed that the complexity of their test is at most the complexity of the Miller primality test (MPT), which is O((logn)^4+o(1))$O\left({(logn)}^{4+o\left(1\right)}\right)$. They conjectured that their test is more effective than the MPT if r$r$ is large.     

On truncated variation,upward truncated variation and downward truncated variation for diffusions

The truncated variation, TV^c${\text{TV}}^c$, is a fairly new concept introduced in ?ochowski (2008) [5]. Roughly speaking, given a càdlàg function f$f$, its truncated variation is “the total variation which does not pay attention to small changes of f$f$, below some threshold c>0$c>0$”. The very basic consequence of such approach is that contrary to the total variation, TV^c${\text{TV}}^c$ is always finite. This is appealing to the stochastic analysis where so-far large classes of processes, like semimartingales or diffusions, could not be studied with the total variation. Recently in ?ochowski (2011) [6], another characterization of TV^c${\text{TV}}^c$ has been found. Namely TV^c${\text{TV}}^c$ is the smallest possible total variation of a function which approximates f$f$ uniformly with accuracy c/2$c/2$. Due to these properties we envisage that TV^c${\text{TV}}^c$ might be a useful concept both in the theory and applications of stochastic processes.