University of Cambridge > Talks.cam > Isaac Newton Institute Seminar Series > Geometric numerical methods for confined Langevin dynamics

Geometric numerical methods for confined Langevin dynamics

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact nobody.

SSD - Stochastic systems for anomalous diffusion

Weak-sense numerical methods for (underdamped) Langevin dynamics in bounded domains are constructed and analysed, with both their finite time convergence and convergence to ergodic limits being proved.  First-order methods are based on an Euler-type scheme interlaced with collisions with the boundary. To achieve second order,  composition schemes are derived based on decomposition of the generator into collisional drift,  impulse, and stochastic momentum evolution.  In a deterministic setting, this approach would typically lead to first-order approximation, even in symmetric compositions, but we find that the stochastic method can provide second-order weak approximation with a single force evaluation, both at finite times and in the ergodic limit.   We provide theoretical and numerical justification for this observation using model problems and compare and contrast the numerical performance of different choices of the ordering of the terms in the splitting scheme. The talk is based on a recent joint work with Ben Leimkuhler (Edinburgh) and Akash Sharma (Gothenburg).

This talk is part of the Isaac Newton Institute Seminar Series series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.

 

© 2006-2024 Talks.cam, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity