Search

The chemistry of divalent lanthanides, Ln 2+ , is a growing sub-field of heavy element chemistry owing to new synthetic approaches. However, some theoretical aspects of these unusual cations are currently underdeveloped, especially as they relate to their dynamic properties in solution. In this work, we address the hydration of two of the classical Ln 2+ cations, Sm 2+ and Eu 2+ , using atomic multipole optimized energetic for biomolecular applications (AMOEBA) force fields. These cations have not been parameterized to date with AMOEBA, and few studies are available because of their instability with respect to oxidation in aqueous media. Coordination numbers (CN's) of 8.2 and 8.1 respectively for Sm 2+ and Eu 2+ , and 8.8 for both Sm 3+ and Eu 3+ have been obtained and are in good agreement with the few available AIMD and X-ray absorption fine structures studies. The decreased CN of Ln 2+ compared with Ln 3+ arises from progressive water exchange events that indicates the gradual stabilization of 8-coordinate structures with respect to 9-coordinate geometries. Moreover, the effects of the chloride counter anions on the coordination of Ln 2+ cations have been studied at different chloride concentrations in this work. Lastly, water exchange times of Ln 2+ cations have been calculated to provide a comprehensive understanding of the behavior of Eu 2+ and Sm 2+ in aqueous chloride media., (© 2022 Wiley Periodicals LLC.)

Scaling factors for atomic charges derived from the RM1 semiempirical quantum mechanical wavefunction in conjunction with CM1 and CM3 charge models have been optimized by minimizing errors in absolute free energies of hydration, ΔG(hyd) , for a set of 40 molecules. Monte Carlo statistical mechanics simulations and free energy perturbation theory were used to annihilate the solutes in gas and in a box of TIP4P water molecules. Lennard-Jones parameters from the optimized potentials for liquid simulations-all atom (OPLS-AA) force field were utilized for the organic compounds. Optimal charge scaling factors have been determined as 1.11 and 1.14 for the CM1R and CM3R methods, respectively, and the corresponding unsigned average errors in ΔG(hyd) relative to experiment were 2.05 and 1.89 kcal/mol. Computed errors in aniline and two derivatives were particularly large for RM1 and their removal from the data set lowered the overall errors to 1.61 and 1.75 kcal/mol for CM1R and CM3R. Comparisons are made to the AM1 method which yielded total errors in ΔG(hyd) of 1.50 and 1.64 kcal/mol for CM1A*1.14 and CM3A*1.15, respectively. This work is motivated by the need for a highly efficient yet accurate quantum mechanical (QM) method to study condensed-phase and enzymatic chemical reactions via mixed QM and molecular mechanical (QM/MM) simulations. As an initial test, the Menshutkin reaction between NH(3) and CH(3) Cl in water was computed using a RM1/TIP4P-Ew/CM3R procedure and the resultant ΔG(‡) , ΔG(rxn) , and geometries were in reasonable accord with other computational methods; however, some potentially serious shortcomings in RM1 are discussed., (Copyright © 2011 Wiley Periodicals, Inc.)