分数布朗运动和反常扩散 - 物理学进展
一个布朗粒子的一维无规运动,粒子每隔τ时间被撞击一次而移动距离l,每次撞
击时向左和向右移动的可能性各占一半。假设粒子从原点出发,在时刻t,粒子已受到了
n=t/τ次撞击。爱因斯坦证得:粒子的平均位移为零,〈x(t)〉=0;方均位移写作〈x2(t)〉=
2Dt,这里D=l2/(2τ)。其实这是扩散方程tP(x,t)=D2
xP(x,t)的高斯分布函数解的前两次
距,也是中心极限定理的直接结果
第25卷 第4期物 理 学 进 展Vol.25,No.4
2005年12月PROGRESSINPHYSICS Dec.,2005
文章编号:10000542(2005)0403599
收稿日期:20051028
基金项目:国家自然科学基金资助项目(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;方均位移写作〈x2(t)〉=
2Dt,这里D=l2/(2τ)。其实这是扩散方程tP(x,t)=D2
xP(x,t)的高斯分布函数解的前两次
距,也是中心极限定理的直接结果。
上世纪40年代以后,布朗运动的动力学研究开始兴起[2,3]。两个等价的方程是Langevin
方程和FokkerPlanck方程[4],前者是关于粒子轨道的随机微分方程;后者是关于分布函数随
时演化的二阶偏微分方程。人们通常认为前者更基本,因为在有些情况下写不出精确的
FokkerPlanck方程,例如色噪声和记忆阻尼。Kubo认为布朗运动可以看成一个系统受到外界
随机扰动的响应,提出了第一和第二涨落耗散定理[5]。站在更广义的角度上,用纯粹的微观
理论来对一个复杂系统的演化过程进行完全的动力学描述并不是适合的,因为哪里包含了大
量的自由度。当今大部分流行的模型是输运理论,即要区分集体自由度(相关变量)和内禀自
由度(非相关变量)。后者在平均意义上被处理成一个热浴,集体和内禀自由度之间的能量转
换的量度就是耗散或摩擦。故许多现象能够类比成布朗运动。
近几年来,在破缺媒介及非大数定理统计下的反常扩散现象引起了人们的极大关注[6~8],
例如在湍流、等离子体、渗透媒介、生长表面和细胞等环境中的系统就表现出偏离布朗运动的
特征。一个自由粒子的方均位移在长时间后正比于时间的分数次幂,称为分数布朗运动:
〈x2(t)〉∝tδ,其中0<δ<1为欠扩散;δ=1是正常扩散;1<δ<2为超扩散;δ=2系弹道扩散。
从数学上看,反常扩散来自于时间和空间上的非局域性,而描写这一运动的主要手段有:连续
时间无规行走(CTRW),分数FokkerPlanck方程,Levy飞行,Tsallis统计,广义Langevin方程
等。现在反常扩散和输运的研究才刚刚起步,仅自由场、线性场和简谐势可以获得精确解。由
于数值求解分数FokkerPlanck方程的困难,人们知道的有势系统的信息还很有限,目前的理
论研究大都是从唯象观点出发,导出分布函数或粒子轨道满足的方程,其中反常指数为一个自
由参数。
玻耳兹曼建立了各态历经理论,即长时间后可观测量的系综平均等于时间平均;各态历
经系统具有唯一的定态分布,分布函数在等能量面上为一常数。这也是统计物理的基础,也
就是微正则系统的等概率法则。该理论的目的在于在相空间理解不可逆的起源。其实这个问
题很重要,构成了平衡态统计物理和大部分非平衡态统计物理的基础,但很少有好的例子表明
在孤立系统(自由粒子)各态历经被破坏。近年来,分数布朗运动和非玻耳兹曼统计为这一问
题的研究打开了话题。众所周知,平衡态统计物理侧重于统计分布,非平衡态统计物理侧重于
动力学。我们感兴趣于两者的结合,也就是在什么动力学规律支配下,系统从不同的初始条件
出发,所达到的渐近稳定态有别于平衡态,而是介于牛顿确定性力学和朗之万随机力学之间。
这不仅需考虑耦合方式的非线性,而且更应该深入探讨热浴的结构。
图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)。不同的
360 物理学进展 25卷
CTRW过程是以平均等待时间T=∫∞
0dtw(t)t和跳跃长度的方均Σ2=∫∞
-∞dxλ(x)x2是否收敛
和发散来分类的。对于标准的布朗运动,两者都是确定的;而对于反常扩散,两者之一必发散。
令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)的FourierLaplace变换
满足以下的代数方程
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(-x2/(4σ2))[Σ2 =2σ2]。对于初始条件W0(x)=δ(x),传播子的
FourierLaplace变换为W(k,u)=(u+K1k2)-1,这里K1=σ2/τ。变回到(x,t)空间,得到著名
的高斯传播子解:W(x,t)=(4πK1t)-1/2exp(-x2/(4K1t))。事实上,只要T和Σ2同时是有限
的,就可以给出相同的结果。
接下来我们考虑平均等待T时间发散,而方均跳跃长度Σ2有限的情况,即一个长尾等待
时间分布。其分布密度函数的渐近表达式为w(t)≈Aα(τ/t)1+α(0<α<1),它的Laplace变换
是w(u)≈1-(uτ)α。那么跳跃密度函数W(x,t)的FourierLaplace变换为
W(k,u)= W0(k)/u
1+Kαu-αk2 (3)
使用分数积分的Laplace变换规则:L{0D-P
t W(x,t)}=u-PW(x,u),p≥0,我们有分数扩散方程
W
t=0D1-α
t Kα 2
x2W(x,t) (4)
式中,RiemannLiouville分数导数(积分)被定义成[11]
0D1-α
t W(x,t)= 1
Γ(α)
t∫t
0dt′W(x,t′)
(t-t′)1-α (5)
进而给出自由粒子方均位移的渐进表示式
〈x2(t)〉= 2Kα Γ(1+α)tα, 0<α<1 (6)
这里广义扩散系数Kα≡σ2/τα。
让我们来考查分数导数和积分的算例[6]与意义。常数C的1/2阶导数是0D1/2
x C=C/
槡πx,而1/2阶积分为0D-1/2
x C=2C x/ 槡π;x的1/2阶导数是0D1/2
x x=2 x/ 槡π,而1/2阶积分
为0D-1/2
x =4x3/2/3槡π。分数速度定义成ν=dμ/dtμx(t),则粒子的轨道写作x(t)=Γ-1(μ)
∫t
0dt′ν(t′)(t-t′)μ-1。扩散粒子的微观运动是完成了一条双绞线和处处非可微的曲线,对于
这个曲线,粒子的轨道是速度的一个功率权重的平均。事实上,它是由粒子的速度记忆效应引
起的,对于分数路径,瞬间速度和位移并不贡献到粒子的宏观运动。
4期包景东:分数布朗运动和反常扩散361
一般阻尼和有势情况下,分数FokkerPlanck方程表为
W(x,ν,t)
t +νW
x+f(x) m
W(x,ν,t)
v =γα0D1-α
t L^
FPW(x,ν,t) (7)
FP算子是L^
FP=
νv+kBT
m
2
ν2。这里分数FokkerPlanck方程仅能描述欠扩散(0<α<1),物理
上适合于长等待系统。但并不依赖于等待时间和跳跃分布的密度函数,它的推导以及线性和
简谐势场的解见文献[9,10]。
另外,从系统分布密度函数随时演化的角度,一些作者[12]引入了非线性媒质中的过阻尼
FokkerPlanck方程,即
W(x,t)
t =
x[U′(x)W(x,t)]+D 2
x2[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)-U0,U0是势的极小值。自
由粒子的方均位移为〈x2(t)〉∝t2/(ν+1),ν<1为超扩散;ν>1为欠扩散;ν=1系正常扩散。
图2 具有相同跳跃步数的布朗运动
(左)与μ=1.5的Levy飞行
(右)的比较
1.2 Levy飞行
与上述分数扩散给出的分数FokkerPlanck方程的情形相反,我们考虑一个长跳跃过程,
即Levy飞行。它具有有限的等待时间但跳跃长度分布的二次距发散,后者一种可能的Fourier
变换是λ(k)=exp(-σμ|k|μ)≈1-σμ|k|μ,这里1<μ<2,其对应的跳跃分布的渐进形式为λ
(x)≈Aμσ-μ|x|-1-μ。自由粒子的扩散方程写作tW=Kμ
-∞Dμx
W(x,t),其解是一个Fox函数[7]
W(x,t)= 1
μ|x|H1,1
2,2
|x|
(Kμt)1/μ
(1,1/μ),(1,1/2)
[ (1,1),(1,1/2 ] )
(10)
自由粒子(甚至简谐势中的粒子)方均位移发散,〈x2(t)〉
→∞。但我们应该关注的是Levy飞行分布的宽度而不是
二次距,那么前者定义为〈x2(t)〉L=∶∫L-1tL/μ1t/μdxx2W(x,t)
≈t2/μ,是一种超扩散[13,14]。这种运动还可以用满足Levy
统计的噪声η(t)驱动的Langevin方程x·
(t)=F(x)/
(γm)+η(t)来描述,其中η(t)分布密度函数的Fourier
变换是
p(k)=∫dηexp(-ikη)p(η)=exp(-D|k|μ)
(11)
式中0<μ
No comments:
Post a Comment