曲酸锌(II)、镉(II)与汞(II)配合物的结构、电子性质与热力学特性：光谱学与密度泛函理论联合解析

《Kuwait Journal of Science》：Structural, electronic, and thermodynamic insights into Zn(II), Cd(II), and Hg(II) kojate systems: combined spectroscopic and DFT analysis

　　

：研究人员合成了曲酸的锌(II)、镉(II)和汞(II)配合物，并通过光谱分析与热分析对其进行了系统研究。红外光谱（FTIR）证实了金属离子与配体羟基氧及羰基氧原子之间存在配位作用。与游离曲酸相比，所有配合物均表现出非电解质特性，且热稳定性显著提高。配合物的电子与磁性数据表明其主要呈现四面体构型。密度泛函理论（DFT）计算为四面体畸变提供了证据，并得出与实验数据一致的热稳定性顺序：锌(II) > 汞(II) > 镉(II)。

论文解读文章

一、研究背景与问题

曲酸（kojic acid），学名5-羟基-2-羟甲基-4H-吡喃-4-酮，是一种γ-吡喃酮衍生物，在生物技术、化妆品和制药领域具有广泛应用。曲酸分子中，5-羟基与4-羰基在羟基去质子化后形成特征性的O,O-供体配位位点，使其能够与多种金属离子发生配位作用。锌(II)、镉(II)和汞(II)同属元素周期表第12族，均具有完整的d¹⁰电子构型，但三者离子尺寸、极化率、生物功能和毒性存在显著差异，这些差异导致其在与配体配位时表现出不同的金属-氧键距离、电荷再分布和配体振动频率。

此前针对曲酸的计算研究主要集中于游离曲酸分子本身，但计算所得分子与实验分离所得固体的组成和结构未必一致。这一矛盾促使研究人员对既有实验数据档案进行重新分析，并将其与在相同密度泛函理论（DFT）水平下计算的1:2型锌(II)和汞(II)曲酸盐模型进行系统比较。研究旨在通过分析闭合度、碳-金属比、表观曲酸根-金属键比、稳定性、配体与金属间距离、缔合能、前线轨道能隙、布居描述符等参数，评估理论模型与实验数据之间的对应关系。

二、研究内容与主要结论

研究人员对锌(II)、镉(II)和汞(II)的曲酸配合物（分别记为Zn-K、Cd-K和Hg-K）进行了合成与表征，并开展了系统的实验测定与DFT计算。实验方面包括元素分析、磁化率测定、漫反射光谱、傅里叶变换红外光谱（FTIR）和热重分析。计算方面采用Psi4软件，在PBE0泛函结合D3色散校正和Becke-Johnson阻尼、def2-SVP基组/有效核势（ECP）水平下，对三个中性[M(koj)₂]模型（M = Zn、Cd、Hg）进行了几何优化和电子结构计算。

实验数据重新分析表明，三种产物的分析闭合度均不完整，Zn-K、Cd-K和Hg-K的未计入部分分别为18.12%、11.21%和8.84%。碳-金属比分别为10.75、17.26和30.29，对应的曲酸根当量与金属之比为1.79、2.88和5.05，均无法用1:2化学计量比解释。理论计算值显示，若为[M(koj)₂]组成，锌、镉、汞的金属含量应分别为18.81%、28.49%和41.55%，而实测值仅为17.33%、20.12%和21.57%，存在明显偏差。因此，实验数据不支持三种固体产物具有统一的[M(koj)₂]分子式。

FTIR光谱分析显示，与游离曲酸相比，三种配合物的O-H和C=O吸收峰均向低波数方向移动，且位移幅度按Zn-K、Cd-K、Hg-K顺序递增，表明金属配位对配体含氧基团产生了不同程度的扰动，但光谱数据本身无法确定质子化状态、供体类型和产物化学计量比。

热重分析重构数据表明，曲酸、Zn-K、Cd-K和Hg-K的T₅₀值（质量保留50%对应的温度）分别为264.3°C、332.5°C、334.4°C和335.7°C，配合物较游离曲酸的T₅₀提高了约68-71°C，热稳定性显著增强。600°C时三种配合物的残余质量分别为22%、23%和24%，但数据无法推断残余物的组成或分解机理。

电子吸收光谱与磁性测量方面，由于锌(II)、镉(II)和汞(II)均为d¹⁰闭壳层构型，配合物预期为抗磁性，常规的d-d跃迁或配体场参数分析不适用。漫反射光谱特征经波数-能量换算后，Zn-K呈现0.806-1.488 eV的宽吸收范围，Cd-K在0.515、0.790和1.638 eV处有吸收特征，Hg-K则在0.533、1.500和2.425 eV处有吸收特征。

DFT计算方面，优化后的几何结构显示锌和镉模型呈现对称的双齿配位模式，而汞模型则表现出两个短Hg-O键（2.082 ?）和两个长Hg-O键（2.464 ?）的不对称配位特征。平均M-O键距离按Zn < Cd < Hg的顺序递增。Mulliken和L?wdin布居分析均显示金属原子上存在显著的电荷再分配，偏离形式M²⁺状态。前线轨道能隙方面，锌和镉模型的HOMO-LUMO能隙接近（4.651 eV和4.660 eV），汞模型则较大（4.910 eV），但这些计算值并不对应任何实验漫反射吸收特征。

气相缔合能计算显示，三个模型的缔合能均为较大的负值，顺序为Zn（-673.36 kcal mol^-1）> Hg（-634.99 kcal mol^-1）> Cd（-617.27 kcal mol^-1），与实验热稳定性顺序一致。但需注意，这些是电子缔合能而非实验生成焓，不能直接与热分解数据比较。

局部C=O力常数与频率计算方面，研究人员通过五点能量扫描获得局部C=O力常数，计算频率与实验频率比较显示：锌模型计算值1629.3 cm^-1对实验值1667 cm^-1（误差-37.7 cm^-1），镉模型1630.3 cm^-1对1656 cm^-1（误差-25.7 cm^-1），汞模型1683.6 cm^-1对1647 cm^-1（误差+36.6 cm^-1）。平均绝对误差（MAE）为33.34 cm^-1，均方根误差（RMSE）为33.77 cm^-1，Pearson相关系数r = -0.844。模型虽能复现频率量级，但未能再现实验上Zn→Cd→Hg的频率变化趋势，尤其是汞模型的预测偏差较大，表明中性1:2模型无法代表Hg-K固体的振动环境。

三、关键技术方法

研究人员采用的主要技术方法包括：傅里叶变换红外光谱（FTIR）用于表征配位方式；热重分析（TGA）用于评估热稳定性；漫反射光谱用于电子吸收特征测定；Gouy磁天平用于磁化率测量；元素分析仪与EDTA滴定用于组成测定；密度泛函理论（DFT）计算采用Psi4软件，在PBE0-D3(BJ)/def2-SVP水平下进行几何优化、布居分析、前线轨道分析和局部C=O力常数五点能量扫描。

四、综合讨论与结论

综合讨论部分指出，实验测量与DFT计算关注的是问题的不同层面。实验数据表明三种分离固体产物均不能用1:2化学计量比描述，而计算工作探讨的是中性[M(koj)₂]模型本身的稳定性与特征。频率量级对锌和镉的吻合表明配体氧供体骨架受到扰动，而汞的偏差则说明单一模型无法解释整个系列化合物的行为。硬软酸碱（HSAB）概念和极化率差异有助于解释第12族元素间键合行为的差异，但不能替代组成与结构数据。

研究结论明确指出：对存档材料的重新分析表明分析闭合度不完整，C/M比各不相同，曲酸根与金属之比为1.79、2.88和5.05，因此未证实制备固体具有[M(koj)₂]类型的统一分子式。计算所得T₅₀值较曲酸高约68-71°C。PBE0-D3(BJ)/def2-SVP水平的计算提供了三个中性1:2体系的几何、能量、电荷、键序、轨道能隙和局部C=O扫描数据，但由于原始计算文件不可获取，无法进行独立验证。计算局部模式频率的MAE为33.34 cm^-1，对锌和镉更接近实验值。总体而言，计算为设计的分子模型提供了有用信息，但无法推导实验固体的分子式或结构。产物应继续以Zn-K、Cd-K和Hg-K命名，直至完成全面的组成与结构测定。