不同晶体结构的金属中应力调控的结构与电子演变 胡书山、黄刚

《Crystals》：Stress-Modulated Structural and Electronic Evolution in Metals with Different Crystal Structures Shushan Hu and Gang Huang

　　**摘要**
尽管宏观应力显著影响金属组件的性能，但其背后的原子-电子耦合机制仍未能得到充分理解，这些机制受不同晶体对称性的支配。为填补这一空白，本研究系统地分析了在施加应力作用下代表性金属材料的结构和电子演变。实验中，X射线衍射（XRD）揭示了冷轧钛合金和硅钢中复杂的宏观残余应力分布。受这些工程观察结果的启发，研究人员采用第一性原理密度泛函理论（DFT）计算来揭示其背后的物理机制。具体而言，研究了面心立方（FCC）铝和铜、体心立方（BCC）铁以及六方密排（HCP）钛晶体在拉伸和压缩条件下的响应。计算了应力依赖的弹性性质、态密度（DOS）、能带结构和声子谱。结果表明，拉伸会软化所有金属（铝在拉伸过程中变得机械不稳定），而压缩则会硬化它们的晶格。从电子层面看，拉伸载荷会使费米能级附近的DOS峰变尖锐，并使导带向其移动；而压缩则会使DOS峰变平缓，导带远离费米能级。声子分析表明，铜和钛在高压下仍保持动态稳定性，而铝和铁则表现出声子模式软化现象。这些由应力引起的变化揭示了关键的原子-电子耦合机制，为通过应力工程调整金属性能提供了理论基础。

**1. 引言**
金属材料因其高强度、延展性以及优异的热电性能而在航空航天、汽车、能源和生物医学领域得到广泛应用[1,2,3]。常见的金属晶体结构包括面心立方（FCC）、体心立方（BCC）和六方密排（HCP），每种结构都具有不同的机械特性。例如，基于FCC的高熵合金（如CoCrFeMnNi）表现出卓越的延展性和韧性[4]，而基于Fe、Mo和Nb的BCC难熔合金在常温和高温下均具有高强度和硬度[5]。HCP金属（如Ti、Zr和Mg）以其高比强度、低密度和耐腐蚀性而著称，非常适合用于航空航天、核能和电子行业[6]。因此，这些金属的微观结构配置对其宏观机械性能和实际应用至关重要[7,8]。

当金属晶体受到外部应力作用时，晶格参数、键长和原子位置会发生变化，从而直接影响其电子结构。近期实验证实了应力对金属电子态的显著影响：施加约0.5%的压缩应变会导致杨氏模量突然下降约10%，在更高压力下再次增加约20%，这归因于新的电子态进入费米面[9]。这凸显了应力调节的电子占据与机械性能之间的紧密联系。正如理论预期的那样，成键和反键电子态决定了金属的强度和凝聚性；因此，改变电子结构可以有效地调节强度和延展性。Borges等人通过第一性原理计算表明，成键态的占据和d轨道DOS的双峰与BCC高熵合金中的位错核心结构和韧性密切相关[10]。此外，电子能级的重新分布可能引发相变或机械不稳定性[11]。总体而言，应力会改变键长、键角和晶格对称性，进而影响费米面附近的能级和DOS。最近的理论研究越来越关注通过可控应变工程来调节多种晶体系统的相稳定性以及电子和光学响应[12,13]。这种原子-电子耦合对于理解应力作用下的金属性能演变至关重要[14,15]。

不同晶体结构对应力的弹性和电子响应存在显著差异。第一性原理计算表明，在应变作用下弹性常数和声子谱会发生变化。例如，Zhang等人结合密度泛函微扰理论（DFPT）和实验数据发现，高达2 GPa的压力会导致铝和铜的声子谱蓝移，减弱电子-声子耦合并降低电阻率[16]。这证明了通过应力调节电子和声子性质的可行性。应力作用下杨氏模量和泊松比等弹性模量的变化也引起了关注[17]。例如，Ho等人的DFT研究在特定应变条件下发现，FCC材料（如铝和铜）沿主轴方向会出现负泊松比，显示出弹性性质的复杂方向依赖性[18]。此外，Singh等人发现费米面的拓扑变化会导致弹性常数的突然软化[19]。同样，不同相的能带结构和DOS显著决定了它们的宏观机械响应。例如，You等人研究了三种Mg-Zn合金相的电子结构和声子特性，观察到键合稳定性的显著差异：立方相几乎各向同性，而六方相则表现出明显的各向异性[20]。这些研究共同强调了应力诱导的改变对带宽、DOS分布和振动模式的影响，所有这些都影响材料的机械响应。

尽管取得了这些进展，但仍存在关键挑战。大多数计算研究关注的是在均匀应变下的理想晶体，而实际材料通常包含缺陷和微观结构异质性。此外，现有的理论模型（如自由电子模型或硬球模型）无法完全捕捉应力诱导的电子-晶格耦合[21]。在实验方面，能够同时测量应力变化和电子结构变化的技术（如费米面演变）仍然有限。此外，大多数研究单独探讨应力对弹性、能带结构或声子行为的影响。因此，亟需一种综合的多尺度方法，能够同时揭示材料在施加机械载荷下的原子、电子和振动响应。

为解决这些问题，本研究结合实验和计算方法，系统地分析了金属的结构和电子演变。这里采用了第一性原理密度泛函理论（DFT），该理论在预测复杂系统的本征电子和磁性质方面已被证明非常可靠[22]。本研究不仅将FCC、BCC和HCP金属视为孤立案例进行研究，还试图建立一个将宏观应力与本征晶格对称性联系起来的通用物理基础。实验中对多相钛合金和硅钢样品进行了宏观XRD残余应力测量，以揭示实际工程材料中的复杂应力状态。然而，由于晶界、沉淀物和次级相的存在，直接从这些多相合金中提取本征的结构-性能关系较为复杂。因此，为了分离这些宏观现象并揭示仅由基本晶体对称性控制的基本物理机制，随后对典型的理想单晶（FCC铝和铜、BCC铁以及HCP钛）在拉伸和压缩条件下进行了第一性原理模拟。计算了应力依赖的晶格常数、杨氏模量、泊松比以及相应的能带结构、DOS和声子谱的变化。这种统一的方法将电子和机械性质联系起来，明确区分了纯几何体积/应变效应与本征电子重构（如轨道杂化），为通过应力工程优化金属性能提供了理论基础。

**2. 实验和计算方法**
**2.1. 实验程序**
本研究中使用的冷轧钛合金和硅钢板由首钢集团（北京）提供，并被切割成矩形样品。样品表面应力测量前先使用120、200、500、1000和2400目的SiC砂纸进行机械研磨，然后使用2.5 μm钻石悬浮液进行精细抛光。接着使用5%高氯酸-乙醇溶液（分析纯度>99%）在30 V和50 mA的条件下进行10秒的电解抛光，以去除机械研磨造成的损伤表面层。样品在无水乙醇中超声清洗并干燥，以获得光滑的镜面表面。

此外，还对冷轧钛合金和硅钢样品进行了残余应力测量，每个表面上预定了九个测量点。使用配备Cr Kα辐射源（用于硅钢）和Cu Kα辐射源（用于钛合金）的X350A X射线应力分析仪（iXpic，加拿大温哥华）在27 kV和8 mA的条件下测量了沿轧制方向（纵向）和横向的残余应力。采用经典的sin2φ方法在约15 μm的表面深度内检测宏观拉伸和压缩应力。在此过程中，记录了足够的2θ范围内的衍射角以确保高精度峰值分析。每个点的测量结果取三次平均值。X射线应力分析仪的仪器误差约为±15 MPa。因此，观察到的广泛残余应力分布并非由于测量限制，而是直接反映了冷轧过程中严重塑性变形引起的微观结构异质性（例如，相界处的局部应力集中、变形纹理和不同的晶粒取向）。由于硅钢板的厚度较薄，仅考虑了平面内的应力（纵向和横向），忽略了垂直于板面的应力。所研究的钛合金和硅钢的详细化学成分分别列在表1和表2中。比较两种钢材的应力分布有助于评估磁性质和晶体结构对机械性能的影响。

**2.2. 计算细节**
进行了第一性原理密度泛函理论（DFT）计算，以研究元素金属的基本机械和电子响应。虽然实验观察主要关注HCP和BCC合金系统，但计算框架扩展到了具有FCC结构的铝（Al）和铜（Cu），以及具有BCC结构的铁（Fe）和具有HCP结构的钛（Ti），从而涵盖了三种最基本的金属晶体家族。这些金属的理想晶体结构分别为：Al和Cu的FCC（空间群Fm3m）、Fe的BCC（空间群Im3m）和Ti的HCP（空间群P63/mmc），如图1所示。所有计算均使用CASTEP代码进行。初始几何优化通过总能量最小化来确定平衡基态结构。收敛标准设定为“精细”，并采用了全细胞松弛。对于变胞计算，使用了固定的基组质量。

随后施加了从-15 GPa到+15 GPa（递增2.5 GPa）的各向同性拉伸和压缩应力。虽然实际制造过程涉及各向异性剪切，但施加各向同性载荷可以为明确区分体积诱导的电子效应和各向异性几何伪影提供基本的理论基线。此外，尽管±15 GPa的应力范围远超过宏观工程屈服极限，但这一极端范围在物理上是必要的，用于探究完美无缺陷晶体的理想理论强度和极限动态稳定性。对于每个应力状态，进行了结构松弛以确定应力依赖的晶格常数。通过向优化后的结构施加小幅度应变来获得弹性常数（Cij）。然后使用标准的Voigt–Reuss–Hill近似方程评估宏观机械性能，包括杨氏模量和泊松比[23,24]。

此外，还计算了总态密度（DOS）、能带结构和声子色散谱。为了区分真实的物理能带结构变化和可能由应力引起的细胞变形的数值伪影，将费米能级严格作为所有电子结构计算的参考零点。交换-相关相互作用采用修正的Perdew–Burke–Ernzerhof（RPBE）泛函在广义梯度近似（GGA）中处理[25,26]。选择RPBE泛函是因为它能够减少标准PBE的过绑定倾向，为过渡金属提供可靠的晶格常数和体积模量。使用超软赝势[27]描述核-价电子相互作用，并采用了330 eV的平面波基组截止能量。

布里渊区采样采用11 × 11 × 11 Monkhorst–Pack k点网格[28]。这种相对密集的高对称性网格对于准确解析费米面的细节和确保DOS收敛至关重要。当总能量差低于1 × 10^-5 eV/原子、最大原子位移小于0.001 ?、残余应力低于0.05 GPa以及每个原子上的Hellmann–Feynman力低于0.03 eV/?时，认为收敛成立。结构优化是使用Broyden–Fletcher–Goldfarb–Shanno（BFGS）算法[29]进行的。所有模拟都考虑了自旋极化，以准确考虑磁性能对机械和电子特性的影响，尤其是在铁磁性的体心立方（BCC）Fe中。在所有材料中保持了相同的计算精度，以确保在施加应力下它们的机械、电子和振动特性的直接可比性。

3. 结果与讨论
3.1 残余应力分布
图2a和图2b分别展示了钛合金和硅钢沿轧制（纵向）和横向方向的残余应力结果。可以看出，钛合金样品的残余应力显著更高，应力分布也更广。这些发现突显了晶体结构和相组成对金属材料残余应力状态和机械行为的重要影响。它们还表明，在解释观察到的机械行为时，应考虑结构和磁性差异（顺磁性与铁磁性）。具体来说，观察到钛合金的最大残余应力明显高于硅钢，这提出了一个关于它们内在结构稳定性极限的基本假设。最终，这些宏观XRD测量的相关性在于为这项研究提供了实证基础：它物理上要求我们后续的原子尺度密度泛函理论（DFT）模拟来揭示为什么面心立方（HCP）结构能够承受更高的应力，而体心立方（BCC）结构可能更容易受到应力诱导的不稳定性。

3.2 第一性原理分析
图3描述了四种金属在拉伸和压缩应力下的晶格常数变化。在拉伸作用下，所有金属的晶格常数都呈现出大致线性的增加，膨胀的幅度顺序为：Al > Cu > Ti > Fe。Al对拉伸应力的敏感性最高，这表明FCC晶格中的Al–Al键更容易伸长，可能是由于Al的低内聚能和各向同性的键合网络。在压缩作用下，所有金属的晶格常数都减小，收缩的幅度顺序为Al > Ti > Cu > Fe。此外，Al的体积收缩最明显，表明其在压缩下的结构稳定性相对较低，这与它的低内在弹性模量一致。Ti的压缩响应比Cu更强，可能是由于HCP晶格固有的各向异性抗压能力。相比之下，Fe的晶格响应最小，反映了BCC结构的高晶格刚度和变形抗力。总体而言，不同的晶格对称性和键合配置导致了各种金属之间的应力响应差异。

图4显示了泊松比随施加应力的变化。在拉伸过程中，Cu和Ti的泊松比相对稳定，表明它们具有较好的横向变形能力。相比之下，Al在约5 GPa的拉伸应力后泊松比急剧增加，表明其横向弹性增强，这可能与FCC晶格的低剪切模量和高角度敏感性有关。Fe在约7.5 GPa时泊松比突然下降，可能是由于BCC晶格中的局部结构重排或电子轨道的突然重组。在压缩作用下，泊松比大多保持稳定，除了Cu在约10 GPa时的短暂跳跃，这可能是由于非线性晶格应变适应或电子密度突变所致。这些变化反映了横向弹性响应的结构依赖性，受晶格各向异性和滑移系特性的控制。

图5展示了在施加应力下杨氏模量的变化。在拉伸作用下，所有金属的模量值都呈下降趋势，表明键合减弱和晶格刚性降低。值得注意的是，Al的模量在约7.5 GPa后变为负值，表明机械不稳定性的开始。这表明FCC Al键已经超过了临界伸长阈值，可能导致内在键断裂和机械完整性的根本丧失。相比之下，Fe和其他金属在拉伸作用下的刚性降低更为渐进，表明它们具有更高的抗拉硬度。在压缩作用下，所有金属的杨氏模量通常增加，反映了由于间距减小而产生的更强的原子间相互作用。有趣的是，Cu在约10 GPa时显示出明显的峰值，随后趋于稳定。这可能表明了暂时的结构重排或中间相的形成，从而在达到新的平衡之前增强了刚性。其他金属的模量增加较为平滑。总体而言，拉伸倾向于软化金属系统，而压缩则导致刚性增加——这与应力诱导的键合重构一致。

图6和图7分别显示了四种金属在拉伸和压缩应力下的总态密度（DOS）曲线。在拉伸载荷下，所有四种金属的DOS峰都在费米能级附近增强且更加尖锐，表明局部电子态的增加和明显的轨道分裂。对投影态密度（PDOS）的分析表明，这些效应在过渡金属Ti和Fe中最为明显。这表明在拉伸作用下，高度局域化的d轨道能级发生了显著的重构，以及自旋极化的强烈调制。相比之下，Al和Cu中的s-p杂化态显示出更宽泛、对应力不那么敏感的响应。这些不同的轨道特定变化直接与观察到的结构软化相关。相反，在压缩应力下，DOS曲线趋于平坦，峰值强度减小，反映了电子态的增强非局域化和金属量的轻微减少。特别是Fe，在应力作用下其DOS的低能区域显示出强烈的不对称性，表明其局域化d态对外部应力具有高敏感性，可能涉及磁交换路径的位移和轨道杂化的变化。相比之下，Al在拉伸和压缩载荷下都保持了最稳定的DOS分布，与其s-p价电子主导和d态贡献最小相一致。

图8和图9分别展示了四种金属在拉伸和压缩应力下的能带结构演变。在拉伸作用下，导带向费米能级移动，靠近带边处出现重叠。这种带收敛增强了费米能级的态密度，提高了载流子迁移率和电导率。这种效应在Ti和Fe中最为明显，表明它们的轨道杂化和传输路径对应力具有强烈的敏感性。在压缩应力下，导带从费米能级处移开，伴随着带重叠区域的减小。这减少了可激发的热态数量，导致电子传输的抑制，略微降低了金属材料的特点。然而，在所有应力条件下，费米能级处的带重叠仍然持续存在，表明这些材料的内在金属特性在所施加的应力范围内得到了保持。

图10和图11展示了四种金属在拉伸和压缩应力下的声子色散曲线。这些谱提供了有关金属晶格在机械扰动下的动态稳定性的直接见解。Ti和Cu在整个应力范围内表现出相对稳定的声子频率，显示出坚固的动态稳定性和对外部机械波动的抵抗。然而，Al在10 GPa的拉伸应力下表现出负频率和声子模式的显著重构。特别是在高对称性边界点附近的横向声学（TA）分支中出现了明显的软化模式。这直接表明了动态不稳定性和应力诱导的结构崩溃的内在倾向。同样，Fe在高拉伸应力下也表现出明显的声子软化。尽管没有出现负频率，但在低至中等频率区域发生了显著变化，反映了由磁弹性耦合驱动的声子模式重新分布。虽然本研究中没有引入明确的点缺陷（例如空位或间隙），但Fe和Al中的异常DOS和能带结构特征表明它们对局部结构扰动具有较高的敏感性。应力诱导的电子结构重构可能降低了缺陷形成能量或扩散屏障，从而影响微观结构的稳定性和动力学行为。这些发现突出了应力调节缺陷动力学的潜力，这对于高性能合金的设计和长期可靠性至关重要。

4. 讨论
相应的应力-结构-性质关系在图12中以示意图的形式展示。实验和计算结果的结合揭示了具有不同晶体对称性的金属系统之间宏观应力与内在物理性质之间的一致性联系。XRD测量证实，钛合金的宏观残余应力更高，应力分布更广。这种宏观差异突出了由它们的主要晶体结构（钛合金的HCP和硅钢的BCC）和磁相组成所控制的不同机械行为。工程合金中的深刻宏观复杂性是我们进行基本第一性原理研究的主要动机，该研究旨在隔离和阐明基本晶体对称性的内在应力响应。我们明确承认，完美的、无缺陷的理论模型无法完全捕捉溶质原子效应（例如Al、V、Si）和支配实际多相合金的相界动力学。然而，这种故意的材料差异具有关键作用：通过解耦复杂的微观结构变量，纯净的单晶模型使我们能够清晰地分离出仅由晶格对称性驱动的纯几何应变效应。

图12示意性地说明了用于研究金属中应力诱导的结构和电子演化的综合实验-计算方法。该示意图突出了如何将XRD衍生的宏观残余应力数据与晶格、弹性、电子和振动性质的第一性原理模拟相结合，从而在施加应力下统一理解原子-电子耦合机制。因此，第一性原理模拟将这些宏观观察与原子尺度上的键合变化联系起来，作为一个通用的理论基础。FCC Al和Cu、BCC Fe以及HCP Ti在等效载荷下表现出明显不同的弹性响应。Al的高拉伸顺应性与其低内聚能和各向同性的键合网络相一致，而Fe较小的晶格敏感性反映了其本质上刚性的BCC框架和强定向键合。Ti的各向异性压缩响应与其实验残余应力分布很好地对应，突出了晶格对称性在调节机械稳定性中的作用。电子结构分析进一步阐明了这些机械行为的起源。重要的是，观察到的定性通用趋势——如拉伸软化、压缩硬化和电子带位移——是由内在轨道重叠控制的稳健物理现象，而不是所选RPBE泛函的假象。这些结果与金属弹性的通用理论基准一致。拉伸应力一致地使DOS峰在费米能级附近变尖锐，表明轨道局域化的增强和可能的电子拓扑转变，特别是在Fe和Ti中。这些变化可能会削弱金属的内聚性，降低弹性模量。相反，压缩倾向于使DOS特征变平，电子非局域化并增强键合。同时的能带结构变化——在拉伸作用下向费米能级移动，在压缩作用下远离费米能级——表明载流子浓度可调节，这是针对应变工程化功能金属的潜在兴趣点。

声子谱提供了对动态稳定性极限的洞察。虽然Cu和Ti在整个应力范围内保持稳定，但Al和Fe在高拉伸载荷下表现出声子软化，表明它们对结构不稳定性或缺陷激活的敏感性。这意味着，在涉及大拉伸应变的服务条件下，这些金属可能需要微观结构强化（例如晶界强化或沉淀硬化）来抑制不稳定性。重要的是，本研究指出了一个原子-电子耦合途径，通过该途径应力修改了原子间距离、改变了轨道重叠，并重新分配了电子态——这些机制反过来决定了观察到的弹性 and 振动响应。这种耦合在部分填充d带的系统中尤为明显，其中磁相互作用和轨道杂化强烈依赖于晶格间距。这些效应可以在高性能合金的设计中加以利用，通过有目的地引入特定的应力状态来定制机械硬度、电子传输或振动阻尼。此外，这些发现对合金设计和工业应用具有重要意义。对于应力与电子结构之间的耦合机制的更深入理解使工程师能够在加工过程中有目的地引入特定的形变状态，从而系统地调整金属性能。例如，应力引起的电导率或磁性能变化可以被用于柔性电子设备或应变传感器中，而控制结构部件中的残余应力可以提高抗疲劳性能。在工业应用中，这里确定的声子模式软化和发展带重构的理论阈值可以直接指导冷轧工艺参数的设定。通过利用这些内在的晶格不稳定性限制，工程师可以优化轧制工艺，使局部微应力保持在阈值以下，从而防止制造过程中的宏观故障。因此，这项研究为航空航天、汽车和能源行业高性能合金的设计以及应力工程化设备的发展提供了指导。

目前的工作结合了实验残余应力表征和第一性原理计算，其中宏观实验观测结果为基本的计算发现提供了现实的经济背景，极大地促进了计算结果的应用。未来的研究应该明确地将微观结构缺陷纳入模拟，并使用原位实验方法（例如基于同步辐射的应变依赖电子结构映射）来捕捉实时的应力-电子相互作用。这些扩展将进一步推进在复杂服役条件下的金属系统的预测性设计。

5. 结论

超越对个别金属的孤立案例研究，这项工作为基本晶体对称性的应力响应机制建立了一个通用的物理基准。通过将实验表征与第一性原理计算相结合，这项研究成功地将复杂的工程微观结构与内在的原子响应区分开来。宏观XRD分析揭示了钛合金和硅钢在残余应力分布方面的显著差异，强调了加工过程中异质工程合金的复杂机械行为。第一性原理结果表明，具有不同晶格对称性的金属（FCC、BCC和HCP）在外加载荷下的响应 fundamentally 不同：拉伸应力通常导致键长延长、模量降低和整体软化，而压缩则会增加晶格刚性。此外，金属的内在磁性能也会影响它们的应力响应。例如，铁这种铁磁材料的应力诱导的d轨道态密度变化表明，在外加载荷下其磁交换路径可能会发生潜在的改变，而在非磁性金属如铝和铜中没有这种现象。

电子结构分析显示，拉伸应力会增强费米能级附近的态密度（DOS），触发轨道重构，并降低导带，从而提高电导率。压缩则会导致DOS变平并减弱金属特性。声子光谱表明铜和钛具有动态稳定性，而铝和铁在高拉伸应力下表现出声子软化现象，这意味着它们的结构稳定性有限。这些结果表明，应力调制从根本上改变了键合特性、电子占据情况以及动态稳定性——为定制金属性能提供了坚实的理论指导。应力诱导的软化可能会在高应变条件下降低屈服强度和延展性，而应力诱导的强化可以提高硬度，这对于控制工程应用中的机械强度至关重要。此外，应力引起的能带结构和磁性的变化为传感器和智能材料的设计提供了机会。