羅曉華 何為

(1.重慶交通大學(xué) 圖書館，重慶 400074;2.重慶大學(xué) 電氣工程學(xué)院，重慶 400044)

其中e是電子電荷，h是諧波數(shù)，V是加速電壓，E是同步粒子能量，β=v/c是粒子無量綱速度。不失一般性，考慮φs=0的穩(wěn)態(tài) (stationary state)情形。當(dāng)φs=0時，方程 (1)退化為

1944年，蘇聯(lián)物理學(xué)家維克斯勒 (次年，麥克米蘭)發(fā)現(xiàn)，在高頻電場中運動的帶電粒子具有保持相位穩(wěn)定的能力，并指出適當(dāng)選擇同步相位，即可實現(xiàn)非同步粒子的能量補償。從此，加速器的能量上限實現(xiàn)了質(zhì)的突破，質(zhì)子能量從當(dāng)初的107eV量級提升到當(dāng)今的1013eV量級。遺憾的是束流能量和束流強度不能兼得。正是這對矛盾激勵著人們不斷對加速器技術(shù)進行探索，并取得了一個又一個的突破，儲存環(huán)技術(shù)和束流冷卻技術(shù)便是典型的例子。直到今天，加速器的束流動力學(xué)問題依然是人們關(guān)注的重要問題之一［1－14］。本文關(guān)心的正是束流冷卻對同步運動的影響。值得注意的是，隨著人們認識能力的不斷提高，先后認識了束流動力學(xué)的線性、非線性和復(fù)雜性，而本文就試圖利用Mathieu方程的穩(wěn)定區(qū)和不穩(wěn)定區(qū) (禁帶)來討論系統(tǒng)的穩(wěn)定性。為此，在經(jīng)典力學(xué)框架內(nèi)，把周期冷卻的同步運動方程化為Mathieu方程，系統(tǒng)呈現(xiàn)出一系列穩(wěn)定區(qū)和不穩(wěn)定區(qū);對可能存在的一階不穩(wěn)定區(qū)進行了討論，用攝動法導(dǎo)出了它的邊界曲線和禁帶寬度，并討論了禁帶穿越和同步運動的穩(wěn)定性，導(dǎo)出了束流冷卻一次諧波振幅的臨界值。

1 運動方程

束流冷卻通常是在環(huán)形加速器 (比如同步加速器、儲存環(huán)等)的一個或幾個直線段上安裝冷卻器來實現(xiàn)的。當(dāng)粒子通過冷卻器時，用一束冷 (發(fā)射度小)的電子同它一道運動，并將整個離子束浸泡其中。如果選擇電子的速度與離子速度相等，離子和電子在運動方向保持同步，同時它們之間還將通過“熱”運動彼交換橫向能量，使得比較熱的離子束逐漸變冷。換句話說，由于束流冷卻，離子的Betatron振蕩將逐漸衰減，而同步振蕩也將逐漸衰減。

在理想情況下，如果選擇相位φ和動量分散系數(shù)δ為正則坐標(biāo)，θ為獨立變數(shù)，粒子同步運動方程可表示為［7，13－14］

其中e是電子電荷，h是諧波數(shù)，V是加速電壓，E是同步粒子能量，β=v/c是粒子無量綱速度。不失一般性，考慮φs=0的穩(wěn)態(tài) (stationary state)情形。當(dāng)φs=0時，方程 (1)退化為

設(shè)系統(tǒng)的冷卻系數(shù)為μ(θ)，對于動量分散為δ的粒子束，經(jīng)過“長度”為Δθ的冷卻后，衰減量Δδ可表示為 Δδ= －μ(θ)δΔθ。于是，正則方程 (2)可化為［7］

或

其中

其中ε是小參數(shù)，僅表示伴隨項是O(ε)級小量。將μ(θ)按Fourier級數(shù)展開

其中Δθ是冷卻器“長度”，而η0=Δθ/2π是冷卻器的相對長度。令

方程 (5)可進一步化為

其中

式中略去了ε2項。這就是只關(guān)心弱冷卻情況，而μ(θ)由式 (6)給出。將式 (6)和 (9)帶入方程(8)，可得粒子同步運動方程

其中

是由于周期冷卻產(chǎn)生的n次諧波振幅。方程 (10)是一般的Hill方程，它的系數(shù)包含了所有Fourier分量，這個系統(tǒng)將存在和共振與差共振。但是由于這些組合共振都比較弱，感興趣的是整數(shù)共振或半整數(shù)共振。由式 (10)可得n次諧波激勵的同步運動方程為

2 系統(tǒng)的穩(wěn)定性

作變數(shù)變換

方程 (12)化為

其中參數(shù)

方程 (14)是一個標(biāo)準(zhǔn)的Mathieu方程。一次項的系數(shù)`是φ的周期函數(shù)。根據(jù)Flogue定理，方程 (14)的解具有如下性質(zhì)

其中d是加速器格子周期，而T是粒子運動周期。粒子在周期中場運動時，它的能量將分裂成帶。事實上，系統(tǒng)的穩(wěn)定性與參數(shù) (q，ε)的分布有關(guān)，在參數(shù) (q，ε)平面上，系統(tǒng)出現(xiàn)了一系列穩(wěn)定和不穩(wěn)定區(qū)。當(dāng)→0時，如果q＞0，系統(tǒng)是穩(wěn)定的;q＜0系統(tǒng)不穩(wěn)定。不過，隨著的增加，即使q＜0，系統(tǒng)也可能是穩(wěn)定的。由于方程 (14)包含振動部分cos 2φ，在q＞0的右半平面內(nèi)出現(xiàn)了一系列穩(wěn)定區(qū)和不穩(wěn)定區(qū)。當(dāng)→0時，這些不穩(wěn)定區(qū)退化為一點，且全都落在ε=0的橫軸上，這些點由方程

給出，并分別稱為一階、二階、三階等不穩(wěn)定區(qū)。對于一次諧波 (n=1)激勵，由式 (15)，(16)可知，這些不穩(wěn)定區(qū)對應(yīng)于νs=1/2，1，3/2，…，共振。這就是說，由于束流冷卻將激發(fā)出不同階次的整數(shù)共振或半整數(shù)共振。圖1給出了數(shù)值結(jié)果。一般說來，在參數(shù) (q，ε)平面上，這種穩(wěn)定和不穩(wěn)定區(qū)有無限多個。對于具體情況只有一個或少數(shù)幾個穩(wěn)定區(qū)和不穩(wěn)定區(qū)起作用。下面對一階不穩(wěn)定區(qū)(νs=1/2)做進一步分析。

圖1 Mathieu方程的穩(wěn)定與不穩(wěn)定區(qū) (陰影區(qū)域)

3 Mathieu方程的一階不穩(wěn)定區(qū)及其禁帶寬度

利用攝動法［15－16］求解方程 (14)。注意到方程 (14)中的 (是一個小參數(shù)，可將x和q按它展開

將式 (17)和 (18)代入方程 (14)，比較 (的同次冪系數(shù)，可得逐次近似方程

方程 (19)的解可表示為

上式表明，當(dāng)k為偶數(shù)時，x0是周期為π的周期函數(shù)，當(dāng)k為奇數(shù)時，x0是周期為2π的周期函數(shù)。我們以k=1為例，求禁帶的邊界方程和相應(yīng)的解。

將式 (22)代入方程 (20)，并令k=1，可得

消去方程中的久期項，要求

于是方程 (23)化為

聯(lián)立求解方程 (24)和 (25)，可得

和

航空發(fā)動機裝配過程數(shù)據(jù)能夠映射為任務(wù)、物料、工藝、質(zhì)量4個視圖上的對應(yīng)信息，分別回答了“由誰在何時來執(zhí)行哪個任務(wù)”、“對哪些物料進行裝配”、“如何裝配”和“記錄哪些質(zhì)量信息”的問題，4個視圖之間存在的邏輯關(guān)聯(lián)關(guān)系如圖2所示。

于是，方程 (26)的解可表示為

或

將式 (27)、(29)代入方程 (21)，可得二級近似方程

消去久期項，可得

精確到二次近似，一階不穩(wěn)定區(qū)的一條邊界曲線和這條邊界曲線上的解可表示為

同樣，將式 (28)、(30)代入方程 (21)，并用類似的方法，可求得一階不穩(wěn)定區(qū)的另一條邊界曲線和這條邊界曲線上的解為

由式 (33)和 (35)，可求得一階不穩(wěn)定區(qū)的禁帶寬度

4 系統(tǒng)的臨界性質(zhì)

注意到當(dāng)諧波數(shù)n=1時，由式 (15)可知，同步運動的Betatron頻率可表示為=q/4。微分該式，由式 (37)可將一階不穩(wěn)定區(qū)的禁帶寬度表示為

再注意到粒子穿越禁帶時振幅呈指數(shù)增長，設(shè)φ0(因而x0)是粒子進入禁帶時的“初始”振幅，穿越禁帶后的振幅φ(因而x)可表示為

其中θ是粒子在禁帶中轉(zhuǎn)過的角度 (或“停留”的時間)。當(dāng)φ滿足條件

時，系統(tǒng)處于臨界狀態(tài)，其中φu和π－φs是同步運動的兩個不穩(wěn)定點，φs是系統(tǒng)的同步相位 (在本文中φs=0)。換句話說，當(dāng)離子穿越禁帶，相位呈指數(shù)增長，如果相位的振幅φ＜φu和φ＜π－φs，系統(tǒng)是穩(wěn)定的，否則系統(tǒng)不穩(wěn)定。由式 (15)、(16)和 (38)可得束流冷卻的一次諧波振幅臨界值

或

其中ωs是同步粒子振蕩頻率，Tc是同步粒子穿越禁帶的臨界時間。當(dāng)δ1＜δ1c系統(tǒng)是穩(wěn)定的，否則系統(tǒng)不穩(wěn)定。

5 結(jié)語

由于冷卻效應(yīng)，同步運動方程化為了參數(shù)激勵的Matheiu方程，由于參數(shù)激勵，系統(tǒng)出現(xiàn)了一系列穩(wěn)定和不穩(wěn)定區(qū)。這就是說，由于束流冷卻，系統(tǒng)的同步運動將在νs=1/2，1，3/2，…，處出現(xiàn)共振。本文從系統(tǒng)的正則方程出發(fā)，在經(jīng)典力學(xué)框架內(nèi)和弱冷卻近似下，把粒子的同步運動方程化為Mathieu方程，利用攝動法導(dǎo)出了一階不穩(wěn)定區(qū) (νs=1/2)的邊界曲線和禁帶寬度。指出了由于冷卻效應(yīng)，當(dāng)粒子穿越禁帶后的相位φ=φ0+Δφ等于不穩(wěn)定相位時，粒子處于臨界狀態(tài)，并由此導(dǎo)出了束流冷卻的一次諧波激勵振幅的臨界值或。當(dāng) δ1＜ δ1c時，系統(tǒng)是穩(wěn)定的，否則系統(tǒng)不穩(wěn)定。

［1］Dumas H S，Ellison J A，Vogt M.First－Order Averaging Principles for Maps with Applications to Accelerator Beam Dynamics［J］.SIAM J Appl Dyn Syst，2004，3(4):409 －432.

［2］羅詩裕，邵明珠，胡西多.二維晶化束的平均場概念與單粒子模型(I)［J］.高能物理與核物理，2004，28(1):96－99.

［3］胡西多，羅詩裕，邵明珠，等.二維晶化束的平均場概念與單粒子模型(II)［J］.高能物理與核物理，2004，28(2):196－199.

［4］Derbenev Y.Six － dimensional muon beam cooling using a homogeneous absorber:Concepts，beam dynamics，cooling decrements，and equilibrium emittances in a helical dipole channel［J］.Phys Rev ST Accel Beams，2008，5:041002.(20 pages)

［5］Qiang J，Lidia S，Ryne R D，et al.Three－dimensional quasistatic model for high brightness beam dynamics simulation［J］.Phys Rev ST Accel Beams，2006，9:044204.(10 pages)

［6］Helmut Wiedemann.Particle accelerator physics［M］.Berlin:Springer，2007:215 －221.

［7］Lee S Y.Accelerator Physics［M］.2nd ed.Shanghai:Fudan University Press，2006:288 －295.

［8］Schroeder C B，Shadwick B A，Esarey E.Radiative damping and electron beam dynamics in plasma－ based accelerators［J］.Phys Rev E，2006，74:026501.(14 pages)

［9］Ferrario M，Alesini D，Bacci A，et al.Measurement of the Double Emittance Minimum in the Beam Dynamics of the Sparc High－Brightness Photoinjector［J］.Phys Rev Lett，2007，99:234801.(5 pages)

［10］Stokes R H，Crandall K R，Stovall J E.Rf quadrupole beam dynamics［J］.IEEE Transactions on Nuclear Science，2007，26(3):3469－3471.

［11］Zhang Z L，Jameson R A，Zhao H W，et al.Beam dynamics design of an RFQ for a planned accelerator，which uses a direct plasma injection scheme［J］.Nuclear Instruments and Methods in Physics Research Section A:Accelerators，Spectrometers，Detectors and Associated Equipment，2008，592(3):197 －200.

［12］Piel C，Dunkel K，Kremer F，Pekeler M.Phase 1 Commissioning Status of the 40 MeV Proton/Deuteron Accelerator SARAF［C］//Proceeding of EPAC08.Genoa:2008.

［13］范麗仙，羅詩裕，邵明珠.正弦平方勢與環(huán)形加速器的彎晶束流引出［J］.原子核物理評論，2011，28(1):63－67.

［14］肖慧娟，羅詩裕，邵明珠.多尺度法與準(zhǔn)等時同步加速器粒子縱向運動的穩(wěn)定性［J］.原子核物理評論，2011，28(3):300－304.

［15］Nayfeh A H.Introduction to Peturbation Techniques［M］.New Youks:John Wiley＆ Sons，1981:226－240.

［16］羅詩裕，邵明珠，羅曉華.正弦平方勢與應(yīng)變超晶格位錯動力學(xué)［J］.中國科學(xué):物理學(xué)，力學(xué)和天文學(xué)，2010，40(2):207－212.