Multi-resolution approaches like QM/MM or, more recently, ML/MM aim to maintain a very high level of accuracy for the description of the solute while solvent molecules are approximated with classical force fields. These approaches are often fruitful when an explicit description of the electronic structure of the solute is desired, for example to study chemical reactions. However, these methods show the same limitations as classical MD simulations and often require vast computational resources.