混沌体系逃逸率的计算研究.pdf

混沌体系逃逸率的计算 赵海军 杜孟利 北京中国科学院理论物理研究所 摘要 本文讨论了混沌体系中逃逸率的计算问题。我们分别用数值模拟和解 析方法计算了Henon．Heiles体系逃逸率随体系能量的变化。在数值模拟时， 我们分别使用了辛算法和非辛算法求解粒子运动的哈密顿方程。我们发现 由数值模拟得到的逃逸率随着模拟时间的增加与解析公式的偏离变大，这 反映了逃逸并不是严格遵守指数规律的。我们还发现虽然单粒子长时间的 轨道依赖算法，但是，由不同算法得到的逃逸率相互之间符合得很好。 1． 引言和模型 H=圭(z+z)+D(x，y)， (1) u(x，y)=三(z2+J，2)+x2y-；y3 x和Y表示粒子的坐标，见．乃表示对应的动量。我们已令粒子的质量 为l。删体系是个著名的混沌模型。研究表明，如果体系能量E<l／12， 粒子运动是规则的：随着体系能量增加，相空间混沌运动的区域逐渐扩 大。 当体系能量接近E=l／6时，相空间九乎全部都是混沌的了[1]。 131 图1给出了u的等值线。我们看到图1中有三个鞍点只0=0，y=1)， 封闭的，运动被限制于封闭线内。能量为1／6的等势线构成三角形 P]BP,。能量高于1／6的等势线不是封闭的，在三个顶角处有开口。能 量E”．=1／6是体系的闽能值。如果粒子能量高于阈能值E。=1／6，就可以 从开口处逃出HH体系的近似三角形区域。对于弹子球体系的研究表明， 混沌体系逃逸过程中势阱中的粒于数应当随时间指数衰减[2]， _Ⅳ0)=N(O)exp(-ta) r21 其中n为逃逸率，其大小依赖于体系能量E。 一 ∥■弋 X 图1．HH体系的等势线。B，置和只为三个鞍点位置 132 2．逃逸率公式 对于粒子的逃逸率与能量的关系，我们可以推导出一个解析公式。假 设体系能量E大于阈能值，我们通过每个鞍点可画一条和逃逸方向垂直的 轴。我们定义由三条不相交的等势线和三条直线所形成的范围为势阱。由 于混沌的影响，给定能量E的运动假定是各态历经的，那么粒子在势阱中相 空间的分 出。分布的实质是这样的： 只hg／(q,p)-雨6和(E肛-H(q地,p))p)) 在空间上是均匀的(二维)，粒子空间位置确定后，动量大小由哈密顿(1) 给出，动量的方向也是均匀分布的。定义AE=E一瓦，0为粒子动量的方 向角，A(,XE)为前面描述的势阱面积，用p(x，Y，口)表示依赖于空间位置和 动量方向的密度分布，则有p(x，少，秒)=—2x—AI一(a—E)。如果势阱中有N个粒子， 单位时间通过Pl点处的开口逃逸的粒子数为 cos(9)，线积分对AlBl线上的经典允许部分 uS,it』dOp(x,y，秒)l卢(x，y)I 进行。体系的三个出1：3具有对称性，因此，总逃逸数是从一个1：3逃逸数的3 倍。N个粒子在势阱中的粒子数随时间的变化可写成(注意P不依赖于坐 标) 警卅咿x扣／2乡瑚王厕 ㈤ =-2万43zlEpN(t) 133 势阱的大小么(缸)没有解析公式。我们用蒙特卡罗方法计算出的结果 如图2中的点所示。我们发现数值计算结果可以用一个二次多项式很好的 表示，即 么(AE)=4+4△E+4△E2(4) 拟合得到 其中心--1．299是AE=0时等势线形成的等边三角形的面积， 4=9．656以及4：一22．6。公式(4)在图2中用实线表示。 A 图2．势阱面积随能量的变化。点是数值计算的