Algorithms for Brownian dynamics across discontinuities

Abstract

The problem of mass diffusion in layered systems has relevance to applications in different scientific disciplines, e.g., chemistry, material science, soil science, and biomedical engineering. The mathematical challenge in these type of model systems is to match the solutions of the time-dependent diffusion equation in each layer, such that the boundary conditions at the interfaces between them are satisfied. As the number of layers increases, the solutions may become increasingly complicated. Here, we describe an alternative computational approach to multi-layer diffusion problems, which is based on the description of the overdamped Brownian motion of particles via the underdamped Langevin equation. In this approach, the probability distribution function is computed from the statistics of an ensemble of independent single particle trajectories. To allow for simulations of Langevin dynamics in layered systems, the numerical integrator must be supplemented with algorithms for the transitions across the discontinuous interfaces. Algorithms for three common types of discontinuities are presented: (i) A discontinuity in the friction coefficient, (ii) a semi-permeable membrane, and (iii) a step-function chemical potential. The general case of an interface where all three discontinuities are present (Kedem-Katchalsky boundary) is also discussed. We demonstrate the validity and accuracy of the derived algorithms by considering a simple two-layer model system and comparing the Langevin dynamics statistics with analytical solutions and alternative computational results.

arXiv commentRevised version. Accepted for publication in the Journal of Computational Physics

References in corpus (13)

• State-dependent diffusion: thermodynamic consistency and its path integral formulation
• A simple and effective Verlet-type algorithm for simulating Langevin dynamics
• Application of the G-JF Discrete-Time Thermostat for Fast and Accurate Molecular Simulations
• Semi-analytical solution of multilayer diffusion problems with time-varying boundary conditions and general interface conditions
• Langevin dynamics in inhomogeneous media: Re-examining the Itô-Stratonovich Dilemma
• Fluctuation-dissipation relation for systems with spatially varying friction
• Isothermal Langevin dynamics in systems with power-law spatially-dependent friction
• Comparison of Modern Langevin Integrators for Simulations of Coarse-Grained Polymer Melts
• Diffusion in heterogeneous discs and spheres: new closed-form expressions for exit times and homogenization formulae
• New homogenization approaches for stochastic transport through heterogeneous media
• The G-JF Thermostat for Accurate Configurational Sampling in Soft-Matter Simulations
• Application of underdamped Langevin dynamics simulations for the study of diffusion from a drug-eluting stent
• Monte Carlo Approaches for Simulating a Particle at a Diffusivity Interface and the "Ito--Stratonovich Dilemma"

Cited by in corpus (7)

• A probabilistic model of diffusion through a semi-permeable barrier
• Renewal equations for single-particle diffusion through a semipermeable interface
• A Langevin dynamics approach for multi-layer mass transfer problems
• Thermodynamics of a Brownian particle in a non-confining potential
• Geometrical optics of large deviations of Brownian motion in inhomogeneous media
• Diffusion with stochastic resetting screened by a semipermeable interface
• Uphill drift in the absence of current in single-file diffusion