基于第一性原理计算的TiAl合金高温氧化机理的理论研究

TiAl合金是轻质高温结构材料，因高温力学性能出色且密度低，在航空航天、汽车制造等领域有很大应用潜力。不过在高温环境中，氧化现象会使材料性能下降，这大大限制了其在工程中的应用。所以要深入探究TiAl合金的高温氧化机理，这既具备理论价值又具有实际意义。

第一性原理计算是基于量子力学的理论工具，它可以从原子层面揭示材料氧化过程的微观机制，能够为优化合金成分、提高抗氧化性能提供理论方面的支持。该计算方法的核心是求解薛定谔方程，通过获取材料体系的电子结构和能量信息来预测其物理化学性质。

在研究TiAl合金高温氧化机理时，第一性原理计算方法通常会涉及几个关键环节，包括构建表面吸附模型、计算氧原子的吸附能和扩散势垒、分析氧化层形成时的电荷转移情况。通过开展这些计算工作，能够明确氧原子在TiAl合金表面的初始吸附行为、氧化层的生长动态，以及合金元素对氧化过程产生影响的具体机制。

在实际应用当中，第一性原理计算这种研究方法不仅能够帮助理解TiAl合金的高温失效机制，还可以为开发新型抗氧化涂层、设计合金成分提供科学方面的依据，从而延长材料在极端环境中的使用时间。当第一性原理计算和实验表征结合在一起之后，能够形成一个理论预测与实验验证的闭环研究体系，这可以显著提高材料研发的效率和准确性。

基于第一性原理计算研究TiAl合金的高温氧化机理，对于推动TiAl合金在高温工程领域的应用有着重要的实践意义。

TiAl合金在高温下氧化，会经历从原子吸附到氧化层逐步生长的多尺度微观变化过程。氧化开始时氧原子在γ - TiAl和α₂ - Ti₃Al两种相表面的吸附情况直接影响后续氧化发展方向。用密度泛函理论（DFT）计算可知，氧原子更倾向附着在配位数较高的Ti位置，这种吸附能力能用公式$E${ads} = E{O/slab} - (E{slab} + \frac{1}{2}E{O_2})来进行定量表示。研究数据表明，γ - TiAl(111)晶面桥位吸附能约为 - 5.2 eV，明显比α₂ - Ti₃Al(0001)晶面的 - 4.6 eV低，这解释了γ相为何更容易生成富含Ti的氧化层。

随着氧含量逐渐升高，最初生成的AlO和TiO两种氧化产物成核所需能量壁垒有明显差异。计算结果显示，AlO成核能垒约为1.8 eV，而TiO成核能垒仅1.2 eV，所以TiO会先在材料表面形成核结构。氧化层逐步生长时，O²⁻、Al³⁺和Ti⁴⁺等离子会发生扩散移动，它们的扩散速度由阿伦尼乌斯方程$D = D_0 \exp\left(-\frac{Q}{RT}\right)$来决定。具体来说，Al³⁺在Al₂O₃晶界处扩散的活化能大约是3.5 eV，而O²⁻在TiO₂内部扩散的活化能高达4.2 eV，这种扩散能力的不同最终使得氧化层呈现出分层结构。

目前相关研究的争议主要集中在合金与氧化层界面的反应过程方面。有学者觉得Ti原子会先穿过界面生成TiO₂过渡层，也有看法认为Al元素的聚集会促使连续Al₂O₃保护层形成。这些微观机制尚不明确的状况，正好为第一性原理计算的进一步研究提供了方向。

第一性原理计算方法基于量子力学基本原理，通过求解体系电子的薛定谔方程来得到材料的微观结构和性能信息。本研究选用密度泛函理论（DFT）处理多电子体系，其基本方程可表示为：

这里的有效势\( V_{\text{eff}}(\mathbf{r}) \)包含电子与原子核的相互作用以及电子间的交换关联效应。计算过程用VASP软件包，交换关联能的描述选用广义梯度近似（GGA）下的PBE泛函，并且引入自旋极化效应，这么做是为了更准确地处理Ti的3d电子。平面波基组的截断能设定为500 eV，布里渊区积分采用Monkhorst - Pack方法，k点网格密度保持在0.04 Å\(^{-1}\)，以此确保能量收敛精度能够达到10\(^{-5}\) eV。
在模型构建方面，搭建γ - TiAl(111)和(001)表面模型，采用周期性边界条件并且添加15 Å的真空层，添加真空层的目的是消除层间的相互作用。O原子吸附模型的覆盖度设置在0.25 - 1.0 ML之间，这样做是为了考察不同氧化程度带来的影响。同时搭建Al\(_2\)O\(_3\)/TiAl界面模型，通过分析界面结合能的方式来探究氧化机制。为了验证模型的可靠性，先计算纯TiAl的晶格常数，得到a = 3.97 Å、c = 4.04 Å，将其与实验值（a = 3.98 Å、c = 4.05 Å）相比较，误差小于0.3%，这表明所建模型具有较高的准确性，计算参数也是比较合理的。这一系列的参数设置工作以及模型构建工作，能够为后续分析氧化机理提供可靠的理论支撑，在后续对于氧化机理的分析中，这些参数设置和模型构建所提供的理论支撑会起到非常重要的作用，使得后续的分析能够基于一个准确且合理的基础之上进行。

### 2.3TiAl合金氧化初期过程的模拟与验证

弄清楚TiAl合金的高温氧化机理，对其氧化初期过程进行模拟和验证是非常重要的一个步骤。这一过程主要依靠第一性原理计算方法，从原子层面去研究氧原子在合金表面的吸附情况、扩散情况以及氧化物成核的具体情况，并且要结合实验数据来验证模拟结果是不是可靠。具体做法是先计算氧原子在TiAl不同晶面上的吸附能，通过计算吸附能来确定最稳定的吸附构型。

吸附能按照下面的公式进行计算：

这里面$E_{\text{slab+O}}$指的是吸附氧原子后体系的总能量，$E_{\text{slab}}$指的是干净表面的能量，$E_{\text{O}_2}$指的是氧气分子的能量。通过比较不同吸附位点的吸附能具体数值，就能够知道氧原子更倾向于占据哪些位置。
之后用NEB（弹性带）方法来计算氧原子从表面向次表面扩散的能垒，计算这个能垒有助于明确氧化过程当中的速率控制步骤。计算扩散能垒的时候需要确定初始态、过渡态和终态的能量差，通常情况下，能垒要是越高，氧原子扩散就越困难。与此同时还要模拟AlO和TiO的成核过程，通过计算成核能垒来仔细分析这两种氧化物形成的难易程度。

为了验证模拟结果是不是准确，需要把计算得到的吸附能、扩散能垒和XPS（X射线光电子能谱）、STM（扫描隧道显微镜）的实验数据，或者是其他的理论结果进行对比。在分析结果一致性和存在差异的原因时，需要考虑泛函选择、温度效应等因素所带来的影响。举例来说，GGA泛函有可能会低估能垒的具体数值，实验条件当中存在的杂质或者缺陷也有可能让理论值和实测值出现偏差。采用理论和实验相结合的方式，能够有效地验证模拟结果是不是可靠，还能够为TiAl合金的高温氧化防护提供相应的理论支持。

深入理解TiAl合金在高温下的氧化表现，关键是掌握氧化层生长动力学以及界面反应的作用机制。氧化层生长的快慢主要由氧离子（O²⁻）和金属离子（Al³⁺、Ti⁴⁺）在氧化层内部的扩散过程决定。按照菲克定律，扩散通量$J$可以用公式$J = -D \nabla C$来描述，其中$D$是扩散系数，$C$为浓度。借助第一性原理计算，把迁移能垒$\Delta E_m$和振动频率$\nu$结合起来，再通过过渡态理论就能够确定扩散系数，具体的公式为：

在这个公式里，\(a\)指的是原子跃迁的距离，\(k_B\)是玻尔兹曼常数，\(T\)代表的是温度。对Al₂O₃和TiO₂中体扩散与晶界扩散的能垒进行计算后可以发现，晶界作为快速扩散通道，会使得氧化层的生长速度加快。
在研究界面反应机理的时候，需要构建合金与氧化层之间的界面模型，并且计算界面处Al原子扩散和氧化反应的能垒。界面反应能垒\(\Delta G_r\)会对氧化层与合金界面的反应快慢产生影响，其表达式为：

对这些能垒进行分析，就能够找出Al发生选择性氧化的关键步骤。之后结合动力学蒙特卡罗（KMC）方法，就能够模拟在1000 K时氧化层厚度$d(t)$随着时间$t$的变化规律，这种变化一般是符合抛物线关系的，其公式为：

这里的$k_p$是抛物线速率常数，它的大小是由氧化层的相组成、晶界缺陷密度以及界面结合强度共同来决定的。晶界缺陷会让扩散能垒降低，而界面结合强度则直接和氧化层是否容易剥落有关系，这两方面的因素会一起对TiAl合金的高温抗氧化能力产生影响。也就是说，晶界缺陷的存在使得扩散能垒降低，这会在一定程度上促进氧化层的生长，而界面结合强度若不足，氧化层就容易剥落，这两者相互作用，综合起来影响着TiAl合金在高温环境下抵抗氧化的能力。

研究合金元素对TiAl合金高温氧化行为的影响，第一性原理计算是有效方法，可开展全面分析。第一步要搭建包含Cr、Nb、Si等合金元素的TiAl晶体模型。掺杂方式有取代型和间隙型两种，取代型如Cr原子取代Ti或Al的晶格位点，间隙型如Si原子占据八面体间隙位置。通过计算纯TiAl和掺杂体系表面O原子的吸附能来评估合金元素对氧化初期过程的影响。

吸附能的计算公式是：

这里面，\(E_{\text{TiAl+O}}\)指的是吸附O原子后体系的总能量，\(E_{\text{TiAl}}\)是纯净表面的能量，\(E_{\text{O}_2}\)是氧分子的能量。从计算结果能够知道，Si掺杂可以明显地降低O原子的吸附能，进而起到抑制氧化初期反应过程的作用。
计算扩散能垒采用爬坡弹性带方法（CI - NEB），通过对比O原子在纯TiAl和掺杂体系中的迁移路径进行分析。不同体系的能垒存在差异，这种差异可以反映出合金元素对氧原子扩散的阻碍效果。就像Nb掺杂之后，O原子的扩散能垒大约提高了0.3 eV，而这会让氧化反应的速率变慢。界面偏聚能的计算公式是：

要是偏聚能为负值的情况出现，那就说明合金元素更愿意在氧化层与基体的界面处富集。Cr和Si呈现出较强的负偏聚能，这就能够解释它们在界面处出现偏聚的现象。

在计算氧化层形成能的时候发现，SiO₂的形成能是 - 9.1 eV，明显低于TiO₂的 - 4.8 eV，这意味着Si更容易优先发生氧化，从而形成致密的保护膜。界面结合能的计算结果还显示出，SiO₂与TiAl的界面结合强度比TiO₂与TiAl的界面结合强度更高，这和在实验当中观察到的Si能降低TiAl氧化速率的现象是相互符合的。把计算得到的数据和实验得到的数据结合起来看，合金元素主要是通过对吸附能、扩散能垒以及氧化层的稳定性进行改变，从而共同对氧化行为产生影响，在这其中Si的“选择性氧化”作用表现得最为突出明显。

这项研究运用第一性原理计算方法来系统地对TiAl合金在高温环境当中的氧化机理展开探究。研究着重对氧化过程里原子扩散、界面反应以及氧化层形成这些微观机制进行剖析。研究人员搭建了TiAl合金表面模型，以此模拟氧原子在合金表面的吸附和扩散情况，得到的结果是发现氧原子存在更倾向于占据Ti原子间隙位置这样的动力学特性。

从计算结果能够知道，TiAl合金初始氧化阶段主要是由Ti元素优先氧化起主导作用，这种优先氧化会生成致密的TiO₂保护层。与此同时Al元素通过扩散作用在氧化层和基体的界面处聚集起来，随着时间的推移逐渐形成连续的Al₂O₃阻挡层。这种双重氧化层结构可以有效地抑制氧原子进一步渗透到合金内部，从而让合金的高温抗氧化能力得到提高。

在研究过程中还发现，温度和氧分压会对氧化层生长速率以及结构稳定性产生显著的影响。在高温环境之中，TiO₂层的晶格缺陷数量会增加，而晶格缺陷数量的增加会加快Al原子的扩散速度，进而推动Al₂O₃层的形成。另外研究通过分析能带结构和态密度，把TiAl合金氧化过程中的电子转移机制明确了出来，这为优化合金成分设计提供了理论方面的支撑。

这项研究具有重要意义，它不仅让人们对TiAl合金高温氧化行为的认识得到了加深，而且为开发新型耐高温合金材料奠定了坚实的基础，所以具备重要的学术价值以及广阔的工程应用前景。