We analyze the random fluctuations of several multiscale algorithms, such as the multiscale finite element method (MsFEM) and the finite element heterogeneous multiscale method (HMM), that have been developed to solve partial differential equations with highly heterogeneous coefficients. Such multiscale algorithms are often shown to correctly capture the homogenization limit when the highly oscillatory random medium is stationary and ergodic. This paper is concerned with the random fluctuations of the solution about the deterministic homogenization limit. We consider the simplified setting of the one-dimensional elliptic equation, where the theory of random fluctuations is well understood. We develop a fluctuation theory for the multiscale algorithms in the presence of random environments with short-range and long-range correlations. For a given mesh size h, we show that the fluctuations converge in distribution in the space of continuous paths to Gaussian processes as the correlation length ε→0. We next derive the limit of such Gaussian processes as h→0 and compare this limit with the distribution of the random fluctuations of the continuous model. When such limits agree, we conclude that the multiscale algorithm captures the random fluctuations accurately and passes the corrector test. This property serves as an interesting benchmark to assess the behavior of the multiscale algorithm in practical situations where the assumptions necessary for the theory of homogenization are not met. What we find is that the computationally more expensive methods MsFEM, and HMM with a choice of parameter δ=h, correctly capture the random fluctuations both for short-range and long-range oscillations in the medium. The less expensive method HMM with δ<h correctly captures the fluctuations for long-range oscillations and strongly amplifies their size in media with short-range oscillations. We present a modified scheme with an intermediate computational cost that captures the random fluctuations in all cases.

We study surface charges on a generic null boundary in three dimensional topological massive gravity (TMG). We construct the solution phase space which involves four independent functions over the two dimensional null boundary. One of these functions corresponds to the massive chiral propagating graviton mode of TMG. The other three correspond to three surface charges of the theory, two of which can always be made integrable, while the last one can become integrable only in the absence of the chiral massive graviton flux through the null boundary. As the null boundary symmetry algebra we obtain Heisenberg ⊕ Virasoro algebra with a central charge proportional to the gravitational Chern-Simons term of TMG. We also discuss that the flux of the chiral massive gravitons appears as the (Bondi) news through the null surface.

We consider the problem of the random fluctuations in the solutions to elliptic PDEs with highly oscillatory random coefficients. In our setting, as the correlation length of the fluctuations tends to zero, the heterogeneous solution converges to a deterministic solution obtained by averaging. When the Green’s function to the unperturbed operator is sufficiently singular (i.e., not square integrable locally), the leading corrector to the averaged solution may be either deterministic or random, or both in a sense we shall explain. Our main application is the solution of an elliptic problem with random Robin boundary condition that may be used to model diffusion of signaling molecules through a layer of cells into a bulk of extracellular medium. The problem is then described by an elliptic pseudo-differential operator (a Dirichlet-to-Neumann operator) on the boundary of the domain with random potential. In the physical setting of a three dimensional extracellular medium on top of a two-dimensional surface of cells forming a layer of epithelium, we show that the approximate corrector to averaging consists of a deterministic correction plus a Gaussian field of amplitude proportional to the correlation length of the random medium. The result is obtained under some assumptions on the four-point correlation function in the medium. We provide examples of such random media based on Gaussian and Poisson statistics.

We construct the boundary phase space in D-dimensional Einstein gravity with a generic given co-dimension one null surface N as the boundary. The associated boundary symmetry algebra is a semi-direct sum of diffeomorphisms of \mathcal{N} and Weyl rescalings. It is generated by D towers of surface charges that are generic functions over \mathcal{N}. These surface charges can be rendered integrable for appropriate slicings of the phase space, provided there is no graviton flux through \mathcal{N}. In one particular slicing of this type, the charge algebra is the direct sum of the Heisenberg algebra and diffeomorphisms of the transverse space, \mathcal{N}_v for any fixed value of the advanced time v. Finally, we introduce null surface expansion- and spin-memories, and discuss associated memory effects that encode the passage of gravitational waves through \mathcal{N}, imprinted in a change of the surface charges.

We consider the theory of correctors to homogenization in stationary transport equations with rapidly oscillating, random coefficients. Let ε << 1 be the ratio of the correlation length in the random medium to the overall distance of propagation. As ε↓0, we show that the heterogeneous transport solution is well-approximated by a homogeneous transport solution. We then show that the rescaled corrector converges in (probability) distribution and weakly in the space and velocity variables, to a Gaussian process as an application of a central limit result. The latter result requires strong assumptions on the statistical structure of randomness and is proved for random processes constructed by means of a Poisson point process.