The ionic Hubbard model is a paradigmatic setup for studying the competition between band and Mott insulating behavior. Within the variationally exact in infinite dimensions Gutzwiller approximation, we derive a compact analytic expression for the phase boundary between the Mott and band insulators. While the method reproduces the expected band–Mott insulator phenomenology, it does not capture the correlated metallic state at finite staggered potential found, for example, in dynamical mean-field theory. This absence highlights that the metallic phase originates from incoherent Hubbard-band physics rather than Fermi-liquid behavior well captured by the Gutzwiller approximation. Our formulation establishes a concise variational framework for the ionic Hubbard model, with natural extensions to nonequilibrium setups and spin-exchange dynamics. Abstract Published by the Jagiellonian University 2026 authors
We investigate paramagnetic metal-insulator transitions in the infinite-dimensional ionic Hubbard model at finite temperatures. By means of the dynamical mean-field theory with an impurity solver of the continuous-time quantum Monte Carlo method, we show that an increase in the interaction strength brings about a crossover from a band insulating phase to a metallic one, followed by a first-order transition to a Mott insulating phase. The first-order transition turns into a crossover above a certain critical temperature, which becomes higher as the staggered lattice potential is increased. Further, analysis of the temperature dependence of the energy density discloses that the intermediate metallic phase is a Fermi liquid. It is also found that the metallic phase is stable against strong staggered potentials even at very low temperatures.
The main flaw of the well-known Brinkman-Rice description, obtained through the Gutzwiller approximation, of the paramagnetic Mott transition in the Hubbard model is in neglecting high-energy virtual processes that generate, for instance, the antiferromagnetic exchange $J\ensuremath{\sim}{t}^{2}/U$. Here, we propose a way to capture those processes by combining the Brinkman-Rice approach with a variational Schrieffer-Wolff transformation, and apply this method to study the single-band metal-to-insulator transition in a Bethe lattice with infinite coordination number, where the Gutzwiller approximation becomes exact. We indeed find for the Mott transition a description very close to the real one provided by the dynamical mean-field theory, an encouraging result in view of possible applications to more involved models.
Introduction In this chapter we study metal–insulator transitions which occur due to correlations between the charge carriers. This quantum phase transitions results from the competition between kinetic energy and Coulomb repulsion where the latter has the tendency to localise the electrons. The simplest model which describes this competition is the one-band Hubbard model with its two energy scales, the bandwidth W and the Coulomb repulsion U between electrons of opposite spins on the same site. For a fixed density of electrons, n = 1; as the ratio U/W increases, we expect to find a zero-temperature phase transition from a metallic state to an insulator where the electrons are frozen at the lattice sites. The method employed here to investigate this transition is the Gutzwiller approach (Gutzwiller, 1965). This is a variational method which played an important role in providing a framework for the ideas of Mott on the problem of the metal–insulator transition (Mott, 1974). On the metallic side of the transition, this approach describes a metal with renormalised Fermi liquid parameters, essentially with enhanced effective mass and uniform susceptibility. This is the so-called highly correlated electron gas (Brinkman and Rice, 1970) for which the metal oxide V 2 O 3 is the prototype system. In the insulating side, Gutzwiller's approximation yields a noninteracting system of localised electrons. This insulating phase is incompressible , i.e. has zero compressibility and infinite zero temperature susceptibility. The latter is not a consequence of interactions but simply due to the divergence of the Curie susceptibility, χ ∝ 1 /T , of non-interacting moments at T = 0. Then, in Gutzwiller's approach the transition is from a paramagnetic metal to a paramagnetic insulator. We will continue using the designation Mott transition to describe this type of phase transition where the insulating phase has no long-range magnetic order. In real systems, most frequently, localisation due to interactions is accompanied by the appearance of antiferromagnetic order.
We investigate the local moment formation in the half-filled
SU(N) Hubbard model under a staggered ionic potential. As the Hubbard
U increases, the charge fluctuations are suppressed and eventually frozen when
U is above a critical value
Uc, marking the development of well-defined local moment with integer
m fermions on the A sublattice and $(N\ensuremath{-}m)$ fermions on the B sublattice, respectively. We obtain an analytical solution for
Uc for the paramagnetic ground state within the variational Gutzwiller approximation and renormalized mean field theory. For large
N,
Uc is found to depend on
N linearly with fixed
m/N, but sublinearly with fixed
m. The local moment formation is accompanied by a peculiar phase transition from the band insulator to the Mott insulator, where the ionic potential and quasiparticle weight are renormalized to zero simultaneously. Inside the Mott phase, the low-energy physics is described by the
SU(N) Heisenberg model with conjugate representations, which is widely studied in the literature.
We review the recently proposed extension of the Gutzwiller approximation, M. Schiro' and M. Fabrizio, Phys. Rev. Lett. 105, 076401 (2010), designed to describe the out-of-equilibrium time-evolution of a Gutzwiller-type variational wave function for correlated electrons. The method, which is strictly variational in the limit of infinite lattice-coordination, is quite general and flexible, and it is applicable to generic non-equilibrium conditions, even far beyond the linear response regime. As an application, we discuss the quench dynamics of a single-band Hubbard model at half-filling, where the method predicts a dynamical phase transition above a critical quench that resembles the sharp crossover observed by time-dependent dynamical mean field theory. We next show that one can actually define in some cases a multi-configurational wave function combination of a whole set of mutually orthogonal Gutzwiller wave functions. The Hamiltonian projected in that subspace can be exactly evaluated and is equivalent to a model of auxiliary spins coupled to non-interacting electrons, closely related to the slave-spin theories for correlated electron models. The Gutzwiller approximation turns out to be nothing but the mean-field approximation applied to that spin-fermion model, which displays, for any number of bands and integer fillings, a spontaneous
Z2 symmetry breaking that can be identified as the Mott insulator-to-metal transition.
A detailed study of electronic phase transitions in the ionic Hubbard model at half filling is presented. Within the dynamical mean field approximation a series of transitions from the band insulator via a metallic state to a Mott--Hubbard insulating phase is found at intermediate values of the one-body potential $\ensuremath{\Delta}$ with increasing the Coulomb interaction
U. We obtain a critical region in which the metallic phase disappears and a coexistence phase between the band and the Mott insulating state sets in. Our results are consistent with those obtained at low dimensions, thus they provide a concrete description for the charge degrees of freedom of the ionic Hubbard model.