Association Free Energies of Metal Cations with Mesylate and Acetate in Brine Calculated via Molecular Dynamics Simulation

TL;DR

This study uses molecular dynamics simulations to calculate and compare the binding free energies of Na+ and Ca2+ ions with sulfonate and acetate groups in brine, revealing differences in binding strength and dynamics.

Contribution

It provides detailed atomistic insights into ion pairing free energies of metal cations with functional groups in brine, including new calculations of binding energies and comparison with experimental data.

Findings

Na+ and Ca2+ bind weakly to sulfonate with contact times of 300-400 ps.

Ca2+ binds strongly to acetate with correlation times of hundreds of nanoseconds.

Binding free energies vary with force field parameters and water models, aligning with experimental data.

Abstract

Here we report results from fully-atomistic molecular dynamics simulations and present the calculated binding free energies of the associations between metal cations and anionic functional groups under brine condition. Specifically the ion pair formations of Na+ and Ca2+ metal cations with sulfonate and carboxylate functional groups were considered. In our simulations, sulfonate and carboxylate functional groups in polyelectrolytes were modeled by mesylate and acetate anions, respectively. The simulations show that the binding strengths of Na+ and Ca2+ with sulfonate group are relatively weak with correlation times for the contact ion pairs being around 300-400 ps. In the case of acetate, the binding strength with Na+ was found to be comparable to the case of sulfonate group, while the binding of acetate with Ca2+ was found to be much stronger with a correlation time of the order of hundreds of nanoseconds. The binding free energies of Na+ and Ca2+ with sulfonate as well as that of Na+ with acetate have been calculated directly from the radial distribution functions. For the case of binding between Ca2+ and acetate constrained MD simulations with umbrella sampling were carried out to improve the sampling of the phase space and the potential of mean force was obtained by the weighted histogram analysis method. The binding free energy of a 2:1 binding event between one Ca2+ ion and two acetates was also calculated. Finally the binding free energies between Ca2+ and acetate were compared for different force field parameters and water models, and also with experimental and calculated values reported in the literature.