A robust numerical algorithm for studying biomolecular transport processes
A robust numerical algorithm for studying biomolecular transport processes
复制标题
DOI:
10.1006/jtbi.2003.3200
复制
发表时间:
2003-04-21
影响因子:
2
通讯作者:
Elston, TC
中科院分区:
文献类型:
--
作者:
Wang, HY;Peskin, CS;Elston, TC
We present a numerical algorithm that is well suited for the study of biomolecular transport processes. In the algorithm a continuous Markov process is discretized as a jump process and the jump rates are derived from local solutions of the continuous system. Consequently, the algorithm has two advantages over standard numerical methods: (1) it preserves detailed balance for equilibrium processes, (2) it is able to handle discontinuous potentials. The formulation of the algorithm also allows us to calculate the effective diffusion coefficient or, equivalently, the randomness parameter. We provide several simple examples of how to implement the algorithm. All the MATLAB functions files needed to reproduce the results presented in the article are available from www.amath.une.edu/Faculty/telston/matlab_functions. (C) 2003 Elsevier Science Ltd. All rights reserved.