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

Numerical calculation of spectral functions of the Bose-Hubbard model using bosonic dynamical mean-field theory

Jaromir Panas, Anna Kauch, J. Kuneš, D. Vollhardt, Krzysztof Byczuk

Published Jul 6, 2015• 13 citations
arXivHTMLPDF

Abstract

We calculate the momentum dependent spectral function of the Bose-Hubbard model on a simple cubic lattice in three dimensions within the bosonic dynamical mean-field theory (B-DMFT). The continuous-time quantum Monte Carlo method is used to solve the self-consistent B-DMFT equations together with the maximum entropy method for the analytic continuation to real frequencies. Results for weak, intermediate, and strong interactions are presented. In the limit of weak and strong interactions very good agreement with results obtained by perturbation theory is found. By contrast, at intermediate interactions the results differ significantly, indicating that in this regime perturbative methods fail to describe the dynamics of interacting bosons.

DOI
10.1103/physrevb.92.045102

Related papers

Strong-coupling solution of the bosonic dynamical mean-field theory

Anna Kauch, Krzysztof Byczuk, D. Vollhardt

We derive an approximate analytical solution of the self-consistency equations of the bosonic dynamical mean-field theory (B-DMFT) in the strong-coupling limit. The approach is based on a linked-cluster expansion in the hybridization function of normal bosons around the atomic limit. The solution is used to compute the phase diagram of the bosonic Hubbard model for different lattices. We compare our results with numerical solutions of the B-DMFT equations and numerically exact methods, respectively. The very good agreement with those numerical results demonstrates that our approach captures the essential physics of correlated bosons both in the Mott insulator and in the superfluid phase. Close to the transition into the superfluid phase the momentum distribution function at zero momentum is found to be strongly enhanced already in the normal phase. The linked-cluster expansion also allows us to compute dynamical properties such as the spectral function of bosons. The evolution of the spectral function across the transition from the normal to the superfluid phase is seen to be characteristically different for the interaction driven and density driven transition, respectively.

arXiv:1203.1236•Physical Review B•May 9, 2012•8 citations

Dynamical mean-field theory for bosons

Peter Anders, Emanuel Gull, Lode Pollet +2 more

We discuss the recently developed bosonic dynamical mean-field theory (B-DMFT) framework, which maps a bosonic lattice model onto the self-consistent solution of a bosonic impurity model with coupling to a reservoir of normal and condensed bosons. The effective impurity action is derived in several ways: (i) as an approximation to the kinetic energy functional of the lattice problem, (ii) using a cavity approach and (iii) using an effective medium approach based on adding a one-loop correction to the self-consistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on the sampling of a perturbation expansion in the hybridization functions and the condensate wave function. As applications of the formalism, we present finite-temperature B-DMFT phase diagrams for the bosonic Hubbard model on a three-dimensional (3D) cubic and a 2D square lattice, the condensate order parameter as a function of chemical potential, critical exponents for the condensate, the approach to the weakly interacting Bose gas regime for weak repulsions and the kinetic energy as a function of temperature.

arXiv:1103.0017•New Journal of Physics•Jul 29, 2011•62 citations

Thermodynamic and spectral properties of quantum many-particle systems

Sebastian Fuchs

The investigation of quantum many-particle systems is a central goal of modern condensed-matter physics. Complex physical phenomena such as high-temperature superconductivity, heavy fermion behavior, and correlation-driven metal-insulator transitions are central aspects of this particular field of science. The Hubbard model is a paradigm of correlated-electron physics. It is one of the simplest quantum-mechanical lattice models capable of capturing relevant physical aspects of strongly correlated electron systems. The Hubbard model also plays an important role in the explanation of many-body phenomena observed in ultra-cold atoms trapped in optical lattices. In this thesis, we study the three-dimensional Hubbard model using quantum Monte Carlo simulations (QMC). Our aim is a precise numerical study of its properties in the thermodynamic limit. However, a direct investigation of lattice models in the thermodynamic limit is generally impossible because QMC simulations cannot treat infinite lattice sizes. Therefore, approximative schemes that allow calculations directly in the thermodynamic limit are advantageous. A widely applied approximation of this kind is the dynamical mean-field theory (DMFT) which maps the infinite-lattice problem onto a single-site impurity model embedded in a mean field. We apply cluster mean-field theories which represent a systematic extension of the DMFT by expanding the single impurity to a finite cluster. The second focus of this thesis is the calculation of spectral properties of the Hubbard model. QMC methods map the quantum-mechanical lattice model on a classical one at the expense of an additional imaginary time dimension. For this reason, QMC algorithms can only provide dynamical data on the imaginary time or frequency axis. The necessary analytic continuation to the physically relevant real axis has proven to be a difficult problem that has to be approached by advanced data analysis tools. The maximum-entropy method (MEM) is the standard technique to handle this problem. We apply the MEM to the problem of extracting momentum-resolved single particle spectra from QMC simulations of the 3D Hubbard model. Furthermore, we discuss and expand algorithms that were recently proposed as alternatives to the standard MEM approach.

•Jan 1, 2011•5 citations

Mean-Field Dynamics of the Bose–Hubbard Model in High Dimension

Shahnaz Farhat, Denis Périce, Sören Petrat

Abstract This paper is a first step toward rigorously establishing the validity of mean-field approximations for the dynamics of quantum systems in a large dimensional limit. Our result is a derivation of the large d mean-field limit of the Bose–Hubbard model on a square lattice. The Bose–Hubbard model effectively describes bosons on a lattice with on-site interactions and nearest-neighbour hopping, serving as a foundational framework for understanding strong particle interactions and the superfluid to Mott insulator phase transition. We prove a trace norm estimate between the one-lattice-site reduced density of the Schrödinger dynamics and the mean-field dynamics in the limit of large dimension. Here, the mean-field approximation is in the hopping amplitude and not in the interaction, leading to a very rich and non-trivial mean-field equation. This mean-field equation does not only describe the condensate, as is the case when the mean-field description comes from a large particle number limit averaging out the interaction, but it allows for a phase transition since it contains the full non-trivial interaction. Our work is a rigorous justification of a simple case of the highly successful dynamical mean-field theory (DMFT) for bosons, which somewhat surprisingly yields many qualitatively correct results in three dimensions.

arXiv:2501.05304•Annales Henri Poincaré•Oct 10, 2026•0 citations

Dynamical Mean Field Solution of the Bose-Hubbard Model

Peter Anders, Emanuel Gull, Lode Pollet +2 more

We present the effective action and self-consistency equations for the bosonic dynamical mean field approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations, we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.

arXiv:1004.0510•Physical Review Letters•Aug 24, 2010•89 citations

Towards out of equilibrium impurity solver for the Dynamical Mean Field Theory

Karpov, A., G. Sultanov, Evgeny Polyakov +1 more

The Dynamical Mean Field Theory (DMFT) is a powerful tool for calculating highly correlated systems (both bosonic and fermionic) in a state of thermodynamic equilibrium. However, in the case of non-equilibrium states, the method has significant limitations that do not allow obtaining correct results. The stumbling block here is the impurity solver: a method for calculating the dynamics of an open system. In this work we present the prototype of a universal impurity solver for the DMFT method which can solve both equilibrium and non-equilibrium problems. We analyse spectral functions of Bose-Hubbard model on the Bethe lattice in different regimes and show that our solver gives correct and physical results.

arXiv:2311.02623•Nov 5, 2023•0 citations
We calculate the momentum dependent spectral function of the Bose-Hubbard model on a simple cubic lattice in three dimensions within the bosonic dynamical mean-field theory (B-DMFT). The continuous-time quantum Monte Carlo method is used to solve the self-consistent B-DMFT equations together with the maximum entropy method for the analytic continuation to real frequencies. Results for weak, intermediate, and strong interactions are presented. In the limit of weak and strong interactions very good agreement with results obtained by perturbation theory is found. By contrast, at intermediate interactions the results differ significantly, indicating that in this regime perturbative methods fail to describe the dynamics of interacting bosons.
We derive an approximate analytical solution of the self-consistency equations of the bosonic dynamical mean-field theory (B-DMFT) in the strong-coupling limit. The approach is based on a linked-cluster expansion in the hybridization function of normal bosons around the atomic limit. The solution is used to compute the phase diagram of the bosonic Hubbard model for different lattices. We compare our results with numerical solutions of the B-DMFT equations and numerically exact methods, respectively. The very good agreement with those numerical results demonstrates that our approach captures the essential physics of correlated bosons both in the Mott insulator and in the superfluid phase. Close to the transition into the superfluid phase the momentum distribution function at zero momentum is found to be strongly enhanced already in the normal phase. The linked-cluster expansion also allows us to compute dynamical properties such as the spectral function of bosons. The evolution of the spectral function across the transition from the normal to the superfluid phase is seen to be characteristically different for the interaction driven and density driven transition, respectively.
We discuss the recently developed bosonic dynamical mean-field theory (B-DMFT) framework, which maps a bosonic lattice model onto the self-consistent solution of a bosonic impurity model with coupling to a reservoir of normal and condensed bosons. The effective impurity action is derived in several ways: (i) as an approximation to the kinetic energy functional of the lattice problem, (ii) using a cavity approach and (iii) using an effective medium approach based on adding a one-loop correction to the self-consistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on the sampling of a perturbation expansion in the hybridization functions and the condensate wave function. As applications of the formalism, we present finite-temperature B-DMFT phase diagrams for the bosonic Hubbard model on a three-dimensional (3D) cubic and a 2D square lattice, the condensate order parameter as a function of chemical potential, critical exponents for the condensate, the approach to the weakly interacting Bose gas regime for weak repulsions and the kinetic energy as a function of temperature.
The investigation of quantum many-particle systems is a central goal of modern condensed-matter physics. Complex physical phenomena such as high-temperature superconductivity, heavy fermion behavior, and correlation-driven metal-insulator transitions are central aspects of this particular field of science. The Hubbard model is a paradigm of correlated-electron physics. It is one of the simplest quantum-mechanical lattice models capable of capturing relevant physical aspects of strongly correlated electron systems. The Hubbard model also plays an important role in the explanation of many-body phenomena observed in ultra-cold atoms trapped in optical lattices. In this thesis, we study the three-dimensional Hubbard model using quantum Monte Carlo simulations (QMC). Our aim is a precise numerical study of its properties in the thermodynamic limit. However, a direct investigation of lattice models in the thermodynamic limit is generally impossible because QMC simulations cannot treat infinite lattice sizes. Therefore, approximative schemes that allow calculations directly in the thermodynamic limit are advantageous. A widely applied approximation of this kind is the dynamical mean-field theory (DMFT) which maps the infinite-lattice problem onto a single-site impurity model embedded in a mean field. We apply cluster mean-field theories which represent a systematic extension of the DMFT by expanding the single impurity to a finite cluster. The second focus of this thesis is the calculation of spectral properties of the Hubbard model. QMC methods map the quantum-mechanical lattice model on a classical one at the expense of an additional imaginary time dimension. For this reason, QMC algorithms can only provide dynamical data on the imaginary time or frequency axis. The necessary analytic continuation to the physically relevant real axis has proven to be a difficult problem that has to be approached by advanced data analysis tools. The maximum-entropy method (MEM) is the standard technique to handle this problem. We apply the MEM to the problem of extracting momentum-resolved single particle spectra from QMC simulations of the 3D Hubbard model. Furthermore, we discuss and expand algorithms that were recently proposed as alternatives to the standard MEM approach.
Abstract This paper is a first step toward rigorously establishing the validity of mean-field approximations for the dynamics of quantum systems in a large dimensional limit. Our result is a derivation of the large d mean-field limit of the Bose–Hubbard model on a square lattice. The Bose–Hubbard model effectively describes bosons on a lattice with on-site interactions and nearest-neighbour hopping, serving as a foundational framework for understanding strong particle interactions and the superfluid to Mott insulator phase transition. We prove a trace norm estimate between the one-lattice-site reduced density of the Schrödinger dynamics and the mean-field dynamics in the limit of large dimension. Here, the mean-field approximation is in the hopping amplitude and not in the interaction, leading to a very rich and non-trivial mean-field equation. This mean-field equation does not only describe the condensate, as is the case when the mean-field description comes from a large particle number limit averaging out the interaction, but it allows for a phase transition since it contains the full non-trivial interaction. Our work is a rigorous justification of a simple case of the highly successful dynamical mean-field theory (DMFT) for bosons, which somewhat surprisingly yields many qualitatively correct results in three dimensions.
We present the effective action and self-consistency equations for the bosonic dynamical mean field approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations, we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.
The Dynamical Mean Field Theory (DMFT) is a powerful tool for calculating highly correlated systems (both bosonic and fermionic) in a state of thermodynamic equilibrium. However, in the case of non-equilibrium states, the method has significant limitations that do not allow obtaining correct results. The stumbling block here is the impurity solver: a method for calculating the dynamics of an open system. In this work we present the prototype of a universal impurity solver for the DMFT method which can solve both equilibrium and non-equilibrium problems. We analyse spectral functions of Bose-Hubbard model on the Bethe lattice in different regimes and show that our solver gives correct and physical results.