“…Note that FMO-MD is of a fully QM scheme where the energy and forces are evaluated at each step. A main purpose of the present study is to investigate the hydration structures related to the gadolinium break [8][9][10][11][12] as well as the lanthanide contraction [1,3], where the FMO-MD simulations were actually applied to t La(III), Eu(III), Gd(III), Tb(III) and Lu(III) ions in a droplet of 64 water molecules. The presence of 4f electrons was implicitly treated to reduce computational times, by using the model core potential (MCP) [37].…”