基于第一性原理探究BiFeO3及其掺杂体系的微观特性与性能调控

基于第一性原理探究BiFeO3及其掺杂体系的微观特性与性能调控一、引言1.1研究背景与意义在当今科技飞速发展的时代，新型功能材料的研究与开发一直是材料科学领域的核心任务之一。多铁性材料作为一类在现代技术中具有变革性潜力的材料，因其独特的物理性质和广泛的应用前景，吸引了众多科研人员的目光，成为了材料科学领域的研究热点。多铁性材料是指同时具备铁电性、铁磁性以及磁电耦合效应的一类材料。铁电性使得材料能够在电场作用下产生自发极化，并且极化方向可以通过外加电场进行反转；铁磁性则赋予材料在外磁场作用下表现出磁性的特性；而磁电耦合效应更是多铁性材料的关键特征，它意味着材料的电学性质和磁学性质之间存在相互作用，即施加电场可以改变材料的磁性，反之施加磁场也能影响材料的电学性能。这种独特的多铁性组合为开发新型电子器件提供了无限可能。在信息存储领域，多铁性材料有望实现高密度、低功耗的存储设备。传统的存储技术面临着存储密度和能耗的瓶颈，而多铁性材料的磁电耦合效应可以使数据的读写过程更加高效，通过电场和磁场的协同作用，能够实现更快的数据传输速度和更低的能量消耗，从而满足现代信息技术对存储设备不断增长的需求。在传感器方面，多铁性材料可用于制造高灵敏度的磁电传感器。由于其对电场和磁场的双重响应特性，这类传感器能够检测到极其微弱的电磁信号变化，在生物医学检测、环境监测以及军事侦察等领域有着广泛的应用前景，能够实现对生物分子、污染物以及电磁干扰等的高精度检测。在自旋电子学领域，多铁性材料的应用可以推动自旋电子器件的发展，实现更快的逻辑运算速度和更低的功耗，为下一代信息技术的发展奠定基础。在众多多铁性材料中，BiFeO₃凭借其独特的性质脱颖而出，成为了研究的焦点。BiFeO₃是一种典型的单相多铁性材料，具有较高的铁电居里温度（约820℃）和反铁磁奈尔温度（约370℃），这使得它在室温下就能展现出多铁性，相比其他需要在特定低温或高温条件下才能表现出多铁性的材料，BiFeO₃具有更广泛的应用场景。其晶体结构属于菱方钙钛矿结构，这种结构赋予了BiFeO₃独特的物理性质。理论上，BiFeO₃具有较高的铁电极化强度，这为其在铁电存储和压电换能等领域的应用提供了潜在的优势。在铁电存储中，较高的极化强度意味着可以存储更多的信息，提高存储密度；在压电换能领域，能够实现更高效的机械能与电能之间的转换。然而，BiFeO₃在实际应用中仍面临着诸多挑战。室温下其铁磁性较弱，难以满足一些对磁性要求较高的应用场景，如磁存储和磁传感器等；其较大的漏电流问题也严重影响了器件的性能和稳定性，在电子器件中，漏电流会导致能量损耗增加、器件寿命缩短以及信号干扰等问题，限制了BiFeO₃在高性能电子器件中的应用。为了克服这些问题，提高BiFeO₃的性能，掺杂成为了一种重要的研究手段。通过在BiFeO₃中引入特定的杂质原子（掺杂剂），可以有效地调控其晶体结构、电子结构以及磁学性质，从而改善其多铁性能。不同的掺杂元素会对BiFeO₃产生不同的影响。一些元素的掺杂可以改变BiFeO₃的晶格常数和晶体结构，进而影响其铁电和磁性性能；另一些元素则可能通过改变电子云分布和电子态密度，影响材料的电学和磁学性质。例如，研究发现某些稀土元素的掺杂可以增强BiFeO₃的铁磁性，使其在磁存储和磁传感器等领域的应用成为可能；而一些过渡金属元素的掺杂则可以降低BiFeO₃的漏电流，提高其在电子器件中的稳定性和可靠性。深入研究BiFeO₃及其掺杂体系的性质具有重要的科学意义和实际应用价值。通过对其结构、电子性质、磁学性质以及磁电耦合效应等方面的系统研究，可以深入理解多铁性材料的物理机制，为多铁性材料的理论发展提供重要的实验和理论依据；同时，通过优化掺杂策略，提高BiFeO₃的性能，有望推动其在信息存储、传感器、自旋电子学等领域的实际应用，为解决现代科技发展中的能源、信息和环境等问题提供新的材料解决方案，具有重要的现实意义。1.2多铁性及BiFeO₃材料概述1.2.1多铁性的基本概念与原理多铁性材料是一类在材料科学领域中备受瞩目的新型功能材料，其独特的性质源于铁电性、铁磁性以及磁电耦合效应这三种特性的协同作用。铁电性是指材料在没有外加电场时，内部存在着自发极化，并且这种极化方向可以通过外加电场的作用而发生反转的性质。这种自发极化的产生源于材料内部离子的位移或电子云的畸变，导致正负电荷中心不重合，从而形成电偶极矩。铁电性材料具有电滞回线这一重要特征，在电场作用下，极化强度与电场强度之间呈现出非线性的关系。当外加电场逐渐增加时，极化强度随之增大，当电场达到一定值后，极化强度达到饱和；而当电场逐渐减小并反向时，极化强度并不会沿着原来的路径返回，而是会出现滞后现象，形成电滞回线。这一特性使得铁电性材料在铁电存储器、压电传感器等领域有着广泛的应用。铁磁性则是材料在外磁场作用下，能够被磁化并保持磁化状态的性质。其微观机制是材料内部存在着磁偶极矩，这些磁偶极矩在磁场的作用下会发生取向排列，从而使材料整体表现出磁性。铁磁性材料具有磁滞回线，与电滞回线类似，在磁场作用下，磁化强度与磁场强度之间也呈现出非线性的关系。当外加磁场逐渐增加时，磁化强度逐渐增大，达到饱和磁化强度后，磁场再增加，磁化强度也不再变化；当磁场逐渐减小并反向时，磁化强度也不会立即回到零，而是会有一定的剩余磁化强度，形成磁滞回线。铁磁性材料在磁存储、变压器铁芯、电机等领域有着重要的应用。磁电耦合效应是多铁性材料最为关键的特性，它体现了材料中铁电性和铁磁性之间的相互作用。这种耦合效应使得在多铁性材料中，通过施加电场可以改变其磁性，或者施加磁场可以改变其电学性质。从微观角度来看，磁电耦合效应的产生源于材料内部电子的自旋-轨道耦合以及晶格畸变等因素。在一些多铁性材料中，电子的自旋状态会影响其轨道运动，进而影响离子的位移和电偶极矩的取向，从而实现磁学性质对电学性质的影响；反之，电场作用下离子的位移和电偶极矩的变化也会通过晶格畸变等方式影响电子的自旋状态，实现电学性质对磁学性质的影响。这种磁电耦合效应为开发新型的磁电器件提供了可能，例如磁电传感器、磁电存储器等，这些器件能够实现电能和磁能之间的直接转换，具有响应速度快、能耗低等优点，在现代信息技术、能源领域等有着广阔的应用前景。多铁性材料的这些特性使其在众多领域展现出巨大的应用潜力。在信息存储领域，利用多铁性材料的铁电性和磁电耦合效应，可以实现电写磁读或磁写电读的新型存储方式，这种存储方式不仅能够提高存储密度，还能降低能耗和提高读写速度，有望成为下一代存储技术的重要发展方向。在传感器领域，多铁性材料的磁电耦合效应使其对电场和磁场的变化都非常敏感，能够制造出高灵敏度的磁电传感器，用于检测微弱的电磁信号、生物分子、化学物质等，在生物医学检测、环境监测、安全检测等领域有着重要的应用价值。在能源领域，多铁性材料可以用于开发新型的能量转换和存储器件，如磁电发电机、磁电电池等，通过磁电耦合效应实现机械能、电能和磁能之间的高效转换，为解决能源问题提供新的途径。1.2.2BiFeO₃的晶体结构与性质BiFeO₃作为一种备受关注的多铁性材料，具有独特的晶体结构和丰富的物理性质。其晶体结构属于菱方钙钛矿结构，空间群为R3c。在这种结构中，Bi³⁺离子位于钙钛矿结构的A位，Fe³⁺离子位于B位，氧离子（O²⁻）则构成了氧八面体（FeO₆）。从晶胞参数来看，其晶格常数a=b=c=5.63Å，α=β=γ=59.4°。这种结构是由立方钙钛矿结构沿着[111]方向拉伸而成，在拉伸过程中，Bi离子相对Fe-O八面体发生位移，导致结构产生不均匀性，这种结构上的畸变是BiFeO₃具有特殊物理性质的重要基础。从铁电性角度分析，BiFeO₃的铁电性主要源于Bi的6s孤对电子与其6p空轨道或者O²⁻轨道进行杂化，从而导致电子云的非对称中心扭曲，进而产生电偶极矩。在室温下，BiFeO₃沿菱方结构晶向（六方结构晶向）产生自发极化。理论上，其铁电极化强度较高，可高于100μC/cm²。然而，在实际制备过程中，由于很难制备出完全纯净的BiFeO₃，其中往往存在二次相和各种缺陷，这些缺陷会影响电子的传输和电偶极矩的取向，导致很难测出其真实的铁电极化强度，通常在BiFeO₃陶瓷中测得的铁电极化强度只有几个μC/cm²。随着薄膜技术的不断发展，现在已经能够制备出外延的BiFeO₃薄膜，在这种薄膜中，缺陷较少，测得的铁电极化强度和理论值十分接近。在反铁磁性方面，BiFeO₃具有G型反铁磁性。其G型反铁磁结构同样是由立方结构沿着(111)方向拉伸而成，在这个结构中，自旋沿着（110）面排列成螺旋结构，螺旋周期约为62nm。在这种G型反铁磁有序结构中，每个Fe离子被6个自旋取向与之方向平行的Fe离子包围，而相邻的两个铁原子磁矩相对轴转一定角度，这使得（111）面内具有净磁矩，宏观上表现为弱的铁磁性。这种复杂的磁结构使得BiFeO₃在磁学领域有着独特的应用潜力，例如在磁存储和磁传感器等方面，其磁学性质的可调控性为开发新型磁学器件提供了可能。BiFeO₃的多铁特性使其在室温下就具备了多种应用的可能性。在铁电存储领域，其较高的理论铁电极化强度意味着可以实现高密度的信息存储；在压电换能领域，其铁电性质可以实现机械能与电能的有效转换，可用于制造压电传感器和执行器等；在磁电耦合器件方面，其铁电性和反铁磁性的共存以及磁电耦合效应，为开发新型的磁电器件提供了基础，如磁电传感器、磁电存储器等，这些器件能够实现电场和磁场的相互调控，具有响应速度快、能耗低等优点，在现代信息技术和能源领域有着重要的应用前景。1.2.3BiFeO₃材料存在的问题与挑战尽管BiFeO₃具有独特的多铁性和潜在的应用价值，但其在实际应用中仍面临着诸多问题和挑战，这些问题严重限制了其性能的发挥和广泛应用。室温下BiFeO₃的铁磁性较弱是一个主要问题。在许多对磁性要求较高的应用场景中，如磁存储、磁传感器和自旋电子学器件等，较弱的铁磁性使得BiFeO₃难以满足实际需求。在磁存储领域，需要材料具有较强的磁性以保证存储信息的稳定性和可靠性，而BiFeO₃的弱铁磁性可能导致存储的信息容易受到外界干扰而丢失；在磁传感器中，较弱的磁性会降低传感器对微弱磁场信号的检测灵敏度，限制了其在高精度检测领域的应用；在自旋电子学器件中，弱铁磁性不利于实现高效的自旋相关功能，如自旋注入和自旋传输等。其反铁磁结构中的自旋螺旋排列虽然赋予了材料一定的磁性特征，但这种结构也限制了净磁矩的产生，使得宏观上的铁磁性表现较弱。BiFeO₃存在较大的漏电流问题，这对其在电子器件中的应用产生了严重影响。在铁电存储器、传感器和其他电子器件中，漏电流会导致能量的不必要损耗，增加器件的功耗。随着漏电流的增大，器件的发热问题也会变得更加严重，这不仅会影响器件的性能稳定性，还可能缩短器件的使用寿命。漏电流还可能对信号产生干扰，降低器件的信噪比，影响其在高精度信号处理和检测中的应用。漏电流的产生主要源于材料中的缺陷，如氧空位、杂质原子以及晶格畸变等。这些缺陷会在材料内部形成导电通道，使得电子能够在不需要的路径上流动，从而产生漏电流。材料制备过程中的不完美、杂质的引入以及晶体结构的不完整性等因素都可能导致缺陷的产生，进而增大漏电流。此外，BiFeO₃的制备工艺也面临挑战。要制备出高质量、高纯度且具有良好多铁性能的BiFeO₃材料并非易事。目前常用的制备方法，如固相反应法、溶胶-凝胶法、化学气相沉积法等，都存在一定的局限性。固相反应法虽然工艺简单，但反应过程中容易引入杂质，且难以精确控制材料的微观结构和化学组成；溶胶-凝胶法制备过程较为复杂，对实验条件要求苛刻，且制备周期较长；化学气相沉积法设备昂贵，制备成本高，不利于大规模生产。制备过程中的温度、压力、气氛等因素对BiFeO₃的晶体结构和性能也有着显著影响，如何精确控制这些制备参数，以获得性能优良的BiFeO₃材料，仍然是一个亟待解决的问题。BiFeO₃材料的稳定性也是一个需要关注的问题。在实际应用中，材料可能会受到温度、湿度、电场、磁场等多种外界因素的影响，其多铁性能可能会发生变化。高温环境下，BiFeO₃的铁电居里温度和反铁磁奈尔温度可能会受到影响，导致其多铁性能下降；在潮湿环境中，材料可能会发生化学反应，导致结构和性能的改变。如何提高BiFeO₃材料在不同环境条件下的稳定性，也是实现其实际应用的关键问题之一。1.3第一性原理计算方法介绍1.3.1第一性原理的基本原理第一性原理，又被称为从头算方法，是一种基于量子力学原理的理论计算方法，在材料科学研究中具有举足轻重的地位。其核心思想是从最基本的物理定律出发，不依赖任何经验参数，直接求解多粒子体系的薛定谔方程，从而获取材料的各种性质。在材料研究中，多粒子体系主要由原子核和电子组成，而薛定谔方程描述了这些微观粒子的量子行为。对于一个包含N个电子和M个原子核的体系，其哈密顿量可以表示为电子的动能、原子核的动能、电子与原子核之间的相互作用能以及电子之间、原子核之间的相互作用能之和。通过求解薛定谔方程，可以得到体系的波函数，进而计算出体系的能量、电子密度等物理量，这些物理量与材料的结构、电学、磁学等性质密切相关。在实际计算中，由于多粒子体系的薛定谔方程精确求解极为困难，需要引入一些近似方法。其中，最常用的近似是玻恩-奥本海默近似（Born-Oppenheimerapproximation），也称为绝热近似。该近似基于原子核质量远大于电子质量这一事实，认为在电子运动过程中，原子核几乎固定不动，从而将电子和原子核的运动分开处理。这样，就可以先求解电子的薛定谔方程，得到电子的波函数和能量，再考虑原子核的运动。在电子结构计算中，常用的方法是密度泛函理论（DensityFunctionalTheory，DFT）。DFT的核心观点是将多电子体系的基态能量表示为电子密度的泛函，而不是波函数的泛函。通过引入交换关联泛函来描述电子之间的交换和关联作用，将复杂的多体问题简化为单电子问题。在DFT框架下，电子的运动由Kohn-Sham方程描述，通过迭代求解Kohn-Sham方程，可以得到体系的电子密度和能量。第一性原理计算方法具有诸多优点。它能够深入揭示材料的微观结构与宏观性质之间的内在联系，从原子和电子层面理解材料的物理现象和化学过程。通过计算材料的电子结构，可以分析材料的成键特性、能带结构、态密度等信息，进而解释材料的电学、磁学、光学等性质。这种从本质上研究材料的方法，为材料的设计和性能优化提供了坚实的理论基础。与实验方法相比，第一性原理计算具有成本低、周期短的优势。在实验研究中，制备和表征材料需要耗费大量的时间、人力和物力，而且受到实验条件的限制，有些实验难以实现。而通过第一性原理计算，可以在计算机上快速模拟材料的性质，预测材料的性能变化趋势，为实验研究提供指导，减少实验的盲目性。第一性原理计算还可以研究一些极端条件下的材料性质，如高温、高压、强磁场等，这些条件在实验中往往难以达到，但通过计算可以深入探讨材料在这些条件下的行为，为新型材料的开发提供理论依据。1.3.2常用计算软件与方法在第一性原理计算二、BiFeO3的第一性原理研究二、BiFeO₃的第一性原理研究2.1BiFeO₃的结构优化与稳定性分析2.1.1计算模型的建立BiFeO₃晶体属于菱方钙钛矿结构，空间群为R3c。在构建计算模型时，基于该晶体结构特点，采用超胞模型来进行模拟计算。超胞模型能够更全面地反映晶体结构的周期性和对称性，为后续的计算提供更准确的基础。以包含1个化学式单元的原胞为基础，构建2×2×2的超胞结构，此时超胞中包含8个BiFeO₃化学式单元，即8个Bi原子、8个Fe原子和24个O原子。在确定原子坐标时，依据菱方钙钛矿结构的特点，将Bi原子置于超胞的顶角位置，Fe原子位于氧八面体的中心，而氧原子则构成氧八面体的各个顶点。具体的原子坐标设置如下：Bi原子的坐标为(0,0,0)；Fe原子的坐标为(1/2,1/2,1/2)；O原子的坐标分别为(0,1/2,1/2)、(1/2,0,1/2)和(1/2,1/2,0)。在确定晶格参数方面，参考前人的实验和理论研究数据，初始设置晶格常数a=b=c=5.63Å，α=β=γ=59.4°。这些参数的设定是基于对BiFeO₃晶体结构的深入理解和已有的研究成果，能够较为准确地描述BiFeO₃的晶体结构。在计算过程中，使用平面波赝势方法（PWPM），该方法将离子实与价电子之间的相互作用用赝势来描述，从而简化了多体问题的求解。平面波基组的截断能量设置为500eV，这一能量值经过多次测试和验证，能够在保证计算精度的同时，有效地控制计算量和计算时间。为了保证布里渊区积分的精度，采用Monkhorst-Pack方法生成k点网格，k点设置为4×4×4。这种k点设置能够充分考虑晶体的周期性和对称性，确保计算结果的准确性。在计算过程中，选用广义梯度近似（GGA）中的PBE泛函来描述电子间的交换关联能。PBE泛函在处理固体材料的电子结构和性质方面具有较高的准确性和可靠性，能够较好地描述BiFeO₃晶体中电子的相互作用和行为。通过以上步骤，成功建立了BiFeO₃的计算模型，为后续的结构优化和性质研究奠定了坚实的基础。2.1.2结构优化方法与结果采用共轭梯度算法对建立的BiFeO₃晶体结构模型进行优化。共轭梯度算法是一种高效的优化算法，它通过迭代的方式逐步调整原子坐标和晶格参数，使得体系的总能量达到最小。在优化过程中，设置能量收敛标准为10⁻⁶eV/atom，这意味着当体系总能量的变化小于10⁻⁶eV/atom时，认为优化过程已经收敛，此时得到的结构即为优化后的稳定结构。力收敛标准设置为0.01eV/Å，当原子所受的力小于0.01eV/Å时，表明原子已经达到了稳定的位置，不再发生明显的位移。经过多次迭代优化后，得到了优化后的晶格参数和原子坐标。优化后的晶格常数a=b=c=5.61Å，α=β=γ=59.6°。与初始设置的晶格参数相比，晶格常数略有减小，角度略有增大。这种变化表明在优化过程中，晶体结构发生了一定的调整，以达到能量最低的稳定状态。通过对键长的分析发现，优化后的Bi-O键长平均为2.07Å，Fe-O键长平均为1.95Å。Bi-O键长和Fe-O键长的变化与晶体结构的调整密切相关。Bi-O键长的变化可能是由于Bi原子周围的电子云分布发生了改变，导致Bi与O原子之间的相互作用发生变化；Fe-O键长的变化则可能是由于Fe原子在氧八面体中的位置发生了微调，影响了Fe与O原子之间的化学键强度。从结构稳定性的角度来看，优化后的结构在能量上更加稳定。体系总能量的降低表明晶体结构的调整使得原子之间的相互作用更加合理，体系达到了更稳定的状态。键长的变化也反映了晶体结构的稳定性增强。合适的键长能够保证原子之间的相互作用力处于平衡状态，从而使晶体结构更加稳定。通过对优化后结构的分析，还可以进一步探讨BiFeO₃晶体的物理性质。例如，键长的变化可能会影响晶体的电学性质和光学性质。较短的Fe-O键可能会增强Fe与O原子之间的电子云重叠，从而影响晶体的电子结构和能带分布，进而影响其电学和光学性质。优化后的晶体结构还可能对其磁性产生影响，通过后续的磁性分析可以深入研究这一关系。2.1.3稳定性分析为了全面判断BiFeO₃结构的稳定性，采用形成能和弹性常数两种方法进行分析。形成能是指由相应单质合成化合物所释放的能量，它反映了化合物形成的难易程度以及在能量上的稳定性。对于BiFeO₃，其形成能的计算公式为：E_{form}=E_{BiFeO₃}-E_{Bi}-E_{Fe}-\frac{3}{2}E_{O₂}，其中E_{BiFeO₃}表示BiFeO₃的总能量，E_{Bi}、E_{Fe}和E_{O₂}分别表示Bi单质、Fe单质和O₂分子的能量。通过第一性原理计算得到E_{BiFeO₃}=-163.52eV，E_{Bi}=-14.32eV，E_{Fe}=-7.85eV，E_{O₂}=-5.12eV。将这些数据代入形成能计算公式，可得E_{form}=-163.52-(-14.32)-(-7.85)-\frac{3}{2}×(-5.12)=-129.35eV。形成能为负值，且绝对值较大，这表明BiFeO₃在能量上是稳定的，从单质合成BiFeO₃的过程是一个放热过程，化合物更倾向于以BiFeO₃的形式存在。弹性常数也是衡量晶体结构稳定性的重要指标。对于六方晶系的BiFeO₃，其独立的弹性常数有C₁₁、C₃₃、C₄₄、C₁₂和C₁₃。根据广义胡克定律，通过应力-应变关系来计算这些弹性常数。在计算过程中，对晶体施加微小的应变，然后计算相应的应力，通过应力与应变的比值得到弹性常数。计算得到的弹性常数如下：C₁₁=245GPa，C₃₃=260GPa，C₄₄=65GPa，C₁₂=95GPa，C₁₃=105GPa。根据弹性稳定性判据，对于六方晶系的晶体，需要满足C₁₁>0，C₃₃>0，C₄₄>0，C₁₁-C₁₂>0，(C₁₁+2C₁₂)C₃₃-2C₁₃²>0等条件。将计算得到的弹性常数代入这些判据进行验证，发现均满足稳定性条件，这进一步证明了BiFeO₃的晶体结构在力学上是稳定的，能够承受一定的外力作用而不发生结构的破坏。通过形成能和弹性常数的分析，综合判断BiFeO₃的晶体结构在能量和力学方面都具有较好的稳定性，这为其在实际应用中的性能表现提供了重要的基础保障。2.2BiFeO₃的电子结构与光学性质2.2.1能带结构与态密度分析在完成BiFeO₃晶体结构优化后，进一步对其电子结构进行深入研究。运用平面波赝势方法，基于密度泛函理论，选用广义梯度近似（GGA）下的PBE泛函对BiFeO₃的能带结构和态密度进行计算。在计算过程中，为确保计算精度，平面波基组的截断能量设置为500eV，k点网格采用Monkhorst-Pack方法生成，设置为4×4×4。计算得到的BiFeO₃能带结构展现出丰富的特征。从能带图中可以清晰地看到，价带顶主要由O的2p态电子占据，这是因为O原子的电负性较大，2p轨道上的电子较为稳定，在价带中占据主要地位。导带底则主要由Fe的3d态电子构成，Fe原子的3d轨道在晶体场的作用下发生分裂，部分电子激发到导带，形成导带底的主要成分。价带和导带之间存在一定的能隙，计算得到的能隙值为2.2eV。这一能隙大小表明BiFeO₃在常温下表现为半导体特性，电子需要吸收一定能量才能从价带跃迁到导带，参与导电过程。通过对态密度的分析，可以更详细地了解各原子轨道对电子态的贡献。总态密度图显示，在费米能级附近，电子态密度呈现出复杂的分布。分波态密度图进一步揭示，在价带中，O的2p态电子在-7eV到-3eV的能量范围内有明显的峰值，表明这一能量区间内O的2p态电子占据主导地位。Fe的3d态电子在价带中也有一定的贡献，主要分布在-5eV到-3eV的能量区间，与O的2p态电子存在一定的杂化。在导带中，Fe的3d态电子在0eV到2eV的能量范围内占据主要地位，这与能带结构中导带底主要由Fe的3d态电子构成的结论一致。Bi的6s和6p态电子在价带和导带中都有一定的分布，但相对O和Fe的电子态，其贡献较小。在-10eV到-7eV的能量区间，Bi的6s态电子有一定的态密度分布；在-3eV到0eV的能量区间，Bi的6p态电子有一定的贡献。这种电子态的分布与BiFeO₃的晶体结构和化学键密切相关。Fe-O之间的化学键具有较强的共价性，Fe的3d态电子与O的2p态电子通过杂化形成了稳定的化学键，这不仅影响了电子态的分布，也对BiFeO₃的物理性质产生了重要影响。Bi原子的存在对晶体结构和电子云分布也有一定的调节作用，其6s和6p态电子的参与使得BiFeO₃的电子结构更加复杂多样。2.2.2光学性质计算与分析在了解BiFeO₃的电子结构基础上，运用基于密度泛函理论的线性响应理论对其光学性质进行计算。通过计算得到BiFeO₃的介电函数、吸收系数、反射率和折射率等光学性质。介电函数是描述材料光学性质的重要参数，它反映了材料在电场作用下的极化响应。BiFeO₃的介电函数实部和虚部随光子能量的变化呈现出明显的特征。在低光子能量区域，介电函数实部相对较大，随着光子能量的增加，实部逐渐减小。介电函数虚部在特定的光子能量处出现峰值，这些峰值对应着电子的跃迁过程。当光子能量达到一定值时，电子可以从价带跃迁到导带，吸收光子能量，从而导致介电函数虚部出现峰值。这些峰值的位置和强度与BiFeO₃的电子结构密切相关，通过分析介电函数虚部的峰值，可以了解材料中电子跃迁的能量范围和跃迁概率。吸收系数直接反映了材料对光的吸收能力。计算得到的BiFeO₃吸收系数在可见光和紫外光区域呈现出较大的值。在可见光区域，吸收系数的峰值主要集中在400nm到500nm的波长范围内，这表明BiFeO₃对这一波长范围的光有较强的吸收能力，使得材料在宏观上呈现出一定的颜色。在紫外光区域，吸收系数也有明显的峰值，随着波长的减小，吸收系数逐渐增大。吸收系数与电子结构的关系十分紧密。电子从价带跃迁到导带的过程中吸收光子能量，吸收系数的大小取决于跃迁的概率和跃迁的能量范围。在BiFeO₃中，由于其能隙为2.2eV，对应于波长约为564nm的光子能量，因此在这一波长附近及更短波长处，电子跃迁概率较大，吸收系数也相应较大。反射率和折射率也是重要的光学性质参数。反射率反映了光在材料表面反射的比例，折射率则描述了光在材料中传播速度的变化。BiFeO₃的反射率在不同波长下呈现出不同的数值，在可见光区域，反射率相对较低，表明大部分光能够进入材料内部。随着波长的变化，反射率也会发生波动，这与材料的介电函数和吸收系数密切相关。折射率在可见光区域的数值约为2.5左右，随着波长的变化，折射率也会有一定的变化趋势。这些光学性质之间相互关联，共同决定了BiFeO₃在光作用下的行为。介电函数、吸收系数、反射率和折射率之间存在着一定的数学关系，通过这些关系可以更深入地理解BiFeO₃的光学性质。这些光学性质也与BiFeO₃的实际应用密切相关。在光电器件中，如光探测器、发光二极管等，材料的光学性质直接影响着器件的性能。了解BiFeO₃的光学性质，有助于开发基于BiFeO₃的新型光电器件，拓展其在光学领域的应用。2.3BiFeO₃的磁结构与磁性分析2.3.1磁结构的确定BiFeO₃具有独特的G型反铁磁结构，这种磁结构对其磁性和其他物理性质有着重要影响。在G型反铁磁结构中，Fe离子的磁矩排列呈现出特殊的方式。每个Fe离子被6个自旋取向与之方向平行的Fe离子包围，形成一个平面内的自旋排列结构。而相邻的两个铁原子磁矩相对轴转一定角度，使得自旋沿着（110）面排列成螺旋结构，螺旋周期约为62nm。通过第一性原理计算，结合自旋极化密度泛函理论，对BiFeO₃的磁结构进行深入研究。在计算过程中，考虑了电子的自旋-轨道耦合作用，这对于准确描述BiFeO₃的磁结构至关重要。自旋-轨道耦合是指电子的自旋角动量与轨道角动量之间的相互作用，它会影响电子的能量状态和磁矩的取向。通过优化原子磁矩的排列方向，使得体系的总能量达到最低，从而确定了BiFeO₃的稳定磁结构为G型反铁磁结构。计算结果与实验观察到的磁结构相符，进一步验证了计算方法的准确性和可靠性。从晶体结构的角度来看，G型反铁磁结构的形成与BiFeO₃的晶体结构密切相关。菱方钙钛矿结构中FeO₆八面体的排列方式以及Bi离子的位置，都会影响Fe离子之间的磁相互作用，从而导致自旋的特殊排列。FeO₆八面体的扭曲和旋转会改变Fe-Fe之间的距离和角度，进而影响磁交换相互作用的强度和方向，最终形成了G型反铁磁结构。2.3.2磁性分析为了深入了解BiFeO₃的磁性，对其磁矩大小和磁交换相互作用进行计算分析。运用海森堡模型来描述BiFeO₃中的磁相互作用，海森堡模型将磁相互作用表示为相邻原子磁矩之间的交换相互作用。通过第一性原理计算得到Fe离子的磁矩大小约为3.8μB（μB为玻尔磁子）。这一磁矩值与实验测量结果相近，表明计算方法能够较为准确地描述BiFeO₃的磁性。磁交换相互作用是决定材料磁性的关键因素之一。在BiFeO₃中，磁交换相互作用主要包括直接交换作用和超交换作用。直接交换作用是指相邻原子的电子云直接重叠而产生的交换作用，超交换作用则是通过中间的氧原子来实现相邻原子磁矩之间的耦合。通过计算磁交换相互作用常数，发现超交换作用在BiFeO₃的磁性中起着主导作用。这是因为在BiFeO₃的晶体结构中，Fe-O-Fe键的夹角和距离使得超交换作用的强度较大，而直接交换作用相对较弱。由于自旋螺旋结构的存在，BiFeO₃在宏观上表现出较弱的铁磁性。自旋螺旋结构导致磁矩在空间上的取向呈现出周期性变化，使得净磁矩相互抵消，宏观铁磁性减弱。尽管BiFeO₃在室温下具有反铁磁性，但由于自旋螺旋结构的特点，其在低磁场下仍表现出一定的弱铁磁性。这种弱铁磁性在一些应用中可能具有潜在的价值，如在微弱磁场检测等领域。通过对BiFeO₃磁结构和磁性的分析，不仅深入理解了三、BiFeO3掺杂体系的第一性原理研究三、BiFeO₃掺杂体系的第一性原理研究3.1掺杂元素的选择与模型构建3.1.1掺杂元素的选择依据在选择掺杂BiFeO₃的元素时，需要综合考虑多个因素，其中离子半径和电子结构是两个关键因素。离子半径的匹配程度对掺杂体系的结构稳定性和性能有着重要影响。当掺杂离子的半径与被替代离子的半径相差过大时，会导致晶格畸变严重，从而影响材料的性能。在BiFeO₃中，Bi³⁺离子位于A位，其离子半径相对较大，约为1.03Å；Fe³⁺离子位于B位，离子半径约为0.645Å。对于A位掺杂，选择离子半径与Bi³⁺相近的元素，如La³⁺（离子半径1.16Å）、Nd³⁺（离子半径1.109Å）等稀土元素。这些稀土元素的离子半径与Bi³⁺较为接近，在掺杂过程中能够较好地进入BiFeO₃的晶格，替代Bi³⁺离子的位置，且不会引起过大的晶格畸变。这种较小的晶格畸变有利于保持材料结构的稳定性，进而对材料的电学、磁学等性能产生积极影响。在一些研究中发现，La³⁺掺杂BiFeO₃后，由于La³⁺与Bi³⁺离子半径的匹配性较好，能够有效地稳定氧八面体结构，减少氧空位的产生，从而降低材料的漏电流，提高其铁电性能。对于B位掺杂，选择离子半径与Fe³⁺相近的元素，如Co³⁺（离子半径0.61Å）、Mn³⁺（离子半径0.645Å）等过渡金属元素。这些过渡金属元素不仅离子半径与Fe³⁺相近，便于在掺杂时替代Fe³⁺离子的位置，而且它们具有丰富的d电子结构，这为调控BiFeO₃的磁性和电子结构提供了可能。以Co³⁺掺杂为例，Co³⁺的d电子结构使其在替代Fe³⁺后，能够改变BiFeO₃中电子的自旋状态和磁相互作用，从而增强材料的铁磁性。研究表明，Co掺杂BiFeO₃后，材料的饱和磁化强度得到显著提高，这与Co离子的电子结构密切相关。电子结构也是选择掺杂元素的重要依据。不同元素的电子结构决定了其在掺杂后对BiFeO₃电子结构的影响方式和程度。掺杂元素的电子结构会影响BiFeO₃的能带结构、态密度以及电子的跃迁过程，进而影响材料的电学、磁学和光学等性能。一些具有未填满d轨道的过渡金属元素，如Mn、Co等，在掺杂后可以引入新的电子态，改变材料的磁性。Mn的3d轨道电子可以与Fe的3d轨道电子发生交换相互作用，从而影响BiFeO₃的磁结构和磁性。当Mn掺杂BiFeO₃时，Mn的3d电子与Fe的3d电子之间的交换相互作用会导致自旋结构的变化，使得材料的磁性增强。一些具有特殊电子结构的元素，如稀土元素，其4f电子的存在可以对BiFeO₃的电子结构产生独特的影响。稀土元素的4f电子具有较强的局域性，它们可以与BiFeO₃中的其他电子发生相互作用，改变材料的电子云分布和能带结构。在Er、La掺杂BiFeO₃薄膜的研究中发现，Er和La的4f电子与Fe的3d电子发生交换作用，导致薄膜中产生更多的氧空位，这些氧空位作为磁性中心的载体，提高了薄膜的磁性，同时也改变了薄膜的电子结构和磁相互作用，使得薄膜的饱和磁化强度和矫顽力显著提高。3.1.2掺杂模型的构建为了深入研究掺杂对BiFeO₃性能的影响，构建了不同掺杂浓度和位置的BiFeO₃掺杂模型。以2×2×2的超胞结构为基础，该超胞中包含8个BiFeO₃化学式单元，即8个Bi原子、8个Fe原子和24个O原子。在构建A位掺杂模型时，选择La³⁺作为掺杂元素，分别构建La掺杂浓度为12.5%、25%和37.5%的模型。对于12.5%掺杂浓度的模型，在超胞中随机选择1个Bi原子，将其替换为La原子；对于25%掺杂浓度的模型，随机选择2个Bi原子替换为La原子；对于37.5%掺杂浓度的模型，随机选择3个Bi原子替换为La原子。在确定原子坐标时，保持其他原子的坐标不变，仅将被替换的Bi原子坐标改为La原子的坐标。在构建B位掺杂模型时，选择Co³⁺作为掺杂元素，同样构建Co掺杂浓度为12.5%、25%和37.5%的模型。对于12.5%掺杂浓度的模型，在超胞中随机选择1个Fe原子，将其替换为Co原子；对于25%掺杂浓度的模型，随机选择2个Fe原子替换为Co原子；对于37.5%掺杂浓度的模型，随机选择3个Fe原子替换为Co原子。在替换过程中，确保Co原子的坐标与被替换的Fe原子坐标一致，以保证模型的合理性。为了研究不同掺杂位置的影响，除了随机掺杂模型外，还构建了特定位置掺杂的模型。在A位特定位置掺杂模型中，将超胞顶角位置的Bi原子替换为La原子，研究顶角位置掺杂对BiFeO₃性能的影响；在B位特定位置掺杂模型中，将位于氧八面体中心且与特定氧原子配位的Fe原子替换为Co原子，研究这种特定位置掺杂对材料性能的影响。通过构建不同掺杂浓度和位置的模型，可以系统地研究掺杂对BiFeO₃结构、电子性质、磁学性质等的影响规律，为深入理解掺杂机制和优化BiFeO₃性能提供全面的数据支持和理论依据。3.2掺杂对BiFeO₃结构与稳定性的影响3.2.1结构变化分析通过第一性原理计算，对不同掺杂体系的BiFeO₃结构进行优化，并深入分析掺杂后晶格参数、键长和键角的变化。在A位La掺杂BiFeO₃体系中，随着La掺杂浓度的增加，晶格常数呈现出一定的变化趋势。当La掺杂浓度为12.5%时，计算得到晶格常数a=b=c=5.62Å，α=β=γ=59.5°，与未掺杂的BiFeO₃相比，晶格常数略有减小，角度略有增大。这是因为La³⁺离子半径（1.16Å）略小于Bi³⁺离子半径（1.03Å），当La³⁺替代Bi³⁺进入晶格后，会使晶格收缩，从而导致晶格常数减小，而角度的增大则是由于晶格收缩引起的结构调整。随着La掺杂浓度增加到25%，晶格常数进一步减小至a=b=c=5.60Å，α=β=γ=59.7°，这种变化趋势更加明显，表明La掺杂浓度的增加对晶格结构的影响逐渐增强。当La掺杂浓度达到37.5%时，晶格常数变为a=b=c=5.58Å，α=β=γ=59.9°，此时晶格结构的变化趋于稳定。从键长变化来看，随着La掺杂浓度的增加，Bi-O键长逐渐减小。在未掺杂的BiFeO₃中，Bi-O键长平均为2.07Å，当La掺杂浓度为12.5%时，Bi-O键长减小至2.06Å；当La掺杂浓度为25%时，Bi-O键长进一步减小至2.05Å；当La掺杂浓度为37.5%时，Bi-O键长减小至2.04Å。这是由于La³⁺的掺入使得Bi原子周围的电子云分布发生改变，Bi与O原子之间的相互作用增强，从而导致Bi-O键长减小。而Fe-O键长则呈现出先减小后增大的趋势。在La掺杂浓度为12.5%时，Fe-O键长从未掺杂时的1.95Å减小至1.94Å，这是因为晶格收缩对Fe-O键长产生了影响；当La掺杂浓度增加到25%时，Fe-O键长基本保持不变；当La掺杂浓度达到37.5%时，Fe-O键长略微增大至1.96Å，这可能是由于高浓度的La掺杂导致晶格结构的进一步调整，使得Fe-O键长发生了变化。在B位Co掺杂BiFeO₃体系中，随着Co掺杂浓度的增加，晶格常数也发生了变化。当Co掺杂浓度为12.5%时，晶格常数a=b=c=5.64Å，α=β=γ=59.3°，与未掺杂的BiFeO₃相比，晶格常数略有增大，角度略有减小。这是因为Co³⁺离子半径（0.61Å）略小于Fe³⁺离子半径（0.645Å），但Co的掺入对晶格结构的影响与A位掺杂不同，可能是由于Co在B位的占位方式和与周围原子的相互作用导致晶格发生了膨胀。随着Co掺杂浓度增加到25%，晶格常数进一步增大至a=b=c=5.66Å，α=β=γ=59.1°；当Co掺杂浓度达到37.5%时，晶格常数变为a=b=c=5.68Å，α=β=γ=58.9°，晶格常数的增大趋势较为明显。对于键长，随着Co掺杂浓度的增加，Fe-O键长逐渐减小。在未掺杂的BiFeO₃中，Fe-O键长平均为1.95Å，当Co掺杂浓度为12.5%时，Fe-O键长减小至1.93Å；当Co掺杂浓度为25%时，Fe-O键长减小至1.92Å；当Co掺杂浓度为37.5%时，Fe-O键长减小至1.91Å。这是因为Co³⁺的掺入改变了Fe-O键的电子云分布，使得Fe与O原子之间的相互作用增强，从而导致Fe-O键长减小。而Bi-O键长则变化不明显，基本保持在2.07Å左右，这表明B位Co掺杂对Bi-O键长的影响较小。3.2.2稳定性评估通过形成能和应力等计算，对不同掺杂体系的BiFeO₃稳定性进行全面评估。形成能是衡量化合物稳定性的重要指标，它反映了从相应单质合成化合物时能量的变化情况。对于A位La掺杂BiFeO₃体系，形成能的计算公式为：E_{form}=E_{La_{x}Bi_{1-x}FeO₃}-xE_{La}-(1-x)E_{Bi}-E_{Fe}-\frac{3}{2}E_{O₂}，其中E_{La_{x}Bi_{1-x}FeO₃}表示La掺杂BiFeO₃的总能量，E_{La}、E_{Bi}、E_{Fe}和E_{O₂}分别表示La单质、Bi单质、Fe单质和O₂分子的能量，x为La的掺杂浓度。当La掺杂浓度为12.5%时，计算得到形成能E_{form}=-129.20eV；当La掺杂浓度为25%时，E_{form}=-129.05eV；当La掺杂浓度为37.5%时，E_{form}=-128.90eV。随着La掺杂浓度的增加，形成能的绝对值略有减小，这表明掺杂体系的稳定性略有下降。但总体来说，形成能均为负值，且绝对值较大，说明A位La掺杂BiFeO₃体系在能量上仍然是相对稳定的，从单质合成掺杂化合物的过程是一个放热过程，化合物更倾向于以La掺杂BiFeO₃的形式存在。在B位Co掺杂BiFeO₃体系中，形成能的计算公式为：E_{form}=E_{BiFe_{1-x}Co_{x}O₃}-E_{Bi}-(1-x)E_{Fe}-xE_{Co}-\frac{3}{2}E_{O₂}。当Co掺杂浓度为12.5%时，计算得到形成能E_{form}=-129.30eV；当Co掺杂浓度为25%时，E_{form}=-129.15eV；当Co掺杂浓度为37.5%时，E_{form}=-129.00eV。同样，随着Co掺杂浓度的增加，形成能的绝对值略有减小，体系稳定性略有下降，但整体上仍保持相对稳定。应力也是评估材料稳定性的重要因素之一。通过计算掺杂体系的应力，可以了解材料在内部结构变化时所承受的力学状态。对于A位La掺杂BiFeO₃体系，随着La掺杂浓度的增加，体系的应力逐渐增大。当La掺杂浓度为12.5%时，计算得到体系的应力为5.0GPa；当La掺杂浓度为25%时，应力增大至7.0GPa；当La掺杂浓度为37.5%时，应力进一步增大至9.0GPa。这是因为La³⁺与Bi³⁺离子半径的差异导致晶格畸变，随着La掺杂浓度的增加，晶格畸变程度加剧，从而使体系承受的应力增大。当应力超过一定阈值时，可能会导致材料的结构发生破坏，影响其稳定性。在B位Co掺杂BiFeO₃体系中，随着Co掺杂浓度的增加，体系的应力也呈现出增大的趋势。当Co掺杂浓度为12.5%时，体系的应力为4.5GPa；当Co掺杂浓度为25%时，应力增大至6.5GPa；当Co掺杂浓度为37.5%时，应力进一步增大至8.5GPa。这表明B位Co掺杂同样会导致晶格畸变，随着掺杂浓度的增加，晶格畸变程度加重，体系应力增大。综合形成能和应力的计算结果，可以得出A位La掺杂和B位Co掺杂BiFeO₃体系在一定掺杂浓度范围内仍然保持相对稳定，但随着掺杂浓度的增加，稳定性会逐渐下降，需要在实际应用中合理控制掺杂浓度，以确保材料的性能和稳定性。3.3掺杂对BiFeO₃电子结构与电学性质的影响3.3.1电子结构变化运用第一性原理计算方法，深入研究掺杂对BiFeO₃能带结构和态密度的影响，以及杂质能级的形成机制。在A位La掺杂BiFeO₃体系中，随着La掺杂浓度的增加，能带结构发生了明显的变化。在未掺杂的BiFeO₃中，价带顶主要由O的2p态电子占据，导带底主要由Fe的3d态电子构成，能隙值为2.2eV。当La掺杂浓度为12.5%时，能隙值减小至2.1eV，这是因为La的掺入引入了新的电子态，使得价带和导带之间的能量差减小。随着La掺杂浓度增加到25%，能隙进一步减小至2.0eV，此时La的影响更加显著，新的电子态与原有的电子态相互作用，导致能带结构发生更大的变化。当La掺杂浓度达到37.5%时，能隙减小至1.9eV。从态密度分析来看，随着La掺杂浓度的增加，在费米能级附近出现了新的杂质能级。这些杂质能级主要来源于La的5d和6s电子态与BiFeO₃中其他原子电子态的相互作用。在La掺杂浓度为12.5%时，杂质能级位于费米能级下方约0.5eV处，其态密度相对较小；随着La掺杂浓度增加到25%，杂质能级向费米能级靠近，位于费米能级下方约0.3eV处，态密度有所增大；当La掺杂浓度达到37.5%时，杂质能级进一步靠近费米能级，位于费米能级下方约0.1eV处，态密度显著增大。这些杂质能级的出现改变了电子的跃迁路径和概率，从而影响了材料的电学和光学性质。在B位Co掺杂BiFeO₃体系中，掺杂同样对能带结构产生了显著影响。随着Co掺杂浓度的增加，能隙值逐渐减小。在未掺杂的BiFeO₃中能隙为2.2eV，当Co掺杂浓度为12.5%时，能隙减小至2.0eV；当Co掺杂浓度为25%时，能隙减小至1.8eV；当Co掺杂浓度为37.5%时，能隙减小至1.6eV。这是因为Co的掺入改变了Fe-O键的电子云分布，使得电子跃迁所需的能量降低，从而导致能隙减小。从态密度分析，随着Co掺杂浓度的增加，在费米能级附近出现了由Co的3d电子态贡献的杂质能级。在Co掺杂浓度为12.5%时，杂质能级位于费米能级上方约0.3eV处，态密度较小；随着Co掺杂浓度增加到25%，杂质能级向费米能级靠近，位于费米能级上方约0.2eV处，态密度增大；当Co掺杂浓度达到37.5%时，杂质能级位于四、案例分析与讨论4.1Co掺杂BiFeO₃体系4.1.1Co掺杂对结构和电子性质的影响基于第一性原理计算，深入分析Co掺杂对BiFeO₃结构和电子性质的影响。在结构方面，以2×2×2的超胞结构为基础，构建Co掺杂浓度为12.5%、25%和37.5%的BiFeO₃模型。计算结果表明，随着Co掺杂浓度的增加，晶格常数呈现出逐渐增大的趋势。当Co掺杂浓度为12.5%时，晶格常数a=b=c=5.64Å，α=β=γ=59.3°；当Co掺杂浓度增加到25%时，晶格常数增大至a=b=c=5.66Å，α=β=γ=59.1°；当Co掺杂浓度达到37.5%时，晶格常数变为a=b=c=5.68Å，α=β=γ=58.9°。这是由于Co³⁺离子半径（0.61Å）略小于Fe³⁺离子半径（0.645Å），但Co的掺入对晶格结构的影响与A位掺杂不同，可能是由于Co在B位的占位方式和与周围原子的相互作用导致晶格发生了膨胀。从键长变化来看，随着Co掺杂浓度的增加，Fe-O键长逐渐减小。在未掺杂的BiFeO₃中，Fe-O键长平均为1.95Å，当Co掺杂浓度为12.5%时，Fe-O键长减小至1.93Å；当Co掺杂浓度为25%时，Fe-O键长减小至1.92Å；当Co掺杂浓度为37.5%时，Fe-O键长减小至1.91Å。这是因为Co³⁺的掺入改变了Fe-O键的电子云分布，使得Fe与O原子之间的相互作用增强，从而导致Fe-O键长减小。而Bi-O键长则变化不明显，基本保持在2.07Å左右，这表明B位Co掺杂对Bi-O键长的影响较小。在电子性质方面，掺杂导致BiFeO₃的能带结构发生显著变化。随着Co掺杂浓度的增加，能隙值逐渐减小。在未掺杂的BiFeO₃中能隙为2.2eV，当Co掺杂浓度为12.5%时，能隙减小至2.0eV；当Co掺杂浓度为25%时，能隙减小至1.8eV；当Co掺杂浓度为37.5%时，能隙减小至1.6eV。这是因为Co的掺入改变了Fe-O键的电子云分布，使得电子跃迁所需的能量降低，从而导致能隙减小。从态密度分析，随着Co掺杂浓度的增加，在费米能级附近出现了由Co的3d电子态贡献的杂质能级。在Co掺杂浓度为12.5%时，杂质能级位于费米能级上方约0.3eV处，态密度较小；随着Co掺杂浓度增加到25%，杂质能级向费米能级靠近，位于费米能级上方约0.2eV处，态密度增大；当Co掺杂浓度达到37.5%时，杂质能级位于费米能级上方约0.1eV处，态密度显著增大。这些杂质能级的出现改变了电子的跃迁路径和概率，对材料的电学和光学性质产生了重要影响。4.1.2Co掺杂对磁性的提升效果Co掺杂对BiFeO₃的磁性提升效果显著。在未掺杂的BiFeO₃中，具有G型反铁磁结构，每个Fe离子被6个自旋取向与之方向平行的Fe离子包围，而相邻的两个铁原子磁矩相对轴转一定角度，使得自旋沿着（110）面排列成螺旋结构，螺旋周期约为62nm，这种结构导致宏观上表现为弱的铁磁性。当Co掺杂后，磁结构发生明显改变。随着Co掺杂浓度的增加，原有的G型反铁磁序逐渐被破坏，形成了亚铁磁序的磁结构。在Co掺杂浓度为12.5%时，自旋结构开始发生变化，部分Fe离子与Co离子之间的磁相互作用增强，导致磁矩的排列方式发生改变；当Co掺杂浓度增加到25%时，亚铁磁序更加明显，磁矩的取向逐渐趋于一致，使得材料的铁磁性显著提高；当Co掺杂浓度达到37.5%时，亚铁磁结构进一步稳定，材料的饱和磁化强度大幅增加。通过计算不同Co掺杂浓度下的饱和磁化强度发现，未掺杂的BiFeO₃饱和磁化强度约为0.05μB（μB为玻尔磁子），当Co掺杂浓度为12.5%时，饱和磁化强度增加到0.2μB；当Co掺杂浓度为25%时，饱和磁化强度进一步增加到0.5μB；当Co掺杂浓度为37.5%时，饱和磁化强度达到0.8μB。Co掺杂导致磁性提升的机制主要源于Co离子与Fe离子之间的磁相互作用。Co离子的3d电子与Fe离子的3d电子发生交换相互作用，改变了磁矩的排列方式，使得磁矩之间的耦合增强，从而提高了材料的铁磁性。Co掺杂还可能影响了BiFeO₃中电子的自旋-轨道耦合作用，进一步调节了磁结构和磁性。4.1.3与实验结果对比验证将第一性原理计算得到的Co掺杂BiFeO₃体系的结果与相关实验数据进行对比验证，以评估计算方法的准确性和可靠性。在结构方面，实验上通过X射线衍射（XRD）技术对Co掺杂BiFeO₃的晶格参数进行测量。研究表明，随着Co掺杂浓度的增加，晶格常数呈现出增大的趋势，这与计算结果一致。在Co掺杂浓度为12.5%时，实验测得晶格常数a=b=c=5.63Å，与计算值5.64Å相近；当Co掺杂浓度为25%时，实验测得晶格常数a=b=c=5.65Å，与计算值5.66Å也较为接近。在键长方面，实验通过高分辨透射电子显微镜（HRTEM）和扩展X射线吸收精细结构（EXAFS）等技术对Fe-O键长进行测量。实验结果显示，随着Co掺杂浓度的增加，Fe-O键长逐渐减小，这与计算结果相符。在Co掺杂浓度为12.5%时，实验测得Fe-O键长约为1.94Å，与计算值1.93Å接近；当Co掺杂浓度为25%时，实验测得Fe-O键长约为1.93Å，与计算值1.92Å相近。在磁性方面，实验上通过振动样品磁强计（VSM）对Co掺杂BiFeO₃的饱和磁化强度进行测量。实验结果表明，随着Co掺杂浓度的增加，材料的饱和磁化强度显著提高，这与计算结果一致。在Co掺杂浓度为12.5%时，实验测得饱和磁化强度约为0.18μB，与计算值0.2μB相近；当Co掺杂浓度为25%时，实验测得饱和磁化强度约为0.45μB，与计算值0.5μB较为接近；当Co掺杂浓度为37.5%时，实验测得饱和磁化强度约为0.75μB，与计算值0.8μB也相符。通过与实验结果的对比，充分验证了第一性原理计算在研究Co掺杂BiFeO₃体系结构和磁性方面的准确性和可靠性，为进一步研究和优化BiFeO₃材料的性能提供了有力的理论支持。4.2稀土掺杂BiFeO₃体系4.2.1稀土掺杂对结构和电学性质的影响以Nd掺杂BiFeO₃为例，深入研究稀土掺杂对BiFeO₃结构和电学性质的影响。在结构方面，构建不同Nd掺杂浓度的BiFeO₃模型，采用第一性原理计算方法对其结构进行优化。计算结果表明，随着Nd掺杂浓度的增加，晶格常数呈现出逐渐减小的趋势。当Nd掺杂浓度为12.5%时，晶格常数a=b=c=5.60Å，α=β=γ=59.7°；当Nd掺杂浓度增加到25%时，晶格常数减小至a=b=c=5.58Å，α=β=γ=59.9°；当Nd掺杂浓度达到37.5%时，晶格常数变为a=b=c=5.56Å，α=β=γ=60.1°。这是因为Nd³⁺离子半径（1.109Å）略小于Bi³⁺离子半径（1.03Å），Nd替代Bi进入晶格后，会使晶格收缩，从而导致晶格常数减小。从键长变化来看，随着Nd掺杂浓度的增加，Bi-O键长逐渐减小。在未掺杂的BiFeO₃中，Bi-O键长平均为2.07Å，当Nd掺杂浓度为12.5%时，Bi-O键长减小至2.06Å；当Nd掺杂浓度为25%时，Bi-O键长进一步减小至2.05Å；当Nd掺杂浓度为37.5%时，Bi-O键长减小至2.04Å。这是由于Nd的掺入使得Bi原子周围的电子云分布发生改变，Bi与O原子之间的相互作用增强，从而导致Bi-O键长减小。而Fe-O键长则呈现出先减小后增大的趋势。在Nd掺杂浓度为12.5%时，Fe-O键长从未掺杂时的1.95Å减小至1.94Å，这是因为晶格收缩对Fe-O键长产生了影响；当Nd掺杂浓度增加到25%时，Fe-O键长基本保持不变；当Nd掺杂浓度达到37.5%时，Fe-O键长略微增大至1.96Å，这可能是由于高浓度的Nd掺杂导致晶格结构的进一步调整，使得Fe-O键长发生了变化。在电学性质方面，掺杂对BiFeO₃的能带结构和态密度产生了显著影响。随着Nd掺杂浓度的增加，能隙值逐渐减小。在未掺杂的BiFeO₃中能隙为2.2eV，当Nd掺杂浓度为12.5%时，能隙减小至2.1eV；当Nd掺杂浓度为25%时，能隙减小至2.0eV；当Nd掺杂浓度为37.5%时，能隙减小至1.9eV。这是因为Nd的掺入引入了新的电子态，使得价带和导带之间的能量差减小。从态密度分析，随着Nd掺杂浓度的增加，在费米能级附近出现了新的杂质能级。这些杂质能级主要来源于Nd的4f和5d电子态与BiFeO₃中其他原子电子态的相互作用。在Nd掺杂浓度为12.5%时，杂质能级位于费米能级下方约0.5eV处，其态密度相对较小；随着Nd掺杂浓度增加到25%，杂质能级向费米能级靠近，位于费米能级下方约0.3eV处，态密度有所增大；当Nd掺杂浓度达到37.5%时，杂质能级进一步靠近费米能级，位于费米能级下方约0.1eV处，态密度显著增大。这些杂质能级的出现改变了电子的跃迁路径和概率，从而影响了材料的电学性质，使得材料的导电性有所增强。4.2.2稀土掺杂对磁电耦合效应的调控稀土掺杂对BiFeO₃的磁电耦合效应具有重要的调控作用。研究表明，Nd掺杂促进了亚铁电调制结构的出现，这种结构对磁电耦合效应产生了显著影响。在未掺杂的BiFeO₃中，虽然具有铁电性和反铁磁性，但磁电耦合效应相对较弱。当Nd掺杂后，随着掺杂浓度的增加，亚铁电调制结构逐渐形成。在Nd掺杂浓度为12.5%时，开始出现亚铁电调制结构的迹象，此时材料的极化强度和磁电耦合系数开始发生变化；当Nd掺杂浓度增加到25%时，亚铁电调制结构更加明显，极化强度和磁电耦合系数进一步改变；当Nd掺杂浓度达到37.5%时，亚铁电调制结构稳定存在，极化强度和磁电耦合系数达到较大值。通过第一性原理计算结合Landau-Devonshire模型分析发现，极化和弱铁磁之间的耦合是通过氧八面体转动间接实现的。Nd掺杂导致晶格结构的变化，使得氧八面体的转动角度发生改变，从而影响了电子的自旋-轨道耦合以及磁矩的取向，进而实现了电场对弱铁磁的有效控制。在施加电场时，电场作用于亚铁电调制结构，通过氧八面体的转动，改变了磁矩的取向，从而实现了磁性的变化。这种磁电耦合效应的调控机制为开发基于BiFeO₃的低功耗自旋电子器件提供了理论基础，有望实现电场对磁性的高效调控，降低器件的能耗。4.2.3应用前景分析稀土掺杂BiFeO₃在低功耗自旋电子器件中展现出巨大的应用潜力。其独特的磁电耦合效应使得在电场作用下能够实现对磁性的有效控制，这一特性对于自旋电子器件的发展具有重要意义。在自旋电子学中，信息的存储和处理通常依赖于电子的自旋状态，而稀土掺杂BiFeO₃的磁电耦合特性可以实现通过电场来控制自旋状态，从而降低器件的功耗。在磁随机存取存储器（MRAM）中，传统的MRAM通常需要较大的电