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

Mott transition in one dimension: Benchmarking dynamical cluster approaches

Matthias Balzer, Werner R. Hanke, Michael Potthoff

Published Jan 29, 2008• 36 citations
arXivHTMLPDF

Abstract

The variational cluster approach (VCA) is applied to the one-dimensional Hubbard model at zero temperature using clusters (chains) of up to ten sites with full diagonalization and the Lanczos method as cluster solver. Within the framework of the self-energy-functional theory (SFT), different cluster reference systems with and without bath degrees of freedom, in different topologies and with different sets of variational parameters, are considered. Static and one-particle dynamical quantities are calculated for half-filling as a function of $U$ as well as for fixed $U$ as a function of the chemical potential to study the interaction- and filling-dependent metal-insulator (Mott) transition. The recently developed $Q$-matrix technique is used to compute the SFT grand potential. For benchmarking purposes we compare the VCA results with exact results available from the Bethe ansatz, with essentially exact dynamical density-matrix renormalization-group data, with (cellular) dynamical mean-field theory and full diagonalization of isolated Hubbard chains. Several issues are discussed including convergence of the results with cluster size, the ability of cluster approaches to access the critical regime of the Mott transition, efficiency in the optimization of correlated-site vs bath-site parameters, and of multidimensional parameter optimization. We also study the role of bath sites for the description of excitation properties and as charge reservoirs for the description of filling dependencies. The VCA turns out to be a computationally cheap method which is competitive with established cluster approaches.

DOI
10.1103/physrevb.77.045133

Related papers

Variational cluster approach to thermodynamic properties of interacting fermions at finite temperatures: A case study of the two-dimensional single-band Hubbard model at half filling

Kazuhiro Seki, Tomonori Shirakawa, Seiji Yunoki

We formulate a finite-temperature scheme of the variational cluster approximation (VCA) particularly suitable for an exact-diagonalization cluster solver. Based on the analytical properties of the single-particle Green's function matrices, we explicitly show the branch-cut structure of logarithm of the complex determinant functions appearing in the self-energy-functional theory (SFT) and whereby construct an efficient scheme for the finite-temperature VCA. We also derive the explicit formulas for entropy and specific heat within the framework of the SFT. We first apply the method to explore the antiferromagnetic order in a half-filled Hubbard model by calculating the entropy, specific heat, and single-particle excitation spectrum for different values of on-site Coulomb repulsion $U$ and temperature $T$. We also calculate the $T$ dependence of the single-particle excitation spectrum in the strong coupling region, and discuss the overall similarities to and the fine differences from the spectrum obtained by the spin-density-wave mean-field theory at low temperatures and the Hubbard-I approximation at high temperatures. Moreover, we show a necessary and sufficient condition for the third law of thermodynamics in the SFT. On the basis of the thermodynamic properties, such as the entropy and the double occupancy, calculated via the $T$ and/or $U$ derivative of the grand potential, we obtain a crossover diagram in the $(U,T)$ plane, which separates a Slater-type insulator and a Mott-type insulator. Next, we demonstrate the finite-temperature scheme in the cluster-dynamical-impurity approximation (CDIA), i.e., the VCA with noninteracting bath orbitals attached to each cluster, and study the paramagnetic Mott metal-insulator transition in the half-filled Hubbard model. Formulating the finite-temperature CDIA, we first address a subtle issue regarding the treatment of the artificially introduced bath degrees of freedom which are absent in the originally considered Hubbard model. We then apply the finite-temperature CDIA to calculate the finite-temperature phase diagram in the $(U,T)$ plane. Metallic, insulating, coexistence, and crossover regions are distinguished from the bath-cluster hybridization-variational-parameter dependence of the grand-potential functional. We find that the Mott transition at low temperatures is discontinuous, and the coexistence region of the metallic and insulating states persists down to zero temperature. The result obtained here by the finite-temperature CDIA is complementary to the previously reported zero-temperature CDIA phase diagram.

arXiv:1808.03807•Physical review. B./Physical review. B•Nov 7, 2018•26 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

Cellular dynamical mean-field theory for the one-dimensional extended Hubbard model

C. J. Bolech, S. S. Kancharla, Gabriel Kotliar

We explore the use of exact diagonalization methods for solving the self-consistent equations of the cellular dynamical mean field theory for the one-dimensional regular and extended Hubbard models. We investigate the nature of the Mott transition and convergence of the method as a function of cluster size as well as the optimal allocation of computational resources between bath and ``cluster-impurity'' sites, with a view to develop a renormalization group method in higher dimensions. We assess the performance of the method by comparing results for the Green's functions in both the spin density wave and charge density wave phases with accurate density matrix renormalization group (DMRG) calculations.

arXiv:cond-mat/0206166•Physical review. B, Condensed matter•Feb 28, 2003•66 citations

Variational cluster approach to ferromagnetism in infinite dimensions and in one-dimensional chains

Matthias Balzer, Michael Potthoff

The variational cluster approach (VCA) is applied to study spontaneous ferromagnetism in the Hubbard model at zero temperature. We discuss several technical improvements of the numerical implementation of the VCA which become necessary for studies of a ferromagnetically ordered phase, e.g., more accurate techniques to evaluate the variational ground-state energy, improved local as well as global algorithms to find stationary points, and different methods to locate the magnetic phase transition. Using the single-site VCA, i.e., the dynamical impurity approximation (DIA), the ferromagnetic phase diagram of the model in infinite dimensions is worked out. The results are compared with previous dynamical mean-field studies for benchmarking purposes. The DIA results provide a unified picture of ferromagnetism in the infinite-dimensional model by interlinking different parameter regimes that are governed by different mechanisms for ferromagnetic order. Using the DIA and the VCA, we then study ferromagnetism in one-dimensional Hubbard chains with nearest- and next-nearest-neighbor hopping ${t}_{2}$. In comparison with previous results from the density-matrix renormalization group, the phase diagram is mapped out as a function of the Hubbard-$U$, the electron filling, and ${t}_{2}$. The stability of the ferromagnetic ground state against local and short-range nonlocal quantum fluctuations is discussed.

arXiv:1007.2517•Physical Review B•Nov 30, 2010•21 citations

Cluster-dynamical mean-field theory of the density-driven Mott transition in the one-dimensional Hubbard model

Massimo Capone, Marcello Civelli, S. S. Kancharla +2 more

The one-dimensional Hubbard model is investigated by means of two different cluster schemes suited to introduce short-range spatial correlations beyond the single-site dynamical mean-field theory, namely, the cellular dynamical mean-field theory, which does not impose the lattice symmetries, and its periodized version in which translational symmetry is recovered. It is shown that both cluster schemes are able to describe with extreme accuracy the evolution of the density as a function of the chemical potential from the Mott insulator to the metallic state. Using exact diagonalization to solve the cluster-impurity model, we discuss the role of the truncation of the Hilbert space of the bath, and propose an algorithm that gives higher weights to the low-frequency hybridization matrix elements and improves the speed of the convergence of the algorithm.

arXiv:cond-mat/0401060•Physical Review B•May 18, 2004•88 citations

Mott transition and pseudogap of the square-lattice Hubbard model: Results from center-focused cellular dynamical mean-field theory

Michael Meixner, Henri Menke, Marcel Klett +4 more

The recently proposed center-focused post-processing procedure [Phys. Rev. Res. 2, 033476 (2020)] of cellular dynamical mean-field theory suggests that central sites of large impurity clusters are closer to the exact solution of the Hubbard model than the edge sites. In this paper, we systematically investigate results in the spirit of this center-focused scheme for several cluster sizes up to 8×8 in and out of particle-hole symmetry. First we analyze the metal-insulator crossovers and transitions of the half-filled Hubbard model on a simple square lattice. We find that the critical interaction of the crossover is reduced with increasing cluster sizes and the critical temperature abruptly drops for the 4×4 cluster. Second, for this cluster size, we apply the center-focused scheme to a system with more realistic tight-binding parameters, investigating its pseudogap regime as a function of temperature and doping, where we find doping dependent metal-insulator crossovers, Lifshitz transitions and a strongly renormalized Fermi-liquid regime. Additionally to diagnosing the real space origin of the suppressed antinodal spectral weight in the pseudogap regime, we can infer hints towards underlying charge ordering tendencies.

arXiv:2310.17302•SciPost Physics•Feb 27, 2024•14 citations
The variational cluster approach (VCA) is applied to the one-dimensional Hubbard model at zero temperature using clusters (chains) of up to ten sites with full diagonalization and the Lanczos method as cluster solver. Within the framework of the self-energy-functional theory (SFT), different cluster reference systems with and without bath degrees of freedom, in different topologies and with different sets of variational parameters, are considered. Static and one-particle dynamical quantities are calculated for half-filling as a function of UUU as well as for fixed UUU as a function of the chemical potential to study the interaction- and filling-dependent metal-insulator (Mott) transition. The recently developed QQQ-matrix technique is used to compute the SFT grand potential. For benchmarking purposes we compare the VCA results with exact results available from the Bethe ansatz, with essentially exact dynamical density-matrix renormalization-group data, with (cellular) dynamical mean-field theory and full diagonalization of isolated Hubbard chains. Several issues are discussed including convergence of the results with cluster size, the ability of cluster approaches to access the critical regime of the Mott transition, efficiency in the optimization of correlated-site vs bath-site parameters, and of multidimensional parameter optimization. We also study the role of bath sites for the description of excitation properties and as charge reservoirs for the description of filling dependencies. The VCA turns out to be a computationally cheap method which is competitive with established cluster approaches.
We formulate a finite-temperature scheme of the variational cluster approximation (VCA) particularly suitable for an exact-diagonalization cluster solver. Based on the analytical properties of the single-particle Green's function matrices, we explicitly show the branch-cut structure of logarithm of the complex determinant functions appearing in the self-energy-functional theory (SFT) and whereby construct an efficient scheme for the finite-temperature VCA. We also derive the explicit formulas for entropy and specific heat within the framework of the SFT. We first apply the method to explore the antiferromagnetic order in a half-filled Hubbard model by calculating the entropy, specific heat, and single-particle excitation spectrum for different values of on-site Coulomb repulsion UUU and temperature TTT. We also calculate the TTT dependence of the single-particle excitation spectrum in the strong coupling region, and discuss the overall similarities to and the fine differences from the spectrum obtained by the spin-density-wave mean-field theory at low temperatures and the Hubbard-I approximation at high temperatures. Moreover, we show a necessary and sufficient condition for the third law of thermodynamics in the SFT. On the basis of the thermodynamic properties, such as the entropy and the double occupancy, calculated via the TTT and/or UUU derivative of the grand potential, we obtain a crossover diagram in the (U,T)(U,T)(U,T) plane, which separates a Slater-type insulator and a Mott-type insulator. Next, we demonstrate the finite-temperature scheme in the cluster-dynamical-impurity approximation (CDIA), i.e., the VCA with noninteracting bath orbitals attached to each cluster, and study the paramagnetic Mott metal-insulator transition in the half-filled Hubbard model. Formulating the finite-temperature CDIA, we first address a subtle issue regarding the treatment of the artificially introduced bath degrees of freedom which are absent in the originally considered Hubbard model. We then apply the finite-temperature CDIA to calculate the finite-temperature phase diagram in the (U,T)(U,T)(U,T) plane. Metallic, insulating, coexistence, and crossover regions are distinguished from the bath-cluster hybridization-variational-parameter dependence of the grand-potential functional. We find that the Mott transition at low temperatures is discontinuous, and the coexistence region of the metallic and insulating states persists down to zero temperature. The result obtained here by the finite-temperature CDIA is complementary to the previously reported zero-temperature CDIA phase diagram.
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.
We explore the use of exact diagonalization methods for solving the self-consistent equations of the cellular dynamical mean field theory for the one-dimensional regular and extended Hubbard models. We investigate the nature of the Mott transition and convergence of the method as a function of cluster size as well as the optimal allocation of computational resources between bath and ``cluster-impurity'' sites, with a view to develop a renormalization group method in higher dimensions. We assess the performance of the method by comparing results for the Green's functions in both the spin density wave and charge density wave phases with accurate density matrix renormalization group (DMRG) calculations.
The variational cluster approach (VCA) is applied to study spontaneous ferromagnetism in the Hubbard model at zero temperature. We discuss several technical improvements of the numerical implementation of the VCA which become necessary for studies of a ferromagnetically ordered phase, e.g., more accurate techniques to evaluate the variational ground-state energy, improved local as well as global algorithms to find stationary points, and different methods to locate the magnetic phase transition. Using the single-site VCA, i.e., the dynamical impurity approximation (DIA), the ferromagnetic phase diagram of the model in infinite dimensions is worked out. The results are compared with previous dynamical mean-field studies for benchmarking purposes. The DIA results provide a unified picture of ferromagnetism in the infinite-dimensional model by interlinking different parameter regimes that are governed by different mechanisms for ferromagnetic order. Using the DIA and the VCA, we then study ferromagnetism in one-dimensional Hubbard chains with nearest- and next-nearest-neighbor hopping t2{t}_{2}t2​. In comparison with previous results from the density-matrix renormalization group, the phase diagram is mapped out as a function of the Hubbard-UUU, the electron filling, and t2{t}_{2}t2​. The stability of the ferromagnetic ground state against local and short-range nonlocal quantum fluctuations is discussed.
The one-dimensional Hubbard model is investigated by means of two different cluster schemes suited to introduce short-range spatial correlations beyond the single-site dynamical mean-field theory, namely, the cellular dynamical mean-field theory, which does not impose the lattice symmetries, and its periodized version in which translational symmetry is recovered. It is shown that both cluster schemes are able to describe with extreme accuracy the evolution of the density as a function of the chemical potential from the Mott insulator to the metallic state. Using exact diagonalization to solve the cluster-impurity model, we discuss the role of the truncation of the Hilbert space of the bath, and propose an algorithm that gives higher weights to the low-frequency hybridization matrix elements and improves the speed of the convergence of the algorithm.
The recently proposed center-focused post-processing procedure [Phys. Rev. Res. 2, 033476 (2020)] of cellular dynamical mean-field theory suggests that central sites of large impurity clusters are closer to the exact solution of the Hubbard model than the edge sites. In this paper, we systematically investigate results in the spirit of this center-focused scheme for several cluster sizes up to 8×8 in and out of particle-hole symmetry. First we analyze the metal-insulator crossovers and transitions of the half-filled Hubbard model on a simple square lattice. We find that the critical interaction of the crossover is reduced with increasing cluster sizes and the critical temperature abruptly drops for the 4×4 cluster. Second, for this cluster size, we apply the center-focused scheme to a system with more realistic tight-binding parameters, investigating its pseudogap regime as a function of temperature and doping, where we find doping dependent metal-insulator crossovers, Lifshitz transitions and a strongly renormalized Fermi-liquid regime. Additionally to diagnosing the real space origin of the suppressed antinodal spectral weight in the pseudogap regime, we can infer hints towards underlying charge ordering tendencies.