-
热电装置是运用Seebeck效应或者Peltier效应将环境中的热量收集起来转化为电能的发电机或者在外加偏压的条件下对某一局部区域进行制冷的冷机. 然而, 目前热电装置的热电转换效率仍然很低. 近来, 人们对开发多端热电器件越来越感兴趣, 因为与传统的两端热电器件相比, 多端热电器件可以通过分开电流和热流来显著提高热力学性能[1].
2011年, Sánchez等[2]首次提出了一种基于库仑阻塞区中两个电容耦合量子点的三端能量转换器. 它可以通过从热库中吸收热能来发电, 这种转换器可以在最佳配置下实现卡诺效率. 后来, Thierschman等[3]通过实验证明了这种能量转换器的工作原理. 作为热机的反向操作, Zhang等[4-6]提出了一种基于两个耦合量子点的三端量子点制冷机, 详细分析了在考虑一些不可逆因素后这种制冷机的热力学性能. 对于基于更多耦合量子点的三端器件, Aniket等[7,8]提出了基于三量子点系统的实际非局部热机和制冷机, 基于主方程, 分析了它们的热力学性能. 一些学者理论上[9-14]和实验上[15-18]研究了基于共振隧穿量子点和电容耦合量子点的许多其他三端热电装置的性能, 取得长足的进展. 除了共振隧穿和耦合量子点之外, 其他纳米结构, 包括量子阱(或半导体超晶格[19-22]、量子霍尔棒[23,24]、纳米线[25,26]、热离子[27]和能量选择性隧道结[28]也被提出用于多端热电器件. 由磁振子、声子和光子等驱动源驱动的三端热电装置也得到了广泛的研究[29-37].
在三端结构热电装置的基础上, Whitney等[38]提出了一种基于三个电容耦合量子点的四端结构热电热机, 该热机即使在两个热库输入的总热功率为零的情况下, 通过热曳拉效应能达到对外输出功率的效果. Su等[39]提出了一种基于麦克斯韦妖驱动的三量子点系统的热力学泵, 他们利用马尔可夫随机热力学分析了驱动泵内热力学过程. Jiang等[40]分析了一种四端介观热电系统中库仑拖曳驱动的热电效应装置. Chen等[41]总结了量子热力学循环, 电子能量选择量子器件, 量子点热管理器件等在理论和实验方面研究的代表性成果和研究进展. 在前人工作的基础上, 本文提出了一种具有三个电容耦合量子点的四端混合驱动制冷机模型. 该模型可以通过高温的热库注入的热能和偏置电压输入的功率共同驱动来实现对低温库的制冷. 相对于三端制冷机, 它可以利用环境中多余的废热来驱动制冷机工作, 提高能源利用. 本文的重点是分析制冷率与制冷系数之间的热力学性能特征, 并对制冷机在最大制冷率条件下的主要性能参数进行了优化, 我们还比较了该制冷机在强电容耦合和弱电容耦合情况下的性能.
-
本文提出的四端热电制冷机的模型如图1(a)所示, 该装置是由三个耦合量子点、两个热库和左右两个电子库组成. 量子点
${\text{Q}}{{\text{D}}_{\text{M}}}$ 通过两个隧道结嵌入到左右两个电子库之间, 且相互之间可以进行粒子和能量的交换, 其中左侧电子库L的温度为$ {T_{\text{L}}} $ 以及化学势为$ {\mu _{\text{L}}} $ , 右侧电子库R的温度为$ {T_{\text{R}}} $ 以及化学势为$ {\mu _{\text{R}}} $ . 量子点$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 与温度为$ {T_{\text{H}}} $ , 化学势为$ {\mu _{\text{H}}} $ 的热库H耦合, 而量子点$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 与温度为$ {T_{\text{C}}} $ , 化学势为$ {\mu _{\text{C}}} $ 的中间温度的热库耦合. 三个量子点${\text{Q}}{{\text{D}}_{\text{M}}}$ ,$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 和$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 相互耦合并且通过长程库仑作用相互影响. 采用$ {U_{{\text{MH}}}} $ ,$ {U_{{\text{MC}}}} $ 和$ {U_{{\text{HC}}}} $ 来表示充电能. 如果量子点彼此距离较远, 它们可以桥接, 以获得强/弱耦合, 同时又能确保电子库和热库之间形成良好的热隔离. 四端结构热电热机的等效电路图如图1(b)所示.图 1 (a) 基于三量子点耦合的四端制冷机的示意图; (b) (a)中装置的等效电路图
Figure 1. (a) Schematic diagram of a four terminal refrigerator based on three coupled quantum dots, and (b) is the equivalent circuit diagram of the device in (a).
三个耦合的量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ ,${\text{Q}}{{\text{D}}_{\text{H}}}$ 和$ {\text{Q}}{{\text{D}}_{\text{C}}} $ , 都有各自的单能级$ {\varepsilon _{\text{M}}} $ ,$ {\varepsilon _{\text{H}}} $ 和$ {\varepsilon _{\text{C}}} $ . 中心系统的量子态可以表示为($ m,h,c $ ), 其中m, h和c分别表示各个量子点能级上的电子占据数. 因为库仑相互作用阻止两个电子同时出现在同一个能级上, 所以三个量子点的单能级$ {\varepsilon _{\text{M}}} $ ,$ {\varepsilon _{\text{H}}} $ 和$ {\varepsilon _{\text{C}}} $ 上各自只能有一个电子或者没有电子. 该系统存在8种量子态, 它们分别为 ($ 0, 0, 0 $ ), ($ 1, 0, 0 $ ), ($ 0, {\text{1}}, 0 $ ), ($ 0, 0, {\text{1}} $ ), ($ {\text{1}}, {\text{1}}, 0 $ ), (${\text{1}}, 0, {\text{1}}$ ), ($ {\text{0}}, {\text{1}}, {\text{1}} $ ), ($ {\text{1}}, {\text{1}}, {\text{1}} $ ), 8种量子态相互转换的动力学过程如图2所示, 其中$\left[ {{\varGamma _\alpha }} \right]_{mhc}^{m'h'c'}$ 表示由于电子隧穿进入库$ \alpha $ 或者隧穿流出库$ \alpha $ ($ \alpha = {\text{S, H, C}} $ 和$\rm S = {\rm{L, R }}$ )而使得该三量子点中心系统从量子态($m{{', }} h{{', }}c{'}$ )跃迁到量子态($m,h,c$ )的速率.在序贯隧穿近似中(
$ \hbar \gamma \ll {k_{\text{B}}}T $ ), 量子点的能级宽度可以忽略, 用主方程来描述量子态占有概率随时间的演化[2]:$ {\dot P_{mhc}} = \displaystyle\sum\limits_{m'h'c{'}} {\left[ M \right]_{mhc}^{m'h'c{'}}{P_{m'h'c{'}}}} , $ 其中
$ \displaystyle\sum\nolimits_{mhc} {{P_{mhc}} = {\text{1}}} $ ,$ {P_{mhc}} $ 代表量子态($ m, h, c $ )出现的概率. 跃迁矩阵所有的非对角元素是通过$\left[ M \right]_{mhc}^{m'h'c'} = \displaystyle\sum\nolimits_{\alpha = {\text{S, H, C}}} {\left[ {{\varGamma _\alpha }} \right]_{mhc}^{m'h'c'}}$ 来确定, 即三量子点系统从量子态($ {\text{m', h', c'}} $ )到量子态($ {\text{m, h, c}} $ )跃迁速率的总和.$\left[ {{\varGamma _\alpha }} \right]_{mhc}^{m'h'c'}{\text{ = }}\left[ {{\gamma _\alpha }} \right]_{mhc}^{m'h'c'}{f_\alpha }\left( {\Delta _{mhc}^{m'h'c'}} \right)$ , 这里的$ \left[ {{\gamma _\alpha }} \right]_{mhc}^{m'h'c'} $ 是三量子点系统从量子态($ m{{', }}h{{', }}c{'} $ )到量子态($ m{\text{, }}h{\text{, }}c $ )时量子点和对应库$ \alpha $ 之间的裸穿速率,$ {f_\alpha }\left( x \right) = 1/\left[ {1 + {{\text{e}}^{(x - {\mu _\alpha })/{k_{\text{B}}}{T_\alpha }}}} \right] $ 是费米函数, 它描述了化学势为$ {\mu _\alpha } $ 、温度为$ {T_\alpha } $ 的库$ \alpha $ 中的电子分布.$ \Delta _{mhc}^{m'h'c'} = {E_{mhc}} - {E_{m'h'c'}} $ ,$ {E_{mhc}} $ 是量子态为($ m{\text{, }}h{\text{, }}c $ )时的能量,$ {k_{\text{B}}} $ 是玻尔兹曼常数. 利用跃迁矩阵的每一列之和为零, 可以直接求出(1)式中的对角元素,$\left[ M \right]_{mhc}^{mhc} = - \displaystyle\sum\nolimits_{mhc \ne m'h'c{'}} {\left[ M \right]_{mhc}^{m'h'c'}}$ .根据主方程(1), 求出在稳态情况下
${\dot P_{mhc}} = 0$ , 每一个量子态($ m{\text{, }}h{\text{, }}c $ )的占有概率${P_{{\text{mhc}}}}$ 的表达式. 从左侧电子库$L$ 进入到量子点$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 中的电子流为$ {I_{\text{L}}} = \left[ {{\varGamma _{\text{L}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{L}}}} \right]_{000}^{100}{P_{100}} + \left[ {{\varGamma _{\text{L}}}} \right]_{110}^{010}{P_{010}} - \left[ {{\varGamma _{\text{L}}}} \right]_{010}^{110}{P_{110}} + \left[ {{\varGamma _{\text{L}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\text{L}}}} \right]_{001}^{101}{P_{101}} + \left[ {{\varGamma _{\text{L}}}} \right]_{111}^{011}{P_{011}} - \left[ {{\varGamma _{\text{L}}}} \right]_{011}^{111}{P_{111}}, $ 从右侧电子库R进入到量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 中的电子流为$ {I_R} = \left[ {{\varGamma _{\text{R}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{R}}}} \right]_{000}^{100}{P_{100}} + \left[ {{\varGamma _{\text{R}}}} \right]_{110}^{010}{P_{010}} - \left[ {{\varGamma _{\text{R}}}} \right]_{010}^{110}{P_{110}} + \left[ {{\varGamma _{\text{R}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\text{R}}}} \right]_{001}^{101}{P_{101}} + \left[ {{\varGamma _{\text{R}}}} \right]_{111}^{011}{P_{011}} - \left[ {{\varGamma _{\text{R}}}} \right]_{011}^{111}{P_{111}}, $ 稳态情况下, 从两侧电子库进入到量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的电荷流满足大小相等, 方向相反, 即$ I = {I_{\text{L}}} = - {I_{\text{R}}} $ , 则制冷机的输入功率为$ P = - I \times e\Delta V , $ 其中
$ \Delta V = {\mu _{\text{R}}} - {\mu _{\text{L}}} $ 为施加的偏置电压.从左侧电子库
$L$ 到量子点$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的热流为$ \begin{split} \;& {J_L} =\left( {{\varepsilon _M} - {\mu _L}} \right)\left( {\left[ {{\varGamma _{\text{L}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{L}}}} \right]_{000}^{100}{P_{100}}} \right) + \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MH}}}} - {\mu _{\text{L}}}} \right)\left( {\left[ {{\varGamma _{\text{L}}}} \right]_{110}^{010}{P_{010}} - \left[ {{\varGamma _{\text{L}}}} \right]_{010}^{110}{P_{110}}} \right) \\&~ + \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MC}}}} - {\mu _{\text{L}}}} \right)\left( {\left[ {{\varGamma _{\text{L}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\text{L}}}} \right]_{001}^{101}{P_{101}}} \right) + \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MH}}}} + {U_{{\text{MC}}}} - {\mu _{\text{L}}}} \right)\left( {\left[ {{\varGamma _{\text{L}}}} \right]_{111}^{011}{P_{011}} - \left[ {{\varGamma _{\text{L}}}} \right]_{011}^{111}{P_{111}}} \right), \end{split} $ 从右侧电子库R到量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的热流为$ \begin{split} \;&{J_{\text{R}}} = \left( {{\varepsilon _{\text{M}}} - {\mu _{\text{R}}}} \right)\left( {\left[ {{\varGamma _{\text{R}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{R}}}} \right]_{000}^{100}{P_{100}}} \right) + \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MH}}}} - {\mu _{\text{R}}}} \right)\left( {\left[ {{\varGamma _R}} \right]_{110}^{010}{P_{010}} - \left[ {{\varGamma _{\text{R}}}} \right]_{010}^{110}{P_{110}}} \right) \\& ~+ \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MC}}}} - {\mu _{\text{R}}}} \right)\left( {\left[ {{\varGamma _{\text{R}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\text{R}}}} \right]_{001}^{101}{P_{101}}} \right) + \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MH}}}} + {U_{{\text{MC}}}} - {\mu _{\text{R}}}} \right)\left( {\left[ {{\varGamma _{\text{R}}}} \right]_{111}^{011}{P_{011}} - \left[ {{\varGamma _{\text{R}}}} \right]_{011}^{111}{P_{111}}} \right), \end{split} $ 则从左右电子库流入量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的总热流为$ {J_{{\text{LR}}}} = {J_{\text{L}}} + {J_{\text{R}}}, $ 从高温热库
$H$ 到量子点$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 的热流为$ \begin{split} \;&{J_{\text{H}}} =\left( {{\varepsilon _{\text{H}}} - {\mu _{\text{H}}}} \right)\left( {\left[ {{\varGamma _{\text{H}}}} \right]_{010}^{000}{P_{000}} - \left[ {{\varGamma _{\text{H}}}} \right]_{000}^{010}{P_{010}}} \right) + \left( {{\varepsilon _{\text{H}}} + {U_{{\text{MH}}}} - {\mu _{\text{H}}}} \right)\left( {\left[ {{\varGamma _{\text{H}}}} \right]_{110}^{100}{P_{100}} - \left[ {{\varGamma _{\text{H}}}} \right]_{100}^{110}{P_{110}}} \right) \\&~ + \left( {{\varepsilon _{\text{H}}} + {U_{{\text{HC}}}} - {\mu _{\text{H}}}} \right)\left( {\left[ {{\varGamma _{\text{H}}}} \right]_{011}^{001}{P_{001}} - \left[ {{\varGamma _{\text{H}}}} \right]_{001}^{011}{P_{011}}} \right) + \left( {{\varepsilon _{\text{H}}} + {U_{{\text{MH}}}} + {U_{{\text{HC}}}} - {\mu _{\text{H}}}} \right)\left( {\left[ {{\varGamma _{\text{H}}}} \right]_{111}^{101}{P_{101}} - \left[ {{\varGamma _{\text{H}}}} \right]_{101}^{111}{P_{111}}} \right), \end{split} $ 从
${\text{C}}$ 库到量子点$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 的热流为$ \begin{split} \;&{J_{\text{C}}} = \left( {{\varepsilon _{\text{C}}} - {\mu _{\text{C}}}} \right)\left( {\left[ {{\varGamma _{\text{C}}}} \right]_{001}^{000}{P_{000}} - \left[ {{\varGamma _{\text{C}}}} \right]_{000}^{001}{P_{001}}} \right) + \left( {{\varepsilon _{\text{C}}} + {U_{{\text{MC}}}} - {\mu _{\text{C}}}} \right)\left( {\left[ {{\varGamma _{\text{C}}}} \right]_{101}^{100}{P_{100}} - \left[ {{\varGamma _{\text{C}}}} \right]_{100}^{101}{P_{101}}} \right) \\& ~ + \left( {{\varepsilon _{\text{C}}} + {U_{{\text{HC}}}} - {\mu _{\text{C}}}} \right)\left( {\left[ {{\varGamma _{\text{C}}}} \right]_{011}^{010}{P_{010}} - \left[ {{\varGamma _{\text{C}}}} \right]_{010}^{011}{P_{011}}} \right) + \left( {{\varepsilon _{\text{C}}} + {U_{{\text{MC}}}} + {U_{{\text{HC}}}} - {\mu _{\text{C}}}} \right)\left( {\left[ {{\varGamma _{\text{C}}}} \right]_{111}^{110}{P_{110}} - \left[ {{\varGamma _{\text{C}}}} \right]_{110}^{111}{P_{111}}} \right), \end{split} $ 将左右两电子库看作一个子系统
$LR$ , 整个制冷机满足热力学第一定律:$ {J_{\text{H}}}{\text{ + }}{J_{\text{C}}}{\text{ + }}{J_{\text{L}}} + {J_{\text{R}}} + P = 0, $ 和热力学第二定律:
$ \dot S = - \frac{{{J_{\text{H}}}}}{{{T_{\text{H}}}}} - \frac{{{J_{\text{C}}}}}{{{T_{\text{C}}}}} - \frac{{{J_{\text{L}}}}}{{{T_{\text{L}}}}} - \frac{{{J_{\text{R}}}}}{{{T_{\text{R}}}}} \geqslant 0. $ 当
$ {T_{\text{H}}} \gt {T_{{\text{L(R)}}}} \gt {T_{\text{C}}} $ 时, 该装置可作为一种利用从高温热库H流出的热流$ {J_{\text{H}}} $ 和外加电压输入的功率共同驱动来实现对低温热库C制冷的混合驱动制冷机, 称之为装置A. 此时的JH, JC, JLR, P四个参数应该满足:$ {J_{\text{H}}} \gt 0 \text{, }~ {J_{\text{C}}} \gt 0 \text{, }~ {J_{{\text{LR}}}} \lt 0 \text{, }~ P \gt 0 , $ 此时, JC可被视为制冷机的制冷率, 此时的制冷系数定义为
$ {\eta _{{\text{COP}}}}{\text{ = }}\frac{{{J_{\text{C}}}}}{{{J_{\text{H}}} + P}}, $ 当输入功率
$ P = 0 $ 时, 该装置可作为一个纯热驱动制冷机, 制冷系数${\eta _{{\text{COP}}}}{\text{ = }}{{{J_{\text{C}}}}}/{{{J_{\text{H}}}}}$ . 当高温热库输入的热流$ {J_{\text{H}}} = 0 $ 时, 该装置可作为一个纯电驱动制冷机, 制冷系数${\eta _{{\text{COP}}}}{\text{ = }}{{{J_{\text{C}}}}}/{P}$ .当
$ {T_{\text{H}}} \gt {T_{\text{C}}} \gt {T_{{\text{L(R)}}}} $ 时, 该装置可作为一种利用高温热库H释放的热流$ {J_{\text{H}}} $ 和外加电压输入的功率共同驱动来实现对低温系统LR制冷的混合驱动制冷机, 称之为装置B. 此时的JH, JC, JLR, P四个参数应该满足:$ {J_{\text{H}}} \gt 0 \text{, }~ {J_{\text{C}}} \lt 0 \text{, }~ {J_{{\text{LR}}}} \gt 0 \text{, }~ P > 0 , $ 其中
$ {J_{{\text{LR}}}} $ 为制冷率, 制冷系数定义为$ {\eta _{{\text{COP}}}}{\text{ = }}\frac{{{J_{{\text{LR}}}}}}{{{J_{\text{H}}} + P}}. $ 采用了一套物理上合理的参数, 三个量子点的能级为
$ {\varepsilon _{\text{H}}} = {\varepsilon _{\text{C}}} = {\varepsilon _{\text{M}}} = 0 $ , 四个库的温度满足${k_{\text{B}}}{T_{\text{L}}}/\gamma = {k_{\text{B}}}{T_{\text{R}}}/\gamma = 12 + {k_{\text{B}}}\Delta T/\gamma$ ,$ {k_{\text{B}}}{T_{\text{C}}}/\gamma = 12 - {k_{\text{B}}}\Delta T/\gamma $ ,$ {k_{\text{B}}}{T_{\text{H}}}{\text{/}}\gamma = 18 $ , 四个库的化学势为$ {\mu _{\text{H}}} = {\mu _{\text{C}}} = {\mu _{\text{L}}} = 0 $ ,$ {\mu _{\text{R}}} = e\Delta V $ , 且充电能为$ {U_{{\text{MH}}}}/\gamma = 50 $ ,$ {U_{{\text{MC}}}}/\gamma = 20 $ ,$ {U_{{\text{HC}}}}/\gamma = 10 $ , 裸穿率为$ \left[ {{\gamma _{\text{H}}}} \right]_{mhc}^{m'h'c'} = \left[ {{\gamma _{\text{C}}}} \right]_{mhc}^{m'h'c'} = \gamma $ ,$ \left[ {{\gamma _{\text{L}}}} \right]_{mhc}^{m'h'c'} = \, \gamma $ 除了$ \left[ {{\gamma _{\text{L}}}} \right]_{100}^{000} = \left[ {{\gamma _{\text{L}}}} \right]_{000}^{100} = \lambda \gamma $ ,$\left[ {{\gamma _{\text{R}}}} \right]_{mhc}^{m'h'c'} = \lambda \gamma$ 除了$ \left[ {{\gamma _{\text{R}}}} \right]_{100}^{000} = \left[ {{\gamma _{\text{R}}}} \right]_{000}^{100} = \, \gamma \, $ .$ \lambda $ ($ {\text{0}} \leqslant \lambda \leqslant {\text{1}} $ )为耗散系数, 当$ \lambda = 0 $ 时, 在两量子点$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 和$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 都为空的情况下, 量子点$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 会与左侧电子库$L$ 解耦. 当$ \lambda \ne 0 $ 时, 在两量子点$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 和$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 均被占据的情况下, 量子点$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 与左右电子库都是耦合的.可以绘制出当耗散系数为
$ \lambda = 0 $ 时, 输入功率$ P $ 以及各量子点与对应库之间的热流JH, JC, JLR随温差$ \Delta T $ 和偏置电压$ e\Delta V $ 变化的三维投影图. 如图3所示, 图中的黑线对应$ P = 0 $ , 红线代表JH = 0, 蓝线表示$ {J_{\text{C}}} = 0 $ , 黄线表示$ {J_{{\text{LR}}}} = 0 $ .图 3 当耗散系数为
$ \lambda = 0 $ 时, (a) 输入功率$ P $ 以及(b)—(d)各量子点与对应库之间的热流JH, JC, JLR随温差$ \Delta T $ 和偏置电压$ e\Delta V $ 变化的三维投影图Figure 3. The three-dimensional projection graphs for (a) input power
$ P $ and (b)–(d) the heat flow JH, JC, JLR varying with the temperature difference$ \Delta T $ and the bias voltage$ e\Delta V $ under the dissipation factor$ \lambda = 0 $ .当耗散系数为
$ \lambda = 0 $ 和$ \lambda = 0.26 $ 时, 装置A和装置B各自对应的工作区域, 如图4所示. 通过数值计算, 发现装置A和装置B对应的工作区域都随着耗散系数的增大而减小. 且当$ \lambda \geqslant 0.015 $ 时, 装置B会消失. 当$ \lambda \leqslant 0.26 $ 时, 装置A只有一种工作模式, 即混合驱动制冷机, 当$ \lambda \geqslant 0.26 $ , 装置A有两种不同的工作模式, 即混合驱动制冷机和纯电驱动制冷机. 与装置A相比, 装置B更为复杂, 因为子系统LR既是电力的驱动端, 又是制冷的目标. 虽然装置B的制冷率明显小于装置A, 但可以在装置A不能工作的区域工作, 进而实现不同的制冷效果.图 4 当耗散系数为
$ \lambda = 0 $ 和$ \lambda = 0.26 $ 时, 装置A和装置B各自对应的工作区域Figure 4. The corresponding working areas of device A and device B when the dissipation factor is (a)
$ \lambda = 0 $ and (b)$ \lambda = 0.26 $ , respectively.本文仅分析装置A即混合驱动制冷机的性能特征. 可以画出在不同耗散系数下, 装置作为混合驱动制冷机时的总熵产率随温差
$ \Delta T $ 和偏置电压$ e\Delta V $ 变化的三维曲线, 如图5所示. 其中绿色面对应$ \lambda = 0 $ , 红色面对应$ \lambda = 0.1 $ , 蓝色面对应$ \lambda = 0.2 $ . 发现总熵产率满足的$ \dot S \gt 0 $ , 并且$ \dot S $ 随着耗散系数的增加而增加.图 5 在不同耗散系数下, 装置作为混合驱动制冷机时的总熵产率随温差
$ \Delta T $ 和偏置电压$ e\Delta V $ 变化的三维曲线Figure 5. The three-dimensional curves of total entropy production rate varying with temperature difference
$ \Delta T $ and bias voltage$ e\Delta V $ when the device is used as a hybrid-driven refrigerator under different dissipation factor.为了简化制冷率和制冷系数并明确各种参数之间的关系, 制冷率和制冷系数可以重新写成
$ {J_{\text{C}}} = {J_{\text{C}}}(eV,\lambda ,\Delta T,{U_{{\text{MH}}}},{U_{{\text{MC}}}}), $ $ \eta = \eta (eV,\lambda ,\Delta T,{U_{{\text{MH}}}},{U_{{\text{MC}}}}). $ 根据以上各公式, 对于给定的
$ \Delta T = 2\gamma /{k_{\text{B}}} $ 和$e\Delta V = 50$ , 绘制出在不同的耗散系数下制冷率和制冷系数随充电能$ {U_{{\text{MH}}}} $ 和$ {U_{{\text{MC}}}} $ 变化的三维图, 如图6所示. 可以看出对于耗散系数为$ \lambda = 0 $ , 当满足${U_{{\text{MH}}}}/\gamma = 50$ 和$ {U_{{\text{MC}}}}/\gamma \simeq 19.5 $ 时, 可以得到最大制冷率, 当满足$ {U_{{\text{MH}}}}/\gamma = 50 $ 和$ {U_{{\text{MC}}}}/\gamma \simeq 23 $ 时, 可以得到最大制冷系数. 且最大制冷率和最大制冷系数均随耗散系数的增大而减小.图 6 在不同耗散系数下, 制冷率和制冷系数随充电能
$ {U_{{\text{MH}}}} $ 和$ {U_{{\text{MC}}}} $ 变化的三维图Figure 6. The three-dimensional diagrams of the cooling rate and the COP varying with charging energy
$ {U_{{\text{MH}}}} $ and$ {U_{{\text{MC}}}} $ under different dissipation factor.在不同的耗散系数下制冷率和制冷系数随充电能
$ {U_{{\text{MH}}}} $ 和偏置电压$ eV $ 变化的三维图, 如图7所示. 可以看出, 制冷率随着偏置电压的增大而增大. -
根据 (1)—(5) 式, (10)和(11)式以及极值条件:
$ \frac{{\partial {J_{\text{C}}}}}{{\partial {U_{{\text{MH}}}}}} = 0 {,~~ }\frac{{\partial {J_{\text{C}}}}}{{\partial {U_{{\text{MC}}}}}} = 0 {,~~ } \frac{{\partial {J_{\text{C}}}}}{{\partial eV}} = 0 . $ 在给定条件
$ \Delta T = 2\gamma /{k_{\text{B}}} $ 下, 优化的制冷率$ {J_{{\text{Copt}}}} $ 和优化制冷率对应的制冷系数随耗散系数$ \lambda $ 的变化曲线, 如图8所示. 优化的制冷率和对应的制冷系数随耗散系数$ \lambda $ 的增大而单调递减, 只有当$ \lambda = 0 $ 时, 才能得到优化的制冷率最大值和对应的制冷系数的最大值. 为了获得最大的制冷率和相应的制冷系数, 应该选择尽可能小的耗散系数. 本文采用了$ {U_{{\text{MH}}}}/\gamma = 50 $ ,$ e\Delta V = 50 $ 进行数值计算.图 8 在给定条件
$ \Delta T = 2\gamma /{k_{\text{B}}} $ 下, (a) 优化的制冷率$ {J_{{\text{Copt}}}} $ 和(b)优化制冷率对应的制冷系数以及(c)对应的充电能$ {U_{{\text{MC}}}} $ 随耗散系数$ \lambda $ 的变化曲线Figure 8. The curves of (a) the optimized cooling rate
$ {J_{{\text{Copt}}}} $ and (b) the COP corresponding to optimized cooling rate and (c) the corresponding charging energy$ {U_{{\text{MC}}}} $ as a function of dissipation factor$ \lambda $ under the given condition$ \Delta T = 2\gamma /{k_{\text{B}}} $ .分析温差
$ \Delta T $ 对混合驱动制冷机最优性能的影响. 利用(16)和(17)式, 只对最大制冷率对应的充电能$ {U_{{\text{MC}}}} $ 进行优化, 而直接取充电能、优化偏压和优化耗散系数为$ {U_{{\text{MH}}}}/\gamma = 50 $ ,$ e\Delta V = 50 $ 和$ \lambda = 0 $ . 可以绘制出最大制冷率$ {J_{{\text{C}}\max }} $ 、最大制冷率对应的COP$ {\eta _{\text{J}}} $ 、最佳充电能量$ {U_{{\text{MC}}}} $ 随温差$ \Delta T $ 的变化曲线, 如图9所示. 从图9(a)可以看出, 最大制冷率$ {J_{{\text{C}}\max }} $ 、最大制冷率对应的COP$ {\eta _{\text{J}}} $ 会随着温差$ \Delta T $ 的增大而逐渐减小. 由图9(b)可知, 最佳充电能量$ {U_{{\text{MC}}}} $ 随温差$ \Delta T $ 的增大而单调减小.图 9 (a) 最大制冷率
$ {J_{{\text{C}}\max }} $ 和最大制冷率对应的COP$ {\eta _{\text{J}}} $ ; (b) 最佳充电能量$ {U_{{\text{MC}}}} $ 随温差$ \Delta T $ 的变化曲线Figure 9. The curves of (a) the maximum cooling rate
$ {J_{{\text{C}}\max }} $ and the COP corresponding to maximum cooling rate$ {\eta _{\text{J}}} $ and (b) the optimal charging energy$ {U_{{\text{MC}}}} $ as a function of temperature difference$ \Delta T $ . -
当量子点
$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 和$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 之间是强耦合$ {U_{{\text{HC}}}} $ 时, 则量子点$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 和$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 同时被占据的可能性就变得很小, 所以就不存在量子态($ {\text{0}}, {\text{1}}, {\text{1}} $ )和($ {\text{1}}, {\text{1}}, {\text{1}} $ ), 此时中心系统仅存在6个量子态. 电荷涨落可分为两个循环, 第一个循环是:$ \begin{split} (0, 0, 0)\to ( 1, 0, 0) \to ( {\text{1}}, {\text{1}}, 0 )\to ( 0, 1, 0 )\to ( 0, 0, 0) \end{split}$ 涉及到系统和高温热库之间的交换能量
$ {U_{MH}} $ . 第二个循环是:$ \begin{split} (0, 0, 0)\to ( 1, 0, 0) \to ( {\text{1}}, 0, 1 )\to ( 0, 0, 1 )\to ( 0, 0, 0) \end{split} $ 涉及到系统和次高温热库之间的交换能量
$ {U_{{\text{MC}}}} $ . 从左侧电子库L到量子点$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的电子流表达式为$ \begin{split} \;&{I_{\text{L}}} = \left[ {{\varGamma _{\text{L}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{L}}}} \right]_{000}^{100}{P_{100}} + \left[ {{\varGamma _{\text{L}}}} \right]_{110}^{010}{P_{010}} \\& ~~ - \left[ {{\varGamma _{\text{L}}}} \right]_{010}^{110}{P_{110}} + \left[ {{\varGamma _{\text{L}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\text{L}}}} \right]_{001}^{101}{P_{101}}, \end{split} $ 从右侧电子库R到量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的电子流表达式为$ \begin{split} \;& {I_{\text{R}}} = \left[ {{\varGamma _{\text{R}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{R}}}} \right]_{000}^{100}{P_{100}} + \left[ {{\varGamma _{\text{R}}}} \right]_{110}^{010}{P_{010}} \\& ~~ - \left[ {{\varGamma _{\text{R}}}} \right]_{010}^{110}{P_{110}} + \left[ {{\varGamma _{\text{R}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\text{R}}}} \right]_{001}^{101}{P_{101}}, \end{split}$ 从左侧电子库L到量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的热流表达式为$ \begin{split} \;& {J_{\text{L}}} = \left( {{\varepsilon _{\text{M}}} - {\mu _{\text{L}}}} \right)\left( {\left[ {{\varGamma _{\text{L}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\text{L}}}} \right]_{000}^{100}{P_{100}}} \right) \\& ~~ + \left( {{\varepsilon _{\text{M}}} + {U_{{\text{MH}}}} - {\mu _{\text{L}}}} \right) \left( {\left[ {{\varGamma _{\text{L}}}} \right]_{110}^{010}{P_{010}} - \left[ {{\varGamma _{\text{L}}}} \right]_{010}^{110}{P_{110}}} \right)\\& ~~ + ( {{\varepsilon _{\text{M}}} + {U_{{\text{MC}}}} - {\mu _{\text{L}}}} )\left( [ {{\varGamma _{\text{L}}}} ]_{101}^{001}{P_{001}} - [ {{\varGamma _{\text{L}}}} ]_{001}^{101}{P_{101}} \right),\\ \end{split} $ 从右侧电子库R到量子点
$ {\text{Q}}{{\text{D}}_{\text{M}}} $ 的热流表达式为$\begin{split} &\\ {J_{\rm{R}}} =\;& \left( {{\varepsilon _{\rm{M}}} - {\mu _{\rm{R}}}} \right)\left( {\left[ {{\varGamma _{\rm{R}}}} \right]_{100}^{000}{P_{000}} - \left[ {{\varGamma _{\rm{R}}}} \right]_{000}^{100}{P_{100}}} \right) \\& + \left( {{\varepsilon _{\rm{M}}} + {U_{\rm MH}} - {\mu _{\rm{R}}}} \right)\\&\times\left( {\left[ {{\varGamma _{\rm{R}}}} \right]_{110}^{010}{P_{010}} - \left[ {{\varGamma _{\rm{R}}}} \right]_{010}^{110}{P_{110}}} \right) \\& + \left( {{\varepsilon _{\rm{M}}} + {U_{\rm MC}} - {\mu _{\rm{R}}}} \right)\\&\times\left( {\left[ {{\varGamma _{\rm{R}}}} \right]_{101}^{001}{P_{001}} - \left[ {{\varGamma _{\rm{R}}}} \right]_{001}^{101}{P_{101}}} \right) , \end{split} $ 从高温热库H到量子点
$ {\text{Q}}{{\text{D}}_{\text{H}}} $ 的热流表达式为$\begin{split} &{J_{\text{H}}} = \left( {{\varepsilon _{\text{H}}} - {\mu _{\text{H}}}} \right)\left( {\left[ {{\varGamma _{\text{H}}}} \right]_{010}^{000}{P_{000}} - \left[ {{\varGamma _{\text{H}}}} \right]_{000}^{010}{P_{010}}} \right) \\& + ( {{\varepsilon _{\text{H}}} + {U_{{\text{MH}}}} - {\mu _{\text{H}}}} ) \Big( [ {\varGamma _{\text{H}}} ]_{110}^{100}{P_{100}} - [ \varGamma _{\text{H}} ]_{100}^{110}{P_{110}} \Big), \end{split} $ 从次高温热库C到量子点
$ {\text{Q}}{{\text{D}}_{\text{C}}} $ 的热流表达式为$ \begin{split} & {J_{\text{C}}} = \left( {{\varepsilon _{\text{C}}} - {\mu _{\text{C}}}} \right)\left( {\left[ {{\varGamma _{\text{C}}}} \right]_{001}^{000}{P_{000}} - \left[ {{\varGamma _{\text{C}}}} \right]_{000}^{001}{P_{001}}} \right) \\& + \left( {{\varepsilon _{\text{C}}} + {U_{{\text{MC}}}} - {\mu _{\text{C}}}} \right) \left( {[ {{\varGamma _{\text{C}}}} ]_{101}^{100}{P_{100}} - [ {{\varGamma _{\text{C}}}} ]_{100}^{101}{P_{101}}} \right). \end{split} $ 通过数值计算, 绘制出该混合驱动制冷机在强耦合
$ {U_{{\text{HC}}}} $ 的情况下, 当$ \lambda = 0 $ 时制冷率$ {J_{\text{C}}} $ 和制冷系数$ {\eta _{{\text{COP}}}} $ 随$ {U_{{\text{MH}}}} $ 和$ {U_{{\text{MC}}}} $ 变化的三维图, 如图10所示. 与弱耦合的情况相比, 强耦合作用下最大制冷率和制冷系数要大得多.图 10 在强耦合
$ {U_{{\text{HC}}}} $ 的情况下, 当$ \lambda = 0 $ 时, (a) 制冷率$ {J_{\text{C}}} $ 和(b)制冷系数$ {\eta _{{\text{COP}}}} $ 随着$ {U_{{\text{MH}}}} $ 和$ {U_{{\text{MC}}}} $ 变化的三维图Figure 10. In the case of strong coupling
$ {U_{{\text{HC}}}} $ , when$ \lambda = 0 $ , the three-dimensional diagrams for (a) cooling rate$ {J_{\text{C}}} $ and (b) the COP varying with$ {U_{{\text{MH}}}} $ and$ {U_{{\text{MC}}}} $ . -
本文研究了利用高温热源释放的热能和外加电压输入电能共同驱动来实现对低温库制冷的四端三量子点制冷机的性能. 主要研究结论: 1)混合驱动制冷机的工作区域、制冷率和制冷系数均随耗散系数的增大而减小; 2)最大制冷率和对应的制冷系数COP随着温差
$ \Delta T $ 的增大而减小; 3)通过合理选择充电能$ {U_{{\text{MC}}}} $ , 使$ \lambda = 0 $ , 并使充电能$ {U_{{\text{MH}}}} $ 和$ e\Delta V $ 尽可能大, 可以使混合驱动制冷机在最大制冷率的最佳状态下运行; 4)与弱耦合的情况相比, 强耦合作用下最大制冷率和制冷系数要大得多.
-
In this paper, a four-terminal hybrid driven refrigerator model with three capacitively coupled quantum dots is proposed, which can be driven by the energy current injected from the highest temperature thermal reservoir and the power input to achieve the refrigeration of the low temperature reservoir. Based on the master equation we derive the expressions for charge current and heat current between three quantum dots and thermal reservoirs in the weak/strong capacitive coupling case, respectively. We numerically analyze the thermodynamic performance characteristics of the refrigerator between the cooling rate and the coefficient of performance, and the main performance parameters of the refrigerator are optimized under the condition of the maximum cooling rate. Finally, we compare the performance of this refrigerator in the strong capacitive coupling case with that in the weak capacitive coupling case.
-
Keywords:
- couple quantum dot /
- four-terminal hybrid refrigerator /
- cooling rate /
- coefficient of performance
[1] Chen L, Ding Z, Sun F 2011 Energy 36 4011
Google Scholar
[2] Sánchez R, Büttiker M 2011 Phys. Rev. B 83 085428
Google Scholar
[3] Thierschmann H, Sánchez R, Sothmann B, Arnold F, Heyn C, Hansen W, Buhmann H, Molenkamp L W 2015 Nat. Nanotechnol 10 854
Google Scholar
[4] Zhang Y C, Lin G X, Chen J C 2015 Phys. Rev. E 91 052118
Google Scholar
[5] Zhang Y C, Huang C K, Lin G X, Chen J C 2015 Energy 85 200
Google Scholar
[6] Zhang Y C, Wang Y, Huang C K, Lin G X, Chen J C 2016 Energy 95 593
Google Scholar
[7] Aniket S. 2020 J. Appl. Phys 127 234903
Google Scholar
[8] Anamika B, Surojit H, Shailendra K, Varshney, Gourab D, Aniket S 2021 Phys. Rev. E 103 012131
Google Scholar
[9] Kano S, Fujii M 2017 Nanotechnology 28 095403
Google Scholar
[10] Lim J S, Sánchez D, López R 2018 N. J. Phys. 20 023038
Google Scholar
[11] Daré A M, Lombardo P 2017 Phys. Rev. B 96 115414
Google Scholar
[12] Su H, Shi Z C, He J Z 2015 Chin. Phys. Lett. 32 100501
Google Scholar
[13] 苏豪, 王家伟, 赵沁园, 何济洲 2016 中国科学: 技术科学 46 1296
Google Scholar
Su H, Wang J W, Zhao Q Y, He J Z 2016 Sci. Sin-Tech. 46 1296
Google Scholar
[14] Shi Z C, Qin W F, He J Z 2016 Mod. Phys. Lett. B 30 1650397
Google Scholar
[15] Roche B, Roulleau P, Jullien T, Jompol, Y, Farrer I, Ritchie D A, Glattli D C. 2015 Nat. Commun 6 6738
Google Scholar
[16] Hartmann F, Pfeffer P, Höfling S, Kamp M, Worschech L. 2015 Phys. Rev. Lett. 114 146805
Google Scholar
[17] Josefsson M, Svilans A, Burke A, Hoffmann E, Fahlvik S, Thelander C, Leijnse M, Linke H 2018 Nat. Nanotechnol 13 920
Google Scholar
[18] Keller A J, Lim J S, Sánchez D, López R, Amasha S, Katine J A 2016 Phys. Rev. Lett. 117 066602
Google Scholar
[19] Sothmann B, Sánchez R, Jordan A N, Büttiker M 2013 N. J. Phys. 15 095021
Google Scholar
[20] Choi Y, Jordan A N 2015 Physica E:Low-dimensional Systems Nanostruct 74 465
Google Scholar
[21] Lin Z B, Yang Y Y, Fu J, Li W, He J Z 2019 Chin. Phys. Lett. 36 060501
Google Scholar
[22] Lin Z B, Li W, Yang Y Y, He J Z 2020 Phys. Rev. E 101 022117
Google Scholar
[23] Sothmann B, Sánchez R, Jordan A N 2014 Europhys. Lett. 107 47003
Google Scholar
[24] Sánchez R, Sothmann B, Jordan A N 2015 Phys. Rev. Lett. 114 146801
Google Scholar
[25] Boukai A I, Bunimovich Y, Tahir-Kheli J, Yu J K, Goddard I W A, Heath J R 2008 Nature 451 168
Google Scholar
[26] Yang Y Y, Xu S, Li W, He J Z 2020 Phys. Scr. 95 095001
Google Scholar
[27] Yang Y Y, Xu S, He J Z 2020 Chin. Phys. Lett. 37 120502
Google Scholar
[28] Su S, Zhang Y, Chen J, Shih T M 2016 Sci. Reports 6 21425
Google Scholar
[29] Shi Z C, Fu J, Qin W F, He J Z 2017 Chin. Phys. Lett. 34 110501
Google Scholar
[30] Jiang J H, Entin-Wohlman O, Imry Y 2013 New Journal of Physics 15 075021
Google Scholar
[31] Li C, Zhang Y, He J 2013 Chin. Phys. Lett. 30 100501
Google Scholar
[32] Rutten B, Esposito M, Cleuren B 2009 Phys. Rev. B 80 235122
Google Scholar
[33] Cleuren B, Rutten B, Van den Broeck C 2012 Phys. Rev. Lett. 108 120603
Google Scholar
[34] 施志诚, 何济洲, 肖宇玲 2015 中国科学: 物理学 力学 天文学 45 050502
Google Scholar
Shi Z C, He J Z, Xiao Y L 2015 Sci. Sin-Phys. Mech. Astron. 45 050502
Google Scholar
[35] Li C, Zhang Y, Wang J, He J 2013 Phys. Rev. E 88 062120
Google Scholar
[36] Wang J H, Lai Y M, Ye Z L, He J Z, Ma Y L, Liang Q H 2015 Phys. Rev. E 91 050102
Google Scholar
[37] 李唯, 符婧, 杨贇贇, 何济洲 2019 物理学报 68 220501
Google Scholar
Li W, Fu J, Yang Y Y, He J Z 2019 Acta Phys. Sin. 68 220501
Google Scholar
[38] Whitney R S, Sánchez R, Haupt F, Splettstoesser J 2016 Phys. E 82 176
Google Scholar
[39] Fu T, Du J Y, Su S H, Su G Z, Chen J C 2021 Eur. Phys. J. Plus 136 1059
Google Scholar
[40] Xi M M, Wang R Q, Lu J C, Chen T Y, Jiang J H 2021 Chin. Phys. Lett 38 088801
Google Scholar
[41] 苏山河, 张艳超, 彭万里, 苏国珍, 陈金灿 2021 中国科学: 物理学 力学 天文学 51 112
Google Scholar
Su S H, Zhang Y C, Peng W L, Su G Z, Chen J C 2021 Sci. Sin-Phys. Mech. Astron. 51 112
Google Scholar
-
图 3 当耗散系数为
$ \lambda = 0 $ 时, (a) 输入功率$ P $ 以及(b)—(d)各量子点与对应库之间的热流JH, JC, JLR随温差$ \Delta T $ 和偏置电压$ e\Delta V $ 变化的三维投影图Fig. 3. The three-dimensional projection graphs for (a) input power
$ P $ and (b)–(d) the heat flow JH, JC, JLR varying with the temperature difference$ \Delta T $ and the bias voltage$ e\Delta V $ under the dissipation factor$ \lambda = 0 $ .图 5 在不同耗散系数下, 装置作为混合驱动制冷机时的总熵产率随温差
$ \Delta T $ 和偏置电压$ e\Delta V $ 变化的三维曲线Fig. 5. The three-dimensional curves of total entropy production rate varying with temperature difference
$ \Delta T $ and bias voltage$ e\Delta V $ when the device is used as a hybrid-driven refrigerator under different dissipation factor.图 8 在给定条件
$ \Delta T = 2\gamma /{k_{\text{B}}} $ 下, (a) 优化的制冷率$ {J_{{\text{Copt}}}} $ 和(b)优化制冷率对应的制冷系数以及(c)对应的充电能$ {U_{{\text{MC}}}} $ 随耗散系数$ \lambda $ 的变化曲线Fig. 8. The curves of (a) the optimized cooling rate
$ {J_{{\text{Copt}}}} $ and (b) the COP corresponding to optimized cooling rate and (c) the corresponding charging energy$ {U_{{\text{MC}}}} $ as a function of dissipation factor$ \lambda $ under the given condition$ \Delta T = 2\gamma /{k_{\text{B}}} $ .图 9 (a) 最大制冷率
$ {J_{{\text{C}}\max }} $ 和最大制冷率对应的COP$ {\eta _{\text{J}}} $ ; (b) 最佳充电能量$ {U_{{\text{MC}}}} $ 随温差$ \Delta T $ 的变化曲线Fig. 9. The curves of (a) the maximum cooling rate
$ {J_{{\text{C}}\max }} $ and the COP corresponding to maximum cooling rate$ {\eta _{\text{J}}} $ and (b) the optimal charging energy$ {U_{{\text{MC}}}} $ as a function of temperature difference$ \Delta T $ .图 10 在强耦合
$ {U_{{\text{HC}}}} $ 的情况下, 当$ \lambda = 0 $ 时, (a) 制冷率$ {J_{\text{C}}} $ 和(b)制冷系数$ {\eta _{{\text{COP}}}} $ 随着$ {U_{{\text{MH}}}} $ 和$ {U_{{\text{MC}}}} $ 变化的三维图Fig. 10. In the case of strong coupling
$ {U_{{\text{HC}}}} $ , when$ \lambda = 0 $ , the three-dimensional diagrams for (a) cooling rate$ {J_{\text{C}}} $ and (b) the COP varying with$ {U_{{\text{MH}}}} $ and$ {U_{{\text{MC}}}} $ . -
[1] Chen L, Ding Z, Sun F 2011 Energy 36 4011
Google Scholar
[2] Sánchez R, Büttiker M 2011 Phys. Rev. B 83 085428
Google Scholar
[3] Thierschmann H, Sánchez R, Sothmann B, Arnold F, Heyn C, Hansen W, Buhmann H, Molenkamp L W 2015 Nat. Nanotechnol 10 854
Google Scholar
[4] Zhang Y C, Lin G X, Chen J C 2015 Phys. Rev. E 91 052118
Google Scholar
[5] Zhang Y C, Huang C K, Lin G X, Chen J C 2015 Energy 85 200
Google Scholar
[6] Zhang Y C, Wang Y, Huang C K, Lin G X, Chen J C 2016 Energy 95 593
Google Scholar
[7] Aniket S. 2020 J. Appl. Phys 127 234903
Google Scholar
[8] Anamika B, Surojit H, Shailendra K, Varshney, Gourab D, Aniket S 2021 Phys. Rev. E 103 012131
Google Scholar
[9] Kano S, Fujii M 2017 Nanotechnology 28 095403
Google Scholar
[10] Lim J S, Sánchez D, López R 2018 N. J. Phys. 20 023038
Google Scholar
[11] Daré A M, Lombardo P 2017 Phys. Rev. B 96 115414
Google Scholar
[12] Su H, Shi Z C, He J Z 2015 Chin. Phys. Lett. 32 100501
Google Scholar
[13] 苏豪, 王家伟, 赵沁园, 何济洲 2016 中国科学: 技术科学 46 1296
Google Scholar
Su H, Wang J W, Zhao Q Y, He J Z 2016 Sci. Sin-Tech. 46 1296
Google Scholar
[14] Shi Z C, Qin W F, He J Z 2016 Mod. Phys. Lett. B 30 1650397
Google Scholar
[15] Roche B, Roulleau P, Jullien T, Jompol, Y, Farrer I, Ritchie D A, Glattli D C. 2015 Nat. Commun 6 6738
Google Scholar
[16] Hartmann F, Pfeffer P, Höfling S, Kamp M, Worschech L. 2015 Phys. Rev. Lett. 114 146805
Google Scholar
[17] Josefsson M, Svilans A, Burke A, Hoffmann E, Fahlvik S, Thelander C, Leijnse M, Linke H 2018 Nat. Nanotechnol 13 920
Google Scholar
[18] Keller A J, Lim J S, Sánchez D, López R, Amasha S, Katine J A 2016 Phys. Rev. Lett. 117 066602
Google Scholar
[19] Sothmann B, Sánchez R, Jordan A N, Büttiker M 2013 N. J. Phys. 15 095021
Google Scholar
[20] Choi Y, Jordan A N 2015 Physica E:Low-dimensional Systems Nanostruct 74 465
Google Scholar
[21] Lin Z B, Yang Y Y, Fu J, Li W, He J Z 2019 Chin. Phys. Lett. 36 060501
Google Scholar
[22] Lin Z B, Li W, Yang Y Y, He J Z 2020 Phys. Rev. E 101 022117
Google Scholar
[23] Sothmann B, Sánchez R, Jordan A N 2014 Europhys. Lett. 107 47003
Google Scholar
[24] Sánchez R, Sothmann B, Jordan A N 2015 Phys. Rev. Lett. 114 146801
Google Scholar
[25] Boukai A I, Bunimovich Y, Tahir-Kheli J, Yu J K, Goddard I W A, Heath J R 2008 Nature 451 168
Google Scholar
[26] Yang Y Y, Xu S, Li W, He J Z 2020 Phys. Scr. 95 095001
Google Scholar
[27] Yang Y Y, Xu S, He J Z 2020 Chin. Phys. Lett. 37 120502
Google Scholar
[28] Su S, Zhang Y, Chen J, Shih T M 2016 Sci. Reports 6 21425
Google Scholar
[29] Shi Z C, Fu J, Qin W F, He J Z 2017 Chin. Phys. Lett. 34 110501
Google Scholar
[30] Jiang J H, Entin-Wohlman O, Imry Y 2013 New Journal of Physics 15 075021
Google Scholar
[31] Li C, Zhang Y, He J 2013 Chin. Phys. Lett. 30 100501
Google Scholar
[32] Rutten B, Esposito M, Cleuren B 2009 Phys. Rev. B 80 235122
Google Scholar
[33] Cleuren B, Rutten B, Van den Broeck C 2012 Phys. Rev. Lett. 108 120603
Google Scholar
[34] 施志诚, 何济洲, 肖宇玲 2015 中国科学: 物理学 力学 天文学 45 050502
Google Scholar
Shi Z C, He J Z, Xiao Y L 2015 Sci. Sin-Phys. Mech. Astron. 45 050502
Google Scholar
[35] Li C, Zhang Y, Wang J, He J 2013 Phys. Rev. E 88 062120
Google Scholar
[36] Wang J H, Lai Y M, Ye Z L, He J Z, Ma Y L, Liang Q H 2015 Phys. Rev. E 91 050102
Google Scholar
[37] 李唯, 符婧, 杨贇贇, 何济洲 2019 物理学报 68 220501
Google Scholar
Li W, Fu J, Yang Y Y, He J Z 2019 Acta Phys. Sin. 68 220501
Google Scholar
[38] Whitney R S, Sánchez R, Haupt F, Splettstoesser J 2016 Phys. E 82 176
Google Scholar
[39] Fu T, Du J Y, Su S H, Su G Z, Chen J C 2021 Eur. Phys. J. Plus 136 1059
Google Scholar
[40] Xi M M, Wang R Q, Lu J C, Chen T Y, Jiang J H 2021 Chin. Phys. Lett 38 088801
Google Scholar
[41] 苏山河, 张艳超, 彭万里, 苏国珍, 陈金灿 2021 中国科学: 物理学 力学 天文学 51 112
Google Scholar
Su S H, Zhang Y C, Peng W L, Su G Z, Chen J C 2021 Sci. Sin-Phys. Mech. Astron. 51 112
Google Scholar
计量
- 文章访问数: 1244
- PDF下载量: 44
- 被引次数: 0