A linear stochastic biharmonic heat equation: hitting probabilities

Consider the linear stochastic biharmonic heat equation on a d–dimen- sional torus ($$d=1,2,3$$ d = 1 , 2 , 3 ), driven by a space-time white noise and with periodic boundary conditions: $$\begin{aligned} \left( \frac{\partial }{\partial t}+(-\varDelta )^2\right) v(t,x)= \sigma \dot{W}(t,x),\ (t,x)\in (0,T]\times {\mathbb {T}}^d, \end{aligned}$$ ∂ ∂ t + ( - Δ ) 2 v ( t , x ) = σ W ˙ ( t , x ) , ( t , x ) ∈ ( 0 , T ] × T d , $$v(0,x)=v_0(x)$$ v ( 0 , x ) = v 0 ( x ) . We find the canonical pseudo-distance corresponding to the random field solution, therefore the precise description of the anisotropies of the process. We see that for $$d=2$$ d = 2 , they include a $$z(\log \tfrac{c}{z})^{1/2}$$ z ( log c z ) 1 / 2 term. Consider D independent copies of the random field solution to (0.1). Applying the criteria proved in Hinojosa-Calleja and Sanz-Solé (Stoch PDE Anal Comp 2021. 10.1007/s40072-021-00190-1), we establish upper and lower bounds for the probabilities that the path process hits bounded Borel sets.This yields results on the polarity of sets and on the Hausdorff dimension of the path process.

Approximation of SDE solutions using local asymptotic expansions

We develop an asymptotic expansion for small time of the density of the solution of a non-degenerate system of stochastic differential equations with smooth coefficients, and apply this to the stepwise approximation of solutions. The asymptotic expansion, which takes the form of a multivariate Edgeworth-type expansion, is obtained from the Kolmogorov forward equation using some standard PDE results. To generate one step of the approximation to the solution, we use a Cornish–Fisher type expansion derived from the Edgeworth expansion. To interpret the approximation generated in this way as a strong approximation we use couplings between the (normal) random variables used and the Brownian path driving the SDE. These couplings are constructed using techniques from optimal transport and Vaserstein metrics. The type of approximation so obtained may be regarded as intermediate between a conventional strong approximation and a weak approximation. In principle approximations of any order can be obtained, though for higher orders the algebra becomes very heavy. In order 1/2 the method gives the usual Euler approximation; in order 1 it gives a variant of the Milstein method, but which needs only normal variables to be generated. However the method is somewhat limited by the non-degeneracy requirement.