, while its metallic phase only occupies a narrow region in the phase diagram. To explore the enlargement of the metallic regime, we extend the ionic Hubbard model to two coupled layers and study the interplay of interlayer hybridization

V

and two types of intralayer staggered potentials

\mathrm{\ensuremath{\Delta}}

: one with the same (in-phase) and the other with a

\ensuremath{\pi}

-phase shift (antiphase) potential between layers. Our determinant quantum Monte Carlo (DQMC) simulations at lowest accessible temperatures demonstrate that the interaction-driven metallic phase between Mott and band insulators expands in the

\mathrm{\ensuremath{\Delta}}\text{\ensuremath{-}}V

phase diagram of bilayer IHM only for in-phase ionic potentials; while antiphase potential always induces an insulator with charge density order. This implies possible further extension of the ionic Hubbard model from the bilayer case here to a realistic three-dimensional model.

DOI: 10.1103/physrevb.93.165146

Related papers

Interaction-driven Band-insulator--to--Metal transition in bilayer ionic Hubbard model

Mi Jiang, T. C. Schulthess

The interaction-driven insulator-to-metal transition has been reported in the ionic Hubbard model (IHM) for moderate interaction $U$ , while its metallic phase only occupies a narrow region in the phase diagram. To explore the enlargement of the metallic regime, we extend the ionic Hubbard model to two coupled layers and study the interplay of interlayer hybridization $V$ and two types of intralayer staggered potentials $Δ$ : one with the same (in-phase) and the other with a $π$ -phase shift (anti-phase) potential between layers. Our determinant Quantum Monte Carlo (DQMC) simulations at lowest accessible temperatures demonstrate that the interaction-driven metallic phase between Mott and band insulators expands in the $Δ-V$ phase diagram of bilayer IHM only for in-phase ionic potentials; while anti-phase potential always induces an insulator with charge density order. This implies possible further extension of the ionic Hubbard model from the bilayer case here to realistic three-dimensional model.The interaction-driven insulator-to-metal transition has been reported in the ionic Hubbard model (IHM) for moderate interaction $U$, while its metallic phase only occupies a narrow region in the phase diagram. To explore the enlargement of the metallic regime, we extend the ionic Hubbard model to two coupled layers and study the interplay of interlayer hybridization $V$ and two types of intralayer staggered potentials $Δ$: one with the same (in-phase) and the other with a $π$-phase shift (anti-phase) potential between layers. Our determinant Quantum Monte Carlo (DQMC) simulations at lowest accessible temperatures demonstrate that the interaction-driven metallic phase between Mott and band insulators expands in the $Δ-V$ phase diagram of bilayer IHM only for in-phase ionic potentials; while anti-phase potential always induces an insulator with charge density order. This implies possible further extension of the ionic Hubbard model from the bilayer case here to realistic three-dimensional model.

arXiv:1604.03271Apr 12, 20161 citation

Quantum Monte Carlo Simulation of the two-dimensional ionic Hubbard model

B. Martinie

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.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.

arXiv:1002.0325Feb 1, 20100 citations

Metallic phase in the two-dimensional ionic Hubbard model

Karim Bouadim, Norman Paris, F. Hébert +2 more

We investigate the phases of the ionic Hubbard model in a two-dimensional square lattice using determinant quantum Monte Carlo. At half-filling, when the interaction strength or the staggered potential dominates, we find Mott and band insulators, respectively. When these two energies are of the same order, we find a metallic region. Charge and magnetic structure factors demonstrate the presence of antiferromagnetism only in the Mott region, although the externally imposed density modulation is present everywhere in the phase diagram. Away from half-filling, other insulating phases are found. Kinetic energy correlations do not give clear signals for the existence of a bond-ordered phase.We investigate the phases of the ionic Hubbard model in a two-dimensional square lattice using determinant quantum Monte Carlo. At half-filling, when the interaction strength or the staggered potential dominates, we find Mott and band insulators, respectively. When these two energies are of the same order, we find a metallic region. Charge and magnetic structure factors demonstrate the presence of antiferromagnetism only in the Mott region, although the externally imposed density modulation is present everywhere in the phase diagram. Away from half-filling, other insulating phases are found. Kinetic energy correlations do not give clear signals for the existence of a bond-ordered phase.

arXiv:0705.0755Physical Review BAug 9, 200770 citations

Quantum Monte Carlo Study of an Interaction-Driven Band-Insulator–to–Metal Transition

Norman Paris, Karim Bouadim, F. Hébert +2 more

We study the transitions from band insulator to metal to Mott insulator in the ionic Hubbard model on a two-dimensional square lattice using determinant quantum Monte Carlo. Evaluation of the temperature dependence of the conductivity demonstrates that the metallic region extends for a finite range of interaction values. The Mott phase at strong coupling is accompanied by antiferromagnetic order. Inclusion of these intersite correlations changes the phase diagram qualitatively compared to dynamical mean field theory.We study the transitions from band insulator to metal to Mott insulator in the ionic Hubbard model on a two-dimensional square lattice using determinant quantum Monte Carlo. Evaluation of the temperature dependence of the conductivity demonstrates that the metallic region extends for a finite range of interaction values. The Mott phase at strong coupling is accompanied by antiferromagnetic order. Inclusion of these intersite correlations changes the phase diagram qualitatively compared to dynamical mean field theory.

arXiv:cond-mat/0607707Physical Review LettersJan 25, 2007102 citations

Dc conductivity of two-dimensional ionic Hubbard model extracted by analytic continuation

Yiqun Huang, Wanpeng Han, Junsong Sun +2 more

The ionic Hubbard model in a two-dimensional square lattice is revisited using determinant quantum Monte Carlo (DQMC) and numerical analytic continuation. From the dc conductivity extracted by stochastic maximum entropy of the DQMC data, we find that the band insulator induced by a staggered potential [Formula: see text] goes through a metallic phase before entering a Mott insulator as the on-site Hubbard interaction U is increased. Moreover, we map out the phase diagram in the [Formula: see text] plane, and reveal a smaller metallic region than that found in the previous studies. Our results provide a more exact location of the phase boundaries, and are an important complement to those obtained by an approximate approach.The ionic Hubbard model in a two-dimensional square lattice is revisited using determinant quantum Monte Carlo (DQMC) and numerical analytic continuation. From the dc conductivity extracted by stochastic maximum entropy of the DQMC data, we find that the band insulator induced by a staggered potential [Formula: see text] goes through a metallic phase before entering a Mott insulator as the on-site Hubbard interaction U is increased. Moreover, we map out the phase diagram in the [Formula: see text] plane, and reveal a smaller metallic region than that found in the previous studies. Our results provide a more exact location of the phase boundaries, and are an important complement to those obtained by an approximate approach.

International Journal of Modern Physics CDec 27, 20230 citations

The phase diagram of the square lattice bilayer Hubbard model: a variational Monte Carlo study

Robert Rüger, Luca F. Tocchio, Roser Valentí +1 more

We investigate the phase diagram of the square lattice bilayer Hubbard model at half-filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of the interlayer hopping and on-site Coulomb repulsion U . With this study we resolve some discrepancies in previous calculations based on the dynamical mean-field theory, and we are able to determine the nature of the phase transitions between metal, Mott insulator and band insulator. In the magnetic case we find only two phases: an antiferromagnetic Mott insulator at small for any value of U and a band insulator at large . At large U values we approach the Heisenberg limit. The paramagnetic phase diagram shows at small a metal to Mott insulator transition at moderate U values and a Mott to band insulator transition at larger U values. We also observe a re-entrant Mott insulator to metal transition and metal to band insulator transition for increasing in the range of . Finally, we discuss the phase diagrams obtained in relation to findings from previous studies based on different many-body approaches.We investigate the phase diagram of the square lattice bilayer Hubbard model at half-filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of the interlayer hopping and on-site Coulomb repulsion U . With this study we resolve some discrepancies in previous calculations based on the dynamical mean-field theory, and we are able to determine the nature of the phase transitions between metal, Mott insulator and band insulator. In the magnetic case we find only two phases: an antiferromagnetic Mott insulator at small for any value of U and a band insulator at large . At large U values we approach the Heisenberg limit. The paramagnetic phase diagram shows at small a metal to Mott insulator transition at moderate U values and a Mott to band insulator transition at larger U values. We also observe a re-entrant Mott insulator to metal transition and metal to band insulator transition for increasing in the range of . Finally, we discuss the phase diagrams obtained in relation to findings from previous studies based on different many-body approaches.

arXiv:1311.6504New Journal of PhysicsMar 7, 201428 citations

The interaction-driven insulator-to-metal transition has been reported in the ionic Hubbard model (IHM) for moderate interaction

U

V

and two types of intralayer staggered potentials

\mathrm{\ensuremath{\Delta}}

: one with the same (in-phase) and the other with a

\ensuremath{\pi}

\mathrm{\ensuremath{\Delta}}\text{\ensuremath{-}}V

The interaction-driven insulator-to-metal transition has been reported in the ionic Hubbard model (IHM) for moderate interaction

U

V

and two types of intralayer staggered potentials

Δ

: one with the same (in-phase) and the other with a

π

-phase shift (anti-phase) potential between layers. Our determinant Quantum Monte Carlo (DQMC) simulations at lowest accessible temperatures demonstrate that the interaction-driven metallic phase between Mott and band insulators expands in the

Δ-V

phase diagram of bilayer IHM only for in-phase ionic potentials; while anti-phase potential always induces an insulator with charge density order. This implies possible further extension of the ionic Hubbard model from the bilayer case here to realistic three-dimensional model.