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:1103.3664 · Physical Review B

Extended self-energy functional approach for strongly correlated lattice bosons in the superfluid phase

Enrico Arrigoni, Michael Knap, Wolfgang von der Linden

Published Jul 29, 2011• 17 citations
arXivHTMLPDF

Abstract

Among the various numerical techniques to study the physics of strongly correlated quantum many-body systems, the self-energy functional approach (SFA) has become increasingly important. In its previous form, however, SFA is not applicable to Bose-Einstein condensation or superfluidity. In this paper, we show how to overcome this shortcoming. To this end, we identify an appropriate quantity, which we term $D$, that represents the correlation correction of the condensate order parameter, as it does the self-energy for Green's function. An appropriate functional is derived, which is stationary at the exact physical realization of $D$ and of the self-energy. Its derivation is based on a functional-integral representation of the grand potential followed by an appropriate sequence of Legendre transformations. The approach is not perturbative and, therefore, applicable to a wide range of models with local interactions. We show that the variational cluster approach based on the extended self-energy functional is equivalent to the ``pseudoparticle'' approach proposed in Phys. Rev. B 83, 134507 (2011). We present results for the superfluid density in the two-dimensional Bose-Hubbard model, which shows a remarkable agreement with those of quantum-Monte-Carlo calculations.

DOI
10.1103/physrevb.84.014535

Related papers

Variational cluster approach for strongly correlated lattice bosons in the superfluid phase

Michael Knap, Enrico Arrigoni, Wolfgang von der Linden

We extend the variational cluster approach to deal with strongly correlated lattice bosons in the superfluid phase. To this end, we reformulate the approach within a pseudoparticle formalism, whereby cluster excitations are described by particlelike excitations. The approximation amounts to solving a multicomponent noninteracting bosonic system by means of a multimode Bogoliubov approximation. A source-and-drain term is introduced in order to break U(1) symmetry at the cluster level. We provide an expression for the grand potential, the single-particle normal and anomalous Green's functions, the condensate density, and other static quantities. As a first nontrivial application of the method we choose the two-dimensional Bose-Hubbard model and evaluate results in both the Mott and the superfluid phases. Our results show an excellent agreement with quantum Monte Carlo calculations.

arXiv:1010.4295•Physical Review B•Apr 7, 2011•29 citations

Self-energy functional theory with symmetry breaking for disordered lattice bosons

Dario Hügel, Hugo U. R. Strand, Lode Pollet

We extend the self-energy functional theory to the case of interacting lattice bosons in the presence of symmetry breaking and quenched disorder. The self-energy functional we derive depends only on the self-energies of the disorder-averaged propagators, allowing for the construction of general non-perturbative approximations. Using a simple single-site reference system with only three variational parameters, we are able to reproduce numerically exact quantum Monte Carlo (QMC) results on local observables of the Bose–Hubbard model with box disorder with high accuracy. At strong interactions, the phase boundaries are reproduced qualitatively but shifted with respect to the ones observed with QMC due to the extremely low condensate fraction in the superfluid phase. Deep in the strongly-disordered weakly-interacting regime, the simple reference system employed is insufficient and no stationary solutions can be found within its restricted variational subspace. By systematically analyzing thermodynamical observables and the spectral function, we find that the strongly interacting Bose glass is characterized by different regimes, depending on which local occupations are activated as a function of the disorder strength. We find that the particles delocalize into isolated superfluid lakes over a strongly localized background around maximally-occupied sites whenever these sites are particularly rare. Our results indicate that the transition from the Bose glass to the superfluid phase around unit filling at strong interactions is driven by the percolation of superfluid lakes which form around doubly occupied sites.

arXiv:1801.07274•Quantum Science and Technology•Apr 25, 2018•2 citations

Self-consistent methods for interacting lattice bosons with U(1)-symmetry-breaking

Dario Hügel

This thesis is dedicated to the derivation and benchmarking of self-consistent numerical methods that can be applied to interacting bosonic lattice models. The central goal is to derive methods with low numerical complexity but high accuracy, to be applied to complex systems which are out-of reach for established methods such as path integral quantum Monte Carlo (QMC) or the density matrix renormalization group. In the first part we derive the self-energy functional theory (SFT) for bosons. Building upon previous works on lattice systems without U(1)-symmetry-breaking, we systematically extend SFT to the possibility of a broken U(1)-symmetry and the presence of disorder. SFT incorporates bosonic dynamical mean-field theory as a certain limit, and represents a general non-perturbative framework, enabling the construction of diagrammatically sound approximations in the thermodynamical limit that are controlled in the number of optimization parameters. Using just three variational parameters, we are able to study the Bose-Hubbard model both in its clean version and in the presence of local disorder, showing excellent agreement with numerically exact QMC results. We systematically analyze the corresponding spectral functions, which cannot be fully captured by QMC. In particular, we find that in the presence of disorder the phase transition from the Bose glass to the superfluid phase at strong interactions is driven by the percolation of superfluid lakes which form around doubly occupied sites, leading to a small condensate fraction over a strongly-localized background. The second part is dedicated to the derivation of reciprocal cluster mean-field theory (RCMF) and its application to the strongly-interacting Harper-Hofstadter-Mott model (HHMm). In RCMF the full lattice in the thermodynamical limit is projected onto finite-size clusters, which are decoupled in reciprocal space through a mean-field decoupling approximation, crucially preserving the symmetries of the non-interacting dispersion. The resulting groundstate phase diagram of the HHMm exhibits band insulating, striped superfluid, and supersolid phases. Furthermore, we observe gapless uncondensed liquid phases at integer fillings, and a metastable competing fractional quantum Hall (fQH) phase. The fQH phase, predicted as the groundstate by other methods, is most likely underestimated by RCMF. We then show how a quasi-one-dimensional geometry stabilizes gapped topologically non-trivial groundstates in the HHMm. We observe quasi-one-dimensional analogues of fQH phases at fillings 1/2 and 3/2, and unconventional gapped non-degenerate groundstates at integer filling with quantized Hall responses. By systematically comparing results computed with RCMF and exact diagonalization (ED), we are able to give conclusive quantitative answers on the phase boundaries of the system, as the two methods approach the thermodynamical limit from opposite sides, since RCMF favours gapless and ED gapped phases.

•Jan 1, 2018•0 citations

Bosonic self-energy functional theory

Dario Hügel, Philipp Werner, Lode Pollet +1 more

We derive the self-energy functional theory for bosonic lattice systems with broken U(1) symmetry by parametrizing the bosonic Baym-Kadanoff effective action in terms of one- and two-point self-energies. The formalism goes beyond other approximate methods such as the pseudoparticle variational cluster approximation, the cluster composite boson mapping, and the Bogoliubov+U theory. It simplifies to bosonic dynamical-mean-field theory when constraining to local fields, whereas when neglecting kinetic contributions of noncondensed bosons, it reduces to the static mean-field approximation. To benchmark the theory, we study the Bose-Hubbard model on the two- and three-dimensional cubic lattice, comparing with exact results from path integral quantum Monte Carlo. We also study the frustrated square lattice with next-nearest-neighbor hopping, which is beyond the reach of Monte Carlo simulations. A reference system comprising a single bosonic state, corresponding to three variational parameters, is sufficient to quantitatively describe phase boundaries and thermodynamical observables, while qualitatively capturing the spectral functions, as well as the enhancement of kinetic fluctuations in the frustrated case. On the basis of these findings, we propose self-energy functional theory as the omnibus framework for treating bosonic lattice models, in particular, in cases where path integral quantum Monte Carlo methods suffer from severe sign problems (e.g., in the presence of nontrivial gauge fields or frustration). Self-energy functional theory enables the construction of diagrammatically sound approximations that are quantitatively precise and controlled in the number of optimization parameters but nevertheless remain computable by modest means.

arXiv:1603.01401•Physical review. B./Physical review. B•Nov 9, 2016•14 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

Benchmarking the variational cluster approach by means of the one-dimensional Bose-Hubbard model

Michael Knap, Enrico Arrigoni, Wolfgang von der Linden

Convergence properties of the variational cluster approach with respect to the variational parameter space, cluster size, and boundary conditions of the reference system are investigated and discussed for bosonic many-body systems. Specifically, the variational cluster approach is applied to the one-dimensional Bose-Hubbard model, which exhibits a quantum phase transition from Mott to superfluid phase. In order to benchmark the variational cluster approach, results for the phase boundary delimiting the first Mott lobe are compared with essentially exact density matrix renormalization group data. Furthermore, static quantities, such as the ground state energy and the one-particle density matrix are compared with high-order strong coupling perturbation theory results. For reference systems with open boundary conditions the variational parameter space is extended by an additional variational parameter which allows for a more uniform particle density on the reference system and thus drastically improves the results. It turns out that the variational cluster approach yields accurate results with relatively low-computational effort for both the phase boundary as well as the static properties of the one-dimensional Bose-Hubbard model, even at the tip of the first Mott lobe where correlation effects are most pronounced.

arXiv:1002.2449•Physical Review B•Jun 17, 2010•15 citations
Among the various numerical techniques to study the physics of strongly correlated quantum many-body systems, the self-energy functional approach (SFA) has become increasingly important. In its previous form, however, SFA is not applicable to Bose-Einstein condensation or superfluidity. In this paper, we show how to overcome this shortcoming. To this end, we identify an appropriate quantity, which we term DDD, that represents the correlation correction of the condensate order parameter, as it does the self-energy for Green's function. An appropriate functional is derived, which is stationary at the exact physical realization of DDD and of the self-energy. Its derivation is based on a functional-integral representation of the grand potential followed by an appropriate sequence of Legendre transformations. The approach is not perturbative and, therefore, applicable to a wide range of models with local interactions. We show that the variational cluster approach based on the extended self-energy functional is equivalent to the ``pseudoparticle'' approach proposed in Phys. Rev. B 83, 134507 (2011). We present results for the superfluid density in the two-dimensional Bose-Hubbard model, which shows a remarkable agreement with those of quantum-Monte-Carlo calculations.
We extend the variational cluster approach to deal with strongly correlated lattice bosons in the superfluid phase. To this end, we reformulate the approach within a pseudoparticle formalism, whereby cluster excitations are described by particlelike excitations. The approximation amounts to solving a multicomponent noninteracting bosonic system by means of a multimode Bogoliubov approximation. A source-and-drain term is introduced in order to break U(1) symmetry at the cluster level. We provide an expression for the grand potential, the single-particle normal and anomalous Green's functions, the condensate density, and other static quantities. As a first nontrivial application of the method we choose the two-dimensional Bose-Hubbard model and evaluate results in both the Mott and the superfluid phases. Our results show an excellent agreement with quantum Monte Carlo calculations.
We extend the self-energy functional theory to the case of interacting lattice bosons in the presence of symmetry breaking and quenched disorder. The self-energy functional we derive depends only on the self-energies of the disorder-averaged propagators, allowing for the construction of general non-perturbative approximations. Using a simple single-site reference system with only three variational parameters, we are able to reproduce numerically exact quantum Monte Carlo (QMC) results on local observables of the Bose–Hubbard model with box disorder with high accuracy. At strong interactions, the phase boundaries are reproduced qualitatively but shifted with respect to the ones observed with QMC due to the extremely low condensate fraction in the superfluid phase. Deep in the strongly-disordered weakly-interacting regime, the simple reference system employed is insufficient and no stationary solutions can be found within its restricted variational subspace. By systematically analyzing thermodynamical observables and the spectral function, we find that the strongly interacting Bose glass is characterized by different regimes, depending on which local occupations are activated as a function of the disorder strength. We find that the particles delocalize into isolated superfluid lakes over a strongly localized background around maximally-occupied sites whenever these sites are particularly rare. Our results indicate that the transition from the Bose glass to the superfluid phase around unit filling at strong interactions is driven by the percolation of superfluid lakes which form around doubly occupied sites.
This thesis is dedicated to the derivation and benchmarking of self-consistent numerical methods that can be applied to interacting bosonic lattice models. The central goal is to derive methods with low numerical complexity but high accuracy, to be applied to complex systems which are out-of reach for established methods such as path integral quantum Monte Carlo (QMC) or the density matrix renormalization group. In the first part we derive the self-energy functional theory (SFT) for bosons. Building upon previous works on lattice systems without U(1)-symmetry-breaking, we systematically extend SFT to the possibility of a broken U(1)-symmetry and the presence of disorder. SFT incorporates bosonic dynamical mean-field theory as a certain limit, and represents a general non-perturbative framework, enabling the construction of diagrammatically sound approximations in the thermodynamical limit that are controlled in the number of optimization parameters. Using just three variational parameters, we are able to study the Bose-Hubbard model both in its clean version and in the presence of local disorder, showing excellent agreement with numerically exact QMC results. We systematically analyze the corresponding spectral functions, which cannot be fully captured by QMC. In particular, we find that in the presence of disorder the phase transition from the Bose glass to the superfluid phase at strong interactions is driven by the percolation of superfluid lakes which form around doubly occupied sites, leading to a small condensate fraction over a strongly-localized background. The second part is dedicated to the derivation of reciprocal cluster mean-field theory (RCMF) and its application to the strongly-interacting Harper-Hofstadter-Mott model (HHMm). In RCMF the full lattice in the thermodynamical limit is projected onto finite-size clusters, which are decoupled in reciprocal space through a mean-field decoupling approximation, crucially preserving the symmetries of the non-interacting dispersion. The resulting groundstate phase diagram of the HHMm exhibits band insulating, striped superfluid, and supersolid phases. Furthermore, we observe gapless uncondensed liquid phases at integer fillings, and a metastable competing fractional quantum Hall (fQH) phase. The fQH phase, predicted as the groundstate by other methods, is most likely underestimated by RCMF. We then show how a quasi-one-dimensional geometry stabilizes gapped topologically non-trivial groundstates in the HHMm. We observe quasi-one-dimensional analogues of fQH phases at fillings 1/2 and 3/2, and unconventional gapped non-degenerate groundstates at integer filling with quantized Hall responses. By systematically comparing results computed with RCMF and exact diagonalization (ED), we are able to give conclusive quantitative answers on the phase boundaries of the system, as the two methods approach the thermodynamical limit from opposite sides, since RCMF favours gapless and ED gapped phases.
We derive the self-energy functional theory for bosonic lattice systems with broken U(1) symmetry by parametrizing the bosonic Baym-Kadanoff effective action in terms of one- and two-point self-energies. The formalism goes beyond other approximate methods such as the pseudoparticle variational cluster approximation, the cluster composite boson mapping, and the Bogoliubov+U theory. It simplifies to bosonic dynamical-mean-field theory when constraining to local fields, whereas when neglecting kinetic contributions of noncondensed bosons, it reduces to the static mean-field approximation. To benchmark the theory, we study the Bose-Hubbard model on the two- and three-dimensional cubic lattice, comparing with exact results from path integral quantum Monte Carlo. We also study the frustrated square lattice with next-nearest-neighbor hopping, which is beyond the reach of Monte Carlo simulations. A reference system comprising a single bosonic state, corresponding to three variational parameters, is sufficient to quantitatively describe phase boundaries and thermodynamical observables, while qualitatively capturing the spectral functions, as well as the enhancement of kinetic fluctuations in the frustrated case. On the basis of these findings, we propose self-energy functional theory as the omnibus framework for treating bosonic lattice models, in particular, in cases where path integral quantum Monte Carlo methods suffer from severe sign problems (e.g., in the presence of nontrivial gauge fields or frustration). Self-energy functional theory enables the construction of diagrammatically sound approximations that are quantitatively precise and controlled in the number of optimization parameters but nevertheless remain computable by modest means.
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.
Convergence properties of the variational cluster approach with respect to the variational parameter space, cluster size, and boundary conditions of the reference system are investigated and discussed for bosonic many-body systems. Specifically, the variational cluster approach is applied to the one-dimensional Bose-Hubbard model, which exhibits a quantum phase transition from Mott to superfluid phase. In order to benchmark the variational cluster approach, results for the phase boundary delimiting the first Mott lobe are compared with essentially exact density matrix renormalization group data. Furthermore, static quantities, such as the ground state energy and the one-particle density matrix are compared with high-order strong coupling perturbation theory results. For reference systems with open boundary conditions the variational parameter space is extended by an additional variational parameter which allows for a more uniform particle density on the reference system and thus drastically improves the results. It turns out that the variational cluster approach yields accurate results with relatively low-computational effort for both the phase boundary as well as the static properties of the one-dimensional Bose-Hubbard model, even at the tip of the first Mott lobe where correlation effects are most pronounced.