[论文解读] Hydration of polyelectrolytes studied by molecular dynamics simulation

Oliver Biermann, E. HAEDICKE+3 more|arXiv (Cornell University)|Jan 9, 2001

Electrostatics and Colloid Interactions参考文献 1被引用 8

一句话总结

本研究采用全原子分子动力学模拟，比较了稀水溶液中羧甲基纤维素钠（CMC）与聚丙烯酸钠（PAA）的溶剂化行为。尽管分子量相近，PAA与水形成的氢键数量（0.083 根/amu）显著多于CMC（0.029–0.036 根/amu），导致其溶剂化作用更强且反离子缔合程度更高；而CMC的构象则高度依赖于取代模式，可形成伸展或致密球状结构。

ABSTRACT

Molecular dynamics simulations of diluted (~2.5 weight percent) aqueous solutions of two polyelectrolytes, namely sodium carboxy methyl cellulose (CMC) and sodium poly(acrylate) (PAA) have been performed. Water and counterions were taken into account explicitly. For CMC the substitution pattern and starting conformation is all-important. Two simulations of CMC oligomers resulted in different structures: One molecule takes a stretched conformation, while the second one keeps a globule-like, toroidal one. PAA is stretched during the whole simulation, with an average characteristic ratio of 8.3. On a local atomistic scale CMC and PAA have different hydrogen-bond properties. The COO- groups of PAA can only act as hydrogen bond acceptors, but due to the high negative charge density there are still more water molecules assembled around PAA than around CMC. There are 0.036 bonds/amu respectively 0.029 bonds/amu to water for the two CMC oligomers, but more than twice as many for PAA: 0.083 bonds/amu. Beside intermolecular hydrogen bonding, there is a significant amount of intramolecular H-bonding for CMC, which is influenced by the COO- groups, which act as strong H-acceptor. In contrast to hydroxy- and carboxylic groups, ether oxygens are hardly involved into hydrogen bonding.

研究动机与目标

• 理解聚合物电解质在水中的原子级溶剂化结构与动力学行为。
• 比较天然聚合物电解质CMC与合成聚合物PAA的溶剂化行为。
• 研究CMC中取代模式对其构象及氢键作用的影响。
• 量化CMC与PAA在氢键作用、反离子缔合及溶剂化壳层结构方面的差异。
• 为未来PAA的粗粒化建模提供原子级数据。

提出的方法

• 对稀水溶液（CMC和PAA浓度约为2.5 wt%）进行了全原子分子动力学模拟。
• 显式模拟了水分子和反离子（Na+），以捕捉离子特异性和溶剂效应。
• 模拟了两种具有不同随机取代模式（DS ≈ 0.7）的CMC低聚物，以评估结构多样性。
• 使用径向分布函数（RDFs）分析聚合物链周围水分子和Na+离子的空间分布。
• 计算氢键统计（单位质量的氢键数）并监测分子内与分子间氢键的动力学行为。
• 通过特征比和随时间变化的构象变化分析链尺寸。

实验结果

研究问题

• RQ1 CMC中的取代模式如何影响其在水溶液中的构象？
• RQ2 CMC与PAA在原子尺度上与水的氢键作用有何差异？
• RQ3 CMC与PAA在电荷密度和离子缔合方面有何不同？
• RQ4 分子内氢键对CMC的溶剂化作用和动力学行为有何影响？
• RQ5 CMC与PAA的溶剂化壳层及Na+的径向分布函数有何差异？

主要发现

• CMC低聚物表现出截然不同的构象：一种保持球状和环状结构，另一种则伸展为伸展构象，具体取决于取代模式。
• PAA在整个模拟过程中保持伸展构象，特征比为8.3，表明其具有高链柔性。
• PAA与水形成的氢键数量为0.083 根/amu，是CMC两种低聚物（0.029–0.036 根/amu）的两倍以上。
• PAA中的COO−基团仅作为氢键受体，但由于电荷密度高，吸引的水分子数量多于CMC。
• CMC II因具有紧凑的球状结构，形成更多分子内氢键，与水形成的分子间氢键更少，相比CMC I。
• 尽管分子量相近，PAA在0.64 nm（屏蔽长度）范围内平均吸引约13个Na+离子，而CMC仅结合0.41个反离子。