李冀 陈亮 冯芒2)‡

1) (广州工业技术研究院,广州 511458)

2) (中国科学院精密测量科学与技术创新研究院,武汉 430071)

1 引言

热传导是物理学中的一个基本现象,它描述了当物质受到加热或冷却时,热量在物质内部以及物质之间传递的机制.对热传导的研究有助于深入了解宏观现象和微观动力学之间的联系.例如,一维模型中的热传导作为一个典型的研究对象,可以帮助我们理解热量传输过程中的非线性和混沌等复杂现象[1].对热传导的深入研究可以揭示热量传输的内在机理和数学模型,为实际应用提供理论依据和指导.

在经典系统中,热传导遵循傅里叶定律,该定律阐述了热流密度与空间中温度场之间的相互关系.假设在三维空间中存在一个宏观系统,在某一特定时刻t,系统空间的温度分布可用T(q,t) 来描述.其中,向量q(x,y,z) 表示系统内的某个空间位置,x,y和z分别表示三维坐标轴上的坐标.根据傅里叶定律,在此情况下,热流密度J(q,t) 可以表示为[2]

其中κ为系统的导热系数.为了理解温度场在时间维度上的变化规律,这里定义系统在三维空间内的能量分布为u(q,t) .将方程(1)代入连续性方程∂u/∂t+∇J=0[3]得

其中c=∂u/∂T是系统的比热容.方程(2)也被称为热扩散方程,它描述了宏观系统中能量转移的规律.傅里叶定律是从大量宏观系统实验中总结出的经验定律,其基本假设是系统中的粒子在进行热量交换时主要依靠的输运过程是扩散.也就是说,粒子通过与周围粒子的随机碰撞来进行能量的交换.从声子的角度来看,这可以理解为数量多的声子向数量少的声子的方向扩散.然而,当系统的尺寸和维度减小到一定程度后,粒子之间的热交换过程可能不再符合扩散输运的规律.例如,热交换过程可能会变为弹道输运[4],此时粒子在进行能量交换时没有发生碰撞,并且可以忽略声子之间的相互作用.在这种情况下,热传导不再符合傅里叶定律.

为了深入了解低维系统中热传导的特性,科学家们针对两种简单模型进行了研究.其中一种模型是谐波晶体,另一种则是FPU (Fermi-Pasta-Ulam)晶体.在这两种模型中,假定晶体包含有限颗粒子,并且晶体两端与两个不同温度且温度恒定的热库相接触.这些粒子在空间中以固定点为中心进行简谐振动,同时它们仅仅与相邻的粒子进行相互作用.当两个粒子之间的相互作用势表示为时,其为FPU 晶体,其中β为非谐参数;i,l表示两个粒子的编号;qi,ql表示粒子在空间上的位置.而当时,称为谐波晶体.综述文献[3,5]详尽地介绍了相关研究,展示了当前的研究现状和未来的发展趋势.人们通常以系统的稳态热流J∼Na作为评判热传导特性的重要标准.其中,N代表模型中的粒子数,而a是一个常数.在傅里叶定律的适用范围内,系统的热流应满足a=–1[6,7].对于这两种晶体,当系统处于三维状态时,其稳态热流遵循傅里叶定律.然而,当系统维度降低至二维或一维时,参数a不再等于1,这表明热传导特性发生了变化.对于一维谐波晶体,当系统处于有序状态,即所有粒子的质量或粒子之间相互作用力等参数均相同时,其稳态热流与系统大小N无关,此时a=0[4,8].给系统中的粒子的质量引入无序因素后,例如系统中粒子质量并不全相等,a的值可能为–3/2[9,10]或者为–1/2[11,12].针对FPU晶体,人们发现了一个有趣的特性: 当FPU 晶体的温度低于某个阈值时,其热传导特性会与谐波晶体极为相似[13].

离子晶体是由囚禁在离子阱中的多个离子组成的晶体,本综述主要探讨离子晶体中的热传导特性.这类晶体因其独特的性质而备受关注: 在真空环境中,囚禁的离子晶体与外部环境没有能量交换,其晶体结构和温度可以通过电场和光场进行精确操控[14–16].这为研究低维晶体在经典或量子状态下的热传导提供了理想的实验平台.在离子晶体中,离子之间的距离相对于其他晶体材料要大得多,这使得离子间的相互作用主要是库仑力.这种情况下,离子会与晶体中所有其他离子进行能量交换.通过调节系统的参数,可以使离子晶体呈现出不同维度的构型,例如直线构型(一维)、之字构型(二维)或螺旋构型(三维)[17],这些不同维度构型的离子晶体为研究热传导提供了丰富的模型.

本综述的结构安排如下: 第2 节详细阐述热传导研究中常用的一些基本原理和方法;第3 节分别采用一维、二维和三维模型来计算离子晶体中的热传导,并展示离子晶体的结构相变与温度分布以及与热流的关系;第4 节在三维离子晶体热传导模型中加入无序度的情况下,探讨离子晶体中热传导的特殊表现和规律;第5 节介绍一些离子晶体中与热传导相关的实验,并总结目前在离子晶体中对热传导研究的主要成果.

2 离子晶体热传导模型

本文的研究对象是两端与不同温度的热库接触的离子晶体模型,期望获得该系统在稳态情况下的热流密度以及温度分布.离子晶体中的热传导研究通常借助郎之万方程来描述系统中每个离子的运动.此方程原用于描述布朗运动,其中摩擦力和随机噪声项在热传导情境中可对应于热库对离子的作用,可以证明,此系统的稳态是唯一的[18].通过对稳态下离子晶体的深入分析,可以得到热流密度、温度分布等关键信息.

该模型涉及两类离子,其中一类离子与热库接触,而另一类离子则不与热库接触.为了计算模型的稳态,有两种常用的方法可供选择: RLL (Rieder Lebowitz Lieb)法和LEGF (Langevin equations and Green’s function)法.RLL 法是由Rieder,Lebowitz 和Lieb 三人首次提出的[4],LEGF 法是由Ford,Kac 和Mazur 三人首次提出[19],并且在后续工作[2,20–23]中得到了发展.这两种方法的主要区别在于如何处理这两类离子.LEGF 法通过哈密顿量进行区分,将总哈密顿量分别按照这两类离子写出各自独立的哈密顿量.相对而言,RLL 法是在运动方程的层面对两者进行区分,即根据不同的离子采用不同的运动方程.

下面基于RLL 方法介绍一些基本原理.首先定义模型中离子的数量为N,离子的质量为m,那么离子晶体的哈密顿量可以写为

其中pi表示第i个离子的动量,V(qi) 表示离子在阱势场中的势能,而U(qi-ql) 则表示离子i与离子l之间的相互作用势能.此时,可以将离子受到的来自阱势场和其他离子的力表示为fi=-∂H/∂qi.假设第1 个离子与高温热库L 接触,第N个离子与低温热库R 相接触.那么系统的运动方程为

其中,γL,R表示热库L 或R 对离子的摩擦系数,ηL,R(t)表示热库L 或R 对离子的随机力.接下来,需要用离散形式定义热流,首先定义由第i个离子流向第l个离子的瞬态的热流密度为ji→l,其内涵是第i个离子向第l个离子做功的功率,即:

定义热库L,R 对第1,N个离子的作用力为fL,fR,那么热库流向离子的热流可以表示为

为了利用连续性方程,定义第i个离子上的能量密度为[3,24]

于是,连续性方程的离散形式可以表示为[3,24]

当处于稳态的时候离子晶体的热流不再变化,即∂ε(qi)/∂t=0 ,使用〈···〉表示物理量的平均值,可得此时的总热流就是离子晶体与热库之间的热流,使用〈J〉代表系统稳态下的总热流,可得

根据温度与能量的关系

3 离子晶体中的热传导

离子阱系统中,产生与离子晶体接触的热库有两种方式.离子晶体被离子阱势场囚禁在真空中,其与外界交换能量的方式通常有两种: 离子阱势场的作用和激光的作用.在与离子阱势场的作用中,通过调整阱势场的频率、幅度以及离子在势场中的位置,可以对离子晶体产生不同程度的加热作用.在与激光的相互作用中,通过调整激光的功率和频率,可以对离子晶体产生不同程度的加热和冷却效果.因此,在实验中可以通过这两种方式在离子晶体的两边产生不同温度的热库: 1)改变离子晶体两端最边缘的若干个离子所处的囚禁势场,这样会使这两类离子受到来自阱的不同的加热作用,此时的囚禁势场相当于两个不同温度的热库;2)让两端最边缘的若干个离子分别受到两束不同参数的激光照射,此时激光就能被视作一个等效的热库.因为使用激光产生热库的方案更容易实现,在接下来的章节中,讨论的模型都以激光作为热库.本节将详细介绍并探讨文献[25–28]针对不同维度模型中的离子晶体热流的研究成果.

3.1 一维模型

考虑一维模型,N个离子沿x轴线性排列,由于离子阱势场的限制,每个离子均以固定位置为中心做简谐振动,用ωx表示离子的宏运动频率.设表示第i个离子简谐振动的中心坐标,而xi和pi表示其相对于运动中心的相对位移和动量.根据海森伯-郎之万方程(Heisenberg-Langevin equation)[29,30],系统的运动方程为

其中ρi表示热库对离子的摩擦系数,ζi(t) 表示热库对离子的随机力.Ail表示离子之间相互作用力的 耦合矩阵,其表达式为[25]

假设离子晶体两端的热库为Markovian 型热库,那么ζi(t) 满足以下条件( ℏ=1):

xi(t) 和pi(t) 为第i个离子在t时刻的位移和动量,I为单位矩阵,Ω为一个大小为 2N×2N的矩阵.

在时刻t,第i个离子的温度可以表示为Ti(t)=当时间t→∞时,可以得到离子的稳态温度=Ti(t→∞) .

在文献[25]中,采用了上述方法对包含100 个离子的模型在不同ρ的取值下的温度分布进行了计算.当两个热库位于离子晶体边缘时,离子晶体中的温度分布如图1 所示.显然,图中的温度分布不符合傅里叶定律.对于符合傅里叶定律的系统,其温度分布应该呈现从高温点逐渐下降到低温点的线性分布,然而在图1 的分布中,与热库接触的离子和未与热库接触的离子间的温度存在显著的跃阶现象,且所有未与热库接触的离子的温度分布几乎完全一致.这种温度分布的出现,可以归因于热库将热量传递到离子晶体的速率远低于离子晶体内部热传导的速度.

图1 ρ 值不同时,由100 个离子组成的离子晶体在稳态下的温度分布[25],其中离子的温度 由平均声子数来表示.第1 个离子和第100 个离子与热库相接触,由图中的红线表示,其他的离子由蓝线表示Fig.1.Temperature distribution of an ion crystal composed of 100 ions under different ρ in stable state,where temperature distribution measured by the mean thermal phonon number.The 1st and the 100th ions are in contact with the thermal bath,indicated by the red line in the figure,while other ions are indicated by the blue line.The figure is taken from Ref.[25].

在另一种情况下,如图2 所示,将高温热库移到离子阱晶体的中央,此时的温度分布将会分为两部分.位于两个热库之间的部分将维持在一个相对较低的温度,而位于高温热库另一侧的部分则展现出较高的温度.值得注意的是,这两个温度之间的差距会随着热导率的降低而逐渐减小并最终趋于一致.从上述计算结果可以看出,相较于符合傅里叶分布的宏观系统,一维离子晶体的温度稳态分布具有多种不同寻常的模式.

图2 由 101 个离子组成的离子晶体在稳态下的温度分布[25],其中离子的温度 由平均声子数来表示.第1 个离子和第51个离子与热库相接触.与热库相接触的离子由图中的红线表示,其他的离子由蓝线表示Fig.2.Temperature distribution of an ion crystal composed of 101 ions under different ρ in stable state,where temperature distribution measured by the mean thermal phonon number.The first ion and the 51st ion are in contact with the thermal bath.Ions in contact with the thermal bath are indicated by the red line in the figure,while other ions are indicated by the blue line.The figure is taken from Ref.[25].

3.2 二维模型

接下来,我们将探讨更为复杂的二维模型.在二维模型中,离子晶体能够发生结构相变,从一维的直线构型转变为二维的之字构型.将模型的两个维度分别定义为x轴和y轴,离子在这两个轴向上的宏运动频率分别为ωx和ωy.定义该模型中的离子数量为N,离子的质量为m,所带电荷为Q.另外定义qi和pi为离子相对坐标零点的位移和动量,qx,i,qy,i和px,i,py,i分别为它 们在两个坐标轴上 的投影.于是系统的哈密顿量可以写为[26]

可以看出,该模型中的离子振动中心位置并不固定,因此当选取不同的离子阱参数时,离子晶体将呈现出不同的构型.在二维模型中,离子的运动方程可以写为[26]

其中µ表示x或y;ηµ,i和Dµ,i表示来自热库的摩擦系数和扩散系数;Wµ,i与热库对离子的随机力εµ,i(t) 有关,并且随机力满足条件〈εµ,i(t)〉=0 和〈εµ,i(t)εµ,i(t′)〉=2Dµ,iδ(t-t′).热库是由激光对离子的作用产生的,在弱光强的条件下可得[26]

假设激光只与x轴上投影坐标最大和最小的几个离子相互作用,分别产生两个不同温度的热库,而其他离子则不与激光发生相互作用.用U(|qi-ql|) 代表离子之间的相互作用势,V(qi) 代表离子阱的势场,离子的能量密度可以表示为[3,24]

定义从第i个离子流向第l个离子的热流密度为

激光与第i个离子作用,使这个离子与热库产生了接触,热库流向第i个离子的热流密度为

定义热流的方向是由编号较小的离子流向编号较大的离子.对于单个离子来说,其自身的能量变化等于流向离子的能量密度减去流出离子的能量密度.因此,第i个离子的连续性方程的离散形式可以写成:

在稳态情况下,流入离子的能量与流出离子的能量相等,因此由(24)式,稳态下的连续性方程可以改写成:

式中的〈jB→i〉可以由Novikov 定律[31]得到

为了得到通过离子晶体的总热流J(t),需要把单个离子的连续性方程扩展成整个离子晶体的连续性方程[3,24]:

其中ε(q,t) 代表离子晶体的能量密度;j(q,t) 代表离子晶体的热流密度;ji(t) 代表经过第i个离子的热流密度,其表达式可以通过傅里叶分析得到[26]

将热流密度j(q,t) 在整个离子晶体范围内进行积分,得到通过离子晶体的总热流J(t) :

利用Novikov 定理写出稳态下离子晶体总热流的表达式[31]:

稳态时离子晶体中的温度分布可以表示为

文献[26]采用数值模拟方法计算了包含N=30 个离子的离子晶体在不同构型下的稳态热流和温度分布,其中编号为i=1,2,3 和i=N– 2,N– 1,N的离子分别受到相同功率不同失谐量的两束激光照射,而其他离子不与激光相互作用.离子晶体的结构相变过程用序参量α来表示,其定义见文献[26].在文献[26]采用的数值模拟条件下,α=13 时,离子晶体为直线构型,随着α逐渐减小,离子晶体相变为之字构型.在离子晶体为直线构型时,离子沿着x轴排列.在直线构型到之字构型的结构相变过程中,离子晶体的温度分布如图3 所示,其中温度分布梯度逐渐增大.在直线构型中,离子晶体温度分布规律与一维模型计算结果一致,数值模拟结果显示离子晶体中间部分的温度为两个热库的平均值.在之字构型中,离子晶体内温度分布更接近符合傅里叶定律的宏观系统.

图3 处于不同构型下的离子晶体在x 方向上的温度分布 [26],其中 α 表示结构相变的序参量.两个颜色不同的区域中的离子和不同参数的激光相互作用.插图表示离子晶体的构型Fig.3.Temperature distribution of ion crystals vs.the distance x in different configurations,where α represents the order parameter of structural phase transition,and the two regions with different colors represent ions interacting with the laser.The insets show the configurations of ion crystals.The figure is taken from Ref.[26].

在二维模型中,离子的振动模式从一维(x方向)拓展到了二维(x和y方向).这两个方向的宏运动频率越接近,离子晶体的构型就越接近之字构型,反之,则更接近直线构型.图4 给出了离子晶体在两种构型之间的相变过程中热流的变化趋势.随着α值的不断增大,离子晶体从之字构型相变到直线构型.在这个过程中,热流逐渐增加.当α=13时,离子晶体完成了从之字构型到直线构型的完全转变,热流达到最大值.然而,随着α值的继续增大,热流量开始逐渐下降,但是下降速度较为缓慢.这种现象出现的原因在于离子晶体从之字构型向直线构型相变的过程中,离子之间的距离减小,导致离子之间的相互作用力增强,这使得离子晶体中的热流增大.另外,在相变点附近,离子热运动产生的噪声增大,这种噪声能够减小离子晶体中的非平衡性影响[32],也使得热流增强.

图4 处于不同构型下的离子晶体的稳态热流 [26],其中α表示结构相变的序参量.插图表示离子晶体的构型Fig.4.Heat flux of ion crystals in different configurations,where α represents the order parameter of structural phase transition.The insets show the configurations of ion crystals.The figure is taken from Ref.[26].

3.3 三维模型

在三维模型中,将研究对象扩展到螺旋构型,为此,需要首先计算出在三维情况下离子晶体内的稳态热流和温度分布的表达式.定义离子晶体中一共有N个离子,且质量都为m,则系统的哈密顿量可以写为

其中q=(qx,i,qy,i,qz,i) 是三维空间中的坐标,p=(px,i,py,i,pz,i) 是三维空间中的动量.定义µ=(x,y,z)代表坐标轴上的3 个方向.两束激光分别与离子晶体两端的离子相互作用,使离子晶体与两个温度不同的热库L 和R 相接触.定义离子晶体中的第1 个到第NL个离子与热库L 相接触,第NR个到第N个离子与热库R 相接触,其他的离子不与热库接触.此时,离子的运动方程可以写为

接下来计算系统的稳态热流,首先定义系统的能量密度为

根据(34)式和(39)式,可以推出系统的热流密度j(q,t) 和热库流向离子的热流密度jB→i分别为[28]

其中G(λ) 为满足以下条件的任意函数[34]:

利用(39)式、(41)式和(42)式,可以写出连续性方程:

通过对j(q,t) 积分,可得系统的热流:

在稳态下,系统的平均热流与热库流向系统的热流相等,于是稳态热流可以写为

其中k代表所有与热库接触的离子的编号,Bk代表与编号为k的离子接触的热库.考虑Novikov定律[31],利用(34)式和(37)式,稳态热流的表达式可化为

文献[28]用上述方法计算了不同构型下离子晶体的热流.在三维模型中,由于运动模式增加了一维,离子晶体的构型需要用另一种方式表达.假设离子晶体处于直线构型时沿着x方向,用ωx,y,z表示离子在3 个方向的宏运动频率,并定义ωy=nyωx,ωz=nzωx,那么可以用 (ny,nz) 表示离子的构型.图5 显示了在文献[28]的计算参数下,ny和nz取值不同时离子晶体结构相变的相图.

图5 在不同 ny 和 nz 的取值下,离子晶体结构相变的相图[28].用L (linear),Z (zigzag)和H (helical)这三个字母代表离子晶体的直线构型、之字构型和螺旋构型.沿HL 线,(ny,nz) 值的变化会产生直线构型和螺旋构型之间的结构相变.沿ZL 线,(ny,nz)值的变化会产生直线构型与之字构型之间的结构相变.沿ZHZ 线,(ny,nz) 值的变化会产生从之字构型到螺旋构型再到之字构型的结构相变.图中的灰色区域代表直线构型.在灰色区域之外,HL 线之上的区域代表处于 x-y 平面上的之字构型,在HL 线之下的区域代表处于x -z 平面上的之字构型,HL 线上为螺旋构型Fig.5.Phase diagram of ion crystal structural phase transitions under different values of ny and nz .L (linear),Z (zigzag),and H(helical) represent the linear,zigzag,and helical configurations of ion crystal.ZL line,HL line and ZHZ line represent the zigzag-linear transition,helical-linear transition and zigzag-helical-zigzag transition.The gray region in the figure represents the linear configuration.Outside of the gray region,the region above the HL line represents the zigzag configuration lying in the x-y plane,while the region below the HL line represents the zigzag configuration lying in the x-z plane.The region on the HL line represents the helical configuration.The figure is taken from Ref.[28].

离子晶体在不同构型下的温度分布如图6(a)所示,其中直线构型和之字构型的温度分布规律与二维模型的计算结果相似.螺旋构型的温度分布相比之字构型更接近宏观状态下符合傅里叶定律的温度分布.从图6(b)可以看出,螺旋构型的离子晶体中温度梯度基本保持不变,而不像之字构型中的温度梯度不断变化.

图6 在离子数 N=30 的离子晶体中,不同构型的离子晶体的稳态温度分布[28],图中L 代表直线构型、Z 代表之字构型、H 代表螺旋构型,图中的参数为 (ny,nz) (a) 3 种不同构型的离子晶体在稳态下的温度分布,两个颜色不同的区域中的离子和不同参数的激光相互作用;(b)在稳态下,3 种不同构型的离子晶体不与激光接触的部分的温度梯度分布,其中 ∆T 代表与系统中心的距离Fig.6.Temperature distribution of ion crystals with 30 ions with different configurations in stable state,where L represents linear configuration,Z represents zigzag configuration,and H represents helical configuration.The labels indicate the values of (ny,nz) :(a) Temperature distribution of three different configurations of ion crystals in stable state,where two different colors represent ions interacting with the laser;(b) distribution of temperature gradients in the parts of ion crystals with three different configurations that are not in contact with the laser in stable state,where ∆T represents the distance from the center of the system.The figure is taken from Ref.[28].

在保持ny=nz的前提下,离子晶体可以直接从直线构型相变到螺旋构型,而这两个构型之间的相变过程中的温度分布和热流与直线构型和之字构型之间的相变相似.图7 给出了当 (ny,nz) 的值沿着图5 中的ZHZ 线变化时离子晶体的温度分布和热流的变化情况.从图7 可以看出,离子晶体无论是处于x-z平面的之字构型还是处于x-y平面的之字构型,其热流均小于处于螺旋构型时的热流.此外,还发现ny和nz的差值越大,热流越小.这是由于y和z两个维度之间的不平衡导致的非线性效应阻碍了热流的传导[28].

图7 在离子数 N=30 的离子晶体中,离子晶体从处于x-z 平面的之字构型相变到螺旋构型再相变到处于x-y 平面的之字构型的过程中的温度分布和热流变化[28] (a)离子晶体在结构相变过程中的3 种典型的温度分布,两个颜色不同的区域中的离子和不同参数的激光相互作用,图中Z 代表之字构型,H 代表螺旋构型,图中的参数为 (ny,nz) ;(b)在 ny 保存不变时,nz 不断增大的过程中的热流变化Fig.7.Temperature distribution and heat flux changes in an ion crystal with 30 ions during the process of zigzag-helical-zigzag transition where Z represents zigzag configuration and H represents helical configuration: (a) Three typical temperature distributions of the ion crystal during the structural phase transition process,where two different colors represent ions interacting with the laser;the labels indicate the values of (ny,nz) ;(b) changes in heat flux during the process of increasing nz while keepingny constant.The figure is taken from Ref.[28].

4 热流与无序度之间的关系

在研究谐波晶体的热流中,无序度是一个重要的因素,通常无序度的变化会导致热流的变化.本节将讨论一种能够解析计算离子晶体稳态热流的方法,并使用该方法在离子晶体中引入无序度.同时,分析无序度对离子晶体中热流的影响.

首先,需要写出系统的哈密顿量H,其表达式为H=HS+HE+Hint,由3 个部分组成.其中,结构哈密顿量HS和环境哈密顿量HE分别表示不与热库接触和与热库接触的两类离子的哈密顿量,相互作用哈密顿量Hint表示前面这两类离子之间相互作用的哈密顿量[27,35].定义系统中有N个不与热库接触的离子,X=(q1,···,q3N)T表示这些离子在x,y,z这3 个方向上的坐标,其中T 表示矩阵转置.P=(p1,···,p3N)T表示离子在3 个方向的动量,用V表示离子之间的耦合矩阵.于是系统的结构哈密顿量可以写为

另一方面,用 (e) 表示与离子接触的热库,每个热库与N(e) 个离子接触.假设离子晶体的两边分别与热库L 和R 接触,那么 (e)=(L) 或 (R) .于是系统的环境哈密顿量可以写为

根据文献[35],热库的性质由ξ和ν决定,表达式分别为

假设热库的类型为Ohmic 热库,此时

其中Λ表示高频截止频率,ξ0表示弛豫率,P(e) 表示与热库(e)接触的离子在坐标轴上的投影.Gˆ(s)是系统格林函数(Green’s function)的拉普拉斯变换(Laplace transform),其表达式为

通过对(54)式和(55)式进行积分,就可以得到稳态下的热流等物理量.文献[27,36,37]提供了一种方法,能够解析地计算积分.首先计算(s)的6N个特征值{sφ}和对应的特征向量{rφ}.根据文献[38],(s) 可以被写成:

并且满足以下条件:

于是(54)式的积分可以被写为[27]

其中A=2kB(T(L)P(L)+T(R)P(R)) ,并且ωφ=-isφ.

(55)式的积分可以被写为[27]

其中∆≠-2ikB(T(L)-T(R)) .

文献[27]采用上述方法计算了离子晶体在不同构型下的热导率κL与无序度间的关系,其中κL的定义为〈J〉=κL(T(L)-T(R))/L,L为离子晶体的长度.上述模型虽然是三维,但是也能计算出一维(直线构型)和二维(之字构型)情况下离子晶体受到无序度的影响.研究者们在耦合矩阵V中引入无序度,具体操作如下: 随机选取离子晶体中的一半离子,将其耦合矩阵对应值进行改变Vii→(1±d)Vii,其中d为表示无序度的参数.在文献[27]中,根据数值计算得到了以下结论: 无序度和离子晶体的长度均会对离子晶体的热导率产生影响.对于一维情况(直线构型),当无序度d=0.005 且离子晶体长度较小时,热导率基本不受无序度的影响.然而,当晶体长度继续增加,热导率出现了明显的下降.对于无序度d=0.02 的情况,热导率受到无序度影响的离子晶体的长度的阈值相对于d=0.005 减小了约一半,并且当晶体长度足够大时,热导率下降到了接近0.对于二维情况(之字构型),当d=0.001 时,热导率就明显受到了无序度的影响,随着离子晶体长度的增加,热导率会下降到接近0.当d的值提高到0.02 时,热导率会在更小的晶体长度处下降到接近0.对于三维情况(螺旋构型),热导率受到无序度的影响也很明显,但是热导率没有出现下降到0 的趋势.

除了离子晶体,研究者也在其他低维系统中对无序度进行了研究.文献[39]在一维FPU 晶体中探讨了无序度与热流之间的关系,他们通过给晶体中的粒子引入不同质量来产生无序度.数值模拟结果显示,在高温或大尺寸的FPU晶体中,无序度对热流的影响较弱.只有在晶体尺寸较小且温度较低的情况下,无序度才会对热流产生显著影响,此时无序度越高,热流越小.文献[6]在二维和三维的谐波晶体中考察了无序度与热流的关系,他们同样通过给粒子赋予不同的质量来产生无序度.对于二维谐波晶体,热流受到无序度和晶体尺寸的共同影响.无序度或晶体尺寸越大,热流越小.当无序度和晶体尺寸达到一定值时,热流会降低至接近零.然而,对于三维谐波晶体,热流受无序度和晶体尺寸的影响较小.尽管无序度和晶体尺寸的增加会使热流减小,但热流并未呈现出趋于零的趋势.

5 总结

本文简要介绍了离子晶体中热传导的研究进展,包括一维、二维和三维模型中温度分布和稳态热流的计算方法,以及在不同维度的离子晶体构型下热流与温度分布的特性.此外,还介绍了无序度对离子晶体热导性的影响.在没有引入无序的情况下,直线构型的离子晶体具有最高的导热性,其次是螺旋构型,最后是之字构型.产生这一现象的原因是: 在离子晶体中运动模式不平衡导致的非线性效应会阻碍热流的传导.直线构型具有最少的运动模式(仅为一维),因此其导热性最佳.之字构型和螺旋构型都包含3 个运动模式,它们之间的区别在于之字构型的3 个运动模式均不相等,而螺旋构型中有两个维度的运动模式是相等的.因此,之字构型的导热性要小于螺旋构型.

在引入无序的情况下,无序会导致离子晶体的热导性下降,其中之字构型受到的影响最大.当离子晶体中的离子数量较少时(25 个以下),直线构型受到影响最小;而当离子数量较多时(120 个以上),螺旋构型受到的影响最小.如果将离子晶体中轴向运动的不平衡也视为一种无序,那么这也能够解释离子晶体中之字构型的导热率最低.

实验上已经成功构造出了离子晶体的多种构型,比如直线构型、之字构型[40]、螺旋构型[41]和带拓扑结的之字构型[42].实验上也实现了离子晶体在不同构型之间的结构相变[43],而且观察到温度对离子晶体构型的影响[44,45].但是到目前为止,还没有实验对上述理论描述的稳态下离子晶体的热传导进行直接验证,主要原因在于: 目前在离子阱系统中,用于测量离子温度的方法(扫描离子荧光线型或测量离子边带谱线等)都不适用于稳态下的离子晶体.比如,在使用扫描离子荧光线型的方法时,需要使用一束不断改变失谐的激光持续照射离子,这相当于离子晶体与一个额外的热库相接触并且破坏了稳态.而扫描离子边带谱线的方法虽然不会破坏稳态,但是只能准确测量运动能量低于20 个声子数的离子体系,相对来说,稳态下离子体系至少会有几百个声子数.目前最有希望对稳态进行无破坏测量的技术是通过使用CCD 拍摄离子照片来识别离子的温度[46].在离子阱系统中已经有研究者在热传导方向进行了相关的实验,比如在文献[47,48]中,用脉冲激光加热离子晶体中的一个离子,然后测量了热量在离子晶体中的传播(扩散、反射等)规律.

目前对有关低维系统热传导的理论已经在实验室中有所应用,比如在离子阱中的协同冷却技术,这种技术能够在只对离子晶体中的一个离子进行冷却的情况下,降低离子晶体中所有离子的温度[49,50].对于离子晶体,存在两种非常有效的冷却方法,分别是EIT (electromagnetically-induced-transparency)冷却[51]和偏振梯度冷却[52].对离子晶体热传导的深入研究有助于更好地理解这几种冷却方法的冷却过程以及所能达到的极限冷却温度.在离子晶体以外的低维系统中,为了测量热流,研究者们发展出了多种实验技术.文献[53]新发展出一种加工技术,能够制造尺寸小于100 nm 的独立三维结构,通过应用这个技术,制造了绝缘热导体和热传感器,并测量了纳米线的热导率.文献[54]发明了一种悬挂装置,这种装置能够用于探测无基底接触的纳米管中的热流.

在宏观系统中,根据变换热学的理论,科学家们已经成功研发出能够控制热量流动方向的超构材料,例如热二极管.同样地,在低维系统中也能够制造出类似功能的微观热二极管.文献[55]在实验中成功使用纳米管制造出了微米级别的热二极管,整流参数约为7%.而文献[56]在理论上模拟了由15 个离子构成的离子晶体组成的热二极管,他们让离子晶体处于直线构型,并调整每个离子所在的势阱,使势阱的频率依次递增.将处于此状态的离子晶体的左右两端与热库接触,理论计算结果显示,当左边的热库温度高于右边的热库时的热流与左边的热库温度低于右边热库时的热流存在明显的差异.理论上,这种热二极管的整流参数最高可达80%.