徐正一,彭欣欣,李連花,周尹敏,裘栩煬,周 敏,徐信業(yè)
(華東師范大學(xué)精密光譜科學(xué)與技術(shù)國家重點實驗室,上海 200062)
核磁共振陀螺儀因其兼具小型化與高精度的特點,可以在一些特殊應(yīng)用中替代較為昂貴的轉(zhuǎn)子陀螺。核磁共振陀螺儀中沒有運動部件,對于振動和過載都不敏感,擁有低功耗和低成本的巨大潛力,并正在向芯片級尺寸和戰(zhàn)略級精度發(fā)展,目前已成為國內(nèi)外研究的熱點領(lǐng)域,且未來有望應(yīng)用于航空航天以及國防戰(zhàn)略領(lǐng)域,同時在民用領(lǐng)域也具有廣泛的應(yīng)用前景。近年來,北京自動化控制設(shè)備研究所[1-2]和國防科技大學(xué)[3-6]等單位也開始了核磁共振陀螺儀的理論和實驗研究。
核磁共振陀螺儀選用極化后的原子核自旋在外磁場中的拉莫爾進動作為轉(zhuǎn)動參考系,因此縱向磁場的不穩(wěn)定會造成拉莫爾進動頻率的變化;而在整個陀螺儀旋轉(zhuǎn)時,由于緩沖氣體的黏滯性,其中的氣體可能并不跟隨陀螺儀同步轉(zhuǎn)動,這種情況可能會使陀螺儀的動態(tài)范圍受限。為了盡可能減少上述問題帶來的影響,核磁共振陀螺儀需要設(shè)計磁屏蔽系統(tǒng)和磁場補償系統(tǒng)以消除外部磁場干擾并維持內(nèi)部磁場的均勻性。
三維磁場鎖定的關(guān)鍵是三維原子磁力計的實現(xiàn)。傳統(tǒng)的基于非線性旋光效應(yīng)的原子磁力計將信號對其他方向上磁場的響應(yīng)視作干擾,因此只能實現(xiàn)一個方向上磁場的探測,而通過對不同方向上的磁場以不同的頻率分別進行調(diào)制,就可以在傳統(tǒng)原子磁力計的基礎(chǔ)上實現(xiàn)三維磁場探測。通過調(diào)制x和z方向磁場,在無外磁場屏蔽的情況下,可實現(xiàn)靈敏度達到1pT·Hz-1/2量級的三維原子磁力計[7];使用雙方向泵浦光頻率調(diào)制的方案,磁力計靈敏度可達到65fT·Hz-1/2[8]?;趦深惒煌脑哟帕τ媽崿F(xiàn)的三維磁力計,可實現(xiàn)x和y方向300fT·Hz-1/2的靈敏度,z方向3pT·Hz-1/2的靈敏度[9]。
考慮到核磁共振陀螺儀小型化的發(fā)展方向,選用兩類不同的原子磁力計實現(xiàn)現(xiàn)有系統(tǒng)中三維磁場的同時測量。下面首先介紹核磁共振陀螺儀的基本原理以及三維磁場鎖定的基本原理,分析三維磁力計誤差信號并給出目前初步實現(xiàn)的三維磁場鎖定情況,最后對現(xiàn)階段實驗的優(yōu)化方向進行分析和討論。
原子核自旋在外加靜磁場中會進行拉莫爾進動。如果對原子施加一定頻率的電磁輻射,且其頻率接近原子核自旋的拉莫爾進動頻率時,原子核自旋會吸收電磁輻射的能量,產(chǎn)生核磁共振現(xiàn)象。這類利用核磁共振過程中的拉莫爾進動作為轉(zhuǎn)動參考系實現(xiàn)的陀螺儀即核磁共振陀螺儀,其工作原理如圖1所示。
圖1 核磁共振陀螺儀的工作原理Fig.1 Operation principle of the nuclearmagnetic resonance gyroscope
核磁共振陀螺儀的基本原理主要包括光泵浦超極化和高靈敏原子磁力計,工作原子一般選用堿金屬原子和惰性氣體原子。首先通過施加同向傳播的圓偏振泵浦光和縱向靜磁場B0以極化堿金屬原子,然后通過堿金屬原子與惰性氣體原子之間的自旋交換碰撞過程實現(xiàn)惰性氣體原子的極化。在垂直于泵浦光傳播方向上施加激勵磁場B1后,惰性氣體原子極化形成的宏觀磁矩進行拉莫爾進動,如果惰性氣體原子的旋磁比為γn,則拉莫爾進動頻率ωL=γnB0。此時在垂直于泵浦光傳播方向上施加探測光后,基于非線性旋光效應(yīng)可以實現(xiàn)高靈敏原子磁力計,用于探測惰性氣體原子極化形成的宏觀磁矩。當(dāng)陀螺儀以角速度ωR繞泵浦光方向旋轉(zhuǎn),則探測到的宏觀磁矩進動角速度ωM=ωL+ωR。故陀螺儀的角速度ωR=ωM-ωL,從而實現(xiàn)了陀螺儀對角速度的測量。可見磁場的穩(wěn)定性對于核磁共振陀螺儀的性能至關(guān)重要,實現(xiàn)三維剩磁的補償有利于提高核磁共振陀螺儀的性能。
如圖1所示,堿金屬原子形成的宏觀磁矩M可以通過布洛赫方程進行描述
(1)
其中,M0為未進動時宏觀磁矩的大小,γe為堿金屬原子的旋磁比,T1和T2分別為堿金屬原子的縱向和橫向弛豫時間。如果令M+=Mx+iMy,B+=Bx+iBy,則關(guān)于M+的布洛赫方程可以寫為
(2)
當(dāng)B0被調(diào)制后,即Bz=B0+Bmcos(ωmt),ωm和Bm分別為靜磁場B0的調(diào)制頻率和幅度,則上述布洛赫方程的解可以寫為[10-11]
(3)
其中,km=γeBm/ωm。當(dāng)滿足|Bx|<<1/γeT2,|By|<<1/γeT2和|Bz-nωm/γe|<<1/γeT2時,可進一步給出δMx和δMy
δMx=γeMzT2Jn(km)×
δMy=γeMzT2Jn(km)×
(4)
由于km與Bm和B0的比值成正比,故通過改變Bm和B0的比值可以改變δMx和δMy的大小。選取n=1可得
(5)
以ωm作為本征頻率,對于y方向上探測到的磁力計信號,如果選用1倍或者2倍頻解調(diào),即p=1或p=2,則通過鎖相放大器解調(diào)得到的X和Y通道的信號SX1和SY1或SX2和SY2分別為
(6)
以ωm的調(diào)制頻率調(diào)制B0后,使用鎖相放大器對原子磁力計探測到的信號進行解調(diào),通過調(diào)節(jié)移相器的相位可以分別獲得平行和垂直于探測光方向上剩磁。而在不同的Bm/B0下,應(yīng)該選用不同的p來優(yōu)化鎖相放大器輸出的Bx和By的誤差信號。當(dāng)Bm/B0≈1.733時,Bx和By的誤差信號分別選用1倍頻和2倍頻解調(diào)能夠獲得較為理想的誤差信號,此時2個誤差信號的幅值較大且大小相仿;當(dāng)Bm/B0<0.5時,Bx和By的誤差信號都選用1倍頻解調(diào)能夠獲得大小相仿的信號,但信號相對較小。
實現(xiàn)Bz的鎖定需要借助于兩種惰性氣體同位素。如果它們的旋磁比分別為γ1和γ2,原子磁力計探測到的拉莫爾進動頻率分別為ωM1和ωM2,則根據(jù)核磁共振陀螺儀的原理,可以得到
(7)
其中,ωR為陀螺儀裝置轉(zhuǎn)動的角速度。故可以獲得Bz的誤差信號Sz為
(8)
故通過測量兩種惰性氣體同位素的拉莫爾進動頻差,其與參考頻率的頻率差即可用于Bz的鎖定。
實驗裝置如圖2所示,實驗中選用充有133CS原子飽和蒸氣、5torr129Xe、45torr131Xe以及200torr N2的1 cm3立方體氣室,并將其加熱至110℃。氣室放置于屏蔽因子達106的五層磁屏蔽筒內(nèi),屏蔽筒內(nèi)通過三維線圈施加磁場。在z方向上施加的靜磁場為10μT,調(diào)制幅度和調(diào)制頻率分別為3.3μT和35kHz,在x方向上施加一個微弱的激勵磁場,驅(qū)動頻率為約118.6Hz,對應(yīng)129Xe在10μT下的拉莫爾進動頻率。使用自制的激光器提供120mW的895nm圓偏振泵浦光以及500μW的852nm線偏振探測光,泵浦光與探測光的傳播方向分別沿+z和-y方向。對探測光的偏振檢測通過沃拉斯頓棱鏡和雙平衡探測器實現(xiàn)。
圖2 實驗裝置示意圖(HP為半波片,QP為四分之一波片,GP為格蘭棱鏡,WP為沃拉斯頓棱鏡,CS為恒流電源,LIA為鎖相放大器)Fig.2 Experimental setup( HP: half-wave plate, QP: quarter-wave plate, GP: Glan-Talyer prism, WP: Wollastonprism, CS: current source, LIA: lock-in amplifier)
雙平衡探測器的信號分成3路送至3臺獨立的鎖相放大器解調(diào),分別獲得3個方向磁場的誤差信號。其中,LIA1和LIA2均以35kHz作為本征頻率,LIA3以118.6Hz作為本征頻率。因為實驗中所使用的Bm/B0=0.17,因此3臺鎖相放大器均輸出1倍頻解調(diào)信號以獲得x和y方向上較為均衡的信號。通過優(yōu)化LIA1和LIA2的相位,可以獲得x和y方向上磁力計信號對Bx和By的響應(yīng)曲線,如圖3和圖4所示。當(dāng)|δBx|和|δBy|均不大于1μT時,x和y方向上的磁力計信號不同時對Bx和By響應(yīng),故可用于分辨剩磁的產(chǎn)生方向。
當(dāng)陀螺儀靜止時,有ωR=0,故δBz=δωM/γe。通過LIA3可給出Bz的響應(yīng)曲線,當(dāng)|δBz|不大于15nT時,LIA3的輸出信號可用于Bz的鎖定,如圖5所示。
圖6中對比了鎖定前后的三維磁力計信號??梢园l(fā)現(xiàn)三維磁場鎖定后,Bx、By和Bz在約5000s內(nèi)的起伏范圍分別不大于6.75nT、5.33nT和6.68nT;而鎖定前Bx、By和Bz在約4000s內(nèi)的起伏范圍分別約15.05nT、39.78nT和16.20nT。
圖3 110℃的氣室溫度下,x和y方向磁力計對Bx的響應(yīng)曲線Fig.3 Magnetometer responses in thexandyaxes toBxat the cell temperature of 110℃
圖4 110℃的氣室溫度下,x和y方向磁力計對By的響應(yīng)曲線Fig.4 Magnetometer responses in thexandyaxes toByat the cell temperature of 110℃
圖5 110℃的氣室溫度下,z方向上磁力計對Bz的響應(yīng)曲線Fig.5 Magnetometer response in thezaxes toBzat the cell temperature of 110℃
圖6 110℃的氣室溫度下,磁場未鎖定和鎖定時Bx、By和Bz的誤差信號Fig.6 Bx,ByandBzerror signals when the magneticfield is locked and unlocked at the cell temperature of 110℃
如圖7、 圖8和圖9所示,計算磁場鎖定前后Bx、By和Bz的誤差信號對應(yīng)的阿倫方差后可見,磁場鎖定后3個方向上的磁場長期穩(wěn)定性都得到了顯著提升,千秒穩(wěn)均提升至少1個量級,但秒穩(wěn)在鎖定前后沒有顯著變化,這主要受限于驅(qū)動三維線圈的電流噪聲。
圖7 110℃的氣室溫度下,磁場鎖定與未鎖定時Bx誤差信號的阿倫方差Fig.7 Allan deviation of theBxerror signal when the magneticfield is locked and unlocked at the cell temperature of 110℃
圖8 110℃的氣室溫度下,磁場鎖定與未鎖定時By誤差信號的阿倫方差Fig.8 Allan deviation of theByerror signal when the magneticfield is locked and unlocked at the cell temperature of 110℃
圖9 不同氣室溫度下磁場鎖定與未鎖定時Bz誤差信號的阿倫方差Fig.9 Allan deviation of theBzerror signals when the magneticfield is locked and unlocked at various cell temperatures
實驗過程中可以觀察到氣室溫度的選擇對于三維磁場鎖定具有較大的影響。相比110℃和98℃,63℃下Bz的響應(yīng)曲線信噪比相對較差,而Bx和By的響應(yīng)曲線信噪比并沒有較大的變化。這是因為在較低的氣室溫度下,CS原子與Xe原子的自旋交換碰撞弛豫較小,造成Xe極化率較小,從而影響測量Bz的磁力計信噪比。由于泵浦光功率較大,因此堿金屬原子的泵浦速率足夠大,不同氣室溫度下堿金屬原子都能保持相對較高的極化率,故對Bx和By的磁力計誤差信號影響較小。
圖9中給出了不同氣室溫度下z方向磁場鎖定后磁力計誤差信號的阿倫方差,可見63℃的氣室溫度下Bz鎖定后在秒穩(wěn)和千秒穩(wěn)均不及98℃以及110℃,而110℃下z方向磁場秒穩(wěn)更優(yōu)于98℃,千秒穩(wěn)則兩者較為接近。這是由于當(dāng)泵浦速率足夠時,較高的氣室溫度下有著更高的堿金屬原子數(shù)密度以及更為劇烈的自旋交換碰撞,因此測量Bz的磁力計有著更好的信噪比,磁場鎖定的效果也越好。
在不同的氣室溫度下,Bz的主要噪聲類型為磁場的白噪聲,其受限于系統(tǒng)內(nèi)的磁場本底噪聲。在較高的氣室溫度下,磁場本底噪聲更小,因而更有利于磁場的鎖定。根據(jù)實驗經(jīng)驗,在較低的氣室溫度下,用于Bz鎖定的伺服系統(tǒng)參數(shù)也更難整定,調(diào)節(jié)不當(dāng)會造成Bz鎖定后出現(xiàn)正弦噪聲。
本文實現(xiàn)了三維磁場的初步鎖定,驗證了基于縱向磁場調(diào)制的三維磁力計的可行性,并對磁場鎖定的性能進行了測量。通過調(diào)制z方向磁場的方式實現(xiàn)了單探測光束對x和y方向剩磁的測量,通過檢測惰性氣體原子的拉莫爾進動頻率實現(xiàn)了對z方向剩磁的測量,從而實現(xiàn)了三維原子磁力計以及三維磁場的初步鎖定,鎖定后x、y和z方向上的磁場起伏不大于6.75nT、5.33nT和6.68nT。
本文探究了氣室溫度對三維磁場鎖定效果的影響。通過觀察不同氣室溫度對z方向磁場鎖定效果的影響,可以認(rèn)為在泵浦速率足夠且氣室本身造成的退極化弛豫較小的情況下,較高的氣室溫度有利于z方向磁場鎖定。根據(jù)不同氣室溫度下z方向磁場鎖定后磁力計誤差信號的阿倫方差,可以認(rèn)為在更高的氣室溫度下具有更小的磁場本底噪聲,因此更加有利于磁場的鎖定。