揭秘磁分子超快退相干机制：自旋-声子相互作用驱动的四阶纯相消过程

《SCIENCE ADVANCES》：A unified ab initio theory of spin-phonon relaxation and decoherence uncovers fast dephasing in magnetic molecules

为了阐明自旋退相干的微观起源，并解决经典磁驰豫理论无法解释实验观测到的自旋相干时间（T2）与自旋弛豫时间（T1）之间的显著差距，研究人员发展了首个统一描述自旋-声子弛豫与退相干过程的第一性原理理论。该工作将量子主方程扩展到密度矩阵的相干项，并完整包含了高达双声子过程的贡献。研究以具有大磁阻塞温度的单分子磁体 DyCp 为模型体系进行计算，揭示出尽管其强大的轴向磁各向异性确保了在 77K 下长达秒量级的缓慢磁弛豫，但由于双声子纯相消机制，其Kramers双重态叠加态的相干时间却短于 10 纳秒。这一过程原则上适用于任何与声子热浴相互作用的量子系统，极大地增进了我们对固态系统中量子退相干的理解。

　　

在量子科学的前沿，一个核心挑战是如何维持量子系统的“相干性”——即其能够同时处于多种状态的“叠加态”特性。这种脆弱的状态是量子计算、精密传感等颠覆性技术的基石。然而，量子系统不可避免地会与环境发生相互作用，导致其相干性迅速衰减，这个过程被称为“退相干”。在固态自旋系统中，例如单分子磁体（SMMs），自旋与环境（主要是晶格振动，即声子）的耦合已知会驱动自旋能量弛豫（即改变自旋的能级布居数），但长期以来，声子对“纯相消”（即不伴随能量交换的相位随机化）的贡献一直被忽视或低估。实验上经常观察到，自旋的相干时间（T₂）远短于由能量弛豫时间（T₁）所设定的理论极限的两倍（2T₁），即使在设法消除了自旋-自旋相互作用等噪声源后，这一差距依然存在。这暗示着，可能有一种未知的、高效的退相干机制在起作用，而现有的理论框架无法对其进行描述。为了填补这一知识空白，并从根本上理解固态量子系统中退相干的微观物理图像，一项结合了理论发展与第一性原理计算的研究在《SCIENCE ADVANCES》上报道了其突破性进展。

该研究主要运用了扩展的量子主方程理论和第一性原理计算相结合的方法。研究人员首先从理论上推导了包含密度矩阵非对角元（相干项）演化的四阶量子主方程，构建了一个能够统一描述自旋-声子弛豫与退相干的微扰理论框架。随后，他们将该理论完全第一性原理地应用于一个具体的、具有记录性长磁弛豫时间的镝基单分子磁体 [DyCp₂^ttt]⁺（简称 DyCp）。计算所需的全部输入参数，包括分子的自旋哈密顿量、晶格声子谱以及自旋-声子耦合系数，均通过结合密度泛函理论和多参考态量子化学方法从头计算获得。

研究结果

量子主方程

为获取全温度范围内的相干时间，研究团队推导了一个四阶生成子，能够描述整个自旋密度矩阵的时间演化。他们将系统哈密顿量分解为自旋、声子及二者相互作用三部分。在弱耦合和马尔可夫近似下，利用正则化T矩阵方法，得到了决定密度矩阵演化的二阶和四阶超算符的显式表达式。这些表达式以广义Lindblad形式写出，并引入了跳跃算符来描述声子浴对自旋的影响。其中，四阶跳跃算符囊括了所有可能影响自旋动力学的双声子过程（双吸收、双发射、吸收-发射各一）。这一理论框架的建立，使得完整模拟自旋相干动力学成为可能。

数值模拟

他们将上述方法数值实现，用于计算DyCp的自旋动力学。DyCp具有J = 15/2的基态和大磁各向异性。模拟准确再现了该分子在整个温度范围内的实验磁弛豫数据。在高温区（>60 K），磁弛豫由二阶理论主导，表现为Orbach机制（通过逐级吸收/发射单声子翻越磁翻转能垒）。在低温区，四阶理论贡献的双声子Raman弛豫机制开始主导。重要的是，模拟显示，尽管该分子的磁化翻转时间在77 K时接近秒量级，但其基态Kramers双重态（KD）量子叠加态的相干时间T₂却异常短暂，在相同温度下短于10纳秒，比磁弛豫时间快了八个数量级。

弛豫与退相干对退相干的贡献

为了理解T₂如此之短的原因，他们分析了其对T₁（弛豫贡献）和T₂^*（纯相消贡献）的分解。结果表明，在二阶水平上，T₂完全由T₁决定，因为单声子过程无法导致零能量的纯相消。然而，在四阶水平上，T₂^*成为限制T₂的主要因素。其根源在于，决定四阶纯相消的跳跃算矩阵元涉及 ΔM_J= ±1 的虚拟跃迁，这类跃迁在具有大磁各向异性的DyCp中矩阵元很大，因此效率极高。相反，决定磁弛豫的四阶矩阵元涉及跨越整个磁翻转能垒的虚拟跃迁，其矩阵元被分子的设计刻意压制，因而过程缓慢。这从机制上解释了为何相干性比磁化强度“翻转”得如此之快。

相干时间的声子分解

进一步分析声子贡献发现，在低温下，退相干主要由能量在25-50 cm^-1的低频声子驱动。随着温度升高，能量在~175-200 cm^-1的声子被热激活，贡献变得同等重要。这些声子模式主要对应于分子中叔丁基（^tBu）基团的振动。计算得到的四阶T₂随温度变化的幂律指数约为2.76，接近但不完全等于标准双声子Raman过程的高温极限（T²），反映了参与过程的声子具有特定的能量分布。

研究结论与讨论

本工作通过将时间局域主方程理论扩展至四阶，首次实现了对自旋弛豫与退相干过程的统一第一性原理描述。该理论揭示了一个在四阶微扰下出现的、纯声子性质的“纯相消”机制。在以DyCp为代表的高磁各向异性单分子磁体中，这一机制极为高效，尽管磁弛豫被工程化地减慢，但它仍能导致超快的量子退相干，使相干时间比磁化翻转时间短八个数量级。这一过程预计具有普适性，原则上适用于任何与玻色子热库耦合的量子系统。

该研究为理解和调控分子量子比特的相干性提供了关键见解。为了减缓退相干，未来策略应聚焦于：1）使用自旋-1/2体系或有机自由基，它们具有高能激发和低自旋-轨道耦合；2）将分子嵌入金属/共价有机框架或刚性笼中，以“硬化”导致退相分的低频分子振动（如叔丁基的扭动）；3）探索高频动力学解耦或利用多能级结构设计对声子退相干具有鲁棒性的量子态。尽管在DyCp中直接测量T₂具有挑战，但本理论对所有输入参数的预测成功复现了磁弛豫实验数据，间接验证了其可靠性。近年来，在高温、动力学解耦实验乃至钟态跃迁处，T₂与2T₁之间的偏差被屡次观测到，本工作所揭示的声子纯相消机制为解释这些现象提供了一个普适的理论候选。总之，这项研究不仅增进了对固态量子退相干微观起源的理解，也为未来设计更长相干时间的量子材料提供了理论蓝图和明确指导。