1. 项目概述当量子涨落遇上应变工程在凝聚态物理和材料科学领域铁电材料一直是一个充满魅力的研究方向。想象一下一块晶体内部的正负电荷中心在不受外电场作用时也能自发地朝一个方向“排队”形成永久的电偶极矩这就是铁电性。这种特性让它们在非易失性存储器、压电传感器、能量收集器乃至未来的神经形态计算器件中扮演着核心角色。然而并非所有具有铁电潜质的材料都能在低温下稳定地展现出这种有序态。有一类特殊的材料被称为量子顺电体它们被强烈的量子涨落“冻结”在了无序的边缘钽酸钾KTaO3就是其中的典型代表。KTaO3在低温下本应发生铁电相变但其晶格中原子核的量子零点运动即量子涨落异常强烈这种无处不在的“微观颤抖”足以破坏任何试图建立起来的铁电有序使其始终保持在一种“欲铁电而未成”的量子顺电态。这就好比试图在波涛汹涌的海面上立起一根静止的旗杆量子涨落就是那永不平静的海浪。传统的第一性原理计算如密度泛函理论DFT在处理这种涉及量子核效应和强非谐性的问题时往往力不从心因为它本质上是零温下的基态理论无法捕捉这些动态涨落。因此长久以来如何精确描述并最终调控KTaO3中的量子涨落与铁电有序的竞争是一个理论上的难题也制约着基于此类材料的功能器件设计。近年来两个关键技术的融合为破解这一难题提供了新思路。一是随机自洽谐波近似SSCHA它是一种能够非微扰地处理晶格非谐性并同时包含核量子效应与热效应的强大框架。二是机器学习力场MLFF它能够从少量的第一性原理计算数据中学习到高精度的原子间相互作用势从而将SSCHA这类需要大量能量和力评估的计算成本降低数个数量级。我们的工作正是将SSCHA与MLFF相结合构建了一个高效且精确的计算工作流系统探究了面内应变如何作为一种“调控旋钮”去影响KTaO3中量子涨落与铁电有序之间的微妙平衡。我们发现即便是微小的应变0.2%-1%也足以打破对称性显著提升铁电相变的能垒让量子涨落“望而却步”从而成功诱导出可观的铁电极化其强度甚至能与经典的铁电体钛酸钡BaTiO3媲美。这项研究不仅深化了对量子顺电体物理的理解更为通过“应变工程”这一实验上可实施的手段设计新型量子功能材料提供了清晰的理论路线图。2. 核心物理图像与计算方法解析2.1 量子顺电体与量子临界点要理解KTaO3首先要建立量子顺电体的物理图像。在经典相变理论中温度降低会抑制热涨落使系统有序化例如铁磁体在居里温度以下出现磁序。对于铁电体其相变通常与一个光学软模的频率软化相关——当温度趋近相变点时某个特定晶格振动模式的频率会趋于零意味着晶格失稳原子发生集体偏移形成自发极化。然而在KTaO3这样的量子顺电体中即便温度降至绝对零度热涨落消失量子涨落依然存在。这些源于海森堡不确定性原理的量子零点运动其能量尺度与铁电相变的能垒相当甚至更大。结果就是软模频率在低温下虽然持续下降但永远不会达到零而是在某个有限值饱和系统被“卡”在有序和无序的边界上。这个状态对应的参数空间位置就被称为量子临界点QCP。处于QCP附近的材料对外界扰动极其敏感微小的压力、电场或——正如我们重点研究的——应变都可能成为压倒骆驼的最后一根稻草将系统推入铁电相。这为操控材料性质提供了一个极具潜力的杠杆。2.2 方法论基石SSCHA与MLFF的强强联合面对量子涨落和非谐性这两个难题我们采用了SSCHAMLFF的组合拳。下面我详细拆解一下这个工作流的设计逻辑和实操要点。2.2.1 随机自洽谐波近似SSCHA超越简谐近似传统的晶格动力学计算基于简谐近似将原子势能面在平衡位置作二阶展开此时振动模式是独立的简谐振子。这对于KTaO3是严重不足的因为它会错误地预测出一个虚频负的振动频率平方意味着结构失稳而实际上量子涨落和非谐相互作用稳定了它。SSCHA的核心思想是自洽地寻找一个最优的“有效谐波”哈密顿量使其在考虑真实非谐势能面的情况下给出的原子分布与真实量子/热分布最接近。具体步骤是一个循环迭代过程假设一个试探性的有效力常数即有效谐波势。基于此有效势通过蒙特卡洛或随机采样方法生成包含量子与热涨落的原子构型集合。对这些采样构型用第一性原理方法计算每个原子所受的真实力。比较真实力与有效谐波势预测的力根据差值更新有效力常数使其预测的力在统计平均意义上最接近真实力。重复步骤2-4直至有效力常数收敛。这个过程能同时处理非谐性和量子效应给出重整化后的声子谱虚频被消除和自由能面。但它的瓶颈在于第3步对每个应变、每个温度条件都需要对超胞中数百个原子构型进行昂贵的DFT力计算计算量巨大。2.2.2 机器学习力场MLFF从数据中学习的“代理模型”这正是MLFF大显身手的地方。我们采用的是一种“在线学习on-the-fly”的MLFF方案。它的工作流程如下初始化与主动学习从一个小的初始构型数据集开始用DFT计算其能量和原子受力。模型训练使用贝叶斯回归方法如高斯过程或神经网络训练一个力场模型。该模型以原子局域环境描述符如原子间距离、角度等为输入预测能量和原子力。不确定性估计MLFF模型的关键优势在于能提供预测的不确定性。在SSCHA的采样过程中对于模型预测不确定性高的新原子构型系统会自动调用DFT进行精确计算并将这个新数据点加入训练集。迭代优化用更新后的数据集重新训练MLFF模型然后继续SSCHA采样。如此循环MLFF模型在“实战”中不断进化其预测范围逐渐覆盖整个相空间感兴趣的区域。最终绝大多数构型的力和能量由快速比DFT快数万倍的MLFF提供只有少数不确定性高的点才需要DFT计算。这样我们在保持第一性原理精度的前提下将SSCHA的计算成本降低了数个数量级使得系统研究应变-温度相图成为可能。实操心得MLFF训练的数据策略在构建MLFF时初始数据集的代表性至关重要。对于KTaO3这类涉及结构相变的体系我们不仅采集了平衡位置的构型还特意通过分子动力学模拟或在软模方向施加位移采样了非平衡态和可能过渡态的构型。这确保了MLFF能准确描述势能面在铁电畸变路径上的形状这是计算能垒和极化的基础。忽略这一点MLFF可能在关键区域预测失准。2.3 计算细节与参数选择我们的第一性原理计算基于VASP软件包采用SCANStrongly Constrained and Appropriately Normed泛函。这是一个meta-GGA级别的泛函相比传统的PBE泛函它能更准确地描述电子局域化效应和范德瓦尔斯相互作用对于钙钛矿氧化物这类强关联电子体系的关键性质如晶格常数、带隙、声子预测通常更可靠。平面波截断能设置为800 eVk点网格密度经过严格测试以确保总能量和力的收敛。所有计算在固定实验晶格常数3.9842 Å的基础上进行只松弛原子位置以隔离应变效应避免引入额外的软模。极化计算采用贝里相位Berry Phase方法这是从第一性原理计算铁电极化矢量的标准方法。它通过计算电子波函数在倒易空间中的几何相位贝里联络来定义极化避免了传统基于离子电荷和位移的模型可能带来的歧义。3. 应变诱导的结构相变与对称性破缺3.1 应变施加与对称性演化我们研究了两种面内应变模式单轴应变沿一个晶轴方向拉伸或压缩垂直方向自由弛豫和双轴应变在两个面内方向施加等量应变。应变范围从0%到1%这个量级在实验上通过衬底失配或机械应力是完全可以实现的。未应变本征的KTaO3在低温下属于中心对称的立方相空间群为Pm-3m或其在计算中常用的P4/mmm二者等价。施加应变后晶格对称性被打破单轴应变立方对称性降低为四方相空间群变为P4mm。这个群是非中心对称的允许沿应变方向c轴产生净极化具备了铁电性的结构基础。双轴应变对称性进一步降低为正交相空间群变为Amm2。此时面内的两个方向a轴和b轴不再等价极化可以出现在面内的某个方向。这种从高对称性到低对称性的转变是应变诱导铁电性的结构序曲。它移除了反演对称性为原子偏离高对称位置、形成电偶极矩扫清了障碍。3.2 键长变化揭示晶格畸变的微观指纹为了定量表征畸变我们固定钾K原子的位置因其质量大、位移小重点分析氧O原子与钽Ta原子之间键长的变化。在钙钛矿ABO3结构中铁电畸变主要表现为B位阳离子Ta相对于氧八面体中心的偏移。我们计算了O1-Ta和O2-Ta键长分别对应a和b方向随温度和应变的变化。结果清晰地显示温度效应对于任一固定应变键长变化量即畸变幅度随温度升高而单调减小。在0 K时畸变最大随着温度升高热涨落增强破坏了部分有序畸变减弱直至在居里温度Tc以上键长变化趋于零系统恢复高对称性P4/mmm。这直观地描绘了从铁电相到顺电相的连续相变过程。应变效应在固定温度下尤其是低温键长变化量随应变增大而单调增加。1%应变下观察到的畸变最显著。这说明应变不仅开启了铁电相变的“开关”还像一个“调幅旋钮”直接控制着极化强度的大小。单轴 vs. 双轴有趣的是在相同应变幅度下单轴和双轴应变诱导的键长变化幅度几乎相同。对于双轴应变a和b方向的畸变也近乎相等。这表明对于KTaO3应变诱导的极化强度主要取决于应变的幅度而非其各向异性程度在研究的范围内。极化方向则由对称性决定单轴应变下沿c轴双轴应变下在ab面内。注意事项理解“固定晶格”的计算设定在本文的计算中我们固定了晶格常数只松弛原子位置。这是一种常用的简化目的是专注于研究“纯”应变效应和离子位移带来的铁电性。在实际薄膜或应变器件中衬底约束会导致面内晶格常数固定相干应变但面外晶格常数可以弛豫。我们的计算设定对应于“面内受夹面外自由”的常见薄膜情况。若考虑完全的体积弛豫定量结果可能会有细微变化但定性趋势应变诱导铁电是稳健的。4. 声子谱、能垒与铁电极化的深度分析4.1 非谐性稳定软模从虚频到实频图4展示了在1%应变、0 K条件下包含与非包含非谐效应的声子谱对比。在传统的简谐近似黑色虚线计算中在布里渊区中心Γ点附近仍然存在明显的虚频分支图中通常表现为负频率或绘图时显示为“下沉”到零以下的模式。这个软模对应着铁电畸变的振动模式虚频意味着在简谐近似下这个畸变方向仍然是失稳的。然而当引入SSCHA方法包含非谐效应后红色实线这个软模频率被重整化到了正的真实频率。其他声子分支则几乎不受影响。这一现象至关重要物理意义它证明了在应变KTaO3中非谐相互作用原子间势能的高阶项与量子涨落共同作用提供了额外的恢复力稳定了晶格阻止了它塌缩到简谐理论预测的另一个能量更低的结构。系统找到了一个量子非谐效应下的新平衡点。计算启示这强烈警示我们对于量子顺电体及其衍生物忽略非谐性会严重误判其动力学稳定性、相变温度乃至是否存在铁电性。简谐近似在此类体系中基本失效。4.2 铁电相变能垒量子涨落与应变的博弈能垒是理解相变的关键。我们将系统从顺电高对称相所有原子在中心位置沿软模方向位移到铁电低对称相计算了这条路径上的能量变化从而得到了铁电相变的能垒。本征KTaO3文献值和我们的计算确认其能垒极低约0.3 meV/原子。这个能量尺度与量子零点运动的能量相当因此量子涨落可以轻易地“跨越”这个能垒使系统在不同极化方向间隧穿宏观上表现为无净极化的量子顺电态。1%单轴应变下的KTaO3我们的计算显示能垒显著提升至~3.1 meV/原子。作为对比典型铁电体BaTiO3的能垒约为4 meV/原子。MLFF验证使用我们训练的MLFF计算的能垒为3.2 meV/原子与DFT直接计算结果高度一致交叉验证了MLFF在描述复杂势能面方面的惊人精度。这个能垒的提升是应变诱导铁电性的核心物理机制。3 meV/原子的能垒已经超出了量子涨落典型的活动范围。量子涨落不再能轻易抹平不同极化态之间的能量差异因此系统可以“定居”在其中一个铁电态中形成长程铁电有序。应变的作用实质上是抬高了铁电相的“能量洼地”或者说压低了顺电相的“能量高原”使得两者之间的能垒相对量子涨落的能量尺度变高了。4.3 极化计算与相图构建基于收敛的SSCHA结构我们使用贝里相位方法计算了每个应变-温度组合下的自发极化强度。极化强度在0 K和1%应变下取得了最大值。单轴应变时极化沿c轴达到0.233 C/m²双轴应变时极化在ab面内沿对角线方向合成为0.293 C/m²。这个数值与室温下BaTiO3的单畴自发极化强度约0.26 C/m²处于同一量级表明应变KTaO3具有实用级别的铁电极化强度。应变-温度相图我们将计算结果整理成相图如图5a示意。相图清晰地显示了一个铁电相和顺电相共存的区域其边界即居里温度线Tc(ε)。随着应变ε增大Tc向高温方向移动。例如在0.5%应变下Tc可能在几十K而在1%应变下Tc可提升至百K量级。这为实验上通过控制应变和温度来操控铁电态提供了精确的“地图”。软模频率的温度演化如图5b所示在铁电相内T Tc软模频率随温度升高而增加。这与传统铁电体如BaTiO3的行为一致温度升高热涨落增强倾向于破坏有序使得恢复力常数增大即原子被“拉”回高对称位置的趋势变强表现为软模频率“硬化”。当T接近Tc时频率仍保持有限正值在Tc处发生相变之后系统进入顺电相。5. 实操启示、潜在挑战与拓展方向5.1 对实验与器件设计的启示应变窗口明确理论预测表明低至0.2%的面内拉伸应变即可在KTaO3中诱导出铁电性。这为实验制备指明了可行的应变范围。现代薄膜外延技术如在不同晶格常数的衬底上生长完全可以实现1%甚至更高的应变。材料选择与优化本研究聚焦纯KTaO3。在实际中掺杂如用Nb部分替代Ta形成KTa1-xNbxO3是调节其量子顺电性的常用手段。我们的方法可以无缝扩展到掺杂体系预测掺杂与应变的协同效应寻找更高的Tc和更大的极化。多功能器件潜力KTaO3是著名的二维电子气2DEG宿主材料其界面具有巨大的Rashba自旋轨道耦合效应。应变诱导的铁电性可以与这些电子性质耦合实现电场控制的自旋输运多铁性耦合或制造非易失性、低功耗的铁电晶体管存储器。5.2 计算实践中的常见问题与排查MLFF训练不收敛或精度差问题力场对某些构型的预测误差大导致SSCHA迭代振荡或给出非物理结果。排查首先检查初始训练集是否覆盖了足够的相空间。对于铁电体系务必包含沿软模方向的一系列位移构型。其次检查描述符描述原子环境的方法是否合适。对于钙钛矿氧化物通常采用原子中心对称函数或矩张量描述符效果较好。最后监控MLFF预测的贝叶斯误差确保在SSCHA采样过程中对高误差构型的DFT回调是有效的。SSCHA计算中声子谱仍有虚频问题即使包含了非谐性在某些q点或某些模式下仍出现小的虚频。排查这可能是由于以下原因a) 超胞尺寸不够大无法准确描述长程相互作用或某些特定的非谐耦合。尝试增大超胞。b) SSCHA迭代未完全收敛。增加蒙特卡洛采样步数降低自洽循环的能量收敛阈值。c) 在非常接近相变点时如T≈Tc涨落极大可能需要更精细的计算或考虑其他序参量的耦合。贝里相位极化计算出现量子化错误问题极化值出现大的、不连续的跳变通常是极化量子的整数倍。排查贝里相位极化依赖于所选的规范其差值才有物理意义。确保计算铁电相和参考顺电相时使用了相同的k点路径和网格。计算时需要沿着铁电畸变的路径选取一系列中间图像结构计算每个结构的极化观察其连续变化。如果直接从顺电相跳到铁电相很可能得到错误的极化量子。5.3 未来拓展方向考虑应变梯度与挠曲电效应在实际纳米结构或畴壁中应变可能是不均匀的存在应变梯度。应变梯度可以直接诱导极化即挠曲电效应。将当前方法与应变梯度建模结合可以研究更复杂的纳米尺度铁电拓扑结构如涡旋、斯格明子。耦合电子-声子计算当前工作主要关注晶格动力学。下一步可以结合GW近似或含时DFT研究应变和铁电极化如何影响KTaO3的电子结构、载流子迁移率以及可能出现的拓扑电子态。高通量筛选与逆向设计将本工作流自动化可以对钙钛矿氧化物家族甚至更广泛的材料进行高通量计算筛选出那些处于量子顺电边缘、可通过微小应变诱导出强铁电性的候选材料实现逆向材料设计。动力学过程模拟利用训练好的MLFF可以进行大规模分子动力学模拟研究铁电畴壁的运动、极化翻转动力学等非平衡过程这对理解器件操作速度至关重要。这项研究通过将前沿的机器学习力场与先进的非谐晶格动力学方法相结合不仅解决了KTaO3这一具体材料体系中的关键物理问题更展示了一套强大的计算范式。这套范式对于研究其他受量子涨落、非谐性支配的复杂功能材料如超导体、热电材料、量子自旋液体等具有普适的参考价值。它让我们看到在“计算显微镜”的视野里我们已能够以前所未有的精度观测并设计量子世界与宏观性能之间的桥梁。