Variational Monte Carlo study of soliton excitations in hard-sphere Bose gases