搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

光子驱动量子点制冷机

李唯 符婧 杨贇贇 何济洲

引用本文:
Citation:

光子驱动量子点制冷机

李唯, 符婧, 杨贇贇, 何济洲

Quantum dot refrigerator driven by photon

Li Wei, Fu Jing, Yang Yun-Yun, He Ji-Zhou
PDF
HTML
导出引用
  • 提出了由两个二能级量子点、一个光子库与两个导体端构成的光子驱动量子点制冷机模型. 基于主方程, 导出了制冷机的制冷率和制冷系数的表达式, 获得了制冷机处于紧耦合时所满足的条件. 接着, 数值模拟出该制冷机处于紧耦合和一般情况下制冷率与制冷系数之间的性能特征图, 确定了制冷机性能的优化范围. 最后, 以最大制冷率、最大制冷率下的制冷系数、最大制冷系数和最大制冷系数下的制冷率作为优化目标, 分析了光子库温度、跃迁系数和温比对制冷机性能的影响.
    A model of quantum dot refrigerator driven by photon, which consists of two two-level quantum dots, a photon reservoir and two leads, is proposed in this paper. Comparing with previous studies, we consider the transitions of electrons between different energy levels in a single quantum dot, which is more practical.Based on the theory of master equation and the assumption of weak coupling, we derive the expression of the cooling rate and the coefficient of performance of the refrigerator and obtain the condition of the tight coupling of the refrigerator operation. Next, we plot numerically the performance characteristic curves between the cooling rate and the coefficient of performance in the case of the tight coupling and in the general case. We find that the curves between the cooling rate and the coefficient of performance are opened loops for tight coupling, but they are closed loops in the general case. And we gain the conclusions that the refrigerator can be reversible under the condition of the tight coupling, while it can be irreversible in the general case. Then the optimally operating range of the refrigerator is determined. Finally, the effect of the temperature of the photon reservoir, transition coefficient, and temperature ratio on the performance of refrigerator under the conditions of the maximum cooling rate are studied, and also the coefficient of performance under the maximum cooling rate, the maximum coefficient of performanceand the cooling rate under the maximum coefficient of performanceare analyzed in detail.
      通信作者: 何济洲, hjzhou@ncu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11365015)资助的课题
      Corresponding author: He Ji-Zhou, hjzhou@ncu.edu.cn
    • Funds: Supported by the National Natural Science Foundations of China (Grant No. 11365015)
    [1]

    Giazotto F, Heikkila T T, Luukanen A, Savin A M, Pekola J P 2006 Rev. Mod. Phys. 78 217Google Scholar

    [2]

    Koumoto K, Mori T 2013 Thermoelectric Nanomaterials: Materials Design and Applications (Vol. 182) (Berlin: Springer Press) pp255-285

    [3]

    Maciá E 2015 Thermoelectric Materials: Advances and Applications (Jenny: Stanford Publishing)

    [4]

    Pichanusakorn P, Bandaru P 2010 Mater. Sci. Eng. R 67 19Google Scholar

    [5]

    Benenti G, Casati G, Saito K, Whitney R S 2017 Phys. Reports 694 1Google Scholar

    [6]

    Sothmann B, Sánchez R, Jordan A N 2014 Nanotechnology 26 032001Google Scholar

    [7]

    Sánchez R, Büttiker M 2011 Phys. Rev. B 83 085428Google Scholar

    [8]

    Dare A M, Lombardo P 2017 Phys. Rev. B 96 115414Google Scholar

    [9]

    Zhang Y C, Gin G X, Chen J C 2015 Phys. Rev. E 91 052118Google Scholar

    [10]

    Sothmann B, Sánchez R, Jordan A N, Büttiker M 2012 Phys. Rev. B 85 205301Google Scholar

    [11]

    Zhang Y C, Zhang X, Ye Z L, Gin G X, Chen J C 2017 Appl. Phys. Lett. 110 153501Google Scholar

    [12]

    Jordan A N, Sothmann B, Sánchez R, Buttiker M 2013 Phys. Rev. B 87 075312Google Scholar

    [13]

    Prance J R, Smith C G, Griffiths J P, Chorley S J, Anderson D, Jones G A C, Farrer I, Ritchie D A 2009 Phys. Rev. Lett. 102 146602Google Scholar

    [14]

    Sothmann B, Sánchez R, Jordan A N, Büttiker M 2013 N. J. Phys. 15 095021Google Scholar

    [15]

    Choi Y J, Jordan A N 2015 Physica E 74 465Google Scholar

    [16]

    Su S H, Zhang Y C, Chen J C, Shih T M 2016 Sci. Rep. 6 21425Google Scholar

    [17]

    Wohlman O E, Imry Y, Aharony A 2015 Phys. Rev. B 91 054302Google Scholar

    [18]

    Szukiewicz B, Eckern U, Wysokinski K I 2016 New J. Phys. 18 023050Google Scholar

    [19]

    Whitney R S, Sánchez R, Haupt F, Splettstoesser J 2016 Physica E 75 257Google Scholar

    [20]

    Lim J S, Sanchez D, Lopez R 2018 New J. Phys. 20 023038Google Scholar

    [21]

    Walldorf N, Jauho A P, Kaasbjerg K 2017 Phys. Rev. B 96 115415

    [22]

    Jiang J H, Entin-Wohlman O, Imry Y 2013 New J. Phys. 15 075021Google Scholar

    [23]

    Su H, Shi Z C, He J Z 2015 Chin. Phys. Lett. 32 100501Google Scholar

    [24]

    Lin Z B, Li W, Fu J, Yang Y Y, He J Z 2019 Chin. Phys. Lett. 36 060501Google Scholar

    [25]

    Josefsson M, Svilans A, Burke A M, Hoffmann E A, Fahlvik S, Thelander C, Leijnse M, Linke H 2018 Nature Nanotechnol. 13 920Google Scholar

    [26]

    Thierschmann H, Sanchez R, Sothmann B, Arnold F, Heyn C, Hansen W, Buhmann H, Molenkamp L W 2015 Nature Nanotechnol. 10 854Google Scholar

    [27]

    Roche B, Roulleau P, Ulien T J, Jompol Y, Farrer I, Ritchie D A, Glattli D C 2015 Nature Commun. 6 6738Google Scholar

    [28]

    Hartmann F, Pfeffer P, Hofling S, Kamp M, Worschech L 2015 Phys. Rev. Lett. 114 146805Google Scholar

    [29]

    Cleuren B, Rutten B, van den Broeck C 2012 Phys. Rev. Lett. 108 120603Google Scholar

    [30]

    Levy A, Alicki R, Kosloff R 2012 Phys. Rev. Lett. 109 248901Google Scholar

    [31]

    Wang J H, Lai Y M, Ye Z L, He J Z, Ma Y L, Liao Q H 2015 Phys. Rev. E 91 050102Google Scholar

    [32]

    van den Broeck C 2005 Phys. Rev. Lett. 95 190602Google Scholar

    [33]

    Yuan Y, Wang R, He J Z, Ma Y L, Wang J H 2014 Phys. Rev. E 90 052151Google Scholar

    [34]

    Sheng S Q, Tu Z C 2014 Phys. Rev. E 89 012129Google Scholar

  • 图 1  光子驱动量子点制冷机模型图

    Fig. 1.  A model of a quantum dot refrigerator driven by photon.

    图 2  紧耦合条件下制冷率与制冷系数在不同温度${T_{\rm{S}}}$下的关系

    Fig. 2.  The relation curves of the cooling rate and the coefficient of performance at different temperature ${T_{\rm{S}}}$ under the condition of tight coupling.

    图 3  紧耦合条件下制冷率与制冷系数在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下的关系

    Fig. 3.  The relation curves of the cooling rate and the coefficient of performance at different transition coefficient ${\varGamma _{n{\rm{r}}}}$ under the condition of tight coupling.

    图 4  一般情况下制冷率与制冷系数在不同温度${T_{\rm{S}}}$下的关系

    Fig. 4.  The relation curves of the cooling rate and the coefficient of performance at different temperature ${T_{\rm{S}}}$ in the general case.

    图 5  一般情况下制冷率与制冷系数在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下的关系

    Fig. 5.  The relation curves of the cooling rate and the coefficient of performanceat different transition coefficient ${\varGamma _{n{\rm{r}}}}$ in the general case.

    图 6  在不同温度${T_{\rm{S}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随温比的变化

    Fig. 6.  The curves of two optimal performance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different temperature ${T_{\rm{S}}}$

    图 7  在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随着温比的变化

    Fig. 7.  The curves of two optimal performance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different transition coefficient ${\varGamma _{n{\rm{r}}}}$

    图 8  在不同温度${T_{\rm{S}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随温比的变化

    Fig. 8.  The curves of two optimalperformance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different temperature ${T_{\rm{S}}}$.

    图 10  在不同温度${T_{\rm{S}}}$下, 两个优化性能参数${\eta ^{\max }}$${\dot Q}_{\rm{R}}^\eta $随温比的变化

    Fig. 10.  The curves of two optimal performance parameters ${\eta ^{\max }}$ and ${\dot Q}_{\rm{R}}^\eta $ changing with the temperature ratio at different temperature ${T_{\rm{S}}}$.

    图 9  在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随着温比的变化

    Fig. 9.  The curves of two optimal performance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different transition coefficient ${\varGamma _{n{\rm{r}}}}$.

  • [1]

    Giazotto F, Heikkila T T, Luukanen A, Savin A M, Pekola J P 2006 Rev. Mod. Phys. 78 217Google Scholar

    [2]

    Koumoto K, Mori T 2013 Thermoelectric Nanomaterials: Materials Design and Applications (Vol. 182) (Berlin: Springer Press) pp255-285

    [3]

    Maciá E 2015 Thermoelectric Materials: Advances and Applications (Jenny: Stanford Publishing)

    [4]

    Pichanusakorn P, Bandaru P 2010 Mater. Sci. Eng. R 67 19Google Scholar

    [5]

    Benenti G, Casati G, Saito K, Whitney R S 2017 Phys. Reports 694 1Google Scholar

    [6]

    Sothmann B, Sánchez R, Jordan A N 2014 Nanotechnology 26 032001Google Scholar

    [7]

    Sánchez R, Büttiker M 2011 Phys. Rev. B 83 085428Google Scholar

    [8]

    Dare A M, Lombardo P 2017 Phys. Rev. B 96 115414Google Scholar

    [9]

    Zhang Y C, Gin G X, Chen J C 2015 Phys. Rev. E 91 052118Google Scholar

    [10]

    Sothmann B, Sánchez R, Jordan A N, Büttiker M 2012 Phys. Rev. B 85 205301Google Scholar

    [11]

    Zhang Y C, Zhang X, Ye Z L, Gin G X, Chen J C 2017 Appl. Phys. Lett. 110 153501Google Scholar

    [12]

    Jordan A N, Sothmann B, Sánchez R, Buttiker M 2013 Phys. Rev. B 87 075312Google Scholar

    [13]

    Prance J R, Smith C G, Griffiths J P, Chorley S J, Anderson D, Jones G A C, Farrer I, Ritchie D A 2009 Phys. Rev. Lett. 102 146602Google Scholar

    [14]

    Sothmann B, Sánchez R, Jordan A N, Büttiker M 2013 N. J. Phys. 15 095021Google Scholar

    [15]

    Choi Y J, Jordan A N 2015 Physica E 74 465Google Scholar

    [16]

    Su S H, Zhang Y C, Chen J C, Shih T M 2016 Sci. Rep. 6 21425Google Scholar

    [17]

    Wohlman O E, Imry Y, Aharony A 2015 Phys. Rev. B 91 054302Google Scholar

    [18]

    Szukiewicz B, Eckern U, Wysokinski K I 2016 New J. Phys. 18 023050Google Scholar

    [19]

    Whitney R S, Sánchez R, Haupt F, Splettstoesser J 2016 Physica E 75 257Google Scholar

    [20]

    Lim J S, Sanchez D, Lopez R 2018 New J. Phys. 20 023038Google Scholar

    [21]

    Walldorf N, Jauho A P, Kaasbjerg K 2017 Phys. Rev. B 96 115415

    [22]

    Jiang J H, Entin-Wohlman O, Imry Y 2013 New J. Phys. 15 075021Google Scholar

    [23]

    Su H, Shi Z C, He J Z 2015 Chin. Phys. Lett. 32 100501Google Scholar

    [24]

    Lin Z B, Li W, Fu J, Yang Y Y, He J Z 2019 Chin. Phys. Lett. 36 060501Google Scholar

    [25]

    Josefsson M, Svilans A, Burke A M, Hoffmann E A, Fahlvik S, Thelander C, Leijnse M, Linke H 2018 Nature Nanotechnol. 13 920Google Scholar

    [26]

    Thierschmann H, Sanchez R, Sothmann B, Arnold F, Heyn C, Hansen W, Buhmann H, Molenkamp L W 2015 Nature Nanotechnol. 10 854Google Scholar

    [27]

    Roche B, Roulleau P, Ulien T J, Jompol Y, Farrer I, Ritchie D A, Glattli D C 2015 Nature Commun. 6 6738Google Scholar

    [28]

    Hartmann F, Pfeffer P, Hofling S, Kamp M, Worschech L 2015 Phys. Rev. Lett. 114 146805Google Scholar

    [29]

    Cleuren B, Rutten B, van den Broeck C 2012 Phys. Rev. Lett. 108 120603Google Scholar

    [30]

    Levy A, Alicki R, Kosloff R 2012 Phys. Rev. Lett. 109 248901Google Scholar

    [31]

    Wang J H, Lai Y M, Ye Z L, He J Z, Ma Y L, Liao Q H 2015 Phys. Rev. E 91 050102Google Scholar

    [32]

    van den Broeck C 2005 Phys. Rev. Lett. 95 190602Google Scholar

    [33]

    Yuan Y, Wang R, He J Z, Ma Y L, Wang J H 2014 Phys. Rev. E 90 052151Google Scholar

    [34]

    Sheng S Q, Tu Z C 2014 Phys. Rev. E 89 012129Google Scholar

  • [1] 周亮亮, 吴宏博, 李学铭, 唐利斌, 郭伟, 梁晶. ZrS2量子点: 制备、结构及光学特性. 物理学报, 2019, 68(14): 148501. doi: 10.7498/aps.68.20190680
    [2] 廖天军, 林比宏, 王宇珲. 新型高效热离子功率器件的性能特性研究. 物理学报, 2019, 68(18): 187901. doi: 10.7498/aps.68.20190882
    [3] 张强强, 胡建勇, 景明勇, 李斌, 秦成兵, 李耀, 肖连团, 贾锁堂. 单光子调制频谱用于量子点荧光寿命动力学的研究. 物理学报, 2019, 68(1): 017803. doi: 10.7498/aps.68.20181797
    [4] 张荣, 卢灿灿, 李倩文, 刘伟, 白龙. 线性不可逆热力学框架下一个无限尺寸热源而有限尺寸冷源的制冷机的性能分析. 物理学报, 2018, 67(4): 040502. doi: 10.7498/aps.67.20172010
    [5] 吴海娜, 孙雪, 公卫江, 易光宇. 电子-声子相互作用对平行双量子点体系热电效应的影响. 物理学报, 2015, 64(7): 077301. doi: 10.7498/aps.64.077301
    [6] 刘志民, 赵谡玲, 徐征, 高松, 杨一帆. 红光量子点掺杂PVK体系的发光特性研究. 物理学报, 2014, 63(9): 097302. doi: 10.7498/aps.63.097302
    [7] 张盼君, 孙慧卿, 郭志友, 王度阳, 谢晓宇, 蔡金鑫, 郑欢, 谢楠, 杨斌. 含有量子点的双波长LED的光谱调控. 物理学报, 2013, 62(11): 117304. doi: 10.7498/aps.62.117304
    [8] 屈俊荣, 郑建邦, 王春锋, 吴广荣, 郝娟. 聚对苯乙炔MOPPV/ZnSe量子点复合材料太阳电池性能研究. 物理学报, 2013, 62(7): 078802. doi: 10.7498/aps.62.078802
    [9] 屈俊荣, 郑建邦, 王春锋, 吴广荣, 王雪艳. 碳纳米管掺杂对聚合物聚(2-甲氧基-5-辛氧基)对苯乙炔-PbSe量子点复合材料性能的影响. 物理学报, 2013, 62(12): 128801. doi: 10.7498/aps.62.128801
    [10] 李霞, 冯东海, 潘贤群, 贾天卿, 单璐繁, 邓莉, 孙真荣. 室温下CdSe胶体量子点超快自旋动力学. 物理学报, 2012, 61(20): 207202. doi: 10.7498/aps.61.207202
    [11] 李霞, 冯东海, 何红燕, 贾天卿, 单璐繁, 孙真荣, 徐至展. CdTe/CdS核壳结构量子点超快载流子动力学 . 物理学报, 2012, 61(19): 197801. doi: 10.7498/aps.61.197801
    [12] 姜冰一, 郑建邦, 王春锋, 郝娟, 曹崇德. 基于GaAs/InAs-GaAs/ZnSe量子点太阳电池结构的优化. 物理学报, 2012, 61(13): 138801. doi: 10.7498/aps.61.138801
    [13] 何弦, 何济洲, 肖宇玲. 四能级量子制冷循环. 物理学报, 2012, 61(15): 150302. doi: 10.7498/aps.61.150302
    [14] 琚鑫, 郭健宏. 点间耦合强度对三耦合量子点系统微分电导的影响. 物理学报, 2011, 60(5): 057302. doi: 10.7498/aps.60.057302
    [15] 陈翔, 米贤武. 量子点腔系统中抽运诱导受激辐射与非谐振腔量子电动力学特性的研究. 物理学报, 2011, 60(4): 044202. doi: 10.7498/aps.60.044202
    [16] 王秀梅, 何济洲, 何弦, 肖宇玲. 非线性二极管系统构成的不可逆热机性能特征分析. 物理学报, 2010, 59(7): 4460-4465. doi: 10.7498/aps.59.4460
    [17] 彭红玲, 韩 勤, 杨晓红, 牛智川. 1.3μm量子点垂直腔面发射激光器高频响应的优化设计. 物理学报, 2007, 56(2): 863-870. doi: 10.7498/aps.56.863
    [18] 邓宇翔, 颜晓红, 唐娜斯. 量子点环的电子输运研究. 物理学报, 2006, 55(4): 2027-2032. doi: 10.7498/aps.55.2027
    [19] 侯春风, 郭汝海. 椭圆柱形量子点的能级结构. 物理学报, 2005, 54(5): 1972-1976. doi: 10.7498/aps.54.1972
    [20] 袁晓利, 施 毅, 杨红官, 卜惠明, 吴 军, 赵 波, 张 荣, 郑有钭. 硅量子点中电子的荷电动力学特征. 物理学报, 2000, 49(10): 2037-2040. doi: 10.7498/aps.49.2037
计量
  • 文章访问数:  3343
  • PDF下载量:  38
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-16
  • 修回日期:  2019-08-29
  • 上网日期:  2019-11-01
  • 刊出日期:  2019-11-20

光子驱动量子点制冷机

    基金项目: 国家自然科学基金(批准号: 11365015)资助的课题

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

English Abstract

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

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

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

      图  1  光子驱动量子点制冷机模型图

      Figure 1.  A model of a quantum dot refrigerator driven by photon.

      在弱耦合的情况下($\varGamma \ll {k_{\rm{B}}}T$)电子的传输过程可以描述为序贯隧穿, 系统量子态的占有几率随时间的变化可以由主方程表示, 即${\dot{ P}} = {{M}}{{{P}}_{{i}}}$$(i = 0, \;{\rm{l}}1, \;{\rm{l2}}, \;{\rm{r1}}, \;{\rm{r2}})$, 其中

      $\begin{split} & {{M}} = \left( {\begin{array}{*{20}{c}} {\mathcal M}_{11} & {{k_{0 \leftarrow {\rm{l2}}}}} & {{k_{0 \leftarrow {\rm{l}}1}}} & {{k_{0 \leftarrow {\rm{r}}2}}} & {{k_{0 \leftarrow {\rm{r1}}}}} \\ {{k_{{\rm{l2}} \leftarrow 0}}} & {\mathcal M}_{22} & {{k_{{\rm{l}}2 \leftarrow {\rm{l}}1}}} & {{k_{{\rm{l}}2 \leftarrow {\rm{r}}2}}} & 0 \\ {{k_{{\rm{l1}} \leftarrow 0}}}&{{k_{{\rm{l1}} \leftarrow {\rm{l}}2}}}& {\mathcal M}_{33} & 0 & {{k_{{\rm{l1}} \leftarrow {\rm{r1}}}}} \\ {{k_{{\rm{r2}} \leftarrow 0}}} & {{k_{{\rm{r2}} \leftarrow {\rm{l2}}}}} & 0 & {\mathcal M}_{44} & {{k_{{\rm{r2}} \leftarrow {\rm{r1}}}}} \\ {{k_{{\rm{r}}1 \leftarrow 0}}} & 0 & {{k_{{\rm{r}}1 \leftarrow {\rm{l1}}}}} & {{k_{{\rm{r1}} \leftarrow {\rm{r2}}}}} & {\mathcal M}_{55} \end{array}} \right), \\ & {\mathcal M}_{11}={{ - {k_{{\rm{l}}2 \leftarrow 0}} - {k_{{\rm{l}}1 \leftarrow 0}} - {k_{{\rm{r}}2 \leftarrow 0}} - {k_{{\rm{r}}1 \leftarrow 0}}}}, \\ & {\mathcal M}_{22}={{ - {k_{0 \leftarrow {\rm{l}}2}} - {k_{{\rm{l1}} \leftarrow {\rm{l}}2}} - {k_{r2 \leftarrow l2}}}}, \\ &{\mathcal M}_{33}={{ - {k_{0 \leftarrow {\rm{l}}1}} - {k_{{\rm{l}}2 \leftarrow {\rm{l}}1}} - {k_{{\rm{r1}} \leftarrow {\rm{l}}1}}}}, \\ & {\mathcal M}_{44}={{ - {k_{0 \leftarrow {\rm{r}}2}} - {k_{{\rm{l}}2 \leftarrow {\rm{r}}2}} - {k_{{\rm{r1}} \leftarrow {\rm{r}}2}}}}, \\ & {\mathcal M}_{55}={{ - {k_{0 \leftarrow {\rm{r1}}}} - {k_{{\rm{l1}} \leftarrow {\rm{r1}}}} - {k_{{\rm{r2}} \leftarrow {\rm{r1}}}}}}. \\[-16pt] \end{split} $

      ${k_{i \leftarrow j}}~(i, j = {\rm{0}}, \;{\rm{l1}}, \;{\rm{l2}}, \;{\rm{r1}}, \;{\rm{r2}})$为从j态到i态的传输几率. 根据费米黄金定则, 传输几率分别为:

      $\begin{split} & {{k_{{\rm{l2}} \leftarrow {\rm{0}}}} = {\varGamma _{{\rm{l2}}}}f\left( {\dfrac{{{\varepsilon _{{\rm{l2}}}} - \mu }}{{{T_{\rm{L}}}}}} \right)},\\ & {{k_{0 \leftarrow {\rm{l}}2}} = } {\varGamma _{l2}}\left[ {1 - f\left( {\dfrac{{{\varepsilon _{{\rm{l2}}}} - \mu }}{{{T_{\rm{L}}}}}} \right)} \right], \end{split}$

      $\begin{split} & {{k_{{\rm{l1}} \leftarrow 0}} = {\varGamma _{{\rm{l1}}}}f\left( {\dfrac{{{\varepsilon _{{\rm{l1}}}} - \mu }}{{{T_{\rm{L}}}}}} \right)},\\ &{{k_{0 \leftarrow {\rm{l1}}}} = } {\varGamma _{{\rm{l1}}}}\left[ {1 - f\left( {\dfrac{{{\varepsilon _{{\rm{l}}1}} - \mu }}{{{T_{\rm{L}}}}}} \right)} \right], \end{split}$

      $\begin{split} & {{k_{{\rm{r}}2 \leftarrow 0}} = {\varGamma _{{\rm{r2}}}}f\left( {\dfrac{{{\varepsilon _{{\rm{r2}}}} - \mu }}{{{T_{\rm{R}}}}}} \right)},\\ &{{k_{0 \leftarrow {\rm{r2}}}} = }{\varGamma _{{\rm{r2}}}}\left[ {1 - f\left( {\dfrac{{{\varepsilon _{{\rm{r2}}}} - \mu }}{{{T_{\rm{R}}}}}} \right)} \right],\end{split}$

      $\begin{split} &{{k_{{\rm{r1}} \leftarrow 0}} = {\varGamma _{{\rm{r1}}}}f\left( {\dfrac{{{\varepsilon _{{\rm{r1}}}} - \mu }}{{{T_{\rm{R}}}}}} \right)},\\ &{{k_{0 \leftarrow {\rm{r}}1}} = } {\varGamma _{{\rm{r1}}}}\left[ {1 - f\left( {\dfrac{{{\varepsilon _{{\rm{r1}}}} - \mu }}{{{T_{\rm{R}}}}}} \right)} \right],\end{split}$

      其中, 比例常数${\varGamma _i}~(i = {\rm{l1}}, \;{\rm{l2}}, \;{\rm{r1}}, \;{\rm{r2}})$是导体库和量子点间的耦合系数, ${k_{i \leftarrow 0}}$表示电子由左(右)导体库传输到左(右)量子点的几率, 而${k_{0 \leftarrow j}}$表示电子由左(右)量子点传输到左(右)导体库的几率.

      $\begin{split} & {{k_{{\rm{r}}1 \leftarrow {\rm{l}}1}} = {\varGamma _{\rm{S}}}n\left( {\dfrac{{{\varepsilon _{\rm{g}}}}}{{{T_{\rm{S}}}}}} \right)},\\ &{{k_{{\rm{l}}1 \leftarrow {\rm{r}}1}} = }{\varGamma _{\rm{S}}}\left[ {1 + n\left( {\dfrac{{{\varepsilon _{\rm{g}}}}}{{{T_{\rm{S}}}}}} \right)} \right],\end{split}$

      $\begin{split} & {{k_{{\rm{l}}2 \leftarrow {\rm{r}}2}} = {\varGamma _{\rm{S}}}n\left( {\dfrac{{{\varepsilon _{\rm{g}}}}}{{{T_S}}}} \right)},\\ &{{k_{{\rm{r}}2 \leftarrow {\rm{l}}2}} = } {\varGamma _{\rm{S}}}\left[ {1 + n\left( {\dfrac{{{\varepsilon _{\rm{g}}}}}{{{T_{\rm{S}}}}}} \right)} \right],\end{split}$

      $\begin{split} &{{k_{{\rm{r}}2 \leftarrow {\rm{r}}1}} = {\varGamma _{n{\rm{r}}}}n\left( {\dfrac{{{\varepsilon _{{\rm{r}}2}} - {\varepsilon _{{\rm{r}}1}}}}{{{T_{\rm{R}}}}}} \right)},\\&{{k_{{\rm{r}}1 \leftarrow {\rm{r}}2}} = } {\varGamma _{n{\rm{r}}}}\left[ {1 + n\left( {\dfrac{{{\varepsilon _{{\rm{r}}2}} - {\varepsilon _{{\rm{r}}1}}}}{{{T_{\rm{R}}}}}} \right)} \right],\end{split}$

      $\begin{split} & {{k_{{\rm{l}}2 \leftarrow {\rm{l}}1}} = {\varGamma _{n{\rm{r}}}}n\left( {\dfrac{{{\varepsilon _{{\rm{l2}}}} - {\varepsilon _{{\rm{l}}1}}}}{{{T_{\rm{L}}}}}} \right)},\\ &{{k_{{\rm{l}}1 \leftarrow {\rm{l}}2}} = } {\varGamma _{n{\rm{r}}}}\left[ {1 + n\left( {\dfrac{{{\varepsilon _{{\rm{l}}2}} - {\varepsilon _{{\rm{l}}1}}}}{{{T_{\rm{L}}}}}} \right)} \right],\end{split}$

      其中, ${\varGamma _{\rm{S}}}$是两个量子点之间的耦合系数, 而${\varGamma _{n{\rm{r}}}}$是单个量子点内电子在两个能级间的跃迁系数; (6)和(7)式表示的是电子在两个量子点之间的传输几率, 而(8)和(9)式表示的是电子在单个量子点内两个能级之间的传输几率; ${\varepsilon _{\rm{g}}} = {\varepsilon _{{\rm{l}}2}} - {\varepsilon _{{\rm{r}}2}} = $${\varepsilon _{{\rm{r}}1}} - {\varepsilon _{{\rm{l}}1}} $为能级差.

      在系统达到稳态时, 存在${{\dot p}_0} = {{\dot p}_{{\rm{l}}1}} = {{\dot p}_{{\rm{l}}2}} = $${{\dot p}_{{\rm{r1}}}} = {{\dot p}_{{\rm{r}}2}} = 0 $, 且占有数几率满足${p_0} + {p_{{\rm{l1}}}} + {p_{{\rm{l}}2}} + $${p_{{\rm{r}}1}} + {p_{{\rm{r}}2}} = 1 $, 结合主方程可以解出各个占有数几率的表达式. 从左导体库通过能级${\varepsilon _{{\rm{l}}1}}$, ${\varepsilon _{{\rm{l}}2}}$流出的电子流${J_{{\rm{l}}1}}$, ${J_{{\rm{l}}2}}$为:

      ${J_{{\rm{l}}1}} = {k_{{\rm{l}}1 \leftarrow 0}}{p_0} - {k_{0 \leftarrow {\rm{l}}1}}{p_{{\rm{l1}}}},$

      ${J_{{\rm{l}}2}} = {k_{{\rm{l}}2 \leftarrow 0}}{p_0} - {k_{0 \leftarrow {\rm{l}}2}}{p_{{\rm{l}}2}}.$

      从能级${\varepsilon _{{\rm{r1}}}}$, ${\varepsilon _{{\rm{r2}}}}$流入右导体库的电子流${J_{{\rm{r1}}}}$, ${J_{{\rm{r}}2}}$为:

      $ {J_{{\rm{r}}1}} = {k_{0 \leftarrow \operatorname{r1}}}{p_{{\rm{r}}1}} - {k_{{\rm{r}}1 \leftarrow 0}}{p_0}, $

      ${J_{{\rm{r}}2}} = {k_{0 \leftarrow {\rm{r}}2}}{p_{{\rm{r}}2}} - {k_{{\rm{r}}2 \leftarrow 0}}{p_0}.$

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

      ${J_{{\rm{S}}1}} = {k_{{\rm{r}}1 \leftarrow {\rm{l}}1}}{p_{{\rm{l}}1}} - {k_{{\rm{l}}1 \leftarrow {\rm{r}}1}}{p_{{\rm{r}}1}},$

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

      ${J_{{\rm{S}}2}} = {k_{{\rm{l}}2 \leftarrow {\rm{r}}2}}{p_{{\rm{r}}2}} - {k_{{\rm{r}}2 \leftarrow {\rm{l}}2}}{p_{{\rm{l}}2}},$

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

      ${J_{n{\rm{r1}}}} = {k_{{\rm{l}}1 \leftarrow {\rm{l}}2}}{p_{{\rm{l}}2}} - {k_{{\rm{l}}2 \leftarrow {\rm{l}}1}}{p_{{\rm{l1}}}},$

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

      ${J_{n{\rm{r2}}}} = {k_{{\rm{r}}2 \leftarrow {\rm{r}}1}}{p_{{\rm{r}}1}} - {k_{{\rm{r}}1 \leftarrow {\rm{r}}2}}{p_{{\rm{r}}2}}.$

      ${J_1} = {J_{{\rm{l}}1}} + {J_{{\rm{l}}2}}$表示从左导体库流出的净电子流, ${J_2} = {J_{{\rm{r}}1}} + {J_{{\rm{r}}2}}$表示流入右导体库的净电子流, 系统遵守电子流守恒, 则${J_1} = {J_2}$. 电子流${J_{{\rm{S}}1}} = {J_{{\rm{l}}1}} + $${J_{n{\rm{r}}1}} $分为两部分, 一部分是来自左导体库的电子流, 另一部分是左量子点高能级向低能级跃迁的电子流. 同理电子流之间的关系为${J_{{\rm{S}}1}} = {J_{{\rm{l1}}}} + {J_{n{\rm{r}}1}} = $${J_{{\rm{r}}1}} + {J_{n{\rm{r2}}}} $, ${J_{{\rm{S}}2}} = - {J_{{\rm{l}}2}} + {J_{n{\rm{r}}1}} = - {J_{{\rm{r}}2}} + {J_{n{\rm{r}}2}}$.

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

      ${{\dot Q}_{\rm{L}}} = ({\varepsilon _{{\rm{l1}}}} - \mu ){J_{{\rm{l1}}}} + ({\varepsilon _{{\rm{l2}}}} - \mu ){J_{{\rm{l2}}}} - ({\varepsilon _{{\rm{l2}}}} - {\varepsilon _{{\rm{l}}1}}){J_{n{\rm{r}}1}},$

      ${{\dot Q}_{\rm{R}}} = - ({\varepsilon _{{\rm{r}}1}} - \mu ){J_{{\rm{r}}1}} - ({\varepsilon _{{\rm{r2}}}} - \mu ){J_{{\rm{r}}2}} + ({\varepsilon _{{\rm{r}}2}} - {\varepsilon _{{\rm{r1}}}}){J_{n{\rm{r}}2}},$

      ${{\dot Q}_{\rm{S}}} = {\varepsilon _{\rm{g}}}({J_{{\rm{S1}}}} + {J_{{\rm{S}}2}}).$

      该制冷机可以理解为: 从光子库吸收能量${{\dot Q}_{\rm{S}}}$, 促使右导体库放出热量${{\dot Q}_{\rm{R}}}$, 释放热量$ - {{\dot Q}_{\rm{L}}}$给左导体库. 其制冷系数定义为

      $\eta = \dfrac{{{{{\dot Q}}_{\rm{R}}}}}{{{{{\dot Q}}_{\rm{S}}}}} = \dfrac{{{\varepsilon _{{\rm{r}}2}}{J_{{\rm{S2}}}} - {\varepsilon _{{\rm{r}}1}}{J_{{\rm{S}}1}} + \mu \left( {{J_{{\rm{S}}1}} - {J_{{\rm{S2}}}}} \right)}}{{{\varepsilon _{\rm{g}}}\left( {{J_{{\rm{S1}}}} + {J_{{\rm{S}}2}}} \right)}},$

      ${{\dot Q}_{\rm{R}}}$为制冷率. 一般地, 制冷率和制冷系数是耦合系数(${\varGamma _{{\rm{l1}}}}$, ${\varGamma _{{\rm{l2}}}}$, ${\varGamma _{{\rm{r1}}}}$, ${\varGamma _{{\rm{r2}}}}$, $ {\varGamma _{\rm{S}}} $)、跃迁系数${\varGamma _{n{\rm{r}}}}$、能级(${\varepsilon _{\rm{g}}}$, ${\varepsilon _{{\rm{r}}1}}$, ${\varepsilon _{{\rm{r}}2}}$)化学势μ和温度(${T_{\rm{L}}}$, ${T_{\rm{R}}}$, ${T_{\rm{S}}}$)的函数. 我们首先使用一般的热力学理论来分析其性能. 熵产生率的表达式为

      ${\dot S} = - \dfrac{{{{{\dot Q}}_{\rm{L}}}}}{{{T_{\rm{L}}}}} - \dfrac{{{{{\dot Q}}_{\rm{R}}}}}{{{T_{\rm{R}}}}} - \dfrac{{{{{\dot Q}}_{\rm{S}}}}}{{{T_{\rm{S}}}}} \geqslant 0,$

      对于模型中任何参数的选择, 都要保证${\dot S}$不小于零. 结合能量守恒定律${{\dot Q}_L} + {{\dot Q}_R} + {{\dot Q}_S} = 0$, 熵产生率可以写成

      ${\dot S} = {{\dot Q}_{\rm{R}}}{F_{\rm{R}}} + {{\dot Q}_{\rm{S}}}{F_{\rm{S}}},$

      热力学力为${F_{\rm{R}}} = 1/{T_{\rm{L}}} - 1/{T_{\rm{R}}}$, ${F_{\rm{S}}} = 1/{T_{\rm{L}}} - 1/{T_{\rm{S}}}$, 且其范围满足${F_{\rm{R}}} \leqslant 0$, ${F_{\rm{S}}} > 0$. 利用(21)和(22)式可以得到制冷系数的表达式为

      $ \eta =\left( {1 - \dfrac{{{T_{\rm{L}}}}}{{{T_{\rm{S}}}}} - \dfrac{{{T_{\rm{L}}}\dot S}}{{{{\dot Q}_{\rm{S}}}}}} \right)\left( {\dfrac{{{T_{\rm{R}}}}}{{{T_{\rm{L}}} - {T_{\rm{R}}}}}} \right), $

      可以理解为: 式中第一部分是热机在温度为${T_{\rm{L}}}$${T_{\rm{S}}}$的热源之间工作时, 热机的卡诺效率被熵产生率相关因子降低, 这个因子是对所涉及过程的不可逆性的衡量; 第二部分是工作在温度为${T_{\rm{L}}}$${T_{\rm{R}}}$的热源之间的制冷机在可逆工作状态下的制冷系数. 当熵产生率为零${\dot S} = 0$, 制冷机达到可逆工作状态, 其制冷系数为可逆值${\eta _{\rm{r}}}=\left( {1 - \dfrac{{{T_{\rm{L}}}}}{{{T_{\rm{S}}}}}} \right)\left( {\dfrac{{{T_{\rm{R}}}}}{{{T_{\rm{L}}} - {T_{\rm{R}}}}}} \right)$.

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

      1)当化学势满足$\mu = ({\varepsilon _{{\rm{r1}}}} + {\varepsilon _{{\rm{r2}}}})/2$时, 根据各个电子流之间的关系, (18)—(20)式可以化简为

      ${{\dot Q}_{\rm{L}}} = - ({\varepsilon _{{\rm{r2}}}} + {\varepsilon _{\rm{g}}} - \mu )({J_{{\rm{S1}}}} + {J_{{\rm{S}}2}}),$

      ${{\dot Q}_{\rm{R}}} = ({\varepsilon _{{\rm{r}}2}} - \mu )({J_{{\rm{S1}}}} + {J_{{\rm{S2}}}}),$

      ${{\dot Q}_{\rm{S}}} = {\varepsilon _{\rm{g}}}({J_{{\rm{S1}}}} + {J_{{\rm{S2}}}}),$

      其中制冷率${{\dot Q}_{\rm{R}}}$与输入热流${{\dot Q}_{\rm{S}}}$成比例关系;

      2)调节参数使得${J_1} = {J_2} = 0$, 此时其他电子流之间还存在${J_{n{\rm{r}}1}} = {J_{n{\rm{r}}2}}$, ${J_{{\rm{S1}}}} = {J_{{\rm{S2}}}}$的关系, 因此${{\dot Q}_{\rm{R}}}$${{\dot Q}_{\rm{S}}}$同样成比例关系.

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

      $\eta = \dfrac{{{{{\dot Q}}_{\rm{R}}}}}{{{{{\dot Q}}_{\rm{S}}}}} = \dfrac{{{\varepsilon _{{\rm{r}}2}} - {\varepsilon _{{\rm{r1}}}}}}{{2{\varepsilon _{\rm{g}}}}},$

      只与能级(${\varepsilon _{\rm{g}}}, {\varepsilon _{{\rm{r1}}}}, {\varepsilon _{{\rm{r2}}}}$)有关. 熵产生率可以表达为${\dot S} = {{\dot Q}_{\rm{S}}}\left( {\eta \,{F_{\rm{R}}}{\rm{ + }}{F_{\rm{S}}}} \right)$, 当$\eta \,{F_{\rm{R}}}{\rm{ + }}{F_{\rm{S}}}\, =\;{\rm{0}}$时, 熵产生率达到${\dot S} = 0$. 因此, 只要模型参数满足关系: $\dfrac{{{\varepsilon _{{\rm{r2}}}} - {\varepsilon _{{\rm{r1}}}}}}{2}\left( {\dfrac{{{T_{\rm{R}}}}}{{{T_{\rm{L}}} - {T_{\rm{R}}}}}} \right)={\varepsilon _{\rm{g}}}\left( {1 - \dfrac{{{T_{\rm{L}}}}}{{{T_{\rm{S}}}}}} \right)$, 制冷机工作在可逆状态. 而紧耦合条件下模型参数的取值范围为

      $ 0 \leqslant \dfrac{{{\varepsilon _{{\rm{r2}}}} - {\varepsilon _{{\rm{r1}}}}}}{{2{\varepsilon _{\rm{g}}}}} \leqslant \left( {1 - \dfrac{{{T_{\rm{L}}}}}{{{T_{\rm{S}}}}}} \right)\left( {\dfrac{{{T_{\rm{R}}}}}{{{T_{\rm{L}}} - {T_{\rm{R}}}}}} \right). $

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

      图  2  紧耦合条件下制冷率与制冷系数在不同温度${T_{\rm{S}}}$下的关系

      Figure 2.  The relation curves of the cooling rate and the coefficient of performance at different temperature ${T_{\rm{S}}}$ under the condition of tight coupling.

      图  3  紧耦合条件下制冷率与制冷系数在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下的关系

      Figure 3.  The relation curves of the cooling rate and the coefficient of performance at different transition coefficient ${\varGamma _{n{\rm{r}}}}$ under the condition of tight coupling.

      图  4  一般情况下制冷率与制冷系数在不同温度${T_{\rm{S}}}$下的关系

      Figure 4.  The relation curves of the cooling rate and the coefficient of performance at different temperature ${T_{\rm{S}}}$ in the general case.

      图  5  一般情况下制冷率与制冷系数在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下的关系

      Figure 5.  The relation curves of the cooling rate and the coefficient of performanceat different transition coefficient ${\varGamma _{n{\rm{r}}}}$ in the general case.

      ${\dot Q}_{\rm{R}}^\eta \leqslant {{\dot Q}_{\rm{R}}} \leqslant {\dot Q}_{\rm{R}}^{\max },$

      ${\eta ^{{Q_{\rm{R}}}}} \leqslant \eta \leqslant {\eta ^{\max }},$

      这四个性能优化参数${\dot Q}_{\rm{R}}^{\max }$, ${\eta ^{{Q_{\rm{R}}}}}$, ${\eta ^{\max }}$, ${\dot Q}_{\rm{R}}^\eta $为制冷机的重要参数决定了制冷机性能的上下界限.

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

      $\dfrac{{\partial {{{\dot Q}}_{\rm{R}}}}}{{\partial {\varepsilon _{{\rm{r1}}}}}} = 0,\dfrac{{\partial {{{\dot Q}}_{\rm{R}}}}}{{\partial {\varepsilon _{{\rm{r2}}}}}} = 0,$

      $\dfrac{{\partial \eta }}{{\partial {\varepsilon _{{\rm{r1}}}}}} = 0,\dfrac{{\partial \eta }}{{\partial {\varepsilon _{{\rm{r2}}}}}} = 0,$

      可以得到这四个性能优化参数${\dot Q}_{\rm{R}}^{\max }$, ${\eta ^{{Q_{\rm{R}}}}}$, ${\eta ^{\max }}$, ${\dot Q}_{\rm{R}}^\eta $的数值解. 进而在紧耦合和一般情况下讨论了光子库温度、跃迁系数、温比$\tau = {T_{\rm{R}}}/{T_{\rm{L}}}$对四个优化性能参数的影响.

      图6图7为紧耦合条件下优化性能参数${\dot Q}_{\rm{R}}^{\max }$, ${\eta ^{{Q_{\rm{R}}}}}$随着温比变化的曲线图, 温度${T_{\rm{L}}}$的变化范围在${T_{\rm{R}}} < {T_{\rm{L}}} < {T_{\rm{S}}}$时, 制冷机才能有制冷作用, 因此图6中温度${T_{\rm{S}}}$不同时, 曲线起点不同. 在做图过程中令化学势$\mu /{k_{\rm{B}}} = 0\;{\rm{K}}$, ${\varepsilon _{{\rm{r1}}}} = - {\varepsilon _{{\rm{r2}}}}$, 其他参数的选取为${\varepsilon _{\rm{g}}}/{k_{\rm{B}}} = 2\;{\rm{K}}$, $ {\varGamma _{{\rm{l1}}}} = {\varGamma _{{\rm{l}}2}} = {\varGamma _{{\rm{r}}1}} = {\varGamma _{{\rm{r2}}}} = {\varGamma _{\rm{S}}} = \varGamma $, ${T_{\rm{R}}} = 4\;{\rm{K}}$. 图6中最大制冷率${\dot Q}_{\rm{R}}^{\max }$和对应的制冷系数${\eta ^{{Q_{\rm{R}}}}}$都随着温比的增大而增大. 图7中, 随跃迁系数${\varGamma _{n{\rm{r}}}}$的增加, 最大制冷率${\dot Q}_{\rm{R}}^{\max }$增加, 但对应的制冷系数${\eta ^{{Q_{\rm{R}}}}}$减少. 因此为了达到提高制冷率的目的, 应尽可能提高跃迁系数${\varGamma _{n{\rm{r}}}}$和温度${T_{\rm{S}}}$.

      图  6  在不同温度${T_{\rm{S}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随温比的变化

      Figure 6.  The curves of two optimal performance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different temperature ${T_{\rm{S}}}$

      图  7  在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随着温比的变化

      Figure 7.  The curves of two optimal performance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different transition coefficient ${\varGamma _{n{\rm{r}}}}$

      图8图10为一般情况下四个优化性能参数${\dot Q}_{\rm{R}}^{\max }$, ${\eta ^{{Q_{\rm{R}}}}}$, ${\eta ^{\max }}$, ${\dot Q}_{\rm{R}}^\eta $随着温比变化的曲线图, 在做图过程中令化学势$\mu /{k_{\rm{B}}} = 0\;{\rm{K}}$, 其他参数的选取为${\varepsilon _{\rm{g}}}/{k_{\rm{B}}} = 2\;{\rm{K}}$, $ {\varGamma _{{\rm{l}}1}} = {\varGamma _{{\rm{l}}2}} = {\varGamma _{{\rm{r}}1}} = {\varGamma _{{\rm{r}}2}} = {\varGamma _{\rm{S}}} = \varGamma $, ${T_{\rm{R}}} = $4 K. 图8图9中最大制冷率${\dot Q}_{\rm{R}}^{\max }$和对应的制冷系数${\eta ^{{Q_{\rm{R}}}}}$都随着温比的增加而增加, 图10中最大制冷系数随着温比的增大而增大, 而对应的制冷率趋于0. 要使得制冷机获得最大制冷率或者最大制冷系数应使左右两导体库的温度相近, 并提高光子库的温度和跃迁系数${\varGamma _{n{\rm{r}}}}$.

      图  8  在不同温度${T_{\rm{S}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随温比的变化

      Figure 8.  The curves of two optimalperformance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different temperature ${T_{\rm{S}}}$.

      图  10  在不同温度${T_{\rm{S}}}$下, 两个优化性能参数${\eta ^{\max }}$${\dot Q}_{\rm{R}}^\eta $随温比的变化

      Figure 10.  The curves of two optimal performance parameters ${\eta ^{\max }}$ and ${\dot Q}_{\rm{R}}^\eta $ changing with the temperature ratio at different temperature ${T_{\rm{S}}}$.

      图  9  在不同跃迁系数${\varGamma _{n{\rm{r}}}}$下, 两个优化性能参数${\dot Q}_{\rm{R}}^{\max }$${\eta ^{{Q_{\rm{R}}}}}$随着温比的变化

      Figure 9.  The curves of two optimal performance parameters ${\dot Q}_{\rm{R}}^{\max }$ and ${\eta ^{{Q_{\rm{R}}}}}$ changing with the temperature ratio at different transition coefficient ${\varGamma _{n{\rm{r}}}}$.

    • 本文主要研究结果如下:

      1)通过分析, 达到紧耦合的情况有两种, 一是使化学势μ与能级${\varepsilon _{{\rm{r}}2}}$, ${\varepsilon _{{\rm{r1}}}}$满足$\mu = ({\varepsilon _{{\rm{r1}}}} + {\varepsilon _{{\rm{r2}}}})/2$, 但其净电子流不为零; 二是使净电子流${J_{{\rm{l}}1}} + {J_{{\rm{l}}2}}$为零;

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

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

参考文献 (34)

目录

    /

    返回文章
    返回