Source author record

Kristin Kirchner

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.PR math.NA math.AP Numerical Analysis math.FA math.ST Methodology Statistics Theory

Catalog footprint

Research overview

5works

8topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2023arXiv

Monte Carlo convergence rates for $k$ th moments in Banach spaces

We formulate standard and multilevel Monte Carlo methods for the $k$ th moment $\mathbb{M}^k_\varepsilon[ξ]$ of a Banach space valued random variable $ξ\colonΩ\to E$ , interpreted as an element of the $k$ -fold injective tensor product space $\otimes^k_\varepsilon E$ . For the standard Monte Carlo estimator of $\mathbb{M}^k_\varepsilon[ξ]$ , we prove the $k$ -independent convergence rate $1-\frac{1}{p}$ in the $L_q(Ω;\otimes^k_\varepsilon E)$ -norm, provided that (i) $ξ\in L_{kq}(Ω;E)$ and (ii) $q\in[p,\infty)$ , where $p\in[1,2]$ is the Rademacher type of $E$ . By using the fact that Rademacher averages are dominated by Gaussian sums combined with a version of Slepian's inequality for Gaussian processes due to Fernique, we moreover derive corresponding results for multilevel Monte Carlo methods, including a rigorous error estimate in the $L_q(Ω;\otimes^k_\varepsilon E)$ -norm and the optimization of the computational cost for a given accuracy. Whenever the type of the Banach space $E$ is $p=2$ , our findings coincide with known results for Hilbert space valued random variables. We illustrate the abstract results by three model problems: second-order elliptic PDEs with random forcing or random coefficient, and stochastic evolution equations. In these cases, the solution processes naturally take values in non-Hilbertian Banach spaces. Further applications, where physical modeling constraints impose a setting in Banach spaces of type $p<2$ , are indicated.

preprint2023arXiv

Multiple and weak Markov properties in Hilbert spaces with applications to fractional stochastic evolution equations

We define various higher-order Markov properties for stochastic processes $(X(t))_{t\in \mathbb{T}}$ , indexed by an interval $\mathbb{T} \subseteq \mathbb{R}$ and taking values in a real and separable Hilbert space $U$ . We furthermore investigate the relations between them. In particular, for solutions to the stochastic evolution equation $\mathcal{L} X = \dot W^Q\!$ , where $\mathcal{L}$ is a linear operator acting on functions mapping from $\mathbb{T}$ to $U$ and $(\dot W^Q(t))_{t\in\mathbb{T}}$ is the formal derivative of a $U$ -valued (cylindrical) $Q$ -Wiener process, we prove necessary and sufficient conditions for the weakest Markov property via locality of the precision operator $\mathcal{L}^*\! \mathcal{L}$ . As an application, we consider the space-time fractional parabolic operator $\mathcal{L} = (\partial_t + A)^γ$ of order $γ\in (1/2,\infty)$ , where $-A$ is a linear operator generating a $C_0$ -semigroup on $U$ . We prove that the resulting solution process satisfies an $N$ th order Markov property if $γ= N \in \mathbb{N}$ and show that a necessary condition for the weakest Markov property is generally not satisfied if $γ\notin \mathbb{N}$ . The relevance of this class of processes is twofold: Firstly, it can be seen as a spatiotemporal generalization of Whittle-Matérn Gaussian random fields if $U = L^2(\mathcal{D})$ for a spatial domain $\mathcal{D}\subseteq\mathbb{R}^d\!$ . Secondly, we show that a $U$ -valued analog to the fractional Brownian motion with Hurst parameter $H \in (0,1)$ can be obtained as the limiting case of $\mathcal{L} = (\partial_t + \varepsilon \, \mathrm{Id}_U)^{H + \frac{1}{2}}$ for $\varepsilon \downarrow 0$ .

preprint2023arXiv

Regularity theory for a new class of fractional parabolic stochastic evolution equations

A new class of fractional-order stochastic evolution equations of the form $(\partial_t + A)^γX(t) = \dot{W}^Q(t)$ , $t\in[0,T]$ , $γ\in (0,\infty)$ , is introduced, where $-A$ generates a $C_0$ -semigroup on a separable Hilbert space $H$ and the spatiotemporal driving noise $\dot{W}^Q$ is the formal time derivative of an $H$ -valued cylindrical $Q$ -Wiener process. Mild and weak solutions are defined; these concepts are shown to be equivalent and to lead to well-posed problems. Temporal and spatial regularity of the solution process $X$ are investigated, the former being measured by mean-square or pathwise smoothness and the latter by using domains of fractional powers of $A$ . In addition, the covariance of $X$ and its long-time behavior are analyzed. These abstract results are applied to the cases when $A := L^β$ and $Q:=\tilde{L}^{-α}$ are fractional powers of symmetric, strongly elliptic second-order differential operators defined on (i) bounded Euclidean domains or (ii) smooth, compact surfaces. In these cases, the Gaussian solution processes can be seen as generalizations of merely spatial (Whittle-)Matérn fields to space-time.

preprint2021arXiv

Multilevel approximation of Gaussian random fields: Covariance compression, estimation and spatial prediction

Centered Gaussian random fields (GRFs) indexed by compacta such as smooth, bounded Euclidean domains or smooth, compact and orientable manifolds are determined by their covariance operators. We consider centered GRFs given as variational solutions to coloring operator equations driven by spatial white noise, with an elliptic self-adjoint pseudodifferential coloring operator from the Hörmander class. This includes the Matérn class of GRFs as a special case. Using biorthogonal multiresolution analyses on the manifold, we prove that the precision and covariance operators, respectively, may be identified with bi-infinite matrices and finite sections may be diagonally preconditioned rendering the condition number independent of the dimension $p$ of this section. We prove that a tapering strategy by thresholding applied on finite sections of the bi-infinite precision and covariance matrices results in optimally numerically sparse approximations. That is, asymptotically only linearly many nonzero matrix entries are sufficient to approximate the original section of the bi-infinite covariance or precision matrix using this tapering strategy to arbitrary precision. The locations of these nonzero matrix entries are known a priori. The tapered covariance or precision matrices may also be optimally diagonally preconditioned. Analysis of the relative size of the entries of the tapered covariance matrices motivates novel, multilevel Monte Carlo (MLMC) oracles for covariance estimation, in sample complexity that scales log-linearly with respect to the number $p$ of parameters. In addition, we propose and analyze a novel compressive algorithm for simulating and kriging of GRFs. The complexity (work and memory vs. accuracy) of these three algorithms scales near-optimally in terms of the number of parameters $p$ of the sample-wise approximation of the GRF in Sobolev scales.

preprint2019arXiv

Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields

We analyze several Galerkin approximations of a Gaussian random field $\mathcal{Z}\colon\mathcal{D}\timesΩ\to\mathbb{R}$ indexed by a Euclidean domain $\mathcal{D}\subset\mathbb{R}^d$ whose covariance structure is determined by a negative fractional power $L^{-2β}$ of a second-order elliptic differential operator $L:= -\nabla\cdot(A\nabla) + κ^2$ . Under minimal assumptions on the domain $\mathcal{D}$ , the coefficients $A\colon\mathcal{D}\to\mathbb{R}^{d\times d}$ , $κ\colon\mathcal{D}\to\mathbb{R}$ , and the fractional exponent $β>0$ , we prove convergence in $L_q(Ω; H^σ(\mathcal{D}))$ and in $L_q(Ω; C^δ(\overline{\mathcal{D}}))$ at (essentially) optimal rates for (i) spectral Galerkin methods and (ii) finite element approximations. Specifically, our analysis is solely based on $H^{1+α}(\mathcal{D})$ -regularity of the differential operator $L$ , where $0<α\leq 1$ . For this setting, we furthermore provide rigorous estimates for the error in the covariance function of these approximations in $L_{\infty}(\mathcal{D}\times\mathcal{D})$ and in the mixed Sobolev space $H^{σ,σ}(\mathcal{D}\times\mathcal{D})$ , showing convergence which is more than twice as fast compared to the corresponding $L_q(Ω; H^σ(\mathcal{D}))$ -rate. For the well-known example of such Gaussian random fields, the original Whittle-Matérn class, where $L=-Δ+ κ^2$ and $κ\equiv \operatorname{const.}$ , we perform several numerical experiments which validate our theoretical results.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Kristin Kirchner

Research overview

Is this your work?

5 publications

Monte Carlo convergence rates for $k$ th moments in Banach spaces

Multiple and weak Markov properties in Hilbert spaces with applications to fractional stochastic evolution equations

Regularity theory for a new class of fractional parabolic stochastic evolution equations

Multilevel approximation of Gaussian random fields: Covariance compression, estimation and spatial prediction

Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields

See the researcher in context

Building this map preview

Kristin Kirchner

Research overview

Is this your work?

5 publications

Monte Carlo convergence rates for kkkth moments in Banach spaces

Multiple and weak Markov properties in Hilbert spaces with applications to fractional stochastic evolution equations

Regularity theory for a new class of fractional parabolic stochastic evolution equations

Multilevel approximation of Gaussian random fields: Covariance compression, estimation and spatial prediction

Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields

See the researcher in context

Building this map preview

Monte Carlo convergence rates for $k$ th moments in Banach spaces