
In this work, we show that uniform integrability is not a necessary condition
for central limit theorems (CLT) to hold for normalized multilevel Monte Carlo
estimators, and we provide near optimal weaker conditions under which the CLT
is achieved. In particular, if the variance decay rate dominates the
computational cost rate (i.e., $\beta> \gamma$), we prove that the CLT always
holds.

Finite volume methods are proposed for computing approximate pathwise
entropy/kinetic solutions to conservation laws with a rough path dependent flux
function. For a convex flux, it is demonstrated that rough path oscillations
may lead to "cancellations" in the solution. Making use of this property, we
show that for $\alpha$H{\"o}lder continuous rough paths the convergence rate
of the numerical methods can improve from $\mathcal{O}(\text{COST}^{\gamma})$,
for some $\gamma \in \left[\alpha/(128\alpha), \alpha/(106\alpha)\right]$,
with $\alpha\in (0, 1)$, to $\mathcal{O}(\text{COST}^{\min(1/4,\alpha/2)})$.
Numerical examples support the theoretical results.

We design and analyse the performance of a multilevel ensemble Kalman filter
method (MLEnKF) for filtering settings where the underlying statespace model
is an infinitedimensional spatiotemporal process. We consider underlying
models that needs to be simulated by numerical methods, with discretization in
both space and time. The multilevel Monte Carlo (MLMC) sampling strategy,
achieving variance reduction through pairwise coupling of ensemble particles on
neighboring resolutions, is used in the samplemoment step of MLEnKF to produce
an efficient hierarchical filtering method for spatiotemporal models. Under
sufficient regularity, MLEnKF is proven to be more efficient for weak
approximations than EnKF, asymptotically in the largeensemble and
finenumericalresolution limit. Numerical examples support our theoretical
findings.

We derive computable error estimates for finite element approximations of
linear elliptic partial differential equations (PDE) with rough stochastic
coefficients. In this setting, the exact solutions contain high frequency
content that standard a posteriori error estimates fail to capture. We propose
goaloriented estimates, based on local error indicators, for the pathwise
Galerkin and expected quadrature errors committed in standard, continuous,
piecewise linear finite element approximations. Derived using easily validated
assumptions, these novel estimates can be computed at a relatively low cost and
have applications to subsurface flow problems in geophysics where the
conductivities are assumed to have lognormal distributions with low regularity.
Our theory is supported by numerical experiments on test problems in one and
two dimensions.

This work embeds a multilevel Monte Carlo sampling strategy into the Monte
Carlo step of the ensemble Kalman filter (EnKF) in the setting of finite
dimensional signal evolution and noisy discretetime observations. The signal
dynamics is assumed to be governed by a stochastic differential equation (SDE),
and a hierarchy of time grids is introduced for multilevel numerical
integration of that SDE. The resulting multilevel EnKF is proved to
asymptotically outperform EnKF in terms of computational cost versus
approximation accuracy. The theoretical results are illustrated numerically.

A formal mean square error expansion (MSE) is derived for EulerMaruyama
numerical solutions of stochastic differential equations (SDE). The error
expansion is used to construct a pathwise a posteriori adaptive time stepping
EulerMaruyama method for numerical solutions of SDE, and the resulting method
is incorporated into a multilevel Monte Carlo (MLMC) method for weak
approximations of SDE. This gives an efficient MSE adaptive MLMC method for
handling a number of lowregularity approximation problems. In lowregularity
numerical example problems, the developed adaptive MLMC method is shown to
outperform the uniform time stepping MLMC method by orders of magnitude,
producing output whose error with high probability is bounded by TOL>0 at the
nearoptimal MLMC cost rate O(TOL^{2}log(TOL)^4).

The difference of the values of observables for the timeindependent
Schroedinger equation, with matrix valued potentials, and the values of
observables for ab initio BornOppenheimer molecular dynamics, of the ground
state, depends on the probability to be in excited states and the
electron/nuclei mass ratio. The paper first proves an error estimate (depending
on the electron/nuclei mass ratio and the probability to be in excited states)
for this difference of microcanonical observables, assuming that molecular
dynamics spacetime averages converge, with a rate related to the maximal
Lyapunov exponent. The error estimate is uniform in the number of particles and
the analysis does not assume a uniform lower bound on the spectral gap of the
electron operator and consequently the probability to be in excited states can
be large. A numerical method to determine the probability to be in excited
states is then presented, based on Ehrenfest molecular dynamics and stability
analysis of a perturbed eigenvalue problem.

BornOppenheimer dynamics is shown to provide an accurate approximation of
timeindependent Schr\"odinger observables for a molecular system with an
electron spectral gap, in the limit of large ratio of nuclei and electron
masses, without assuming that the nuclei are localized to vanishing domains.
The derivation, based on a Hamiltonian system interpretation of the
Schr\"odinger equation and stability of the corresponding HamiltonJacobi
equation, bypasses the usual separation of nuclei and electron wave functions,
includes caustic states and gives a different perspective on the
BornOppenheimer approximation, Schr\"odinger Hamiltonian systems and numerical
simulation in molecular dynamics modeling at constant energy microcanonical
ensembles.