Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2212.13678 · Physical Review Research

HubbardNet: Efficient predictions of the Bose-Hubbard model spectrum with deep neural networks

Ziyan Zhu, Marios Mattheakis, Weiwei Pan, Efthimios Kaxiras

Published Oct 26, 2023• 9 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevresearch.5.043084

Related papers

Putting machine learning to the test in a quantum many-body system

Yilun Gao, Alberto Rodríguez, Rudolf A. Römer

Abstract Quantum many-body systems pose a formidable computational challenge due to the exponential growth of their Hilbert space. While deep learning (DL) has shown promise as an alternative paradigm, many applications in these systems are still to mature further, beyond the proof-of-concept stages of energy estimation at the lower end of the spectrum. Here, we extend DL beyond this frontier by extensively testing HubbardNet , a deep neural network architecture for the Bose–Hubbard model. Pushing improvements in the optimizer and learning rates, and introducing physics-informed output activations that can resolve extremely small wave-function amplitudes, we achieve ground-state energy errors reduced by orders of magnitude and wave-function fidelities exceeding 99%. We further assess physical relevance by analyzing generalized inverse participation ratios and multifractal dimensions for ground and excited states in one and two dimensions, demonstrating that optimized DL models reproduce localization, delocalization, and multifractality trends across the spectrum. Crucially, these qualitative predictions remain robust across four decades of the interaction strength, e.g. spanning across superfluid, Mott-insulating, as well as quantum chaotic regimes. Together, these results suggest DL as a viable qualitative predictor of many-body structure, complementing the quantitative strengths of exact diagonalization and tensor-network methods.

arXiv:2602.01981•New Journal of Physics•Sep 1, 2026•0 citations

Machine learning density functional theory for the Hubbard model

James Nelson, Rajarshi Tiwari, Stefano Sanvito

The solution of complex many-body lattice models can often be found by defining an energy functional of the relevant density of the problem. For instance, in the case of the Hubbard model the spin-resolved site occupation is enough to describe the system's total energy. Similarly to standard density functional theory, however, the exact functional is unknown, and suitable approximations need to be formulated. By using a deep-learning neural network trained on exact-diagonalization results, we demonstrate that one can construct an exact functional for the Hubbard model. In particular, we show that the neural network returns a ground-state energy numerically indistinguishable from that obtained by exact diagonalization and, most importantly, that the functional satisfies the two Hohenberg-Kohn theorems: for a given ground-state density it yields the external potential, and it is fully variational in the site occupation.

arXiv:1810.12700•Physical review. B./Physical review. B•Feb 15, 2019•44 citations

Ground state phase diagram of the one-dimensional Bose-Hubbard model from restricted Boltzmann machines

Kristopher McBrian, Giuseppe Carleo, Ehsan Khatami

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.

arXiv:1903.03076•Journal of Physics Conference Series•Oct 1, 2019•23 citations

Symmetries and Many-Body Excitations with Neural-Network Quantum States

Kenny Choo, Giuseppe Carleo, Nicolas Regnault +1 more

Artificial neural networks have been recently introduced as a general ansatz to represent many-body wave functions. In conjunction with variational Monte Carlo calculations, this ansatz has been applied to find Hamiltonian ground states and their energies. Here, we provide extensions of this method to study excited states, a central task in several many-body quantum calculations. First, we give a prescription that allows us to target eigenstates of a (nonlocal) symmetry of the Hamiltonian. Second, we give an algorithm to compute low-lying excited states without symmetries. We demonstrate our approach with both restricted Boltzmann machines and feed-forward neural networks. Results are shown for the one-dimensional spin-1/2 Heisenberg model, and for the one-dimensional Bose-Hubbard model. When comparing to exact results, we obtain good agreement for a large range of excited-states energies. Interestingly, we find that deep networks typically outperform shallow architectures for high-energy states.

arXiv:1807.03325•Physical Review Letters•Oct 19, 2018•246 citations

Accurate neural quantum states for interacting lattice bosons

Zakari Denis, Giuseppe Carleo

In recent years, neural quantum states have emerged as a powerful variational approach, achieving state-of-the-art accuracy when representing the ground-state wave function of a great variety of quantum many-body systems, including spin lattices, interacting fermions or continuous-variable systems. However, accurate neural representations of the ground state of interacting bosons on a lattice have remained elusive. We introduce a neural backflow Jastrow Ansatz, in which occupation factors are dressed with translationally equivariant many-body features generated by a deep neural network. We show that this neural quantum state is able to faithfully represent the ground state of the 2D Bose-Hubbard Hamiltonian across all values of the interaction strength. We scale our simulations to lattices of dimension up to 20 × 20 while achieving the best variational energies reported for this model. This enables us to investigate the scaling of the entanglement entropy across the superfluid-to-Mott quantum phase transition, a quantity hard to extract with non-variational approaches.

arXiv:2404.07869•Quantum•Jun 17, 2025•12 citations

Symmetries and Many-Body Excitations with Neural-Network Quantum States

Kenny Choo, Giuseppe Carleo, Nicolas Regnault +1 more

Artificial neural networks have been recently introduced as a general ansatz to represent many-body wave functions. In conjunction with variational Monte Carlo calculations, this ansatz has been applied to find Hamiltonian ground states and their energies. Here, we provide extensions of this method to study excited states, a central task in several many-body quantum calculations. First, we give a prescription that allows us to target eigenstates of a (nonlocal) symmetry of the Hamiltonian. Second, we give an algorithm to compute low-lying excited states without symmetries. We demonstrate our approach with both restricted Boltzmann machines and feed-forward neural networks. Results are shown for the one-dimensional spin-1/2 Heisenberg model, and for the one-dimensional Bose-Hubbard model. When comparing to exact results, we obtain good agreement for a large range of excited-states energies. Interestingly, we find that deep networks typically outperform shallow architectures for high-energy states.

•Oct 19, 2018•0 citations
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 MMM sites, we obtain the energy spectrum as an analytical function of the on-site Coulomb repulsion, UUU, and the total number of particles, NNN, from a single training. This approach bypasses the need to solve a new Hamiltonian for each different set of values (U,N)(U,N)(U,N) and generalizes well even for out-of-distribution (U,N)(U,N)(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 Quantum many-body systems pose a formidable computational challenge due to the exponential growth of their Hilbert space. While deep learning (DL) has shown promise as an alternative paradigm, many applications in these systems are still to mature further, beyond the proof-of-concept stages of energy estimation at the lower end of the spectrum. Here, we extend DL beyond this frontier by extensively testing HubbardNet , a deep neural network architecture for the Bose–Hubbard model. Pushing improvements in the optimizer and learning rates, and introducing physics-informed output activations that can resolve extremely small wave-function amplitudes, we achieve ground-state energy errors reduced by orders of magnitude and wave-function fidelities exceeding 99%. We further assess physical relevance by analyzing generalized inverse participation ratios and multifractal dimensions for ground and excited states in one and two dimensions, demonstrating that optimized DL models reproduce localization, delocalization, and multifractality trends across the spectrum. Crucially, these qualitative predictions remain robust across four decades of the interaction strength, e.g. spanning across superfluid, Mott-insulating, as well as quantum chaotic regimes. Together, these results suggest DL as a viable qualitative predictor of many-body structure, complementing the quantitative strengths of exact diagonalization and tensor-network methods.
The solution of complex many-body lattice models can often be found by defining an energy functional of the relevant density of the problem. For instance, in the case of the Hubbard model the spin-resolved site occupation is enough to describe the system's total energy. Similarly to standard density functional theory, however, the exact functional is unknown, and suitable approximations need to be formulated. By using a deep-learning neural network trained on exact-diagonalization results, we demonstrate that one can construct an exact functional for the Hubbard model. In particular, we show that the neural network returns a ground-state energy numerically indistinguishable from that obtained by exact diagonalization and, most importantly, that the functional satisfies the two Hohenberg-Kohn theorems: for a given ground-state density it yields the external potential, and it is fully variational in the site occupation.
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.
Artificial neural networks have been recently introduced as a general ansatz to represent many-body wave functions. In conjunction with variational Monte Carlo calculations, this ansatz has been applied to find Hamiltonian ground states and their energies. Here, we provide extensions of this method to study excited states, a central task in several many-body quantum calculations. First, we give a prescription that allows us to target eigenstates of a (nonlocal) symmetry of the Hamiltonian. Second, we give an algorithm to compute low-lying excited states without symmetries. We demonstrate our approach with both restricted Boltzmann machines and feed-forward neural networks. Results are shown for the one-dimensional spin-1/2 Heisenberg model, and for the one-dimensional Bose-Hubbard model. When comparing to exact results, we obtain good agreement for a large range of excited-states energies. Interestingly, we find that deep networks typically outperform shallow architectures for high-energy states.
In recent years, neural quantum states have emerged as a powerful variational approach, achieving state-of-the-art accuracy when representing the ground-state wave function of a great variety of quantum many-body systems, including spin lattices, interacting fermions or continuous-variable systems. However, accurate neural representations of the ground state of interacting bosons on a lattice have remained elusive. We introduce a neural backflow Jastrow Ansatz, in which occupation factors are dressed with translationally equivariant many-body features generated by a deep neural network. We show that this neural quantum state is able to faithfully represent the ground state of the 2D Bose-Hubbard Hamiltonian across all values of the interaction strength. We scale our simulations to lattices of dimension up to 20 × 20 while achieving the best variational energies reported for this model. This enables us to investigate the scaling of the entanglement entropy across the superfluid-to-Mott quantum phase transition, a quantity hard to extract with non-variational approaches.
Artificial neural networks have been recently introduced as a general ansatz to represent many-body wave functions. In conjunction with variational Monte Carlo calculations, this ansatz has been applied to find Hamiltonian ground states and their energies. Here, we provide extensions of this method to study excited states, a central task in several many-body quantum calculations. First, we give a prescription that allows us to target eigenstates of a (nonlocal) symmetry of the Hamiltonian. Second, we give an algorithm to compute low-lying excited states without symmetries. We demonstrate our approach with both restricted Boltzmann machines and feed-forward neural networks. Results are shown for the one-dimensional spin-1/2 Heisenberg model, and for the one-dimensional Bose-Hubbard model. When comparing to exact results, we obtain good agreement for a large range of excited-states energies. Interestingly, we find that deep networks typically outperform shallow architectures for high-energy states.