
Dense kernel matrices $\Theta \in \mathbb{R}^{N \times N}$ obtained from
point evaluations of a covariance function $G$ at locations $\{ x_{i} \}_{1
\leq i \leq N}$ arise in statistics, machine learning, and numerical analysis.
For covariance functions that are Green's functions of elliptic boundary value
problems and homogeneouslydistributed sampling points, we show how to identify
a subset $S \subset \{ 1 , \dots , N \}^2$, with $\# S = O ( N \log (N)
\log^{d} ( N /\epsilon ) )$, such that the zero fillin incomplete Cholesky
factorisation of the sparse matrix $\Theta_{ij} 1_{( i, j ) \in S}$ is an
$\epsilon$approximation of $\Theta$. This factorisation can provably be
obtained in complexity $O ( N \log( N ) \log^{d}( N /\epsilon) )$ in space and
$O ( N \log^{2}( N ) \log^{2d}( N /\epsilon) )$ in time; we further present
numerical evidence that $d$ can be taken to be the intrinsic dimension of the
data set rather than that of the ambient space. The algorithm only needs to
know the spatial configuration of the $x_{i}$ and does not require an analytic
representation of $G$. Furthermore, this factorization straightforwardly
provides an approximate sparse PCA with optimal rate of convergence in the
operator norm. Hence, by using only subsampling and the incomplete Cholesky
factorization, we obtain, at nearly linear complexity, the compression,
inversion and approximate PCA of a large class of covariance matrices. By
inverting the order of the Cholesky factorization we also obtain a solver for
elliptic PDE with complexity $O ( N \log^{d}( N /\epsilon) )$ in space and $O (
N \log^{2d}( N /\epsilon) )$ in time.

Probabilistic integration of a continuous dynamical system is a way of
systematically introducing model error, at scales no larger than errors
introduced by standard numerical discretisation, in order to enable thorough
exploration of possible responses of the system to inputs. It is thus a
potentially useful approach in a number of applications such as forward
uncertainty quantification, inverse problems, and data assimilation. We extend
the convergence analysis of probabilistic integrators for deterministic
ordinary differential equations, as proposed by Conrad et al.\ (\textit{Stat.\
Comput.}, 2017), to establish meansquare convergence in the uniform norm on
discrete or continuoustime solutions under relaxed regularity assumptions on
the driving vector fields and their induced flows. Specifically, we show that
randomised highorder integrators for globally Lipschitz flows and randomised
Euler integrators for dissipative vector fields with polynomiallybounded local
Lipschitz constants all have the same meansquare convergence rate as their
deterministic counterparts, provided that the variance of the integration noise
is not of higher order than the corresponding deterministic integrator. These
and similar results are proven for probabilistic integrators where the random
perturbations may be statedependent, nonGaussian, or noncentred random
variables.

A strong mode of a probability measure on a normed space $X$ can be defined
as a point $u$ such that the mass of the ball centred at $u$ uniformly
dominates the mass of all other balls in the smallradius limit. Helin and
Burger weakened this definition by considering only pairwise comparisons with
balls whose centres differ by vectors in a dense, proper linear subspace $E$ of
$X$, and posed the question of when these two types of modes coincide. We show
that, in a more general setting of metrisable vector spaces equipped with
measures that are finite on bounded sets, the density of $E$ and a uniformity
condition suffice for the equivalence of these two types of modes. We
accomplish this by introducing a new, intermediate type of mode. We also show
that these modes can be inequivalent if the uniformity condition fails. Our
results shed light on the relationships between among various notions of
maximum a posteriori estimator in nonparametric Bayesian inference.

We consider the use of randomised forward models and loglikelihoods within
the Bayesian approach to inverse problems. Such random approximations to the
exact forward model or loglikelihood arise naturally when a computationally
expensive model is approximated using a cheaper stochastic surrogate, as in
Gaussian process emulation (kriging), or in the field of probabilistic
numerical methods. We show that the Hellinger distance between the exact and
approximate Bayesian posteriors is bounded by moments of the difference between
the true and approximate loglikelihoods. Example applications of these
stability results are given for randomised misfit models in large data
applications and the probabilistic solution of ordinary differential equations.

Consider an infinite sequence $(U_n)_{n\in\mathbb{N}}$ of independent Cauchy
random variables, defined by a sequence $(\delta_n)_{n\in\mathbb{N}}$ of
location parameters and a sequence $(\gamma_n)_{n\in\mathbb{N}}$ of scale
parameters. Let $(W_n)_{n\in\mathbb{N}}$ be another infinite sequence of
independent Cauchy random variables defined by the same sequence of location
parameters and the sequence $(\sigma_n\gamma_n)_{n\in\mathbb{N}}$ of scale
parameters, with $\sigma_n\neq 0$ for all $n\in\mathbb{N}$. Using a result of
Kakutani on equivalence of countably infinite product measures, we show that
the laws of $(U_n)_{n\in\mathbb{N}}$ and $(W_n)_{n\in\mathbb{N}}$ are
equivalent if and only if the sequence $(\vert
\sigma_n\vert1)_{n\in\mathbb{N}}$ is squaresummable.

The Bayesian perspective on inverse problems has attracted much mathematical
attention in recent years. Particular attention has been paid to Bayesian
inverse problems (BIPs) in which the parameter to be inferred lies in an
infinitedimensional space, a typical example being a scalar or tensor field
coupled to some observed data via an ODE or PDE. This article gives an
introduction to the framework of wellposed BIPs in infinitedimensional
parameter spaces, as advocated by Stuart (Acta Numer. 19:451559, 2010) and
others. This framework has the advantage of ensuring uniformly wellposed
inference problems independently of the finitedimensional discretisation used
for numerical solution. Recently, this framework has been extended to the case
of a heavytailed prior measure in the family of stable distributions, such as
an infinitedimensional Cauchy distribution, for which polynomial moments are
infinite or undefined. It is shown that analogues of the KarhunenLo\`eve
expansion for squareintegrable random variables can be used to sample such
measures on quasiBanach spaces. Furthermore, under weaker regularity
assumptions than those used to date, the Bayesian posterior measure is shown to
depend Lipschitz continuously in the Hellinger and total variation metrics upon
perturbations of the misfit function and observed data.

This article extends the framework of Bayesian inverse problems in
infinitedimensional parameter spaces, as advocated by Stuart (Acta Numer.
19:451559, 2010) and others, to the case of a heavytailed prior measure in
the family of stable distributions, such as an infinitedimensional Cauchy
distribution, for which polynomial moments are infinite or undefined. It is
shown that analogues of the KarhunenLo\`eve expansion for squareintegrable
random variables can be used to sample such measures on quasiBanach spaces.
Furthermore, under weaker regularity assumptions than those used to date, the
Bayesian posterior measure is shown to depend Lipschitz continuously in the
Hellinger metric upon perturbations of the misfit function and observed data.

Given a sequence of Cauchydistributed random variables defined by a sequence
of location parameters and a sequence of scale parameters, we consider another
sequence of random variables that is obtained by perturbing the location or
scale parameter sequences. Using a result of Kakutani on equivalence of
infinite product measures, we provide sufficient conditions for the equivalence
of laws of the two sequences.

We consider the problem of providing optimal uncertainty quantification (UQ)
 and hence rigorous certification  for partiallyobserved functions. We
present a UQ framework within which the observations may be small or large in
number, and need not carry information about the probability distribution of
the system in operation. The UQ objectives are posed as optimization problems,
the solutions of which are optimal bounds on the quantities of interest; we
consider two typical settings, namely parameter sensitivities (McDiarmid
diameters) and output deviation (or failure) probabilities. The solutions of
these optimization problems depend nontrivially (even nonmonotonically and
discontinuously) upon the specified legacy data. Furthermore, the extreme
values are often determined by only a few members of the data set; in our
principal physicallymotivated example, the bounds are determined by just 2 out
of 32 data points, and the remainder carry no information and could be
neglected without changing the final answer. We propose an analogue of the
simplex algorithm from linear programming that uses these observations to offer
efficient and rigorous UQ for highdimensional systems with highcardinality
legacy data. These findings suggest natural methods for selecting optimal
(maximally informative) next experiments.

We propose a composite layered structure for tunable, lowloss plasmon
resonances, which con sists of a noblemetal thin film coated in graphene and
supported on a hexagonal boron nitride (hBN) substrate. We calculate electron
energy loss spectra (EELS) for these structures, and nu merically demonstrate
that bulk plasmon losses in noblemetal films can be significantly reduced, and
surface coupling enhanced, through the addition of a graphene coating and the
widebandgap hBN substrate. Silver films with a trilayer graphene coating and
hBN substrate demonstrated sur face plasmondominant spectral profiles for
metallic layers as thick as 34 nm. A continuedfraction expression for the
effective dielectric function, based on a specular reflection model which
includes boundary interactions, is used to systematically demonstrate plasmon
peak tunability for a variety of configurations. Variations include substrate,
plasmonic metal, and individual layer thickness for each material. Mesoscale
calculation of EELS is performed with individual layer dielectric functions as
input to the effective dielectric function calculation, from which the loss
spectra are directly determined.

We consider the effective behaviour of a rateindependent process when it is
placed in contact with a heat bath. The method used to "thermalize" the process
is an interiorpoint entropic regularization of the MoreauYosida incremental
formulation of the unperturbed process. It is shown that the heat bath destroys
the rate independence in a controlled and deterministic way, and that the
effective dynamics are those of a nonlinear gradient descent in the original
energetic potential with respect to a different and nontrivial effective
dissipation potential.