The origin of degeneracies and crossings in the 1d Hubbard model

The paper is devoted to the connection between integrability of a finite quantum system and degeneracies of its energy levels. In particular, we analyse in detail the energy spectra of finite Hubbard chains. Heilmann and Lieb demonstrated that in these systems there are crossings of levels of the same parameter-independent symmetry. We show that this apparent violation of the Wigner-von Neumann noncrossing rule follows directly from the existence of nontrivial conservation laws and is a characteristic signature of quantum integrability. The energy spectra of Hubbard chains display many instances of permanent (at all values of the coupling) twofold degeneracies that cannot be explained by parameter-independent symmetries. We relate these degeneracies to the different transformation properties of the conserved currents under spatial reflections and the particle-hole transformation and estimate the fraction of doubly degenerate states. We also discuss multiply degenerate eigenstates of the Hubbard Hamiltonian. The wavefunctions of many of these states do not depend on the coupling, which suggests the existence of an additional parameter-independent symmetry.

The inhomogeneous extended Bose-Hubbard model: A mean-field theory

We develop an inhomogeneous mean-field theory for the extended Bose-Hubbard model with a quadratic, confining potential. In the absence of this potential, our mean-field theory yields the phase diagram of the homogeneous extended Bose-Hubbard model. This phase diagram shows a superfluid (SF) phase and lobes of Mott-insulator (MI), density-wave (DW), and supersolid (SS) phases in the plane of the chemical potential mu and on-site repulsion U; we present phase diagrams for representative values of V, the repulsive energy for bosons on nearest-neighbor sites. We demonstrate that, when the confining potential is present, superfluid and density-wave order parameters are nonuniform; in particular, we obtain, for a few representative values of parameters, spherical shells of SF, MI, DW, and SS phases. We explore the implications of our study for experiments on cold-atom dipolar condensates in optical lattices in a confining potential.

Bose-Hubbard model in a strong effective magnetic field: Emergence of a chiral Mott insulator ground state

Motivated by experiments on Josephson junction arrays, and cold atoms in an optical lattice in a synthetic magnetic field, we study the ``fully frustrated'' Bose-Hubbard model with half a magnetic flux quantum per plaquette. We obtain the phase diagram of this model on a two-leg ladder at integer filling via the density matrix renormalization group approach, complemented by Monte Carlo simulations on an effective classical XY model. The ground state at intermediate correlations is consistently shown to be a chiral Mott insulator (CMI) with a gap to all excitations and staggered loop currents which spontaneously break time-reversal symmetry. We characterize the CMI state as a vortex supersolid or an indirect exciton condensate, and discuss various experimental implications.

Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium

We derive exact expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electromagnetic field. We also evaluate these expressions for the homogeneous case in equilibrium, where all time dependence and external fields vanish. Unlike similar sum rules for the Fermi-Hubbard model, in the Bose-Hubbard model case, the sum rules often depend on expectation values that cannot be determined simply from parameters in the Hamiltonian like the interaction strength and chemical potential but require knowledge of equal-time many-body expectation values from some other source. We show how one can approximately evaluate these expectation values for the Mott-insulating phase in a systematic strong-coupling expansion in powers of the hopping divided by the interaction. We compare the exact moment relations to the calculated moments of spectral functions determined from a variety of different numerical approximations and use them to benchmark their accuracy. DOI: 10.1103/PhysRevA.87.013628

Chiral Mott insulator with staggered loop currents in the fully frustrated Bose-Hubbard model

Motivated by experiments on Josephson junction arrays in a magnetic field and ultracold interacting atoms in an optical lattice in the presence of a ``synthetic'' orbital magnetic field, we study the ``fully frustrated'' Bose-Hubbard model and quantum XY model with half a flux quantum per lattice plaquette. Using Monte Carlo simulations and the density matrix renormalization group method, we show that these kinetically frustrated boson models admit three phases at integer filling: a weakly interacting chiral superfluid phase with staggered loop currents which spontaneously break time-reversal symmetry, a conventional Mott insulator at strong coupling, and a remarkable ``chiral Mott insulator'' (CMI) with staggered loop currents sandwiched between them at intermediate correlation. We discuss how the CMI state may be viewed as an exciton condensate or a vortex supersolid, study a Jastrow variational wave function which captures its correlations, present results for the boson momentum distribution across the phase diagram, and consider various experimental implications of our phase diagram. Finally, we consider generalizations to a staggered flux Bose-Hubbard model and a two-dimensional (2D) version of the CMI in weakly coupled ladders.

Phase diagram of the half-filled ionic Hubbard model

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 Delta and the on-site Hubbard U. We find that for a finite Delta and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold U = U-AF via a first-order phase transition. For U smaller than U-AF 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,Delta) or, equivalently, on decreasing U below U-AF(T,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 U-AF. At U-N > U-AF, when the Neel temperature T-N for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second-order transition. For U >> Delta, T-N similar to t(2)/U(1 - x(2)), where x = 2 Delta/U and thus T-N increases with increase in Delta/U. In the three-dimensional parameter space of (U/t, T/t, and Delta/t), as T increases, the surface of first-order transition at U-AF(T,Delta) and that of the second-order transition at U-N(T,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.

Feshbach modulation spectroscopy of the Fermi-Hubbard model

In the vicinity of a Feshbach resonance, a system of ultracold atoms in an optical lattice undergoes rich physical transformations which involve molecule formation and hopping of molecules on the lattice and thus goes beyond a single-band Hubbard model description. We explore theoretically the response of this system to a harmonic modulation of the magnetic field, and thus of the scattering length, across the Feshbach resonance. In the regime in which the single-band Hubbard model is still valid, we provide results for the doublon production as a function of the various parameters, such as frequency, amplitude, etc., that characterize the field modulation, as well as the lattice depth. The method may uncover a route towards the efficient creation of ultracold molecules and also provide an alternative to conventional lattice-depth-modulation spectroscopy.

Classical two-site fidelity and quantum phase transitions in the Ising model and the extended Hubbard model

In this Letter, the classical two-site-ground-state fidelity (CTGF) is exploited to identify quantum phase transitions (QPTs) for the transverse field Ising model (TFIM) and the one-dimensional extended Hubbard model (EHM). Our results show that the CTGF exhibits an abrupt change around the regions of criticality and can be used to identify QPTs in spin and fermionic systems. The method is especially convenient when it is connected with the density-matrix renormalization group (DMRG) algorithm. (C) 2008 Elsevier B.V. All rights reserved.

Magnetic orderings and phase separations in the zero-bandwidth limit of the extended Hubbard model with intersite magnetic interactions

This is an author-created, un-copyedited version of an article accepted for publication in Acta Physica Polonica A. The Version of Record is available online at http://przyrbwn.icm.edu.pl/APP/PDF/118/a118z2p31.pdf

Possible Exotic Phases in the One-Dimensional Extended Hubbard Model

We investigate numerically the ground state phase diagram of the one-dimensional extended Hubbard model, including an on--site interaction U and a nearest--neighbor interaction V. We focus on the ground state phases of the model in the V >> U region, where previous studies have suggested the possibility of dominant superconducting pairing fluctuations before the system phase separates at a critical value V=V_PS. Using quantum Monte Carlo methods on lattices much larger than in previous Lanczos diagonalization studies, we determine the boundary of phase separation, the Luttinger Liquid correlation exponent K_rho, and other correlation functions in this region. We find that phase separation occurs for V significantly smaller than previously reported. In addition, for negative U, we find that a uniform state re-enters from phase separation as the electron density is increased towards half filling. For V < V_PS, our results show that superconducting fluctuations are not dominant. The system behaves asymptotically as a Luttinger Liquid with K_rho < 1, but we also find strong low-energy (but gapped) charge-density fluctuations at a momentum not expected for a standard Luttinger Liquid.

Kondo-attractive-Hubbard model for the ordering of local magnetic moments in superconductors

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Spin, charge, and orbital correlations in the one-dimensional t(2g)-orbital Hubbard model

We present the zero-temperature phase diagram of the one-dimensional t(2g)-orbital Hubbard model, obtained using the density-matrix renormalization group and Lanczos techniques. Emphasis is given to the case of the electron density n=5 corresponding to five electrons per site, while several other cases for electron densities between n=3 and 6 are also studied. At n=5, our results indicate a first-order transition between a paramagnetic (PM) insulator phase, with power-law slowly decaying correlations, and a fully polarized ferromagnetic (FM) state by tuning the Hund's coupling. The results also suggest a transition from the n=5 PM insulator phase to a metallic regime by changing the electron density, either via hole or electron doping. The behavior of the spin, charge, and orbital correlation functions in the FM and PM states are also described in the text and discussed. The robustness of these two states against varying parameters suggests that they may be of relevance in quasi-one-dimensional Co-oxide materials, or even in higher dimensional cobaltite systems as well.

Competition between spin, charge, and bond waves in a Peierls-Hubbard model

We study a one-dimensional extended Peierls-Hubbard model coupled to intracell and intercell phonons for a half-filled band. The calculations are made using the Hartree-Fock and adiabatic approximations for arbitrary temperature. In addition to static spin, charge, and bond density waves, we predict intermediate phases that lack inversion symmetry, and phase transitions that reduce symmetry on increasing temperature.