不同尺寸金纳米簇（GNCs）对水分子结构与动力学的影响：分子动力学（MD）模拟研究

《Applied Nano》：Structure and Dynamics of Water Molecules with Gold Nanoclusters as Seen from MD Simulations

　　

研究人员采用分子动力学（MD）模拟研究了分散于不同尺寸金纳米簇（GNCs）之间的水分子的结构性质与动力学。将五种不同直径的球形GNCs置于固定边长的立方模拟盒中，其余体积由足够的水分子填充，以在300 K温度下重现这些Au/水体系的实验密度。MD模拟表明，Au纳米簇尺寸对Au/水体系中水分子的动力学具有显著影响。通过计算Au/水体系中水分子质心的均方位移?Δr²(t)?_i的斜率并采用爱因斯坦关系，获得了这些分子的平均扩散系数D。研究发现，随着Au纳米簇尺寸增大，D值降低；而当Au纳米簇直径增大至1.5 nm时，水的剪切黏度值显著增大。相反，GNCs尺寸对径向分布函数中对应最近邻第一配位壳层的第一峰位置影响较弱。随着距离r增大，第二峰和第三峰的高度降低并消失，这可能表明H₂O分子正在远离羟基。

论文解读：不同尺寸金纳米簇水分散体系中水分子的结构与动力学

一、研究背景与问题
近年来，纳米医学与基于纳米颗粒的递送体系取得显著进展。金纳米簇（gold nanoclusters, GNCs）及其水分散体系在放射治疗、能量产生等方面受到关注，但准确、实时监测有机纳米簇在生物系统中的转归仍面临挑战。检测方法较多依赖荧光，而高分辨率无标记成像常受有机材料光学对比度限制，影响对纳米簇命运的整体认识。金属在水中的分散（metal dispersion in water, MDW）可结合多种成像与分析方法同时检测载体，但将这些复杂工具常规化仍受限。实验研究水含金属分散体系的结构与动力学，可与计算机模拟互补；MD模拟能在全原子分辨率下提供实验难以获得的细节，但其可靠性需要通过与既有实验数据比较来验证。与此同时，水分子氢键具有共价特征，不能仅由库仑和Lennard-Jones相互作用描述，已有水模型（如SPCE、TIP3P、TIP4P、TIP5P）在势参数或几何结构上不同，均不能定量反映水的全部性质。因此，需要选择更合适的水模型，并在统一模型下考察“裸”Au纳米簇尺寸对水结构性质与动力学的影响。

二、研究开展、结论与意义
研究人员开展MD模拟，构建固定边长立方模拟盒中放置近似球形GNCs、其余体积填充水分子的Au/水体系，并在300 K下使水含量对应实验密度。研究比较纯水体系不同水模型的径向分布函数（radial distribution function, RDF）、均方位移（mean square displacement, MSD）、扩散系数D和剪切黏度（shear viscosity），选定TIP3P作为后续Au/水体系的基础水模型。随后考察直径0.6、1.0、1.5、2.0和2.5 nm的球形Au纳米簇（分别含Au₁₃、Au₄₃、Au₈₇、Au₂₄₉、Au₄₅₉）对水结构及动力学的影响。结果表明，水扩散系数随金纳米球尺寸增大而降低；当Au纳米簇直径增大至1.5 nm时，水的剪切黏度明显增大，归因于纳米簇表面积增大使水分子与更大摩擦表面相互作用、流动性降低。GNCs尺寸对RDF中最近邻第一配位壳层第一峰位置影响较弱；随距离r增大，第二、第三峰降低并消失，表明H₂O分子远离OH基团。该研究在分子尺度上给出了Au/水体系中尺寸相关的相对趋势，并提示经典模型与计算资源方面的限制，对理解含金纳米颗粒水体系的结构与动力学具有一定启示。

三、关键技术方法（不超过250字）
采用MD模拟，以GROMACS 5.0.2和OpenMD生成近似球形GNCs；使用优化势液体模拟全原子力场（Optimized Potentials for Liquid Simulations-All Atom, OPLS-AA），金原子采用12–6 Lennard–Jones（LJ）势且电荷设为零；在正则系综（NVT）300 K、v-rescale恒温器、粒子网格Ewald（Particle Mesh Ewald, PME）静电及1 nm截断下模拟。先以SPCE、TIP3P、TIP4P、TIP5P纯水体系的RDF、MSD、D和剪切黏度与实验比较，选定TIP3P。纯水体系平衡20 ns、生产5 ns；Au/水体系平衡200 ns、分析5 ns，并用250 ps短模拟估计黏度；由MSD斜率和爱因斯坦关系求D，由剪切应力张量相关函数求剪切黏度。无样本队列来源。

四、研究结果
2.1 纯水体系模拟
通过计算纯水体系中原子对RDF、MSD及剪切黏度，并比较SPCE、TIP3P、TIP4P和TIP5P模型与实验数据，研究人员发现TIP3P预测的扩散系数与实验值偏差较大，但其黏度最接近实验值，且RDF与实验分布更吻合，因此选择TIP3P作为研究Au纳米簇影响的基础模型。

2.2 Au/水体系模拟
通过构建五种不同直径球形GNCs的Au/水体系并计算RDF、MSD、D和剪切黏度，研究人员发现：GNCs尺寸对RDF第一峰位置影响较弱；随r增大，第二、第三峰高度下降并消失，表明H₂O分子远离OH基团；Au-O第一峰位置随纳米簇尺寸变化几乎不变。MSD斜率随纳米球直径增大而降低，由爱因斯坦关系得到的D随Au纳米簇尺寸增大而降低；最小Au纳米簇使自扩散系数较纯水增加，可能与四面体氢键网络破坏有关。剪切黏度随Au纳米簇直径增大至1.5 nm而显著增大，与纳米簇表面积增大、摩擦表面增大和流动性降低相关。由于模拟盒体积固定，增大纳米簇尺寸同时提高Au体积分数、减少水分子数，最大纳米簇的相邻周期镜像表面间距减小也可能影响输运性质，因此结果主要按相对趋势解释。

五、讨论部分总结
讨论指出，所用模型为Au–水界面的简化经典描述：采用零部分电荷的非极化LJ参数，未显式考虑电子极化、镜像电荷效应以及极小金簇电子结构的尺寸依赖变化；这些限制对亚纳米尺度Au₁₃等可能最显著。因此，结果主要应解释为同一相互作用模型下不同纳米簇尺寸的相对趋势，而非Au–水界面的定量精确描述；更严格处理需极化力场或从头算MD。固定盒子体积下，纳米簇尺寸变化伴随Au体积分数和界面区域相对贡献变化，最大纳米簇的周期镜像表面间距减小也可能影响计算输运性质。研究还指出大规模MD模拟需要显著计算资源和高性能并行计算；通过将含Au₂₄₉和Au₄₅₉两个复杂体系的Green–Kubo积分计算范围扩展至超过50 ps，可为进一步确认结果可靠性奠定基础。

六、研究结论翻译
研究含金纳米簇（GNCs）水的结构性质的一个关键挑战，是分子性质与相组织之间的关系。分子动力学（MD）模拟为这类Au纳米簇/水体系中微观性质与宏观性质关系问题提供了数值途径。研究人员的MD模拟用于研究Au/水体系中水分子的结构性质与动力学。MD模拟结果分析表明，水的扩散系数随金纳米球尺寸增大而降低；这种行为可由纳米球生长伴随的水分子迁移率下降来解释。研究人员的计算还表明，当Au纳米簇直径增大至1.5 nm时，水的黏度明显增大。Au纳米簇直径增大对应于纳米簇表面积增大，因此水分子与更大的摩擦表面相互作用，导致流动性降低。对于最大的纳米簇，相邻周期镜像之间表面间距减小可增强水分子在其间的保持，从而影响计算得到的输运性质。通过将含Au₂₄₉和Au₄₅₉的两个最复杂体系的Green–Kubo积分计算范围扩展至超过50 ps，可为进一步确认结果可靠性奠定基础。模拟还表明，Au纳米簇尺寸对相应函数中对应最近邻第一配位壳层的第一峰位置影响较弱；随着距离r进一步增大，第二峰和第三峰高度下降并消失，表明H₂O分子远离OH基团。将MD模拟整合到溶于水体系的Au纳米颗粒研究中具有显著潜力，但也面临若干挑战。开展大规模MD模拟需要显著计算能力，常需高性能处理资源和图形处理器并行计算以提高效率；缺乏这些资源时，模拟时间更短、体系尺寸更小，限制其覆盖相关过程的能力。研究人员认为，该研究可为含金纳米颗粒水体系的结构与动力学问题提供启示。