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:1611.10185 · Physical review. B./Physical review. B

Infinite occupation number basis of bosons: Solving a numerical challenge

Andreas Geißler, Walter Hofstetter

Published Jun 29, 2017• 1 citation
arXivHTMLPDF

Abstract

In any bosonic lattice system, which is not dominated by local interactions and thus ``frozen'' in a Mott-type state, numerical methods have to cope with the infinite size of the corresponding Hilbert space even for finite lattice sizes. While it is common practice to restrict the local occupation number basis to ${N}_{c}$ lowest occupied states, the presence of a finite condensate fraction requires the complete number basis for an exact representation of the many-body ground state. In this work we present a truncation scheme to account for contributions from higher number states. By simply adding a single coherent-tail state to this common truncation, we demonstrate increased numerical accuracy and the possible increase in numerical efficiency of this method for the Gutzwiller variational wave function and within dynamical mean-field theory.

DOI
10.1103/physrevb.95.224516

Related papers

Delocalization and superfluidity of ultracold bosonic atoms in a ring lattice

Fernanda Pinheiro, A. F. R. de Toledo Piza

Properties of bosonic atoms in small systems with a periodic quasi one-dimensional circular toroidal lattice potential subjected to rotation are examined by performing exact diagonalization in a truncated many body space. The expansion of the many-body Hamiltonian is considered in terms of the first band Bloch functions, and no assumption regarding restriction to nearest-neighbor hopping (tight-binding approximation) is involved. A finite size version of the zero temperature phase diagrams of Fisher et al. \cite{Fisher} is obtained and the results, in remarkable quantitative correspondence with the results available for larger systems, discussed. Ground state properties relating to superfluidity are examined in the context of two-fluid phenomenology. The basic tool, consisting of the intrinsic inertia associated with small rotation angular velocities in the lab frame, is used to obtain ground state `superfluid fractions' numerically. They are analytically associated with one-body, uniform solenoidal currents in the case of the adopted geometry. These currents are in general incoherent superpositions of contributions from each eigenstate of the associated reduced one-body densities, with the corresponding occupation numbers as weights. Full coherence occurs therefore only when only one eigenstate is occupied by all bosons. The obtained numerical values for the superfluid fractions remain small throughout the parameter region corresponding to the `Mott insulator to superfluid' transition, and saturate at unity only as the lattice is completely smoothed out.

arXiv:1301.4672•Journal of Physics B Atomic Molecular and Optical Physics•Oct 15, 2013•7 citations

Exact solution of the Bose-Hubbard model on the Bethe lattice

Guilhem Semerjian, Marco Tarzia, Francesco Zamponi

The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.

arXiv:0904.3075•Physical Review B•Jul 31, 2009•42 citations

THE BOSONIC STRUCTURE OF FERMIONS

J. RUAN, R. J. CREWTHER

We bosonize fermions by identifying their occupation numbers as the binary digits of a Bose occupation number. Unlike other schemes, our method allows infinitely many fermionic oscillators to be constructed from just one bosonic oscillator.

arXiv:hep-th/9407178•Modern Physics Letters A•Oct 30, 1994•2 citations

The Bosonic Structure of Fermions

J. Ruan, R. J. Crewther

We bosonize fermions by identifying their occupation numbers as the binary digits of a Bose occupation number. Unlike other schemes, our method allows infinitely many fermionic oscillators to be constructed from just one bosonic oscillator.

•Jan 1, 1994•2 citations

Bosonic field digitization for quantum computers

Alexandru Macridin, Andy C. Y. Li, S. Mrenna +1 more

Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods that can accurately predict error scaling as a function of the resolution and parameters of the model and that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation number cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite-size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wave function sampling, and the finite representation size. Numerical simulations of small size ${\mathrm{\ensuremath{\Phi}}}^{4}$ problems demonstrate that the boson mass optimizing the sampling of the ground state wave function is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wave functions does not necessarily result in accurate representation. We develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.

arXiv:2108.10793•Physical Review A•May 4, 2022•23 citations

Hermitian boson mapping and finite truncation

Calvin W. Johnson, Joseph N. Ginocchio

Starting from a general, microscopic fermion-to-boson mapping that preserves Hermitian conjugation, we discuss truncations of the boson Fock space basis. We give conditions under which the exact boson images of finite fermion operators are also finite [e.g., a (1+2)-body fermion Hamiltonian is mapped to a (1+2)-body boson Hamiltonian] in the truncated basis. For the most general case, where the image is not necessarily exactly finite, we discuss how to make practical and controlled approximations.

arXiv:nucl-th/9405026•Physical Review C•Aug 1, 1994•7 citations
In any bosonic lattice system, which is not dominated by local interactions and thus ``frozen'' in a Mott-type state, numerical methods have to cope with the infinite size of the corresponding Hilbert space even for finite lattice sizes. While it is common practice to restrict the local occupation number basis to Nc{N}_{c}Nc​ lowest occupied states, the presence of a finite condensate fraction requires the complete number basis for an exact representation of the many-body ground state. In this work we present a truncation scheme to account for contributions from higher number states. By simply adding a single coherent-tail state to this common truncation, we demonstrate increased numerical accuracy and the possible increase in numerical efficiency of this method for the Gutzwiller variational wave function and within dynamical mean-field theory.
Properties of bosonic atoms in small systems with a periodic quasi one-dimensional circular toroidal lattice potential subjected to rotation are examined by performing exact diagonalization in a truncated many body space. The expansion of the many-body Hamiltonian is considered in terms of the first band Bloch functions, and no assumption regarding restriction to nearest-neighbor hopping (tight-binding approximation) is involved. A finite size version of the zero temperature phase diagrams of Fisher et al. \cite{Fisher} is obtained and the results, in remarkable quantitative correspondence with the results available for larger systems, discussed. Ground state properties relating to superfluidity are examined in the context of two-fluid phenomenology. The basic tool, consisting of the intrinsic inertia associated with small rotation angular velocities in the lab frame, is used to obtain ground state `superfluid fractions' numerically. They are analytically associated with one-body, uniform solenoidal currents in the case of the adopted geometry. These currents are in general incoherent superpositions of contributions from each eigenstate of the associated reduced one-body densities, with the corresponding occupation numbers as weights. Full coherence occurs therefore only when only one eigenstate is occupied by all bosons. The obtained numerical values for the superfluid fractions remain small throughout the parameter region corresponding to the `Mott insulator to superfluid' transition, and saturate at unity only as the lattice is completely smoothed out.
The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.
Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods that can accurately predict error scaling as a function of the resolution and parameters of the model and that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation number cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite-size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wave function sampling, and the finite representation size. Numerical simulations of small size ${\mathrm{\ensuremath{\Phi}}}^{4}$ problems demonstrate that the boson mass optimizing the sampling of the ground state wave function is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wave functions does not necessarily result in accurate representation. We develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.
Starting from a general, microscopic fermion-to-boson mapping that preserves Hermitian conjugation, we discuss truncations of the boson Fock space basis. We give conditions under which the exact boson images of finite fermion operators are also finite [e.g., a (1+2)-body fermion Hamiltonian is mapped to a (1+2)-body boson Hamiltonian] in the truncated basis. For the most general case, where the image is not necessarily exactly finite, we discuss how to make practical and controlled approximations.