Table of Contents
Fetching ...

Neural network approaches for variance reduction in fluctuation formulas

Grigorios Pavliotis, Renato Spacek, Gabriel Stoltz, Urbain Vaes

TL;DR

The paper tackles variance reduction in transport-coefficient estimation, recasting estimators as inner products with Poisson-solution operators $G$ of the form $-\mathcal{L}G=g$. It introduces physics-informed neural networks (PINNs) to approximate the Poisson solutions and builds three control variate frameworks (forward, adjoint, and combined) within Green--Kubo and Half-Einstein estimators, plus a zero-cost adjoint variant. It provides rigorous bias and variance analyses and demonstrates substantial variance reduction in underdamped Langevin dynamics and multiscale SDEs, achieving favorable cost-variance trade-offs even when approximations are rough. The results suggest PINNs as a practical tool for variance reduction in moderately high-dimensional stochastic systems where exact deterministic Poisson solves are infeasible, enabling efficient computation of mobility and other transport coefficients via fluctuation formulas.

Abstract

We propose a method utilizing physics-informed neural networks (PINNs) to solve Poisson equations that serve as control variates in the computation of transport coefficients via fluctuation formulas, such as the Green--Kubo and generalized Einstein-like formulas. By leveraging approximate solutions to the Poisson equation constructed through neural networks, our approach significantly reduces the variance of the estimator at hand. We provide an extensive numerical analysis of the estimators and detail a methodology for training neural networks to solve these Poisson equations. The approximate solutions are then incorporated into Monte Carlo simulations as effective control variates, demonstrating the suitability of the method for moderately high-dimensional problems where fully deterministic solutions are computationally infeasible.

Neural network approaches for variance reduction in fluctuation formulas

TL;DR

The paper tackles variance reduction in transport-coefficient estimation, recasting estimators as inner products with Poisson-solution operators of the form . It introduces physics-informed neural networks (PINNs) to approximate the Poisson solutions and builds three control variate frameworks (forward, adjoint, and combined) within Green--Kubo and Half-Einstein estimators, plus a zero-cost adjoint variant. It provides rigorous bias and variance analyses and demonstrates substantial variance reduction in underdamped Langevin dynamics and multiscale SDEs, achieving favorable cost-variance trade-offs even when approximations are rough. The results suggest PINNs as a practical tool for variance reduction in moderately high-dimensional stochastic systems where exact deterministic Poisson solves are infeasible, enabling efficient computation of mobility and other transport coefficients via fluctuation formulas.

Abstract

We propose a method utilizing physics-informed neural networks (PINNs) to solve Poisson equations that serve as control variates in the computation of transport coefficients via fluctuation formulas, such as the Green--Kubo and generalized Einstein-like formulas. By leveraging approximate solutions to the Poisson equation constructed through neural networks, our approach significantly reduces the variance of the estimator at hand. We provide an extensive numerical analysis of the estimators and detail a methodology for training neural networks to solve these Poisson equations. The approximate solutions are then incorporated into Monte Carlo simulations as effective control variates, demonstrating the suitability of the method for moderately high-dimensional problems where fully deterministic solutions are computationally infeasible.
Paper Structure (43 sections, 16 theorems, 167 equations, 7 figures, 7 tables)

This paper contains 43 sections, 16 theorems, 167 equations, 7 figures, 7 tables.

Key Result

Corollary 2.1

Suppose that as:semigroup_decay holds. Then $\mathcal{L}^{-1}$ is a well-defined bounded operator on $L^2_0(\mu)$, which satisfies the following operator identity: Thus, Poisson equations of the form $-\mathcal{L} G = g$ for $g \in L^2_0(\mu)$ are well-posed and admit a unique solution in $L^2_0(\mu)$.

Figures (7)

  • Figure 4.1: Training iteration steps vs training loss for the 2D Langevin dynamics for both $\psi_g$ and $\psi_f^*$ approximate solutions. Top: linear scale. Bottom: logarithmic y-axis.
  • Figure 4.2: Asymptotic variance over time for all estimators for the two-dimensional Langevin dynamics. Left: Green--Kubo estimators. Right: half-Einstein estimators. For Green--Kubo, the asymptotic variance is defined as the variance divided by time, while the two coincide for half-Einstein.
  • Figure 4.3: Variances of the Green--Kubo estimators presented in \ref{['sec:cv_methodology']}, when these are used to calculate the coefficients of the homogenized equation associated to \ref{['eq:multiscale']}. These values were calculated based on $K = 1000$ realizations of each estimator.
  • Figure 4.4: Variances of the half-Einstein estimators presented in \ref{['sec:cv_methodology']}, when these are used to calculate the coefficients of the homogenized equation associated to \ref{['eq:multiscale']}. These values were calculated based on $K = 1000$ realizations of each estimator.
  • Figure C.1: Absolute bias and variance of the Green--Kubo estimator \ref{['eq:GK_estimator']}, in the simple setting of \ref{['appendix:weight_num_ill']}.
  • ...and 2 more figures

Theorems & Definitions (30)

  • Corollary 2.1: Well-posedness of Poisson equations
  • Proposition 2.2: Bounds on bias for the standard Green--Kubo estimator
  • proof
  • Proposition 2.3: Bounds on variance for the standard Green--Kubo estimator
  • proof
  • Proposition 2.4: Bounds on bias for the general half-Einstein estimator
  • Remark 2.5: Consistency conditions on $w$
  • proof : Proof of \ref{['prop:bias_gen_standard_HE']}
  • Proposition 2.6: Bounds on variance for the generalized half-Einstein estimator
  • proof
  • ...and 20 more