The application of CMD to organometallic complexes containing lanthanide-series atoms is not popular because of the difficulty of describing the interactions with explicit force-field functions. A modeling and development protocol for the force field of Gd 3+ complexes in aqueous solution was recently published for relatively small ligands, and the water dynamics of the hydration shell around a Gd 3+ complex in a solvent was successfully simulated.