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

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

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

定义

m s m_{s} ms:代表不同物种的个数
m r m_{r} mr:代表不同反应的个数,物种与反应共同组成一个反应系统。
S j S_{j} Sj:代表不同的物种, j j j的取值从1到 m s m_{s} ms
n j n_{j} nj:物种 S j S_{j} Sj对应的粒子数
n = ( n 1 , ⋯   , n m s ) T \mathbf{n}=\left(n_{1}, \cdots, n_{m_{s}}\right)^{\mathrm{T}} n=(n1,,nms)T:粒子数目的状态向量
r i j ∈ N ( p i j ∈ N ) r_{i j} \in \mathbb{N}\left(p_{i j} \in \mathbb{N}\right) rijN(pijN):经过反应 i i i n j n_{j} nj的变化量。
s j ˉ = p i j − r i j s_{\bar{j}}=p_{i j}-r_{i j} sjˉ=pijrij:化学计量系数,代表物种 j j j经过反应 i i i后的粒子数目变化
根据以上定义,我们可以把反应 i i i对于状态空间的影响表示为 ∑ j r i j S j → ∑ j p i j S j \sum_{j} r_{i j} S_{j} \rightarrow \sum_{j} p_{i j} S_{j} jrijSjjpijSj

内部反应等待时间

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值