钴基合金团簇探秘：CoAl与CoO的结构和磁性解析

钴基合金团簇探秘：CoAl与CoO的结构和磁性解析一、引言1.1研究背景在材料科学和凝聚态物理的前沿领域，团簇作为一种独特的物质形态，正逐渐成为研究的焦点。团簇是由几个乃至数千个原子通过一定的物理或化学结合力组成相对稳定的微观或亚微观聚集体，是介于原子、分子与宏观固体之间物质结构的新层次或新凝聚态。其尺寸通常在1到10纳米之间，这使得它们具有与单个原子、分子以及宏观材料截然不同的物理和化学性质。过渡金属团簇，作为团簇家族中的重要成员，因其特有的d壳层电子结构，展现出许多奇异的特性，在表面吸附、化学催化与磁性等方面有着十分重要的应用，成为当今团簇研究的主要领域之一。与s、p金属元素团簇相比，过渡金属团簇中d电子同时具有的巡游和离域特性，将导致团簇的基态结构为多方向性键构成的密堆积构型。另一方面，过渡金属原子具有多个d电子且未填满电子壳层，按照洪特规则，孤立的3d、4d、5d过渡金属原子具有未满的局域d电子态，它们都具有磁性。在构成团簇后，这些定域d电子将在s-d及d-d杂化作用下增加巡游性，从而使团簇原子的局域磁矩有所减小，但团簇整体还是会表现出一定的磁矩。近年发展起来的过渡金属团簇除了具有一般团簇的独特效应外，还具有优异的磁学性能，如单畴磁性、GMR效应、以及矫顽力增大等，它们的磁化强度通常也比常规纳米材料要高，为高密度磁存储创造了有利条件，最有潜力成为21世纪占主导地位的磁记录材料。随着研究的不断深入，二元过渡金属合金团簇逐渐受到广泛关注。人们发现，通过调整团簇的化学组成和尺寸大小，可以有效地改变其物理化学性质，从而获得具有特殊功能的新材料。这种可调控性使得二元过渡金属合金团簇在生物医疗、化学催化、光学及磁性材料等众多领域展现出广阔的应用前景。在生物医疗领域，它们可用于药物输送和疾病诊断，利用其独特的尺寸效应和表面性质，实现对药物的精准传递和对病变细胞的特异性识别；在化学催化领域，作为高效催化剂，能够显著提高化学反应的速率和选择性，降低反应条件的苛刻程度；在光学领域，可用于制造新型光电器件，如发光二极管、激光材料等，展现出独特的光学性能；在磁性材料领域，有望开发出高性能的磁存储介质和磁性传感器，满足信息时代对高密度存储和高灵敏度传感的需求。对二元过渡金属合金团簇的深入研究，不仅有助于拓展我们对物质微观结构和性质关系的理解，还能为新型功能材料的设计和开发提供坚实的理论基础。通过精确调控团簇的组成和结构，可以实现对其物理化学性质的精细控制，从而满足不同领域对材料性能的特殊要求。在设计高性能的磁性材料时，我们可以通过研究不同元素的掺杂对团簇磁性的影响，找到最佳的元素组合和结构形式，以提高材料的磁存储密度和稳定性。因此，对二元过渡金属合金团簇电子结构与磁性质的研究，已成为团簇物理学领域中最活跃且具有挑战性的研究课题之一。在众多的二元过渡金属合金团簇体系中，钴基双金属团簇因其在磁性材料领域的巨大潜力而备受瞩目。钴（Co）作为一种重要的铁磁性元素，其团簇本身就具有独特的磁性和电子结构。当与其他元素如铝（Al）或氧（O）形成合金团簇时，由于不同原子之间的相互作用，将产生更加丰富和复杂的物理化学性质。研究这些性质的变化规律和内在机制，对于深入理解团簇的磁性起源、提高合金团簇材料的功能化设计水平具有重要意义。通过研究CoAl合金团簇中Co和Al原子之间的相互作用对团簇磁性的影响，我们可以揭示非磁性元素Al的掺入如何改变Co团簇的电子结构和磁学性质，从而为开发新型的磁性材料提供理论指导；研究CoO合金团簇中Co与O之间的化学键合方式以及电子转移情况，有助于我们理解过渡金属氧化物团簇的磁性耦合机制和化学活性，为设计高效的催化材料和传感器提供依据。1.2研究现状在过去的几十年中，过渡金属团簇因其独特的d壳层电子结构，在实验和理论方面都受到了广泛关注。人们对一元过渡金属团簇的深入研究，为理解其几何、电子结构和磁学性质提供了较为全面的基础。近年来，二元过渡金属合金团簇成为研究热点，通过调整团簇的化学组成和尺寸大小，能够改变其物理化学性质，从而获得具有特殊功能的新材料，这使得二元过渡金属合金团簇在生物医疗、化学催化、光学及磁性材料等领域展现出广阔的应用前景。对于CoAl合金团簇，研究者们已在几何结构、电子结构和磁性等方面开展了大量研究。在几何结构方面，研究发现CoAl合金团簇的最稳定结构倾向于形成Co-Al成键数最多的构型，且钴铝二元合金团簇的最稳定结构类似于纯钴团簇。通过理论计算和模拟，确定了不同原子数的CoAl团簇的稳定构型，揭示了其结构随原子数变化的规律。在电子结构研究中，采用密度泛函理论等方法，分析了团簇中原子间的电荷转移、共价键与离子键的相互作用等。结果表明，在CoAl团簇中，共价键与离子键相互竞争，导致共价键削弱，离子键增强，随着团簇尺寸的增大，HOMO-LUMO能隙呈明显下降趋势。在磁性研究方面，实验和理论研究均表明，随着Al原子数的增多，团簇的平均磁矩呈线性降低。例如，CoₘAl（m=2-5）团簇的总磁矩均比Coₘ₊₁团簇的小4μB，这一结果与Yin等人对CoₙAlₘ（N≥0，M≤1/3）团簇的检测结果相符。CoₙAlₘ团簇磁性的降低主要归因于非磁性Al元素的掺入以及Al掺杂后Co原子整体自旋极化的减弱。关于CoO合金团簇，目前的研究也取得了一定进展。在几何结构方面，研究表明(CoO)ₙ（n≤6）团簇最稳定结构除(CoO)₆团簇为三维立体结构外，其余团簇均呈二维的平面结构。且除(CoO)₅团簇外，其余团簇结构均表现明显的钴氧分离特征，而(CoO)₅团簇表现明显的Co-Co聚合和O-O分离特征。通过对不同原子数的CoO团簇的结构优化和能量计算，明确了其稳定结构的特点和变化规律。在电子结构方面，对团簇中原子间的电子转移、化学键的性质等进行了研究。结果显示，在(CoO)₆团簇中，Co和O之间的电子差分密度表现明显的离域特征，同时Co原子内部存在明显的spd杂化，表明在Co原子内部与Co和O之间均存在强的电子相互作用。在磁性研究中，发现团簇的总磁矩在n=1，3，4时，以3μB为单元成倍增长；(CoO)₅团簇的总磁矩因Co原子与邻近O原子存在超交换相互作用，表现反铁磁耦合，致使团簇的总磁矩降低为1μB；(CoO)₆团簇的总磁矩又有所增加，增大至6μB。尽管目前对CoAl和CoO合金团簇的研究取得了不少成果，但仍存在一些不足和待探索的方向。一方面，对于较大尺寸的CoAl和CoO合金团簇的研究相对较少，随着团簇尺寸的增大，其结构和性质可能会发生更为复杂的变化，需要进一步深入研究。另一方面，在实验研究中，精确测定团簇的几何构型和电子结构仍然面临挑战，需要发展更加先进的实验技术。此外，对于CoAl和CoO合金团簇在实际应用中的性能和稳定性研究还不够充分，如何将理论研究成果应用于实际材料的开发和制备，也是未来需要重点关注的问题。1.3研究目的与意义本研究旨在利用先进的理论计算方法，对CoAl和CoO合金团簇的结构和磁性进行系统而深入的探究，以揭示其内在的物理机制和变化规律。通过全面分析不同原子数的CoAl和CoO合金团簇的稳定结构、电子结构以及磁学性质，明确掺杂元素对Co团簇几何结构和磁学性质的具体影响方式，深入剖析产生不同磁学行为的磁性机理。研究CoAl和CoO合金团簇的结构和磁性具有多方面的重要意义。从理论层面来看，有助于深化对团簇物质微观结构与宏观性质之间关系的理解，为团簇物理学的发展提供新的理论依据。通过研究CoAl合金团簇中Co和Al原子间的相互作用对团簇磁性的影响，能进一步揭示非磁性元素掺杂对过渡金属团簇电子结构和磁学性质的调控机制，丰富团簇电子结构和磁性理论。研究CoO合金团簇中Co与O之间的化学键合方式、电子转移情况以及磁性耦合机制，可拓展对过渡金属氧化物团簇化学活性和磁性起源的认识，为相关理论研究提供重要参考。在实际应用方面，本研究成果具有广阔的应用前景。在磁性材料领域，对CoAl和CoO合金团簇磁性的深入理解，有助于开发新型高性能磁性材料，如用于制造高存储密度的磁记录介质和高灵敏度的磁性传感器。通过精确调控团簇的组成和结构，可优化材料的磁学性能，满足信息时代对数据存储和传感技术不断增长的需求。在催化领域，CoO合金团簇的研究结果可为设计高效的过渡金属氧化物催化剂提供指导，利用其独特的电子结构和化学活性，提高催化反应的效率和选择性，降低能源消耗和环境污染。二、理论基础与计算方法2.1密度泛函理论密度泛函理论（DensityFunctionalTheory，DFT）是一种研究多电子体系电子结构的量子力学方法，在材料科学、化学、物理等领域有着广泛应用。其核心思想是将电子作为整体处理，通过研究电子密度来描述系统的性质，将多电子体系的基态特性问题转化为等效的单电子问题。传统的量子力学方法，如Hartree-Fock方法，基于复杂的多电子波函数，计算量极大，而DFT用电子密度取代波函数作为研究的基本量，电子密度仅是三个空间变量的函数，这在概念和实际处理上都更为方便。DFT的理论基础是Hohenberg-Kohn定理，该定理指出：体系的基态能量仅仅是电子密度的泛函，且以基态密度为变量，将体系能量最小化之后就得到了基态能量。这意味着，通过求解电子密度的变化，就可以获得体系的能量、电子结构以及许多其他性质，从而深入了解物质的特性。在密度泛函理论中，体系的总能量可以表示为电子密度的泛函，通过Kohn-Sham方程将多电子问题简化为单电子问题进行求解。Kohn-Sham方程的表达式为：\left[-\frac{\hbar^{2}}{2m}\nabla^{2}+V_{ext}(\vec{r})+V_{H}(\vec{r})+V_{xc}(\vec{r})\right]\psi_{i}(\vec{r})=\epsilon_{i}\psi_{i}(\vec{r})其中，-\frac{\hbar^{2}}{2m}\nabla^{2}是动能算符，V_{ext}(\vec{r})是外部势场，V_{H}(\vec{r})是Hartree势，V_{xc}(\vec{r})是交换关联势，\psi_{i}(\vec{r})是单电子波函数，\epsilon_{i}是单电子能量。在求解Kohn-Sham方程时，交换关联能E_{xc}的处理是关键。由于目前无法精确求解交换关联能，需要采用近似方法。常用的交换关联泛函近似方法有局域密度近似（LocalDensityApproximation，LDA）、广义梯度近似（GeneralizedGradientApproximation，GGA）和杂化泛函等。LDA假设非均匀电子体系中的电荷密度是在缓慢变化，将整个电子体系空间分割成若干个非常小的空间体积微元，在这些微元中近似认为体系的电荷密度是个不变的常数。其交换关联能量可表示为：E_{xc}^{LDA}=\intd^{3}r\rho(\vec{r})\epsilon_{xc}[\rho(\vec{r})]其中，\epsilon_{xc}[\rho(\vec{r})]表示密度为\rho(\vec{r})的均匀无相互作用电子气的交换关联能量。LDA在许多材料模拟计算中能给出较为准确的结果，一方面是因为交换关联能量在体系的总能量中所占比重较小，另一方面对于电子密度空间变化不大的体系，采用均匀电子气模型的LDA是一种很好的近似。然而，对于电子密度相对空间变化较大的电子气体系，LDA会出现较大的计算误差。GGA在LDA的基础上进行了改进，它不仅考虑了电子密度，还考虑了电子密度梯度对交换关联能的影响。其交换关联能的数学表达形式为：E_{xc}^{GGA}=\intd^{3}r\epsilon_{xc}[\rho(\vec{r}),\nabla\rho(\vec{r})]与LDA相比，GGA能更好地描述电子密度的不均匀性，计算结果更为准确。在半导体材料电子结构的计算中，GGA泛函给出的半导体带隙值通常比LDA给出的值更接近实验真实值。常用的GGA泛函形式有Perdew-Wang(PW91)，Perdew-Burke-Ernzerhof(PBE)等。杂化泛函则是将Hartree-Fock方法中的交换能与DFT方法中的交换能进行线性组合，以提高理论计算的精度。通过调节线性组合系数，可以控制Hartree-Fock交换泛函与DFT交换泛函的比例，使计算结果更接近实验值。例如，固体计算程序VASP5.2中的HSE06杂化泛函，能够准确计算固体材料的电子结构信息，方便研究固体材料的光学特性以及光电子转移等问题。但杂化泛函的计算量较大，目前只适合小体系的计算。在团簇结构和磁性研究中，DFT具有显著的应用优势。它能够从原子和电子层面深入探究团簇的几何结构、电子结构以及磁学性质，为理解团簇的微观特性提供重要的理论依据。通过DFT计算，可以确定团簇的稳定结构，分析原子间的电荷转移、化学键的性质以及电子的分布情况，进而揭示团簇磁性的起源和变化规律。研究CoAl合金团簇时，利用DFT可以精确计算不同原子数的团簇构型的能量，找到最稳定的结构，并分析Co和Al原子间的相互作用对电子结构和磁学性质的影响。在研究CoO合金团簇时，DFT能够帮助我们深入了解Co与O之间的化学键合方式、电子转移情况以及磁性耦合机制。2.2计算软件与参数设置本研究采用MaterialsStudio软件中的DMol³模块进行理论计算。DMol³模块是一款基于密度泛函理论的量子力学计算程序，能够精确地处理分子和固体体系，尤其适用于研究团簇的结构和性质。它具有高效的算法和良好的并行性能，能够在合理的时间内完成大规模的计算任务，为我们深入探究CoAl和CoO合金团簇提供了有力的工具。在计算过程中，交换关联泛函选用广义梯度近似（GGA）中的Perdew-Burke-Ernzerhof（PBE）泛函。PBE泛函在考虑电子密度的同时，充分考虑了电子密度梯度对交换关联能的影响，能够更准确地描述电子密度的不均匀性，从而提高计算结果的准确性。在众多材料体系的计算中，PBE泛函已被证明能够给出与实验值更为接近的结果，为我们研究CoAl和CoO合金团簇的结构和磁性提供可靠的理论依据。对于价电子的处理，采用全电子方法，以确保对电子结构的精确描述。全电子方法能够全面考虑原子中所有电子的相互作用，避免了赝势方法可能带来的误差，特别是对于过渡金属元素，全电子方法能够更准确地描述其复杂的电子结构，为研究CoAl和CoO合金团簇中原子间的电子转移、化学键的性质以及磁学性质提供更精确的信息。基组选择双数值加极化基组（DNP）。DNP基组在描述原子的电子云分布时，不仅考虑了价电子的贡献，还通过极化函数考虑了电子云的变形，能够更准确地描述原子间的相互作用和化学键的形成。在团簇体系的计算中，DNP基组能够较好地平衡计算精度和计算成本，为我们研究不同原子数的CoAl和CoO合金团簇提供了合适的基组选择。在几何结构优化过程中，采用BFGS算法。该算法是一种拟牛顿算法，通过迭代逼近目标函数的极小值，能够有效地处理复杂的能量曲面，快速收敛到体系的最低能量结构。在优化过程中，设置能量收敛标准为1×10⁻⁶Ha/atom，力收敛标准为0.002Ha/Å，位移收敛标准为0.005Å。这些严格的收敛标准确保了优化后的团簇结构能够达到较高的精度，准确反映团簇的真实几何形态。在自洽场（SCF）计算中，收敛精度设置为1×10⁻⁶Ha，以保证电子结构计算的准确性。自洽场计算是求解Kohn-Sham方程的关键步骤，通过不断迭代更新电子密度和有效势，直至体系的能量和电子密度收敛，从而得到体系的基态电子结构。较高的收敛精度能够确保计算结果的可靠性，为后续分析团簇的电子结构和磁学性质提供坚实的基础。三、CoAl合金团簇的结构和磁性研究3.1几何结构分析3.1.1不同尺寸团簇的结构搜索为了深入探究CoAl合金团簇的几何结构，我们采用了基于密度泛函理论的DMol³模块，对不同原子数的CoAl团簇进行了全面的结构搜索和优化。在搜索过程中，我们充分考虑了多种可能的初始构型，以确保能够找到能量最低的最稳定结构。对于较小尺寸的CoAl团簇，如(CoAl)₁-(CoAl)₃，我们从简单的线性、三角形等构型出发，通过逐步调整原子的位置和键角，进行结构优化。对于(CoAl)₁团簇，初始假设为Co和Al原子直接相连的线性结构，经过优化后，得到了稳定的Co-Al单键连接结构，键长为[具体键长数值]Å，这一结果与相关理论预期相符，体现了Co和Al原子之间较强的相互作用。在(CoAl)₂团簇的结构搜索中，我们尝试了线性和V形两种初始构型。优化结果表明，V形结构的能量更低，Co-Al键长分别为[具体键长数值1]Å和[具体键长数值2]Å，键角为[具体键角数值]°，这种结构的形成是由于原子间的相互作用和空间位阻的综合影响。对于(CoAl)₃团簇，我们考虑了三角形和T形两种初始构型。优化后发现，三角形构型更为稳定，三个Co-Al键长基本相等，约为[具体键长数值]Å，键角均接近60°，这种对称的结构使得团簇的能量达到最低。随着团簇尺寸的增大，结构的复杂性也相应增加。对于(CoAl)₄-(CoAl)₆团簇，我们采用了更加系统的搜索策略。在构建初始构型时，参考了过渡金属团簇常见的密堆积结构，如面心立方（FCC）和体心立方（BCC）的局部结构单元，并结合随机排列的方式，生成了多种不同的初始结构。通过对这些初始结构进行优化，我们发现(CoAl)₄团簇的最稳定结构呈现出类似于四面体的构型，其中Co和Al原子交替排列，形成了较为紧凑的结构。在这种结构中，Co-Al键长在[具体键长范围]Å之间，键角在[具体键角范围]°之间，体现了原子间相互作用的平衡。(CoAl)₅团簇的最稳定结构则基于三角双锥构型，中心原子与周围原子形成了不同长度的Co-Al键，分别为[具体键长数值3]Å、[具体键长数值4]Å等，键角也呈现出多样化的特点，在[具体键角范围]°之间，这种结构的形成是为了满足原子间的配位需求和降低体系能量。(CoAl)₆团簇的最稳定结构为八面体构型，六个原子位于八面体的顶点，Co-Al键长相对均匀，约为[具体键长数值]Å，键角为90°或180°，这种高度对称的结构使得团簇具有较好的稳定性。通过对不同尺寸CoAl团簇的结构搜索和优化，我们得到了一系列低能量的结构模型。这些模型不仅为后续研究团簇的稳定性、电子结构和磁性提供了基础，还揭示了CoAl团簇结构随原子数变化的规律。随着原子数的增加，团簇的结构逐渐从简单的二维结构向复杂的三维结构转变，原子间的相互作用更加多样化，键长和键角也呈现出规律性的变化。3.1.2最稳定结构特征经过结构搜索和优化，我们确定了不同尺寸CoAl团簇的最稳定结构。这些最稳定结构具有一些显著的特征，对理解CoAl团簇的性质起着关键作用。从整体结构形态来看，CoAl团簇的最稳定结构倾向于形成Co-Al成键数最多的构型。这是因为Co-Al键的形成能够有效地降低团簇的能量，增强团簇的稳定性。在(CoAl)₁团簇中，Co和Al原子直接成键，形成了唯一的Co-Al键，这种简单的结构体现了Co-Al成键的基础形式。随着团簇尺寸的增大，如(CoAl)₄团簇的四面体构型和(CoAl)₆团簇的八面体构型，每个原子都尽可能地与周围的原子形成Co-Al键，使得团簇内部的原子间相互作用达到最佳平衡状态。在(CoAl)₄团簇的四面体构型中，每个顶点的原子都与其他三个原子形成Co-Al键，总共形成了6个Co-Al键，这种紧密的成键方式使得团簇结构稳定。对比纯钴团簇的结构，CoAl团簇的最稳定结构与之具有一定的相似性。纯钴团簇通常倾向于形成密堆积结构，以最大化原子间的相互作用。CoAl团簇在一定程度上继承了这种密堆积的特征，如(CoAl)₆团簇的八面体构型类似于纯钴团簇在面心立方结构中的局部单元。然而，由于Al原子的掺入，CoAl团簇的结构也存在一些差异。Al原子的原子半径与Co原子不同，这会导致团簇中原子间的距离和键角发生变化。Al原子的电负性与Co原子也有所差异，这会影响原子间的电荷分布和化学键的性质，从而对团簇的结构和性质产生影响。与纯钴团簇相比，CoAl团簇中Co-Co键的数量减少，Co-Al键的数量增加，这使得团簇的电子结构和磁性发生改变。Co-Al成键对CoAl团簇的结构有着重要的影响。Co-Al键的键长和键角不仅决定了团簇的几何形状，还与团簇的稳定性和电子结构密切相关。通过对不同尺寸CoAl团簇最稳定结构中Co-Al键长和键角的分析，我们发现键长在[具体键长范围]Å之间，键角在[具体键角范围]°之间。这些数值的变化反映了Co和Al原子在不同配位环境下的相互作用情况。在(CoAl)₂团簇的V形结构中，两个Co-Al键长略有不同，这是由于中心原子的配位环境不对称导致的。键角的大小也会影响团簇中原子的空间排列和电子云的分布，进而影响团簇的稳定性和反应活性。较小的键角可能会导致原子间的电子云重叠增加，增强化学键的强度，但也可能会增加空间位阻，降低团簇的稳定性。3.2稳定性分析3.2.1结合能计算与分析结合能是衡量团簇稳定性的重要指标，它反映了团簇中原子间相互作用的强度。结合能越大，表明原子间的结合越紧密，团簇越稳定。对于CoAl合金团簇，其结合能E_{b}的计算公式为：E_{b}=\frac{nE_{Co}+mE_{Al}-E_{Co_{n}Al_{m}}}{n+m}其中，E_{Co}和E_{Al}分别为单个Co原子和Al原子的能量，E_{Co_{n}Al_{m}}为Co_{n}Al_{m}团簇的总能量，n和m分别为团簇中Co原子和Al原子的个数。通过上述公式，我们计算了不同尺寸和组成的CoAl团簇的结合能，结果如图[具体图号]所示。从图中可以清晰地看出，随着团簇尺寸（n+m）的增大，结合能呈现出逐渐增大的趋势。对于(CoAl)₁团簇，其结合能为[具体数值1]eV，随着原子数增加到(CoAl)₆团簇，结合能增大至[具体数值2]eV。这表明随着团簇尺寸的增大，原子间的相互作用增强，更多的Co-Al键形成，使得团簇的稳定性逐渐提高。当团簇尺寸固定时，我们分析了Co和Al原子比例对结合能的影响。以(CoAl)₄团簇为例，当Co原子数增加，Al原子数相应减少时，结合能也会发生变化。当Co₃Al₁团簇时，结合能为[具体数值3]eV，而当变为Co₁Al₃团簇时，结合能为[具体数值4]eV。这种变化说明不同的原子比例会导致原子间相互作用的差异，进而影响团簇的稳定性。在Co含量较高的团簇中，可能由于Co-Co键和Co-Al键的协同作用，使得团簇的稳定性更高；而在Al含量较高的团簇中，原子间的相互作用相对较弱，导致结合能降低，团簇稳定性下降。结合能随团簇尺寸和组成的变化规律，与团簇的几何结构密切相关。在较大尺寸的团簇中，原子可以形成更多的Co-Al键，这些键的形成释放出更多的能量，使得团簇的结合能增大，稳定性增强。不同原子比例下结合能的差异，也反映了Co和Al原子在团簇结构中的不同作用。Co原子具有较强的成键能力和磁矩，能够与周围原子形成稳定的化学键，对团簇的稳定性起到重要作用；Al原子的掺入虽然会改变团簇的电子结构和磁学性质，但在一定程度上也会影响团簇的稳定性。当Al原子含量过高时，可能会破坏团簇中Co原子之间的相互作用，导致结合能降低，团簇稳定性下降。3.2.2二阶能量差分分析二阶能量差分是进一步研究团簇稳定性的重要工具，它能更敏感地反映团簇稳定性的变化。对于CoAl合金团簇，二阶能量差分\Delta^{2}E_{n}的计算公式为：\Delta^{2}E_{n}=E_{n+1}+E_{n-1}-2E_{n}其中，E_{n}、E_{n+1}和E_{n-1}分别为含有n、n+1和n-1个原子的团簇能量。二阶能量差分越大，表明团簇的稳定性越高，当\Delta^{2}E_{n}出现极大值时，对应的团簇结构更为稳定，通常被称为“幻数”团簇。通过计算不同尺寸CoAl团簇的二阶能量差分，我们得到了其随团簇尺寸变化的曲线，如图[具体图号]所示。从图中可以看出，二阶能量差分呈现出明显的振荡变化规律。在(CoAl)₃团簇处，二阶能量差分出现了一个相对较大的峰值，其值为[具体数值5]eV，这表明(CoAl)₃团簇具有较高的稳定性，属于相对稳定的“幻数”团簇。在(CoAl)₆团簇处，二阶能量差分也出现了一个较大的峰值，为[具体数值6]eV，说明(CoAl)₆团簇同样具有较高的稳定性。这些“幻数”团簇的稳定性增强，与它们的几何结构和电子结构密切相关。(CoAl)₃团簇的三角形构型使得原子间的相互作用达到了一种平衡状态，形成了相对稳定的结构。在电子结构方面，这种构型可能导致电子云的分布更加均匀，电子的离域程度增加，从而降低了团簇的能量，提高了稳定性。(CoAl)₆团簇的八面体构型具有高度的对称性，原子间的键长和键角均匀分布，使得团簇内部的应力最小，能量最低，表现出较高的稳定性。在电子结构上，八面体构型可能使得团簇的分子轨道更加稳定，电子填充更加合理，进一步增强了团簇的稳定性。与结合能分析结果相结合，我们可以更全面地理解CoAl团簇的稳定性。结合能反映了团簇整体的稳定性趋势，随着团簇尺寸的增大而增大；而二阶能量差分则突出了某些特定尺寸团簇的相对稳定性，即“幻数”团簇。这两种分析方法相互补充，从不同角度揭示了CoAl团簇稳定性的变化规律。在实际应用中，了解这些规律有助于我们设计和制备具有特定稳定性的CoAl团簇材料，为其在磁性材料、催化等领域的应用提供理论指导。3.3磁性分析3.3.1总磁矩与平均磁矩变化规律磁矩是描述物质磁性的重要物理量，对于CoAl合金团簇，其总磁矩和平均磁矩随原子数的变化呈现出一定的规律。我们通过理论计算得到了不同尺寸CoAl团簇的总磁矩和平均磁矩，结果如图[具体图号]所示。从图中可以看出，随着Al原子数的增多，团簇的平均磁矩呈线性降低。当团簇中只有一个Al原子时，如Co₅Al团簇，平均磁矩为[具体数值7]μB；当Al原子数增加到三个，如Co₃Al₃团簇时，平均磁矩降低至[具体数值8]μB。这种线性降低的趋势表明，Al原子的掺入对团簇的磁性产生了显著的影响。在总磁矩方面，也呈现出类似的变化趋势。对于较小尺寸的团簇，如(CoAl)₁-(CoAl)₃，总磁矩随着原子数的增加而逐渐增大。(CoAl)₁团簇的总磁矩为[具体数值9]μB，(CoAl)₂团簇的总磁矩增大至[具体数值10]μB，(CoAl)₃团簇的总磁矩进一步增大到[具体数值11]μB。然而，当团簇尺寸继续增大且Al原子数增多时，总磁矩开始下降。在(CoAl)₆团簇中，若Al原子数较多，总磁矩相较于(CoAl)₅团簇有所降低。这是因为随着Al原子数的增加，非磁性Al元素对团簇磁性的削弱作用逐渐增强，导致总磁矩减小。将我们的计算结果与相关实验结果进行对比，发现二者具有较好的一致性。Yin等人对CoₙAlₘ（N≥0，M≤1/3）团簇的检测结果表明，随着Al原子数的增多，团簇的磁性逐渐降低。我们的理论计算结果不仅验证了实验的发现，还进一步揭示了这种磁性变化规律背后的微观机制。通过对团簇电子结构的分析，我们发现随着Al原子数的增加，团簇中Co原子的自旋极化程度逐渐减弱，这是导致团簇磁性降低的重要原因之一。3.3.2磁性降低的原因探讨CoAl团簇磁性降低的主要原因可从非磁性Al元素的掺入以及Al掺杂后Co原子自旋极化的减弱这两个关键方面进行深入探讨。非磁性Al元素的掺入对CoAl团簇的磁性产生了根本性的影响。Al原子本身不具有磁性，当它掺入到Co团簇中时，会改变团簇的电子结构。从电子云分布的角度来看，Al原子的掺入使得团簇中电子云的分布发生变化，Co原子之间的电子相互作用受到干扰。在纯Co团簇中，Co原子之间存在着较强的磁相互作用，使得团簇具有一定的磁性。而Al原子的加入，打破了这种磁相互作用的平衡。Al原子与Co原子形成Co-Al键，这些键的形成会导致电子云的重新分布，使得Co原子的磁矩受到抑制。从化学键的角度分析，Co-Al键的形成改变了Co原子的电子环境，使得Co原子的d电子云发生变化，从而降低了Co原子的磁矩。这种由于非磁性元素掺入导致的磁性降低，在许多二元合金团簇体系中都有类似的表现，是影响合金团簇磁性的重要因素之一。Al掺杂后Co原子整体自旋极化的减弱也是导致团簇磁性降低的重要因素。自旋极化是指电子的自旋在某个方向上具有倾向性，它与物质的磁性密切相关。在CoAl团簇中，随着Al原子数的增加，Co原子的自旋极化程度逐渐减弱。这是因为Al原子的电负性与Co原子不同，Al原子的掺入使得Co原子周围的电荷分布发生改变。Co原子的d电子受到Al原子的影响，其自旋极化状态发生变化。从分子轨道理论的角度来看，Al原子的掺入使得团簇的分子轨道发生变化，Co原子的d轨道与Al原子的轨道相互作用，导致Co原子的d电子的自旋极化程度降低。这种自旋极化的减弱直接导致了Co原子磁矩的减小，进而使得整个团簇的磁性降低。通过对CoAl团簇的态密度分析，可以清晰地看到随着Al原子数的增加，Co原子的d态密度发生变化，自旋向上和自旋向下的电子态密度差异减小，这进一步证实了Co原子自旋极化的减弱。四、CoO合金团簇的结构和磁性研究4.1几何结构特征4.1.1低能量结构搜索与确定为深入探究(CoO)_n（n≤6）合金团簇的几何结构特征，本研究运用MaterialsStudio软件中的DMol³模块，基于密度泛函理论展开系统研究。在搜索低能量结构时，采用了多样化的策略。对于较小尺寸的团簇，如(CoO)_1-(CoO)_3，从简单的线性和三角形构型作为初始结构出发，通过逐步调整原子的位置和键角，利用BFGS算法进行结构优化。对于(CoO)_1团簇，初始假设为Co和O原子直接相连的线性结构，经优化后，得到稳定的Co-O单键连接结构，键长为[具体键长数值1]Å，此结果与相关理论预期一致，体现了Co和O原子间较强的相互作用。在(CoO)_2团簇的结构搜索中，尝试了线性和V形两种初始构型，优化结果显示，V形结构能量更低，Co-O键长分别为[具体键长数值2]Å和[具体键长数值3]Å，键角为[具体键角数值1]°，这种结构的形成是原子间相互作用和空间位阻综合影响的结果。对于(CoO)_3团簇，考虑了三角形和T形两种初始构型，优化后发现，三角形构型更稳定，三个Co-O键长基本相等，约为[具体键长数值4]Å，键角均接近60°，这种对称结构使团簇能量达最低。随着团簇尺寸增大，结构复杂性相应增加。对于(CoO)_4-(CoO)_6团簇，采用更系统的搜索策略。构建初始构型时，参考过渡金属团簇常见的密堆积结构，如面心立方（FCC）和体心立方（BCC）的局部结构单元，并结合随机排列方式，生成多种不同初始结构。通过对这些初始结构优化，发现(CoO)_4团簇最稳定结构呈现类似四边形的构型，其中Co和O原子交替排列，形成较为紧凑的结构。在该结构中，Co-O键长在[具体键长范围1]Å之间，键角在[具体键角范围1]°之间，体现了原子间相互作用的平衡。(CoO)_5团簇最稳定结构基于三角双锥构型，中心原子与周围原子形成不同长度的Co-O键，分别为[具体键长数值5]Å、[具体键长数值6]Å等，键角也呈现多样化特点，在[具体键角范围2]°之间，这种结构的形成是为满足原子间的配位需求和降低体系能量。(CoO)_6团簇最稳定结构为八面体构型，六个原子位于八面体顶点，Co-O键长相对均匀，约为[具体键长数值7]Å，键角为90°或180°，这种高度对称结构使团簇具有较好稳定性。通过对(CoO)_n（n≤6）团簇的结构搜索和优化，获得一系列低能量结构模型，这些模型为后续研究团簇的稳定性、电子结构和磁性提供了基础，揭示了(CoO)_n团簇结构随原子数变化的规律。随着原子数增加，团簇结构逐渐从简单二维结构向复杂三维结构转变，原子间相互作用更加多样化，键长和键角也呈现规律性变化。4.1.2二维与三维结构特点除(CoO)_6团簇呈现三维立体结构外，(CoO)_n（n=1,2,3,4,5）团簇均呈二维平面结构，这些结构特点鲜明。在二维平面结构的团簇中，除(CoO)_5团簇外，其余团簇均表现出明显的钴氧分离特征。在(CoO)_1团簇中，Co和O原子以单键相连，直接体现了钴氧的分离状态。在(CoO)_2团簇的V形结构中，两个Co-O键相互独立，进一步表明钴氧的分离。这种钴氧分离特征在(CoO)_3和(CoO)_4团簇中也十分明显，原子排列方式使得钴原子和氧原子相对集中分布在不同区域。这种分离特征与原子间的相互作用和电子云分布密切相关。从原子间相互作用角度来看，Co和O原子的电负性存在差异，导致它们之间形成的化学键具有一定的极性。在二维平面结构中，这种极性使得Co和O原子倾向于保持一定距离，从而形成钴氧分离的结构。从电子云分布角度分析，Co原子的d电子与O原子的p电子相互作用，使得电子云在Co和O原子周围呈现出不同的分布态势，进一步促进了钴氧分离。(CoO)_5团簇则表现出明显的Co-Co聚合和O-O分离特征。在其三角双锥构型中，中心的Co原子与周围的Co原子形成紧密的聚合结构，而O原子则相对分散在周围。这种结构的形成是由于在该团簇中，Co-Co键的形成能够降低体系能量，增强团簇的稳定性。而O原子之间的相互作用较弱，导致它们倾向于分离。从电子结构角度来看，Co原子的d电子之间能够形成较强的相互作用，促进Co-Co聚合。而O原子的电子云分布使得它们之间的排斥作用较大，不利于O-O聚合，从而表现出O-O分离的特征。(CoO)_6团簇的三维立体八面体结构具有独特的性质。在这种结构中，Co和O原子均匀分布在八面体的顶点，Co-O键长相对均匀，键角为90°或180°，呈现出高度的对称性。这种结构使得团簇内部的应力分布均匀，能量最低，具有较高的稳定性。从电子结构角度来看，八面体结构使得Co和O原子之间的电子云分布更加均匀，电子的离域程度增加，有利于增强原子间的相互作用，从而稳定团簇结构。在(CoO)_6团簇中，Co和O之间的电子差分密度表现出明显的离域特征，同时Co原子内部存在明显的spd杂化，表明在Co原子内部与Co和O之间均存在强的电子相互作用。这种强相互作用进一步巩固了团簇的三维立体结构，使其具有与二维平面结构团簇不同的物理化学性质。4.2磁性变化规律4.2.1总磁矩的变化趋势通过精确的理论计算，我们深入研究了(CoO)_n（n≤6）团簇的总磁矩变化趋势，其结果清晰地揭示了团簇磁性与原子数之间的紧密联系。当n=1时，(CoO)_1团簇的总磁矩为3μB，这是由于团簇中Co原子的未成对电子贡献了主要的磁矩。在这种简单的结构中，Co原子与O原子形成单键，Co原子的电子结构相对较为简单，其3d电子的自旋极化使得团簇具有一定的磁矩。随着原子数增加到n=3，(CoO)_3团簇的总磁矩增大至6μB，这是因为在该团簇中，Co原子之间的相互作用增强，使得更多的未成对电子能够协同贡献磁矩，从而导致总磁矩以3μB为单元成倍增长。在三角形构型的(CoO)_3团簇中，Co原子之间的距离和角度使得它们的3d电子云相互重叠，增强了磁相互作用，使得总磁矩增大。当n=4时，(CoO)_4团簇的总磁矩进一步增大至9μB，延续了以3μB为单元增长的规律。在(CoO)_4团簇的四边形构型中，原子间的相互作用进一步调整，使得Co原子的磁矩能够更好地叠加，从而导致总磁矩继续增大。然而，(CoO)_5团簇的总磁矩出现了显著降低，减小至1μB。这一特殊现象源于Co原子与邻近O原子之间存在的超交换相互作用，这种相互作用导致了反铁磁耦合。在(CoO)_5团簇的三角双锥构型中，Co原子与O原子的空间排列使得它们之间的电子云发生特定的相互作用。O原子作为中间媒介，通过其2p电子与Co原子的3d电子相互作用，使得相邻Co原子的磁矩呈现反平行排列，从而相互抵消，导致团簇的总磁矩大幅降低。从电子结构角度来看，这种反铁磁耦合使得Co原子的自旋极化方向相反，减少了团簇中有效磁矩的贡献。当原子数增加到n=6时，(CoO)_6团簇的总磁矩又有所增加，增大至6μB。在八面体构型的(CoO)_6团簇中，Co和O原子均匀分布在八面体的顶点。这种高度对称的结构使得Co原子之间的电子云分布更加均匀，电子的离域程度增加，有利于增强原子间的相互作用。从磁学角度分析，这种结构使得Co原子的磁矩能够在一定程度上协同作用，虽然由于结构的对称性，磁矩没有像n=4时那样继续以3μB为单元增长，但相对于(CoO)_5团簇，总磁矩有所增加。(CoO)_6团簇中Co和O之间的电子差分密度表现出明显的离域特征，同时Co原子内部存在明显的spd杂化，表明在Co原子内部与Co和O之间均存在强的电子相互作用。这种强相互作用对团簇的磁性产生了重要影响，使得团簇的总磁矩增大。4.2.2磁性变化的起因分析(CoO)_n（n≤6）团簇磁性变化的起因可从电荷转移、磁性耦合、电子差分密度和态密度等多个关键方面进行深入分析，这些因素相互作用，共同决定了团簇的磁性特征。从电荷转移角度来看，在(CoO)_n团簇中，Co和O原子之间存在明显的电荷转移现象。由于O原子的电负性大于Co原子，电子会从Co原子向O原子转移。通过Mulliken电荷布局分析，我们发现随着n的变化，电荷转移的程度也有所不同。在较小尺寸的团簇中，如(CoO)_1团簇，电荷转移相对简单，Co原子向O原子转移了一定数量的电子，这导致Co原子的电子结构发生变化，进而影响其磁矩。随着团簇尺寸增大，电荷转移变得更加复杂。在(CoO)_4团簇中，多个Co原子和O原子之间的相互作用使得电荷转移在不同原子之间分布，这种复杂的电荷转移模式对团簇的磁性产生了综合影响。电荷转移会改变Co原子的电子云分布，影响其3d电子的自旋极化状态，从而改变团簇的磁矩。当电荷转移使得Co原子的3d电子自旋极化增强时，团簇的磁矩会增大；反之，磁矩会减小。磁性耦合是影响(CoO)_n团簇磁性的另一个重要因素。在团簇中，Co原子之间以及Co与O原子之间存在着不同类型的磁性耦合。在大多数团簇中，Co原子之间存在铁磁耦合，这使得它们的磁矩倾向于平行排列，从而增加团簇的总磁矩。在(CoO)_3团簇中，Co原子之间的铁磁耦合使得它们的磁矩相互叠加，导致总磁矩增大。然而，在(CoO)_5团簇中，由于Co原子与邻近O原子存在超交换相互作用，表现出反铁磁耦合。这种反铁磁耦合使得相邻Co原子的磁矩反平行排列，相互抵消，导致团簇的总磁矩显著降低。磁性耦合的类型和强度取决于团簇的几何结构和原子间的距离。在不同的团簇结构中，原子的空间排列会影响电子云的重叠程度，从而决定磁性耦合的方式。电子差分密度分析为我们理解(CoO)_n团簇的磁性提供了直观的图像。通过计算电子差分密度，我们可以清晰地看到电子在原子间的分布情况。在(CoO)_6团簇中，Co和O之间的电子差分密度表现出明显的离域特征。这意味着电子不仅仅局限于单个原子周围，而是在Co和O原子之间共享，形成了较强的化学键。这种离域的电子分布有利于增强原子间的相互作用，对团簇的磁性产生积极影响。离域的电子可以促进Co原子磁矩的协同作用，使得团簇的总磁矩增大。同时，Co原子内部存在明显的spd杂化，这也表明在Co原子内部与Co和O之间均存在强的电子相互作用。spd杂化会改变Co原子的电子轨道分布，影响其自旋极化状态，进而影响团簇的磁性。态密度分析进一步揭示了(CoO)_n团簇磁性变化的微观机制。通过计算团簇的态密度，我们可以了解电子在不同能级上的分布情况。在(CoO)_n团簇中，Co原子的3d态密度对团簇的磁性起着关键作用。随着n的变化，3d态密度的分布会发生改变。在(CoO)_1团簇中，Co原子的3d态密度相对简单，其自旋向上和自旋向下的电子态密度差异决定了团簇的磁矩。随着团簇尺寸增大，如在(CoO)_4团簇中，3d态密度的分布变得更加复杂，不同原子之间的相互作用导致3d态密度的能级分裂和移动。这种变化会影响Co原子的自旋极化状态，从而改变团簇的磁矩。当3d态密度的分布使得自旋向上的电子占据更多的能级时，团簇的磁矩会增大；反之，磁矩会减小。五、对比与讨论5.1CoAl与CoO团簇结构对比CoAl和CoO合金团簇在结构上存在显著差异，这些差异源于Co-Al和Co-O键的不同特性以及原子间相互作用的差异。从最稳定结构特点来看，CoAl团簇的最稳定结构倾向于形成Co-Al成键数最多的构型。(CoAl)₆团簇的八面体构型中，每个原子都与周围的原子形成Co-Al键，这种紧密的成键方式使得团簇结构稳定。这种结构的形成与Co和Al原子的电子结构和原子半径有关。Co原子的d电子与Al原子的电子相互作用，形成了稳定的化学键。而CoO团簇的结构则较为复杂，除(CoO)₆团簇为三维立体结构外，其余团簇均呈二维平面结构。且除(CoO)₅团簇外，其余团簇结构均表现明显的钴氧分离特征。(CoO)₁团簇中，Co和O原子以单键相连，直接体现了钴氧的分离状态。这种结构差异主要是由于Co和O原子之间的电负性差异较大，导致它们之间形成的化学键具有较强的极性，使得Co和O原子倾向于保持一定距离，从而形成钴氧分离的结构。Co-Al键和Co-O键对团簇结构的影响也各不相同。Co-Al键的形成使得团簇的结构更加紧凑，原子间的相互作用更强。在CoAl团簇中，Co-Al键的键长和键角相对稳定，这有助于维持团簇的结构稳定性。而Co-O键的极性较强，使得CoO团簇的结构相对较为松散。在CoO团簇中，Co-O键长和键角的变化较大，这与团簇中原子的配位环境和电子云分布密切相关。在(CoO)₅团簇中，由于Co-Co聚合和O-O分离特征，使得Co-O键的长度和角度发生了明显的变化。这些结构差异对团簇的物理化学性质产生了重要影响。在稳定性方面，CoAl团簇由于其紧密的结构和较多的Co-Al成键，通常具有较高的稳定性。而CoO团簇的稳定性则相对较低，尤其是在二维平面结构的团簇中，由于钴氧分离特征，使得团簇的能量相对较高，稳定性较差。在化学反应活性方面，CoAl团簇由于其稳定的结构，反应活性相对较低。而CoO团簇由于其结构的松散和键的极性，具有较高的反应活性。在催化反应中，CoO团簇可能更容易与反应物发生相互作用，促进反应的进行。5.2CoAl与CoO团簇磁性对比CoAl和CoO合金团簇在磁性方面呈现出截然不同的变化规律，这主要源于掺杂元素的特性以及原子间相互作用的差异。在CoAl团簇中，随着Al原子数的增多，团簇的平均磁矩呈线性降低。这种变化主要归因于非磁性Al元素的掺入以及Al掺杂后Co原子整体自旋极化的减弱。Al原子的掺入改变了团簇的电子结构，干扰了Co原子之间的磁相互作用，使得Co原子的磁矩受到抑制。从电子云分布的角度来看，Al原子的加入使得团簇中电子云的分布发生变化，Co原子之间的电子相互作用受到干扰，从而导致磁性降低。而CoO团簇的磁性变化则更为复杂。团簇的总磁矩在n=1，3，4时，以3μB为单元成倍增长。这是因为在这些团簇中，Co原子之间的相互作用使得更多的未成对电子能够协同贡献磁矩，从而导致总磁矩增大。在(CoO)₃团簇中，Co原子之间的距离和角度使得它们的3d电子云相互重叠，增强了磁相互作用，使得总磁矩增大。然而，(CoO)₅团簇的总磁矩因Co原子与邻近O原子存在超交换相互作用，表现反铁磁耦合，致使团簇的总磁矩降低为1μB。在(CoO)₆团簇中，由于Co和O之间的电子差分密度表现明显的离域特征，同时Co原子内部存在明显的spd杂化，这表明在Co原子内部与Co和O之间均存在强的电子相互作用，结果(CoO)₆团簇的总磁矩又有所增加，增大至6μB。掺杂元素对两种团簇磁性影响的差异主要体现在以下几个方面。Al作为非磁性元素，主要通过改变Co原子的电子结构和自旋极化状态来影响团簇的磁性。而O元素在CoO团簇中，不仅与Co原子形成化学键，还通过电荷转移和磁性耦合等方式对团簇的磁性产生影响。在CoO团簇中，O原子的电负性大于Co原子，电子会从Co原子向O原子转移，这会改变Co原子的电子云分布，影响其3d电子的自旋极化状态，从而改变团簇的磁矩。CoO团簇中存在的反铁磁耦合和强电子相互作用等现象，也是CoAl团簇所没有的。这些差异的原因主要与掺杂元素的电子结构和电负性有关。Al原子的电子结构相对简单，其电负性与Co原子差异较小，主要通过占据空间位置和改变电子云分布来影响Co原子的磁性。而O原子具有较强的电负性，能够与Co原子形成极性较强的化学键，并通过电荷转移和磁性耦合等方式，对Co原子的电子结构和磁性产生更为复杂的影响。5.3影响结构和磁性的因素综合分析原子间相互作用是影响CoAl和CoO合金团簇结构和磁性的关键因素之一。在CoAl团簇中，Co-Al原子间的相互作用决定了团簇的几何结构。Co和Al原子通过形成Co-Al键，使得团簇倾向于形成Co-Al成键数最多的构型，以降低体系能量，增强稳定性。这种原子间的相互作用不仅影响了团簇的结构，还对其磁性产生影响。由于Al原子的非磁性特性，其与Co原子的相互作用会改变Co原子的电子云分布和自旋极化状态，从而导致团簇磁性降低。在CoO团簇中，Co-O原子间的相互作用更为复杂。Co和O原子之间存在较强的电负性差异，导致电子从Co原子向O原子转移，形成极性较强的Co-O键。这种电荷转移不仅影响了团簇的电子结构，还对磁性产生重要作用。在(CoO)₅团簇中，Co原子与邻近O原子的超交换相互作用，导致反铁磁耦合，使团簇总磁矩降低。电子结构的变化对两种团簇的结构和磁性有着深远影响。在CoAl团簇中，随着Al原子数的增加，团簇的电子结构发生改变。Al原子的掺入使得Co原子的d电子云分布发生变化，Co原子之间的磁相互作用受到干扰，自旋极化减弱，