光子驱动量子点制冷机*

李唯 符婧 杨贇贇 何济洲

(南昌大学物理系,南昌 330031)

提出了由两个二能级量子点、一个光子库与两个导体端构成的光子驱动量子点制冷机模型.基于主方程,导出了制冷机的制冷率和制冷系数的表达式,获得了制冷机处于紧耦合时所满足的条件.接着,数值模拟出该制冷机处于紧耦合和一般情况下制冷率与制冷系数之间的性能特征图,确定了制冷机性能的优化范围.最后,以最大制冷率、最大制冷率下的制冷系数、最大制冷系数和最大制冷系数下的制冷率作为优化目标,分析了光子库温度、跃迁系数和温比对制冷机性能的影响.

1 引 言

相对于传统的热机而言,热电器件具有结构轻便、体积小、无运动部件等优点,在工业和汽车尾气废热回收、太阳能电池开发、计算机芯片散热与制冷等方面有广泛的应用前景[1].但是,相对低的热电转换效率或输出功率又限制了热电器件的商业化应用.因此,人们一直在致力于寻找高热电优值的热电材料(高电传导、低热传导),如量子阱、超晶格、量子点、硅纳米线、碳纳米管、分子结等纳米热电材料[2-4].近年来,三端结构纳米热电效应和热电装置性能的研究,在理论[5-24]与实验[23-28]方面引起了研究人员的极大兴趣,成为纳米能量转换器件研究的热点.

例如Cleuren等[29]提出了一种利用光子直接进行制冷的光子驱动纳米制冷机.这种装置没有运动部件,也没有净电流,制冷是“冷”电子取代“热”电子的结果.然而,Amikam 等[30]认为 Cleuren 等人在处理量子点间的电子传输时,忽略了量子点之间的电子跃迁.Wang等[31]研究了一种光子驱动四能级制冷机装置,该装置由两个玻色库与一个四能级系统连接而成,由光子驱动进行制冷.基于前人的工作,我们考虑了单量子点内不同能级间的电子跃迁,提出了光子驱动量子点制冷机新模型,分析了该光子驱动纳米制冷机的热力学性能,并讨论了主要参数对制冷机的优化性能影响.

2 模型与理论分析

光子驱动量子点制冷机模型如图1所示.制冷机由两个导体端,一个光子库和两个量子点系统组成.两个导体库的温度分别是TL,TR与光子库的温度TS之间的关系满足TS＞TL＞TR,且两个导体库的化学势相同µl=µr=µ.左边量子点只与左导体库耦合,右边量子点只与右导体库耦合.每个量子点都具有较低和较高的两个能级,且两个量子点的高能级的能级差和低能级的能级差相同.在这个模型中,量子点之间的电子跃迁是由光子驱动的.在高温光子的影响下,“冷”电子通过两个量子点的较低能级从左导体库向右导体库传输,同时“热”电子通过两个量子点的较高能级向相反的方向传输.对右导体库来说,“冷”电子替换了“热”电子,对其起到制冷的效果.这里的“冷”、“热”电子分别代表着能级低于费米能级(或化学势)εF=µ,和能级高于费米能级εF的电子.同时,考虑单个量子点不同能级间电子跃迁是由声子驱动的.左量子点上电子从高能级向低能级跃迁放出热量给左库,同时右量子点上电子从低能级向高能级跃迁从右库吸收热量,同样对右库有制冷的效果.图中黑色双向箭头表示电子的隧穿,棕色波浪线表示由光子驱动引起的电子跃迁,红色波浪线表示由声子驱动引起的电子跃迁.

图1 光子驱动量子点制冷机模型图Fig.1.A model of a quantum dot refrigerator driven by photon.

在弱耦合的情况下(Γ ≪kBT )电子的传输过程可以描述为序贯隧穿,系统量子态的占有几率随时间的变化可以由主方程表示,即P=MPi(i=0,l1,l2,r1,r2),其中

ki←j(i,j=0,l1,l2,r1,r2)为从j态到i态的传输几率.根据费米黄金定则,传输几率分别为:

其中,比例常数 Γi(i=l1,l2,r1,r2) 是导体库和量子点间的耦合系数,ki←0表示电子由左(右)导体库传输到左(右)量子点的几率,而 k0←j表示电子由左(右)量子点传输到左(右)导体库的几率.

其中,ΓS是两个量子点之间的耦合系数,而 Γnr是单个量子点内电子在两个能级间的跃迁系数;(6)和(7)式表示的是电子在两个量子点之间的传输几率,而(8)和(9)式表示的是电子在单个量子点内两个能级之间的传输几率;εg=εl2- εr2=εr1- εl1为能级差.

从左量子点低能级向右量子点低能级跃迁的电子流为

从右量子点高能级向左量子点高能级跃迁的电子流为

从左量子点内高能级向低能级跃迁的电子流为

从右量子点内低能级向高能级跃迁的电子流为

J1=Jl1+Jl2表示从左导体库流出的净电子流,J2=Jr1+Jr2表示流入右导体库的净电子流,系统遵守电子流守恒,则 J1=J2.电子流JS1=Jl1+Jnr1分为两部分,一部分是来自左导体库的电子流,另一部分是左量子点高能级向低能级跃迁的电子流.同理电子流之间的关系为JS1=Jl1+Jnr1=

利用上述稳态时电子流的表达式,可以获得电子在传输过程中引起的热量交换.从左、右导体库流出的热流和从光子库吸收的热流的具体表达式为:

热力学力为 FR=1/TL-1/TR,FS=1/TL-1/TS,且其范围满足 FR≤0,FS＞0.利用 (21)和 (22)式可以得到制冷系数的表达式为

可以理解为: 式中第一部分是热机在温度为 TL和TS的热源之间工作时,热机的卡诺效率被熵产生率相关因子降低,这个因子是对所涉及过程的不可逆性的衡量;第二部分是工作在温度为 TL和TR的热源之间的制冷机在可逆工作状态下的制冷系数.当熵产生率为零=0,制冷机达到可逆工作状态,其制冷系数为可逆值

当两个热力学力 FR,FS同时为零时,熵产生率为零=0,制冷机达到可逆状态,但这个条件意味着系统达到完全平衡 TL=TR=TS,失去了物理意义.另外,通过对参数进行适当的调整,使得制冷机内各个热流之间互成比例关系,也就是紧耦合条件[32-34].经分析可得以下两种情况,使制冷机处于紧耦合条件下:

1) 当化学势满足 µ=(εr1+εr2)/2 时,根据各个电子流之间的关系,(18)—(20)式可以化简为

2)调节参数使得 J1=J2=0,此时其他电子流之间还存在的关系,因此同样成比例关系.

在这两个条件下制冷系数都为

只与能级(εg,εr1,εr2)有关.熵产生率可以表达为当 η FR+FS=0 时,熵产生率达到=0.因此,只要模型参数满足关系:制冷机工作在可逆状态.而紧耦合条件下模型参数的取值范围为

3 性能特征分析

利用 (26)和 (28)式,通过数值模拟,可以分析制冷机的性能特征.图2和图3分别是紧耦合条件下在不同光子库温度TS和不同跃迁系数Γnr下,改变能级 εr2得到的制冷率和制冷系数之间的关系特征曲线图.做图中,紧耦合条件利用的是µ=(εr1+εr2)/2 的 情况.为方便起见,我们取µ/kB=0K 0 K,则 εr1=-εr2.其他参数的选取为 εg/kB=2K,Γl1=Γl2=Γr1=Γr2=ΓS=Γ,TL=8K,TR=4K.能级 εr2的取值范围可以利用条件=0 来确定.从图形来看,特征曲线是类似抛物线的开型曲线,制冷机的性能系数可以达到可逆值 ηr.温度 TS增大时,最大制冷率和最大制冷系数也增大;增大跃迁 系数 Γnr,最大制 冷率也增大.一般 地,利 用(19)和(21)式,可以做出不同光子库温度 TS和不同跃迁系数 Γnr下制冷率与制冷系数的特征曲线,即图4和5.同样改变能级 εr2的值,且利用制冷率为零的条件确定它的范围,化学势取 µ/kB=0K,能级 εr1/kB=-1K,其他参数不变.从图形可以看到,特征曲线是封闭的闭合曲线,起点和终点都在原点上,制冷机不能达到可逆工作状态.跃迁系数Γnr越小,最大制冷率下的制冷系数 ηQR与最大制冷系数 ηmax的值越接近,即可以在得到最大制冷率的同时得到较大的制冷系数.最大制冷率下的制冷系数 ηQR,和最大制冷系数 ηmax下的制冷率都可以通过数值模拟的方法计算.作为制冷机应该尽可能获得较大的制冷率和制冷系数,因此,制冷机的优化范围为:

图2 紧耦合条件下制冷率与制冷系数在不同温度 TS 下的关系Fig.2.The relation curves of the cooling rate and the coefficient of performance at different temperature TS under the condition of tight coupling.

图3 紧耦合条件下制冷率与制冷系数在不同跃迁系数Γnr下的关系Fig.3.The relation curves of the cooling rate and the coefficient of performance at different transition coefficientΓnr under the condition of tight coupling.

图4 一般情况下制冷率与制冷系数在不同温度 TS 下的关系Fig.4.The relation curves of the cooling rate and the coefficient of performance at different temperature TS in the general case.

图5 一般情况下制冷率与制冷系数在不同跃迁系数Γnr下的关系Fig.5.The relation curves of the cooling rate and the coefficient of performanceat different transition coefficientΓnr in the general case.

4 优化性能分析

利用(19)—(21)式和极值条件

图6 在不同温度 TS 下,两个优化性能参数 和 ηQR随温比的变化Fig.6.The curves of two optimal performance parameters and ηQR changing with the temperature ratio at different temperature TS.

图7 在不同跃迁系数 Γnr 下,两个优化性能参数和ηQR随着温比的变化Fig.7.The curves of two optimal performance parameters and ηQR changing with the temperature ratio at different transition coefficient Γnr.

图8 在不同温度 TS 下,两个优化性能参数 和 ηQR随温比的变化Fig.8.The curves of two optimalperformance parametersand ηQR changing with the temperature ratio at different temperature TS.

图9 在不同跃迁系数 Γnr 下,两个优化性能参数和ηQR随着温比的变化Fig.9.The curves of two optimal performance parameters and ηQR changing with the temperature ratio at different transition coefficient Γnr.

图10 在不同温度TS下,两个优化性能参数ηmax和随温比的变化Fig.10.The curves of two optimal performance parameters ηmaxand changing with the temperature ratio at different temperature TS.

5 结 论

本文主要研究结果如下:

1)通过分析,达到紧耦合的情况有两种,一是使化学势µ与能级εr2,εr1满足µ=(εr1+εr2)/2,但其净电子流不为零;二是使净电子流Jl1+Jl2为零;

2)在紧耦合条件下,制冷率和制冷系数之间的性能特征图为开型,制冷机可以达到可逆状态;最大制冷率和对应的制冷系数都随光子库温度的增大而增大;最大制冷率随跃迁系数的增大而增大,而对应的制冷系数随跃迁系数的增大而减小;因此合理选择光子库温度、跃迁系数和温比可以提高制冷机的性能;

3)在一般情况下,制冷率和制冷系数之间的性能特征图为闭型,制冷机不能达到可逆状态;优化性能参数有类似于紧耦合条件下的结论.