Molecular Dynamics Investigation of Alkali Metal Ions in Liquid and Aqueous Ammonia

详细信息    查看全文

• 作者：Esam A. Orabi ; Guillaume Lamoureux
• 刊名：Journal of Chemical Theory and Computation
• 出版年：2013
• 出版时间：May 14, 2013
• 年：2013
• 卷：9
• 期：5
• 页码：2324-2338
• 全文大小：715K
• 年卷期：v.9,no.5(May 14, 2013)
• ISSN：1549-9626

文摘

A polarizable potential model for M⁺鈥揘H₃ interactions (M⁺ = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺) is optimized based on the ab initio properties of the ion鈥揳mmonia dimers calculated at the MP2 level of theory. The optimized model reproduces the ab initio binding energies of M⁺(NH₃)_n (n = 2鈥?) and M⁺(NH₃)_n(H₂O)_m (n, m = 1鈥? and n + m 鈮?4) clusters and gives relative solvation free energies in liquid ammonia in good agreement with experimental data, without further adjustments. It also reproduces binding cooperativity in ion鈥揳mmonia and ion鈥揳mmonia鈥搘ater clusters. The model is used in molecular dynamics simulations of isolated ions in liquid ammonia and in aqueous ammonia solutions with various ammonia molar fractions (0.0 鈮?x_NH3 鈮?1.0). Simulations in liquid ammonia show coordination numbers of 4.0 for Li⁺, 5.3 for Na⁺, 6.1 for K⁺, 6.7 for Rb⁺, and 7.7 for Cs⁺, in very good agreement with available experimental results. Simulations of ions in aqueous ammonia show preferential solvation by water in their first solvation shells and preferential solvation by ammonia in their second shells. Potentials of mean force are calculated between each ion and NH₃ in liquid water, and between each ion and H₂O in liquid ammonia. The results suggest that, in liquid water, Li⁺ and Na⁺ bind NH₃ in their second solvation shells only, while Cs⁺ binds NH₃ in its first solvation shell only (K⁺ and Rb⁺ ions show only weak affinity for NH₃ in water). In liquid ammonia, the ions bind H₂O in their first solvation shells with an affinity following the trend Li⁺ > Na⁺ > K⁺ 鈮?Rb⁺ > Cs⁺.