Combining the Fluctuating Charge Method, Non-Periodic Boundary Conditions and Meta-Dynamics: Aqua Ions as case...

We present the current status of development of our code for performing Molecular Dynamics (MD) simulations exploiting a polarizable force field based on the Fluctuating Charge (FQ) method and non-Periodic Boundary Conditions (NPBC). Continuing on the path set in a previous work, we increased the capabilities of the code by implementing a number of new features, including: a non-iterative algorithm for rigid trigonal molecule simulations; two additional temperature coupling schemes; a meta-dynamics based approach for effective free energy evaluations. Although these are well known algorithms, each present in one or more widely used MD packages, they have now been tested, for the first time, in the context of the FQ model coupled with NPBC. As case studies, we considered three aqueous ions of increasing charge, namely Na⁺, Ca²⁺ and La³⁺, at infinite dilution. In particular, by exploiting a computational approach recently proposed by our group and based on the metadynamics technique, we focused on the important role played by solvent polarization on ionic hydration structures, also investigating the free energy landscapes of ion coordination and the water exchange rates. Such an approach, previously tested with standard non-polarizable models, was applied here to evaluate the effects of explicit polarization on water exchange barriers between different solvent coordination structures. Moreover, we have analyzed and discussed in some detail non-linear electrostatic effects arising from solvent polarization while going from a mono- to a di- and trivalent ion.