层状羟基碳酸盐KAl(CO3)(OH)2的负线性压缩性及结构起源：力学与光学响应的第一性原理研究

《Materials Advances》：Anomalous mechanical and optical response of layered dawsonite-type KAl(CO3)(OH)2: negative linear compressibility and its structural originOpen Access

　　

层状羟基碳酸盐作为CO₂封存、催化及耐压应用的功能材料正受到日益关注，但其在压缩条件下的力学和光学响应仍认识不足。研究人员报告在钾道森石KAl(CO₃)(OH)₂——一种层状羟基碳酸盐——中发现负线性压缩性（negative linear compressibility, NLC），采用第一性原理密度泛函理论（density functional theory, DFT）计算，压力范围为0–10 GPa。NLC在机械稳定区（0–5 GPa）内稳健，该区间内所有Born稳定性判据均满足。约5 GPa以上结构失去机械稳定性（C₅₅和C₆₆变为负值），7–10 GPa数据仅为完整性而报告，因其对应于假设的、机械不稳定的状态。在静水压缩下，晶胞体积减小，而晶格b参数膨胀，从而产生沿b轴的NLC。结合原子间距离分析与晶体轨道重叠布居/晶体轨道哈密顿布居（crystal orbital overlap/Hamilton populations, COOP/COHP）分析表明，异常变形源于O–H?O氢键的各向异性增强以及碳酸根离子（CO₃^2?）的刚性。晶格动力学计算显示，具有负Grüneisen参数的模式定域于O–H伸缩振动，其氢键与b轴成一定角度，从而建立压力诱导氢键增强与异常压力响应之间的关联。弹性常数揭示强各向异性及约5 GPa以上的机械不稳定性，而计算得到的红外（infrared, IR）和拉曼（Raman）光谱、纵光学–横光学（longitudinal optical–transverse optical, LO–TO）分裂及介电函数与道森石族已报道的实验范围一致，并揭示与NLC相关的压力诱导双折射（birefringence）增加。这些结果建立了层状羟基碳酸盐中化学键、晶格动力学与异常力学和光学行为之间的直接联系，并为设计具有定制压力响应性质的材料开辟了前景。

研究背景：道森石族化合物通式为MAl(CO₃)(OH)₂（M = Na⁺、K⁺或NH₄⁺），其结构由沿c轴延伸的[Al₂(OH)₄(CO₃)₂]^2?碳铝酸盐链构成，链间通过O–H?O氢键和K⁺等阳离子相互作用连接。该类材料因可固定CO₂、用作催化前驱体、阻燃添加剂及耐压功能材料而受到关注。然而，道森石族在压力下的力学和光学响应研究十分有限；Na-道森石仅有少量压缩和热膨胀研究，K-道森石的系统压力性质缺失，其负线性压缩性（NLC）此前未见报道。因此，研究人员有必要在宽压力范围内阐明KAl(CO₃)(OH)₂的结构、化学键、弹性、振动和介电行为，并揭示NLC的结构起源。

研究人员开展了0–10 GPa范围内的第一性原理研究，发现K-道森石沿b轴具有NLC，且该效应存在于机械稳定窗口0–5 GPa内；约5 GPa以上C₅₅和C₆₆变为负值，表明机械稳定性丧失，故7–10 GPa数据仅作为趋势报告。NLC来源于碳铝酸盐链的特定架构、CO₃^2?基团的刚性以及压力下O–H?O氢键的方向性增强；沿a、c轴压缩伴随链沿b轴分离，形成铰链式变形。该研究将化学键、晶格动力学、弹性各向异性与压力诱导双折射联系起来，说明NLC并非孤立力学异常，而是具有可测光学特征的结构–动力学现象。其意义在于深化层状羟基碳酸盐结构–性质关系认识，并为CO₂封存、传感及光机械应用中的压力响应材料设计提供依据。

关键技术方法（不超过250字）：研究人员基于CRYSTAL软件实施全电子DFT计算，采用双ζ价极化基组（double-zeta valence with polarization, DZVP），以PBEsol为主，并与PBE、PBE-D3(BJ)、B3LYP比较；k点取Monkhorst–Pack 4×4×4，固定晶胞体积并弛豫晶格和原子坐标模拟静水压缩，用BM2/BM3/BM4状态方程（equation of state, EoS）拟合；以有限变形法求弹性常数并检验Born判据，经Voigt–Reuss–Hill平均求多晶模量；计算Γ点声子、LO–TO分裂、IR/Raman强度、介电张量、COOP/COHP与Mulliken布居及模式Grüneisen参数。研究未涉及实验样本队列，均为计算模型。

研究结果：
3. 晶体结构与化学键：通过结构优化、形变密度、Mulliken布居和COOP/COHP分析，确认KAl(CO₃)(OH)₂属Cmcm空间群，由[Al₂(OH)₄(CO₃)₂]^2?链、K⁺通道和O–H?O₃氢键构成；化学键呈层级特征，包括CO₃^2?和O–H的强共价键、Al–O的离子–共价混合键、K–O的纯离子键及中等强度氢键。PBEsol对晶胞体积和结构参数描述最接近实验，PBE-D3(BJ)表明色散作用较小。

4. 状态方程及压力对结构和化学键的影响：通过固定体积静水压缩和Birch–Murnaghan状态方程拟合，PBEsol给出B₀约47 GPa；随压力增加晶胞体积和a、c参数下降，而b参数由11.8661 ?增至11.9584 ?（0–10 GPa，7–10 GPa仅为不稳定趋势），产生沿b轴NLC。H?O₃距离在10 GPa内缩短8.6%，O₁–H?O₃角由165.8°增至168.3°（7 GPa），ICOOP/ICOHP显示氢键共价成分增强；Al–O₁削弱、Al–O₂增强，C–O键几乎不变，带隙由5.35 eV略降至5.29 eV，说明电子结构重排与异常变形一致。

5. 压力下弹性和力学性质：通过弹性常数、Born判据和Voigt–Reuss–Hill平均，发现0–5 GPa内机械稳定，7 GPa时C₅₅和C₆₆为负；弹性响应强各向异性，C₁₂ > C₁₁导致S₁₂负值近抵消S₂₂，使β_b约?1.66 TPa^?1；B_a = 54.6 GPa、B_c = 213.7 GPa而B_b为负，E_a、E_b、E_c分别约30.05、67.84、164.09 GPa；多晶B_H随压力增加，而G_H、E_H、维氏硬度（Vickers hardness, H_V）下降，弹性各向异性参数增大，声速和德拜温度（Debye temperature, Θ_D）下降，最小热导率（minimum thermal conductivity, k_min）略降；C₅₅在零压仅3.51 GPa，提示xz剪切模式接近软不稳定。

6. 压力下振动和介电性质：
6.1. 振动谱的对称性和一般特征：通过Γ点声子计算和对称性分解，20原子原胞给出60个模式；Al不参与Raman活性模式，K主要贡献低频晶格振动，高频模定域于CO₃^2?和OH；[CO₃]与[OH]子晶格计算支持谱峰归属，并显示CO₃^2?对称性降低导致简并解除和ν₁红外激活。

6.2. 红外光谱：通过解析IR强度，计算得ν₃约1442/1527 cm^?1、ν₁约1096 cm^?1、ν₂约859 cm^?1、ν₄约744/726 cm^?1，ν_OH约3496/3524 cm^?1；结果与道森石族实验范围一致，ν_OH高度定域，反映中等强度且不等价氢键。

6.3. 拉曼光谱：通过解析Raman强度，最强线为CO₃^2?的ν₁（1096 cm^?1）和ν₃（1499 cm^?1），ν_OH的Ag模约3527 cm^?1、B_1g模约3487 cm^?1；低频区包含晶格模及[CO₃]、[OH] librational振动。

6.4. LO–TO分裂和反射光谱：通过非解析修正和Lyddane–Sachs–Teller关系，发现1527、1442、491 cm^?1等模具有较大LO–TO分裂，而744、635 cm^?1模分裂很小；B_3u模的LST乘积2.04与ε_0,xx/ε_∞,xx = 1.98在3%内一致，其他偏振因指认不唯一未作定量检验；反射谱在1000–1500 cm^?1显示由ν₃、ν₂主导的Reststrahlen带和各向异性。

6.5. 压力对振动谱的影响：通过频率–压力拟合模式Grüneisen参数，C–O内模γ约0.1–0.2，含K和Al–O变形的晶格模γ约0.2–2.6；O–H伸缩振动具有负γ（约?0.19至?0.25），是压力下氢键增强的谱学指纹，但其与NLC为相关而非直接因果；NLC因果链为H?O₃缩短、刚性CO₃^2?与OH铰链式旋转、沿b轴各向异性膨胀。作者指出同类铰链机制可能伴随沿b轴负热膨胀，但完整准谐分析未在本工作完成。

6.6. 介电性质：通过Drude–Lorentz模型和CPHF方法计算介电张量，ε_∞主值在零压为1.996、2.433、2.288（x、y、z），随压力略增；ε_0,xx随压力显著下降，ε_0,yy和ε_0,zz变化较小；双折射Δn = n_y ? n_x由0 GPa的0.147增至7 GPa的0.159，dΔn/dP约0.0017 GPa^?1，与b轴伸长及NLC相关。

6.7. 热力学函数：通过声子热力学计算，零点能E₀由265.6增至274.2 kJ mol^?1（7 GPa），而热能量E_T、熵S和定容热容C_V随压力下降；低频晶格模对压力敏感，其“冻结”解释熵和热容降低，正γ模主导总体趋势，负γ模仅部分补偿。

讨论部分总结：原文未单设讨论部分，结论段将NLC归纳为层状碳铝酸盐框架中由氢键各向异性增强、CO₃^2?刚性和铰链式旋转共同导致的结构–动力学现象；弹性各向异性、剪切模量和硬度随压力下降、OH伸缩模负Grüneisen参数及双折射增加分别构成其力学、振动和光学关联。作者强调预测NLC仅存在于0–5 GPa机械稳定窗口，7–10 GPa数据不应解释为物理可实现状态；PBEsol因最接近实验结构而作为分析基础。该工作为层状羟基碳酸盐结构–性质关系及压力响应材料设计提供依据，并指出需实验验证NLC及进一步研究更高压力下声子和弹性异常。

研究结论部分翻译：研究人员在0–10 GPa压力范围内对钾道森石KAl(CO₃)(OH)₂的晶体结构、化学键、弹性、振动和介电性质进行了系统第一性原理研究。需要强调，预测的NLC存在于机械稳定压力窗口（0–5 GPa）内；超出该范围，弹性常数C₅₅和C₆₆变为负值，表明机械稳定性丧失，因此7–10 GPa结果仅作为趋势给出，不应解释为晶体物理可实现状态。采用PBEsol、PBE-D3和杂化B3LYP等多种交换–相关泛函评估了预测可靠性并证明计算模型选择的合理性；PBEsol最能重现实验结构参数，因而构成本分析基础。K-道森石中的化学键呈复杂层级：CO₃^2?和羟基内为强共价键，Al–O为中间离子–共价键，K–O为纯离子相互作用，H?O为中等强度氢键；这种组合决定结构显著层状特征及其高度各向异性压力响应。本工作的中心结果是发现沿晶格b轴的NLC：静水压缩下晶胞体积减小而b参数增加，这是道森石族中首次报道的NLC。该异常源于碳铝酸盐链的特定架构、CO₃^2?基团刚性及压力下氢键方向性增强；沿a、c轴压缩伴随链沿b轴分离，形成异常变形模式。弹性常数和多晶模量显示强各向异性：沿c轴单轴压缩模量大于沿a轴的三倍以上，且杨氏模量E_b随压力增加快于其他方向。机械稳定性保持至约5 GPa，超过后出现不稳定迹象，与声子谱演化一致；维氏硬度随压力下降，反映压力诱导化学键重排。振动光谱包括IR和Raman强度可在孤立[CO₃]和[OH]子晶格模型内较好描述。压力下多数模式频率增加，但OH伸缩振动和若干框架变形模表现负Grüneisen参数；这些模式定域于OH伸缩振动并反映压力诱导氢键增强，其负Grüneisen参数为NLC铰链机制所依赖的氢键重排提供谱学指纹。基于Lyddane–Sachs–Teller框架分析LO–TO分裂揭示压力下介电响应异常重分布，包括与异常晶格动力学相关模式的贡献；电子介电张量逐渐更具各向异性，导致伴随b轴压力诱导伸长的双折射增加。综合结构、弹性、振动和介电数据表明，K-道森石中的NLC是由氢键各向异性增强和碳铝酸盐框架特定架构引起的结构–动力学现象。这些结果深化了对层状羟基碳酸盐结构–性质关系的理解，并为设计具有定制异常力学和光学响应的功能材料提供设计原则。研究结果直接关系到道森石在地质条件下（包括CO₂封存）的行为，以及用于传感和光机械应用的压力响应材料开发。进一步研究的一个有前景方向是对预测NLC进行实验验证，并详细研究更高压力下的声子和弹性异常。