BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Kristin Kirchner

Kristin Kirchner contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.PRmath.NAmath.APNumerical Analysismath.FAmath.STMethodologyStatistics Theory

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access