Abstract Motivated by recent advances in the representation of ground state wavefunctions of quantum many-body systems using restricted Boltzmann machines as variational ansatz, we utilize an open-source platform for constructing such ansatz called NetKet to explore the extent of applicability of restricted Boltzmann machines to bosonic lattice models. Within NetKet, we design and train these machines for the one-dimensional Bose-Hubbard model through a Monte Carlo sampling of the Fock space. We vary parameters such as the strength of the onsite repulsion, the chemical potential, the system size and the maximum site occupancy and use converged equations of state to identify phase boundaries between the Mott insulating and superfluid phases. We compare the average density and the energy to results from exact diagonalization and map out the ground state phase diagram, which agrees qualitatively with previous finding obtained through conventional means.
Recently, the use of neural quantum states for describing the ground state of many- and few-body problems has been gaining popularity because of their high expressivity and ability to handle intractably large Hilbert spaces. In particular, methods based on variational Monte Carlo have proven to be successful in describing the physics of bosonic systems such as the Bose-Hubbard (BH) model. However, this technique has not been systematically tested on the parameter space of the BH model, particularly at the boundary between the Mott insulator and superfluid phases. In this work, we evaluate the capabilities of variational Monte Carlo with a trial wavefunction given by a Restricted Boltzmann Machine to reproduce the quantum ground state of the BH model on several points of its parameter space. To benchmark the technique, we compare its results to the ground state found through exact diagonalization for small one-dimensional chains. In general, we find that the learned ground state correctly estimates many observables, reproducing to a high degree the phase diagram for the first Mott lobe and part of the second one. However, we find that the technique is challenged whenever the system transitions between excitation manifolds, as the ground state is not learned correctly at these boundaries. We improve the quality of the results produced by the technique by proposing a method to discard noisy probabilities learned in the ground state.
We present a deep neural network (DNN) -based model (HubbardNet) to variationally find the ground-state and excited-state wave functions of the one-dimensional and two-dimensional Bose-Hubbard model. Using this model for a square lattice with
M sites, we obtain the energy spectrum as an analytical function of the on-site Coulomb repulsion,
U, and the total number of particles,
N, from a single training. This approach bypasses the need to solve a new Hamiltonian for each different set of values
(U,N) and generalizes well even for out-of-distribution
(U,N). Using HubbardNet, we identify the two ground-state phases of the Bose-Hubbard model (Mott insulator and superfluid). We show that the DNN-parametrized solutions are in excellent agreement with results from the exact diagonalization of the Hamiltonian, and it outperforms exact diagonalization in terms of computational scaling. These advantages suggest that our model is promising for efficient and accurate computation of exact phase diagrams of many-body lattice Hamiltonians.
Abstract Projected variational wavefunctions such as the Gutzwiller, many-body correlator and Jastrow ansatzes have provided crucial insight into the nature of superfluid-Mott insulator transition in the Bose Hubbard model (BHM) in two or more spatial dimensions. However, these ansatzes have no obvious tractable and systematic way of being improved. A promising alternative is to use Neural-network quantum states (NQS) based on Restricted Boltzmann Machines (RBMs). With binary visible and hidden units NQS have proven to be a highly effective at describing quantum states of interacting spin- 1 2 lattice systems. The application of NQS to bosonic systems has so far been based on one-hot encoding from machine learning where the multi-valued site occupation is distributed across several binary-valued visible units of an RBM. Compared to spin- 1 2 systems one-hot encoding greatly increases the number of variational parameters whilst also making their physical interpretation opaque. Here we revisit the construction of NQS for bosonic systems by reformulating a one-hot encoded RBM into a correlation operator applied to a reference state, analogous to the structure of the projected variational ansatzes. In this form we then propose a number of specialisations of the RBM motivated by the physics of the BHM and the ability to capture exactly the projected variational ansatzes. We analyse in detail the variational performance of these new RBM variants for a 10 × 10 BHM, using both a standard Bose condensate state and a pre-optimised Jastrow + many-body correlator state as the reference state of the calculation. Several of our new ansatzes give robust results as nearly good as one-hot encoding across the regimes of the BHM, but at a substantially reduced cost. Such specialised NQS are thus primed tackle bosonic lattice problems beyond the accuracy of classic variational wavefunctions.
Benchmarking a restricted Boltzmann machine on the
Z2 Bose-Hubbard chain in the adiabatic hard-core regime
We study the ground state of the
Z2 Bose-Hubbard chain in the adiabatic hard-core limit at half filling using variational Monte Carlo with a shallow restricted Boltzmann machine as the variational ansatz. In this context, the neural quantum state is compared with the established adiabatic description of the model. The variational results reproduce the overall structure of the phase diagram obtained from magnetization observables, distinguish the polarized and Néel-ordered regions, and capture representative spin patterns and site occupations for the staggered insulating configurations selected by a weak symmetry-breaking field. Taken together, these results show that a shallow restricted Boltzmann machine reproduces the main adiabatic phase structure of the one-dimensional
Z2 Bose-Hubbard chain and captures the selected symmetry-broken insulating configurations at half filling.
The superfluid-insulator transition in systems of lattice bosons is usually analyzed in the framework of the Bose-Hubbard model, and has been extensively studied by theory and simulations. Less attention has been paid to the remnants of the transition in truncated lattices, with or without periodic boundary conditions. Here we consider the hard-core limit of the extended Bose-Hubbard model on small square and triangular grids -- i.e., sections of the square and triangular lattices containing up to 13 sites. By mapping out the zero-temperature phase diagram through exact diagonalization, we find ground-state characteristics that are markedly different from those emerging in the thermodynamic limit, together with similarities. The dichotomy between superfluid-like and insulating-like behavior is then investigated in two-dimensional systems of a few interacting bosons in the continuum, subject to confining and optical-lattice potentials mimicking the
3×3 square grid. Using path-integral Monte Carlo simulations, we compute kinetic and potential energies, as well as superfluidity and exchange-cycle statistics, finding hints of Bose-Hubbard behavior even in systems of just five particles.
We use a convolutional restricted Boltzmann machine (CRBM) neural network to construct a variational wave function (WF) for the Hubbard model on a square lattice and study it using the variational Monte Carlo (VMC) method. In the wave function, the CRBM acts as a correlation factor to a mean-field BCS state. The number of variational parameters in the WF does not grow automatically with the lattice size and it is computationally much more efficient compared to other neural network based WFs. We find that in the intermediate to strong coupling regime of the model at half-filling, the wave function outperforms even the highly accurate long range backflow-Jastrow correlated wave function. Using the WF, we study the ground state of the half-filled model as a function of onsite Coulomb repulsion
U. We consider two cases for the next-nearest-neighbor hopping parameter, e.g.,
t′=0 as well as a frustrated model case with
t′=0. By examining several quantities, e.g., double occupancy, charge gap, momentum distribution, and spin-spin correlations, we find that the weekly correlated phase in both cases is paramagnetic metallic (PM). As
U is increased, the system undergoes a first-order Mott transition to an insulating state at a critical
Uc, the value of which depends upon
t′. The Mott state in both cases is spin gapped with long range antiferromagnetic (AF) order. Remarkably, the AF order emerges spontaneously from the wave function which does not have any explicitly broken symmetry in it. Apart from some quantitative differences in the results for the two values of
t′, we find some interesting qualitative differences in the way the Mott transition takes place in the two cases.