用于弹道声子动力学的离散实空间密度矩阵方法

《Journal of Computational Physics》：Discrete Real-Space Density Matrix Approach for Ballistic Phonon Dynamics

　　•提出了一种用于声子的离散实空间密度矩阵方法。•对非均匀晶格系统中的相干弹道输运进行了建模。•实空间传播避免了使用显式的位置-动量相空间网格。•带限能量算符提升了大规模模拟的可扩展性。

引言
声子输运决定了纳米尺度器件的热行为，对于高密度集成电路以及电子和声子异质结构而言愈发重要[1]、[11]、[12]、[13]。随着器件尺寸的缩小，界面效应、量子化效应和非局部效应变得越来越显著，这使得传统的宏观热传导模型无法用来描述热输运现象[11]、[12]。这一情况在薄膜系统中尤为突出，因为在这些系统中弹道输运和扩散输运机制共存，且热导率会受到边界散射、改变后的态密度以及量子化声子模式的影响[14]、[15]。

近年来，基于单粒子密度矩阵及其维格纳-魏尔变换来描述热输运，已被证明是一种极为有效的方法，能够全面捕捉固体中的相干和非局部输运现象。在这一领域具有里程碑意义的成果是Simoncelli[1]以及Marzari和Mauri的研究，他们提出了固体热输运的维格纳表述方法。他们证明，无论是来自佩尔斯-玻尔兹曼理论中的粒子状输运，还是对热传导有贡献的波动效应，都可以在统一的框架内得到描述[1]、[16]。此外，他们还指出，在复杂晶体以及低热导率材料中，声子描述的相干性及非对角元素起着至关重要的作用[1]、[16]。

Camiola、Romano和Vitanza[2]则采用了另一种互补的方法。他们从密度算符的量子力学刘维尔动力学出发，明确推导出了声子维格纳方程。这种方法清晰地展示了基于算符的微观动力学、维格纳变换与宏观可解释的输运量之间的联系[2]、[17]。这一点非常重要，因为它表明维格纳表示并非仅仅是一种替代符号，而是直接源于基于密度算符的描述[2]、[17]。

从数学上讲，晶格系统的维格纳函数可以被视为一种自然的动力学量。在这方面，Spohn从根本上阐明了晶格动力学、维格纳函数以及弱非谐系统的声子玻尔兹曼方程之间的关系[18]。具体而言，维格纳函数起到了连接微观波动动力学与动力学描述之间的桥梁作用[18]。

本文所提出的方法正是建立在这些基础研究之上的。不过，其主要目的并非重复阐述关于声子输运的通用维格纳理论，而是将其针对现代应用中一类特别重要的问题进行优化：集成电路中的空间非均匀异质结构。现有的主流文献大多关注周期性、类块体系统或大量利用对称性，而本文则侧重于离散的材料界面、局部质量变化、空间依赖的耦合关系以及有限几何结构。对于这类结构，需要一种能够在统一框架内整合实空间分辨率、量子力学相干性以及局部非均匀性的表示方法[1]、[5]、[15]。

虽然像格林函数方法这样的标准方法在处理静态谱传输和开接触问题时效率很高，但本文所采用的时域密度矩阵公式可以直接处理瞬态密度和空间相关性，并给出它们的位置-波矢表示。

基于此，以下的论述将从实空间中相互耦合的谐振子系统开始。采用质量加权公式、过渡到第二量子化，再通过单粒子密度矩阵进行描述，这不仅保证了形式上的一致性，更重要的是，能够将局部结构信息直接纳入哈密顿量中[3]、[5]、[13]。在这种方法中，材料界面或异质层序列自然而然地表现为能量算符的空间依赖矩阵元素，无需将其视为均匀参考系统的扰动。相比那些更偏向于倒易空间或连续介质模型的表示方法，这体现了该方法的主要优势。

对于谐振晶格中的纯相干初始状态，本文的密度矩阵公式与传统的谐振声子波包动力学是一致的，例如Schelling等人的方法[19]。第5.3节中的分析验证表明，该方法能够再现相应的色散关系、传播速度以及界面反射和透射现象。此外，本文公式还有一个优点，那就是能够直接传播完整的单粒子关联矩阵及其离散维格纳表示，从而无需进行全局模态分解，就能同时获取局部种群信息和非局部空间相关性。

随后的离散维格纳-魏尔变换将这种实空间描述转换为位置-动量相空间，同时不会丢失底层的原子级信息[4]、[5]、[11]、[15]。这对于异质结构尤为重要，因为在界面、薄膜内部或材料突变处存在的微观相干性可以直接决定输运行为[2]、[4]。因此，本文采用的方法旨在在已有文献的基础上，进一步发展上述的维格纳和密度矩阵方法，但重点在于适用于非均匀器件结构的离散实空间表述方法。对于这类特定的系统，目前可用且易于应用的表述方法相对较少。

因此，本文在现有主流文献中的定位如下：它基于这样一个核心观点，即相干声子输运最有效地可以通过密度矩阵和维格纳形式来描述，然后有针对性地将这一框架拓展为集成电路中异质结构所需的原子级离散描述。后续章节将在借鉴通用声子维格纳理论及其源自密度算符动力学的推导的基础上[1]、[2]，详细阐述一种能够系统处理局部非均匀性和空间离散耦合结构的表述方法。

片段摘录
离散实空间力常数模型
为描述声子的弹道输运，需在第二量子化框架下考虑晶格动力学。由于本研究聚焦于纯粹的弹道输运状态，该模型中忽略了非谐效应。通过明确考虑各个晶格点之间的相关性，可以灵活地表征空间非均匀性。诸如材料界面或原子质量的局部变化之类的空间差异，可以直接融入模型中。

能量算符的物理意义
对于均匀晶格，能量算符D具有平移不变性。因此，两个晶格点之间的耦合并不取决于它们的绝对位置，而只取决于它们之间的相对距离。假设晶格是无限的，那么在其离散表示中，D具有托普利茨矩阵的结构，即Dmn=dm?n。该矩阵的每一行都具有相同的耦合模式，对应于由相应晶格索引偏移后的模板。

算符D的作用
在位置-动量相空间中描述声子输运时，需要对单粒子密度矩阵进行离散维格纳-魏尔变换[1]、[2]、[4]、[5]。为此，首先引入质心坐标和相对坐标：χ=l+m2，ξ=l?m。这样，密度矩阵就可以表示为ρlm=ρ(χ,ξ)的形式。从密度矩阵的运动方程?tρlm=i∑kdklρkm?i∑ndmnρln出发，用坐标χ和ξ重新表示后，可得?tρ(χ,ξ)=i∑s[d(χ+12s,s?ξ)ρ(χ?ξ?

结果
通过一个模型系统，对用于描述声子弹道输运的所推导出的数值方法进行了验证和系统分析。测试系统为一个一维量子线，其总长度为L=64.63纳米，且在χ=45.08纳米处存在从硅到Si0.9Ge0.1的材料界面。这个数值测试系统是一个简化的一维模型，其中声子输运仅限于Si/SiGe纳米线的轴向。这种近似方式能够捕捉到主要的输运方向。

结论
本文开发了一种用于描述空间非均匀异质结构中相干弹道声子输运的实空间密度矩阵形式，并通过量子线中的动力学过程对其进行了验证。通过将空间依赖的耦合参数直接纳入离散哈密顿量中，该模型能够精确、原子级地呈现尖锐材料界面处的量子力学耦合和干涉过程。本研究的一个核心物理成果是

数据与代码获取
支持本研究结果的数据以及用于数值模拟的MATLAB代码，均可在TUDOData[28]网站上免费获取。

CRediT作者贡献说明
D. Mai：概念构建、方法论、软件实现、初稿撰写、可视化。D. Schulz：监督指导、方法论、验证、修改润色。

利益冲突声明
作者声明没有已知的财务利益或个人关系可能影响本文所述的工作。

D. Mai | D. Schulz