Boson-Gutzwiller Ground state configuration of the Boson-Gutzwiller wavefunction using imaginary time evolution (gradient-descent) method.