透闪石中的缺陷物理与铁替代：复杂硅酸盐的第一性原理研究

《Acta Materialia》：Defect physics and iron substitution in tremolite: A first-principles study of a complex silicate

　　

复杂硅酸盐矿物具有化学成分多样且结构框架复杂的特点，然而其缺陷形成及缺陷介导性能的控制机制仍知之甚少。研究人员利用第一性原理计算，系统研究了代表性钙镁闪石——透闪石中的本征点缺陷形成与Fe替代。研究发现缺陷形成存在显著的亚晶格依赖性：骨架阳离子空位形成能高且位点依赖性弱，而氧空位对位点高度敏感并在p型条件下具有热力学可及性。核心发现是电荷补偿普遍主导缺陷稳定性，决定了本征缺陷与外征缺陷截然不同的静电通路。本征反位缺陷通过与本征氧空位形成施主-受主配对而获得稳定化，而外征Fe掺杂则触发相反的补偿逻辑，优先与本征受主配对形成高度稳定的中性缺陷团簇。至关重要的是，这些外征补偿通路的激活受局域水合与氧化还原条件控制。这种环境依赖的通路选择调控着Fe2?/Fe3?比例，为天然闪石中观察到的价间电荷转移和系统颜色变化（从浅绿色到深棕色）提供了原子尺度起源。更广泛地，这些结果确立了电荷补偿作为复杂硅酸盐中缺陷形成的统一原理，并证明通过本征补偿网络进行过渡金属掺杂是比调控孤立本征缺陷更有效的性能调控途径。

透闪石中的缺陷物理与铁替代：复杂硅酸盐的第一性原理研究——论文解读

**研究背景与意义**

硅酸盐矿物是地球上最丰富的无机固体类别，构成了广泛地质、环境及先进工业系统的结构基础。传统上被视为化学惰性和电子惰性的硅酸盐，已被大量实验和理论研究表明其能够承载点缺陷、杂质掺入及辐照诱导中心，这些缺陷会改变局域电子结构并产生光学活性、电荷俘获及缺陷介导的化学活性。从材料化学视角看，硅酸盐在宽带隙氧化物中占据独特地位：与简单二元氧化物不同，其延展的Si–O框架施加严格的拓扑约束，同时容纳相当大的局部结构灵活性。闪石类矿物是其中特别复杂且探索不足的一类，其晶体结构由角共享SiO?四面体构成的一维双链组成，通过二价阳离子（如Mg2?和Ca2?）占据的边共享八面体带连接，常伴有结构羟基。透闪石（Ca?Mg?Si?O??(OH)?）是典型的闪石矿物，也是天然软玉（jade）的主要成分。尽管已有广泛的矿物学和光谱学表征，但透闪石中本征点缺陷及缺陷-掺杂剂相互作用的系统性原子尺度理解仍然缺乏。过渡金属掺杂进一步丰富了硅酸盐的缺陷景观。在天然透闪石-阳起石固溶体系列中，观察到从浅绿色到深棕黑色的系统颜色演变，由Fe2?和Fe3?物种的不同比例驱动，调控着价间电荷转移光学跃迁的强度。地质研究表明，体相Fe2?/Fe3?比例对水合和氧化还原环境等外部条件敏感响应，但这种依赖性的原子尺度热力学起源仍不清楚。本研究旨在建立连接复杂硅酸盐宿主中缺陷形成能、电荷补偿和热力学转变水平的统一微观框架。

**研究开展与主要结论**

研究人员对透闪石中所有晶体学不等价位点的空位、反位和间隙缺陷进行了系统考察，评估了其形成能、电荷态稳定性和转变水平。研究发现缺陷形成存在明显的亚晶格依赖性：骨架阳离子空位形成能高且弱位点依赖，而氧空位强烈位点敏感且在p型条件下热力学可及。核心发现是电荷补偿普遍主导缺陷稳定化：本征反位缺陷通过与本征氧空位形成施主-受主配对获得稳定，而外征Fe掺杂触发相反的补偿逻辑，与本征受主配对形成高度稳定的中性缺陷团簇。这种环境依赖的通路选择控制着Fe2?/Fe3?比例，为天然闪石中观察到的价间电荷转移和颜色变化提供了原子尺度起源。论文发表在《Acta Materialia》。

**关键技术方法**

研究人员采用基于密度泛函理论（DFT）的第一性原理计算方法，使用维也纳第一性原理模拟软件包（VASP）实现。交换关联作用采用Perdew-Burke-Ernzerhof（PBE）广义梯度近似，电子-离子相互作用采用投影缀加平面波（PAW）方法描述。计算采用550 eV平面波动能截断，并考虑DFT-D3色散校正。缺陷计算使用含82个原子的透闪石常规晶胞，采用3×1×5的Γ中心k点网格采样。缺陷形成能计算采用标准超胞形式，其中有限尺寸静电校正采用Kumagai和Oba的各向异性扩展Freysoldt（eFNV）方案。元素化学势由多维热力学相空间约束确定，选取透闪石稳定性窗口内的代表性顶点条件（A点）作为基准。

**研究结果**

**3.1 晶体结构与电子性质**：透闪石结晶于单斜C2/m空间群，优化晶格常数与实验晶体学数据吻合良好（偏差约0.5%以内）。计算的电子带隙为4.96 eV，呈现绝缘体特性。价带顶（VBM）由O 2p态主导，导带底（CBM）主要源于Si–O反键杂化态。由于带边由硅酸盐框架而非阳离子衍生态控制，氧相关缺陷预计引入深转变水平，而大多数阳离子中心缺陷应保持电子良性。

**3.2 稳定化学势区域**：构建了Ca–Mg–Si–O–H五元体系的多维热力学相图，考虑的竞争相包括二元氧化物（MgO、CaO）、氢氧化物（Mg(OH)?、Ca(OH)?）、无水硅酸盐（SiO?、镁橄榄石Mg?SiO?、斜硅钙石Ca?SiO?、透辉石CaMg(SiO?)?）以及含水镁硅酸盐（滑石Mg?Si?O??(OH)?和蛇纹石Mg?Si?O?(OH)?）。在Δμ_O=?1.50 eV条件下，选取Mg(OH)?和Ca(OH)?同时饱和析出的顶点作为缺陷计算基准。

**3.3 缺陷形成与自补偿**：碱土阳离子空位（V_Ca和V_Mg）形成能为3.80–4.58 eV（VBM处），在n型条件下降至约?5.0 eV，作为强电子陷阱。硅空位（V_Si）形成能高达7.94–8.48 eV，反映切断四个Si–O共价键和破坏刚性双链骨架的高昂代价。氧空位（V_O）在p型条件下形成能可降至?0.70 eV（+2电荷态），呈现强位点依赖性；O-1至O-4位点暴露于结构空腔，具有较大弛豫自由度，而O-5和O-6受刚性拓扑约束。氢空位（V_H）作为深受主，形成能2.44 eV；羟基空位（V_OH）偏好中性态，形成能恒为4.36 eV。同价反位缺陷（Ca_Mg、Mg_Ca）形成能低（0.73–0.89 eV），而异价反位（Ca_Si、Mg_Si）形成能高（5.05–5.92 eV），但在n型条件下被强烈稳定化。间隙缺陷中，碱土间隙（Ca_i、Mg_i）为强施主，形成能分别为?2.28和?0.63 eV。

**3.4 施主-受主补偿**：氧空位（V_O••）作为双施主，可与受主型反位缺陷Ca_Si″和Mg_Si″通过库仑吸引形成电中性缺陷团簇(Ca_Si-V_O)×和(Mg_Si-V_O)×，形成能分别为3.49 eV和3.40 eV，相比孤立反位缺陷降低超过2 eV。这种稳定化不仅源于静电补偿，还涉及局域电子转移和晶格畸变释放。

**3.5 缺陷能量学的稳健性与热力学敏感性**：在化学势窗口的三个极端点（A、B、C）评估发现，V_Si形成能始终超过7.5 eV，确立骨架破坏在所有平衡条件下热力学不可及。V_O对化学势不敏感但强烈依赖费米能级：VBM处为放热双施主，CBM处升至禁止值。受主型缺陷（V_Ca、V_Mg、Ca_Si、Mg_Si）呈现相反行为。异价替代显示化学势敏感性差异：Ca相关缺陷的能量位移大于Mg相关缺陷，反映了化学势窗口宽度的差异。

**3.6 缺陷转变水平**：V_Ca引入VBM上方0.32和0.95 eV两个受主水平。V_Mg根据位点不同引入浅受主（0.04–0.33 eV）和深受主（0.87–1.00 eV）水平。V_Si表现出多电子俘获特征。V_O位点1–4表现为深双施主，具有非相邻(+2/0)转变（2.47–2.82 eV）；位点5–6具有稳定的+1中间态。V_H表现出两性行为，施主水平0.11 eV，受主水平1.32 eV。由于V_O和V_H的转变水平深在带隙中，室温下不能被热电离提供自由载流子，主要作为深载流子俘获中心而非有效掺杂剂。

**3.7 Fe替代与氧化还原驱动的电荷补偿**：Fe优先占据Mg位点（Fe_Mg）。在p型条件下（氧化环境），+1有效电荷态（Fe_Mg•，Fe3?）被优先稳定；在n型条件下（还原环境），中性态（Fe_Mg×，Fe2?）占优。中性Fe_Mg×缺陷形成能低至0.72–0.78 eV，表明Fe易于替代进入Mg亚晶格。孤立Fe_Mg引入深转变水平；Fe_Mg×（Fe2?态）引入4.00 μ_B磁矩，主要由Fe 3d轨道贡献。由于Fe3?施主与V_O••施主之间存在施主-施主排斥，氧化态Fe3?掺杂剂需要本征受主（V_H′、V_Ca″、V_Mg″）实现电荷补偿。(Fe_Mg-V_H)×团簇形成能与孤立替代相当（0.76 eV），而(2Fe_Mg-V_Ca)×和(2Fe_Mg-V_Mg)×团簇具有负形成能（?0.18和?0.19 eV），远低于孤立空位成本。磁矩分析确认团簇中每个Fe原子为高自旋3d?构型（Fe3?态），每个Fe捐赠一个价电子补偿邻近受主空位。

**3.8 Hubbard U效应**：GGA+U计算（U_eff=1–5 eV）显示，孤立Fe_Mg×（Fe2?）形成能对U弱依赖（0.18–0.31 eV）。两个三聚体Fe3?团簇随U增大而进一步稳定，U=5 eV时(2Fe_Mg-V_Mg)×达?1.10 eV，(2Fe_Mg-V_Ca)×达?0.87 eV。(Fe_Mg-V_H)×对U几乎不敏感（约+0.5 eV），归因于半满壳层能量对在位库仑修正的不变性。化学势分析表明，随着脱水进行(Fe_Mg-V_H)×逐渐稳定化，最终形成能变为负值，激活脱氢-氧化路径。

**讨论总结**

本研究表明透闪石的缺陷化学兼具稳健与脆弱特征：阳离子空位形成能位点不敏感，源于修饰阳离子的离子键合特性和硅酸盐框架的延展连通性；而阴离子相关缺陷揭示位点敏感性，氧空位在p型条件下因完全电离为+2态且静电未屏蔽，对局部几何灵活性高度敏感。电荷补偿被确立为复杂硅酸盐中缺陷稳定化的普遍驱动力：本征受主型反位与氧空位施主配对，而外征过渡金属Fe3?施主与本征受主配对，遵循严格静电选择规则。环境依赖的缺陷团簇化机制为天然闪石的光学性质和颜色变化提供了微观起源：富水环境对应非补偿Fe2?替代（浅绿色），部分脱水激活氢空位配对，无水条件促进骨架空位团簇化（高Fe3?浓度，深棕黑色）。这些发现为复杂硅酸盐和先进陶瓷材料的缺陷工程提供了实用指导：通过引入多价过渡金属激活本征补偿网络，是比操纵孤立本征缺陷更有效的性能调控策略。

**研究结论翻译**

研究人员对透闪石中的本征点缺陷和Fe替代进行了第一性原理研究。基于各亚晶格的结构角色和相关的电荷不平衡，缺陷形成呈现出清晰的层级结构。阳离子空位在p型条件下具有高且弱位点依赖的形成能（图3a），但在n型条件下作为强稳定化的本征电子陷阱。氧空位能量上更可及但表现出巨大的位点敏感性，尤其在p型条件下，空腔暴露位点与结构约束位点之间的结构弛豫差异显著（图3b）。Mg位点上的Fe替代在整个费米能级范围内能量上高度有利（图7a）。本征和外征缺陷主要在宽带隙内引入深转变水平（图6），充当深俘获中心。本工作的核心成果是确定电荷补偿为复杂硅酸盐中缺陷稳定化的主导机制。施主-受主补偿显著降低缺陷形成能，表明电荷补偿缺陷团簇而非孤立带电点缺陷支配缺陷种群。优选的团簇化构型由严格静电选择规则动态决定：本征受主型反位仅与本征氧空位施主结合，而外征多价Fe施主触发反转的补偿逻辑，仅与本征受主（如质子或骨架阳离子空位）配对，形成异常稳定的中性缺陷团簇。最终，这种环境驱动的缺陷-掺杂剂耦合控制了材料的宏观价态比和光学性质。更广泛地说，有意识地利用这些本征补偿网络的受控过渡金属替代，为调控复杂硅酸盐和先进陶瓷材料的功能性质提供了现实可行的能量学途径。