In this thesis, several aspects of correlation effects in band insulators that were hitherto not addressed, both in equilibrium and non-equilibrium, are studied using the Ionic Hubbard model (IHM) at half-filling and Dynamical mean-field theory (DMFT). The IHM is an extension of the Hubbard model on a bipartite lattice with a staggered on-site “ionic potential Δ added in. The Hubbard model by itself is a tight-binding model of electrons hopping between (orbitals localized on) lattice sites, with amplitudes t and t’ for first neighbor and second neighbor hopping respectively, together with a coulomb repulsion energy cost of U for double occupancy (i.e., with both up and down spin electrons) at any lattice site. The studies in this thesis address strong correlation effects in this model using DMFT techniques which map the lattice problem to a self consistently embedded quantum impurity problem. The impurity problem itself is solved using state of the art continuous time quantum monte-carlo (CTQMC) techniques (Ch.s 3 and 4) as well as the somewhat simpler Iterated Perturbation Theory (IPT) (Ch.s 3, 5 and 6).
Chapter 3 of the thesis is devoted to a detailed analysis of the properties of the half-filled t-U-Δ IHM (i.e., with t’=0) on a Bethe lattice of infinite connectivity. It is shown conclusively here that for a finite Δ and at zero temperature (T=0), the IHM has a first order transition from the Band Insulator (BI) phase to an antiferromagnetic insulating (AFI) phase at a threshold U = UAF (which increases with Δ) such that the BI to AFI transition preempts the transition from the BI to the correlation induced paramagnetic metallic (PM) phase. Evidence is also presented for a quantum transition to a sliver of half-metal (HM) phase just after the AF order turns on, followed by a return to the AFI phase on further increasing U. With increasing T, the AF order is lost via a first-order transition for weak to intermediate U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, the transition from the AFI to the paramagnetic Mott insulating phase (adiabatically connected to the BI phase) is of second order. As T increases, the range of U over which the AF order is stable shrinks, collapsing eventually to a line of tricritical points that separates the surfaces of first- and second-order phase transitions in the T-U-Δ space.
In Chapter 4 of the thesis the t-t’-U-Δ IHM on a square lattice is studied using DMFT+CTQMC. It is shown here that the presence of t’ frustrates the AF order for weak to intermediate correlations, resulting in the stabilization of the correlation induced PM phase over a broad region of the parameter space, although the AFI is always the stable phase for large enough U. In the intermediate coupling regimes, two other interesting phases show up as a result of the competition among U, Δ and t’ : a Ferrimagnetic metallic (FM) phase which has non-zero values of the uniform as well as staggered magnetization; and even more interesting, an…
We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity using dynamical mean-field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered ionic potential $\mathrm{\ensuremath{\Delta}}$ and the on-site Hubbard
U. We find that for a finite $\mathrm{\ensuremath{\Delta}}$ and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold
U=UAF via a first-order phase transition. For
U smaller than
UAF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing
U. We show that the results obtained within both methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature ${T}_{AF}(U,\mathrm{\ensuremath{\Delta}})$ [or, equivalently, on decreasing
U below ${U}_{AF}(T,\mathrm{\ensuremath{\Delta}})$], within both methods, for weak to intermediate values of
U/t. In the strongly correlated regime, where the effective low-energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite temperature
T,DMFT+CTQMC shows a second phase transition (not seen within
DMFT+IPT) on increasing
U beyond
UAF. At
UN>UAF, when the Neel temperature
TN for the effective Heisenberg model becomes lower than
T, the AFM order is lost via a second-order transition. For $U\ensuremath{\gg}\mathrm{\ensuremath{\Delta}}, {T}_{N}\ensuremath{\sim}{t}^{2}/U(1\ensuremath{-}{x}^{2})$, where $x=2\mathrm{\ensuremath{\Delta}}/U$ and thus
TN increases with increase in $\mathrm{\ensuremath{\Delta}}/U$. In the three-dimensional parameter space of $(U/t,\phantom{\rule{0.28em}{0ex}}T/t,\phantom{\rule{4.pt}{0ex}}\text{and}\phantom{\rule{4.pt}{0ex}}\mathrm{\ensuremath{\Delta}}/t)$, as
T increases, the surface of first-order transition at ${U}_{AF}(T,\mathrm{\ensuremath{\Delta}})$ and that of the second-order transition at ${U}_{N}(T,\mathrm{\ensuremath{\Delta}})$ approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.
We investigate the ionic Hubbard model (IHM) at half-filling in the limit of strong correlations and large ionic potential. The low-energy effective Hamiltonian in this limit, obtained by a similarity transformation, is a modified $t\ensuremath{-}J$ model with effective second-neighbor hopping terms. We explore the possibilities of
d-wave pairing and extended
s-wave pairing superconducting (SC) phases on a two-dimensional square lattice at zero temperature within a Gutzwiller projected renormalized mean-field theory. In the sector of solutions that forbid spin-ordering, the system shows a finite nonzero
d-wave as well as extended
s-wave pairing amplitude for $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}U\ensuremath{\gg}t$. The width of the superconducting phase in $U\ensuremath{-}\mathrm{\ensuremath{\Delta}}$ regime shrinks with increase in
U and $\mathrm{\ensuremath{\Delta}}$, though the extended
s-wave pairing phase is higher in energy than the
d-wave pairing superconducting phase. But in a spin-resolved renormalized mean-field calculation, which allows for an antiferromagnetic (AF) order along with the
d-wave or extended
s-wave pairing, the SC phase is no longer viable and the system shows a direct transition from an AF ordered phase to a paramagnetic band insulator. Except for a thin sliver of a half-metallic AF phase close to the AF transition point, most of the AF ordered phase is a Mott insulator. We benchmarked the AF Mott insulator to band insulator transition within the Gutzwiller projected renormalized mean-field theory against the dynamical mean-field theory solved using continuous time quantum Monte Carlo. Our work suggests that the ground-state phase diagram of the IHM at half-filling in the limit of extreme correlations does not have any SC phase. The SC phase seen in the paramagnetic sector is a metastable phase, being higher in energy than the AF Mott insulator phase.
The Quantum Monte Carlo simulations of the ionic Hubbard model on a two-dimensional square lattice at half filling were performed. The method based on the direct-space, proposed by Suzuki and al., Hirsch and al., was used. Cycles of increasing and decreasing values of the Coulomb interaction $ U $ were performed for fixed temperature ($ kT=0.01 $). Results indicate that, at low temperature, the two insulator phases are separated by a metallic phase for weak to intermediate values of the staggered potential $ Δ$. For large Coulomb repulsion the system is in a Mott insulator with an antiferromagnetism order. On increasing and decreasing the Coulomb interaction $ U $ the metal-Mott insulator transition shows an hysteresis phenomenon while the metal-band insulator transition is continue. For large $ Δ$ it seems that the metallic region shrinks to a single metallic point. However, the band insulator to the Mott insulator transition is not direct for the studied model. A phase diagram is drawn for the temperature $ kT=0.01 $. For $ Δ=0.5 $ cycles of increasing and decreasing temperature were programmed for different values of the Coulomb interaction $ U $ . A behaviour change appears for $ U\simeq 1.75 $. This suggests that a crossover line divides the metallic region of the phase diagram.
Strongly correlated electron systems are one of the most fascinating problems in current physics. The strong electron-electron interaction in these materials leads to the emergence of nontrivial elementary excitations (quasiparticles, QPs) above the ground state ranging from fractional spins in quasi-one-dimensional materials to magnetic monopole in the pyrochlore lattice. The condensation of these quasiparticles upon changing some external parameters may stabilize new exotic states of matter. Experimental measurements such as inelastic neutron scattering provide us with valuable information about the excitation spectrum of such systems which require microscopic models to be described. This thesis is devoted to a detailed analysis of the excitation spectrum and of the quantum phase transitions in the one-dimensional (1D) ionic Hubbard model (IHM). The IHM consists of a nearest-neighbor (n.n.) hopping, onsite Hubbard interaction, and an ionic (staggered) potential separating the odd and even sites energetically. The model exhibits two continuous phase transitions on increasing the Hubbard interaction identified by a low-energy effective field theory and confirmed by a rigorous density matrix renormalization group (DMRG) analysis after several attempts. The first transition occurs from band insulator (BI) phase to the 2-fold degenerate spontaneously dimerized insulator (SDI) phase. The transition is in the Ising universality class as is plausible from symmetry considerations. The SDI phase becomes unstable towards a quasi-long-range order Mott insulator (MI) phase at a second transition point resembling the Kosterlitz-Thouless (KT) transition in the frustrated Heisenberg chain. We employ continuous unitary transformations (CUT) to systematically map the IHM to effective Hamiltonians (almost) conserving the number of QPs in the system. Using an analysis in the BI regime where electrons and holes define QPs, the low-energy excitation spectrum of the model is quantitatively determined in the BI phase almost up to the first transition point. The transition from the BI to the SDI phase is signaled by the vanishing of an S=0 exciton mode at the total momentum K=\pi. The condensation of these excitons beyond the first transition point is described by a BCS-type-theory showing the stabilization of the SDI phase. The mean-field solution indicates no second phase transition to the quasi-long-range order MI phase. This is interpreted as the effect of strong quantum fluctuations in one dimension. We consider the IHM in the dimer limit where the uniform chain is separated into independent dimers. The different phases of the IHM are studied by increasing the interdimer hopping and reaching the uniform limit. This dimer limit satisfactorily produces the excitation spectrum of the BI phase confirming the vanishing of an S=0 exciton mode at the first transition point. It is found that the SDI-to-MI transition takes place by softening of a magnetic S=1 excitation, i.e., a triplon. We report rigorous results for the gapless triplon dispersion in the MI phase and discuss the binding effects in the 2-triplon sector.
We investigate the effects of Coulomb disorder on the half-filled ionic Hubbard model within the framework of dynamical mean-field theory (DMFT), using the equation-of-motion method as the impurity solver, combined with typical medium theory (TMT). In contrast to conventional disordered Hubbard models, the ionic Hubbard model contains a staggered ionic potential that introduces additional competition between band-insulating and correlation-driven insulating tendencies. We analyze how this ionic potential modifies the electronic phases induced by Coulomb disorder through the arithmetic and geometric averages of the local density of states. In the absence of disorder, the system exhibits metallic (M), Mott-insulating (MI), and band-insulating (BI) phases. When Coulomb disorder is introduced, Anderson-insulating (AI) behavior and localized states (LS) inside the Mott gap are also observed. In particular, the LS regime appears between the metallic and Mott-insulating phases at small ionic potentials, whereas the AI regime develops for sufficiently large values of the Coulomb interaction and ionic potential. Our results show that the interplay between ionic potential, electron correlation, and disorder substantially influences the localization behavior and insulating characteristics of the ionic Hubbard model.
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.