This is a simple python/cython code implementing the homogeneous Gutzwiller variational wave function for the Bose-Hubbard model. The search for the optimal wave-function parameters is performed through Simulated Annealing, a Monte Carlo method for stochastic optimization.
This program was used to produce the Gutzwiller phase diagram reported in: Density-dependent hopping for ultracold atoms immersed in a Bose-Einstein-condensate vortex lattice [arXiv:1711.10234 cond-mat.quant-gas], by R. H. Chaviguri, T. Comparin, M. Di Liberto, and M. A. Caracanhas.
New in v1.0.2:
- Adds reference to paper defining the notation for the Bose-Hubbard model.
New in v1.0.1:
- Adds basic tests.
- Adds dependence on future library (for the tests only).