Related papers: Monte Carlo PDE Solvers for Nonlinear Radiative Bo…
Monte Carlo is famous for accepting model extensions and model refinements up to infinite dimension. However, this powerful incremental design is based on a premise which has severely limited its application so far: a state-variable can…
This work demonstrates algorithms to accurately compute solutions to thermal radiation transport problems using a reduced floating-point precision implementation of the Implicit Monte Carlo method. Several techniques falling into the…
We introduce a new family of numerical algorithms for approximating solutions of general high-dimensional semilinear parabolic partial differential equations at single space-time points. The algorithm is obtained through a delicate…
In recent work, Li et al.\ (Comm.\ Math.\ Sci., 7:81-107, 2009) developed a diffuse-domain method (DDM) for solving partial differential equations in complex, dynamic geometries with Dirichlet, Neumann, and Robin boundary conditions. The…
A method is proposed to handle the sign problem in the simulation of systems having indefinite or complex-valued measures. In general, this new approach, which is based on renormalisation blocking, is shown to yield statistical errors…
The Monte Carlo implementation of different approaches for diffractive scattering in $e - p$ collisions (resolved $\PO$, pQCD, soft color interactions) is described, with emphasis on the construction of the hadronic final state. Simple…
The multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to…
We develop efficient and accurate numerical methods to solve a class of shallow shell problems of the von Karman type. The governing equations form a fourth-order coupled system of nonlinear biharnomic equations for the transverse…
We introduce a neural-preconditioned iterative solver for Poisson equations with mixed boundary conditions. Typical Poisson discretizations yield large, ill-conditioned linear systems. Iterative solvers can be effective for these problems,…
We extend Monte Carlo samplers based on piecewise deterministic Markov processes (PDMP samplers) by formally defining different boundary conditions such as sticky floors, soft and hard walls and teleportation portals. This allows PDMP…
We study the continuum limit of a family of kinetic Monte Carlo models of crystal surface relaxation that includes both the solid-on-solid and discrete Gaussian models. With computational experiments and theoretical arguments we are able to…
We study reaction-diffusion equations in cylinders with possibly nonlinear diffusion and possibly nonlinear Neumann boundary conditions. We provide a geometric Poincar\'e-type inequality and classification results for stable solutions, and…
We provide a probabilistic representations of the solution of some semilinear hyperbolicand high-order PDEs based on branching diffusions. These representations pave theway for a Monte-Carlo approximation of the solution, thus bypassing the…
The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC…
Discrete time stochastic optimal control problems and Markov decision processes (MDPs), respectively, serve as fundamental models for problems that involve sequential decision making under uncertainty and as such constitute the theoretical…
Recently, a class of efficient spectral Monte-Carlo methods was developed in \cite{Feng2025ExponentiallyAS} for solving fractional Poisson equations. These methods fully consider the low regularity of the solution near boundaries and…
We present a general framework for performing "bold-line" diagrammatic Monte Carlo calculations based on using an analytical partial resummation as a starting point for stochastically summing all diagrams. As an example we solve the one…
We develop several deep learning algorithms for approximating families of parametric PDE solutions. The proposed algorithms approximate solutions together with their gradients, which in the context of mathematical finance means that the…
We consider a quasilinear equation given in the half-space, i.e. a so called boundary reaction problem. Our concerns are a geometric Poincar\'e inequality and, as a byproduct of this inequality, a result on the symmetry of low-dimensional…
The Poisson-Boltzmann equation (PBE) is an implicit solvent continuum model for calculating the electrostatic potential and energies of ionic solvated biomolecules. However, its numerical solution remains a significant challenge due strong…