Models of boundary behavior of particles diffusing between two concentrations

Abstract

Flux between regions of different concentration occurs in nearly every device involving diffusion, whether an electrochemical cell, a bipolar transistor, or a protein channel in a biological membrane. Diffusion theory has calculated that flux since the time of Fick (1855), and the flux has been known to arise from the stochastic behavior of Brownian trajectories since the time of Einstein (1905), yet the mathematical description of the behavior of trajectories corresponding to different types of boundaries is not complete. We consider the trajectories of non-interacting particles diffusing in a finite region connecting two baths of fixed concentrations. Inside the region, the trajectories of diffusing particles are governed by the Langevin equation. At the interface between the region and the baths, trajectories are set by a control mechanism that modifies dynamics so the concentration of particles remains (nearly) constant. We analyze different models of controllers and derive equations for the time evolution and spatial distribution of particles inside the domain. Our analysis shows a distinct difference between the time evolution and the steady state concentrations. While the time evolution of the density is governed by an integral operator, the spatial distribution is governed by the familiar Fokker-Planck operator. The boundary conditions for the time dependent density depend on the model of the controller; however, this dependence disappears in the steady state, if the controller is of a renewal type. Renewal-type controllers, however, produce spurious boundary layers that can be catastrophic in simulations of charged particles, because even a tiny net charge can have global effects. The design of a non-renewal controller that maintains concentrations of non-interacting particles without creating spurious boundary layers at the interface requires the solution of the time-dependent Fokker-Planck equation with absorption of outgoing trajectories and a source of ingoing trajectories on the boundary (the so called albedo problem).

Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids

“…A stochastic analysis of the trajectories of atoms is possible 152,153 and a multidimensional ͑nonequilibrium͒ Fokker-Planck equation can be derived by analysis 59 or steepest descent arguments. 154 The multidimensional Fokker-Planck equation can be reduced 155 to the PNP equations by a closure procedure 59,156,157 with no more ͑than the considerable͒ arbitrariness of closures of equilibrium systems 32,158 ͑which do not form obviously convergent or uniformly convergent series, for example, and thus have unknown errors͒. The Fokker-Planck equation includes the continuity equation and so does its derivates, the drift diffusion or PNP equation.…”

“…If different scales are merged as in Eq. ͑18͒, the variables in the different partial differential equations may not be comparable or even have the same units ͑e.g., the concentrations of species and the distribution function of locations of the atoms of that species 59,157,166 21,22,76,78 as written, without change. The cost function is the macroscopic ͑hydrodynamic͒ part of energy.…”

Poisson-Nernst-Planck-Fermi theory for modeling biological ion channels

“…If we extend our theory by introducing local probabilities p i (r) = v i C i (r) that depend on location, in effect allowing probabilities to depend on location as in the theory of stochastic processes 47 (applied for example to ionic channels 11,48,49 ), the electrochemical potential can be generalized locally to…”

Stochastic chemical reactions in microdomains

“…To compute the mean and the variance of the number of molecules N(t) surviving in Ω at time t, we use here the renewal equation [4,3,20],…”