带电测试粒子在磁化史瓦西黑洞中的混沌运动.pdf
《带电测试粒子在磁化史瓦西黑洞中的混沌运动.pdf》由会员分享,可在线阅读,更多相关《带电测试粒子在磁化史瓦西黑洞中的混沌运动.pdf(13页珍藏版)》请在咨信网上搜索。
1、当史瓦西黑洞周围存在渐近均匀的外部磁场时,描述带电粒子在史瓦西黑洞附近运动的哈密顿系统会变为不可积系统.类似于这样的相对论哈密顿系统不存在有显式分析解的2 部分分离形式,给显式辛算法的构建和应用带来困难。近一年以来的系列工作提出将相对论哈密顿系统分解为具有显式分析解的2 个以上分离部分形式,成功解决了许多相对论时空构建显式辛算法的难题.最近的工作回答了哈密顿系统显式可积分离数目对长期数值积分精度有何影响、哪种显式辛算法有最佳长期数值性能这两个问题,指出哈密顿有最小可积分离数目即3 部分分裂解形式并且应用于优化的4 阶分段龙格库塔显式辛算法可取得最好精度.由此选择上述数值积分方法并利用庞加莱截面
2、、最大李雅普诺夫指数和快速李雅普诺夫指标研究在磁化史瓦西黑洞附近运动的带电粒子轨道动力学.结果显示:针对某特定的粒子能量和角动量,较小的外部磁场很难形成混沌轨道;较大的正磁场参数容易使轨道产生混沌,并且随着磁场的增大,轨道的混沌程度也随之加强;粒子能量适当变大也可以加剧混沌程度,但负磁场参数和粒子角动量变大都会减弱混沌.关键词天体力学,黑洞,磁场,混沌,算法中图分类号:P138;文献标识码:A1引言近几年来,引力波 1-2 的数十起成功探测以及M87星系中心超大质量黑洞(M87*)照片 3-4)和银河系中心黑洞(SgrA*)照片 5-6 陆续证实了爱因斯坦广义相对论引力理论关于黑洞存在的预言.
3、黑洞物理性质和引力检验已成为研究热点问题史瓦西黑洞是爱因斯坦场方程的一个解.由于史瓦西黑洞存在4 个独立的运动常数,即粒子能量、角动量、粒子的4 速度关系和类似的卡特常数(亦称为方位角运动)7 ,所以是可积系统尽管该系统理论上存在解析解,但分析解很难用初等函数表示,仅仅可以用椭圆积分形式表达而已.大多数可观测2022-04-08收到原稿,2 0 2 2-0 5-3 1收到修改稿*国家自然科学基金项目(119 7 3 0 2 0)资助tM到的黑洞都有等离子体吸积盘,可以在黑洞外产生磁场.这样的磁场一般比较弱小,虽然对黑洞时空没有什么影响,但对荷质比(电荷与质量之比)较大的粒子产生洛伦兹力作用对粒
4、子运动影响不可忽略 8-12 .当史瓦西时空有外部电磁场浸入时,哈密顿雅可比方程很可能会变为不可分离变量的情形,即不存在类似的卡特常数,从而相应的哈密顿系统变为不可积系统.带电粒子在磁化黑洞周围的运动于一定条件下会产生混沌现象 13-2 0 判定试验粒子轨道的混沌性方法有庞加莱截面、最大李雅普诺夫指数 2 1 和快速李雅普诺夫指标 2 2 等.数值方法是研究混沌系统最常用的方法.对于39-1(2)64卷哈密顿系统长期数值积分方法来说,保持哈密顿相流的辛积分器 2 3 自然成为优先选择,可使能量积分误差无长期增长趋势辛算法分为显式辛算法 2 3-2 4 和隐式辛算法 2 5;显隐混合算法 2 6
5、-3 1 也常被使用显式算法比隐式方法具有计算效率优势.相对论时空对应的哈密顿一般没有显式分析解的两部分分离形式,导致显式辛算法的构建和应用有极大困难求解这些相对论时空问题,如采用辛算法的话,过去只能是隐式 2 5 或者显隐混合 2 6-3 1 的形式.当然,结合中点置换的扩大相空间显式类辛方法 14,3 2-3 7 可以适合用于求解不能变量分离的相对论时空问题.近一年以来伍教授团队在天体物理杂志等期刊上发表了系列工作 17,2 0,3 8-4 4 ,,提出将相对论时空对应的哈密顿或时间变换的哈密顿分解为具有显式分析解的2 个以上分离部分形式,成功解决了许多相对论时空构建显式辛算法的难题.这一
6、思路还可以应用于Yoshida高阶显式辛算法 4 5-4 6 问题是哈密顿系统显式可积分离数目是否影响长期数值积分精度还有一个问题是何种显式辛算法有最佳长期数值性能作者最近发表在天体物理杂志上的一个工作 4 4 回答了这两个问题,指出哈密顿有最小可积分离数目即3 部分分裂解形式并且应用于优化的4 阶分段龙格库塔显式辛算法(Partitioned-Runge-Kutta,PRK64)47可取得最好精度.注意这种4 阶PRK方法与同阶Yoshida算法相比,包含更多额外的时间系数和子哈密顿解的组合,故优化后的PRK算法在计算上比Yoshida算法更耗时但额外增加的时间不太多,并且拥有更小的截断误差
7、,因而这种方法值得推荐应用.本文主要目的是沿着近期工作 4 4 的思路,采用4 阶优化的PRK显式辛方法应用于具有显式分析解的3 部分哈密顿分解形式并结合庞加莱截面、最大李雅普诺夫指数和快速李雅普诺夫指标研究在磁化史瓦西黑洞附近运动的带电粒子轨道动力学.2牛物理模型本节先给出描述带电粒子在附带渐近均匀磁场的史瓦西黑洞周围运动的动力学模型,再给出适合于求解磁化史瓦西时空哈密顿系统的辛算法.天文学报2.1史瓦西时空在球坐标(t,r,)中,t表示坐标时间,r、表示3 维空间,类似于空间直角坐标系中的c、y、么.描述一个带电荷q的试验粒子在磁化史瓦西黑洞附近运动的动力学哈密顿函数方程为:T=zgl(p
8、-qA)(pu-qAv),H2其中,q是粒子的电荷,p是由一组标准哈密顿正则方程多一兴确定的产义动量即力表示协变广义动量,A、A,表示电磁四向量势,上下标、v表示不同坐标符号,代表球坐标(t,r,)中任意坐标符号,史瓦西黑洞度规的非零逆变分量gu为:gt=1/gt=-(grr=1/grr=1这里gt、9 r、9 o o、9 s s 是史瓦西黑洞度规的协变分量.4速度iu是坐标u关于原时的导数.两个常数协变动量分量为:pt=gtut=-(1P=gosp+qAp=r?sin262Pt、P分别表示在t与方向上的广义动量,B表示磁场强度,E是粒子的恒定能量,L是粒子的恒定角动量.(4)式中电磁场势的唯
9、一非零协变分量采用 4 9 的形式:B将(3)-(5)式代入(1)式,磁化史瓦西哈密顿可以简化为2 自由度4 维相空间系统,即:2HP-22P1(L_Brsin?e)2个22r2 sin?04期(1)212t=-E,qB2BsinH2E22(3)L.(4)(5)(6)39-264卷Pr、p e 分别表示在r与方向上的广义动量,这里=qB.(6)式实际是经过无量纲化处理后所得到的.光速c和引力常数G单位化处理:C=G=1,无量纲处理就是以黑洞质量M和粒子质量m来进行标度变换:ttM,r rM,B B/M,E mE,pr mpr,L mML,pe mMpe,q mq,H mH.其中m为试验粒子质量
10、,M引力源质量变为几何化单位,即:M=1.哈密顿函数(6)式除含两个常数E和L外,还有哈密顿本身总是一个守恒量.对于类时时空,这个常数是:H2当史瓦西黑洞外无渐近均匀磁场时,哈密顿函数(6)式的第4 个运动积分存在,故而是严格可积的;但当存在渐近均匀磁场时,第4 个运动积分不存在,从而导致系统不可积,进而可能会产生混沌.2.2算法实现文献 3 9 将哈密顿函数(6)式分解为4 个显式可积求解部分构建显式辛算法.作者近期的工作 4 4 指出哈密顿函数(6)式还可分解为3、5个显式可积求解部分来构建显式辛算法,发现3 部分分解方法精度最高.按照这一思路,3 个部分裂解方法如下:H=Hi+H2+H3
11、,其中每个部分的哈密顿为:2H1=1一2(12r2 sin?01P2H2721H很明显,3 个部分哈密顿函数H1、H 2、H 3 都是可积且可以写出解析解,其解析解是时间t的显式函数.求解3 部分的算子依次定义为礼H设h为积分步长,对哈密顿H提供2 个一阶辛算子:S1=HH HR,周娜英等:带电测试粒子在磁化史瓦西黑洞中的混沌运动hS2H=S1H(2)S1H*(2)3个2 阶对称乘积可以构造4 阶Yoshida显式辛积分器 4 5;S4H=S2H(cih)S2H(czh)S2H(cih),(15)1(7)(8)E2+2L gBr2 sin24期S1H=HH HN H.(13)其中S1H*是S1
12、H的共轭辛算子.于是,这2 个一阶辛算子的积是一个对称组合,也就是二阶显式辛方法:(14)式中:C1=1/(22 1/3),C2=12 C1.进一步可以建立优化的4 阶分段龙格库塔PRK64显式辛算法 4 4,4 7 ;PRKg4=S1Ha12 S1Ha1 S1Ha1o S1Ha。S1Has S1Ha,S1Hag S1Ha,S1Haa S1Ha,SIHa S1Ha1:从第(n1)步到第n步,具体离散差分方案如下:S1He2=sin on-/Hs-pfa cos On-/po.m-1,(9)fi=arctan 2(e2,e1),f2=tan(en-1-fi),gH2=fi+arctan(ei+e
13、2)hpo,n-1+fal,(10)(11)(12)hrH2=1/(e1 cos On-1+e2 sin On-1).Hl:=pa+L-(ra)sin2 Ha/2)2/(r-H2)sin2 gHla-L-(rH2)?sin2 g2/2/rH2 _E2/rH2 2)2),(16)39-364卷S1H*p:H=po,n-1+hl-cos On-1(L-(r-1Posin 0n-1/2)2)/(rm-1 sin 0n-1)-B cos n-(L-r-1sin On-1/2)/sin On-1.天文学报sin2 gH2/21)/(rH2)3 sin2 gHa-(cos 0H2 L-(rH2)2sin?
14、gH2/2)/sin gH2.sin 0n-1/2)/(r-1 sin?On-1)-(L-rn-1 sin?n-1/2)/n-1-E2/(rn-1-2),4期f#=arctan 2(e,e),f=tan(0n-1-f),a=f+aretan(e1?+e2)hpi+f2,r*Ha=1/(ei cos On-1+eg sin On-1).式中计算出的1,2,12 列举在表1中.所谓优化算法意味着系数ai、b,之间的自由系数在5阶(优化的4 阶显式辛算法)截断误差项中系数平方和最小化与非优化算法相比,优化算法可以减小截断误差(12)-(16)式都是求解哈密顿函数(6)式或(8)式的显式辛算法.表1P
15、RK64算法相关系数Table 1 Correlation coefficient of PRKs4 integratorsIntegratorsCoefficient3类数值分析先评估算法优劣,再挑选好的算法并结合混沌指标来研究带电粒子轨道的动力学,特别探讨3 个动力学参数对轨道混沌的影响.3.1算算法对比取时间步长为h=1以及参数为磁场参数=8.910-4,粒子能量E=0.995和角动量L=4.6.初始条件为pr=0、=/2 一旦初始半径r给定,pe初值取正值并且由(7)式确定图1为两种4 阶算法的哈密顿误差图图1(a)为两种4 阶算法积分初始半径r=11的轨道所得哈密顿误差H=H+1/2
16、图,其中红色的曲线表示Yoshida4阶显式辛算法S4,绿色曲线表示PRK4阶显式PRK641=12=0.0792036964311962=11=0.130311410182166Q3=Q10=0.2228614958676084=g=-0.3667132680474265=8=0.324648188689706Q6=7=0.109688477876750辛算法PRK64可以看到两条曲线都不随时间的增大而增大,在长期积分中保持平稳显式辛算法S4的哈密顿误差阶数值最后稳定在10 的-9 阶到-8 阶次之间,而优化后的显式辛算法PRK64算法的哈密顿误差阶数值稳定在10 的-13 阶到-12 阶次
17、之间,比S4算法精度高4 个数量级.图1(b)为初始半径r=88的误差图两种4 阶算法的精度非常高大轨道初始半径比小轨道初始半径精度高的原因在于前者平均轨道周期要小于后者,但积分到=10 4 后两条曲线全部上移,原因是截断误差远小于计算机的舍入误差导致的.即使误差有长期增长趋势,但精度仍可达到10 的-13 阶次左右.还要注意到优化后的显式辛算法需要使用更多的计算时间.表2 中列出两种算法的中央处理器(Central39-464卷ProcessingUnit,CPU)时间,表明优化后的PRK方法计算时间比S4算法稍长一些,但总的来说,计算时间依旧很短,最长CPU时间不超过3 min.因此,考虑
18、到精度和效率,PRK64算法值得在实算中使用.-6-7-8-9-1050-12-13-14-15-160-10(b)r=88-11-12HV-13014-15-1650图1两种4 阶算法的哈密顿误差图。磁场参数=8.910-4,粒子能量E=0.995,角动量L=4.6.Fig.1 Hamiltonian errors for two fourth-order integrators.Magnetic parameter =8.910-4,particle energyE=0.995,angular momentum L=4.6.表2 图1中两种算法的CPU时间Table 2 CPU times
19、 for the two integrators inFig.1IntegratorsS4PRK64r=11122r=881483.2轨道动力学图1中的两个轨道对同一算法的哈密顿误差的长期表现不同是由于两轨道有不同的动力学性质.实际上,初始半径r=11对应有序轨道,即轨道是拟周期的;而初始半径r=88对应混沌轨道,即轨周娜英等:带电测试粒子在磁化史瓦西黑洞中的混沌运动(a)r=11S4PRK.412345671gTPRKWMV124期道演化对初始条件微小改变有指数式的敏感依赖性.轨道的有序和混沌动力学属性可借助轨道的3 维空间构型和平面投影图来粗略地观察选取初始轨道半径r=25,粒子能量E=0
20、.995和角动量L=4图2(a)-(c)画出了3 个不同磁场参数的3 维空间轨道构型,图中红色为轨道在赤道面0=元/2 上的投影在图2(a)中,=0显示试验粒子在史瓦西黑洞时空中的运动轨道限于赤道平面上,并且围绕黑洞做圆轨道周期运动在图2(b)中,=1.210-4对应3 维空间轨道不在一个平面上,其投影是一组相交环曲线在图2(c)中,=0.01对应3 维空间轨道也不在一个平面上,其投影由一些互相交错环曲线组成.图2(c)与图2 (b)相比,轨迹更杂乱、更没有规则.这两个运动轨道是混沌还是有序难以判断.轨道有序或混沌的准确判定还需要借助其他方法.3.2.1混沌指标判断带电粒子在史瓦西黑洞周围运动
21、的轨道是否混沌有许多方法.例如庞加莱截面法、李雅普345671gt1532*24诺夫指数和快速李雅普诺夫指标等.庞加莱截面适合研究2 自由度的4 维保守系统,可用于判断系统的运动状态.若以0=元/2 为曲面,上方为0 元/2,下方为0 元/2,可以按照线性插值来求曲面上的点.当曲面上只得到几个点或者这些点构成一封闭曲线时,系统做拟周期运动;当截面上得到的点是杂乱无章随机分布在一个区域时,系统做混沌运动.李雅普诺夫指数(Lyapunov indicator)是衡量两邻近轨道随时间平均指数分离比的指标,能够反映轨道混沌的强度.最大李雅普诺夫指数判断混沌较为准确,计算方法主要有变分法和两粒子法 2
22、1.本文使用变分法计算:(17)T80 T13(0)1式中:S(T)和S(0)分别表示T时刻与初始时刻的切向量.入值趋于稳定到一个正值,说明有界轨道处在混沌状态;入值趋于0,说明有序轨道处于有序状态.39-564卷入是定义的一个值,用来表示轨道的混沌状态.快速李雅普诺夫指标(Fast Lyapunov indica-tor,FLI)与李雅普诺夫指数相比,可以更快地看到轨道的混沌,得到混沌指数,其定义为 50 :FLI=lg S()I.若FLI随时间增长很快,则有界轨道混沌;若FLI随时间增长非常慢,则有界轨道有序.(a)=00=元/2(b)=1.210-49=元/2(c)=0.019=元/2图
23、2 3 维轨道构型及其在赤道面0=元/2 上的投影Fig.2 Three-dimensional trajectories and their projectionsat the plane=/23.2.2粒子运动轨道的动力学图3 为赤道面=元/2 上的庞加莱截面.图3(a)为图2(b)中粒子轨道在赤道面=/2上的庞加莱截面图,截面上的轨道是闭合的环轨道,说明图2(b)中粒子做有序的拟周期运动.图3(b)为图2(c)中粒子轨迹的庞加莱截面图,截面上的轨道是随机离散的点分布在一个面区域上,表明粒子做混天文学报沌运动显然,相比3 维空间构型和平面投影图,庞加莱截面区分轨道有序或混沌要更清楚、准确带
24、电试验粒子在史瓦西黑洞周围运动时,磁场破坏了第4 个运动积分的存在,导致系统不可积,是粒子运动混沌的根本原因.(18)0.20.10.0-0.1-0.200.40.30.20.10.0-0.1-0.2-0.3-0.4图3 粒子轨道在赤道面=元/2 上的庞加莱截面.图(a):图2(b)中粒子的轨道.图(b):图2(c)中粒子的轨道.Fig.3 Poincar cross-section of the particle orbit at theplane =/2.Panel(a):the particle trajectory is the sameas that in Fig.2(b).Pane
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 带电 测试 粒子 磁化 史瓦西 黑洞 中的 混沌 运动
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。