We present a novel formulation based on quantum Monte Carlo techniques for the treatment of volume polarization due to quantum mechanical penetration of the solute charge density in the solvent domain. The method allows to accurately solve Poisson's equation of the solvation model coupled with the Schrodinger equation for the solute. We demonstrate the performance of the approach on a representative set of solutes in water solvent and give a detailed analysis of the dependence of the volume polarization on the solute cavity and the treatment of electron correlation.