基于密度泛函理论（DFT）研究碳掺杂与外电场对二氧化钛（001）面电子结构及光学性能的协同影响

一、 研究作者、机构及发表信息

本研究的主要作者包括Cuihua Zhao、Dewei Huang和Jianhua Chen*（通讯作者）。所有作者均来自中国广西大学，具体隶属于广西大学资源、环境与材料学院，以及广西有色金属及特色材料加工重点实验室。此外，Jianhua Chen 也隶属于广西大学金属资源利用与环境保护创新中心。该项研究成果发表于Elsevier旗下的期刊《Journal of Materiomics》（材料学报），文章于2018年5月16日在线发表，最终收录于该期刊2018年的第4卷，第247至255页。文章的标题为“DFT study for combined influence of C-doping and external electric field on electronic structure and optical properties of TiO2 (001) surface”。

二、 学术背景与研究目的

本研究属于计算材料科学和半导体光催化领域。二氧化钛（TiO2）因其在光催化净化水与空气、自清洁表面、抗菌防雾等诸多领域的广泛应用，数十年来一直是研究热点。然而，纯TiO2作为光催化剂存在主要缺陷：其光响应范围主要在紫外区（因其带隙较宽，约3.2 eV），对太阳光中占大部分的可见光吸收率低，且光生电子与空穴复合率高，导致太阳能转换效率低下。为克服这些缺陷，研究者们尝试了多种方法，其中非金属掺杂（如氮、碳、硫掺杂）被认为是最有吸引力的策略之一，因其有望将TiO2的光吸收边红移至可见光区。已有研究表明，碳掺杂TiO2（C-TiO2）在可见光下表现出相对较好的光催化活性。

除了材料改性，外部电场也被认为是一种潜在的增效手段。电场可以促进光生电子与空穴的分离，增加半导体载流子浓度，并且基于Franz-Keldysh效应，电场还能导致光吸收边红移，从而增强对太阳光的利用。电场具有清洁、易得、强度和方向可调等优点。然而，将碳掺杂与外部电场两种手段结合，系统研究其对TiO2，特别是高活性（001）晶面电子结构和光学性能的协同影响，在微观机理层面的探索尚不充分。

因此，本研究旨在通过第一性原理计算，深入探究碳掺杂与外部电场对TiO2（001）表面的协同作用机制。具体目标包括：1) 构建稳定的C-TiO2（001）表面模型；2) 分析碳掺杂本身对TiO2（001）电子结构（能带、态密度、电荷布居）和光学性质的影响；3) 系统研究不同强度外电场（沿特定方向）对C-TiO2（001）电子结构和光学性质的调控作用；4) 从电子结构层面解释这种协同作用如何增强TiO2在可见光区的光催化活性，为开发高效光电催化材料提供理论依据。

三、 详细研究流程与方法

本研究是一项纯粹的理论计算工作，其核心流程基于密度泛函理论（Density Functional Theory, DFT）展开，不涉及实验样品制备与测试。所有计算均使用Materials Studio软件中的CASTEP（Cambridge Serial Total Energy Package）模块完成。

1. 计算方法与参数设定： * 交换关联泛函：采用广义梯度近似（GGA）下的PW91泛函。 * 电子关联修正：为更准确地描述TiO2中Ti的3d电子强关联效应，采用了DFT+U方法，U值设定为3.5 eV。 * 赝势：使用超软赝势（Ultrasoft Pseudopotentials）来描述价电子与离子实之间的相互作用。 * 截断能与K点网格：经过测试，平面波截断能设定为340 eV，布里渊区K点采样采用3×3×1的Monkhorst-Pack网格，以确保体系处于最稳定状态。 * 收敛标准：几何优化计算设置了严格的收敛容差：最大作用力小于0.08 eV/Å，最大位移小于0.002 Å，最大应力小于0.1 GPa，最大能量变化小于2.0×10⁻⁵ eV/atom。自洽场（SCF）收敛容差设定为2.0×10⁻⁶ eV/atom。

2. 计算模型构建： * 基底选择：选择锐钛矿相TiO2中光催化活性更高的（001）晶面作为研究对象。 * 表面模型：从优化的TiO2单胞出发，构建了（3×3×1）的超晶胞表面模型。为消除周期性边界条件中表面间的相互作用，在垂直于表面方向设置了20 Å的真空层。最终模型包含12层原子，其顶部和底部六层原子结构对称。 * 掺杂模型：研究仅考虑碳原子替代氧位点（O substitutional sites）的情况。通过计算不同氧位点（O12c, O33c, O43c, O63c）的替代形成能，确定最稳定的掺杂构型。计算表明，替代表面第一层二配位氧原子（O12c）的形成能最低（0.190 eV），因此选择此构型作为后续研究的C-TiO2（001）模型。在该模型中，碳原子浓度约为0.93%。

3. 研究步骤与流程： 整个研究遵循清晰的递进式流程： * 步骤一：纯TiO2（001）与C-TiO2（001）的基态性质计算。 首先对纯TiO2（001）表面模型进行几何优化和电子结构计算，作为基准。然后，对构建的C-TiO2（001）模型（C替代O12c位）进行几何优化，分析其结构弛豫（如键长变化）。随后，计算两者的能带结构、分波态密度（PDOS）、电子密度差以及Mulliken电荷布居，对比碳掺杂引入的影响。 * 步骤二：C-TiO2（001）在外电场下的电子结构计算。 在优化后的C-TiO2（001）模型上，沿表面法线方向（z方向）施加一系列不同强度的外部静电场（0.1 eV, 0.3 eV, 0.5 eV, 0.7 eV, 0.9 eV, 1.0 eV）。对每个电场强度下的体系进行自洽计算，分析其能带结构（区分自旋向上和自旋向下）和分波态密度（PDOS）的变化。作为对比，也计算了纯TiO2（001）在相同电场下的能带结构。 * 步骤三：光学性质计算。 基于上述各体系（纯TiO2、C-TiO2、不同电场下的C-TiO2）的基态电子结构，利用CASTEP模块计算它们的光学吸收谱。重点观察吸收边位置和可见光区吸收强度的变化。 * 步骤四：数据分析与机理阐释。 综合所有计算结果，分析碳掺杂和外电场各自及协同作用对TiO2（001）表面电子结构（带隙、态密度分布、电荷转移）和光学吸收特性的影响规律，并从物理机制（如杂质能级、能带弯曲、Franz-Keldysh效应、电子自旋极化等）角度进行解释，最终关联到光催化性能的提升。

四、 主要研究结果

1. 碳掺杂TiO2（001）的电子结构与光学性质： * 结构弛豫：碳原子替代O12c位后，两个C-Ti键的键长（2.061 Å 和 2.335 Å）比掺杂前的O-Ti键（~1.995-2.000 Å）更长。这归因于碳原子共价半径（0.77 Å）大于氧原子（0.73 Å），以及C-Ti键的离子性弱于O-Ti键。 * 电子结构：能带计算表明，C-TiO2（001）的带隙为2.29 eV，比纯TiO2（001）的2.83 eV显著减小。态密度分析揭示了原因：碳掺杂在费米能级附近引入了由C 2p轨道贡献的杂质能级，同时使导带整体向费米能级方向移动，共同导致带隙变窄。 * 成键与电荷：电子密度差和Mulliken电荷布居分析显示，C-Ti键的离子性弱于O-Ti键，共价性更强。碳原子获得的电荷（-0.43 e）少于氧原子（-0.72 e），钛原子失去的电荷也略少（1.31 e vs 1.42 e），这与氧的电负性（3.44）强于碳（2.55）相符。 * 光学性质：光学吸收谱显示，碳掺杂使TiO2（001）的光吸收边发生红移，向长波方向移动，表明其对可见光的吸收能力增强。

2. 外部电场对C-TiO2（001）电子结构的调控： * 带隙变化：当施加沿z方向的外部电场时，C-TiO2（001）的带隙进一步减小。在0.1 eV至0.5 eV的电场范围内，带隙随电场增强而持续减小（从2.28 eV降至0.76 eV），这主要源于导带的下移和展宽。 * 电子自旋极化：当电场强度超过0.5 eV（达到0.7 eV及以上）时，体系出现显著的电子自旋极化现象。此时，自旋向上和自旋向下的能带结构发生分离：自旋向上能带仍显示带隙（但数值继续变化），而自旋向下能带则表现为金属态（穿过费米能级）。随着电场从0.7 eV增至1.0 eV，自旋向上能带的带隙从2.15 eV减小至1.76 eV，同时自旋向下能带整体向高能方向移动。 * 态密度与成键变化：态密度分析表明，在低电场（≤0.5 eV）下，C 2s与Ti的s-p轨道在-3 eV至0 eV范围内的杂化较强。电场强度增至0.7 eV后，该能量区间的C 2p态密度急剧下降，C与Ti的杂化减弱。Mulliken电荷分析进一步证实，电场强度超过0.5 eV后，电荷从C原子向Ti原子转移的速度加快，C原子所带负电荷减少（从-0.43 e变为-0.20 e @ 1.0 eV），Ti原子所带正电荷增加（从1.31 e变为1.47 e @ 1.0 eV），表明C-Ti键的强度随强电场而减弱。 * 与纯TiO2的对比：纯TiO2（001）在相同电场下也呈现带隙减小和自旋极化现象，但电场对C-TiO2（001）的影响幅度更大，表明碳掺杂与电场的结合产生了更强的协同效应。

3. 外部电场对C-TiO2（001）光学性质的调控： * 光学吸收计算显示，随着外电场的增强，C-TiO2（001）在紫外区的吸收峰强度减弱，但吸收峰变得更尖锐且明显红移，导致吸收边向长波方向延伸。 * 在低电场（0.1-0.5 eV）下，吸收峰的变化相对较小。当电场达到并超过0.7 eV（即发生电子自旋极化时），吸收边的红移变得非常显著，在可见光区（400-800 nm）的光吸收强度大大增强。电场在0.7 eV至1.0 eV范围内引起的红移增幅最大。 * 这种红移现象可以用Franz-Keldysh效应解释：在外电场作用下，电子穿透势垒的几率增加，使得能量小于带隙的光子也能引发电子跃迁，从而导致吸收边向低能方向移动。同时，电场引起的能带倾斜有助于光生电子和空穴的分离，提高了载流子浓度。

五、 研究结论与价值

本研究通过系统的DFT计算，得出以下核心结论：碳掺杂与外加电场的协同作用能够显著调控并优化TiO2（001）表面的电子结构和光学性质，从而有效增强其在可见光区的光催化活性。具体而言： 1. 碳掺杂的单独作用：在TiO2（001）中引入碳杂质，通过引入杂质能级和使导带下移，将带隙从2.83 eV减小至2.29 eV，实现了光吸收边的初步红移。 2. 外电场的单独与协同作用：外部电场能进一步减小C-TiO2的带隙，并在强度超过阈值（0.5 eV）时诱导产生电子自旋极化，形成独特的半金属态电子结构。更重要的是，电场通过Franz-Keldysh效应和促进电荷分离，使C-TiO2的光吸收边发生更大幅度的红移，并显著增强其在可见光区的光吸收强度，这种效应在0.7-1.0 eV的强电场区间尤为突出。 3. 科学价值：该研究从原子和电子尺度深入揭示了碳掺杂与外部电场对TiO2光催化材料性能的协同调控机理，特别是明确了电场诱导电子自旋极化这一关键现象及其对光学性质的显著影响，为理解“光电催化”这一新兴领域的微观过程提供了重要的理论见解。 4. 应用价值：研究结果为设计和发展高效的光电催化材料体系提供了明确的理论指导。它提示，通过合理的元素掺杂（如碳掺杂）并结合可调的外部电场，可以动态、高效地拓宽TiO2等宽禁带半导体的光响应范围并提升其光生载流子分离效率，这对于开发高性能的太阳能转化与环境净化技术具有潜在的应用前景。

六、 研究亮点

1. 研究视角新颖：首次系统地从第一性原理层面，研究了碳掺杂与外部电场协同作用对TiO2（001）面电子结构及光学性能的影响，而非单独研究其中一种改性手段。
2. 发现了关键的电子结构转变：明确指出了外部电场存在一个强度阈值（~0.5 eV），超过此阈值后，C-TiO2体系会发生电子自旋极化，自旋向下通道呈现金属态，这是一种新颖的物态调控，对光生载流子的输运行为可能产生重要影响。
3. 揭示了清晰的构效关系与机制：研究建立了从原子掺杂→电子结构（带隙、态密度、电荷分布）变化→光学性质（吸收边红移）→光催化活性提升的完整逻辑链条，并用量子力学计算数据给予了详细佐证。明确将可见光吸收的增强归因于带隙减小和Franz-Keldysh效应的共同作用。
4. 方法严谨：研究采用了DFT+U方法修正强关联效应，测试并选择了合理的计算参数（截断能、K点、真空层厚度），并通过对比不同掺杂位点的形成能确定了最稳定的模型，确保了计算结果的可靠性。
5. 对比研究深入：不仅研究了C-TiO2在电场下的行为，还与纯TiO2在相同电场下的响应进行了对比，有力地证明了协同效应的存在和优势。

七、 其他有价值内容

研究中还指出，他们仅考虑了碳替代氧位点的情况，而前人研究（Muhich等）表明某些间隙碳位点可能更稳定。这提示未来研究可以进一步探索不同掺杂类型（间隙或共掺杂）与电场协同的效应。此外，研究模型基于固定的碳浓度（~0.93%），不同掺杂浓度的影响也是一个值得拓展的方向。这些点为后续更深入的理论与实验研究提供了线索。