Related papers: Perturbative density functional methods for choles…
We consider the application of multilevel Monte Carlo methods to steady state Darcy flow in a random porous medium, described mathematically by elliptic partial differential equations with random coefficients. The levels in the multilevel…
We prove the existence of weak solutions to a viscoelastic phase separation problem in two space dimensions. The mathematical model consists of a Cahn-Hilliard-type equation for two-phase flows and the Peterlin-Navier-Stokes equations for…
McKean-Vlasov stochastic differential equations (MVSDEs) describe systems whose dynamics depend on both individual states and the population distribution, and they arise widely in neuroscience, finance, and epidemiology. In many…
We introduce a new path integral Monte Carlo method for investigating nonadiabatic systems in thermal equilibrium and demonstrate an approach to reducing stochastic error. We derive a general path integral expression for the partition…
There has been substantial interest in developing Markov chain Monte Carlo algorithms based on piecewise-deterministic Markov processes. However existing algorithms can only be used if the target distribution of interest is differentiable…
This article describes Monte-Carlo algorithms for charged systems using constrained updates for the electric field. The method is generalized to treat inhomogeneous dielectric media, electrolytes via the Poisson-Boltzmann equation and…
We consider numerical methods for the Poisson-Nernst-Planck-Cahn-Hilliard (PNPCH) equations with steric interactions. We propose a novel energy stable numerical scheme that respects mass conservation and positivity at the discrete level.…
We use the Multi Level Monte Carlo method to estimate uncertainties in a Henry-like salt water intrusion problem with a fracture. The flow is induced by the variation of the density of the fluid phase, which depends on the mass fraction of…
A technique for reducing the number of integrals in a Monte Carlo calculation is introduced. For integrations relying on classical or mean-field trajectories with local weighting functions, it is possible to integrate analytically at least…
We present a reduced-scaling auxiliary-field quantum Monte Carlo (AFQMC) framework designed for large molecular systems and ensembles, with or without coupling to optical cavities. Our approach leverages the natural block sparsity of…
Using thermodynamic and variational principles we examine a basic phase field model for a mixture of two incompressible fluids in strongly perforated domains. With the help of the multiple scale method with drift and our recently introduced…
In the present paper we consider the initial data, external force, viscosity coefficients, and heat conductivity coefficient as random data for the compressible Navier--Stokes--Fourier system. The Monte Carlo method, which is frequently…
Particle Markov Chain Monte Carlo (PMCMC) is a general computational approach to Bayesian inference for general state space models. Our article scales up PMCMC in terms of the number of observations and parameters by generating the…
Monte Carlo methods are widely used in particle physics to integrate and sample probability distributions (differential cross sections or decay rates) on multi-dimensional phase spaces. We present a Neural Network (NN) algorithm optimized…
Macroscopically heterogeneous materials, characterised mostly by comparable heterogeneity lengthscale and structural sizes, can no longer be modelled by deterministic approach instead. It is convenient to introduce stochastic approach with…
We study the static properties of cubic blue phases by numerically minimising the three-dimensional, Landau-de Gennes free energy for a cholesteric liquid crystal close to the isotropic-cholesteric phase transition. Thus we are able to…
We observe stationary random tessellations $X=\{\Xi_n\}_{n\ge1}$ in $\mathbb{R}^d$ through a convex sampling window $W$ that expands unboundedly and we determine the total $(k-1)$-volume of those $(k-1)$-dimensional manifold processes which…
Stochastic expansion-based methods of uncertainty quantification, such as polynomial chaos and separated representations, require basis functions orthogonal with respect to the density of random inputs. Many modern engineering problems…
Soft condensed matter structures often challenge us with complex many-body phenomena governed by collective modes spanning wide spatial and temporal domains. In order to successfully tackle such problems mesoscopic coarse-grained (CG)…
Doubts have been raised on the representation of chiral spin liquids exhibiting topological order in terms of projected entangled pair states (PEPSs). Here, starting from a simple spin-1/2 chiral frustrated Heisenberg model, we show that a…