金剛石氮空位色心耦合機(jī)械振子和腔場系統(tǒng)中方差壓縮研究?

廖慶洪 葉楊 李紅珍 周南潤

1)(南昌大學(xué)電子信息工程系,南昌 330031)

2)(清華大學(xué),低維量子物理國家重點(diǎn)實(shí)驗(yàn)室,北京 100084)

3)(西安交通大學(xué)應(yīng)用物理系,西安 710049)

1 引 言

隨著納米技術(shù)的不斷進(jìn)步,機(jī)械振子的尺寸已經(jīng)可以做到微米甚至納米量級.納米機(jī)械振子有著較高的振動頻率和品質(zhì)因子、極小的質(zhì)量以及較低的耗散,在高精度位移測量、力測量、質(zhì)量測量和精確的量子計算等方面有著重要的應(yīng)用[1?4].此外,機(jī)械振子可以用來研究宏觀系統(tǒng)的量子特性,比如宏觀量子疊加態(tài)、量子糾纏態(tài)、以及用于量子信息處理等.要讓納米機(jī)械振子在如此廣泛的應(yīng)用中都發(fā)揮作用,人們希望能更加方便地操控其狀態(tài).因此,如何把機(jī)械振子與一個兩能級系統(tǒng)耦合起來,實(shí)現(xiàn)強(qiáng)耦合就很值得研究[5].另一方面,在固態(tài)系統(tǒng)中,金剛石氮空位中心(NV色心)在室溫下具有毫秒量級的長相干時間,也很容易被微波操控,并被激光高效率地讀出其狀態(tài),是實(shí)現(xiàn)量子計算的優(yōu)良載體,被認(rèn)為有望實(shí)現(xiàn)室溫下的量子計算和高靈敏度量子測量以及在量子信息處理、納米尺度空間中的磁場與電場、溫度的靈敏探測等相關(guān)領(lǐng)域具有巨大的應(yīng)用價值[6].把金剛石氮空位中心與納米機(jī)械振子耦合起來,用于精密測量,量子信息處理等成為了一個熱門的研究題目[5].利用梯度磁場誘導(dǎo)的金剛石NV色心與微納米機(jī)械振子運(yùn)動之間的耦合,實(shí)現(xiàn)了百納米尺度的金剛石質(zhì)心振動與其中的色心電子自旋之間的強(qiáng)耦合,并且制備了金剛石NV色心的薛定諤貓態(tài)和實(shí)現(xiàn)雙縫物質(zhì)波干涉[7].文獻(xiàn)[8]提出了基于金剛石NV色心與機(jī)械振子的耦合,利用動力學(xué)解耦技術(shù)實(shí)現(xiàn)超靈敏質(zhì)量探測方案.基于Stark效應(yīng)的高靈敏度電場探測[9],基于NV色心的光檢測磁共振譜峰值隨著溫度的變化而變化,可實(shí)現(xiàn)納米尺度的高靈敏度溫度探測[10],NV色心與外部機(jī)械振子以及所在的金剛石機(jī)械振子耦合分別可以實(shí)現(xiàn)髙靈敏度的機(jī)械振子運(yùn)動[11]以及金剛石內(nèi)部的應(yīng)力等[12]物理量探測.Li等[13]理論上提出了利用通電的碳納米管機(jī)械振子和單個NV色心自旋耦合的混合型量子器件,實(shí)現(xiàn)了NV色心單自旋體系與納米機(jī)械振子強(qiáng)耦合機(jī)制及相干操控.文獻(xiàn)[14]研究了借助兩能級發(fā)射器由真空漲落引起的能級移動探索振動石墨烯機(jī)械振子位移在量子水平的壓縮和測量方案.文獻(xiàn)[15,16]研究了超導(dǎo)量子比特耦合NV色心系綜系統(tǒng)中的量子態(tài)的傳輸和存儲以及存在內(nèi)相互作用的多體開放量子系統(tǒng)中的糾纏動力學(xué)特性.

本文研究金剛石NV色心同時耦合腔場和機(jī)械振子系統(tǒng)中機(jī)械振子的方差壓縮,主要討論金剛石NV色心初態(tài)和NV色心與機(jī)械振子耦合強(qiáng)度對聲子場方差壓縮的影響.結(jié)果發(fā)現(xiàn):控制和操縱系統(tǒng)參數(shù)可以制備壓縮時間長、壓縮幅度大的聲子場壓縮態(tài),并且通過調(diào)控NV色心初態(tài)以及磁場梯度可以實(shí)現(xiàn)對機(jī)械振子方差壓縮非經(jīng)典特性的操控.

2 物理模型及其基本方程

本文研究光懸浮的納米金剛石(內(nèi)含單個NV色心)系統(tǒng)[7],物理模型如圖1所示,納米金剛石囚禁在諧振勢中,囚禁納米金剛石振動頻率為ωm,質(zhì)量為m,納米金剛石質(zhì)心運(yùn)動哈密頓量可以描述為若外加磁場具有梯度,則導(dǎo)致金剛石NV色心與納米金剛石質(zhì)心運(yùn)動的耦合,耦合哈密頓量為考慮金剛石NV色心和腔場相互作用,其相互作用哈密頓量為同時金剛石NV色心由外加激光場驅(qū)動,則系統(tǒng)的哈密頓量為

其中,ω為腔場頻率,ω0為金剛石NV色心躍遷頻率,ωc為經(jīng)典驅(qū)動場頻率;g為金剛石NV色心與腔場的耦合系數(shù);λ為NV色心與經(jīng)典驅(qū)動場的耦合系數(shù);G=gsμBGma0/?為NV色心與機(jī)械振子的耦合系數(shù);Gm為磁場梯度;b+和b分別為腔場的產(chǎn)生和湮滅算符;|0〉和|1〉構(gòu)成量子比特,分別對應(yīng)σz=?1和σz=+1的本征矢[8],令|0〉=|g〉,|1〉=|e〉,則σ+=|e〉〈g|和σ?=|g〉〈e|分別為金剛石NV色心的上升和下降算符.

圖1 金剛石NV色心耦合腔場和機(jī)械振子模型示意圖Fig.1.Schematic diagram of nitrogen-vacancy center in diamond coupled to cavity field and mechanical resonator.

因?yàn)殓壅儞Q不會改變系統(tǒng)的本征值,引入一個幺正變換U1=exp(?iωcσzt/2),將(1)式的哈密頓量變換到旋轉(zhuǎn)坐標(biāo)表象中,則

其中

式中Δ=ω0?ωc為金剛石NV色心躍遷頻率與NV色心驅(qū)動場頻率的失諧量.

由綴飾態(tài)基|+〉,|?〉和|e〉,|g〉之間的關(guān)系可得

則系統(tǒng)在綴飾態(tài)表象下的哈密頓量為

引入幺正操作U2=exp(iωcSzt/2),系統(tǒng)哈密頓量(8)式變?yōu)?/p>

其中

本文主要研究NV色心與機(jī)械振子、腔場之間的相互作用,則考慮系統(tǒng)哈密頓量的相互作用部分,在相互作用繪景下,令U4=exp(?iH(0)4t),可得系統(tǒng)相互作用繪景下的哈密頓量為

其中Δ1=?+ωc? ω,Δ2=? ? ωm.考慮特殊情況Δ1=0,Δ2=0,金剛石NV色心初始處于cosθ1|+〉+sinθ1|?〉,機(jī)械振子處于基態(tài)|0m〉,腔場處于單光子態(tài)|1f〉,則系統(tǒng)初態(tài)為

在哈密頓量(12)式作用下,系統(tǒng)任意時刻的態(tài)矢量為

求解薛定諤方程,可得

金剛石機(jī)械振子的振動模式可以處理成聲子振動模,研究聲子場的壓縮,引入正交算符

根據(jù)方差的定義

其中

由任意時刻態(tài)矢量(14)式,可以得到金剛石機(jī)械振子的約化密度矩陣為

其中

所以

當(dāng)聲子場處于壓縮態(tài)時,其中一個正交算符的量子漲落小于1/2,減小了聲子場某一部分的量子起伏,使測量的精確度提高,由此壓縮態(tài)在精密測量、引力波探測、光通信、量子信息處理等方面具有重要應(yīng)用.

為了更好地闡明金剛石機(jī)械振子方差壓縮的物理機(jī)制和原因,給出機(jī)械振子的相干性的定義和計算表達(dá)式為[19,20]

由(19)和(27)式可得機(jī)械振子的相干性為

利用以上結(jié)果討論金剛石NV色心初態(tài)和NV色心與機(jī)械振子耦合強(qiáng)度對聲子場方差壓縮影響.

3 結(jié)果與討論

圖2展示了聲子場方差壓縮因子ΔX2和相干性隨時間的動力學(xué)演化特性,圖中實(shí)線和虛線分別對應(yīng)聲子場方差壓縮因子ΔX2和相干性演化規(guī)律,其中金剛石NV色心和腔場耦合系數(shù)g=1,θ1=π/3,金剛石NV色心和機(jī)械振子耦合強(qiáng)度G=0.1,機(jī)械振子相干性整體向上平移了0.3數(shù)值,由圖2可以看出,聲子場方差壓縮因子ΔX2和相干性隨時間的演化呈現(xiàn)振蕩行為,t=0初始時刻,聲子場方差壓縮因子ΔX2為0.5,機(jī)械振子相干性為零,不存在壓縮,隨著金剛石NV色心和機(jī)械振子以及腔場的相互作用進(jìn)行,金剛石NV色心和機(jī)械振子以及腔場之間信息和能量不斷地交換和轉(zhuǎn)移,聲子場方差壓縮因子ΔX2小于0.5,聲子場出現(xiàn)壓縮,此時相干性C大于零呈現(xiàn)相干性,機(jī)械振子出現(xiàn)非經(jīng)典性,當(dāng)相互作用演化到t≈24時,聲子場方差壓縮因子ΔX2呈現(xiàn)最大壓縮振幅,此時出現(xiàn)最大方差壓縮,其物理機(jī)制為此時機(jī)械振子約化密度矩陣具有最大相干性(由圖2中虛線可知),機(jī)械振子非經(jīng)典性越強(qiáng),方差壓縮越大.

圖2 聲子場方差壓縮因子ΔX2和相干性隨時間的演化(參數(shù)g=1,θ1= π/3,G=0.1)Fig.2.The time evolution of the quadrature squeezing of the phonon field ΔX2and coherence.The parameters are g=1,θ1= π/3,G=0.1.

圖3給出了聲子場方差壓縮因子ΔX1和ΔX2隨時間的演化規(guī)律,圖3(a)和圖3(b)分別對應(yīng)聲子場方差壓縮因子ΔX1和ΔX2的演化特性,其中參數(shù)g=1,θ1=0,G=0.1.由(15)式可知當(dāng)θ1=0導(dǎo)致M2(t)=M5(t)=M6(t)=0,從而N1=N2=0,此時由(26)式可以看到聲子場方差壓縮因子ΔX1和ΔX2隨時間具有相同的演化規(guī)律.由圖3(a)和圖3(b)中可以看出,聲子場方差壓縮因子ΔX1和ΔX2隨時間演化取值都大于等于0.5,即聲子場在正交算符X1和X2方向上都不存在方差壓縮,不存在非經(jīng)典性,主要原因?yàn)榇藭r機(jī)械振子約化密度矩陣不存在相干性,機(jī)械振子處于混態(tài)ρm=ρ00|0〉〈0|+ρ11|1〉〈1|+ρ22|2〉〈2|,此時機(jī)械振子不存在非經(jīng)典特性,X1和X2方向不存在方差壓縮.

由于可以通過控制外加經(jīng)典場實(shí)現(xiàn)金剛石NV色心初始態(tài)的調(diào)控,圖4給出了聲子場方差壓縮因子ΔX2隨時間的演化動力學(xué)特性,圖中實(shí)線和虛線分別對應(yīng)金剛石NV色心不同的初態(tài)θ1=π/3和θ1=π/3.5,其他參數(shù)與圖2相同.由圖中實(shí)線可知,聲子場方差壓縮因子ΔX2初始時不存在方差壓縮,隨著金剛石NV色心和機(jī)械振子以及腔場的相互作用進(jìn)行,聲子場方差壓縮因子ΔX2小于0.5,聲子場出現(xiàn)壓縮,由圖4中插圖可知當(dāng)相互作用演化到t≈24時,此時出現(xiàn)最大方差壓縮,對比圖3中實(shí)現(xiàn)和虛線可以看到,隨著θ1的減小,聲子場方差壓縮越來越小,壓縮時間越來越短,由此結(jié)論通過調(diào)控金剛石NV色心初態(tài),可以制備壓縮時間長,壓縮幅度大的聲子場壓縮態(tài),實(shí)現(xiàn)對機(jī)械振子方差壓縮非經(jīng)典特性的操控,理論上提供了一種調(diào)控聲子場方差壓縮的方式.

圖3 聲子場方差壓縮因子ΔX1和ΔX2隨時間的演化(參數(shù)g=1,θ1=0,G=0.1)Fig.3.The time evolution of the quadrature squeezing of the phonon field ΔX1and ΔX2.The parameters are g=1,θ1=0,G=0.1.

為了探討金剛石NV色心與機(jī)械振子耦合強(qiáng)度對系統(tǒng)聲子場方差壓縮的影響,圖5給出了金剛石NV色心與機(jī)械振子耦合強(qiáng)度取不同值時聲子場方差壓縮因子ΔX2隨時間的演化行為,圖中實(shí)線和虛線分別對應(yīng)金剛石NV色心與機(jī)械振子耦合強(qiáng)度G=0.05和G=0.1,其他參數(shù)與圖2相同.由圖5可以看到,金剛石NV色心與機(jī)械振子耦合強(qiáng)度影響聲子場方差壓縮因子ΔX2的振蕩幅度以及壓縮大小.對比插圖中的實(shí)線和虛線可知,金剛石NV色心與機(jī)械振子耦合強(qiáng)度越大,正交算符X2上漲落振蕩的幅度越大,聲子場方差壓縮幅度越大,機(jī)械振子方差壓縮非經(jīng)典特性越強(qiáng).由金剛石NV色心與機(jī)械振子的耦合系數(shù)表達(dá)式G=gsμBGma0/?可知,實(shí)驗(yàn)上可以通過調(diào)控磁場梯度Gm實(shí)現(xiàn)金剛石NV色心與機(jī)械振子耦合強(qiáng)度的控制,從而實(shí)現(xiàn)對機(jī)械振子方差壓縮非經(jīng)典特性的操控.

圖4 聲子場方差壓縮因子ΔX2隨時間的演化(實(shí)線和虛線分別對應(yīng)θ1=π/3和θ1=π/3.5,其他參數(shù)與圖2相同)Fig.4.The time evolution of the quadrature squeezing of the phonon field ΔX2.The solid line and the dashed line correspond to θ1= π/3, θ1= π/3.5,respectively.The other parameters are the same as Fig.2.

考慮耗散對系統(tǒng)聲子場方差壓縮的影響,系統(tǒng)的主方程為

κ和γ分別為腔場衰減率和NV色心自發(fā)輻射衰減率.通過數(shù)值求解主方程,圖6給出了系統(tǒng)存在腔場衰減和NV色心的自發(fā)輻射對機(jī)械振子方差壓縮的影響,圖中實(shí)線、虛線和點(diǎn)虛線分別對應(yīng)參數(shù)κ=γ=0,κ=γ=0.001和κ=γ=0.005.對比插圖中的實(shí)線和虛線可以看到:考慮系統(tǒng)的腔場衰減和NV色心自發(fā)輻射耗散時,由于存在能量耗散過程,聲子場方差壓縮幅度越來越小,壓縮時間越來越短,機(jī)械振子方差壓縮非經(jīng)典特性越來越弱;隨著耗散強(qiáng)度的繼續(xù)增大,由插圖中的點(diǎn)虛線(對應(yīng)參數(shù)κ=γ=0.005)發(fā)現(xiàn),聲子場方差壓縮效應(yīng)消失,耗散破壞了壓縮效應(yīng).

圖5 聲子場方差壓縮因子ΔX2隨時間的演化(實(shí)線和虛線分別對應(yīng)G=0.05和G=0.1,其他參數(shù)與圖2相同)Fig.5.The time evolution of the quadrature squeezing of the phonon field ΔX2.The solid line and the dashed line correspond to G=0.05,G=0.1,respectively.The other parameters are the same as Fig.2.

圖6 聲子場方差壓縮因子ΔX2隨時間的演化(實(shí)線、虛線和點(diǎn)虛線分別對應(yīng)κ=γ=0,κ=γ=0.001和κ=γ=0.005,其他參數(shù)與圖2相同)Fig.6.The time evolution of the quadrature squeezing of the phonon field ΔX2.The solid line,dashed line and dotted-dashed line correspond to κ = γ =0,κ=γ=0.001 and κ=γ=0.005,respectively.The other parameters are the same as Fig.2.

這里簡單地給出該方案在實(shí)驗(yàn)上實(shí)現(xiàn)的可行性.文獻(xiàn)[21,22]利用光鑷技術(shù)在液體中實(shí)現(xiàn)了內(nèi)含NV色心的納米金剛石的囚禁,Neukirch等[23]在實(shí)驗(yàn)上實(shí)現(xiàn)了內(nèi)含NV色心的納米金剛石在自由空間中的光學(xué)囚禁,文獻(xiàn)[24]利用相似的技術(shù)實(shí)驗(yàn)上也實(shí)現(xiàn)了納米粒子在真空中的囚禁以及質(zhì)心運(yùn)動的探測.McCutcheon和Loncar[25]設(shè)計和分析了NV色心放在光子晶體微腔中心或腔表面位置均能實(shí)現(xiàn)NV色心與腔場的強(qiáng)耦合方案,隨后文獻(xiàn)[26]實(shí)驗(yàn)上實(shí)現(xiàn)了單個NV色心與光子晶體諧振腔的確定性耦合.文獻(xiàn)[27]研究了兩能級發(fā)射器與腔場耦合強(qiáng)度是發(fā)射器與機(jī)械振子耦合強(qiáng)度的20倍時系統(tǒng)中量子關(guān)聯(lián)特性,同時在NV色心下面放置一個磁針,可以產(chǎn)生很強(qiáng)的梯度磁場[28],實(shí)現(xiàn)NV色心電子自旋與納米金剛石振動的耦合,控制磁場梯度可以很方便地實(shí)現(xiàn)金剛石NV色心與機(jī)械振子的耦合強(qiáng)度大小的操控.因此,本文研究的金剛石NV色心同時耦合腔場和機(jī)械振子系統(tǒng)中聲子場的方差壓縮系統(tǒng)具備實(shí)現(xiàn)的可行性.

4 結(jié) 論

本文充分利用金剛石NV色心在室溫下具有毫秒量級的長相干時間以及容易被微波操控等優(yōu)勢,研究了驅(qū)動金剛石NV色心同時耦合腔場和機(jī)械振子系統(tǒng)中聲子場的方差壓縮動力學(xué)特性,引入機(jī)械振子的相干性闡明了聲子場方差壓縮的物理機(jī)制,討論了金剛石NV色心初態(tài)和NV色心與機(jī)械振子耦合強(qiáng)度對聲子場方差壓縮的影響,結(jié)果發(fā)現(xiàn)控制和操縱系統(tǒng)參數(shù)可以制備壓縮時間長、壓縮幅度大的聲子場壓縮態(tài),并且通過調(diào)控NV色心初態(tài)以及磁場梯度可以實(shí)現(xiàn)對機(jī)械振子方差壓縮非經(jīng)典特性的操控,從而在理論上提供了一種調(diào)控聲子場方差壓縮的方式.本文的研究結(jié)果將在量子信息處理以及精密測量等研究領(lǐng)域具有重要的意義和應(yīng)用前景.

感謝清華大學(xué)交叉信息研究院助理研究員尹璋琦博士在本文修改過程中提供的有益討論和幫助.

[1]Carr D W,Evoy S,Sekaric L,Craighead H G,Parpia J M 1999Appl.Phys.Lett.75 920

[2]Blick R H,Roukes M L,Wegscheider W,Bichler M 1998Phys.B:Condensed Matter249 784

[3]Caves C M,Thorne K S,Drever R W P,Sandberg V D,Zimmermann M 1980Rev.Mod.Phys.52 341

[4]Sekaric L,Parpia J M,Craighead H G,Feygelson T,Houston B H,Butler J E 2002Appl.Phys.Lett.81 4455

[5]Xiang Z L,Ashhab S,You J Q,Nori F 2013Rev.Mod.Phys.85 623

[6]Doherty M W,Manson N B,Delaney P,Jelezko F,Wrachtrup J,Hollenberg L C L 2013Phys.Rep.528 1

[7]Yin Z,Li T,Zhang X,Duan L M 2013Phys.Rev.A88 033614

[8]Zhao N,Yin Z Q 2014Phys.Rev.A90 042118

[9]Dolde F,Fedder H,Doherty M W,N?bauer T,Rempp F,Balasubramanian G 2011Nat.Phys.7 459

[10]Toyli D M,de las Casas C F,Christle D J,Dobrovitski V V,Awschalom D D 2013Proc.Natl.Acad.Sci.USA110 8417

[11]Kolkowitz S,Jayich A C,Unterreithmeier Q P,Bennett S D,Rabl P,Harris J G,Lukin M D 2012Science335 1603

[12]Ovartchaiyapong P,Lee K W,Myers B A,Jayich A C 2011Nat.Commun.5 4429

[13]Li P B,Xiang Z L,Rabl P,Nori F 2016Phys.Rev.Lett.117 015502

[14]Muschik C A,Moulieras S,Bachtold A,Koppens F H,Lewenstein M,Chang D E 2014Phys.Rev.Lett.112 223601

[15]Liu B Y,Cui W,Dai H Y,Chen X,Zhang M 2017Chin.Phys.B26 090303

[16]Liu B Y,Dai H Y,Chen X,Zhang M 2015Eur.Phys.J.D69 104

[17]Rabl P,Cappellaro P,Dutt M V G,Jiang L,Maze J R,Lukin M D 2009Phys.Rev.B79 041302

[18]Liu Y X,Sun C P,Nori F 2006Phys.Rev.A74 052321

[19]Walls D F,Milburn G J,Garrison J C 1994Quantum Optics(Berlin:Springer-Verlag)pp297–303

[20]Yu C S,Song H S 2009Phys.Rev.A80 022324

[21]Horowitz V R,Alemán B J,Christle D J,Cleland A N,Awschalom D D 2012Proc.Natl.Acad.Sci.USA109 13493

[22]Geiselmann M,Juan M L,Renger J,Say J M,Brown L J,de Abajo F J,Koppens F,Quidant R 2013Nat.Nanotechnol.8 175

[23]Neukirch L P,Gieseler J,Quidant R,Novotny L,Nick V A 2013Opt.Lett.38 2976

[24]Gieseler J,Deutsch B M,Quidant R,Novotny L 2012Phys.Rev.Lett.109 103603

[25]Mccutcheon M W,Loncar M 2008Opt.Express16 19136

[26]Englund D,Shields B,Rivoire K,Hatami F,Vu?kovi? J,Park H,Lukin M D 2010Nano Lett.10 3922

[27]Restrepo J,Favero I,Ciuti C 2017Phys.Rev.A95 023832

[28]Mamin H J,Poggio M,Degen C L,Rugar D 2007Nat.Nanotechnol.2 301