楊秦男 張延惠? 蔡祥吉 蔣國(guó)輝 徐學(xué)友

1)(山東師范大學(xué)物理與電子科學(xué)學(xué)院,濟(jì)南 250014)

2)(山東省科學(xué)院情報(bào)研究所,濟(jì)南 250014)

(2012年10月10日收到;2012年12月21日收到修改稿)

1 引言

混沌及分形理論是近些年來(lái)發(fā)展起來(lái)的新型非線性理論,在很多方面有著重要研究意義[1].混沌簡(jiǎn)單來(lái)說(shuō)就是研究對(duì)象系統(tǒng)無(wú)規(guī)行為中的規(guī)律性[2],而分形是混沌性質(zhì)中存在的規(guī)律性中的一種.混沌分形與氣象學(xué)、地質(zhì)學(xué)以及物理學(xué)的眾多分支學(xué)科有著廣泛的聯(lián)系,在解決非線性問(wèn)題中起著重要作用.二維介觀器件體系中粒子在輸運(yùn)過(guò)程中的諸多現(xiàn)象和行為都存在著混沌及分形性質(zhì),通過(guò)改變器件參數(shù)或環(huán)境(磁場(chǎng)等)來(lái)探究它們對(duì)粒子逃逸曲線混沌性質(zhì)的影響規(guī)律,同時(shí)發(fā)現(xiàn)體系中存在的分形自相似結(jié)構(gòu)[3-6].本文以逃逸曲線定性比較和分形維數(shù)定量計(jì)算兩種手段來(lái)研究器件中混沌吸引子的變化規(guī)律;首次運(yùn)用“眼式分析”和相似比比較的方法來(lái)驗(yàn)證逃逸曲線中存在的分形自相似結(jié)構(gòu).本文的研究對(duì)分析混沌逃逸規(guī)律、分形自相似結(jié)構(gòu)以及器件輸運(yùn)的穩(wěn)定性控制等具有理論和實(shí)際應(yīng)用意義.

2 器件介紹及理論模型

本文采用的RIKENSide-wall Sinai Billard器件(圖1)由日本理化研究所(RIKEN)納米電子材料實(shí)驗(yàn)室的Cooper和Bird在晶圓上制得,并由新南威爾士大學(xué)的Micolich等[7]進(jìn)行了低溫下電子量子干涉效應(yīng)等實(shí)驗(yàn).為敘述簡(jiǎn)便,本文的研究中將該器件簡(jiǎn)稱為RIKEN器件.

RIKEN臺(tái)球器件可理解為如圖2所示的二維器件模型,圖中標(biāo)度的量級(jí)為微米級(jí).器件底部的兩個(gè)開口很小,粒子無(wú)法通過(guò),因此可以忽略其對(duì)粒子逃逸的影響.器件腔內(nèi)部區(qū)域可視為勢(shì)阱,直線和底部圓弧部分為勢(shì)壘[8].粒子在器件中運(yùn)動(dòng),到達(dá)開口位置時(shí)逃逸.當(dāng)粒子能量增大時(shí),粒子會(huì)從勢(shì)壘上勢(shì)能較低的區(qū)域逃逸,從而導(dǎo)致勢(shì)阱面積增大(但勢(shì)阱的形狀基本不變),粒子的逃逸率增加,相同情況下的粒子逃逸時(shí)間也變短[9-11].本文的研究主要基于統(tǒng)計(jì)單個(gè)粒子逃逸前的碰撞次數(shù),為控制變量以及簡(jiǎn)化計(jì)算過(guò)程,認(rèn)為粒子只從開口處逃逸而無(wú)法從勢(shì)壘逃逸,即二維器件腔內(nèi)部為無(wú)限深勢(shì)阱,邊界勢(shì)壘的勢(shì)能無(wú)限大.

粒子在二維器件腔內(nèi)與器件壁的碰撞視為完全彈性碰撞,碰撞前后的運(yùn)動(dòng)軌跡遵從反射定律[12].由于器件的對(duì)稱性,可定義兩側(cè)頂端的任意一個(gè)開口為粒子的入射口.本文以下研究定義器件左側(cè)開口為粒子入射位置,定義粒子到達(dá)器件兩側(cè)頂端任一開口位置時(shí)逃逸器件腔.

圖1 RIKEN器件的電鏡掃描圖

圖2 RIKEN器件的理論模型

3 粒子初始位置及器件參數(shù)對(duì)輸運(yùn)混沌性質(zhì)的影響

粒子從入射到出射的過(guò)程中會(huì)與勢(shì)壘發(fā)生若干次碰撞,在器件參數(shù)(主要是開口寬度、腔長(zhǎng)以及圓弧半徑)確定的情況下,粒子逃逸前的碰撞次數(shù)由入射角度決定.當(dāng)入射角度微小改變時(shí),碰撞次數(shù)可能不會(huì)明顯改變,也可能會(huì)發(fā)生劇烈變化[13].圖3(a)為入射角分別為21°和22°時(shí)粒子逃逸前的碰撞情況比較,入射角的微小改變沒(méi)有引起碰撞軌跡和次數(shù)的明顯變化;圖3(b)為入射角為41°和41.3°時(shí)粒子逃逸前的碰撞情況比較,入射角的微小改變引起了碰撞軌跡和次數(shù)的劇烈變化.這種產(chǎn)生碰撞次數(shù)劇烈變化時(shí)對(duì)應(yīng)的角度稱作奇異點(diǎn),奇異點(diǎn)的集合就是混沌吸引子.產(chǎn)生這種劇烈變化的原因主要是粒子與器件底部弧線區(qū)域碰撞,當(dāng)初始入射角近似相同的兩個(gè)粒子到達(dá)圓弧部分時(shí),它們反射軌跡之間的夾角及距離差異會(huì)增大,因此圓弧對(duì)微小差異起到了放大的作用(以下簡(jiǎn)稱為圓弧的放大作用).隨著反彈次數(shù)的增加這種差距會(huì)越來(lái)越大,最終導(dǎo)致了入射條件近似的粒子逃逸前碰撞次數(shù)的明顯差異.當(dāng)然,器件的開口邊緣位置、拐角附近區(qū)域、直線區(qū)域等也會(huì)產(chǎn)生這種效果,下文將分別予以討論.器件中粒子逃逸的混沌現(xiàn)象由此產(chǎn)生,這個(gè)過(guò)程稱為“混沌逃逸”[3],其逃逸時(shí)間和碰撞次數(shù)等敏感依賴于粒子初始角度的改變.

圖3 (a)入射角21°(實(shí)線)與22°(虛線)比較;(b)入射角41°(實(shí)線)與 41.3°(虛線)比較

圖4為計(jì)算得到的隨粒子入射角連續(xù)微小改變,粒子碰撞次數(shù)的變化曲線(以下簡(jiǎn)稱為粒子的逃逸曲線).圖中的黑色區(qū)域是粒子碰撞次數(shù)連續(xù)劇烈變化所導(dǎo)致,也就是奇異點(diǎn)所在的混沌區(qū)域;白色區(qū)域是碰撞次數(shù)變化不明顯的區(qū)域.由于碰撞次數(shù)在60次以上的入射情況(圖4(a))很少,混沌區(qū)域的碰撞次數(shù)多在60次以內(nèi),故而可以忽略掉60次以上的碰撞次數(shù),只考慮60次以內(nèi)的碰撞(圖4(b)),本文以下研究均采用此方法.

3.1 粒子初始位置

粒子從器件的左側(cè)開口入射,初始位置改變對(duì)應(yīng)的粒子逃逸曲線變化情況如圖5所示.粒子在不同初始位置下入射,其逃逸曲線中的黑色區(qū)域和白色區(qū)域分布情況基本一致,而且黑色區(qū)域所占比例不變.因此粒子初始位置改變對(duì)逃逸曲線混沌性質(zhì)的影響基本可以忽略,因?yàn)樘右萸€混沌的產(chǎn)生主要是由于入射角的改變,而不是粒子的初始位置.在以下各節(jié)討論中,為控制參數(shù)單一變化,我們均默認(rèn)粒子的初始入射位置為器件左側(cè)開口的中間位置.

3.2 開口寬度

圖4 RIKEN器件碰撞次數(shù)隨入射角變化關(guān)系 (a)碰撞次數(shù)＞60次;(b)碰撞次數(shù)＜60次

圖5 粒子初始位置對(duì)逃逸前碰撞次數(shù)的影響,右側(cè)數(shù)據(jù)為粒子到開口頂端的距離

器件開口對(duì)粒子的逃逸具有選擇作用,尤其是在開口邊緣位置這種選擇作用更明顯.圖3(b)中,當(dāng)軌跡相近的兩粒子分別到達(dá)左端開口邊緣位置時(shí),一個(gè)粒子逃逸而另一個(gè)粒子反射回器件腔無(wú)法逃逸.器件開口寬度的增大使粒子的逃逸概率增大,原逃逸過(guò)程中在開口邊緣無(wú)法逃逸的粒子現(xiàn)在會(huì)被開口選擇,逃逸器件腔.開口邊緣對(duì)粒子逃逸的選擇作用明顯減弱,由此逃逸過(guò)程中碰撞次數(shù)的奇變點(diǎn)減少,逃逸曲線中的混沌區(qū)域會(huì)逐漸被白色區(qū)域取代.圖6為開口寬度分別為0.05—0.30μm時(shí)的逃逸曲線比較,可以看到隨開口寬度增加,白色區(qū)域逐漸增寬增多,器件中粒子逃逸曲線的混沌性減弱.

圖6 開口寬度對(duì)逃逸前碰撞次數(shù)的影響,右側(cè)數(shù)據(jù)為對(duì)應(yīng)的開口寬度

3.3 器件腔長(zhǎng)

圖7 為不同腔長(zhǎng)中粒子的逃逸曲線比較.當(dāng)腔長(zhǎng)增加時(shí)粒子逃逸前在腔內(nèi)的碰撞次數(shù)增加,器件勢(shì)壘的放大作用越來(lái)越明顯.當(dāng)粒子到達(dá)器件的開口位置時(shí),放大的軌道差異最終體現(xiàn)在逃逸前的碰撞次數(shù)差異上,器件腔長(zhǎng)的增加導(dǎo)致逃逸曲線混沌性的增強(qiáng).

3.4 圓弧半徑

器件底部圓弧半徑的變化體現(xiàn)在兩個(gè)方面：一是R≤0.35μm時(shí)弧度等于180°的半圓(如圖8(a));二是R＞0.35μm時(shí)弧度小于180°的圓弧(如圖8(b)).

3.4.1 R≤0.35μm

圖9為R≤0.35μm時(shí)器件底部圓弧半徑變化時(shí)的逃逸曲線比較.在R=0時(shí)(此時(shí)器件的底部為一條直線)的粒子逃逸曲線中,入射角在0°和180°附近的碰撞次數(shù)整體較多且變化較大,存在一定的混沌性質(zhì);其他大部分角度碰撞次數(shù)沒(méi)有明顯的奇變,不存在混沌性質(zhì).當(dāng)圓弧半徑增大時(shí),粒子到達(dá)圓弧的概率逐漸增大,圓弧區(qū)域的放大作用越來(lái)越明顯,導(dǎo)致了粒子逃逸曲線中奇變點(diǎn)個(gè)數(shù)增多.因此圖中的黑色混沌區(qū)域越來(lái)越寬,逃逸曲線的混沌性越來(lái)越強(qiáng).

圖7 器件腔長(zhǎng)對(duì)逃逸前碰撞次數(shù)的影響,右側(cè)數(shù)據(jù)為對(duì)應(yīng)的腔長(zhǎng)

圖8 (a)R≤0.35μm時(shí)器件形狀示意圖;(b)R＞0.35μm時(shí)器件形狀示意圖

但是圖9中的混沌區(qū)域在R=0.20μm之后并沒(méi)有明顯的增多,反而在R接近0.35μm時(shí)有了明顯的減少.如圖8(a)所示,拐角①和拐角②對(duì)粒子入射軌道的微小差距均有放大作用.R較小時(shí)兩個(gè)拐角的影響是分別獨(dú)立的;R較大(接近0.35μm)時(shí),將很少有粒子在兩個(gè)拐角處連續(xù)碰撞,兩個(gè)拐角可視為一個(gè)拐角,其放大作用減弱.由此說(shuō)明在該器件中拐角對(duì)混沌性質(zhì)的影響是明顯的：拐角越多越分散,器件的混沌性越強(qiáng);拐角數(shù)量減少或相鄰拐角距離很近時(shí),器件的混沌性有所減弱.

3.4.2 R＞0.35μm

圖10為R＞0.35μm時(shí)器件底部圓弧半徑變化時(shí)的逃逸曲線比較.很明顯,隨著圓弧半徑R的逐漸增大,圖中的黑色區(qū)域逐漸減少,奇異點(diǎn)的個(gè)數(shù)下降,器件的混沌性減弱.

當(dāng)R＞0.35μm時(shí),器件中拐角的數(shù)量及拐角之間的距離不再改變.圓弧的半徑(或者弧度)是影響輸運(yùn)過(guò)程混沌性質(zhì)的主要原因,半徑越大(弧度越小),圓弧部分的放大效果越弱,粒子逃逸過(guò)程中碰撞次數(shù)的奇變點(diǎn)減少,混沌性減弱.可以預(yù)見(jiàn)當(dāng)R→∞時(shí),原來(lái)的弧線將逐漸變成直線,其放大作用等效于直線,其逃逸曲線將與圖9中R=0時(shí)的逃逸曲線相同,因?yàn)樵跇O限情況下兩者所構(gòu)造的器件形狀是完全相同的.

圖9 R≤0.35μm時(shí)圓弧半徑對(duì)逃逸前碰撞次數(shù)的影響,右側(cè)數(shù)據(jù)為對(duì)應(yīng)的圓弧半徑

4 逃逸曲線中分形自相似結(jié)構(gòu)的討論

4.1 逃逸曲線中的分形自相似結(jié)構(gòu)

總結(jié)以上各種情況下的混沌逃逸曲線,奇異點(diǎn)較多的黑色區(qū)域和對(duì)入射條件的微小改變不敏感的白色區(qū)域相互交錯(cuò),形成一類復(fù)雜的曲線,這樣的體系顯示出混沌動(dòng)力學(xué)的標(biāo)志——分形結(jié)構(gòu)[3].當(dāng)然該混沌系統(tǒng)中可能存在“黑帶”自相似的分形結(jié)構(gòu)[14,15],也就是說(shuō)在不同的放大級(jí)別上,逃逸曲線的黑白區(qū)域分布是相似的[16].

回到原始參數(shù)下RIKEN器件中粒子逃逸曲線中來(lái)：由于該逃逸曲線為非對(duì)稱結(jié)構(gòu),無(wú)法從中央?yún)^(qū)域逐級(jí)放大來(lái)尋找其自相似結(jié)構(gòu),所以我們采用將曲線中的混沌區(qū)域逐段放大的方法來(lái)尋找.圖11為混沌區(qū)域存在的分形自相似結(jié)構(gòu),可以看到各級(jí)放大之后的圖形與上一級(jí)圖形在黑白區(qū)域的分布上是一致的,各級(jí)圖形之間有著很好的相似度.

圖10 R＞0.35μm時(shí)圓弧半徑對(duì)逃逸前碰撞次數(shù)的影響,右側(cè)數(shù)據(jù)為對(duì)應(yīng)的圓弧半徑及圓心角角度

圖11 原RIKEN器件中混沌區(qū)域的分形自相似結(jié)構(gòu) 自上而下各級(jí)圖形的取值范圍分別為(49.8°,62.2°),(51.9°,54.6°),(52.35°,52.95°),(52.449°,52.579°),(52.4694°,52.4977°),(52.47395°,52.48003°);右側(cè)數(shù)據(jù)為放大級(jí)數(shù)

4.2 自相似區(qū)域的討論

我們把自相似區(qū)域的位置用圓(或者橢圓)逐級(jí)別標(biāo)示出來(lái),以探討各放大級(jí)的位置與分形自相似性的關(guān)系.圓的直徑(或者橢圓的長(zhǎng)軸)Di為圖11自相似結(jié)構(gòu)中對(duì)應(yīng)的第i條曲線的入射角范圍,得到如圖12(a)所示的嵌套圓結(jié)構(gòu),我們將之稱為“眼式結(jié)構(gòu)”,該結(jié)構(gòu)也是自相似結(jié)構(gòu),放大后的局部與整體是相似的.圖中的橫軸為各級(jí)自相似圖形的入射角度;由于圓(或者橢圓)表示的是一維長(zhǎng)度量,因此不需要考慮圖像縱坐標(biāo)的標(biāo)度.自相似放大區(qū)域占各級(jí)圖形的比例較低(約占22%),從第三層圓開始其位置已無(wú)法仔細(xì)分辨.

為使圖形易于觀察,在圖12(a)的基礎(chǔ)上將圓的直徑各級(jí)放大(圓心位置不變)得到新的結(jié)構(gòu)D′i=Di×ai-1(4.6≥ a≥1).圖 12(b)為 a=3時(shí)的“眼式結(jié)構(gòu)”,圖中的1,2兩個(gè)圓已經(jīng)相交,第2層及以后的圓仍為嵌套關(guān)系,由此說(shuō)明圖12(a)中1,2兩層的自相似性并不好.結(jié)合圖11的自相似圖形可以發(fā)現(xiàn),其中的前兩個(gè)圖形黑白區(qū)域分布有較明顯差異,而其他圖形之間無(wú)明顯差異.由此得到結(jié)論：逃逸曲線的分形自相似結(jié)構(gòu)中,各級(jí)放大區(qū)域在位置標(biāo)度上也是自相似的,其相似度的好壞與逃逸曲線的相似性一致,各級(jí)放大區(qū)域在位置標(biāo)度上的相似度越好逃逸曲線的自相似性越好.

統(tǒng)計(jì)圖11中各級(jí)放大區(qū)域以及相鄰放大區(qū)域間的比值(簡(jiǎn)稱相似比),得到表1.表中給出了各級(jí)相似結(jié)構(gòu)的放大區(qū)域以及相鄰√放大級(jí)之間的相似比比值.利用標(biāo)準(zhǔn)差公式算得表1中相似比的標(biāo)準(zhǔn)差為0.0062,因此可以認(rèn)為相鄰放大區(qū)域間的比值是基本相等的.由此得到結(jié)論：分形自相似結(jié)構(gòu)中,相鄰放大級(jí)之間的相似比基本一致,在本結(jié)構(gòu)中為21.71%.此方法可以用來(lái)檢驗(yàn)自相似結(jié)構(gòu)中放大級(jí)是否存在缺失,對(duì)于放大級(jí)的準(zhǔn)確尋找具有一定的意義.

圖12 (a)原始尺度的“眼式結(jié)構(gòu)”;(b)a=3時(shí)的“眼式結(jié)構(gòu)”

表1 各放大級(jí)區(qū)域及相鄰放大級(jí)之間的比值

5 逃逸曲線的分形維數(shù)及討論

分形幾何與混沌動(dòng)力學(xué)具有密切的聯(lián)系,前文得到的混沌逃逸曲線需要分形幾何來(lái)描述.分形維數(shù)是衡量一個(gè)體系混沌性質(zhì)的重要參量,可以定量描述混沌性質(zhì)的強(qiáng)弱.對(duì)于不規(guī)則形狀的曲線一般采用盒計(jì)數(shù)法來(lái)計(jì)算其分形維數(shù),將混沌逃逸曲線入射角平均細(xì)分為ε份,找出混沌區(qū)域所占的份數(shù)N(ε),則分形維數(shù)在包含混沌區(qū)域的N(ε)個(gè)盒子中,會(huì)出現(xiàn)全部混沌區(qū)域或者部分混沌區(qū)域兩種情況,當(dāng)ε改變時(shí)N(ε)也隨之改變,因此計(jì)算得到的分形維數(shù)Df是變化的.從長(zhǎng)遠(yuǎn)角度考慮,ε越大N(ε)中全部混沌區(qū)域所占的比例也就越多,當(dāng)ε→∞時(shí)所有N(ε)個(gè)盒子均為全部混沌區(qū)域,體系的分形維數(shù)Df趨于穩(wěn)定值[17].

圖13 采用盒計(jì)數(shù)法計(jì)算原器件的分形維數(shù)

圖13 為采用盒計(jì)數(shù)法算得的原器件的分形維數(shù),逐步增大盒子的數(shù)量ε,曲線趨于穩(wěn)定的斜率k=0.734即為體系的分形維數(shù)Df.由于盒計(jì)數(shù)法中的混沌盒子數(shù)量N(ε)不會(huì)超過(guò)盒子的總數(shù)量ε,所以采用這種方法計(jì)算得到的分形維數(shù)Df≤1.采用該方法計(jì)算其他參數(shù)下器件中粒子逃逸曲線的分形維數(shù),進(jìn)而可以得到分形維數(shù)與器件開口、腔長(zhǎng)、圓弧半徑等的關(guān)系曲線,如圖14所示.

圖14(a)中開口寬度d→0時(shí)分形維數(shù)下降非常快,這主要是由于開口寬度很小時(shí),粒子在各角度入射時(shí)均難以逃逸,碰撞次數(shù)奇變少所導(dǎo)致.由于開口很小時(shí)粒子無(wú)法入射和逃逸器件腔,所以此段曲線在實(shí)際應(yīng)用中可以忽略.此后隨著開口寬度的增加體系混沌性減弱,分形維數(shù)Df逐漸降低趨向于0.圖14(b)中腔長(zhǎng)的增加使粒子的混沌性增強(qiáng),分形維數(shù)Df逐漸增加趨向于1.圖14(c)中在R=(0,0.20)區(qū)間分形維數(shù)隨圓弧半徑增加(弧度不變)而增加,這主要是圓弧的放大作用增強(qiáng)所導(dǎo)致.在R=(0.20,0.35)區(qū)間,器件拐角區(qū)域合并導(dǎo)致混沌性減弱,而圓弧半徑增加導(dǎo)致混沌性增強(qiáng),在二者共同影響下分形維數(shù)并沒(méi)有明顯增加.圖14(d)中圓弧半徑增加(弧度減小)導(dǎo)致其放大作用下降以及混沌性減弱,分形維數(shù)隨之減小.可以看到不同器件參數(shù)下的分形維數(shù)變化情況與第3節(jié)中討論的逃逸曲線的混沌性質(zhì)是一致的.

圖14 分形維數(shù)與器件參數(shù)的關(guān)系曲線

6 結(jié)論

1)RIKEN器件的開口寬度、腔長(zhǎng)、圓弧、拐角等均會(huì)影響粒子逃逸曲線的混沌性質(zhì),圓弧區(qū)域(主要是弧度)對(duì)混沌性質(zhì)的影響比直線區(qū)域更顯著,器件拐角的數(shù)量及距離都會(huì)對(duì)混沌性質(zhì)產(chǎn)生影響.

2)通過(guò)對(duì)器件參數(shù)的控制改變器件的混沌性,對(duì)器件中粒子輸運(yùn)過(guò)程的穩(wěn)定性控制有一定意義.

3)采用逃逸曲線定性比較和分形維數(shù)定量計(jì)算兩種方法均可以描述體系的混沌性質(zhì),兩種方法得到的混沌性質(zhì)變化規(guī)律是一致的.

4)在RIKEN器件的粒子逃逸曲線中發(fā)現(xiàn)了很好的分形自相似結(jié)構(gòu),首次運(yùn)用“眼式結(jié)構(gòu)”分析法及相似比比較的方法對(duì)自相似結(jié)構(gòu)進(jìn)行了驗(yàn)證,得到了自相似結(jié)構(gòu)與其位置和相似比的關(guān)聯(lián)性.

[1]Hu G 1994 Random Force and Nonlinear Systems(Shanghai：Shanghai Technology and Education Press)p114(in Chinese)[胡崗 1994 隨機(jī)力與非線性系統(tǒng)(上海：上?？萍冀逃霭嫔?第114頁(yè)]

[2]Zhong Y X 2010 Discussion on Chaos and Fractals(Beijing：Peking University Press)p5(in Chinese)[鐘云霄2010混沌與分形淺談(北京：北京大學(xué)出版社)第5頁(yè)]

[3]Hansen P,Mitchell K A,Delos JB 2006 Phys.Rev.E 73 066226

[4]Ree S,Reichl L E 2002 Phys.Rev.E 65 055205

[5]Micolich A P,Taylor RP,Newbury R,Fromhold T M,Tench CR 2000 Europhys.Lett.49 417

[6]Jiang GH,Zhang Y H,Cai X J,Yang QN 2011 Shandong Science24 22(in Chinese)[蔣國(guó)輝,張延惠,蔡祥吉,楊秦男2011山東科學(xué)24 22]

[7]Micolich A P 2000 Ph.D.Dissertation(Australia：The University of New South Wales)

[8]Fromhold TM,Tench CR,Taylor RP,Micolich A P,Newbury R 1998 Physica B 249 334

[9]Zhao H J,Du M L 2007 Acta Phys.Sin.56 3827(in Chinese)[趙海軍,杜孟利2007物理學(xué)報(bào)56 3827]

[10]Zhao H J,Du M L 2007 Phys.Rev.E 76 027201

[11]Song X F,Du M L,Zhao H J 2012 Sci.Sin.Phys.Mech.Astron.42 127(in Chinese)[宋新芳,杜孟利,趙海軍2012中國(guó)科學(xué)：物理學(xué)力學(xué) 天文學(xué)42 127]

[12]Xu X Y,Zhang Y H,Huang F Z,Lin SL,Du M L 2005 Acta Phys.Sin.54 4538(in Chinese)[徐學(xué)友,張延惠,黃發(fā)忠,林圣路,杜孟利2005物理學(xué)報(bào)54 4538]

[13]Jiang G H,Zhang Y H,Bian H T,Xu X Y 2011 Chin.Phys.Lett.28 120507

[14]Mitchell K A,Handley JP,Tighe B,Knudson SK,Delos JB 2003 Chaos13 880

[15]Mitchell K A,Handley JP,Knudson S K,Delos JB 2003 Chaos 13 892

[16]Huang R S 2000 Chaos and Its Application(Wuhan：Wuhan University Press)p179(in Chinese)[黃潤(rùn)生2000混沌及其應(yīng)用(武漢：武漢大學(xué)出版社)第179頁(yè)]

[17]Taylor RP,Micolich A P,Jonas D 1999 Nature 399 422