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

Journal of the Korean Physical Society

Exact diagonalization and quantum Monte Carlo study of an ionic Hubbard model in two dimensions

Jongweon Cho, Ji-Woo Lee

Published Mar 1, 2017• 2 citations
Publisher

Abstract

We study quantum phase transitions of an ionic Hubbard model in two dimensions. The ionic Hubbard model explains the quantum states of strongly correlated electrons under the influence of checkerboard-type alternating chemical potentials. For a given amplitude of the alternating potentials Δ, we obtain quantum ground states as we tune the local repulsive energy U between a spin-up electron and a spin-down electron by using an exact diagonalization method of a modified Lanczos algorithm. The system undergoes a quantum phase transition from a band insulator to a Mott insulator as U increases at half-filling. We find the signature of a quantum phase transition by investigating the behavior of ground-state energies and that of double occupancies for the size of L × L = 4 × 4, which was the largest possible lattice in this work. We compare our results with those of quantum Monte Carlo simulations employing the Hirsch-Fye algorithm.

DOI
10.3938/jkps.70.494

Related papers

Quantum Monte Carlo study on an ionic Hubbard model in one dimension

Sung Moon Kim, Ji-Woo Lee

We study an ionic Hubbard model in one dimension by using a quantum Monte Carlo method. The ionic Hubbard model we study is a strongly correlated model for electrons in alternating potentials with a lattice period of 2. By modifying the Hirsch-Fye algorithm, we incorporate the local variation of the potential in the weight of the partition function at each imaginary-time segment. For a given alternating potential of strength Δ = 1.0, as we change the local Coulomb repulsion, U , between spin-up and spin-down electrons, we calculate various physical quantities, including kinetic energies, correlation energies, and band energies. We try to find the signature of the phase transitions from the behavior of the kinetic energy and that of double occupancies at finite temperatures.

Journal of the Korean Physical Society•Feb 1, 2012•2 citations

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.

arXiv:1002.0325•Feb 1, 2010•0 citations

Quantum Monte Carlo study of the one-dimensional ionic Hubbard model

Tim Wilkens, Richard M. Martin

Quantum Monte Carlo methods are used to study a quantum phase transition in a one-dimensional Hubbard model with on-site interaction U and a staggered ionic potential $(\ensuremath{\Delta}).$ Using recently formulated methods, the electronic polarization and localization are determined directly from the correlated ground state wave function and compared to results of previous work using exact diagonalization and Hartree-Fock. We find a transition from a symmetric band insulator to a broken-symmetry bond ordered (BO) phase as the ratio of $U/\ensuremath{\Delta}$ is increased, with a metallic point at the transition. Since it is known that at $\ensuremath{\Delta}=0$ the usual Hubbard model is a Mott insulator (MI) with no long-range order, we have searched for a second transition to this state by (i) increasing U at fixed $\ensuremath{\Delta}$ and (ii) decreasing $\ensuremath{\Delta}$ at fixed U. We find no transition from the BO to MI state, and we propose that the symmetric MI state exists only at $\ensuremath{\Delta}\ensuremath{\equiv}0.$

arXiv:cond-mat/0007472•Physical review. B, Condensed matter•May 16, 2001•70 citations

Nature of the insulating phases in the half-filled ionic Hubbard model

A. P. Kampf, Michael Sekania, G. I. Japaridze +1 more

We investigate the ground-state phase diagram of the one-dimensional 'ionic' Hubbard model with an alternating periodic potential at half-filling by numerical diagonalization of finite systems with the Lanczos and density matrix renormalization group methods. We identify an insulator–insulator phase transition from a band to a correlated insulator with simultaneous charge and bond-charge order. The transition point is characterized by the vanishing of the optical excitation gap while simultaneously the charge and spin gaps remain finite and equal. Indications for a possible second transition into a Mott-insulator phase are discussed.

arXiv:cond-mat/0304697•Journal of Physics Condensed Matter•Aug 15, 2003•76 citations

Ground-state phase diagram of the one-dimensional Hubbard model with an alternating chemical potential

Hiromi Otsuka, Masaaki Nakamura

We investigate the ground-state phase diagram of the one-dimensional half-filled Hubbard model with an alternating potential---a model for the charge-transfer organic materials and the ferroelectric perovskites. We numerically determine the global phase diagram of this model using the level-crossing and the phenomenological renormalization-group methods based on the exact diagonalization calculations. Our results support the mechanism of the double phase transitions between Mott and band insulators pointed out by Fabrizio, Gogolin, and Nersesyan [Phys. Rev. Lett. 83, 2014 (1999)]: We confirm the existence of the spontaneously dimerized phase as an intermediate state. Further we provide numerical evidence to check the criticalities on the phase boundaries. Especially, we perform the finite-size-scaling analysis of the excitation gap to show the two-dimensional Ising transition in the charge part. On the other hand, we confirm that the dimerized phase survives in the strong-coupling limit, which is one of the resultants of competition between the ionicity and correlation effects.

arXiv:cond-mat/0403630•Physical Review B•Apr 8, 2005•34 citations

Exact Diagonalization Study of Mott Transition in the Hubbard Model on an Anisotropic Triangular Lattice

Takashi Koretsune, Yukitoshi Motome, Akira Furusaki

We study Mott transition in the two-dimensional Hubbard model on an anisotropic triangular lattice. We use the Lanczos exact diagonalization of finite-size clusters up to eighteen sites, and calculate Drude weight, charge gap, double occupancy and spin structure factor. We average these physical quantities over twisted boundary conditions in order to reduce finite-size effects. We find a signature of the Mott transition in the dependence of the Drude weight and/or charge gap on the system size. We also examine the possibility of antiferromagnetic order from the spin structure factor. Combining these information, we propose a ground-state phase diagram which has a nonmagnetic insulating phase between a metallic phase and an insulating phase with antiferromagnetic order. Finally, we compare our results with those reported in the previous theoretical studies, and discuss the possibility of an unconventional insulating state.

arXiv:cond-mat/0703111•Journal of the Physical Society of Japan•Jul 9, 2007•75 citations
We study quantum phase transitions of an ionic Hubbard model in two dimensions. The ionic Hubbard model explains the quantum states of strongly correlated electrons under the influence of checkerboard-type alternating chemical potentials. For a given amplitude of the alternating potentials Δ, we obtain quantum ground states as we tune the local repulsive energy U between a spin-up electron and a spin-down electron by using an exact diagonalization method of a modified Lanczos algorithm. The system undergoes a quantum phase transition from a band insulator to a Mott insulator as U increases at half-filling. We find the signature of a quantum phase transition by investigating the behavior of ground-state energies and that of double occupancies for the size of L × L = 4 × 4, which was the largest possible lattice in this work. We compare our results with those of quantum Monte Carlo simulations employing the Hirsch-Fye algorithm.
We study an ionic Hubbard model in one dimension by using a quantum Monte Carlo method. The ionic Hubbard model we study is a strongly correlated model for electrons in alternating potentials with a lattice period of 2. By modifying the Hirsch-Fye algorithm, we incorporate the local variation of the potential in the weight of the partition function at each imaginary-time segment. For a given alternating potential of strength Δ = 1.0, as we change the local Coulomb repulsion, U , between spin-up and spin-down electrons, we calculate various physical quantities, including kinetic energies, correlation energies, and band energies. We try to find the signature of the phase transitions from the behavior of the kinetic energy and that of double occupancies at finite temperatures.
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.
Quantum Monte Carlo methods are used to study a quantum phase transition in a one-dimensional Hubbard model with on-site interaction U and a staggered ionic potential $(\ensuremath{\Delta}).$ Using recently formulated methods, the electronic polarization and localization are determined directly from the correlated ground state wave function and compared to results of previous work using exact diagonalization and Hartree-Fock. We find a transition from a symmetric band insulator to a broken-symmetry bond ordered (BO) phase as the ratio of $U/\ensuremath{\Delta}$ is increased, with a metallic point at the transition. Since it is known that at $\ensuremath{\Delta}=0$ the usual Hubbard model is a Mott insulator (MI) with no long-range order, we have searched for a second transition to this state by (i) increasing U at fixed $\ensuremath{\Delta}$ and (ii) decreasing $\ensuremath{\Delta}$ at fixed U. We find no transition from the BO to MI state, and we propose that the symmetric MI state exists only at $\ensuremath{\Delta}\ensuremath{\equiv}0.$
We investigate the ground-state phase diagram of the one-dimensional 'ionic' Hubbard model with an alternating periodic potential at half-filling by numerical diagonalization of finite systems with the Lanczos and density matrix renormalization group methods. We identify an insulator–insulator phase transition from a band to a correlated insulator with simultaneous charge and bond-charge order. The transition point is characterized by the vanishing of the optical excitation gap while simultaneously the charge and spin gaps remain finite and equal. Indications for a possible second transition into a Mott-insulator phase are discussed.
We investigate the ground-state phase diagram of the one-dimensional half-filled Hubbard model with an alternating potential---a model for the charge-transfer organic materials and the ferroelectric perovskites. We numerically determine the global phase diagram of this model using the level-crossing and the phenomenological renormalization-group methods based on the exact diagonalization calculations. Our results support the mechanism of the double phase transitions between Mott and band insulators pointed out by Fabrizio, Gogolin, and Nersesyan [Phys. Rev. Lett. 83, 2014 (1999)]: We confirm the existence of the spontaneously dimerized phase as an intermediate state. Further we provide numerical evidence to check the criticalities on the phase boundaries. Especially, we perform the finite-size-scaling analysis of the excitation gap to show the two-dimensional Ising transition in the charge part. On the other hand, we confirm that the dimerized phase survives in the strong-coupling limit, which is one of the resultants of competition between the ionicity and correlation effects.
We study Mott transition in the two-dimensional Hubbard model on an anisotropic triangular lattice. We use the Lanczos exact diagonalization of finite-size clusters up to eighteen sites, and calculate Drude weight, charge gap, double occupancy and spin structure factor. We average these physical quantities over twisted boundary conditions in order to reduce finite-size effects. We find a signature of the Mott transition in the dependence of the Drude weight and/or charge gap on the system size. We also examine the possibility of antiferromagnetic order from the spin structure factor. Combining these information, we propose a ground-state phase diagram which has a nonmagnetic insulating phase between a metallic phase and an insulating phase with antiferromagnetic order. Finally, we compare our results with those reported in the previous theoretical studies, and discuss the possibility of an unconventional insulating state.