基于分子动力学模拟的金属玻璃塑性变形原子尺度机制研究

金属玻璃作为一种通过快速冷凝制备而成的非晶态合金材料，因其内部原子排列呈现长程无序而短程有序的独特结构，展现出了远超传统晶态金属的优异性能。这种材料具备极高的强度、较大的弹性极限以及良好的耐磨耐腐蚀性，这些特性使其在航空航天、精密仪器及电子信息等前沿高端制造领域拥有广阔的应用前景。然而，金属玻璃在室温下往往缺乏宏观塑性变形能力，表现出明显的脆性断裂特征。一旦发生屈服，其变形通常会极快地局域化在极窄的剪切带内，导致材料在几乎没有预兆的情况下发生灾难性失效。这种室温脆性严重制约了金属玻璃作为结构材料的广泛应用，因此深入探究其塑性变形的内在微观机制，成为材料科学领域亟待解决的关键问题。

在实验研究手段方面，尽管现有的高分辨率电子显微镜等观测技术能够提供一定的微观结构信息，但由于金属玻璃塑性变形过程具有瞬时性和高度局域化的特点，直接捕捉原子尺度的动态演化轨迹依然存在极大的技术难度。分子动力学模拟技术作为一种基于经典牛顿力学的计算方法，能够从原子尺度上追踪粒子随时间的运动轨迹，从而为解析材料微观机制提供了强有力的工具。该技术通过建立原子间相互作用势函数，构建合理的物理模型，并对系统进行能量最小化与弛豫处理，进而模拟材料在特定载荷或温度条件下的动态响应过程。利用分子动力学模拟，研究人员能够精确地获得应力应变数据、原子位移矢量以及配位数变化等关键信息，进而从原子层面上揭示剪切带萌生、扩展及相互作用的详细物理图像。这种方法不仅克服了实验条件的限制，还能直观地展示变形过程中微观结构的演变规律。通过系统研究金属玻璃塑性变形的原子尺度机制，能够为设计高塑性的金属玻璃成分体系以及制定相应的加工工艺规范提供坚实的理论依据，这对推动新型高性能非晶合金材料的实际工程化应用具有重要的指导价值。

针对金属玻璃塑性变形的原子尺度模拟研究，构建科学合理的模型与精确设置计算参数是获得可靠物理结果的前提。模拟的核心任务在于从原子层面重现非晶态合金的结构特征与力学响应，这要求初始原子构型必须具备真实的非晶短程有序与长程无序特征。模型构建通常采用“熔融-急冷”法，即首先建立具有简单立方或面心立方结构的金属晶体模型，随后在远高于熔点的温度下进行充分弛豫以消除晶体记忆，最后通过极快的冷却速率将液态金属“冻结”至室温，从而形成金属玻璃结构。在确定模拟体系尺寸时，需兼顾计算效率与统计代表性，通常选取包含数万至数十万个原子的立方体盒子，确保模型尺寸远大于金属玻璃典型的剪切变形区特征尺度，以规避尺寸效应带来的误差。

原子间相互作用的精确描述依赖于势函数的选择，本研究选用嵌入原子势（EAM）或修正嵌入原子势（MEAM），该势函数能够准确描述金属键的多体相互作用，尤其适用于多元合金体系的能量与力计算。系统总能量计算公式为：

$$E_{total} = \frac{1}{2} \sum_{i} \sum_{j \neq i} \phi_{ij}(r_{ij}) + \sum_{i} F_i(\rho_i)$$

式中 $E_{total}$ 为体系总能量，$\phi_{ij}(r_{ij})$ 为原子 $i$ 与 $j$ 之间的对势能，$r_{ij}$ 为原子间距，$F_i$ 为嵌入能，$\rho_i$ 为其他原子在原子 $i$ 处产生的局域电子密度。模拟过程采用周期性边界条件以消除表面效应，模拟系综选择等温等压（NPT）系综以实现体积与压力的动态调节。温度与压力控制分别采用Nosé-Hoover热浴与Parrinello-Rahman压力浴算法。在力学加载阶段，通过改变模拟盒子尺寸施加恒定的单向应变或剪切应变，加载速率设定在 $10^8$ 至 $10^9 \text{s}^{-1}$ 量级，以平衡原子热运动与变形速率的关系。为确保结果可靠性，需对时间步长进行调试，通常取 $1$ 至 $2 \text{fs}$，并监控能量守恒与结构稳定性，直至模型达到平衡状态。

在金属玻璃塑性变形的研究中，追踪原子演化行为是揭示其微观变形机制的关键环节。本次研究依托分子动力学模拟方法，构建了能够精确反映原子运动的物理模型，通过对体系在不同应变阶段进行长时间序列的观测，实现了对原子尺度动态过程的捕捉与分析。这一过程的核心在于利用模拟输出的原子坐标与速度信息，结合可视化技术，将无序体系的结构演变转化为可观测的物理图像。

为了深入剖析原子层面的演化细节，本研究采取了多维度综合分析的策略。在原子位移层面，通过计算不同时刻原子相对于参考构型的位置矢量，能够直观地识别出体系中参与剧烈运动的活跃原子区域，这些区域通常与剪切带的萌生和扩展密切相关。针对原子键取向的分析，则侧重于考察邻近原子间排列顺序的改变，利用键对分析技术或键角分布函数，量化表征在塑性流动过程中原子成键模式的重构特征，从而判断短程有序结构在变形中的稳定性。

局部原子堆垛结构的变化是考察金属玻璃塑性起源的重要指标。研究通过引入 Voronoi 多面体指数等结构表征方法，统计分析了不同应变水平下局域原子团簇的类型与数量分布。这一分析能够有效揭示变形过程中诸如二十面体等高致密团簇的解体与重组规律，阐明自由体积是如何通过原子堆垛的调整而产生并聚集的。通过对不同加载条件下的模拟结果进行对比，能够发现原子演化行为既存在普遍共性，也表现出显著的差异性。这种共性通常体现在原子重排的局部化特征上，而差异则反映了应力状态或温度对原子迁移能力的不同影响。综上所述，对原子演化行为的系统追踪，为从原子尺度深入理解金属玻璃的塑性变形机制提供了坚实而直观的模拟观测基础。

在基于分子动力学模拟所获取的原子演化追踪数据基础上，对剪切转变区的识别是深入理解金属玻璃塑性变形机制的首要环节。本研究采用局部原子应变标准差作为关键判定参数，通过设定合理的阈值参数，精确提取出模拟体系中原子排列发生显著重排的集团区域，从而实现对剪切转变区的有效捕捉。通过对不同变形阶段剪切转变区的数量进行统计，结果显示在弹性阶段及塑性变形初期，其数量维持在较低水平且分布相对离散。随着外部加载应力的持续增加，体系内部能量逐渐积累，剪切转变区的数量呈现非线性增长态势，且其在空间分布上表现出从均匀弥散向局部聚集转变的特征。

针对剪切转变区的尺寸变化分析表明，单个剪切转变区在形核初期尺寸较小，主要局限于数个至数十个原子的范畴。在后续加载过程中，这些微小的激活区域通过吸纳周围原子的协同运动实现体积膨胀，呈现出明显的长大趋势。这一过程并非孤立进行，相邻的剪切转变区在空间扩展中逐渐相互接近并发生重叠与贯通，最终形成贯穿整个模拟模型的原子尺度剪切带。这一从独立形核到相互作用、再到连通成带的动态演化路径，完整再现了金属玻璃内部结构失稳与重排的物理图景。

剪切转变区的上述演化过程与宏观塑性变形的发生发展存在严格的对应关系。当剪切转变区处于弥散分布且相互独立时，材料宏观表现为均匀的弹性变形或微小的塑性屈服；而当大量剪切转变区在局部区域密集丛集并相互贯通时，体系内部产生不可逆的结构重排，宏观上则对应着显著的塑性流动及剪切带的形成。这充分证实了剪切转变区作为金属玻璃塑性变形基本载体的核心地位，其形核与长大的微观机制直接决定了材料宏观力学响应的强弱与变形模式。

在金属玻璃塑性变形的分子动力学模拟研究中，自由体积被视为描述原子尺度结构演化最关键的参量之一。自由体积是指原子堆垛过程中存在的、未被原子占据的微小间隙空间。通过对模拟数据的深入分析，能够清晰观测到在外加应力作用下，自由体积并非均匀分布，而是呈现出明显的空间聚集特征。随着剪切带的形成与扩展，高自由体积区域逐渐在剪切带内部连通，导致局部原子体积发生显著膨胀。这种空间分布的变化直接反映了材料内部结构的软化过程，是塑性变形在微观层面的直接体现。

在塑性变形过程中，自由体积的迁移规律表现出高度的动态性。原子通过位置交换，使得自由体积在三维空间内发生扩散与流动。模拟数据显示，自由体积的迁移速率与局部原子重排的强度呈现正相关关系。当自由体积浓度超过某一临界值时，原子间的约束力减弱，原子团簇更易于发生协同重排。这种重排并非单个原子的独立跳跃，而是涉及数十个原子的集体运动，从而构成了剪切转变区的核心结构单元。

自由体积的迁移通过促进局部原子重排，直接调控着剪切转变区的形核与扩展机制。在形核阶段，自由体积的富集降低了局部剪切模量，为剪切转变区的启动提供了结构基础。而在扩展阶段，自由体积的持续流动使得剪切转变区能够克服周围基体的弹性阻力，诱发邻近区域的原子重排，进而导致剪切带的快速贯通。这一过程表明，自由体积不仅是塑性变形的产物，更是驱动变形持续发展的动力源。

综上所述，自由体积与原子重排的协同作用深刻影响着金属玻璃的宏观塑性变形行为。高密度的自由体积迁移与剧烈的原子重排对应着材料良好的塑性，而自由体积的匮乏则会导致脆性断裂。明确这一调控机制，对于通过成分设计改善金属玻璃的塑性具有重要的指导意义。

本文通过对金属玻璃塑性变形过程的分子动力学模拟，深入揭示了其在原子尺度下的变形机制与微观结构演化规律。研究表明，金属玻璃作为一种长程无序而短程有序的非晶态材料，其塑性变形行为与内部原子团的拓扑结构及能量状态密切相关。在模拟过程中，通过构建合理的非晶模型并施加剪切载荷，能够清晰观察到应力诱发下的自由体积产生与湮灭过程，这直接对应着材料宏观上的屈服与硬化现象。核心原理在于，金属玻璃不具备传统晶体材料中的位错滑移机制，其塑性流动主要依赖于局部原子团簇的剪切转变。这些微小的剪切转变事件在应力作用下随机形核并逐渐扩展，形成贯穿样品的剪切带，进而导致材料发生宏观失效。

操作层面上，分子动力学模拟通过求解牛顿运动方程，精确追踪了成千上万个原子的实时位置与速度变化，从而计算出体系的应力-应变响应及原子位移矢量。这种从原子层面出发的研究方法，克服了实验手段在观测瞬时微观结构变化方面的局限性，验证了剪切带内部原子具有显著的高流动性特征，而带外基体则保持相对刚性的弹性变形。这一发现对于理解金属玻璃的室温脆性断裂问题具有重要的理论指导意义。在实际应用中，明确塑性变形的原子尺度机制，有助于材料科研人员通过成分设计或工艺调控来优化材料的自由体积分布，进而改善金属玻璃的塑性加工能力与抗断裂性能。研究结果不仅丰富了非晶态物理的基础理论，也为开发高性能、高韧性的金属玻璃结构材料提供了科学依据，推动了该先进材料在精密器件及航空航天等关键领域的工程化应用进程。