关于碘化铅中层间滑动驱动的有序-无序相变的研究报告

一、 研究团队与发表信息 本研究的主要作者包括Seyeong Cha、Giyeok Lee、Sol Lee、Sae Hee Ryu、Yeongsup Sohn、Gijeong An、Changmo Kang、Minsu Kim、Kwanpyo Kim、Aloysius Soon与Keun Su Kim。他们分别来自韩国延世大学物理系（Yonsei University, Department of Physics）、延世大学材料科学与工程系（Yonsei University, Department of Materials Science and Engineering）以及美国劳伦斯伯克利国家实验室先进光源（E. O. Lawrence Berkeley National Laboratory, Advanced Light Source）。该研究成果以题为“Order-disorder phase transition driven by interlayer sliding in lead iodides”的论文形式，于2023年发表在Nature Communications期刊上。

二、 学术背景与目的 本研究属于凝聚态物理、二维层状材料科学领域。二维层状材料因其在纳米电子学和光电子学方面的应用潜力而备受关注，同时也是研究多种基础量子现象（如相变）的平台。在碘化铅（PbI₂）这类材料中，已观察到多种相变，例如电荷密度波、磁性和结构相变。然而，对于其中一些相变的原子尺度微观机制，特别是层间堆垛和弱耦合在其中扮演的角色，尚不完全清晰。PbI₂本身作为一种重要的层状材料，被广泛用于合成钙钛矿太阳能电池的关键前驱体——卤化铅钙钛矿。

以往，PbI₂的低温晶体结构被普遍认为是中心对称的2H多型体。然而，最近的一些拉曼光谱和X射线衍射研究显示出一些令人困惑的结果，暗示可能存在其他多型体。此外，尽管卤化铅钙钛矿已知存在随温度变化的两步相变，但其二维对应物PbI₂中却未报道过此类相变。

因此，本研究旨在澄清PbI₂低温下的真实晶体结构，探究其是否以及存在何种相变，揭示其微观机制，并阐明该机制如何影响材料的光电与光学特性。具体目标是：1）精确确定PbI₂在低温下的晶体结构；2）发现并表征可能存在的相变；3）揭示相变的微观物理机制；4）分析相变对材料电子结构、光电导和光致发光性质的影响。

三、 详细研究流程 本研究结合了先进的材料表征技术与第一性原理计算，是一个多技术联用的系统性研究。主要流程可分为以下五个相互关联的部分：

1. 低温晶体结构确定：

   • 研究对象：高质量的块状PbI₂单晶（纯度99.9999%）。
   • 实验方法：首先，利用扫描透射电子显微镜（Scanning Transmission Electron Microscopy, STEM）对样品进行原子尺度成像。特别是采用高角度环形暗场成像（HAADF-STEM）技术，它可以提供原子序数对度（Z-对度）图像，从而清晰分辨出铅（Pb）原子和碘（I）原子的位置。STEM图像沿特定晶轴方向采集，以观察层内原子排列和层间堆垛顺序。
   • 同步分析：为了从电子结构层面验证晶体结构，研究团队采用了角分辨光电子能谱（Angle-Resolved Photoemission Spectroscopy, ARPES）。ARPES测量使用同步辐射光源（Advanced Light Source, beamline 7.0.2），光子能量范围为60-250 eV。样品在超高真空环境中原位解理以获得清洁表面。ARPES通过测量光电子的能量和动量，可以直接绘制材料的能带结构。对于PbI₂，其能带数目和色散特征对晶体结构的对称性（如是否为中心对称）和周期性（如原胞大小）非常敏感。
   • 理论计算：通过密度泛函理论（Density-Functional-Theory, DFT）计算了不同堆垛结构（如2H和4H）的电子能带结构。计算采用了VASP软件包，使用PBE泛函处理交换关联能，并考虑了范德华相互作用修正（DFT-D2）。为了更准确地描述能带，还使用了HSE06杂化泛函，并计入了自旋轨道耦合效应。

2. 相变的发现与表征：

   • 研究过程：在确定了低温结构后，研究团队系统地进行了变温ARPES实验。样品温度（Ts）从40 K变化到300 K。在多个温度点采集了沿高对称方向的ARPES能谱图。
   • 数据分析：为了定量分析电子结构随温度的变化，研究者从ARPES数据中提取了价带顶附近能量分布曲线（Energy-Distribution Curves, EDCs）的半高全宽（Full Width at Half Maximum, FWHM）。FWHM的变化反映了电子态寿命或能带展宽的程度，可以指示电子-声子耦合强度或无序度的变化。

3. 光电与光学性质关联测量：

   • 实验设计：在ARPES测量的同时，研究团队同步原位测量了两种关键物性：光电导（光电流）和光致发光（Photoluminescence, PL）。这是本实验设计的巧妙之处，确保了在不同温度下对同一区域、同一样品的电子结构、光电响应和光学发射进行关联性测量，避免了样品差异和测量条件不同带来的误差。
   • 测量方法：光电流通过测量样品与地之间的电流获得，光源为用于ARPES的同步辐射光。光致发光信号则通过探测样品受激后发出的绿色荧光强度来获取（使用96 eV光子激发）。

4. 相变微观机制的揭示：

   • TEM动态观测：为了直接观察相变前后原子排列的变化，研究团队分别在低于和高于相变温度（100 K和300 K）下进行了HAADF-STEM成像。通过对比两个温度下的STEM图像，可以直接观察Pb原子位置的分布变化。
   • 数据分析技术：为了更精确地分析原子位置，研究人员对STEM图像中的亮斑（对应Pb原子）进行了二维高斯函数拟合，从而提取出Pb原子的空间分布图。此外，还对STEM图像进行了傅里叶变换，通过分析倒易空间衍射点的形状（尖锐的点状 vs. 沿某个方向拉长的条纹状）来评估长程有序性。
   • 第一性原理计算验证：为了理解实验观察到的现象背后的物理原因，研究者构建了双层PbI₂模型，采用DFT计算了层间相对滑移（interlayer sliding）的能量势垒。他们系统地计算了上层相对于下层在不同滑动方向（如<110>和<100>方向）上位移时的能量变化（δE），绘制了能量变化图。

5. 数据整合与机制阐释：

   • 工作流程：将来自ARPES（电子结构展宽）、光电导（电流异常）、PL（发光强度斜率变化）、STEM（原子位置无序化）和DFT（滑动势垒各向异性）的所有数据和结果进行交叉比对和综合分析，构建了一个完整的逻辑链条，以阐明“层间滑动”是驱动有序-无序相变的根本机制。

四、 主要研究结果 1. 确认低温结构为4H多型体，而非传统认为的2H结构： * STEM结果：HAADF-STEM图像（图1d）清晰显示，相邻层中的Pb原子并非垂直堆叠，并且相邻层中的PbI₆八面体呈现镜像对称的交替倾斜（用红色和蓝色菱形表示），这与2H结构预期的一致性倾斜（灰色菱形）不符。这直接支持了4H多型体结构（图1c），其特点是具有镜面对称的八面体对沿堆叠方向交替排列，导致面外方向周期加倍，并且结构是非中心对称的。 * ARPES与理论验证：ARPES测量的能带结构（图1g, h）显示，在价带顶到-4 eV的能量范围内存在大量能带。DFT计算表明，2H-PbI₂在此能量区间应有6条价带，而4H-PbI₂由于原胞加倍（kz方向能带折叠）以及非中心对称性引起的Rashba型自旋劈裂，预计有多达24条价带。实验观测到的复杂能带结构与4H结构的预测高度一致。此外，kz方向的色散测量（图1i）显示周期为0.45 Å⁻¹，对应于4H结构的周期。非破坏性的XRD实验（补充信息）也支持这一结论。

1. 发现并表征了在120 K处的有序-无序相变：

   • ARPES证据：变温ARPES数据显示（图2a-d），在40 K到100 K之间，能谱仅有轻微的热展宽。然而，当温度升至200 K和300 K时，价带谱线出现显著的、突然的展宽。对FWHM的定量分析（图2e）表明，在40-120 K区间，FWHM的增加与电子-声子耦合导致的kBT量级展宽一致。但在温度高于120 K后，FWHM出现急剧增加，其能量尺度高达约50倍kBT，这无法用电子-声子耦合解释，而强烈暗示了无序度的突然增加，即发生了有序-无序相变。通过曲线拟合，确定其临界温度Tc约为120 K。作为对比，对已知具有有序-有序相变的1T-TiSe₂进行的类似分析显示，其FWHM的变化在Tc前后均大致遵循kBT趋势，从而凸显了本研究中相变类型的特殊性。

2. 相变伴随显著的光电与光学性质变化：

   • 光电导异常：同步测量的光电流随温度变化曲线（图2f）显示，在温度升至约160 K时，光电流出现一个明显的“凹陷”（dip）。这种在相变点附近出现的电导异常现象，通常归因于相变区域共存相和涨落导致的载流子散射增强或局域化。
   • 光致发光斜率变化：同步测量的PL强度（绿色发光）随温度升高而衰减，但其衰减速率在约130 K处发生突变（图2g中红线的斜率变化）。这表明在无序相中，电子和空穴的散射机制从主要由本征过程（如声子散射、辐射复合）主导，转变为由无序势等外部因素主导，从而抑制了辐射复合效率。

3. 揭示相变机制为层间滑动：

   • TEM直接观测：对比100 K和300 K的HAADF-STEM图像（图3a, b）发现，虽然镜像对称八面体对交替排列的特征在两个温度下都存在，但Pb原子的位置分布发生了变化。通过二维高斯拟合提取的Pb原子位置图（图3c, d）清晰显示，在100 K时，绝大多数Pb原子位于或接近4H结构预期的位置（红色圆圈）；而在300 K时，多个层中的Pb原子相对于4H结构的位置发生了横向偏移。
   • 傅里叶分析：对STEM图像进行傅里叶变换（图3e, f），100 K下的衍射图样呈现清晰、尖锐的点阵，对应良好的长程有序；而300 K下的衍射点则在z轴方向展宽为椭圆形（甚至在大面积选区衍射中呈现条纹），表明层间堆垛在面内方向出现无序，即层与层之间发生了不规则的相对滑动。
   • DFT计算阐明各向异性与低势垒：DFT计算的双层滑动能量图（图4b）揭示了层间滑动势垒具有明显的三度各向异性。沿着<110>方向的一个特定路径（图中红色箭头，从A位点滑动到C位点），滑动势垒极低，仅为10.6 meV/原子，这个值与相变临界温度Tc (~120 K, ~10.3 meV) 惊人地接近。而沿其他方向（蓝色和黄色箭头）的滑动势垒则高出2.4-3.4倍（图4c, d）。这完美解释了为什么在高温无序相中，层倾向于沿着这个低能量路径发生不规则滑动，同时也说明了滑动在面内是各向异性的。这种滑动机制不同于传统的“堆垛无序”，后者是指几种明确定义的堆垛结构（如立方和六方）的非周期性排列；而本研究中的滑动是沿着一个近乎连续的二维低能量路径发生的，更类似于二维熔化。

五、 结论与意义 本研究的结论是：碘化铅（PbI₂）在低温下的晶体结构是非中心对称的4H多型体，而非此前公认的中心对称2H多型体。该材料在约120 K处经历一个由层间滑动驱动的有序-无序相变。这一相变在电子结构上表现为价带谱线的急剧展宽，并伴随光电导异常和光致发光效率变化率的突变。层间滑动具有显著的面内各向异性，其最低能量势垒仅为10.6 meV/原子，与相变临界温度相对应。

本研究的科学价值在于： 1. 机制创新：明确地将“层间滑动”确立为二维范德华层状材料中一种关键的、此前可能被忽视的相变驱动机制。这为理解其他具有弱层间耦合的层状材料中的结构相变和物性变化提供了新的视角。 2. 关联性突破：成功地将原子尺度的结构变化（通过STEM观测）、电子结构演化（通过ARPES）、宏观光电与光学性质（通过同步光电导和PL测量）以及理论计算（DFT势垒）紧密联系起来，完整地揭示了从微观机制到宏观物性的因果链条。 3. 修正认知：修正了关于PbI₂低温晶体结构的长期认知，并首次在PbI₂中报告了这种独特的相变，丰富了对其物相行为的理解。

应用价值方面，PbI₂作为制备高效钙钛矿太阳能电池的关键前驱体材料，其相变和相关的层间滑动行为可能会影响钙钛矿薄膜的结晶过程、界面特性以及最终器件的性能与稳定性。理解这一基础相变机制，可能为优化钙钛矿材料的制备工艺和性能提供新的思路。

六、 研究亮点 1. 多技术深度融合与同步测量：本研究最突出的亮点是将ARPES、STEM、光电导、PL和DFT计算多种尖端技术有机结合，并创新性地在ARPES测量过程中同步原位获取光电导和PL信号，确保了数据在时间和空间上的一致性，极大地增强了论证的说服力。 2. 原子尺度直接观测与机制确认：利用高分辨率STEM在相变温度上下直接“看到”了Pb原子位置的分布从有序变为无序，并通过傅里叶分析证实了层间滑动导致的堆垛无序，为相变机制提供了最直接的实验证据。 3. 低能量势垒的精确计算与完美关联：DFT计算不仅预测了层间滑动势垒的面内各向异性，更重要的是计算出的最低滑动势垒（10.6 meV/原子）与实验测得的相变临界温度在能量尺度上高度吻合，为“热涨落克服滑动势垒引发相变”的物理图像提供了坚实的理论支撑。 4. 清晰的物理图像：研究构建了一个非常清晰的物理图像：低温下，层状PbI₂以有序的4H方式堆叠；温度升高至Tc时，热涨落足以克服沿特定<110>方向的极低层间滑动势垒，导致层开始不规则滑动，系统进入无序相；这种原子位置的无序化引起电子态展宽（ARPES结果），进而导致载流子散射机制改变（光电导异常）和辐射复合效率变化（PL斜率改变）。

七、 其他有价值的发现 1. 非中心对称性的确认：确定了低温4H相的非中心对称性，这可能导致新的物理效应，如可能存在的压电性或非线性光学响应，为未来探索PbI₂的其他功能特性提供了线索。 2. 与1T-TMDs的对比与联系：研究过程中与1T-TiSe₂（一种典型的过渡金属硫族化合物）的有序-有序相变进行了对比，突出了本研究发现的相变类型（有序-无序）的特殊性。同时，研究也指出层间滑动机制在理解其他层状材料（如1T-TaS₂）相变中的潜在重要性，建立了更广泛的联系。