Models of boundary behavior of particles diffusing between two concentrations

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).

[1]  W. Saslow.  Voltaic cells for physicists: Two surface pumps and an internal resistance , 1999 .

[2]  M. Muir.  Physical Chemistry , 1888, Nature.

[3]  A. Kainz,et al.  The analytic structure of the stationary kinetic boundary layer for Brownian particles near an absorbing wall , 1991 .

[4]  C. D. Levermore,et al.  the distribution of exit times for weakly colored noise , 1989 .

[5]  A. Hodgkin.  The conduction of the nervous impulse , 1964 .

[6]  U. M. Titulaer,et al.  The kinetic boundary layer around an absorbing sphere and the growth of small droplets , 1989 .

[7]  George H. Weiss,et al.  Equilibration in two chambers connected by a capillary , 2003 .

[8]  A. Hodgkin.  Ionic Currents Underlying Activity in the Giant Axon of the Squid , 1949 .

[9]  W. Im,et al.  A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels. , 2000, Biophysical journal.

[10]  J. Mccammon,et al.  Molecular dynamics with stochastic boundary conditions , 1982 .

[11]  BioMOCA: A Transport Monte Carlo Model for Ion Channels , 2003 .

[12]  B. Roux,et al.  Statistical mechanical equilibrium theory of selective ion channels. , 1999, Biophysical journal.

[13]  W. Im,et al.  Ions and counterions in a biological channel: a molecular dynamics simulation of OmpF porin from Escherichia coli in an explicit membrane with 1 M KCl aqueous salt solution. , 2002, Journal of molecular biology.

[14]  W. Saslow.  Consider a spherical battery... , 1994 .

[15]  G. Uhlenbeck,et al.  On the Theory of the Brownian Motion II , 1945 .

[16]  R. Eisenberg,et al.  Diffusion as a chemical reaction: Stochastic trajectories between fixed concentrations , 1995 .

[17]  R. W. Kelsall,et al.  The Monte Carlo method for semiconductor device simulation , 1995 .

[18]  M. Saraniti,et al.  Self-consistent particle-based simulations of three-dimensional ionic solutions , 2003 .

[19]  A. Berezhkovskii,et al.  Kinetics of escape through a small hole , 2002 .

[20]  Zeev Schuss,et al.  Theory and Applications of Stochastic Differential Equations , 1980 .

[21]  B. Eisenberg.  Ionic Channels in Biological Membranes: Natural , 1998 .

[22]  Jonathan V Selinger,et al.  The kinetic boundary layer for the Klein-Kramers equation; A new numerical approach , 1984 .

[23]  Martin H. Levinson.  A Thread across the Ocean: The Heroic Story of the Transatlantic Cable , 2003 .

[24]  Mark F. Schumaker.  Boundary conditions and trajectories of diffusion processes , 2002 .

[25]  S. Selberherr.  Analysis and simulation of semiconductor devices , 1984 .

[26]  Shin-Ho Chung,et al.  Reservoir boundaries in Brownian dynamics simulations of ion channels. , 2002, Biophysical Journal.

[27]  M. Karplus,et al.  Deformable stochastic boundaries in molecular dynamics , 1983 .

[28]  S. Orey,et al.  ANALYTICAL TREATMENT OF ONE‐DIMENSIONAL MARKOV PROCESSES , 1971 .

[29]  W. Ebeling,et al.  Generalizations of Onsagers Semiphenomenological Theory of Electrolytic Conductance , 1978 .

[30]  S. Karlin,et al.  A second course in stochastic processes , 1981 .

[31]  U. M. Titulaer.  A systematic solution procedure for the Fokker-Planck equation of a Brownian particle in the high-friction case , 1978 .

[32]  Shin-Ho Chung,et al.  Tests of continuum theories as models of ion channels. II. Poisson-Nernst-Planck theory versus brownian dynamics. , 2000, Biophysical journal.

[33]  M. A. Burschka,et al.  The kinetic boundary layer for the Fokker-Planck equation with absorbing boundary , 1981 .

[34]  B. Hille.  Ionic channels of excitable membranes , 2001 .

[35]  S. Chung,et al.  Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics. , 2000, Biophysical journal.

[36]  B. Eisenberg,et al.  Ion permeation and glutamate residues linked by Poisson-Nernst-Planck theory in L-type calcium channels. , 1998, Biophysical journal.

[37]  U. M. Titulaer,et al.  The moment method for boundary layer problems in Brownian motion theory , 1989 .

[38]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[39]  A. Hux.  MEASUREMENT OF CURRENT-VOLTAGE RELATIONS IN THE MEMBRANE OF THE GIANT AXON OF LOLIGO , 2006 .

[40]  C. Doering,et al.  Mean exit times for particles driven by weakly colored noise , 1989 .

[41]  Boaz Nadler,et al.  The Stationary Arrival Process of Independent Diffusers from a Continuum to an Absorbing Boundary Is Poissonian , 2001, SIAM J. Appl. Math..

[42]  Boaz Nadler,et al.  Connecting a Discrete Ionic Simulation to a Continuum , 2003, SIAM J. Appl. Math..

[43]  R. Eisenberg,et al.  Three-Dimensional Continuum Simulations of Ion Transport Through Biological Ion Channels: Effect of Charge Distribution in the Constriction Region of Porin , 2002 .

[44]  F. W. Wiegel,et al.  BROWNIAN MOTION IN A FIELD OF FORCE , 1986 .

[45]  E. J. Watson,et al.  A drop of ink falls from my pen...It comes to earth, I know not when , 1985 .