K-feldspar中Si/Al-order是否影响其水化层结构？——orthoclase与microcline–水界面的2D/3D AFM与MD研究

《Physical Chemistry Chemical Physics》：Does the Si/Al-order in K-feldspars influence their hydration structure?Open Access

　　

云特性影响天气和气候，例如反照率或降水可能性，强烈取决于云滴是液态还是冻结。在冰核粒子（ice nucleating particles, INP）存在下的异质冰核化在冰云形成中占主导。INP的重要例子是K-feldspar矿物尘粒。K-feldspar一词描述化学组成相同（KSi₃AlO₈）但晶体结构、Si/Al-order和（表面）微观结构不同的一组矿物。发现具有高Si/Al-order的K-feldspar（microcline）比具有较低Si/Al-order的K-feldspar（orthoclase和sanidine）具有更高冰核效率。然而，这种冰核效率差异的起源尚未阐明。在此，研究人员展示orthoclase–水和microcline–水界面水化层结构的实验与理论比较。水化层结构强烈影响浸没冻结中的异质冰核化机制。在orthoclase (001)–和(010)–水界面的二维和三维原子力显微镜（atomic force microscopy, AFM）实验以及分子动力学（molecular dynamics, MD）模拟均揭示出距表面0.5 nm以上的类似棋盘状结构，水层间距为(0.27 ± 0.03) nm。研究人员在orthoclase和microcline上未观察到Si/Al-order对水化结构的显著影响，但模拟表明水化层动力学发生变化。

研究背景：云滴为液态或冻结会显著影响云反照率和降水，纯水均质冰核化约发生在?38 °C，而冰核粒子（ice nucleating particles, INP）可在更高温度触发异质冰核化（heterogeneous ice nucleation）。K-feldspar（钾长石）是高效矿物INP。长石为架状硅酸盐，(Si,Al)O₄四面体共氧连接，框架负电荷由K⁺、Na⁺、Ca²⁺补偿，并易沿(001)和(010)解理。Atkinson et al. 2013将长石尘识别为混合相云主要INP；碱长石通常比斜长石更活跃，K-feldspar成核温度尤其高，但不同样品活性差异及机制仍不清楚。显微与液滴冻结实验将IN活性与微观表面形貌联系；perthitic microstructure（条纹长石显微结构，富Na和富K域出溶）样品更活跃，出溶产生晶格应变并可形成裂纹，冰核化常始于裂纹或孔等特征。电子显微显示六方冰（hexagonal ice, I_h）在不同长石晶面同向生长，提示外延生长，解释为冰I_h棱柱面与长石高能(100)面晶格匹配，仅见于台阶或裂纹。Si/Al-order（硅/铝有序度）也影响活性：有序microcline（微斜长石）高于无序sanidine（透长石）；原子力显微镜（atomic force microscopy, AFM）沉积冰核化显示microcline台阶和裸无缺陷(001)台面均可成核，sanidine主要沿台阶，提示Si/Al-order可能影响浸没冻结（immersion freezing），而浸没冻结原子尺度机制由界面平衡水化层结构决定。问题在于化学组成、Si/Al-order和表面形貌内在关联，perthitic microstructure与Si/Al-order难分离；同时缺乏原子尺度表面结构实验；环境长石表面接触水膜，故长石–水界面研究重要。Dickbreder and Sabath et al.用2D/3D AFM研究microcline(001)–水界面，发现0.5 nm内两层水及横向周期结构；Fenter et al.用X射线反射率（X-ray reflectivity, XRR）发现orthoclase(001)两层水间距(0.26 ± 0.03) nm，(010)仅一层，但缺横向信息；Nakouzi et al.和Umeda et al.用AFM与分子动力学（molecular dynamics, MD）研究orthoclase和albite界面。理论多聚焦microcline，但未解释其优异性，Kiselev et al.、MD和DFT结果不一。Si/Al-order对界面水结构和冰核活性的影响仍缺实验与理论比较。因此，研究人员开展orthoclase–和microcline–水界面的深入实验与理论比较：室温AFM，300 K MD；结论为orthoclase(001)/(010)均有两层棋盘状水化结构，间距(0.27 ± 0.03) nm；与microcline实验无显著差异，MD显示动力学变化；意义在于分离Si/Al-order影响，提示动力学异质性可能贡献冰核活性差异，下一步应关注表面形貌。

方法上，研究人员采用二维/三维AFM在固–液界面成像orthoclase (001)/(010)–水界面，并用unDrift校正漂移；以X射线衍射（X-ray diffraction, XRD）和化学分析确定单晶组成及Al在T1/T2位占据；以MD在300 K模拟orthoclase与microcline (001)/(010)–水界面，采用LAMMPS、CLAYFF力场、SPC/Fw水模型、PPPM长程静电和150 ns轨迹（前50 ns平衡）。orthoclase单晶购自SurfaceNet GmbH（德国），来自马达加斯加同一地区；microcline样品来自Dickbreder and Sabath et al.前期研究；化学分析由Mikroanalytisches Laboratorium Kolbe（德国）完成。

2.1 Hydration layer structure on K-feldspar (001)：样品表征显示两个orthoclase单晶组成和Si/Al-order相似，XRD确定Al在T1/T2位占据在orthoclase范围；所用microcline更有序（三斜，约77% Al占T1(o)），组成相近，故水化差异可归因于Si/Al-order。MD显示orthoclase和microcline(001)水分子配位K⁺、硅醇基（silanol）和铝醇基（aluminol）；microcline的12个aluminol完全有序占T1(o)，orthoclase的8个aluminol占5个T1(o)和3个T2(o)，破坏对称；1D密度均见第一水化层（first hydration layer, HL1）和第二水化层（second hydration layer, HL2）两峰，0.5–1 nm达体相；水氧/氢密度相似；2D密度显示HL1受表面K⁺、silanol、aluminol模板化，HL2具表面晶胞周期；Al无序仅扰动HL1，HL2影响小。2D AFM见点状周期结构，晶胞0.84 × 1.35 nm²和89.9°，与XRD一致；模拟HL2与实验点状结构吻合，说明成像HL2。3D AFM见约0.6和0.4 nm棋盘状振荡；基于溶剂探针近似（solvent tip approximation, STA），亮区为高水密度；位置A/B的Δν_exc vs. z_p+q_s剖面峰间距(0.27 ± 0.03) nm，与XRR和MD一致；两剖面峰间距(0.12 ± 0.04) nm，略大于microcline(001)的0.07 nm。模拟垂直切片和局部水密度剖面也见棋盘状，每位置两峰，平均间距0.27 nm。结论：orthoclase与microcline(001)水化结构相似；HL2有序无明显差异，orthoclase的Al无序仅轻微扰动HL1。

2.2 Hydration layer structure on K-feldspar (010)：MD显示orthoclase/microcline(010)需移除面内一半阳离子以保持电中性；microcline(010) α终止的24个羟基均与Si键合；orthoclase(010)有5个aluminol和19个silanol，破坏对称；1D水氧密度峰多于(001)，第一小峰与表面K⁺重合，随后HL1、HL2和第三小峰，0.7–1 nm达体相；orthoclase/microcline(010)水氧/氢密度相似；2D图显示HL1沿[001]有序，HL2横向有序弱于(001)；orthoclase(010)的HL1密度峰更明显，尤其aluminol周围；未占据K⁺位与邻近羟基HL1密度最大相关、HL2密度较低。2D AFM见沿[001]条纹，晶胞0.71 × 0.83 nm²和66°，与XRD一致；模拟HL2弱可见。3D AFM见约0.5 nm棋盘状振荡；沿[001]每剖面两峰，沿[100]位置A三峰、B一峰；A峰与B谷重叠；同位置峰间距(0.27 ± 0.03) nm，不同位置(0.12 ± 0.01) nm。模拟垂直切片和局部密度剖面一致（沿[001]每位置两峰，沿[100] A三峰B两峰，同剖面平均0.27 nm）。比较(001)/(010)：2D实验和模拟分别为点状与条纹；3D和MD垂直结构相似，均为0.5 nm棋盘状、两位置，同位置峰距水分子大小，不同位置半距。与microcline比较：实验精度内无差异；orthoclase–水界面更模糊；模拟显示orthoclase(001)/(010)第一水化层与microcline有定量差异，由Al替代无序造成，Si/Al-order对HL2影响很小。需研究Al无序对水化层动力学影响。

2.3 Hydration layer dynamics from simulation：MD显示四个表面K⁺在表面与溶液间慢交换（数十ns），需100 ns轨迹、前50 ns平衡；部分K⁺在数十ns于平衡位间移动但不完全入液；部分溶剂化阳离子在(001)可占表面与HL1间亚稳态位，不在(010)。径向分布函数和配位数显示orthoclase/microcline相似晶面配位数相近，除microcline(010)无aluminol；aluminol配位数高于silanol且不同长石面相似；K⁺在(010)第一溶剂化壳比(001)多约一个水分子。水驻留时间最短在K⁺（10¹–10² ps），最长在aluminol（10²–10³ ps），silanol为10¹–10³ ps；orthoclase(001)的silanol/aluminol对称性破坏，orthoclase(010)含aluminol而microcline(010)无，使orthoclase化学相似位点驻留时间分布更宽；microcline和orthoclase(010)最多半数K⁺位占据，邻近K⁺位全程占据的羟基驻留最短，水分子取代邻近空K⁺位者驻留最长。检查四个不同Al分布orthoclase slab，发现表面阳离子动力学和局部水化结构有定量差异，主要归因局部aluminol密度，较小程度归因表面层净电荷。结论：orthoclase的Al无序相比microcline使化学相似位点水驻留时间分布变宽；第一水化层动力学变化不影响HL2结构；驻留时间跨度可超三个数量级但仍快于AFM针尖振荡，故AFM成像平均水密度。

讨论与结论总结：化学组成、Si/Al-order和表面组成是长石冰核效率的潜在因素；浸没冻结机制依赖界面水化层结构。AFM未见不同Si/Al-order诱导的显著水化结构差异；模拟显示aluminol分布对平衡水密度分布有较小定量影响，HL2在microcline和orthoclase相似；但orthoclase第一水化层水交换动力学更异质。较慢水动力学是过冷水中有序冰核出现的前兆；若“动力学异质性”在中等过冷度持续，可能贡献二者冰核活性差异。因化学组成和Si/Al-order对平坦低能解理面水化结构无显著影响，下一步应研究特定表面形貌（如perthitic microstructure）对水化层结构和动力学的影响；Si/Al-order仍可能通过复杂表面特征或间接影响微结构、出溶、缺陷群而影响冰核活性。

研究结论翻译：化学组成、Si/Al-order和表面组成已被识别为潜在影响长石冰核效率的因素。在浸没冻结中，冰核机制强烈依赖界面水化层结构。研究人员的AFM实验未观察到orthoclase和microcline表面不同Si/Al-order诱导的显著水化结构差异，但模拟揭示不同aluminol分布对界面平衡水密度分布有相对较小的定量影响。特别是第二水化层结构——可能在二维AFM中成像——在microcline和orthoclase表面相当相似。然而，第一水化层中的水交换动力学在orthoclase上比microcline更异质。较慢的水动力学已被识别为过冷水中有序冰核出现的前兆。如果这种“动力学异质性”在中等过冷度下持续，它可能贡献microcline和orthoclase之间冰核活性的差异。由于化学组成和Si/Al-order似乎对orthoclase和microcline平坦、低能解理面的水化结构无显著影响，下一个问题是特定表面形貌——例如在perthitic microstructure背景下——对水化层结构和动力学有何影响。Si/Al-order仍可能影响更复杂表面特征上的冰核，或通过间接效应，例如影响微观结构、出溶或缺陷群，这些可能在决定冰核活性中起重要作用。