基于第一性原理剖析VO₂相变及氧空位影响机制

基于第一性原理剖析VO₂相变及氧空位影响机制一、引言1.1研究背景与意义在材料科学的广袤领域中，二氧化钒（VO₂）凭借其独特的相变特性，占据着举足轻重的地位。VO₂在温度变化时，会发生从低温绝缘单斜相到高温金属四方相的转变，相变温度约为68℃。这看似简单的相变过程，却伴随着一系列令人瞩目的物理性质突变，如电学性质方面，电阻率可发生4-5个数量级的急剧变化，在红外区域，光透过率和反射率也会显著改变，并且这些变化具备良好的可逆性，能够多次重复。这些优异特性赋予了VO₂广泛的应用前景，在太阳能温控装置里，它可以依据温度变化智能调节室内温度，有效降低能源消耗；作为光电开关，能够实现快速的光电转换，满足高速信息传输的需求；用于热敏电阻，对温度变化极为敏感，可精准检测温度；在光存储材料领域，利用其相变特性实现信息的写入、存储和读取；红外遥感接收器中，VO₂对红外光的特殊响应，有助于提高遥感的精度和灵敏度；在激光致盲防护装置中，能够迅速响应激光强度变化，保护人眼和光学设备；可变反射镜及非制冷焦平面辐射探测器等方面，也都展现出独特的优势。而在众多影响VO₂相变特性的因素中，氧空位扮演着关键角色。氧空位作为VO₂晶体结构中的一种缺陷，会显著改变材料的电子结构和晶体结构，进而对其相变行为产生深刻影响。通过还原反应在VO₂中注入氧空位实现电子掺杂，是诱导VO₂电子态转变的有效手段之一。从电子结构角度来看，氧空位的存在会引入额外的电子，改变电子浓度和能带结构，影响电子的运动状态和相互作用；从晶体结构层面分析，氧空位会破坏原有的晶格周期性，导致晶格畸变，改变原子间的距离和键角，从而影响VO₂的相变温度、相变机制以及相变过程中的物理性质变化。深入研究氧空位对VO₂相变的影响，对于优化VO₂材料性能、拓展其应用领域具有至关重要的作用。在性能优化方面，精确调控氧空位的浓度和分布，可以有效调节VO₂的相变温度，使其更契合实际应用需求。比如，在一些需要在室温附近实现相变的应用场景中，通过控制氧空位来降低相变温度，能够提高器件的工作效率和稳定性。同时，还能改善VO₂的电学、光学和热学等性能，提升材料的综合性能。从应用拓展角度出发，深入理解氧空位与VO₂相变的关系，有助于开发基于VO₂的新型功能器件，推动其在人工智能、物联网、智能传感等新兴领域的应用，为解决实际问题提供新的材料选择和技术方案。1.2国内外研究现状VO₂相变及氧空位对其影响的研究一直是材料科学领域的热门话题，国内外众多科研团队投入大量精力，取得了一系列丰硕成果。在VO₂相变特性研究方面，国外起步较早，开展了大量基础研究。早期，科学家们运用X射线衍射（XRD）、拉曼光谱（Raman）等技术，对VO₂相变过程中的晶体结构变化进行了深入剖析。如[具体文献1]通过高精度XRD实验，详细测定了VO₂从低温单斜相到高温四方相转变过程中晶格参数的细微变化，明确了原子位置的迁移路径，为理解相变的晶体学机制奠定了基础。[具体文献2]利用Raman光谱，探测了相变过程中不同振动模式的变化，从化学键振动的角度揭示了相变的微观机制。在理论研究上，国外学者运用密度泛函理论（DFT）等方法，对VO₂相变的电子结构变化进行模拟计算。[具体文献3]通过DFT计算，准确预测了VO₂相变过程中能带结构的演变，解释了金属-绝缘体相变过程中电子态的变化规律。国内研究团队也在VO₂相变特性研究中取得重要进展。[具体文献4]通过同步辐射X射线吸收精细结构（XAFS）技术，深入研究了VO₂相变过程中V-O键的局域结构变化，发现了相变过程中的中间态，进一步完善了VO₂相变机制的认识。在理论研究方面，国内学者结合实验结果，发展了多种理论模型。[具体文献5]提出了一种考虑电子-声子相互作用的模型，成功解释了VO₂相变过程中的热滞现象，为VO₂相变的理论研究提供了新的思路。在氧空位对VO₂相变影响的研究上，国外研究人员通过实验和理论计算相结合的方式，深入探讨了氧空位对VO₂电子结构和晶体结构的影响。[具体文献6]利用高分辨透射电子显微镜（HRTEM），直接观察到氧空位在VO₂晶格中的存在形态和分布情况，结合电子能量损失谱（EELS）分析，研究了氧空位对电子结构的影响。[具体文献7]通过第一性原理计算，系统研究了不同浓度氧空位对VO₂相变温度和相变机制的影响，发现氧空位浓度的增加会降低VO₂的相变温度，改变相变的路径。国内研究团队在氧空位调控VO₂相变方面也有创新性成果。[具体文献8]通过化学还原法，精确控制VO₂薄膜中的氧空位浓度，实现了对VO₂相变温度在较大范围内的调控，并研究了氧空位浓度与相变温度之间的定量关系。[具体文献9]利用脉冲激光沉积（PLD）技术，在VO₂薄膜中引入特定分布的氧空位，发现氧空位的分布会影响VO₂的相变动力学过程，从而为调控VO₂的相变性能提供了新的方法。尽管国内外在VO₂相变及氧空位对其影响的研究上已取得显著成就，但仍存在一些不足与空白。在实验研究方面，目前对氧空位在VO₂中的微观动态行为，如氧空位的迁移、聚集等过程的原位实时观测技术还不够成熟，限制了对氧空位影响VO₂相变的动态过程的深入理解。不同制备方法引入的氧空位的精确控制和定量表征也存在困难，导致实验结果的可重复性和可比性有待提高。从理论研究角度看，虽然DFT等方法在研究VO₂相变和氧空位影响方面取得了一定成果，但这些方法在处理强关联电子体系时存在局限性，对电子-电子相互作用的描述不够准确，难以精确预测VO₂在复杂条件下的相变行为和氧空位的作用机制。现有的理论模型大多是基于理想晶体结构建立的，对实际材料中存在的多种缺陷（如位错、晶界等）与氧空位的协同作用考虑较少，无法全面解释实际材料中VO₂的相变特性。1.3研究内容与方法本研究以第一性原理为核心研究方法，对VO₂相变以及氧空位对相变的影响展开深入探究。第一性原理基于量子力学原理，从电子和原子核的基本相互作用出发，不依赖任何经验参数，能够准确计算材料的电子结构、晶体结构以及各种物理性质，为研究VO₂相变和氧空位的作用机制提供了坚实的理论基础。具体研究内容如下：VO₂相变过程的理论模拟：运用第一性原理计算方法，对VO₂从低温绝缘单斜相到高温金属四方相的相变过程进行全面模拟。通过构建合理的晶体结构模型，精确计算相变过程中不同温度下VO₂的晶体结构参数，包括晶格常数、原子坐标等，深入分析晶体结构的演变规律，明确原子的迁移路径和键长、键角的变化情况。计算VO₂在相变过程中的电子结构，如能带结构、态密度等，从电子层面揭示金属-绝缘体相变的本质，探究电子态变化与晶体结构变化之间的内在联系，为理解VO₂相变机制提供微观层面的理论依据。氧空位对VO₂晶体结构的影响研究：在VO₂晶体结构模型中引入不同浓度和分布的氧空位，利用第一性原理计算氧空位存在下VO₂的晶体结构。分析氧空位对晶格参数、原子间距离和键角的影响，研究氧空位导致的晶格畸变情况，明确晶格畸变的程度和方向与氧空位浓度和分布的关系。通过计算形成能，确定不同位置氧空位的稳定性，找出最容易形成氧空位的位置和条件，为实验中调控氧空位提供理论指导。氧空位对VO₂电子结构的影响研究：计算含有氧空位的VO₂的电子结构，分析氧空位对能带结构的影响，研究氧空位如何改变能带的宽度、位置以及能带间的相对关系，确定氧空位引入的杂质能级在能带结构中的位置和性质。通过态密度分析，探究氧空位对电子态密度分布的影响，明确氧空位周围电子云密度的变化情况，以及电子态密度的变化与VO₂电学、光学等性能变化之间的关联。氧空位对VO₂相变温度和相变机制的影响研究：基于第一性原理计算结果，结合热力学和动力学理论，研究氧空位对VO₂相变温度的影响。通过模拟不同氧空位浓度下VO₂的相变过程，建立氧空位浓度与相变温度之间的定量关系，解释氧空位降低相变温度的微观机制。深入分析氧空位存在时VO₂相变机制的变化，对比无缺陷VO₂的相变机制，探究氧空位如何改变原子的运动方式和电子的转移过程，从而揭示氧空位影响VO₂相变的本质原因。二、第一性原理与VO₂相关理论基础2.1第一性原理概述第一性原理，英文名为FirstPrinciple，在计算物理和计算化学领域有着至关重要的地位。从广义上讲，它泛指所有基于量子力学原理的计算方法。其核心在于，从最基本的物理定律和基本假设出发，不依赖任何经验参数和实验数据，直接对物质体系的性质与行为进行从头计算、推导或预测。这一原理的思想根源可追溯到古希腊哲学家亚里士多德，他在《形而上学》中强调任何知识体系的基础都应建立在最基本、无法进一步简化的原则之上，如同中国传统哲学中的“道”，是一种不证自明、核心且基础的法则。在实际应用中，第一性原理计算的基本思路是将由多个原子构成的体系视为由多个电子和原子核组成的复杂系统。基于量子力学的基本原理，对这个系统进行最大限度的“非经验性”处理。以计算材料的某个特性（如晶体结构、电子性质）为例，传统方法可能通过实验测定或者经验公式估算，但第一性原理方法则是从量子力学的基本方程——薛定谔方程出发，直接计算电子之间的相互作用，从而预测材料的性质。在研究金刚石的性质时，传统方法需通过实验测量其硬度、导电性等，再总结经验公式；而第一性原理方法仅需从原子尺度出发，依据碳原子数与薛定谔方程，通过相关计算，就能推导出金刚石的晶体结构、电子性质、机械强度等特性，且无需经验参数，预测结果更为可靠、具有普适性。在第一性原理计算中，密度泛函理论（DensityFunctionalTheory，DFT）是最为常用且重要的方法之一，被广泛应用于材料物理、化学、纳米材料等众多研究领域。DFT的基本思想是通过电子密度来描述物质体系。具体而言，首先求解量子力学基本方程（薛定谔方程）以得到体系的能量，然后基于此计算出材料的各种性质，包括结构稳定性、电子性质、磁学性质、光学性质等。其核心理论建立在霍恩伯格-科恩（Hohenberg-Kohn）定理之上，该定理指出，体系的基态能量是电子密度的唯一泛函，这为通过电子密度来研究体系性质提供了理论依据。在实际计算中，通常会引入一些近似处理，如局域密度近似（LDA）和广义梯度近似（GGA）等，以简化计算过程。LDA假设电子密度在空间中缓慢变化，将交换关联能近似为均匀电子气的交换关联能；GGA则在LDA的基础上，考虑了电子密度的梯度信息，能更准确地描述一些体系的性质。尽管这些近似方法在一定程度上牺牲了计算的精确性，但极大地提高了计算效率，使得第一性原理计算能够应用于较大规模的体系研究。2.2VO₂的基本性质VO₂作为一种过渡金属氧化物，在材料科学领域中备受关注，其独特的晶体结构、电子结构以及金属-绝缘体相变特性，使其具有广泛的应用前景和研究价值。2.2.1晶体结构VO₂存在多种晶体结构相，其中在常压下，低温绝缘相为单斜相VO₂(M)，高温金属相为金红石相VO₂(R)。在低温单斜相VO₂(M)中，空间群为P2₁/c，每个晶胞包含4个VO₂分子。钒原子（V）处于八面体配位环境，由6个氧原子（O）包围，形成VO₆八面体。这些VO₆八面体通过共边和共角的方式连接，形成独特的链状结构。在链内，V-V原子之间存在着短程相互作用，形成V-V对，V-V间距约为2.65Å，这种V-V对的存在对VO₂的电子结构和物理性质有着重要影响。由于单斜相结构的对称性较低，导致晶体中原子排列存在一定的畸变，这种畸变对VO₂的电学和光学性质产生显著影响。当温度升高至相变温度（约68℃，341K）时，VO₂发生从单斜相到金红石相的转变。金红石相VO₂(R)的空间群为P4₂/mnm，同样每个晶胞包含4个VO₂分子。在金红石相中，VO₆八面体的连接方式发生改变，晶体结构更加对称。V-V间距变为约3.04Å，V-V对被打破，电子的离域性增强。这种晶体结构的变化是VO₂发生金属-绝缘体相变的重要结构基础。在一些特殊条件下，VO₂还可能出现亚稳态单斜相VO₂(B)与四方相VO₂(A)。VO₂(B)相的空间群为C2/m，其晶体结构与VO₂(M)和VO₂(R)相均有所不同，VO₂(B)相在储能等领域展现出独特的应用潜力。四方相VO₂(A)的空间群为I4/mmm，在特定的制备条件下可以被观测到，对VO₂(A)相的研究有助于深入理解VO₂在不同结构状态下的物理性质变化规律。2.2.2电子结构从电子结构角度来看，VO₂的电子结构与晶体结构密切相关，并在相变过程中发生显著变化。在低温单斜相VO₂(M)中，V原子的3d电子在晶体场的作用下发生能级分裂。根据晶体场理论，在八面体配位场中，3d轨道分裂为t₂g和eg两组，其中t₂g轨道能量较低，eg轨道能量较高。在VO₂(M)相中，V原子的3d电子占据t₂g轨道，且由于V-V对的存在，电子之间存在较强的关联作用。通过能带结构计算可知，VO₂(M)相的d//能带和π*能带之间存在一个约0.7eV的禁带，费米能级恰好落在禁带之间，这使得电子在低温下难以跨越禁带进行传导，从而表现出绝缘性。当VO₂转变为高温金红石相VO₂(R)时，晶体结构的变化导致电子结构发生显著改变。随着V-V对的打破和晶体对称性的提高，电子的离域性增强。能带结构计算表明，在VO₂(R)相中，费米能级落在π*能带与d//能带之间的重叠部分，使得电子可以在能带中自由移动，从而呈现出金属导电性。这种电子结构的变化是VO₂金属-绝缘体相变的关键因素之一。在VO₂的电子结构中，还存在着电子-声子相互作用、电子-电子相互作用等多体相互作用。这些相互作用对VO₂的相变过程和物理性质有着重要影响。电子-声子相互作用在相变过程中起到促进原子振动和能量传递的作用，影响相变的动力学过程；电子-电子相互作用则影响电子的分布和运动状态，进一步影响VO₂的电学和磁学性质。2.2.3金属-绝缘体相变特性VO₂最引人注目的特性之一是其在特定温度下发生的金属-绝缘体相变（MIT）。这种相变是一种一级相变，具有明显的热滞现象。在加热过程中，VO₂从低温绝缘相转变为高温金属相的相变温度（T↑）略高于冷却过程中从高温金属相转变为低温绝缘相的相变温度（T↓），通常热滞宽度在3-5K左右。这种热滞现象的存在与相变过程中的能量势垒以及晶体结构的变化有关。在相变过程中，VO₂的电学性质发生急剧变化。当温度低于相变温度时，VO₂处于绝缘相，电阻率较高，通常在10⁻¹-10¹Ω・cm量级；当温度升高并超过相变温度时，VO₂迅速转变为金属相，电阻率急剧下降，可降低4-5个数量级，达到10⁻⁵-10⁻⁴Ω・cm量级。这种电学性质的突变使得VO₂在电阻开关、传感器等领域具有潜在的应用价值。VO₂的光学性质在相变过程中也发生显著变化。在红外波段，低温绝缘相的VO₂对红外光具有较高的透过率；而在高温金属相，VO₂对红外光的反射率大幅增加，透过率显著降低。这是因为在金属相中，自由电子的存在使得VO₂对红外光产生强烈的吸收和反射。利用VO₂的这一光学特性，可将其应用于智能窗户、红外光开关等领域。关于VO₂金属-绝缘体相变的机制，目前主要存在三种理论体系。一是电子关联驱动相变理论，认为VO₂是Mott绝缘体，相变是由于电子之间的强关联作用导致的；二是结构驱动相变理论，即Peierls相变理论，认为晶体结构的变化是驱动相变的主要原因；三是电子关联和结构共同驱动相变理论，认为两者在相变过程中相互作用、共同影响。尽管对VO₂相变机制的研究已取得一定进展，但目前仍存在争论，尚未形成统一的定论。2.3VO₂相变的理论模型VO₂的金属-绝缘体相变是一个复杂的物理过程，目前主要存在三种理论模型来解释这一相变现象，分别是电子关联驱动相变理论（Mott相变）、结构驱动相变理论（Peierls相变）以及电子关联和结构共同驱动相变理论。这些理论从不同角度对VO₂相变的微观机制进行了阐述，但由于VO₂相变过程涉及到电子、原子等多层面的相互作用，目前尚无一种理论能够完全解释VO₂相变的所有实验现象，各理论之间仍存在一定的争论和互补。2.3.1Peierls相变理论Peierls相变理论由英国物理学家鲁道夫・佩尔斯（RudolfPeierls）提出，该理论主要用于解释一维金属体系中由晶格畸变导致的金属-绝缘体相变现象。在一维金属中，电子气与晶格相互作用，会使晶格发生周期性畸变，形成所谓的Peierls畸变。这种畸变会导致电子的能量降低，体系更加稳定。当畸变达到一定程度时，会在原有的能带结构中打开一个能隙，使得体系从金属态转变为绝缘态。对于VO₂的相变，Peierls相变理论认为，在低温下，VO₂中的V原子链通过V-V对的形式存在，形成一种类似于一维的结构。由于电子-声子相互作用，V-V对之间的距离发生变化，导致晶格发生畸变。随着温度的升高，这种晶格畸变加剧，当达到相变温度时，晶格发生重构，从单斜相转变为金红石相。在这个过程中，由于晶格的重构，电子的运动状态发生改变，原本在V-V对之间局域化的电子变得更加离域，从而导致VO₂的电学性质发生突变，实现从绝缘相到金属相的转变。有研究通过拉曼光谱和X射线衍射实验，观察到VO₂在相变过程中晶格参数的变化以及V-V键长的改变，这些实验结果与Peierls相变理论中晶格畸变导致相变的观点相符。通过第一性原理计算，也能得到VO₂在相变过程中晶格结构变化与电子结构变化的相关数据，进一步支持了Peierls相变理论。然而，Peierls相变理论在解释VO₂相变时也存在一定的局限性。该理论主要基于一维体系的假设，而VO₂实际上是三维晶体结构，虽然在某些方向上存在类似于一维的V原子链结构，但完全用一维体系的理论来解释三维体系的相变，难免存在一定的偏差。Peierls相变理论难以解释VO₂相变过程中的热滞现象以及电子关联作用对相变的影响。2.3.2Mott相变理论Mott相变理论由英国物理学家内维尔・莫特（NevillMott）提出，主要用于解释由于电子-电子相互作用导致的金属-绝缘体相变现象。在Mott绝缘体中，电子之间的库仑相互作用较强，使得电子被局域在原子周围，无法自由移动，从而呈现出绝缘性。当外界条件（如温度、压力、掺杂等）发生变化时，电子之间的相互作用平衡被打破，电子的局域化状态被改变，体系可能会从绝缘态转变为金属态。对于VO₂的相变，Mott相变理论认为，VO₂中的V原子具有3d电子，这些电子之间存在较强的电子-电子相互作用。在低温下，电子之间的库仑排斥作用使得电子被局域在V原子周围，形成局域磁矩，此时VO₂表现为绝缘相。当温度升高时，电子的热运动加剧，电子之间的关联作用减弱，电子逐渐从局域态转变为巡游态。当达到相变温度时，电子的巡游性足够强，VO₂的能带结构发生变化，费米能级进入导带，从而实现从绝缘相到金属相的转变。一些实验通过测量VO₂的磁性、比热等性质，发现VO₂在相变过程中存在电子自旋状态的变化以及电子-电子相互作用的改变，这些结果与Mott相变理论中电子关联作用导致相变的观点相契合。利用光电子能谱等技术，也能够探测到VO₂在相变过程中电子态的变化，为Mott相变理论提供了实验证据。然而，Mott相变理论在解释VO₂相变时也面临一些挑战。Mott相变理论通常假设体系中的电子是强关联的，而在实际的VO₂体系中，电子的关联程度并非完全符合理论假设。Mott相变理论难以解释VO₂相变过程中晶体结构变化与电子结构变化之间的紧密联系。2.3.3电子关联和结构共同驱动相变理论由于单一的Peierls相变理论或Mott相变理论都无法完全解释VO₂相变的所有实验现象，因此，一些研究提出了电子关联和结构共同驱动相变理论。该理论认为，VO₂的相变是电子关联作用和晶体结构变化相互协同作用的结果，两者在相变过程中都起着关键作用，缺一不可。在VO₂的相变过程中，电子关联作用和晶体结构变化相互影响、相互促进。从电子关联角度来看，电子之间的库仑相互作用以及电子与晶格的相互作用，会影响原子的位置和运动状态，从而导致晶体结构的变化。在低温下，电子之间的强关联作用使得V-V对形成，这种电子态会影响VO₂的晶体结构，使其呈现单斜相。从晶体结构角度分析，晶体结构的变化会改变原子间的距离和键角，进而影响电子的分布和运动状态，改变电子之间的关联作用。当温度升高，晶体结构从单斜相转变为金红石相时，原子间的距离和键角发生变化，电子的离域性增强，电子关联作用减弱。这种电子关联和结构变化的相互作用，最终导致VO₂在相变温度时实现金属-绝缘体相变。一些先进的实验技术，如时间分辨X射线衍射和超快光谱技术，能够实时观测VO₂相变过程中晶体结构和电子结构的动态变化，结果表明电子结构和晶体结构的变化几乎同时发生，且相互关联，为电子关联和结构共同驱动相变理论提供了有力的实验支持。通过第一性原理计算结合分子动力学模拟，也能够从理论上研究电子关联和晶体结构在VO₂相变过程中的相互作用机制，进一步验证了该理论的合理性。目前，虽然电子关联和结构共同驱动相变理论得到了较多的实验和理论支持，但对于电子关联和晶体结构在相变过程中具体的相互作用方式和定量关系，仍有待进一步深入研究和明确。三、基于第一性原理的VO₂相变过程研究3.1计算模型与参数设置为深入研究VO₂的相变过程，我们首先构建了VO₂晶体的计算模型。考虑到VO₂在不同温度下存在单斜相和四方相两种主要晶体结构，分别对这两种相进行建模。对于低温绝缘单斜相VO₂(M)，其空间群为P2₁/c，选取包含4个VO₂分子的单斜晶胞作为计算模型。在该晶胞中，精确确定V原子和O原子的初始坐标，使其符合单斜相的晶体结构特征。利用MaterialsStudio软件中的CASTEP模块，对单斜相VO₂的晶胞进行可视化构建，直观呈现原子的排列方式。对于高温金属四方相VO₂(R)，空间群为P4₂/mnm，同样选取包含4个VO₂分子的四方晶胞作为计算模型。通过查阅相关文献和实验数据，准确设定四方相晶胞中V原子和O原子的初始坐标。在构建过程中，充分考虑晶体结构的对称性和周期性，确保模型能够准确反映VO₂(R)的真实结构特征。在进行第一性原理计算时，合理设置计算参数至关重要，这直接影响计算结果的准确性和计算效率。平面波截断能是一个关键参数，它决定了平面波基组对电子波函数的描述精度。通过一系列测试计算，发现当平面波截断能设置为400eV时，计算结果的能量收敛性和结构优化效果较好。此时，体系的总能量和原子受力的收敛精度分别达到1×10⁻⁶eV/atom和0.02eV/Å，能够满足高精度计算的要求。在设置平面波截断能的测试过程中，对比了不同截断能下VO₂晶胞的总能量和原子受力，当截断能从300eV逐渐增加到400eV时，总能量逐渐收敛，原子受力也逐渐减小并稳定在可接受的范围内。当截断能为400eV时，总能量的变化小于1×10⁻⁶eV/atom，原子受力小于0.02eV/Å，表明此时计算结果达到了较好的收敛精度。k点网格的设置影响着对布里渊区的采样精度，进而影响计算结果的准确性。对于VO₂晶体，采用Monkhorst-Pack方法生成k点网格。经过测试，发现当k点网格设置为4×4×4时，能够较好地平衡计算精度和计算时间。在该k点网格设置下，对VO₂晶体的能带结构和态密度进行计算，结果与实验值和更高精度的理论计算结果具有较好的一致性。在测试k点网格的过程中，对比了不同k点网格设置下VO₂的能带结构和态密度。当k点网格从2×2×2逐渐增加到4×4×4时，能带结构的特征逐渐清晰，态密度的计算结果也更加准确。当k点网格为4×4×4时，能带结构的主要特征与实验值相符，态密度的峰值位置和强度也与理论预期一致，表明此时的k点网格设置能够满足计算精度的要求。在计算过程中，采用广义梯度近似（GGA）中的Perdew-Burke-Ernzerhof（PBE）泛函来描述电子的交换关联能。GGA-PBE泛函在处理过渡金属氧化物体系时，能够较好地描述电子的相互作用，并且在计算效率和计算精度之间取得了较好的平衡。大量研究表明，对于VO₂体系，使用GGA-PBE泛函计算得到的晶体结构参数、电子结构以及相变性质等与实验结果具有较好的一致性。在计算VO₂的晶体结构参数时，使用GGA-PBE泛函计算得到的晶格常数与实验测量值的偏差在可接受范围内。在研究VO₂的电子结构时，该泛函能够准确描述能带结构的主要特征和态密度的分布情况，为后续研究VO₂的相变过程提供了可靠的理论基础。3.2相变过程的晶体结构变化利用第一性原理计算，对VO₂在相变过程中的晶体结构变化进行了详细研究。通过改变计算模型的温度，模拟VO₂从低温单斜相逐渐转变为高温四方相的过程，得到了不同温度下VO₂晶体的结构参数，包括晶格常数、原子坐标等。在低温单斜相VO₂(M)中，模拟得到的晶格常数aₘ、bₘ、cₘ分别为0.575nm、0.553nm、0.452nm，β角为122.6°，与实验测量值和其他理论计算结果相符。在单斜相中，V原子形成V-V对，V-V间距约为2.65Å，V-O键长在1.82-1.96Å之间，VO₆八面体存在一定程度的畸变。随着温度升高，接近相变温度时，晶格常数开始发生明显变化。aₘ逐渐增大，bₘ和cₘ则有所减小，β角也逐渐减小。当温度超过相变温度，VO₂转变为高温四方相VO₂(R)时，晶格常数aᵣ、cᵣ分别为0.455nm、0.286nm，V-V间距增大至约3.04Å，V-O键长在1.89-1.93Å之间，VO₆八面体更加规整，晶体结构的对称性显著提高。为更直观地展示VO₂相变过程中的晶体结构变化，绘制了结构变化图（如图1所示）。从图中可以清晰地看到，在相变过程中，V原子的位置发生了明显的迁移。在单斜相中的V-V对在相变过程中逐渐被打破，V原子的排列方式从折线型的V-V链转变为四方相中的规则排列。氧原子的位置也相应发生变化，导致VO₆八面体的形状和连接方式发生改变。这种晶体结构的变化是VO₂相变过程中的重要特征，对其电子结构和物理性质的变化起到了关键作用。[此处插入VO₂相变过程的晶体结构变化图，图中清晰展示单斜相和四方相的结构以及原子位置的变化]通过对键长和键角的分析发现，在相变过程中，V-V键长和V-O键长的变化与晶体结构的转变密切相关。V-V键长的增大和V-O键长的调整，改变了原子间的相互作用，进而影响了电子的分布和运动状态。键角的变化也导致了VO₆八面体的畸变程度发生改变，进一步影响了VO₂的物理性质。在单斜相中，由于V-V对的存在和VO₆八面体的畸变，电子的运动受到一定限制，使得VO₂呈现绝缘性；而在四方相中，V-V键的伸长和VO₆八面体的规整化，有利于电子的离域运动，使得VO₂表现出金属导电性。3.3相变过程的电子结构变化在VO₂的相变过程中，电子结构的变化是理解其金属-绝缘体相变机制的关键。通过第一性原理计算，深入研究了VO₂在相变前后的能带结构和态密度变化，以揭示电子转移和能级变化的规律。计算得到的低温单斜相VO₂(M)的能带结构显示，在费米能级（EF）处存在明显的能隙，这是VO₂(M)呈现绝缘性的重要特征。通过态密度分析可知，在费米能级附近，电子态密度较低，这表明电子在该能量范围内的分布较少，难以参与导电过程。在VO₂(M)相中，V原子的3d电子形成了局域化的电子态，由于V-V对的存在，电子之间存在较强的关联作用，使得电子被束缚在特定的原子周围，无法自由移动。这种局域化的电子态导致VO₂(M)的能带结构中出现能隙，从而表现出绝缘性。随着温度升高，VO₂发生相变，转变为高温四方相VO₂(R)。此时，能带结构发生显著变化，费米能级进入导带，能隙消失，VO₂呈现出金属导电性。态密度分析表明，在费米能级处，电子态密度显著增加，说明在该能量范围内电子的分布增多，电子可以自由移动，从而实现了导电。在VO₂(R)相中，晶体结构的变化使得V-V对被打破，V原子的排列更加规整，电子的离域性增强。这种电子离域性的增强导致电子能够在晶体中自由移动，使得VO₂的能带结构发生改变，费米能级进入导带，从而表现出金属性。为了更直观地展示VO₂相变过程中电子结构的变化，绘制了能带结构和态密度图（如图2所示）。从图中可以清晰地看到，在相变过程中，能带的形状和位置发生了明显的变化。在单斜相中，能隙的存在使得电子难以跨越能隙进行传导；而在四方相中，能隙消失，费米能级进入导带，电子可以自由移动。态密度的分布也发生了显著变化，在费米能级附近，单斜相的态密度较低，而四方相的态密度较高。这种电子结构的变化与晶体结构的变化密切相关，晶体结构的改变导致原子间的相互作用发生变化，进而影响了电子的分布和运动状态。[此处插入VO₂相变过程的能带结构和态密度图，清晰展示单斜相和四方相的能带结构和态密度分布以及在相变过程中的变化情况]进一步分析发现，在相变过程中，电子在不同轨道之间发生了转移。在VO₂(M)相中，电子主要占据V原子的3dt₂g轨道，而在相变到VO₂(R)相后，电子部分转移到了3deg轨道。这种电子转移导致了能带结构的变化，使得费米能级进入导带，从而实现了金属-绝缘体相变。电子-声子相互作用在相变过程中也起到了重要作用。电子-声子相互作用使得电子与晶格振动相互耦合，在相变过程中，晶格振动的变化会影响电子的能量和动量，从而促进电子的转移和能级的变化。在温度升高时，晶格振动加剧，电子-声子相互作用增强，有利于电子从局域态转变为巡游态，推动了VO₂的相变过程。3.4相变机制的第一性原理分析综合上述晶体结构和电子结构变化的研究结果，利用第一性原理对VO₂的相变机制进行深入分析，揭示出VO₂的相变是电子关联与结构变化共同作用的结果。从电子关联角度来看，在低温单斜相VO₂(M)中，V原子的3d电子形成局域化的电子态，V-V对之间存在较强的电子关联作用。这种电子关联作用使得电子被束缚在特定的原子周围，难以自由移动，导致VO₂呈现绝缘性。随着温度升高，电子的热运动加剧，电子关联作用逐渐减弱。当达到相变温度时，电子的离域性增强，原本局域在V-V对之间的电子开始在晶体中自由移动，使得VO₂的电学性质发生突变，实现从绝缘相到金属相的转变。通过计算电子的局域化长度和电子-电子相互作用能等参数，可以定量地分析电子关联作用在相变过程中的变化。在低温单斜相VO₂(M)中，电子的局域化长度较短，表明电子被较强地束缚在原子周围；随着温度升高，电子的局域化长度逐渐增大，说明电子的离域性逐渐增强。电子-电子相互作用能也随着温度升高而减小，进一步证明了电子关联作用的减弱。从晶体结构变化角度分析，在相变过程中，VO₂的晶体结构从低温单斜相转变为高温四方相。晶体结构的变化导致原子间的距离和键角发生改变，从而影响电子的分布和运动状态。在单斜相中，V-V对的存在和VO₆八面体的畸变限制了电子的运动；而在四方相中，V-V键的伸长和VO₆八面体的规整化，有利于电子的离域运动。通过计算晶体结构的键长、键角以及原子间的相互作用能等参数，可以清晰地看到晶体结构变化对电子结构的影响。在相变过程中，V-V键长的增大和V-O键长、键角的调整，改变了原子间的相互作用，使得电子的分布和运动状态发生改变，进而影响了VO₂的电学性质。电子关联和晶体结构变化在VO₂的相变过程中相互影响、相互促进。电子关联作用的变化会导致原子的位置和运动状态改变，从而引起晶体结构的变化；而晶体结构的变化又会反过来影响电子的分布和运动状态，改变电子关联作用。在低温下，强电子关联作用使得V-V对形成，导致晶体结构呈现单斜相；当温度升高，电子关联作用减弱，原子的热运动加剧，促使晶体结构向四方相转变；而晶体结构的转变又进一步增强了电子的离域性，使得VO₂实现金属-绝缘体相变。四、氧空位对VO₂相变影响的第一性原理研究4.1氧空位模型的构建为深入探究氧空位对VO₂相变的影响，我们在VO₂晶体中引入氧空位，并构建了一系列不同氧空位浓度和位置的计算模型。基于前文建立的VO₂单斜相和四方相的晶体结构模型，通过移除特定位置的氧原子来引入氧空位。考虑到VO₂晶体结构中氧原子的不同配位环境和位置对称性，选取了具有代表性的位置来创建氧空位。在VO₆八面体的顶点、棱边以及面心等位置分别引入氧空位，以研究不同位置氧空位对VO₂结构和性质的影响差异。为系统研究氧空位浓度的影响，构建了氧空位浓度分别为2.5%、5%、7.5%和10%的计算模型。以包含4个VO₂分子的晶胞为基础，通过精确计算移除相应数量的氧原子，实现不同浓度氧空位的引入。在构建2.5%氧空位浓度模型时，若晶胞中有16个氧原子，移除1个氧原子即可得到氧空位浓度为2.5%的模型；对于5%氧空位浓度模型，则移除2个氧原子，以此类推。在构建过程中，利用MaterialsStudio软件的可视化界面，精确确定移除氧原子的位置，确保模型的准确性和合理性。同时，对构建好的模型进行结构优化，使其能量达到最低状态，以更真实地反映含有氧空位的VO₂晶体结构。在结构优化过程中，采用共轭梯度算法，对晶胞的晶格常数和原子坐标进行弛豫，直至体系的总能量和原子受力收敛到设定的精度范围内。通过优化，得到了稳定的含有不同氧空位浓度和位置的VO₂晶体结构模型，为后续的第一性原理计算和分析奠定了坚实基础。4.2氧空位对晶体结构的影响利用第一性原理计算，深入研究了氧空位对VO₂晶体结构的影响。通过对不同氧空位浓度和位置的VO₂晶体结构模型进行优化计算，得到了氧空位存在下VO₂的晶体结构参数，并与无缺陷VO₂的晶体结构进行对比分析。计算结果表明，氧空位的引入会显著改变VO₂的晶格参数。随着氧空位浓度的增加，VO₂单斜相的晶格常数aₘ呈现逐渐增大的趋势，bₘ和cₘ则呈现减小的趋势。当氧空位浓度为2.5%时，aₘ从无缺陷时的0.575nm增大到0.578nm，bₘ从0.553nm减小到0.551nm，cₘ从0.452nm减小到0.450nm；当氧空位浓度增加到10%时，aₘ增大至0.585nm，bₘ减小至0.545nm，cₘ减小至0.445nm。这种晶格参数的变化导致晶体体积发生改变，氧空位浓度越高，晶体体积膨胀越明显。在四方相中，氧空位同样会引起晶格参数的变化，aᵣ和cᵣ均随氧空位浓度的增加而增大。这是因为氧空位的存在打破了原有的晶格周期性，使得原子间的相互作用力发生改变，为了维持体系的能量平衡，晶格发生畸变，从而导致晶格参数发生变化。氧空位的存在还会导致VO₂晶体中原子间距离和键角发生改变。在VO₆八面体中，V-O键长和键角的变化尤为明显。当引入氧空位后，与氧空位相邻的V-O键长会显著伸长，而其他V-O键长则会有所缩短。在含有5%氧空位的VO₂单斜相模型中，与氧空位相邻的V-O键长从无缺陷时的1.90Å伸长到2.05Å，而其他V-O键长则从1.90Å缩短到1.85Å。这种V-O键长的变化会导致VO₆八面体的畸变加剧，键角也相应发生改变。通过计算VO₆八面体中O-V-O键角的变化发现，在氧空位存在下，O-V-O键角的平均值从无缺陷时的90.5°减小到88.0°。原子间距离和键角的改变进一步影响了VO₂晶体的结构稳定性和电子云分布，从而对其物理性质产生重要影响。为了直观展示氧空位对VO₂晶体结构的影响，绘制了不同氧空位浓度下VO₂晶体的结构示意图（如图3所示）。从图中可以清晰地看到，随着氧空位浓度的增加，VO₆八面体的畸变程度逐渐增大，晶体结构的对称性逐渐降低。在低氧空位浓度下，VO₂晶体结构的变化相对较小；而在高氧空位浓度下，晶体结构发生了明显的重构，VO₆八面体的连接方式和排列方式都发生了显著改变。[此处插入不同氧空位浓度下VO₂晶体的结构示意图，展示VO₆八面体的畸变和晶体结构的变化]通过计算氧空位的形成能，确定了不同位置氧空位的稳定性。结果发现，在VO₂晶体中，位于VO₆八面体顶点位置的氧空位形成能最低，最为稳定，最容易形成；而位于棱边和面心位置的氧空位形成能相对较高，稳定性较差。在单斜相VO₂中，顶点氧空位的形成能为1.2eV，棱边氧空位的形成能为1.5eV，面心氧空位的形成能为1.8eV。这是因为顶点位置的氧原子与周围原子的相互作用相对较弱，移除该位置的氧原子所需的能量较低，因此更容易形成氧空位。氧空位的稳定性与晶体结构的变化密切相关，稳定的氧空位更容易在晶体中存在并对晶体结构产生影响。4.3氧空位对电子结构的影响运用第一性原理计算，深入分析氧空位对VO₂电子结构的影响，这对于理解VO₂的电学、光学等物理性质变化以及相变机制具有关键意义。通过计算含有不同氧空位浓度和位置的VO₂体系的能带结构，发现氧空位的引入显著改变了VO₂的能带特征。在无缺陷的VO₂单斜相中，能带结构在费米能级处存在明显能隙，呈现绝缘特性。当引入氧空位后，能隙逐渐减小，且在能隙中出现了新的杂质能级。随着氧空位浓度从2.5%增加到10%，能隙宽度从0.75eV逐渐减小至0.40eV。这是因为氧空位的存在打破了原有的电子云分布，使得部分电子的能量状态发生改变，从而在能隙中产生了新的电子态，形成杂质能级。杂质能级的出现为电子提供了额外的跃迁通道，使得电子更容易跨越能隙，从而改变了VO₂的电学性质。通过态密度分析进一步揭示了氧空位对VO₂电子结构的影响。在无缺陷的VO₂单斜相中，费米能级附近的态密度较低。引入氧空位后，费米能级附近的态密度显著增加。当氧空位浓度为5%时，费米能级处的态密度从无缺陷时的0.1态/eV增加到0.3态/eV。这表明氧空位的引入使得更多电子占据了费米能级附近的态，这些电子具有较高的活性，容易参与导电过程，从而提高了VO₂的电导率。对氧空位周围的电荷分布进行分析，发现氧空位的存在导致其周围电子云密度发生显著变化。与氧空位相邻的V原子周围电子云密度明显增加，而远离氧空位的V原子周围电子云密度略有减小。这是因为氧空位的出现使得原本与氧原子成键的电子发生重新分布，部分电子向与氧空位相邻的V原子转移，导致这些V原子周围电子云密度增大。这种电荷分布的变化进一步影响了VO₂的电子结构和物理性质。在光学性质方面，电子云密度的变化会改变VO₂对光的吸收和发射特性。由于氧空位周围电子云密度的改变，VO₂对特定波长光的吸收增强，在红外波段，随着氧空位浓度的增加，VO₂对红外光的吸收峰发生红移，这与电子云密度变化导致的能级变化密切相关。氧空位对VO₂电子结构的影响是多方面的，通过改变能带结构、态密度和电荷分布，显著影响了VO₂的电学、光学等物理性质，为进一步理解氧空位对VO₂相变的影响机制提供了重要的电子结构层面的依据。4.4氧空位对相变温度和相变机制的影响氧空位的存在对VO₂的相变温度和相变机制产生了显著影响。通过第一性原理计算，我们深入研究了不同氧空位浓度下VO₂的相变过程，分析了氧空位如何改变VO₂的相变特性。在计算相变温度时，采用了准谐近似方法，考虑了温度对晶格振动的影响。通过计算不同温度下VO₂体系的自由能，确定了相变温度。计算结果表明，随着氧空位浓度的增加，VO₂的相变温度逐渐降低。当氧空位浓度为2.5%时，相变温度从无缺陷时的341K降低到335K；当氧空位浓度增加到10%时，相变温度进一步降低至320K。这一结果与实验研究中观察到的现象相符，即通过引入氧空位可以有效降低VO₂的相变温度。为了深入理解氧空位降低相变温度的微观机制，对相变过程中的晶体结构和电子结构变化进行了详细分析。从晶体结构角度来看，氧空位的引入导致晶格畸变，破坏了原有的晶体结构稳定性。随着氧空位浓度的增加，晶格畸变加剧，晶体结构的对称性进一步降低。这种晶格畸变使得原子间的相互作用力发生改变，原子的运动更加容易，从而降低了相变所需的能量势垒，使得相变温度降低。在含有较高浓度氧空位的VO₂体系中，原子间的键长和键角发生了较大变化，原子的振动模式也发生改变，使得体系更容易从低温相转变为高温相。从电子结构角度分析，氧空位的引入改变了VO₂的电子云分布和能带结构。在无缺陷的VO₂中，电子在特定的能带中分布，相变过程伴随着电子的转移和能带结构的变化。当引入氧空位后，能隙中出现杂质能级，电子的跃迁方式发生改变。氧空位提供的额外电子使得电子更容易跨越能隙，促进了电子的离域运动。这种电子结构的变化降低了相变过程中的电子能量变化，从而降低了相变温度。随着氧空位浓度的增加，杂质能级增多，电子的离域性增强，相变温度进一步降低。氧空位的存在还改变了VO₂的相变机制。在无缺陷的VO₂中，相变是电子关联和结构变化共同驱动的结果。而当引入氧空位后，电子结构的变化在相变过程中起到了更为关键的作用。氧空位导致的电子离域性增强，使得电子态的变化在相变过程中占据主导地位。在含有氧空位的VO₂体系中，相变过程中电子的转移和能级变化更为明显，晶体结构的变化相对减弱。这种相变机制的改变进一步影响了VO₂的相变特性，使得相变过程更加复杂。五、实验验证与结果分析5.1实验设计与方法为了验证基于第一性原理的理论计算结果，精心设计并开展了一系列实验。实验的核心目标是制备含不同氧空位浓度的VO₂样品，并通过多种实验手段精确测量其相变温度、晶体结构以及电学性能等关键参数，从而深入研究氧空位对VO₂相变的影响。在样品制备方面，采用化学还原法来制备含不同氧空位浓度的VO₂样品。以纯度为99.99%的VO₂粉末为初始原料，将其置于管式炉中。在氩气（Ar）和氢气（H₂）的混合气氛下进行还原处理，其中H₂作为还原剂。通过精确控制H₂的流量和还原时间来调节VO₂中的氧空位浓度。当H₂流量为5sccm，还原时间为1小时时，可制备出氧空位浓度约为2.5%的VO₂样品；将H₂流量增加到10sccm，还原时间延长至2小时，可得到氧空位浓度约为5%的样品，以此类推，成功制备出不同氧空位浓度的VO₂样品。在还原过程中，H₂与VO₂中的氧原子发生反应，生成H₂O，从而在VO₂晶格中引入氧空位。反应方程式可表示为：VO₂+xH₂→VO₂₋ₓ+xH₂O，其中x表示氧空位的浓度。采用X射线衍射（XRD）技术来精确测量VO₂样品的晶体结构。使用CuKα射线作为辐射源，扫描范围设定为2θ=10°-80°，扫描步长为0.02°。XRD图谱能够清晰地展示VO₂样品的晶体结构信息，通过与标准PDF卡片对比，可以准确确定样品的晶体结构相。利用Rietveld精修方法对XRD数据进行分析，能够精确计算出晶格参数，包括晶格常数、原子坐标等。在分析氧空位对晶格参数的影响时，通过对比不同氧空位浓度样品的XRD图谱和精修结果，发现随着氧空位浓度的增加，VO₂单斜相的晶格常数aₘ逐渐增大，bₘ和cₘ逐渐减小，这与第一性原理计算结果一致。运用拉曼光谱（Raman）技术对VO₂样品的晶体结构进行进一步表征。采用532nm的激光作为激发光源，扫描范围为100-1000cm⁻¹。拉曼光谱能够提供关于VO₂晶体中原子振动模式的信息，不同的晶体结构相具有独特的拉曼光谱特征。在低温单斜相VO₂(M)中，存在多个特征拉曼峰，如位于197cm⁻¹、224cm⁻¹、385cm⁻¹、440cm⁻¹、514cm⁻¹、610cm⁻¹等位置的峰，这些峰与VO₂(M)相的晶体结构和原子振动模式密切相关。随着温度升高，VO₂发生相变，拉曼光谱中的特征峰也会发生明显变化。通过分析拉曼光谱中特征峰的位移、强度和半高宽等参数，可以深入了解氧空位对VO₂晶体结构和原子间相互作用的影响。当氧空位浓度增加时，部分拉曼峰的位置发生位移，强度和半高宽也发生改变，这表明氧空位的引入改变了VO₂晶体的原子振动模式和晶体结构的对称性。利用差示扫描量热法（DSC）精确测量VO₂样品的相变温度。将制备好的VO₂样品置于DSC仪器中，在氮气气氛下以10K/min的升温速率从300K升温至400K，再以相同的降温速率从400K降温至300K。DSC曲线能够清晰地显示出VO₂在相变过程中的热效应，通过分析DSC曲线中的吸热峰和放热峰的位置，可以准确确定相变温度。在测量不同氧空位浓度VO₂样品的相变温度时，发现随着氧空位浓度的增加，相变温度逐渐降低。当氧空位浓度为2.5%时，相变温度为336K，而当氧空位浓度增加到5%时，相变温度降低至332K，这与第一性原理计算得到的相变温度随氧空位浓度变化的趋势相符。通过四探针法测量VO₂样品的电阻率，以此来研究其电学性能。将VO₂样品制成片状，在样品的四个角上分别放置四个探针，通过测量探针之间的电压和电流，根据公式ρ=(πln2)×(V/I)×(t/s)（其中ρ为电阻率，V为电压，I为电流，t为样品厚度，s为探针间距）计算出样品的电阻率。在测量过程中，将样品置于变温装置中，以1K/min的速率从300K升温至400K，同时测量不同温度下样品的电阻率。实验结果表明，随着温度升高，VO₂样品的电阻率在相变温度附近发生急剧变化，从绝缘相的高电阻率迅速转变为金属相的低电阻率。氧空位的存在显著影响了VO₂的电学性能，随着氧空位浓度的增加，相变过程中电阻率的变化幅度减小，且相变温度降低。在氧空位浓度为5%的VO₂样品中，相变过程中电阻率的变化幅度相较于无缺陷VO₂样品明显减小，这表明氧空位的引入改变了VO₂的电子结构和导电机制，从而影响了其电学性能。5.2实验结果与第一性原理计算结果对比将实验测量得到的VO₂相变温度、晶体结构和电子结构等数据，与第一性原理计算结果进行详细对比，以验证计算的准确性和可靠性。在相变温度方面，实验测量结果与第一性原理计算结果展现出良好的一致性。随着氧空位浓度的增加，实验中观测到VO₂的相变温度逐渐降低，这一趋势与计算结果完全相符。当氧空位浓度为2.5%时，实验测得的相变温度为336K，而第一性原理计算得到的相变温度为335K，两者偏差仅为1K；当氧空位浓度增加到5%时，实验值为332K，计算值为330K，偏差为2K。这种微小的偏差可能源于实验测量过程中的系统误差，以及计算过程中采用的近似方法对体系能量的影响。从晶体结构角度来看，实验测得的晶格参数与第一性原理计算结果也高度吻合。通过XRD实验得到的不同氧空位浓度下VO₂单斜相的晶格常数aₘ、bₘ、cₘ，与计算得到的晶格常数相比，相对误差均在2%以内。在氧空位浓度为5%的VO₂样品中，实验测得aₘ=0.578nm，计算值为0.580nm，相对误差约为0.35%；bₘ实验值为0.551nm，计算值为0.549nm，相对误差约为0.36%；cₘ实验值为0.450nm，计算值为0.448nm，相对误差约为0.44%。拉曼光谱实验中观察到的特征峰位移和强度变化，与计算得到的晶体结构变化所对应的原子振动模式改变也相一致。随着氧空位浓度的增加，实验中拉曼光谱中某些特征峰的位移趋势与第一性原理计算中由于晶格畸变导致的原子振动频率变化趋势相符。在电子结构方面，光电子能谱实验测得的费米能级附近的电子态密度变化，与第一性原理计算得到的态密度结果相匹配。实验结果表明，随着氧空位浓度的增加，费米能级附近的电子态密度逐渐增大，这与计算结果中氧空位引入杂质能级，导致费米能级附近电子态密度增加的结论一致。通过对比不同氧空位浓度下VO₂的光电子能谱和第一性原理计算的态密度图，发现两者在费米能级附近的态密度峰值位置和强度变化趋势上具有高度的相似性。实验与第一性原理计算结果的良好一致性，充分验证了基于第一性原理计算方法的可靠性，为深入理解VO₂相变以及氧空位对相变的影响提供了有力的支持。5.3结果讨论实验结果与第一性原理计算结果的高度一致性，充分验证了基于第一性原理研究VO₂相变及氧空位影响的可靠性和有效性。这表明第一性原理计算能够准确地预测VO₂在不同条件下的相变行为和结构、电子性质变化，为进一步研究VO₂材料提供了有力的理论工具。尽管实验与计算结果总体相符，但仍存在一些细微差异。在相变温度的测量和计算中，实验值与计算值之间存在一定偏差，这可能是由于实验过程中样品的不均匀性、测量误差以及计算过程中对体系的简化处理等因素导致。在实际样品中，可能存在杂质、缺陷等，这些因素会影响VO₂的相变温度。在计算过程中，虽然采用了较为精确的理论方法，但仍无法完全考虑到所有的相互作用和实际因素，从而导致计算结果与实验值存在一定差异。在晶体结构和电子结构的分析中，也存在一些差异。实验中通过XRD和拉曼光谱等技术得到的晶体结构信息，与计算结果在某些细节上不完全一致。这可能是由于实验技术本身的分辨率限制，无法准确测量到一些微小的结构变化。在电子结构方面，实验测量的电子态密度与计算结果在某些能量区域的差异，可能是由于实验测量过程中的表面效应、电子弛豫等因素的影响。这些差异为进一步研究提供了方向。在后续研究中，可以进一步优化实验条件，提高样品的质量和测量的精度，以减少实验误差。在计算方面，可以采用更精确的理论方法，考虑更多的相互作用和实际因素，如电子-声子耦合、电子-电子关联等，以提高计算结果的准确性。氧空位对VO₂相变的影响研究具有重要的实际应用潜力。通过引入氧空位来调控VO₂的相变温度，能够使其更适合在室温附近工作，满足智能窗户、传感器等领域的应用需求。在智能窗户中，VO₂薄膜的相变温度接近室温时，能够根据环境温度自动调节对红外光的透过率，实现室内温度的智能调节，从而降低能源消耗。在传感器领域，通过调控相变温度，可以提高传感器对温度、气体等的灵敏度和响应速度。氧空位对VO₂电学和光学性能的影响，为开发新型的光电器件提供了可能。利用氧空位导致的VO₂电学性能变化，可以制备高性能的电阻开关、场效应晶体管等电子器件。在光学领域，氧空位对VO₂光学性能的调控，有望应用于光调制器、光探测器等光电器件的研发，提高器件的性能和功能多样性。六、结论与展望6.1研究成果总结本研究基于第一性原理，深入系统地探究了VO₂的相变过程以及氧空位对相变的影响，取得了一系列具有重要理论和实际意义的研究成果。在VO₂相变过程的研究中，通过精心构建VO₂单斜相和四方相的晶体结构模型，并合理设置平面波截断能为400eV、k点网格为4×4×4以及采用GGA-PBE泛函等关键计算参数，确保了计算结果的高精度和可靠性。计算结果清晰地揭示了VO₂在相变过程中晶体结构和电子结构的变化规律。从晶体结构方面来看，随着温度升高，VO₂从低温单斜相逐渐转变为高温四方相，晶格常数、原子坐标以及键长和键角都发生了显著变化。V-V对在相变过程中被打破，V原子的排列方式从折线型的V-V链转变为四方相中的规则排列，VO₆八面体的形状和连接方式也相应改变，晶体结构的对称性显著提高。从电子结构角度分析，低温单斜相VO₂(M)在费米能级处存在明显能隙，呈现绝缘性；而高温四方相VO₂(R)的费米能级进入导带，能隙消失，表现出金属导电性。在相变过程中，电子在不同轨道之间发生转移，从主要占据V原子的3dt₂g轨道部分转移到3deg轨道，电子-声子相互作用在这个过程中也起到了重要的促进作用。通过对这些变化的深入分析，明确了VO₂的相变是电子关联与结构变化共同作用的结果，电子关联作用的变化导致原子位置和运动状态改变，进而引起晶体结构变化，而晶体结构的变化又反过来影响电子的分布和运动状态，两者相互影响、相互促进，最终导致VO₂实现金属-绝缘体相变。在氧空位对VO₂相变影响的研究中，成功构建了不同氧空位浓度和位置的计算模型，全面深入地研究了氧空位对VO₂晶体结构、电子结构、相变温度和相变机制的影响。研究发现，氧空位的引入显著改变了VO₂的晶格参数，随着氧空位浓度的增加，VO₂单斜相的晶格常数aₘ逐渐增大，bₘ和cₘ逐渐减小，晶体体积膨胀，四方相的晶格参数也发生相应变化。氧空位还导致VO₂晶体中原子间距离和键角改变，VO₆八面体的畸变加剧，原子间的相互作用力发生改变。通过计算氧空位的形成能，确定了位于VO₆八面体顶点位置的氧空位形成能最低，最为稳定，最容易形成。在电子结构方面，氧空位的引入使VO₂的能带结构发生显著变化，能隙逐渐减小，且在能隙中出现新的杂质能级。随着氧空位浓度的增加，费米能级附近的态密度显著增加，氧空位周围的电子云密度发生变化，与氧空位相邻的V原子周围电子云密度明显增加，这些变化导致VO₂的电学和光学等物理性质发生改变。在相变温度和相变机制方面，随着氧空位浓度的增加，VO₂的相变温度逐渐降低。这是因为氧空位的引入导致晶格畸变，破坏了原有的晶体结构稳定性，降低了相变所需的能量势垒；同时，氧空位改变了VO₂的电子云分布和能带结构，促进了电子的离域运动，降低了相变过程中的电子能量变化。氧空位的存在还改变了VO₂的相变机制，使得电子结构的变化在相变过程中起到更为关键的作用，相变过程中电子的转移和能级变化更为明显，晶体结构的变化相对减弱。为了验证第一性原理计算结果的准确性和可靠性，精心设计并开展了实验研究。采用化学还原法成功制备了含不同氧空位浓度的VO₂样品，并运用XRD、