Comparative CHARMM and AMOEBA Simulations of Lanthanide Hydration Energetics and Experimental Aqueous-Solution Structures
The accurate understanding of metal ion hydration in solutions is a prerequisite for predicting stability, reactivity, and solubility. Herein, additive CHARMM force field parameters were developed to enable molecular dynamics simulations of lanthanide (Ln) speciation in water. Quantitatively similar to the much more resource-intensive polarizable AMOEBA potential, the CHARMM simulations reproduce the experimental hydration free energies and correlations in the first shell (Ln-oxygen distance and hydration number). Comparisons of difference pair-distribution functions obtained from the two simulation approaches with those from high-energy X-ray scattering experiments reveal good agreement of first coordination sphere correlations for the Lu3+ ion (CHARMM only), but further improvement to both approaches is required to reproduce the broad, non-Gaussian distribution seen from the La3+ experiment. Second-coordination sphere comparisons demonstrate the importance of explicitly including an anion in the simulation. This work describes the usefulness of less resource-intensive additive potentials in some complex chemical systems such as solution environments where multiple interactions have similar energetics. In addition, 3-dimensional descriptions of the La3+ and Lu3+ coordination geometries are extracted from the CHARMM simulations and generally discussed in terms of potential improvements to solute-structure modeling within solution environments.
- Research Organization:
- Argonne National Lab. (ANL), Argonne, IL (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Basic Energy Sciences (BES)
- DOE Contract Number:
- AC02-06CH11357
- OSTI ID:
- 1461312
- Journal Information:
- Journal of Chemical Theory and Computation, Vol. 14, Issue 3; ISSN 1549-9618
- Publisher:
- American Chemical Society
- Country of Publication:
- United States
- Language:
- English
Similar Records
Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La3+ Interaction with Bioinspired Ligands
CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field