-
Multilevel Monte Carlo methods for the Dean-Kawasaki equation from Fluctuating Hydrodynamics
Authors:
Federico Cornalba,
Julian Fischer
Abstract:
Stochastic PDEs of Fluctuating Hydrodynamics are a powerful tool for the description of fluctuations in many-particle systems. In this paper, we develop and analyze a Multilevel Monte Carlo (MLMC) scheme for the Dean--Kawasaki equation, a pivotal representative of this class of SPDEs. We prove analytically and demonstrate numerically that our MLMC scheme provides a significant reduction in computa…
▽ More
Stochastic PDEs of Fluctuating Hydrodynamics are a powerful tool for the description of fluctuations in many-particle systems. In this paper, we develop and analyze a Multilevel Monte Carlo (MLMC) scheme for the Dean--Kawasaki equation, a pivotal representative of this class of SPDEs. We prove analytically and demonstrate numerically that our MLMC scheme provides a significant reduction in computational cost (with respect to a standard Monte Carlo method) in the simulation of the Dean--Kawasaki equation. Specifically, we link this reduction in cost to having a sufficiently large average particle density, and show that sizeable cost reductions can be obtained even when we have solutions with regions of low density. Numerical simulations are provided in the two-dimensional case, confirming our theoretical predictions.
Our results are formulated entirely in terms of the law of distributions rather than in terms of strong spatial norms: this crucially allows for MLMC speed-ups altogether despite the Dean--Kawasaki equation being highly singular.
△ Less
Submitted 8 May, 2024; v1 submitted 15 November, 2023;
originally announced November 2023.
-
Density fluctuations in weakly interacting particle systems via the Dean-Kawasaki equation
Authors:
Federico Cornalba,
Julian Fischer,
Jonas Ingmanns,
Claudia Raithel
Abstract:
The Dean-Kawasaki equation - one of the most fundamental SPDEs of fluctuating hydrodynamics - has been proposed as a model for density fluctuations in weakly interacting particle systems. In its original form it is highly singular and fails to be renormalizable even by approaches such as regularity structures and paracontrolled distributions, hindering mathematical approaches to its rigorous justi…
▽ More
The Dean-Kawasaki equation - one of the most fundamental SPDEs of fluctuating hydrodynamics - has been proposed as a model for density fluctuations in weakly interacting particle systems. In its original form it is highly singular and fails to be renormalizable even by approaches such as regularity structures and paracontrolled distributions, hindering mathematical approaches to its rigorous justification. It has been understood recently that it is natural to introduce a suitable regularization, e.g., by applying a formal spatial discretization or by truncating high-frequency noise.
In the present work, we prove that a regularization in form of a formal discretization of the Dean-Kawasaki equation indeed accurately describes density fluctuations in systems of weakly interacting diffusing particles: We show that in suitable weak metrics, the law of fluctuations as predicted by the discretized Dean-Kawasaki SPDE approximates the law of fluctuations of the original particle system, up to an error that is of arbitrarily high order in the inverse particle number and a discretization error. In particular, the Dean-Kawasaki equation provides a means for efficient and accurate simulations of density fluctuations in weakly interacting particle systems.
△ Less
Submitted 21 November, 2023; v1 submitted 1 March, 2023;
originally announced March 2023.
-
The Regularised Inertial Dean-Kawasaki equation: discontinuous Galerkin approximation and modelling for low-density regime
Authors:
Federico Cornalba,
Tony Shardlow
Abstract:
The Regularised Inertial Dean-Kawasaki model (RIDK) -- introduced by the authors and J. Zimmer in earlier works -- is a nonlinear stochastic PDE capturing fluctuations around the mean-field limit for large-scale particle systems in both particle density and momentum density. We focus on the following two aspects. Firstly, we set up a Discontinuous Galerkin (DG) discretisation scheme for the RIDK m…
▽ More
The Regularised Inertial Dean-Kawasaki model (RIDK) -- introduced by the authors and J. Zimmer in earlier works -- is a nonlinear stochastic PDE capturing fluctuations around the mean-field limit for large-scale particle systems in both particle density and momentum density. We focus on the following two aspects. Firstly, we set up a Discontinuous Galerkin (DG) discretisation scheme for the RIDK model: we provide suitable definitions of numerical fluxes at the interface of the mesh elements which are consistent with the wave-type nature of the RIDK model and grant stability of the simulations, and we quantify the rate of convergence in mean square to the continuous RIDK model. Secondly, we introduce modifications of the RIDK model in order to preserve positivity of the density (such a feature only holds in a ''high-probability sense'' for the original RIDK model). By means of numerical simulations, we show that the modifications lead to physically realistic and positive density profiles. In one case, subject to additional regularity constraints, we also prove positivity. Finally, we present an application of our methodology to a system of diffusing and reacting particles. Our Python code is available in open-source format.
△ Less
Submitted 18 September, 2023; v1 submitted 20 July, 2022;
originally announced July 2022.
-
Multi-Objective reward generalization: Improving performance of Deep Reinforcement Learning for applications in single-asset trading
Authors:
Federico Cornalba,
Constantin Disselkamp,
Davide Scassola,
Christopher Helf
Abstract:
We investigate the potential of Multi-Objective, Deep Reinforcement Learning for stock and cryptocurrency single-asset trading: in particular, we consider a Multi-Objective algorithm which generalizes the reward functions and discount factor (i.e., these components are not specified a priori, but incorporated in the learning process). Firstly, using several important assets (cryptocurrency pairs B…
▽ More
We investigate the potential of Multi-Objective, Deep Reinforcement Learning for stock and cryptocurrency single-asset trading: in particular, we consider a Multi-Objective algorithm which generalizes the reward functions and discount factor (i.e., these components are not specified a priori, but incorporated in the learning process). Firstly, using several important assets (cryptocurrency pairs BTCUSD, ETHUSDT, XRPUSDT, and stock indexes AAPL, SPY, NIFTY50), we verify the reward generalization property of the proposed Multi-Objective algorithm, and provide preliminary statistical evidence showing increased predictive stability over the corresponding Single-Objective strategy. Secondly, we show that the Multi-Objective algorithm has a clear edge over the corresponding Single-Objective strategy when the reward mechanism is sparse (i.e., when non-null feedback is infrequent over time). Finally, we discuss the generalization properties with respect to the discount factor. The entirety of our code is provided in open source format.
△ Less
Submitted 17 February, 2023; v1 submitted 9 March, 2022;
originally announced March 2022.
-
The Dean-Kawasaki equation and the structure of density fluctuations in systems of diffusing particles
Authors:
Federico Cornalba,
Julian Fischer
Abstract:
The Dean-Kawasaki equation - a strongly singular SPDE - is a basic equation of fluctuating hydrodynamics; it has been proposed in the physics literature to describe the fluctuations of the density of $N$ independent diffusing particles in the regime of large particle numbers $N\gg 1$. The singular nature of the Dean-Kawasaki equation presents a substantial challenge for both its analysis and its r…
▽ More
The Dean-Kawasaki equation - a strongly singular SPDE - is a basic equation of fluctuating hydrodynamics; it has been proposed in the physics literature to describe the fluctuations of the density of $N$ independent diffusing particles in the regime of large particle numbers $N\gg 1$. The singular nature of the Dean-Kawasaki equation presents a substantial challenge for both its analysis and its rigorous mathematical justification. Besides being non-renormalisable by the theory of regularity structures by Hairer et al., it has recently been shown to not even admit nontrivial martingale solutions.
In the present work, we give a rigorous and fully quantitative justification of the Dean-Kawasaki equation by considering the natural regularisation provided by standard numerical discretisations: We show that structure-preserving discretisations of the Dean-Kawasaki equation may approximate the density fluctuations of $N$ non-interacting diffusing particles to arbitrary order in $N^{-1}$ (in suitable weak metrics). In other words, the Dean-Kawasaki equation may be interpreted as a "recipe" for accurate and efficient numerical simulations of the density fluctuations of independent diffusing particles.
△ Less
Submitted 5 July, 2023; v1 submitted 14 September, 2021;
originally announced September 2021.
-
Well-posedness for a regularised inertial Dean-Kawasaki model for slender particles in several space dimensions
Authors:
Federico Cornalba,
Tony Shardlow,
Johannes Zimmer
Abstract:
A stochastic PDE, describing mesoscopic fluctuations in systems of weakly interacting inertial particles of finite volume, is proposed and analysed in any finite dimension $d\in\mathbb{N}$. It is a regularised and inertial version of the Dean-Kawasaki model. A high-probability well-posedness theory for this model is developed. This theory improves significantly on the spatial scaling restrictions…
▽ More
A stochastic PDE, describing mesoscopic fluctuations in systems of weakly interacting inertial particles of finite volume, is proposed and analysed in any finite dimension $d\in\mathbb{N}$. It is a regularised and inertial version of the Dean-Kawasaki model. A high-probability well-posedness theory for this model is developed. This theory improves significantly on the spatial scaling restrictions imposed in an earlier work of the same authors, which applied only to significantly larger particles in one dimension. The well-posedness theory now applies in $d$-dimensions when the particle-width $ε$ is proportional to $N^{-1/θ}$ for $θ>2d$ and $N$ is the number of particles. This scaling is optimal in a certain Sobolev norm. Key tools of the analysis are fractional Sobolev spaces, sharp bounds on Bessel functions, separability of the regularisation in the $d$-spatial dimensions, and use of the Faà di Bruno's formula.
△ Less
Submitted 8 February, 2021; v1 submitted 12 May, 2020;
originally announced May 2020.
-
A priori positivity of solutions to a non-conservative stochastic thin-film equation
Authors:
Federico Cornalba
Abstract:
Stochastic conservation laws are often challenging when it comes to proving existence of non-negative solutions. In a recent work by J. Fischer and G. Grün (2018, Existence of positive solutions to stochastic thin-film equations, SIAM J. Math. Anal.), existence of positive martingale solutions to a conservative stochastic thin-film equation is established in the case of quadratic mobility. In this…
▽ More
Stochastic conservation laws are often challenging when it comes to proving existence of non-negative solutions. In a recent work by J. Fischer and G. Grün (2018, Existence of positive solutions to stochastic thin-film equations, SIAM J. Math. Anal.), existence of positive martingale solutions to a conservative stochastic thin-film equation is established in the case of quadratic mobility. In this work, we focus on a larger class of mobilities (including the linear one) for the thin-film model. In order to do so, we need to introduce nonlinear source potentials, thus obtaining a non-conservative version of the thin-film equation. For this model, we assume the existence of a sufficiently regular local solution (i.e., defined up to a stopping time $τ$) and, by providing suitable conditions on the source potentials and the noise, we prove that such solution can be extended up to any $T>0$ and that it is positive with probability one. A thorough comparison with the aforementioned reference work is provided.
△ Less
Submitted 19 November, 2018;
originally announced November 2018.
-
From weakly interacting particles to a regularised Dean--Kawasaki model
Authors:
Federico Cornalba,
Tony Shardlow,
Johannes Zimmer
Abstract:
The evolution of finitely many particles obeying Langevin dynamics is described by Dean-Kawasaki equations, a class of stochastic equations featuring a non-Lipschitz multiplicative noise in divergence form. We derive a regularised Dean-Kawasaki model based on second order Langevin dynamics by analysing a system of particles interacting via a pairwise potential. Key tools of our analysis are the pr…
▽ More
The evolution of finitely many particles obeying Langevin dynamics is described by Dean-Kawasaki equations, a class of stochastic equations featuring a non-Lipschitz multiplicative noise in divergence form. We derive a regularised Dean-Kawasaki model based on second order Langevin dynamics by analysing a system of particles interacting via a pairwise potential. Key tools of our analysis are the propagation of chaos and Simon's compactness criterion. The model we obtain is a small-noise stochastic perturbation of the undamped McKean-Vlasov equation. We also provide a high-probability result for existence and uniqueness for our model.
△ Less
Submitted 9 September, 2019; v1 submitted 15 November, 2018;
originally announced November 2018.
-
A regularised Dean-Kawasaki model: derivation and analysis
Authors:
Federico Cornalba,
Tony Shardlow,
Johannes Zimmer
Abstract:
The Dean-Kawasaki model consists of a nonlinear stochastic partial differential equation featuring a conservative, multiplicative, stochastic term with non-Lipschitz coefficient, and driven by space-time white noise; this equation describes the evolution of the density function for a system of finitely many particles governed by Langevin dynamics. Well-posedness for the Dean-Kawasaki model is open…
▽ More
The Dean-Kawasaki model consists of a nonlinear stochastic partial differential equation featuring a conservative, multiplicative, stochastic term with non-Lipschitz coefficient, and driven by space-time white noise; this equation describes the evolution of the density function for a system of finitely many particles governed by Langevin dynamics. Well-posedness for the Dean-Kawasaki model is open except for specific diffusive cases, corresponding to overdamped Langevin dynamics. There, it was recently shown by Lehmann, Konarovskyi, and von Renesse that no regular (non-atomic) solutions exist. We derive and analyse a suitably regularised Dean-Kawasaki model of wave equation type driven by coloured noise, corresponding to second order Langevin dynamics, in one space dimension. The regularisation can be interpreted as considering particles of finite size rather than describing them by atomic measures. We establish existence and uniqueness of a solution. Specifically, we prove a high-probability result for the existence and uniqueness of mild solutions to this regularised Dean-Kawasaki model.
△ Less
Submitted 21 January, 2019; v1 submitted 5 February, 2018;
originally announced February 2018.
-
A nonlocal stochastic Cahn-Hilliard equation
Authors:
Federico Cornalba
Abstract:
We consider a stochastic extension of the nonlocal convective Cahn-Hilliard equation containing an additive Wiener process noise. We first introduce a suitable analytical setting and make some mathematical and physical assumptions. We then establish, in a variational context, the existence of a weak statistical solution for this problem. Finally we prove existence and uniqueness of a strong soluti…
▽ More
We consider a stochastic extension of the nonlocal convective Cahn-Hilliard equation containing an additive Wiener process noise. We first introduce a suitable analytical setting and make some mathematical and physical assumptions. We then establish, in a variational context, the existence of a weak statistical solution for this problem. Finally we prove existence and uniqueness of a strong solution.
△ Less
Submitted 7 March, 2016; v1 submitted 27 October, 2015;
originally announced October 2015.