BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

G. A. Pavliotis

G. A. Pavliotis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPmath.NAcond-mat.stat-mechNumerical Analysismath.APmath.PRphysics.comp-phmath.SPmath.STnlin.CDStatistics Theory

Catalog footprint

What is connected

21works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Mobility estimation for Langevin dynamics using control variates

The scaling of the mobility of two-dimensional Langevin dynamics in a periodic potential as the friction vanishes is not well understood for non-separable potentials. Theoretical results are lacking, and numerical calculation of the mobility in the underdamped regime is challenging because the computational cost of standard Monte Carlo methods is inversely proportional to the friction coefficient, while deterministic methods are ill-conditioned. In this work, we propose a new variance-reduction method based on control variates for efficiently estimating the mobility of Langevin-type dynamics. We provide bounds on the bias and variance of the proposed estimator, and illustrate its efficacy through numerical experiments, first in simple one-dimensional settings and then for two-dimensional Langevin dynamics. Our results corroborate previous numerical evidence on the scaling of the mobility in the low friction regime for a simple non-separable potential.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Mean Field Fokker-Planck Equation with Nonlinear No-flux Boundary Conditions

We consider the mean field Fokker-Planck equation subject to nonlinear no-flux boundary conditions, which necessarily arise when subjecting a system of Brownian particles interacting via a pair potential in a bounded domain. With the additional presence of an external potential $V_1$, we show, by analysing the linearised Fokker-Planck operator, that the spectral properties of the equilibrium densities can differ considerably when compared with previous studies, e.g., with periodic boundary conditions. Amongst other mean field models of complex many-body particle systems, we present numerical experiments encompassing in a wide range of physical applications, including: generalised exponential models (Gaussian, Morse); a Kuramoto model, for noisy coupled oscillators; and an Onsager model for liquid crystals. We showcase our results by using the numerical methods developed in the pseudospectral collocation scheme 2DChebClass.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mean-field limits for interacting diffusions with colored noise: phase transitions and spectral numerical methods

In this paper we consider systems of weakly interacting particles driven by colored noise in a bistable potential, and we study the effect of the correlation time of the noise on the bifurcation diagram for the equilibrium states. We accomplish this by solving the corresponding McKean-Vlasov equation using a Hermite spectral method, and we verify our findings using Monte Carlo simulations of the particle system. We consider both Gaussian and non-Gaussian noise processes, and for each model of the noise we also study the behavior of the system in the small correlation time regime using perturbation theory. The spectral method that we develop in this paper can be used for solving linear and nonlinear, local and nonlocal (mean-field) Fokker-Planck equations, without requiring that they have a gradient structure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scaling limits for the generalized Langevin equation

In this paper, we study the diffusive limit of solutions to the generalized Langevin equation (GLE) in a periodic potential. Under the assumption of quasi-Markovianity, we obtain sharp longtime equilibration estimates for the GLE using techniques from the theory of hypocoercivity. We then prove asymptotic results for the effective diffusion coefficient in three limiting regimes: the short memory, the overdamped and the underdamped limits. Finally, we employ a recently developed spectral numerical method in order to calculate the effective diffusion coefficient for a wide range of (effective) friction coefficients, confirming our asymptotic results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Long-time behaviour and phase transitions for the McKean--Vlasov equation on the torus

We study the McKean-Vlasov equation \[ \partial_t \varrho= β^{-1} Δ\varrho + κ\nabla \cdot (\varrho \nabla (W \star \varrho)) \, , \] with periodic boundary conditions on the torus. We first study the global asymptotic stability of the homogeneous steady state. We then focus our attention on the stationary system, and prove the existence of nontrivial solutions branching from the homogeneous steady state, through possibly infinitely many bifurcations, under appropriate assumptions on the interaction potential. We also provide sufficient conditions for the existence of continuous and discontinuous phase transitions. Finally, we showcase these results by applying them to several examples of interaction potentials such as the noisy Kuramoto model for synchronisation, the Keller--Segel model for bacterial chemotaxis, and the noisy Hegselmann--Krausse model for opinion dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Brownian Motion in an N-scale periodic Potential

We study the problem of Brownian motion in a multiscale potential. The potential is assumed to have N+1 scales (i.e. N small scales and one macroscale) and to depend periodically on all the small scales. We show that for nonseparable potentials, i.e. potentials in which the microscales and the macroscale are fully coupled, the homogenized equation is an overdamped Langevin equation with multiplicative noise driven by the free energy, for which the detailed balance condition still holds. The calculation of the effective diffusion tensor requires the solution of a system of N coupled Poisson equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Noise-induced transitions in rugged energy landscapes

We consider the problem of an overdamped Brownian particle moving in multiscale potential with N + 1 characteristic length scales: the macroscale and N separated microscales. We show that the coarse-grained dynamics is given by an overdamped Langevin equation with respect to the free energy and with a space dependent diffusion tensor, the calculation of which requires the solution of N fully coupled Poisson equations. We study in detail the structure of the bifurcation diagram for one-dimensional problems and we show that the multiscale structure in the potential leads to hysteresis effects and to noise-induced transitions. Furthermore, we obtain an explicit formula for the effective diffusion coefficient for a self-similar separable potential and we investigate the limit of infinitely many small scales.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spectral methods for multiscale stochastic differential equations

This paper presents a new method for the solution of multiscale stochastic differential equations at the diffusive time scale. In contrast to averaging-based methods, e.g., the heterogeneous multiscale method (HMM) or the equation-free method, which rely on Monte Carlo simulations, in this paper we introduce a new numerical methodology that is based on a spectral method. In particular, we use an expansion in Hermite functions to approximate the solution of an appropriate Poisson equation, which is used in order to calculate the coefficients of the homogenized equation. Spectral convergence is proved under suitable assumptions. Numerical experiments corroborate the theory and illustrate the performance of the method. A comparison with the HMM and an application to singularly perturbed stochastic PDEs are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient numerical calculation of drift and diffusion coefficients in the diffusion approximation of kinetic equations

In this paper we study the diffusion approximation of a swarming model given by a system of interacting Langevin equations with nonlinear friction. The diffusion approximation requires the calculation of the drift and diffusion coefficients that are given as averages of solutions to appropriate Poisson equations. We present a new numerical method for computing these coefficients that is based on the calculation of the eigenvalues and eigenfunctions of a Schrödinger operator. These theoretical results are supported by numerical simulations showcasing the efficiency of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variance Reduction using Nonreversible Langevin Samplers

A standard approach to computing expectations with respect to a given target measure is to introduce an overdamped Langevin equation which is reversible with respect to the target distribution, and to approximate the expectation by a time-averaging estimator. As has been noted in recent papers, introducing an appropriately chosen nonreversible component to the dynamics is beneficial, both in terms of reducing the asymptotic variance and of speeding up convergence to the target distribution. In this paper we present a detailed study of the dependence of the asymptotic variance on the deviation from reversibility. Our theoretical findings are supported by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A multiscale analysis of diffusions on rapidly varying surfaces

Lateral diffusion of molecules on surfaces plays a very important role in various biological processes, including lipid transport across the cell membrane, synaptic transmission and other phenomena such as exo- and endocytosis, signal transduction, chemotaxis and cell growth. In many cases, the surfaces can possess spatial inhomogeneities and/or be rapidly changing shape. Using a generalisation of the model for a thermally excited Helfrich elastic membrane, we consider the problem of lateral diffusion on quasi-planar surfaces, possessing both spatial and temporal fluctuations. Using results from homogenisation theory, we show that, under the assumption of scale separation between the characteristic length and time scales of the membrane fluctuations and the characteristic scale of the diffusing particle, the lateral diffusion process can be well approximated by a Brownian motion on the plane with constant diffusion tensor $D$ which depends in a highly nonlinear way on the detailed properties of the surface. The effective diffusion tensor will depend on the relative scales of the spatial and temporal fluctuations and, for different scaling regimes, we prove the existence of a macroscopic limit in each case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A New Mode Reduction Strategy for the Generalized Kuramoto-Sivashinsky Equation

Consider the generalized Kuramoto-Sivashinsky (gKS) equation. It is a model prototype for a wide variety of physical systems, from flame-front propagation, and more general front propagation in reaction-diffusion systems, to interface motion of viscous film flows. Our aim is to develop a systematic and rigorous low-dimensional representation of the gKS equation. For this purpose, we approximate it by a renormalization group (RG) equation which is qualitatively characterized by rigorous error bounds. This formulation allows for a new stochastic mode reduction guaranteeing optimality in the sense of maximal information entropy. Herewith, noise is systematically added to the reduced gKS equation and gives a rigorous and analytical explanation for its origin. These new results would allow to reliably perform low-dimensional numerical computations by accounting for the neglected degrees of freedom in a systematic way. Moreover, the presented reduction strategy might also be useful in other applications where classical mode reduction approaches fail or are too complicated to be implemented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A new stochastic mode reduction strategy for dissipative systems

We present a new methodology for studying non-Hamiltonian nonlinear systems based on an information theoretic extension of a renormalization group technique using a modified maximum entropy principle. We obtain a rigorous dimensionally reduced description for such systems. The neglected degrees of freedom by this reduction are replaced by a systematically deefined stochastic process under a constraint on the second moment. This then forms the basis of a computationally efficient method. Numerical computations for the generalized Kuramoto-Sivashinsky equation sup- port our method and reveal that the long-time underlying stochastic process of the fast (unresolved) modes obeys a universal distribution which does not depend on the initial conditions and which we rigorously derive by the maximum entropy principle.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Additive noise effects in active nonlinear spatially extended systems

We examine the effects of pure additive noise on spatially extended systems with quadratic nonlinearities. We develop a general multiscale theory for such systems and apply it to the Kuramoto-Sivashinsky equation as a case study. We first focus on a regime close to the instability onset (primary bifurcation), where the system can be described by a single dominant mode. We show analytically that the resulting noise in the equation describing the amplitude of the dominant mode largely depends on the nature of the stochastic forcing. For a highly degenerate noise, in the sense that it is acting on the first stable mode only, the amplitude equation is dominated by a pure multiplicative noise, which in turn induces the dominant mode to undergo several critical state transitions and complex phenomena, including intermittency and stabilisation, as the noise strength is increased. The intermittent behaviour is characterised by a power-law probability density and the corresponding critical exponent is calculated rigorously by making use of the first-passage properties of the amplitude equation. On the other hand, when the noise is acting on the whole subspace of stable modes, the multiplicative noise is corrected by an additive-like term, with the eventual loss of any stabilised state. We also show that the stochastic forcing has no effect on the dominant mode dynamics when it is acting on the second stable mode. Finally, in a regime which is relatively far from the instability onset, so that there are two unstable modes, we observe numerically that when the noise is acting on the first stable mode, both dominant modes show noise-induced complex phenomena similar to the single-mode case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Corrections to Einstein's relation for Brownian motion in a tilted periodic potential

In this paper we revisit the problem of Brownian motion in a tilted periodic potential. We use homogenization theory to derive general formulas for the effective velocity and the effective diffusion tensor that are valid for arbitrary tilts. Furthermore, we obtain power series expansions for the velocity and the diffusion coefficient as functions of the external forcing. Thus, we provide systematic corrections to Einstein's formula and to linear response theory. Our theoretical results are supported by extensive numerical simulations. For our numerical experiments we use a novel spectral numerical method that leads to a very efficient and accurate calculation of the effective velocity and the effective diffusion tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exponential return to equilibrium for hypoelliptic quadratic systems

We study the problem of convergence to equilibrium for evolution equations associated to general quadratic operators. Quadratic operators are non-selfadjoint differential operators with complex-valued quadratic symbols. Under appropriate assumptions, a complete description of the spectrum of such operators is given and the exponential return to equilibrium with sharp estimates on the rate of convergence is proven. Some applications to the study of chains of oscillators and the generalized Langevin equation are given.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Noise induced state transitions, intermittency and universality in the noisy Kuramoto-Sivashinsky equation

We analyze the effect of pure additive noise on the long-time dynamics of the noisy Kuramoto-Sivashinsky (KS) equation in a regime close to the instability onset. We show that when the noise is highly degenerate, in the sense that it acts only on the first stable mode, the solution of the KS equation undergoes several transitions between different states, including a critical on-off intermittent state that is eventually stabilized as the noise strength is increased. Such noise-induced transitions can be completely characterized through critical exponents, obtaining that both the KS and the noisy Burgers equation belong to the same universality class. The results of our numerical investigations are explained rigorously using multiscale techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical methods for stochastic partial differential equations with multiples scales

A new method for solving numerically stochastic partial differential equations (SPDEs) with multiple scales is presented. The method combines a spectral method with the heterogeneous multiscale method (HMM) presented in [W. E, D. Liu, and E. Vanden-Eijnden, Comm. Pure Appl. Math., 58(11):1544--1585, 2005]. The class of problems that we consider are SPDEs with quadratic nonlinearities that were studied in [D. Blomker, M. Hairer, and G.A. Pavliotis, Nonlinearity, 20(7):1721--1744, 2007.] For such SPDEs an amplitude equation which describes the effective dynamics at long time scales can be rigorously derived for both advective and diffusive time scales. Our method, based on micro and macro solvers, allows to capture numerically the amplitude equation accurately at a cost independent of the small scales in the problem. Numerical experiments illustrate the behavior of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymptotic analysis for the generalized langevin equation

Various qualitative properties of solutions to the generalized Langevin equation (GLE) in a periodic or a confining potential are studied in this paper. We consider a class of quasi-Markovian GLEs, similar to the model that was introduced in \cite{EPR99}. Geometric ergodicity, a homogenization theorem (invariance principle), short time asymptotics and the white noise limit are studied. Our proofs are based on a careful analysis of a hypoelliptic operator which is the generator of an auxiliary Markov process. Systematic use of the recently developed theory of hypocoercivity \cite{Vil04HPI} is made.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymptotic Analysis of the Green-Kubo Formula

A detailed study of various distinguished limits of the Green-Kubo formula for the self-diffusion coefficient is presented in this paper. First, an alternative representation of the Green-Kubo formula in terms of the solution of a Poisson equation is derived when the microscopic dynamics is Markovian. Then, the techniques developed in \cite{golden2, AvelMajda91} are used to obtain a Stieltjes integral representation formula for the symmetric and antisymmetric parts of the diffusion tensor. The effect of irreversible microscopic dynamics on the diffusion coefficient is analyzed and various asymptotic limits of physical interest are studied. Several examples are presented that confirm the findings of our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiscale Modelling and Inverse Problems

The need to blend observational data and mathematical models arises in many applications and leads naturally to inverse problems. Parameters appearing in the model, such as constitutive tensors, initial conditions, boundary conditions, and forcing can be estimated on the basis of observed data. The resulting inverse problems are often ill-posed and some form of regularization is required. These notes discuss parameter estimation in situations where the unknown parameters vary across multiple scales. We illustrate the main ideas using a simple model for groundwater flow. We will highlight various approaches to regularization for inverse problems, including Tikhonov and Bayesian methods. We illustrate three ideas that arise when considering inverse problems in the multiscale context. The first idea is that the choice of space or set in which to seek the solution to the inverse problem is intimately related to whether a homogenized or full multiscale solution is required. This is a choice of regularization. The second idea is that, if a homogenized solution to the inverse problem is what is desired, then this can be recovered from carefully designed observations of the full multiscale system. The third idea is that the theory of homogenization can be used to improve the estimation of homogenized coefficients from multiscale data.

0 citations0 reviewsNo code linkedOpen access