While implicit solvent models are typically sufficient to describe bulk solvent properties, it has long been appreciated that short-range effects cannot be captured accurately with these continuous models. In an attempt to recover short-range solvent effects in QM calculations, the concept of microsolvation has been used, where solvent molecules from the first solvation shell are represented explicitly in the computation while molecules further away are approximated by continuum models. The major challenge for microsolvation approaches is that the positioning of solvent molecules around a central solute is not trivial, i.e., a possibly multitude of solvent configurations must be approximated by a single placement.