Efficient computation of free energy surfaces of Diels–Alder reactions in explicit solvent at ab initio QM/MM level