趙麗梅 張國鋒
(北京航空航天大學(xué)物理科學(xué)與核能工程學(xué)院,北京 100191)(2017年5月24日收到;2017年8月22日收到修改稿)
隨著量子信息科學(xué)的蓬勃發(fā)展,量子力學(xué)和熱力學(xué)的交叉領(lǐng)域逐漸引起了越來越多物理學(xué)家的關(guān)注[1,2].自從量子熱力學(xué)這一概念被Kieu首次提出以來,物理學(xué)者們在這一領(lǐng)域做了大量的研究工作,使得量子熱力學(xué)循環(huán)性能成為近年來研究的熱點(diǎn)[3,4].目前研究比較多的體系有自旋系統(tǒng)[5?27]、諧振子系統(tǒng)[28?31]二能級(jí)或多能級(jí)系統(tǒng)[32,33]以及腔量子電動(dòng)力學(xué)系統(tǒng)[34,35]等.
Dzyaloshinski-Mariya(DM)相互作用[36,37]是一種具有反對(duì)稱的各向異性超交換相互作用,可用來解釋反鐵磁材料中的弱鐵磁性以及其他的新奇性質(zhì)[38].近年來,人們也廣泛討論和研究了在各向異性交換作用下的海森伯模型,發(fā)現(xiàn)DM相互作用會(huì)對(duì)研究系統(tǒng)的性質(zhì)產(chǎn)生重要的影響.例如,Zhang[39]探討了具有DM相互作用的兩個(gè)自旋各向異性系統(tǒng)糾纏量子熱機(jī).Wang等[40]研究了具有DM相互作用的各向異性海森伯XXZ模型的熱糾纏等,這些研究工作都取得了豐富的成果.但是到目前為止,尚未見具體的DM相互作用對(duì)量子熱機(jī)影響的探究.
為了進(jìn)一步豐富量子熱力學(xué)循環(huán)理論,本文以帶有DM相互作用的兩比特模型作為工作物質(zhì),構(gòu)建一個(gè)四能級(jí)糾纏量子Otto熱機(jī)和量子Stirling熱機(jī),通過保持其他參數(shù)不變,僅改變DM相互作用參數(shù)來研究DM相互作用對(duì)量子熱機(jī)的影響.而這種操作在具體的實(shí)際過程中也很容易實(shí)現(xiàn)[41,42],例如,Niu等[43]研究了當(dāng)保持外磁場不變,只有粒子間的自旋耦合系數(shù)J發(fā)生變化的四能級(jí)量子Otto熱機(jī).
通過數(shù)值計(jì)算,探究了循環(huán)過程中在不同的DM相互作用參數(shù)下對(duì)熱機(jī)的基本熱力學(xué)參數(shù)的影響,并且利用等值線圖描繪和分析了糾纏與熱量、功和機(jī)械效率之間的關(guān)系.
量子Otto熱機(jī)是以具有DM相互作用各向異性的兩比特海森伯模型為工作物質(zhì),其哈密頓量可寫為
式中J表示自旋耦合系數(shù),反鐵磁性物質(zhì)對(duì)應(yīng)J>0,鐵磁性物質(zhì)對(duì)應(yīng)J<0,為了計(jì)算方便,本文只取J>0;B是沿Z方向的連續(xù)外磁場;D表示DM相互作用矢量,σi(i=1,2)是泡利算符.為了簡便,我們只考慮D=Dz時(shí)的情況,此時(shí)的哈密頓量可化簡為
選|00〉,|11〉,|10〉,|01〉為基矢,則HDM的本征態(tài)為
其中θ=arctanD,相應(yīng)的四個(gè)本征值為
當(dāng)系統(tǒng)處于熱平衡態(tài)時(shí),其狀態(tài)密度算符可寫為ρ=∑ipi|ψi〉〈ψi|. 其中,pi=exp(?βEi)/Z,配分函數(shù)Z=Tr[exp(?βEi)],β=1/(kBT).
圖1 (網(wǎng)刊彩色)當(dāng) T1=2,T2=1,J=1和B=4時(shí),(a)QH,(b)QL,(c)W 和(d)η隨D1和D2變化的等值線圖Fig.1.(color online)Variations of heat transferred(a)QH,(b)QL,(c)net work input W and(d)the efficiency η of the QHE with variables D1and D2in isoline map for parameters T1=2,T2=1,J=1 and B=4.
量子Otto熱機(jī)循環(huán)由兩個(gè)量子等容過程和兩個(gè)量子絕熱過程構(gòu)成,下面是量子Otto熱機(jī)循環(huán)的四個(gè)階段.
階段一:初始時(shí),體系的四個(gè)本征態(tài)的概率為pi,DM耦合相互作用為D=D1,外部磁場B保持不變,將體系與溫度為TH的高溫?zé)嵩唇佑|,經(jīng)過一段時(shí)間,達(dá)到熱力學(xué)平衡,在這個(gè)過程中系統(tǒng)本征能量Ei保持不變,體系的四個(gè)本征態(tài)的概率變?yōu)?在這個(gè)階段中,只有熱量的交換,不做功.
階段二:移除熱源,體系經(jīng)歷一個(gè)量子絕熱過程,DM耦合作用從D1變到D2,這個(gè)過程要足夠緩慢,使得每個(gè)量子態(tài)的概率不變,依舊為,體系各本征能級(jí)Ei變?yōu)?因此在這個(gè)過程中,只有做功而無熱量的傳遞,指出要使該過程保持絕熱,耦合作用的變化需要滿足傳統(tǒng)的量子絕熱定理[16].
階段三:與第一階段相似,體系與溫度為TL的低溫?zé)嵩闯掷m(xù)接觸,經(jīng)過一段時(shí)間的耦合作用后達(dá)到熱平衡狀態(tài),工作物質(zhì)的每個(gè)態(tài)的占有概率由變?yōu)閜i,體系各本征能級(jí)保持不變.在此過程中,由于占有概率的改變使體系與熱源之間存在熱傳遞,沒有做功.
圖2 (網(wǎng)刊彩色)當(dāng)T1=2,T2=1,J=1和B=6時(shí),(a)QH,(b)QL,(c)W 和(d)η隨D1和D2變化的等值線圖Fig.2.(color online)Variations of heat transferred(a)QH,(b)QL,(c)net work input W and(d)the efficiency η of the QHE with variables D1and D2in isoline map for parameters T1=2,T2=1,J=1 and B=6.
階段四:與第二階段類似,移除熱源,體系再次經(jīng)歷量子絕熱過程,這一過程,DM耦合作用從D2變到D1,每個(gè)本征態(tài)的占有概率pi均保持不變,體系各本征能級(jí)由又回到了Ei.此過程中體系只做功,沒有熱量的交換.
根據(jù)上述過程,可以得到在階段一和階段三的兩個(gè)等容過程中,轉(zhuǎn)移的能量分別為
式中,QH>0和QL<0分別代表系統(tǒng)吸熱和放熱.因此整個(gè)循環(huán)的輸出功為
量子熱機(jī)的機(jī)械效率為
當(dāng)T1=2,T2=1,J=1,磁場B分別為4和6時(shí),我們給出了以D1,D2為變量的函數(shù),通過數(shù)值分析,繪出了基本熱力學(xué)量QH,QL,W,η隨D1,D2變化的等值線圖(如圖1和圖2所示).
從圖1和圖2可以得到4個(gè)結(jié)論:1)DM相互作用不僅能夠提高系統(tǒng)的熱糾纏度,而且會(huì)對(duì)Otto熱機(jī)的基本量子熱力學(xué)量產(chǎn)生很大的影響,既能影響工作熱力學(xué)量的大小,也能影響其變化趨勢;2)當(dāng)系統(tǒng)做正功時(shí),滿足條件D1<D2,這個(gè)結(jié)論也可以從(7)式中推理得到,當(dāng)W>0,有化簡后得到B/T1<B/T2;3)當(dāng)W>0時(shí),QH>?QL>0始終成立,且在不同的外磁場參數(shù)下該糾纏熱機(jī)機(jī)械效率都小于卡諾熱機(jī)的效率,因此熱力學(xué)第二定律始終成立;4)隨著外磁場B的增大,做功和效率的取值范圍也會(huì)越來越大.
在本節(jié)中,考慮了同樣的耦合系統(tǒng)作為工作物質(zhì)構(gòu)建了量子Stirling循環(huán).量子Stirling循環(huán)是經(jīng)典Stirling循環(huán)的量子推廣,包括兩個(gè)量子等溫過程和兩個(gè)量子等容過程[44].該熱機(jī)的四個(gè)過程可以做如下表述.
階段一:是量子等溫過程.在該過程中,系統(tǒng)與處于TH的高溫?zé)嵩唇佑|.DM耦合作用緩慢的從D1變到D2,外部磁場B保持不變,并且在此過程中,系統(tǒng)與熱源在任何一個(gè)時(shí)刻都保持熱平衡.系統(tǒng)從高溫?zé)嵩次鼰釣镼1.
階段二:是量子等容過程.在該過程中,DM耦合作用為D2保持不變,系統(tǒng)的溫度從TH減小到TL,該過程系統(tǒng)不對(duì)外做功,但釋放熱量為Q23.
階段三:是另外一個(gè)量子等溫過程.在該過程中,系統(tǒng)與溫度為TL的低溫?zé)嵩唇佑|,DM耦合作用從D2變到D1.系統(tǒng)的溫度保持不變.系統(tǒng)向低溫?zé)嵩瘁尫艧崃?記為Q2).
階段四:是另一個(gè)量子等容過程.在該過程中,DM耦合作用為D1保持不變,系統(tǒng)吸收熱量(記為Q41),但不對(duì)外做功.整個(gè)系統(tǒng)循環(huán)過程如圖3所示.
圖3 (網(wǎng)刊彩色)量子Stirling熱機(jī)的D-U(DM作用:內(nèi)能)Fig.3.(color online)Scheme illustration of a quantum Stirling heat engine cycle based on DM interaction-energy D-U diagram.
在兩個(gè)等容過程中,由于系統(tǒng)對(duì)外不做功,因此其熱傳遞等于內(nèi)能U的改變.則該系統(tǒng)內(nèi)能
從(9)式可以得到系統(tǒng)與回?zé)崞髦g的熱交換Q23和Q41分別為
兩個(gè)等溫過程中,系統(tǒng)與外界環(huán)境之間的熱交換為dQ=TdS,這里的S是工作物質(zhì)的熵.則該系統(tǒng)熵為
代入到dQ=TdS中,得到在兩個(gè)等溫過程中系統(tǒng)從高溫?zé)嵩次盏臒崃恳约跋虻蜏責(zé)嵩瘁尫诺臒崃縌1和Q2分別為
與量子Otto循環(huán)很不相同,量子Stirling循環(huán)可以使用回?zé)崞?通過回?zé)崞鞅憧梢允寡h(huán)過程中的階段二釋放的熱量會(huì)在階段四中重新被系統(tǒng)吸收,從而達(dá)到改善熱機(jī)循環(huán)的性能的目的.典理想氣體為工作物質(zhì)的循環(huán)而言,系統(tǒng)達(dá)到最佳回?zé)岬臈l件是Q23+Q41=0.然而,以量子系統(tǒng)為工作物質(zhì)的循環(huán)卻并非如此.事實(shí)上,對(duì)于量子循環(huán)而言,系統(tǒng)與回?zé)崞髟趦蓚€(gè)等容過程中的凈熱傳遞量需要滿足:ΔQ=Q23+Q41,因此系統(tǒng)在整個(gè)循環(huán)過程中,從高溫?zé)嵩?溫度為T1)吸收的熱量和從低溫?zé)嵩?溫度為T2)釋放的熱量為
(13)式中,如果ΔQ<0,δ=0;ΔQ>0,δ=1;(14)式中,如果ΔQ>0,δ=0;ΔQ<0,δ=1.以這些熱交換為基礎(chǔ),根據(jù)熱力學(xué)第一定律可以得出系統(tǒng)對(duì)外做功為
W>0代表系統(tǒng)對(duì)外界做功.量子熱機(jī)的機(jī)械效率為
令T1=2,T2=1,J=1,B=4,通過數(shù)值軟件,就可以繪制出量子Stirling熱機(jī)在整個(gè)循環(huán)過程中ΔQ不同情況下的基本熱力學(xué)量QH,QL,W,η隨D1,D2的變化的等值線圖,結(jié)果如圖4–圖6所示.
從所示的等值線圖中可以得出3個(gè)結(jié)論:1)DM相互作用會(huì)對(duì)量子Stirling循環(huán)熱機(jī)的基本量子熱力學(xué)產(chǎn)生重要影響;2)當(dāng)W>0時(shí),QH>?QL>0,滿足條件D1<D2,這個(gè)結(jié)論也可以從(14)式中推理得到,最終化簡后得到B/T1<B/T2;3)在量子Stirling熱機(jī)循環(huán)中,熱機(jī)由于回?zé)崞鞯氖褂?其循環(huán)效率甚至?xí)^Carnot效率,熱力學(xué)第二定律在整個(gè)循環(huán)中都始終成立.
圖4 (網(wǎng)刊彩色)在Q23+Q41>0情況下,T1=2,T2=1,J=1和 B=4時(shí),(a)QH,(b)QL,(c)W 和 (d)η隨D1和D2變化的等值線圖Fig.4.(color online)In the case of Q23+Q41>0,variations of heat transferred(a)QH,(b)QL,(c)net work input W and(d)the efficiency η of the QHE with variables D1and D2in isoline map for parameters T1=2,T2=1,J=1 and B=4.
圖5 (網(wǎng)刊彩色)在Q23+Q41=0情況下,T1=2,T2=1,J=1和B=4時(shí),(a)QH,(b)QL,(c)W 和 (d)η隨D1和D2變化的等值線圖Fig.5.(color online)In the case of Q23+Q41=0,variations of heat transferred(a)QH,(b)QL,(c)net work input W and(d)the efficiency η of the QHE with variables D1and D2in isoline map for parameters T1=2,T2=1,J=1 and B=4.
以含有DM相互作用的兩比特模型為工作物質(zhì)構(gòu)建了一個(gè)糾纏量子熱機(jī).分別介紹了兩種不同的量子熱機(jī)模型:量子Otto熱機(jī)和量子Stirling熱機(jī).通過保持其他參量不變,只有DM相互作用發(fā)生改變,利用數(shù)值運(yùn)算,分析出了熱機(jī)循環(huán)中DM相互作用與熱傳遞、做功以及效率等熱力學(xué)量之間的關(guān)系.結(jié)果表明:DM相互作用對(duì)兩種熱機(jī)的基本熱力學(xué)量都會(huì)產(chǎn)生重要影響;發(fā)現(xiàn)量子Otto熱機(jī)的機(jī)械效率始終小于卡諾熱機(jī)的效率,但是量子Stirling熱機(jī)的循環(huán)效率大于量子Otto熱機(jī),甚至可以超過Carnot效率.在這兩個(gè)糾纏體系中,熱力學(xué)第二定律都依然成立.
[1]Scovil H E D,Schulz-Dubois E O 1959Phys.Rev.Lett.2 262
[2]Geusic J E,Schulz-Dubois E O,Scovil H E D 1967Phys.Rev.156 343
[3]Kieu T D 2004Phys.Rev.Lett.93 140403
[4]Kieu T D 2006Eur.Phys.J.D39 115
[5]Altintas F,Hardal A U C,Mustecaplioglu O E 2015Phys.Rev.A91 023816
[6]Wang X G 2001Phys.Rev.A64 012313
[7]Thomas G,Johal R S 2011Phys.Rev.E83 031135
[8]Huang X L,Wang L C,Yi X X 2013Phys.Rev.E87 012144
[9]Zhou Y,Zhang G F,Li S S 2009Europhys.Lett.86 50004
[10]Zhang G F 2007Phys.Rev.A75 034304
[11]Feldmann T,KosloffR 2004Phys.Rev.E70 046110
[12]Feldmann T,KosloffR 2003Phys.Rev.E68 016101
[13]KosloffR,Feldmann T 2002Phys.Rev.E65 055102
[14]Henrich M J,Mahler G,Michel M 2007Phys.Rev.E75 051118
[15]Zhang T,Liu W T,Chen P X,Li Z 2007Phys.Rev.A75 062102
[16]Thomas G,Johal R S 2014Eur.Phys.J.B87 166
[17]Huang X L,Wang T,Yi X X 2012Phys.Rev.E86 051105
[18]Huang X L,Liu Y,Wang Z,Niu X Y 2014Eur.Phys.J.Plus129 4
[19]Wu F,Chen L,Sun F,Wu C,Li Q 2006Phys.Rev.E73 016103
[20]Ivanchenko E A 2015Phys.Rev.E92 032124
[21]Altintas F,Müstecaplio?lu ? E 2015Phys.Rev.E92 022142
[22]He X,He J,Zheng J 2012Physica A391 6594
[23]Cakmak S,Altintas F,Müstecaplio?lu ? E 2016Eur.Phys.J.Plus131 197
[24]Wang H,Liu S,He J 2009Phys.Rev.E79 041113
[25]Hubner W,Lefkidis G,Dong C D,Chaudhuri D 2014Phys.Rev.B90 024401
[26]Azimi M,Chotorlishvili L,Mishra S K,Vekua T,Hubner W,Berakdar J 2014New J.Phys.16 063018
[27]Albayrak E 2013Int.J.Quantum.Inform.11 1350021
[28]Dillenschneider R,Lutz E 2009Europhys.Lett.88 50003
[29]Woo C H,Wen H,Semenov A A,Dudarev S L,Ma P W 2015Phys.Rev.B91 104306
[30]Ro?nagel J,Abah O,Schmidt-Kaler F,Singer K,Lutz E 2014Phys.Rev.Lett.112 030602
[31]Zhang X Y,Huang X L,Yi X X 2014J.Phys.A:Math.Theor.47 455002.
[32]Wang R,Wang J,He J,Ma Y 2013Phys.Rev.E87 042119
[33]Uzdin R,KosloffR 2014Europhys.Lett.108 40001
[34]Altintas F,Hardal A U C,Mustecaplioglu O E 2015Phys.Rev.A91 023816
[35]Quan H T,Zhang P,Sun C P 2006Phys.Rev.E73 036122
[36]Dzyaloshinskii I 1958J.Phys.Chem.Sol.4 241
[37]Moriya T 1960Phys.Rev.Lett.4 228
[38]Sun Q F,Xie X C,Wang J 2007Phys.Rev.Lett.98 196801
[39]Zhang G F 2008Eur.Phys.J.D49 123
[40]Li D C,Wang X P,Cao Z L 2008J.Phys.Condens.Matter20 325229
[41]Zhong X M,Nguyen B A,Yun J X 2016Phys.Rev.E94 042135
[42]RoSSnagel J,Dawkins S T,Tolazzi K N 2016Science352 325
[43]Niu X Y,Huang X L,Shang Y F,Wang X Y 2015Int.J.Mod.Phys.B29 1550086
[44]Huang X L,Niu X Y,Xiu X M,Yi X X 2014Eur.Phys.J.D68 32