重新审视晶体和非晶TiO2中的氮掺杂

《physica status solidi (b)–– basic solid state physics》：Revisiting the Nitrogen Doping in Crystalline and Amorphous TiO2

　　

在这项工作中，研究人员提出了一项关于二氧化钛（TiO₂）中一个众所周知缺陷——即替代氧位的氮（N_O）——的综合方法论研究，主要关注点在于对晶体和非晶相中三种不同电子结构方法的比较评估。尽管文献中已广泛研究了氮掺杂TiO₂，但缺乏同时针对主要晶体多晶型（锐钛矿和金红石）以及非晶二氧化钛（titania）的理论方法系统比较，尤其是对于后者，定量数据仍然稀少。从标准密度泛函理论（DFT）出发，研究人员将分析扩展到已充分验证的Lichtsteiner的DFT+U方法，最终引入Teles等人最近提出的DFT-1/2方法，该方法具有完美再现TiO₂及许多其他半导体实验带隙的优点。由于缺陷态的描述强烈依赖于准确的电子结构，这相对于传统的DFT（作为一种基态理论，系统性地低估带隙）是一个实质性的改进。详细讨论了各方法之间的相似性和差异。

**论文解读：重新审视晶体与非晶TiO₂中的氮掺杂**

**研究背景与问题**

二氧化钛（TiO₂）作为一种重要的光催化材料，在分解水制氢、污染物降解等领域具有广泛应用。然而，其宽带隙（约3.0–3.2 eV）限制了可见光吸收，仅能利用约4%的太阳光谱（紫外区）。通过掺杂非金属元素（如氮）引入缺陷态，可在带隙中创建局域能级，从而增强可见光响应。尽管已有大量理论和实验研究，但针对替代氧位的氮（N_O）缺陷，系统比较不同电子结构方法（特别是同时涵盖晶体多晶型——锐钛矿和金红石——以及非晶相）的研究仍然缺乏。标准密度泛函理论（DFT）因自相互作用误差而严重低估带隙，而更精确的杂化泛函或格林函数方法计算成本高昂。因此，需要一种低成本且准确的从头算方法，能够处理大超胞缺陷计算。本研究旨在回答以下方法论问题：电子结构方法的选择如何影响缺陷态的结构和电子性质预测？能否为有序和无序TiO₂系统建立统一可靠的计算协议？

**研究内容与意义**

研究人员以固定浓度（~0.78%）的替代氮（N_O）缺陷为研究对象，系统比较了标准DFT、DFT+U（采用Lichtsteiner线性响应方法确定U值）以及近期发展的DFT-1/2方法（Teles等人提出，可完美再现实验带隙）在锐钛矿、金红石和非晶TiO₂中的表现。研究发现，DFT-1/2方法能准确描述带隙和缺陷态位置，与高成本准粒子方法结果一致，且计算效率高。论文发表在《physica status solidi (b)–– basic solid state physics》。该工作填补了方法学对比的空白，尤其为非晶TiO₂提供了定量数据，为后续缺陷工程和光催化性能优化奠定了可靠的理论基础。

**主要关键技术方法**

1. **第一性原理计算框架**：基于自旋极化DFT，采用VASP代码、PBE交换关联泛函和PAW赝势，平面波截断能600 eV。晶体相（锐钛矿、金红石）和非晶相（192原子超胞，熔融淬火法生成，来自Prasai等人工作）均在相同条件下处理。
2. **DFT+U校正**：采用Cococcioni线性响应方法，从第一性原理确定Ti 3d电子的Hubbard U参数（锐钛矿5.67 eV，金红石5.88 eV，非晶相按配位环境加权平均为5.55 eV）。
3. **DFT-1/2校正**：基于Slater半占据技术，通过优化氧、钛和氮的PAW/PBE赝势截断半径，最大化带隙。对于氮杂质，开发了双参数优化协议，同时最小化基本带隙并最大化导带底（CBM）与氮深能级之间的能量差，确定最优截断半径（金红石3.2 Bohr，锐钛矿3.1 Bohr，非晶相3.2 Bohr）。

**研究结果**

**3.1 DFT与DFT+U分析晶体TiO₂特征**

通过DFT+U优化，锐钛矿和金红石的晶格参数与实验吻合良好，带隙分别提升至2.767 eV（间接）和2.456 eV（间接），显著缩小了与实验值的差距。对于氮掺杂锐钛矿，DFT和DFT+U均显示两个N诱导态：一个深中隙态（约1.3 eV高于价带顶VBM）和一个浅态紧邻VBM上沿。DFT+U增强了该深态的局域化（N 2p成分从42%增至50%）。对于氮掺杂金红石，仅出现一个单深的杂化态（N 2p与O 2p混合），DFT+U同样增强其局域化（N 2p成分从38%增至48%）。两种方法一致表明，可见光活性增强源于N 2p孤立态位于VBM之上，而非N与O态混合导致的带隙窄化。

**3.2 DFT-1/2计算能带结构与氮掺杂效应**

DFT-1/2完美再现了纯相锐钛矿（3.61 eV）和金红石（3.15 eV）的实验带隙，与高成本准粒子方法结果一致。对于氮掺杂体系，通过双参数优化协议确定了氮的截断半径。在锐钛矿中，两个N诱导态（深态和浅态）的深态N 2p成分在DFT-1/2下为49.6%，与DFT+U（49.8%）接近。基本带隙（间接）从DFT的0.915 eV增至DFT-1/2的1.194 eV。在金红石中，深态N 2p成分在DFT-1/2下达到最大值51.5%，基本带隙从DFT的0.720 eV增至DFT-1/2的1.035 eV。三种方法给出的能带色散和分布趋势一致，但DFT-1/2提供了最接近实验的带隙值。

**3.3 非晶TiO₂中替代氮掺杂**

非晶模型（192原子）中Ti原子主要为六配位，但存在少量四、五、七配位，DFT+U采用加权平均U值（5.55 eV）。纯非晶相带隙：DFT 2.430 eV，DFT+U 2.615 eV，DFT-1/2 3.521 eV。氮掺杂在三种氧位（二、三、四配位）中，三配位氧位形成能最低（4.70 eV，富氧条件），且仅比二配位高0.051 eV，表明多种构型均可实现。氮掺杂引入多个深中隙态，类似锐钛矿，其N 2p成分随方法从DFT的32.9%增至DFT-1/2的42.9%。基本带隙（DFT 1.02 eV，DFT+U 1.173 eV，DFT-1/2 1.268 eV）趋势与晶体相一致，但DFT-1/2对非晶带隙的打开幅度（33%）小于晶体（约50%），符合实验趋势。

**总结与结论**

本研究重新审视了TiO₂中常见的替代氮（N_O）缺陷，在固定浓度（~0.78%）下，系统比较了标准DFT、DFT+U和DFT-1/2三种方法对锐钛矿、金红石和非晶相的影响。DFT-1/2方法能准确打开纯相带隙至接近实验值，且计算成本远低于准粒子方法。对于氮掺杂缺陷，研究人员开发了一种双参数优化协议，通过结合带边与隙态能量最大化以及基本带隙最小化，唯一确定了氮杂质的势函数截断半径。三种方法在缺陷态位置和组成上趋势一致，但DFT-1/2提供了最接近实验的定量描述。该工作为利用低成本的DFT-1/2方法进行缺陷工程计算提供了可靠协议，尤其适用于非晶TiO₂这一新兴光催化体系。