Wednesday, August 20, 2014

brownian01 方均位移写作〈x 2(t)〉= 2Dt,这里D=l 2/ (2τ)。高斯分布函数解的前两次 距 欠扩散: 长的等待时间,分数 FokkerPlanck方程; 超扩散,


在爱因斯坦100年前创立的无规行走扩散理论中,行走是一个不连续过程,即假设在每个 节拍振荡时间内,粒子随机跳跃一次。而连续时间无规行走引入了两次成功跳跃之间有一个 等待时间分布,且等待时间分布的一次距和跳跃长度分布的二次距中有一个是发散的。其揭 示欠扩散对应于长的等待时间,由分数 FokkerPlanck方程来描写;而超扩散则意味着粒子存 在长跳跃的物理机制,Levy飞行就属于此类,这种运动也可以通过长尾分布的非高斯噪声激 励来产生。布朗运动的高斯解是中心极限定理的必然结果。便于数值模拟有势中的反常输运 的当属含非Ohmic记忆摩擦的广义 Langevin方程。系统长时间动力学行为取决于噪声功率 谱的低频分布,因此当热色噪声的低频被滤掉,导致系统的有效摩擦消失,自由粒子呈现热扩 散的极限情形:弹道扩散。在这一过程中,Kubo第二涨落耗散定理被满足,但系统的稳定态部 分依赖于它的初始速度分布,各态历经性和平衡态被破坏。

分数布朗运动和反常扩散
包景东 (北京师范大学物理系,北京 100875)




下面以各向同性情形为例,讨论黑体辐射场。
1,Stefan-Boltzmann定律:(总)辐射通量B(T)为
40
()()d4cBTTνTρνσ+∞
=∫=,                        (2)
其中σ =5.67×10-5 erg cm-2 s-1 K-4称为Stefan-Boltzmann常数。易于得知,辐射场(总)能量密度ρ(T)=aT 4,其中a=4σ/c=7.57×10-15 erg cm-3 K-4为辐射密度常数。
2,Wien位移定律:令dρλ(T)/dλ=0,可以求得谱能量密度最高的光子的所对应的波长λmax。它与T的关系为
λmax T = 0.29 cm K,                                (3)
其中式右的常数称为Wien常数。量ρλ(T) = ρν (T)|ν=c/λ是以波长表征的能量密度。另外,我们也可以根据dρν(T)/dν=0得到谱能量密度最高的光子的所对应的频率νmax。请问:是否存在νmax = c/λmax?(答案是否定的。提示:|ρν dν| = |ρλdλ|,dν = -(c/λ2)dλ。)
摘 要:  本文评述了分数布朗运动和反常扩散现象及描写它们的几种数学方式。报告了我们 在弹道扩散的产生条件、起源和长时间效应方面的工作。 关键词:分数布朗运动;反常扩散;连续时间无规行走;广义Langevin方程;弹道扩散 中图分类号:O415.6       文献标识码:A 0 引言 历史上,布朗运动现象追源于JanIngenhausz在1785年观测木炭粉末在酒表面的运动;后 来RobertBrown在1828年观测花粉、尘埃、烟灰在水表面的扩散;1905年AlbertEinstein用随 机行走解释流体分子的热运动;1926年 JeanPerrin因为测量 Avogadro常数获 Nobel物理奖。 整整一百年前,伟大的物理学家爱因斯坦在一篇开创性的文章[1],建立了布朗运动的扩散理 论。他考虑一个布朗粒子的一维无规运动,粒子每隔τ时间被撞击一次而移动距离 l,每次撞 击时向左和向右移动的可能性各占一半。假设粒子从原点出发,在时刻 t,粒子已受到了 n=t/ τ次撞击。爱因斯坦证得:粒子的平均位移为零,〈x(t)〉=0;方均位移写作〈x 2(t)〉= 2Dt,这里D=l 2/ (2τ)。其实这是扩散方程 t P(x ,t )=D 2 xP(x ,t)的高斯分布函数解的前两次 距,也是中心极限定理的直接结果


上世纪40年代以后,布朗运动的动力学研究开始兴起[2, 3]。两个等价的方程是 Langevin 方程和FokkerPlanck方程 [4],前者是关于粒子轨道的随机微分方程;后者是关于分布函数随 时演化的二阶偏微分方程。人们通常认为前者更基本,因为在有些情况下写不出精确的 FokkerPlanck方程,例如色噪声和记忆阻尼。Kubo认为布朗运动可以看成一个系统受到外界 随机扰动的响应,提出了第一和第二涨落耗散定理[5]。站在更广义的角度上,用纯粹的微观 理论来对一个复杂系统的演化过程进行完全的动力学描述并不是适合的,因为哪里包含了大 量的自由度。当今大部分流行的模型是输运理论,即要区分集体自由度(相关变量)和内禀自
由度(非相关变量)。后者在平均意义上被处理成一个热浴,集体和内禀自由度之间的能量转 换的量度就是耗散或摩擦。故许多现象能够类比成布朗运动


 反常扩散的描述方法
1.1 连续时间无规行走(CTRW) CTRW模型[7~10]中有两个基本要素:一次跳跃的长 度和两次跳跃之间的等待时间。两者均是随机变量,分 别由跳跃分布密度函数(pdf )φ(x ,t)所定义。那么跳 跃长度的pdf为λ(x )=∫ ∞ 0dt φ(x ,t),等待时间的 pdf为 w(t )=∫ ∞ -∞dx φ(x ,t)。其物理意义分别是:λ(x)dx表 示在间隔(x ,x+dx)跳跃长度的几率;w(t )dt表示在间 隔(t ,t+dt)一次等待时间的几率。如果跳跃长度和等 待时间是独立的,那么 φ(x ,t )=w(t )λ(x)。不同的
603 物 理 学 进 展     25卷
CTRW过程是以平均等待时间T=∫ ∞ 0dtw(t )t和跳跃长度的方均∑2=∫ ∞ -∞dx λ(x )x 2是否收敛 和发散来分类的。对于标准的布朗运动,两者都是确定的;而对于反常扩散,两者之一必发散。 令W(x ,t)代表发现无规行走者在t时刻处于x位置的分布密度函数,它满足一个广义主 方程 W(x ,t )=∫ ∞ -∞dx′∫ ∞ 0 dt′W(x′ ,t′ )φ(x-x′ ,t-t′ )+Ψ(t )δ (x ) (1) 式中Ψ(t )=1-∫ t 0dt′w(t′)为(0,t)时间段内粒子等待不跳的几率。 对方程(1)进行空间Fourier变换和时间 Laplace变换,则 W(x ,t)的 FourierLaplace变换 满足以下的代数方程 W(k ,u)=1-w(u) u W0(k ) 1-φ(k ,u) (2) 这里W0(k)是初始条件W0(x)的Fourier变换。 让我们考虑一个泊松等待时间分布w(t )=τ-1exp(-t/ τ)[T=τ]和高斯跳跃长度分布λ (x )=(4πσ 2)-1/ 2exp(-x 2/ (4σ 2))[∑2=2σ 2]。对于初始条件 W0(x)=δ(x),传播子的 FourierLaplace变换为W(k ,u)=(u+K 1k 2)-1,这里 K 1=σ 2/ τ。变回到(x ,t)空间,得到著名 的高斯传播子解:W(x ,t )=(4πK 1t )-1/ 2exp(-x 2/ (4K 1t ))。事实上,只要T和∑2同时是有限 的,就可以给出相同的结果。 接下来我们考虑平均等待T时间发散,而方均跳跃长度∑2有限的情况,即一个长尾等待 时间分布


接下来我们考虑平均等待T时间发散,而方均跳跃长度∑2有限的情况,即一个长尾等待 时间分布。其分布密度函数的渐近表达式为w(t )≈A α(τ/t )1+α(0<α<1),它的Laplace变换 是w(u)≈1-(uτ) α。那么跳跃密度函数W(x ,t)的FourierLaplace变换为 W(k ,u)= W0(k )/u 1+K αu -αk 2 (3) 使用分数积分的Laplace变换规则:L{ 0D-P t W(x ,t)}=u -PW(x ,u),p ≥0,我们有分数扩散方程  W  t=0D 1-α t K α  2  x 2W(x ,t ) (4) 式中,RiemannLiouville分数导数(积分)被定义成[11] 0D 1-α t W(x ,t )= 1 Γ(α)   t ∫ t 0dt′W(x ,t′ ) (t-t′ ) 1-α (5) 进而给出自由粒子方均位移的渐进表示式 〈x 2(t)〉= 2K αΓ (1+α)t α, 0<α<1 (6) 这里广义扩散系数K α≡σ 2/ τ α。 让我们来考查分数导数和积分的算例[6]与意义。常数 C的 1/2阶导数是0D 1/ 2 x C=C/ π槡 x,而1/2阶积分为0D-1/ 2 x C=2C x/ 槡 π;x的1/2阶导数是0D 1/ 2 x x=2 x/ 槡 π,而1/2阶积分 为0D-1/ 2 x =4x 3/ 2/ 3槡π。分数速度定义成 ν=d μ/ dt μx(t),则粒子的轨道写作 x(t )=Γ-1(μ) ∫ t 0dt′ ν (t′)(t-t′ ) μ-1。扩散粒子的微观运动是完成了一条双绞线和处处非可微的曲线,对于 这个曲线,粒子的轨道是速度的一个功率权重的平均。事实上,它是由粒子的速度记忆效应引 起的,对于分数路径,瞬间速度和位移并不贡献到粒子的宏观运动




接下来我们考虑平均等待T时间发散,而方均跳跃长度∑2有限的情况,即一个长尾等待 时间分布。其分布密度函数的渐近表达式为w(t )≈A α(τ/t )1+α(0<α<1),它的Laplace变换 是w(u)≈1-(uτ) α。那么跳跃密度函数W(x ,t)的FourierLaplace变换为 W(k ,u)= W0(k )/u 1+K αu -αk 2 (3) 使用分数积分的Laplace变换规则:L{ 0D-P t W(x ,t)}=u -PW(x ,u),p ≥0,我们有分数扩散方程  W  t=0D 1-α t K α  2  x 2W(x ,t ) (4) 式中,RiemannLiouville分数导数(积分)被定义成[11] 0D 1-α t W(x ,t )= 1 Γ(α)   t ∫ t 0dt′W(x ,t′ ) (t-t′ ) 1-α (5) 进而给出自由粒子方均位移的渐进表示式 〈x 2(t)〉= 2K αΓ (1+α)t α, 0<α<1 (6) 这里广义扩散系数K α≡σ 2/ τ α。 让我们来考查分数导数和积分的算例[6]与意义。常数 C的 1/2阶导数是0D 1/ 2 x C=C/ π槡 x,而1/2阶积分为0D-1/ 2 x C=2C x/ 槡 π;x的1/2阶导数是0D 1/ 2 x x=2 x/ 槡 π,而1/2阶积分 为0D-1/ 2 x =4x 3/ 2/ 3槡π。分数速度定义成 ν=d μ/ dt μx(t),则粒子的轨道写作 x(t )=Γ-1(μ) ∫ t 0dt′ ν (t′)(t-t′ ) μ-1。扩散粒子的微观运动是完成了一条双绞线和处处非可微的曲线,对于 这个曲线,粒子的轨道是速度的一个功率权重的平均。事实上,它是由粒子的速度记忆效应引 起的,对于分数路径,瞬间速度和位移并不贡献到粒子的宏观运动。
613 4 期 包景东:分数布朗运动和反常扩散
一般阻尼和有势情况下,分数FokkerPlanck方程表为  W(x ,ν ,t )  t +ν W  x+f (x )m  W(x ,ν ,t )  v =γ α0D 1-α t L ^ FPW(x ,ν ,t ) (7) FP算子是L ^ FP= νv+k BT m  2 ν 2。这里分数FokkerPlanck方程仅能描述欠扩散(0<α<1),物理 上适合于长等待系统。但并不依赖于等待时间和跳跃分布的密度函数,它的推导以及线性和 简谐势场的解见文献[9,10]


2 Levy飞行 与上述分数扩散给出的分数 FokkerPlanck方程的情形相反,我们考虑一个长跳跃过程, 即Levy飞行。它具有有限的等待时间但跳跃长度分布的二次距发散,后者一种可能的Fourier 变换是λ(k )=exp(-σ μ| k| μ)≈1-σ μ| k| μ,这里1<μ<2,其对应的跳跃分布的渐进形式为λ (x )≈A μσ-μ| x|-1-μ。自由粒子的扩散方程写作 t W=K μ -∞D μ xW(x ,t),其解是一个Fox函数[7] W(x ,t )= 1 μ|x|H 1, 1 2, 2 |x| (K μt ) 1/ μ (1, 1/ μ),(1, 1/ 2) (1, 1),(1, 1/ 2[ ] ) (10) 自由粒子(甚至简谐势中的粒子)方均位移发散,〈x 2(t)〉 →∞。但我们应该关注的是Levy飞行分布的宽度而不是 二次距,那么前者定义为〈x 2(t)〉 L=∶∫ L 1/ μt -L 1/ μt dxx 2W(x ,t ) ≈t 2/ μ,是一种超扩散[13,14]。这种运动还可以用满足 Levy 统计的噪声 η(t)驱动的 Langevin方程 x · (t )=F(x)/ (γm)+η(t)来描述,其中 η(t)分布密度函数的 Fourier 变换是 p (k )=∫ dηexp(-ik η)p (η)=exp(-D|k| μ) (11) 式中0<μ<2。方程(11)所对应的描写超扩散的分数 FokkerPlanck方程为
623 物 理 学 进 展     25卷
  t W(x ,t )=-  x
F(x )W(x ,t ) γ( ) m +DμW(x ,t ) (12) 这里Riesz分数算符μ由它的Fourier变换来确定,即μ=-(2π)-1∫dkexp(ik



第25卷 第4期 物 理 学 进 展 Vol.25,No.4  2005年12月 PROGRESSINPHYSICS Dec.,2005 文章编号:10000542(2005)0403599
收稿日期:20051028 基金项目:国家自然科学基金资助项目(10235020,10075007);教育部跨世纪优秀人才基金项目;高等学校博 士点专项科研基金(20050027001)
分数布朗运动和反常扩散
包景东 (北京师范大学物理系,北京 100875)
摘 要:  本文评述了分数布朗运动和反常扩散现象及描写它们的几种数学方式。报告了我们 在弹道扩散的产生条件、起源和长时间效应方面的工作。 关键词:分数布朗运动;反常扩散;连续时间无规行走;广义Langevin方程;弹道扩散 中图分类号:O415.6       文献标识码:A 0 引言 历史上,布朗运动现象追源于JanIngenhausz在1785年观测木炭粉末在酒表面的运动;后 来RobertBrown在1828年观测花粉、尘埃、烟灰在水表面的扩散;1905年AlbertEinstein用随 机行走解释流体分子的热运动;1926年 JeanPerrin因为测量 Avogadro常数获 Nobel物理奖。 整整一百年前,伟大的物理学家爱因斯坦在一篇开创性的文章[1],建立了布朗运动的扩散理 论。他考虑一个布朗粒子的一维无规运动,粒子每隔τ时间被撞击一次而移动距离 l,每次撞 击时向左和向右移动的可能性各占一半。假设粒子从原点出发,在时刻 t,粒子已受到了 n=t/ τ次撞击。爱因斯坦证得:粒子的平均位移为零,〈x(t)〉=0;方均位移写作〈x 2(t)〉= 2Dt,这里D=l 2/ (2τ)。其实这是扩散方程 t P(x ,t )=D 2 xP(x ,t)的高斯分布函数解的前两次 距,也是中心极限定理的直接结果。 上世纪40年代以后,布朗运动的动力学研究开始兴起[2, 3]。两个等价的方程是 Langevin 方程和FokkerPlanck方程 [4],前者是关于粒子轨道的随机微分方程;后者是关于分布函数随 时演化的二阶偏微分方程。人们通常认为前者更基本,因为在有些情况下写不出精确的 FokkerPlanck方程,例如色噪声和记忆阻尼。Kubo认为布朗运动可以看成一个系统受到外界 随机扰动的响应,提出了第一和第二涨落耗散定理[5]。站在更广义的角度上,用纯粹的微观 理论来对一个复杂系统的演化过程进行完全的动力学描述并不是适合的,因为哪里包含了大 量的自由度。当今大部分流行的模型是输运理论,即要区分集体自由度(相关变量)和内禀自
由度(非相关变量)。后者在平均意义上被处理成一个热浴,集体和内禀自由度之间的能量转 换的量度就是耗散或摩擦。故许多现象能够类比成布朗运动。 近几年来,在破缺媒介及非大数定理统计下的反常扩散现象引起了人们的极大关注[6~8], 例如在湍流、等离子体、渗透媒介、生长表面和细胞等环境中的系统就表现出偏离布朗运动的 特征。一个自由粒子的方均位移在长时间后正比于时间的分数次幂,称为分数布朗运动: 〈x 2(t)〉∝t δ,其中0<δ<1为欠扩散;δ=1是正常扩散;1<δ<2为超扩散;δ=2系弹道扩散。 从数学上看,反常扩散来自于时间和空间上的非局域性,而描写这一运动的主要手段有:连续 时间无规行走(CTRW),分数FokkerPlanck方程,Levy飞行,Tsallis统计,广义 Langevin方程 等。现在反常扩散和输运的研究才刚刚起步,仅自由场、线性场和简谐势可以获得精确解。由 于数值求解分数FokkerPlanck方程的困难,人们知道的有势系统的信息还很有限,目前的理 论研究大都是从唯象观点出发,导出分布函数或粒子轨道满足的方程,其中反常指数为一个自 由参数。 玻耳兹曼建立了各态历经理论,即长时间后可观测量的系综平均等于时间平均;各态历 经系统具有唯一的定态分布,分布函数在等能量面上为一常数。这也是统计物理的基础,也 就是微正则系统的等概率法则。该理论的目的在于在相空间理解不可逆的起源。其实这个问 题很重要,构成了平衡态统计物理和大部分非平衡态统计物理的基础,但很少有好的例子表明 在孤立系统(自由粒子)各态历经被破坏。近年来,分数布朗运动和非玻耳兹曼统计为这一问 题的研究打开了话题。众所周知,平衡态统计物理侧重于统计分布,非平衡态统计物理侧重于 动力学。我们感兴趣于两者的结合,也就是在什么动力学规律支配下,系统从不同的初始条件 出发,所达到的渐近稳定态有别于平衡态,而是介于牛顿确定性力学和朗之万随机力学之间。 这不仅需考虑耦合方式的非线性,而且更应该深入探讨热浴的结构。
图1 连续时间无规行走(CTRW)模型的 图示。粒子在某一晶格位置上的 等待时间长短用一个圆圈表征,其 中等待时间正比于圆的直径
本文第1节评述了几种描写分数布朗运动和反常扩散现象的数学手段,并与正常的布朗 运动做了比较;第2节报告了我们最近关于弹道扩散的工作,包括产生的条件,可能的物理起 源以及在这一过程中Kubo第一和第二涨落耗散定理的 适用性;小结被最后写在第3节。 1 反常扩散的描述方法
1.1 连续时间无规行走(CTRW) CTRW模型[7~10]中有两个基本要素:一次跳跃的长 度和两次跳跃之间的等待时间。两者均是随机变量,分 别由跳跃分布密度函数(pdf )φ(x ,t)所定义。那么跳 跃长度的pdf为λ(x )=∫ ∞ 0dt φ(x ,t),等待时间的 pdf为 w(t )=∫ ∞ -∞dx φ(x ,t)。其物理意义分别是:λ(x)dx表 示在间隔(x ,x+dx)跳跃长度的几率;w(t )dt表示在间 隔(t ,t+dt)一次等待时间的几率。如果跳跃长度和等 待时间是独立的,那么 φ(x ,t )=w(t )λ(x)。不同的
603 物 理 学 进 展     25卷
CTRW过程是以平均等待时间T=∫ ∞ 0dtw(t )t和跳跃长度的方均∑2=∫ ∞ -∞dx λ(x )x 2是否收敛 和发散来分类的。对于标准的布朗运动,两者都是确定的;而对于反常扩散,两者之一必发散。 令W(x ,t)代表发现无规行走者在t时刻处于x位置的分布密度函数,它满足一个广义主 方程 W(x ,t )=∫ ∞ -∞dx′∫ ∞ 0 dt′W(x′ ,t′ )φ(x-x′ ,t-t′ )+Ψ(t )δ (x ) (1) 式中Ψ(t )=1-∫ t 0dt′w(t′)为(0,t)时间段内粒子等待不跳的几率。 对方程(1)进行空间Fourier变换和时间 Laplace变换,则 W(x ,t)的 FourierLaplace变换 满足以下的代数方程 W(k ,u)=1-w(u) u W0(k ) 1-φ(k ,u) (2) 这里W0(k)是初始条件W0(x)的Fourier变换。 让我们考虑一个泊松等待时间分布w(t )=τ-1exp(-t/ τ)[T=τ]和高斯跳跃长度分布λ (x )=(4πσ 2)-1/ 2exp(-x 2/ (4σ 2))[∑2=2σ 2]。对于初始条件 W0(x)=δ(x),传播子的 FourierLaplace变换为W(k ,u)=(u+K 1k 2)-1,这里 K 1=σ 2/ τ。变回到(x ,t)空间,得到著名 的高斯传播子解:W(x ,t )=(4πK 1t )-1/ 2exp(-x 2/ (4K 1t ))。事实上,只要T和∑2同时是有限 的,就可以给出相同的结果。 接下来我们考虑平均等待T时间发散,而方均跳跃长度∑2有限的情况,即一个长尾等待 时间分布。其分布密度函数的渐近表达式为w(t )≈A α(τ/t )1+α(0<α<1),它的Laplace变换 是w(u)≈1-(uτ) α。那么跳跃密度函数W(x ,t)的FourierLaplace变换为 W(k ,u)= W0(k )/u 1+K αu -αk 2 (3) 使用分数积分的Laplace变换规则:L{ 0D-P t W(x ,t)}=u -PW(x ,u),p ≥0,我们有分数扩散方程  W  t=0D 1-α t K α  2  x 2W(x ,t ) (4) 式中,RiemannLiouville分数导数(积分)被定义成[11] 0D 1-α t W(x ,t )= 1 Γ(α)   t ∫ t 0dt′W(x ,t′ ) (t-t′ ) 1-α (5) 进而给出自由粒子方均位移的渐进表示式 〈x 2(t)〉= 2K αΓ (1+α)t α, 0<α<1 (6) 这里广义扩散系数K α≡σ 2/ τ α。 让我们来考查分数导数和积分的算例[6]与意义。常数 C的 1/2阶导数是0D 1/ 2 x C=C/ π槡 x,而1/2阶积分为0D-1/ 2 x C=2C x/ 槡 π;x的1/2阶导数是0D 1/ 2 x x=2 x/ 槡 π,而1/2阶积分 为0D-1/ 2 x =4x 3/ 2/ 3槡π。分数速度定义成 ν=d μ/ dt μx(t),则粒子的轨道写作 x(t )=Γ-1(μ) ∫ t 0dt′ ν (t′)(t-t′ ) μ-1。扩散粒子的微观运动是完成了一条双绞线和处处非可微的曲线,对于 这个曲线,粒子的轨道是速度的一个功率权重的平均。事实上,它是由粒子的速度记忆效应引 起的,对于分数路径,瞬间速度和位移并不贡献到粒子的宏观运动。 613 4 期 包景东:分数布朗运动和反常扩散
一般阻尼和有势情况下,分数FokkerPlanck方程表为  W(x ,ν ,t )  t +ν W  x+f (x )m  W(x ,ν ,t )  v =γ α0D 1-α t L ^ FPW(x ,ν ,t ) (7) FP算子是L ^ FP= νv+k BT m  2 ν 2。这里分数FokkerPlanck方程仅能描述欠扩散(0<α<1),物理 上适合于长等待系统。但并不依赖于等待时间和跳跃分布的密度函数,它的推导以及线性和 简谐势场的解见文献[9,10]。 另外,从系统分布密度函数随时演化的角度,一些作者[12]引入了非线性媒质中的过阻尼 FokkerPlanck方程,即  W(x ,t )  t =  x[U′ (x )W(x ,t)]+D 2  x 2[W(x ,t)] ν (8) 它对应一个在Ito意义下的态有关的乘性Langevin方程 x · (t )=-U′ (x )+ 2 槡D[W(x ,t)] (ν-1)/ 2η(t ) (9) 这里η(t)是一个均值为零,方差是1的高斯白噪声。 方程(8)已被应用到多孔媒介中的气体渗透问题,薄液体膜中的传播问题,Marshak波的 热传播 以 及 表 面 生 长。其 稳 定 解 为 Wst(x)=[1-(ν-1)βV(x)] 1/ (ν-1) + /Z,这 里 [f ]+=max {f , 0},Z是归一化常数,β=Z ν-1/ (ν D),V(x)=U(x)-U 0,U 0是势的极小值。自 由粒子的方均位移为〈x 2(t)〉∝t 2/ (ν+1),ν<1为超扩散;ν>1为欠扩散;ν=1系正常扩散。
图2 具有相同跳跃步数的布朗运动 (左)与 μ=1.5的 Levy飞 行 (右)的比较
1.2 Levy飞行 与上述分数扩散给出的分数 FokkerPlanck方程的情形相反,我们考虑一个长跳跃过程, 即Levy飞行。它具有有限的等待时间但跳跃长度分布的二次距发散,后者一种可能的Fourier 变换是λ(k )=exp(-σ μ| k| μ)≈1-σ μ| k| μ,这里1<μ<2,其对应的跳跃分布的渐进形式为λ (x )≈A μσ-μ| x|-1-μ。自由粒子的扩散方程写作 t W=K μ -∞D μ xW(x ,t),其解是一个Fox函数[7] W(x ,t )= 1 μ|x|H 1, 1 2, 2 |x| (K μt ) 1/ μ (1, 1/ μ),(1, 1/ 2) (1, 1),(1, 1/ 2[ ] ) (10) 自由粒子(甚至简谐势中的粒子)方均位移发散,〈x 2(t)〉 →∞。但我们应该关注的是Levy飞行分布的宽度而不是 二次距,那么前者定义为〈x 2(t)〉 L=∶∫ L 1/ μt -L 1/ μt dxx 2W(x ,t ) ≈t 2/ μ,是一种超扩散[13,14]。这种运动还可以用满足 Levy 统计的噪声 η(t)驱动的 Langevin方程 x · (t )=F(x)/ (γm)+η(t)来描述,其中 η(t)分布密度函数的 Fourier 变换是 p (k )=∫ dηexp(-ik η)p (η)=exp(-D|k| μ) (11) 式中0<μ<2。方程(11)所对应的描写超扩散的分数 FokkerPlanck方程为
623 物 理 学 进 展     25卷
  t W(x ,t )=-  x
F(x )W(x ,t ) γ( ) m +DμW(x ,t ) (12) 这里Riesz分数算符μ由它的Fourier变换来确定,即μ=-(2π)-1∫dkexp(ik μ)| k| μ[13]。 1.3 Tsallis统计
在标准的统计物理中,可以假设一个泛函使它在一定条件下的极值函数正好是一个满足 中心极限定理的分布密度,这个泛函就是熵。例如一个粒子的 BoltzmannGibbsShannon熵写 作S[q ]=-k B∫dxp(x )ln[σp(x)],这里 σ是一个特征长度。在条件∫dxp(x)=1和〈x 2〉= ∫dxx 2p (x )=σ 2约束下,通过引入乘子,对熵泛函求极值,可得到分布密度函数:p (x )=exp(- βx 2)/Z,且能确定乘子β=1/ (2σ 2)=1/k BT。 Tsallis [15]最近引入了一个如下的广义q熵 S q[p (x ,v)]= 1-∫ p q(x ,v )dx dv q-1 (13) 约束条件是 ∫ p (x ,v )dx dv=1,  ∫ p q(x ,v )E(x ,v )dx dv=U (14) 当q →1,S q恢复到通常的Boltzmann熵。这里p (x ,v)和E(x ,v)分别是粒子分布密度函数和能 量。通过引入Lagrange乘子,对熵泛函求极值,得到粒子在简谐势中的稳定分布 p q(x ,v )≈ 1-(1-q )β 2[ω 2x 2+mv 2( ) ] 1/ (1-q ) (15)   应该指出的是,Tsallis熵引发了关于非广延性统计的广泛研究,特别地,它被应用到折叠 蛋白质的重组、宇宙射线流、湍流、正负电子湮灭以及金融领域。但是,正如我们熟知的那样, 具有相同渐进形式的两个分布并不一定系同一分布。Tsallis分布下的噪声就与 Levy分布下 的噪声不同,并且产生的机制也不一样,前者是由于依赖分布密度函数有关的扩散系数导致 的,而后者的长尾分布实际上是来自于一个无穷的二次距造成的。 1.4 非Ohmic摩擦 系统的运动方程由如下的广义Langevin方程所描述: mx ¨(t )+m∫ t 0γ(t-s )x · (s )ds+U′ (x )=ξ (t ) (16) 式中的无规力ξ (t)服从高斯分布,其系综平均为零,关联满足Kubo第二涨落耗散定理 〈ξ (t )ξ (t′)〉=mk BT γ(t-t′ )=2mk BT∫ ∞ 0 dω πJ (ω)cos ω(t-t′ ) (17) 这里k B是 Boltzmann常数,T是热浴温度,热浴振子的谱密度 J(ω)为非 Ohmic形式[16],即 J (ω)=γ δ(ω/ ω r ) δ-1,其中ω r是一个使mγ δ具有时间倒数量纲的参考频率。 对于不存在势的情况,方程(16)的解能用Laplace变换法获得,粒子速度的平均和关联为 〈v (t)〉=v 0E 2-α[-(ω αt ) 2-α] (18) 和
633 4 期 包景东:分数布朗运动和反常扩散
〈v (t )v (t′)〉=k BT mE 2-α[-(ω α|t-t′| ) 2-α]+ v 2 0-k BT( )m E 2-α[-(ω αt ) 2-α]E 2-α[-(ω αt′ ) 2-α] ( 19) 式中ω2-α α =γ α槇 ω1-αsin-1(απ/ 2),E α(x)是 MittagLeffler函数,特别地,当 α=1,E 1(x)退化为 指数函数e x。从方程(19)可以看出,二次速度关联函数包含了稳定和老化部分。对于等时间 情况,粒子的均方位移写作 〈x 2(t)〉=2k BT mt 2E 2-α, 3[-(ω αt ) 2-α] (20) 在长时间极限下,〈x 2(t)〉≈2k BT mω-2 α (ω αt ) α Γ(1+α)。
图3 100Mo+100Mo反应核系统的熔合几 率随质心能量的变化,其中黑点是 实验结果
1.5 应用 我们现在应用非 Ohmic摩擦模型到一个典型的重 核熔合系统:100Mo+100Mo [17],计算其熔合几率 Pfus作为 质心能量Ecm的函数,结果见图3。势被选成一个倒谐振 子势:U(x )=-1 2mω 2 bx 2,我们用双中心核表面模型和 墙加窗一体耗散机制获得形变势能、惯性和摩擦张量, 核在鞍点的温度定义为aT 2=Ecm-U b,其中能级密度参 数a=A/10(A是复合核核子数),U b为鞍点势能。从图 中可见,正常扩散给出的熔合几率随质心能量变化较 快,而用 δ=0.7的欠扩散所计算的结果与实验值较好 地符合。这是因为复合核大形变状态,其颈部窗口近似 闭合,那么混沌修正核单体耗散使其大为降低,导致能量耗散变弱。我们认为需考虑欠扩散机 制,以符合重核熔合几率随质心能量缓慢变化这一实验结果。 总之,连续时间无轨行走(CTRW)模型和 Tsallis熵与统计模型在理论上极大地扩展了标 准的统计物理,在观念上非常吸引人;但从考虑复杂的有势层面上来看,非 Ohmic摩擦框架下 的广义Langevin方程更容易获得反常输运的数值结果。 2 推广的热布朗运动:弹道扩散 采用Laplace变换法,我们可以给出广义Langevin方程(16)在无势情况的解 x (t )=x 0+v 0H(t )+1 m∫ t 0H(t-t′ )ξ (t′ )dt′ v (t )=v 0H · (t )+1 m∫ t 0H · (t-t′ )ξ (t′ )dt′ (21) 这里x 0和v 0是粒子的初始位置和速度,响应函数H(t)为 ^ H(s )=[s 2+s ^ γ(s)]-1的逆 Laplace 变换,^ γ为阻尼核函数γ(t)的Laplace变换。响应函数导数的一般形式可写为H · (t )=b+∑ jres
643 物 理 学 进 展     25卷
[ H ^ (s)]exp(s j ,t),其中b=(1+lims →0^ γ(s )/s )-1,s j是特征方程s+^ γ(s )=0的非零根。 在长时间极限下,我们有 {〈v 2(t→∞)〉}=k BT m +b 2 {v 2 0}-k BT( )m {〈x 2(t→∞)〉}= bk BT m +b 2 {v 2 0}-k BT( )[ ] m t 2 (22) 式中〈…〉和{…}分别代表对系综和初始态的平均。这表明:粒子方均位移正比于时间的平 方,呈现出弹道扩散,这是超扩散的极限[18]。若粒子的初始分布不是在平衡态,则可引入一个 系统的有效温度Teff=T+b 2(T 0-T),这里T 0=∶{v 2 0}m/k B,Teff=∶〈v 2(t →∞)〉m/k B。这也 意味着系统具有非唯一的定态,也就是说对如此的非马尔科夫布朗动力学,各态历经性被破 坏,Kubo第一涨落耗散定理被破坏。其原因在于 ^ γ(s )=∫ ∞ 0γ(t )dt=0,即系统的有效阻尼消 失。产生弹道扩散的一个等价条件是热色噪声的低频部分被滤掉[18]。 为了探讨弹道扩散的物理起源,我们推广熟知的系统加热库(环境)模型,前者处理为布 朗粒子,两者的耦合提供了布朗粒子所受到的无规力和阻尼。总的Hamiltonian量写作 H=1 2m x 2+U(x )+∑ N j=1 mj 2( q 2 j+ω 2 jq 2 j)+g(x , x ,q j , q j[ ] ) (23) 这里x和q j分别是系统和环境振子的坐标。通常人们仅考虑坐标坐标耦合,即 g=-c j xq j+ c 2 jx 2/ (2mjω 2 j),而认为速度有关的耦合通过变换可以等价于坐标坐标耦合。但我们发现:虽 然绝热积分环境振子的坐标后所得到的关于系统坐标的广义 Langevin方程形式上相同,但两 者的热噪声谱有很大的不同[19,20]。速度有关的耦合可以在量子超导干涉仪,带电粒子在辐射 场运动和磁场中涡旋扩散中发现。例如,如果假设热库振子的频率具有简谐分布(噪声的阻 尼和频率参数分别是Γ和Ω),则系统坐标和热库坐标耦合所产生的热色噪声(简谐噪声)谱 密度函数为 S(ω)= 2γ 0Ω 4k BT (ω 2-Ω 2) 2+Γ2ω 2 (24) 式中γ 0是系统的阻尼强度。而对于系统坐标(速度)和热库速度(坐标)耦合情况,耦合项为g =-d j x q j (-d j xq j ),其所产生的热色噪声(简谐速度噪声[21])谱密度函数是 S(ω)= 2γ 0Γ2ω 2k BT (ω 2-Ω 2) 2+Γ2ω 2 (25) 显然,对(25)式而言,S(ω=0)=0,则 b=(1+γ 0ΓΩ-2)-1;而对(24)式而言 S(0)=2γ 0k BT, 则b=0。 现对无势系统添加上一个常数力F,则长时间后,粒子出现一个耗散加速度 a=F mb,这里 0≤b ≤1。对通常的布朗运动、欠扩散和超扩散情况,b=0,即长时间后外力被摩擦阻力平衡; 但对弹道扩散过程,由于系统的有效阻尼为零,则外力仅被部分平衡。此外,粒子的渐进平均 速度为〈v (t →∞)〉=v 0b,粒子对其的初始状态保留着部分记忆。所以我们称这一结果为介于 牛顿力学和朗之万随机力学之间的动力学[22]。另外,我们还研究了在热简谐速度噪声驱动下 的粒子在无界势和不对称周期势中运动的行为[23]。结果发现,一个闪烁棘轮整流弹道扩散所 653 4 期 包景东:分数布朗运动和反常扩散
给出的定向流远大于整流正常扩散的。
图4 在热简谐速度噪声作用下被广义Langevin方程描述的自由粒子的速度(a)和坐 标(b)的二次距。其中热浴温度T=1. 0 3 结语 在爱因斯坦100年前创立的无规行走扩散理论中,行走是一个不连续过程,即假设在每个 节拍振荡时间内,粒子随机跳跃一次。而连续时间无规行走引入了两次成功跳跃之间有一个 等待时间分布,且等待时间分布的一次距和跳跃长度分布的二次距中有一个是发散的。其揭 示欠扩散对应于长的等待时间,由分数 FokkerPlanck方程来描写;而超扩散则意味着粒子存 在长跳跃的物理机制,Levy飞行就属于此类,这种运动也可以通过长尾分布的非高斯噪声激 励来产生。布朗运动的高斯解是中心极限定理的必然结果。便于数值模拟有势中的反常输运 的当属含非Ohmic记忆摩擦的广义 Langevin方程。系统长时间动力学行为取决于噪声功率 谱的低频分布,因此当热色噪声的低频被滤掉,导致系统的有效摩擦消失,自由粒子呈现热扩 散的极限情形:弹道扩散。在这一过程中,Kubo第二涨落耗散定理被满足,但系统的稳定态部 分依赖于它的初始速度分布,各态历经性和平衡态被破坏。 作者感谢与德国Augsburg大学的 PeterHanggi教授,中国原子能科学研究院的卓益忠研 究员和中国科学院力学研究所的赵江林博士进行过的有益讨论和交流。
参考文献 [1]EinsteinA.AnnPhys(Leipzig),1905,17:549. [2]KramersHA.Physica ,1940,7:284. [3]ChandrasekharS.RevModPhys ,1943,15: 1;WangMC,UhlenbeckGE.RevModPhys ,1945,17:323. [4]RiskenH.TheFokkerplanckEquation(Springer ,Berlin,1989). [5]TodaM,KuboR,SaitoN.StatisticalPhysicsⅡ(SpringerVerlag ,SecondCorrectedPrinting ,1995). [6]KlafterJ ,ShlesingerMF,ZumofenG.Phys.Today ,1996,February :33;SokolovI ,KlafterJ ,BlumenA. Phys.Today ,2002,November :48. [7]MetzlerR,KlafterJ.PhysRep,2000,339:1. [8]MetzlerR,KlafterJ.JPhys ,2004,A37:R161.
663 物 理 学 进 展     25卷
[9]BarkaiE.ChemPhys ,2002,284:13. [10]BarkaiE.PhysRev ,2001,E63:046118;BarkaiE,SilbeyRJ.JPhysChem,2000,B104:3866. [11]PodlubnyI.FractionalDifferentialEquations(AcademicPress ,1999). [12]LenziEK,AnteneodoC,BorlandL.Phys.Rev ,2001,E63:051109;BorlandL.PhysRev ,1998,E57: 6634. [13]JespersenS,MetzlerR,FogedbyHC.PhysRev ,1999,E59:2736. [14]BaoJD,WangHY,JiaY,ZhuoYZ.PhysRev ,2005,E72:051105. [15]TsallisT,MendesRS,PlastinoAR.Physica ,1998,A261:534. [16]WeissU.QuantumDissipativeSystems(WorldScientific ,Singapore ,1999). [17]BaoJD,ZhuoYZ.PhysRev ,2003,C67:064606;ZhaoJL,BaoJD.Physica ,2005,A356:517. [18]BaoJD,ZhuoYZ.PhysRevLett ,2003,91:138104. [19]BaoJD,ZhuoYZ.PhysRev ,2005,E71:010102(R). [20]BaiZW,BaoJD,SongYL.PhysRev ,2005,E72:06xxxx. [21]BaoJD,SongYL,JiQ,ZhuoYZ.PhysRev2005,E72:011113. [22]BaoJD,HanggiP,ZhuoYZ.PhysRev ,2005,E72;06xxxx. [23]BaoJD.PhysRev ,2004,E69:016124. FRACTIONALBROWNIANMOTION ANDANOMALOUSDIFFUSION BAOJingdong (DepartmentofPhysics ,BeijingNormalUniversity ,Beijing100875,China) Abstract : ThephenomenaoffractionalBrownianmotionandanomalousdiffusionandtheir descriptivewaysinthemathematicsarereviewedanddiscussed.Ourworksforthecondition,origin anddynamicaleffectsofballisticdiffusionarereported. Keywords :  fractionalbrownianmotion;anomalousdiffusion;generalizedlangevin equation;nonohmicfriction;ballisticdiffusion

No comments:

Post a Comment