第二章 主方程(Master equation)

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,第二章 主方程,(Master equation),这里我们研究概率分布随时间的演化。,随机过程：与时间有关的随机变量,(time-dependent random variable),我们只考虑仅有短程记忆的过程,马尔科夫过程,(Markov process),，该过程的时间演化方程就是主方程。,主方程是统计物理里最重要的方程之一，它几乎是普遍适用的，并广泛地被应用于化学，生物学，人口动力学，布朗运动，流体，半导体，金融等问题。,2.1,主方程的推导,（,I,）一般情形,对于随机变量,Y,的概率密度，将采用以下的记号来表示：,（随机变量,Y,在 时刻取 值的概率）,;,（随机变量,Y,在 时刻取 值，在 时刻取 值的联合概率）,;,（随机变量,Y,在 时刻取 值，在,时刻取 值，在 时刻取 值的联合概率）。,联合概率密度是正的,:,它们可以被约化,:,并且是归一化的,:,一般性质,:,不同时刻概率密度之间的关系（上式对,积分）：,随机变量与时间有关的矩,(,表征随机变量在不同时刻的值之间的相关,):,平稳过程：如果一个过程对一切,n,与,都有：,在平衡时，所有物理过程都是平稳的。对一个平稳过程，有：,条件概率：,而 只依赖于,-,时间差的绝对值。,（在 时刻取 值的随机变量,Y,，在 时刻取 值的概率）,;,它由如下恒等式来定义：,（,II,）马尔科夫过程,(Markov process),=,（固定 时，随机变量,Y,具有值 的联合概率密度）,联合条件概率密度：,对马尔科夫过程我们有,(,其中,),：,即,t,n,时刻取,y,n,的条件概率完全由,t,n-1,时刻,y,n-1,的值确定。马尔科夫过程完全由 和,转移概率,两个函数确定。例如：,对,y,2,积分，容易得到（,Chapman-Kolmogorov,方程）：,Chapman-Kolmogor,ov,方程的重要性：告诉我们对马尔科夫过程来说两个相继步骤的转移概率是两个单个步骤转移概率的,乘积,，而且相继的步骤是统计,独立,的。,(III),主方程,(Master equation),计算（*,1,）的时间导数我们必须考虑：,这里我们定义 是系统在时间间隔 内，从态,y,1,变到态,y,2,的,单位时间,的转变概率密度（转移率）。因此,在时间 内，从态,y,1,转变,到态,y,2,的概率密度为 ；,在时间,内,不转变,的概率密度为 。所以有：,（*,3),的时间导数为：,在时刻,t,1,+,（,是一个非常小的正数），由定义我们有（以连续变量为例）：,（*,1,）,当,=0,时，由（*,1,）得：（*,2,）,由（*,1-3,）我们发现：,（*,4,）,这就是主方程。,(IV),细致平衡和,Monte Carlo,模拟,为简单记这里我们考虑离散的情形，这时主方程可写为：,这和我们以前学过的统计物理里的刘维尔定理很相似。对平稳过程，我们有,因此对不同的平衡态有（这里我们略去了时间）：,这就是,细致平衡（,detailed balance,）,。对统计物理研究的很多系统而言，转移概率一般是不含时的，即与系统是否处于平衡态无关，因此我们一般有：。可以证明（见下），即使系统初始处于非平衡态时（这时概率密度函数与时间有关），经过足够长的时间后系统将逐渐进入平衡态，这是我们对系统进行,Monte Carlo,模拟,的理论基础。,练习：,考虑相对熵：。这里 是系统处于非平衡态的概率,密度函数，则是系统处于平衡态的概率密度函数。证明 和 ，其中等号仅当系统处于平衡态时成立。由此有：,(V),福克,-,普朗克,(Fokker-Planck),方程,由于转移概率 将随,的增大而迅速减小，我们把,WP,1,按,的幂次展开：,当,y,是一个,连续变量,，而且,y,的改变以,小跳跃,的方式发生时，我们可导出 的偏微分方程,-,福克,-,普朗克方程。先做变量代换,:,类似地,这里 是跳跃的大小。于是主方程变为,上式右边第一项和最后一项可消去，因此得到：,（,+,）,这就是,福克,-,普朗克,(Fokker-Planck),方程,。,其中 是,第,n,级跃变矩,：,2.2,马尔科夫链,(Markov chain),马尔科夫,链：是马尔科夫过程的一个例子，是在,离散,时刻出现的,离散,随机变量,Y,取值之间的转移。,设,Y,可取值,，基本时间间隔为,1,，从,t=0,到,t=1,我们有：,引入 我们可把上式改写为矩阵方程：,在,s,时刻，我们有：,P,(s),在,s,很大时的行为依赖于转移矩阵的结构。,若,Q,的某个幂次的全部元素都是正的（正则矩阵），则,P,(s),趋向唯一的确定的与初态无关的定态 ：,且易证明：,一个例子（雷克书,P.173,）：考虑两个罐子,A,和,B,，有三个红球和两个白球分配给它们，并总使得,A,中有两个球。共有下面三种位形：,位形间的转移为：无规则地从,A,和,B,中各取一个球进行交换。,转移矩阵为 且易知 是正则的,:,令 表示定态，由方程：可解出定态，结果为：,，与初态无关。,2.3,无规行走和扩散方程,考虑一个粒子在,x,轴上运动，且各步行走是统计独立的。设步长为,l,，步间时间为,，,n=0,，,1,，,2,为粒子的绝对位置，则有：,若粒子向左向右运动的概率均为,1/2,，则,原方程可简化为：,把上式写为求导的形式，我们有：,令 并在 为有限的条件下取极限 便可得到扩散方程,:,假定初始时刻 并引入,P,1,(x,t),对,x,的傅里叶变换（特征函数），,扩散方程可变为：,该方程的解为：，再取逆变换，可得：,这是粒子在,t=0,从,x=0,出发，到,t,时刻于,x,点找到它的概率。,一阶矩和二阶矩,：,一阶矩,：把扩散方程两边乘以,并对位置积分后，我们发现：,因此粒子的平均距离不随时间改变；,二阶矩,：把扩散方程两边乘以,并对位置积分后，我们发现：,这正是扩散过程的特征。,2.4,生灭过程，主方程的求解,生灭过程：在一个时刻只能进行一步转移。,我们这里处理一个可用生成函数严格求解的情形。,再假定生和灭的概率正比于现存细菌数，则有 和 ，上式两边同乘以 并化简，即得线性生灭过程的主方程：,考虑,t,时刻有,m,个细菌的一个群体：,在时间 内死亡一个细菌的概率为,在时间 内出生一个细菌的概率为,在时间 内细菌数目不变的概率为,在时间 内出生或死亡数超过,1,的概率为零,.,于是有,:,对依赖于,离散,随机变量的主方程的求解：生成函数法,生成函数（,characteristic function,）可写为：,对,z,求导后令,z-1,可得到随机变量,n,的各阶距：,一般地，我们有：,因此对一般的多项式函数,r(n),我们有：,由此有：,把以上表达式带入到主方程中我们有：,（*）,因此方程,(*),和主方程是等价的，我们只需解方程,(*),。,容易发现,(*),可由方程组,及 得到。,从,dF=0,我们发现,F(z,t)=C,2,，由,一般解为,设,t=0,时，细菌数目为,m,，则,故若 则,结果我们求得,2.5,离散平稳马尔科夫过程的普遍解,对离散平稳马尔科夫过程，,Chapman-Kolmogorov,方程变为：,这里 和,由概率和条件概率的定义我们还有,:,即,和,考虑离散随机变量 和离散时间 ，其中,n,和,是整数。这时我们得到了一个马尔科夫链，我们有：,这里 是系统处于,k,态时下一步跳到,n,态的条件概率，它包含了系统转移机制的一切必要信息。组成了矩阵,Q,的分量：,由（,2.2,）节我们并有：,转移矩阵,Q,lxl,转移矩阵,Q,一般不是对称阵，因而其左，右本征矢量不同。其左本征矢量问题可写为：,右本征矢量问题可写为：,其中,是方程,det|Q-I,|=0,的解。由以上两式可以证明：,正交归一性：即,Q,可以用其左，右本征矢展开：因此我们有,Q,至少有一本征值为,1,，且,若所有,由上可知,则对足够大的,s,方程 不成立，这不可能。,再由,P,=,PQ,及,i,X,i,=,X,i,Q,和,i,Y,i,=,QY,i,，易得,因,Y,i,构成完备本征矢，若所有,i,1,则对所有,i,均有,PY,i,=0,，这不可能。故存在,i,使得,最后由右本征矢方程两边取绝对值得：,对所有的,m,求和并考虑,Q,的归一性即得：,（注意雷克书中的证明有一步是错误的。）,对正则转移矩阵，若,Q,只有一个本征值 则,2.6,近似方法,-,展开,（,I,）简单例子：一维无规行走,考虑一个有边界条件的一维无规行走，其主方程为：,这里,-LnL,而且,L1,，因而系统大小,=2L+11,。我们引入：,x=n/L,并记,(x,t)=P,1,(n,t),主方程可改写为：,由于,1/L,是小量，我们可以把,对,1/L,展开：,情形,1,：,=,:,这时上式右边第一项,1/L,项消失。为简单记我们令,=1,并记,这样,重新标度后我们有：,这是扩散方程,(Fokker-Planck,方程,),。,情形,2,：,:,这时只用考虑主方程右边第一项。令,=t/L,我们有：,这是一个有向无规行走且,x(,),满足：,(II),一般情形,这里考虑连续时间和离散随机变量的主方程并假定转移率,W,与时间无关，这样主方程可写为：,类似于连续随机变量的情形我们可以定义跃变矩：,跃变矩是随机变量,n,的方程。我们一般感兴趣的是随机变量,n,及其各级矩的运动方程，这些已知的话系统的性质就基本清楚了。,的运动方程：在主方程两边乘以,n,并对,n,求和，在对右边第一项作交换,nm,后我们获得：,的运动方程：在主方程两边乘以,并对,n,求和，在对右边第一项作交换,nm,后我们获得：,因此不用解主方程，通过转移率,W(n,m),我们就可以得到系统的大量信息,。,近似：,W,对系统参量,的展开,对大系统，我们可以把,W,对表征系统大小的参量,做展开（因,1/,是一个小量），并将其带入到主方程中，获得一个近似的主方程，这个方程的解可能对系统的性质做出较好的描述。,在转移率中重要的参量是密度,m/,和步长,n=n-m,。因此我们把,W(m,n),展开为：,这里,f(,),是,的任意函数。对大,我们略去上式中的高阶项并带入到主方程中，得：,对大量独立客体的行为，根据中心极限定理我们知道 ，宽度,n,正比于 于是我们可以把,n,在其平均值附近展开：,其中 是,n,对其平均值 的偏移。,我们可以把主方程用,x,来表示，在,n,取值,n-n,+,n,内，,我们定义（这样,(t),显式地依赖于,t,）：,其中,于是我们有：和,主方程随之变为：,把上式右边第一项在 附近作泰勒展开,并重新标定时间,f(,)t=,后，主方程最终变为：,其中 和,在主方程里保留到 ，可得：,要满足上式，只须取,这里,在主方程里保留到,的零级项，可得关于概率密度,的,Fokker-Planck,方程：,由上式即可得扰动,x,的平均值和矩等的运动方程。对平稳过程，上述方程右端的系数与时间无关。如在,=0,有,，,并定义,及 和,，,上述,Fokker-Planck,方程可变为一个广义扩散方程：,主方程展开到,的 级项时，,W,的泰勒展开的第二项,1,将开始有贡献。,2.7,非线性生灭过程,-,马尔萨斯方程,对,线性,生灭过程：,对,非线性,生灭过程：我们假设社会成员间的竞争使得死亡率加大，因此死亡率中还有一个正比于其它个体密度的项 贡献，这里,是系统的大小。这样转移率变为：,由于在时间,内不转变的概率为 ，我们发现主方程可写为：,考虑,t,时刻有,m,个人的一个社会，我们有：,容易发现转移率为：,在,t,时刻个体平均数的方程为：,这个方程中一级矩的演化依赖于二级矩。为此把,n,在其平均值附近展开：,并保留到,级的项，我们发现：,此方程的解为,特点：,当 则 人口消亡；,当 则 人口趋于一个稳定值（,时趋于无穷大）。即当转移率改变时，存在一个从一个状态到另一状态的“