CTRW的化学主方程推导(上)

最新推荐文章于 2022-09-21 00:35:55 发布

IceelfLuo

最新推荐文章于 2022-09-21 00:35:55 发布

阅读量574

点赞数

分类专栏：算法文章标签：算法概率论

本文链接：https://blog.csdn.net/IceelfLuo/article/details/107891535

版权

本文详细探讨了连续时间随机行走(CTRW)模型下化学主方程的推导过程，涉及定义、内部反应等待时间和推广的化学主方程。通过对反应等待时间的分析，阐述了如何构建非同质的多维CTRW过程，并引入了概率密度函数，为后续的化学主方程建立奠定基础。

摘要由CSDN通过智能技术生成

连续时间随机行走(continuous time random walk)是一种常见的扩散模型，化学主方程也和扩散有关，所以思考CTRW模型下的化学主方程的推导。公式比较多，分两部分来写。

定义

$m_{s}$ ：代表不同物种的个数
$m_{r}$ ：代表不同反应的个数，物种与反应共同组成一个反应系统。
$S_{j}$ ：代表不同的物种， $j$ 的取值从1到 $m_{s}$
$n_{j}$ ：物种 $S_{j}$ 对应的粒子数
$\mathbf{n}=\left(n_{1}, \cdots, n_{m_{s}}\right)^{\mathrm{T}}$ ：粒子数目的状态向量
$r_{i j} \in \mathbb{N}\left(p_{i j} \in \mathbb{N}\right)$ ：经过反应 $i$ 后 $n_{j}$ 的变化量。
$s_{\bar{j}}=p_{i j}-r_{i j}$ ：化学计量系数，代表物种 $j$ 经过反应 $i$ 后的粒子数目变化
根据以上定义，我们可以把反应 $i$ 对于状态空间的影响表示为 $\sum_{j} r_{i j} S_{j} \rightarrow \sum_{j} p_{i j} S_{j}$

内部反应等待时间

单一反应 $i$ 可以由它的反应等待时间 $\tau_{i}^{r}$ 来刻画（其中上标 $r$ 是为了与全局的等待时间区分，特指反应的等待时间)， $\tau_{i}^{r}$ 的 PDF，我们记为 $\psi_{i}^{r}$ ，考虑 $\psi_{i}^{r}$ 的表达式：
定义： $h_{i}(\mathbf{n})$ 为反应 $i$ 可能发生的反应事件数目，依赖于状态 $n$
那么会有 $h_{i}(\mathbf{n})=\prod_{j} \frac{n_{j} !}{r_{i j} !\left(n_{j}-r_{i j}\right) !} \circ$
这 $h_{i}$ 个可能事件中，任意单一反应事件可以由内部反应等待时间 $\theta_{l}^{(i)}\left(i i d, l=1, \cdots, h_{i}(\mathbf{n})\right)$ 来刻画，其中 $\theta_{l}^{(i)}$ 的 PDF记为 $p_{i}$