Ion-Specific Anomalous Water Diffusion in Aqueous Electrolytes: A Machine-Learned Many-Body Force Field Study with MACE
3M 碘化銫溶液使水擴散加速 23%，MACE 模型成功破解此異常機制。
- MACE 模型克服傳統古典力場限制，成功重現水分子異常擴散現象。
- 碘離子水合層缺乏能量障壁，促使水分子高速交換並帶動擴散加速。
- 鈉離子的 1.2 kcal 能壘強烈束縛水分子，導致溶液整體流動遲滯。
在 3M 濃度的碘化銫（CsI）水溶液中，水分子的擴散速度竟然比純水快上 23%；反觀在同濃度的氯化鈉（NaCl）溶液中，水分子的移動卻減速了 19%。這個看似違反直覺的離子特異性「異常擴散」現象，超過六十年來始終是傳統分子動力學模擬無法準確重現的物理謎題。義大利拉奎拉大學的研究團隊近期透過 MACE 機器學習力場，成功在原子尺度下解開了這個加速與減速的微觀機制。
傳統模型難解的 23% 水分子加速擴散現象
早在半個多世紀前，化學家就提出了「結構形成（kosmotropic）」與「結構破壞（chaotropic）」離子的概念，藉此解釋不同鹽類如何干擾液態水的氫鍵網路。這些概念後來被總結為著名的霍夫邁斯特序列（Hofmeister series），廣泛用於解釋電解質溶液中的熱力學與傳輸現象。在巨觀表現上，最顯著的證據就是水分子擴散係數的濃度依賴性。核磁共振（NMR）實驗明確指出，當鹽濃度低於 3M 時，添加破壞性離子的 CsI 溶液會讓水分子擴散得更快，而添加形成性離子的 NaCl 溶液則會讓水分子變慢。
傳統的古典力場在處理這類離子溶液時，遭遇了系統性的挫敗。幾乎所有非極化力場（例如 Madrid-2019 模型），甚至部分極化力場，都會得出「只要加鹽，水分子擴散就會變慢」的錯誤結論。無論鹽類的化學性質為何，這些模型都無法捕捉到 CsI 溶液中水分子的加速動態。研究人員過去曾指出，必須建立具備動態電荷轉移的明確波動電荷模型，才有可能捕捉到輕微破壞性溶液（如 KCl）中的水分子加速現象，但這依然無法泛用至整個霍夫邁斯特序列。
為了解決這個運算瓶頸，過去曾有研究引入第一性原理分子動力學（AIMD），直接在密度泛函理論（DFT，一種量子力學計算方法）層級處理電子自由度。雖然 AIMD 能夠重現定性的實驗趨勢，並指出離子並未劇烈破壞氫鍵網路，而是改變了動態分佈的尾端特徵，但受限於龐大的運算成本，這類模擬通常只能處理約 100 個原子，且時間尺度僅有數十皮秒（picoseconds），難以取得具統計意義的黏度與擴散係數數據。後續發展的多體勢能函數雖然極度精確，同樣因為算力消耗過高，無法常規應用於高濃度鹽類與長時程動態研究。
涵蓋最高 3.56 m 濃度的水分子擴散率對比
本次研究採用了新一代基於等變圖神經網路的 MACE 框架，建構出專屬的機器學習力場（MLFF）。團隊利用 DFT 計算出的能量、作用力與應力數據來訓練 MACE 模型。這種架構透過高階等變訊息傳遞機制，在保留量子力學級別準確度的同時，將運算成本降低了數個數量級，使得奈秒（nanosecond）等級、包含數千原子的長時間軌跡模擬成為可能。
在常溫常壓的設定下，研究團隊針對濃度範圍介於 0.89 至 3.56 m（mol/kg） 的 NaCl 與 CsI 水溶液進行了長時間的分子動力學模擬。結果顯示，MACE 模型精準重現了實驗中觀察到的異常擴散現象。在 NaCl 溶液中，水分子的相對擴散係數持續低於 1，且隨濃度增加而穩定下降；而在 CsI 溶液中，相對擴散係數則穩定大於 1，明確展現出有別於純水的加速特徵。
相較於其他團隊先前開發的 DeePMD 模型，本次的 MACE 模型在 NaCl 溶液的模擬上展現出顯著的量化精準度提升。這歸功於 MACE 架構能夠更準確地捕捉第一水合層中 Na+ 與水分子之間強烈的交互作用。此外，MACE 模型也準確預測了剪切黏度（shear viscosity）的變化趨勢，對於 CsI 給出了低於純水的相對黏度，對於 NaCl 則給出了高於純水的預測，全面克服了傳統非極化力場的系統性偏差。
半徑 3.2 埃與 4.3 埃揭示的水合層結構差異
要釐清異常擴散的成因，必須從離子周圍的水合層（hydration shell）結構著手。研究團隊透過徑向分佈函數（RDFs）分析發現，在 NaCl 溶液中，鈉離子（Na+）與氯離子（Cl-）與周圍水分子建立了高度結構化的關聯性。特別是 Na+，其第一水合層呈現出極高的峰值與極深的低谷，代表水分子被緊緊束縛在離子周圍，形成延伸至 1 奈米的多層穩定結構。
對比之下，CsI 溶液展現出截然不同的幾何分布。銫離子（Cs+）與碘離子（I-）的徑向分佈函數僅有較低的第一峰值與寬闊平淺的低谷，這意味著它們的周圍缺乏定義明確的邊界。這種鬆散、彌漫的水合層結構，讓水分子能夠保有高度的排列彈性。儘管兩者的結構差異顯著，但整體水網路的氫鍵數量分析顯示，隨著鹽濃度增加，兩種溶液的平均氫鍵數都呈現相似的下降趨勢。這證明單純觀察「氫鍵數量」這個指標，不足以解釋兩者截然不同的動態行為。
在精確界定水合層的影響範圍上，測量團隊計算了離子與氧原子配位數的反曲點，將 Na+ 的第一水合層截斷半徑設定為 3.2 埃（Å），而 I- 的截斷半徑則延伸至 4.3 埃（Å）。透過將水分子依照其所處的特定水合層進行分類追蹤，研究人員進一步將第一水合層分解為「專屬單一離子」的純水分子，以及「被多個離子共享」的重疊水分子，藉此分離出個別離子對水分子動態的具體貢獻。
由 1.2 kcal 能量障壁驅動的微觀加減速機制
動力學分析揭示了關鍵的微觀加減速機制。在 NaCl 溶液中，專屬於 Na+ 第一水合層內的純水分子，在所有測試濃度下都呈現出顯著的擴散遲緩；甚至在較遠的第二水合層中，Na+ 依然發揮著可觀的遲滯效應。Cl- 雖然也有減速作用，但在低濃度時其周圍水分子的行為仍相當接近體相純水（bulk water），這證實了溶液整體的黏滯度攀升主要歸咎於 Na+ 強大的束縛力。
情境切換到 CsI 溶液，推動水分子加速的引擎主要來自 I- 離子。專屬於 I- 第一水合層的水分子展現出比體相純水更活躍的移動能力。這項觀察與實驗上的黏度 B 係數完美吻合：I- 具有最大的負值（-0.068），象徵其具備最強的流動促進能力；而 Na+ 則具有極大的正值（0.086），象徵著強烈的遲滯效應。
計算團隊進一步繪製了離子與氧原子的平均力勢能（PMF）地貌圖，為這兩個對立現象提供了熱力學上的終極解答。在低濃度下，氧原子與 Na+ 之間存在高達 1.2 kcal/mol 的能量障壁，阻擋了第一與第二水合層之間的水分子交換。這道障壁遠高於純水內部的 0.6 kcal/mol，導致水分子被死死「鎖定」在 Na+ 旁。
相反地，氧原子與 I- 或 Cs+ 之間的勢能圖上完全找不到明確的障壁，只存在一個寬廣的淺谷；水分子跨越約 0.6 kcal/mol 的微小爬升後，就能進入平坦的自由能地貌，輕易與外部體相水進行高速交換。這些數據勾勒出一幅連貫的物理圖像：異常擴散的本質並非整個水網路被徹底摧毀，而是由第一及第二水合層內的局部微觀動態所主導。結構鬆散的 I- 成為了促進水分子的交換中樞，而緊密咬合的 Na+ 則扮演了拖慢整體流動的煞車皮。
機器學習力場填補了量子力學與長時間動態之間的鴻溝，首度在原子層級證實，離子特異性的擴散異常取決於微觀的自由能地貌。