One-dimensional Hubbard Model Research Articles

Disordered impenetrable two-component fermions in one dimension

We study the one-dimensional Hubbard model for two-component fermions with infinitely strong on-site repulsion ($t\ensuremath{-}0$ model) in the presence of disorder. Our analytical treatment demonstrates that the type of disorder drastically changes the nature of the emerging phases. The case of spin-independent disorder can be treated as a single-particle problem with Anderson localization. On the contrary, recent numerical findings show that spin-dependent disorder, which can be realized as a random magnetic field, leads to the many-body localization-delocalization transition. We find an explicit analytic expression for the matrix elements of the random magnetic field between the eigenstates of the $t\ensuremath{-}0$ model with potential disorder on a finite lattice. Analysis of the matrix elements supports the existence of the many-body localization-delocalization transition in this system and provides an extended physical picture of the random magnetic field.

Spin-incoherent liquid and interaction-driven criticality in the one-dimensional Hubbard model

Although the one dimensional (1D) repulsive Fermi-Hubbard model has been intensively studied over many decades, a rigorous understanding of many aspects of the model is still lacking. In this work, based on the solutions to the thermodynamic Bethe ansatz equations, we provide a rigorous study on the following: (1) We calculate the fractional excitations of the system in various phases, from which we identify the parameter regime featuring the spin incoherent Luttinger liquid (SILL). We investigate the universal properties and the asymprotic of correlation functions of the SILL. (2) We study the interaction-driven phase transition and the associated criticality, and build up an essential connection between the Contact susceptibilities and the variations of density, magnetization and entropy with respect to the interaction strength. As an application of these concepts, which hold true for higher dimensional systems, we propose a quantum cooling scheme based on the interaction-driven refrigeration cycle.

Controlling inversion and time-reversal symmetries by subcycle pulses in the one-dimensional extended Hubbard model

Owning to their high controllability, laser pulses have contributed greatly to our understanding of strongly correlated electron systems. However, typical multicycle pulses do not control the symmetry of systems, which plays an important role in the emergence of novel quantum phases. Here, we demonstrate that subcycle pulses whose oscillation is less than one period within a pulse envelope can control inversion and time-reversal symmetries in the electronic states of the one-dimensional extended Hubbard model. Using an ultrashort subcycle pulse, one can generate a steady electric current (SEC) in a photoexcited state due to an Aharonov-Bohm flux instantaneously introduced through the phase of an electric field. Consequently, time-reversal symmetry is broken. In contrast, a broad subcycle pulse does not induce SEC but instead generates electric polarization, thus breaking inversion symmetry. Both symmetry breakings in a photoexcited state can be monitored by second-harmonic generation. These findings provide a new methodology for designing the symmetries of electronic states and open up a new field of subcycle-pulse engineering.

The one-dimensional Hubbard model in a magnetic field: Density profiles and ground-state phase diagram

We study the single-band one-dimensional Hubbard model in an arbitrary magnetic field using the cumulant Green's functions method (CGFM) with clusters containing 5, 6, 7, and 8 sites. We calculate the occupation numbers as functions of the chemical potential in a positive magnetic field and identify a partially filled region with negative magnetization. We obtain the ground-state phase diagram in chemical potential vs. magnetic field coordinates. We include a straightforward application of the CGFM by using a cluster as a correlated quantum dot connected to correlated leads to construct a single-electron transistor that presents a spin-polarised conductance.

Thermofield Theory for Finite-Temperature Electronic Structure.

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. They also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

Drinfeld realization of the centrally extended psl(2|2) Yangian algebra with the manifest coproducts

The Lie superalgebra psl(2|2) is recognized as a quite special algebra. In mathematics, it has the vanishing Killing form and allows for the three-dimensional central extension. In physics, it shows up as the symmetry of the one-dimensional Hubbard model and the asymptotic S-matrix of the anti-de Sitter/conformal field theory (AdS/CFT) correspondence. In this paper, we present the Drinfeld realization of the Yangian algebra associated with the centrally extended Lie superalgebra psl(2|2). Furthermore, we show that it possesses Hopf algebra structures, particularly the coproducts. The idea to prove the existence of manifest coproducts is as follows. First, we shall introduce them to Levendorskii’s realization, a system of a finite truncation of Drinfeld generators. Second, we show that Levendorskii’s realization is isomorphic to the Drinfeld realization by induction.

Time-resolved single-particle spectrum of the one-dimensional extended Hubbard model after interaction quenches

We investigate the non-equilibrium dynamics of the one-dimensional extended Hubbard model after interaction quenches. In strong-coupling regime with large on-site interaction, the ground states of this model with small and large nearest-neighbor interactions are in spin-density-wave and charge-density-wave phases, respectively. Combining twisted boundary conditions with the time-dependent Lanczos method, we obtain snapshots of the time-dependent single-particle spectrum after quenches. We find that for quench within the same phase, the single-particle spectrum becomes close to that of the quenched Hamiltonian immediately after the quench. While for quench across the critical point, the afterward evolution process depends mainly on the distribution of the initial state among the eigenstates of the quenched Hamiltonian. Our finding may serve as a way to detect the phase transition in ultracold atom systems with interactions.

Relation of entanglement entropy and particle number fluctuations in one-dimensional Hubbard model

Relation of entanglement entropy and particle number fluctuations in one-dimensional Hubbard model

Doping effects in high-harmonic generation from correlated systems

Using the one-dimensional Hubbard model, which is commonly used for describing, e.g., high-${T}_{c}$ superconducting cuprates, we study high-harmonic generation (HHG) from doped, correlated materials. Doping is modeled by changing the number of electrons in the lattice from the conventional half-filling case. For relatively small Hubbard $U$, i.e., small electron-electron correlation, we find little to no effect of doping on the dynamics and the HHG spectra. For increasing $U$ the degree of doping has a marked effect on the dynamics and spectra. We explain these findings through the quasiparticle-based doublon-holon picture. The dynamics are separated into two types, first, doublon and holon movement, and, second, doublon-holon pair creation and annihilation. Doping results in all configurations containing doublons or holons. Those quasiparticles can move at no extra cost in energy regardless of the correlation level. This motion at no energy cost increases the high-harmonic gain for low- and medium-harmonic orders. We discuss that in the high-$U$ limit, antiferromagnetic ordering becomes increasingly unlikely with increasing doping rates and explain an associated drop in the high-order harmonics relative to the case of half-filling.

Quantum phase transitions in the matrix product states of the one-dimensional boson Hubbard model

We study quantum phase transitions in the matrix product states of the one-dimensional boson Hubbard model, whose amount of entanglement is limited by the size of the matrices used in the representation of the states. By measuring entanglement entropy and other physical properties, we observe that the Mott-insulator-to-superfluid transitions occur sharp and continuous, accompanied by shifting transition points that are not blurred by finite entanglement effects. This strongly suggests that the transition always occurs between the Mott insulator and a mean-field-like compressible state, followed by the more entangled superfluid state. Both $O(2)$ and the commensurate-incommensurate transitions are studied, whose properties can be characterized by entanglement spectra and critical exponents.

Normal and abnormal thermalization indicators in a one-dimensional low-density Jaynes-Cummings Hubbard model with and without dipole-dipole interaction.

In the one-dimensional low-density Jaynes-Cummings Hubbard (JCH) model, we find that when the hopping strength is much smaller than the coupling strength, the average restricted energy gap ratio exhibits an abnormal statistical behavior that is neither a Poisson nor a Gaussian orthogonal ensemble. But the average half-chain entanglement entropy exhibits ergodicity, and the eigenstate thermalization hypothesis (ETH) is valid for the observable. These results are quite different from those of the standard JCH model. In addition, when the hopping and the coupling strengths are of the same order, quantum chaos still appears in the low-density JCH model, which is in contrast to the integrability of the one-dimensional hard-core bosons. Finally, the dipole-dipole interaction breaks the particle-hole symmetry and leads the abnormal statistical properties to be closer to those of the integrable system at the weak hopping strength limit, but the quantum chaos properties cannot be affected when the hopping strength is of the same order as the coupling strength. Our results demonstrate the counterintuitive behavior in the low-density JCH model and explain the physics behind them from the perspective of the energy spectrum.

Resurgence and renormalons in the one-dimensional Hubbard model

We use resurgent analysis to study non-perturbative aspects of the one-dimensional, multicomponent Hubbard model with an attractive interaction and arbitrary filling. In the two-component case, we show that the leading Borel singularity of the perturbative series for the ground-state energy is determined by the energy gap, as expected for superconducting systems. This singularity turns out to be of the renormalon type, and we identify a class of diagrams leading to the correct factorial growth. As a consequence of our analysis, we propose an explicit expression for the energy gap at weak coupling in the multi-component Hubbard model, at next-to-leading order in the coupling constant. In the two-component, half-filled case, we use the Bethe ansatz solution to determine the full trans-series for the ground state energy, and the exact form of its Stokes discontinuity.

Variational counterdiabatic driving of the Hubbard model for ground-state preparation

Counterdiabatic (CD) protocols enable fast driving of quantum states by invoking an auxiliary adiabatic gauge potential (AGP) that suppresses transitions to excited states throughout the driving process. Usually, the full spectrum of the original unassisted Hamiltonian is a prerequisite for constructing the exact AGP, which implies that CD protocols are extremely difficult for many-body systems. Here, we apply a variational CD protocol recently proposed by P. W. Claeys et al. [Phys. Rev. Lett. 123, 090602 (2019)] to a two-component fermionic Hubbard model in one spatial dimension. This protocol engages an approximated AGP expressed as a series of nested commutators. We show that the optimal variational parameters in the approximated AGP satisfy a set of linear equations whose coefficients are given by the squared Frobenius norms of these commutators. We devise an exact algorithm that escapes the formidable iterative matrix-vector multiplications and evaluates the nested commutators and the CD Hamiltonian in analytic representations. We then examine the CD driving of the one-dimensional Hubbard model up to $L = 14$ sites with driving order $l \leqslant 3$. Our results demonstrate the usefulness of the variational CD protocol to the Hubbard model and permit a possible route towards fast ground-state preparation for many-body systems.

Excitonic Correlations, Spin-State Ordering, and Magnetic-Field Effects in One-Dimensional Two-Orbital Hubbard Model for Spin-Crossover Region

The electronic properties of excitonic insulators have been examined precisely in recent years. Pictures of exciton condensation may be applied to the spin-state transition observed in perovskite cobalt oxides. We examine the crystal-field and magnetic-field dependences of spatial spin structures on the basis of the density matrix renormalization group method using an effective model for the one-dimensional two-orbital Hubbard model in strong-coupling limit. We find an excitonic insulating (EI) phase and a spin-state ordering (SSO) phase in the intermediate region between low-spin and high-spin phases. In the EI phase, spin-triplet excitons are spatially fluctuating due to quantum effects, and an incommensurate spin correlation realizes. The analyses of a spin gap and degeneracy of entanglement spectra suggest the realization of the Haldane-like edge state in the EI phase. In the SSO phase, 3-fold or incommensurate SSO structures realize depending on the crystal-field splitting. These structures are stabilized as a result of the competition of exchange interactions between spin states.

Integrable quenches in the Hubbard model

We study the quench dynamics of the one-dimensional Hubbard model through the quench action formalism. We introduce a class of integrable initial states—expressed as product states over two sites—for which we can provide an exact characterisation of the late-time regime. This is achieved by finding a closed-form expression for the overlaps between our states and the Bethe ansatz eigenstates, which we check explicitly in the limits of low densities and infinite repulsion. Our solution gives access to the stationary values attained by local observables (we show the explicit example of the density of doubly occupied sites) and the asymptotic entanglement dynamics directly in the thermodynamic limit. Interestingly, we find that for intermediate interaction strength Rényi entropies display a double-slope structure.

Level statistics of the one-dimensional ionic Hubbard model

In this paper we analyze the spectral level statistics of the one-dimensional ionic Hubbard model, the Hubbard model with an alternating on-site potential. In particular, we focus on the statistics of the gap ratios between consecutive energy levels. This quantity is often used in order to signal whether a many-body system is integrable or chaotic. A chaotic system has typically the statistics of a Gaussian ensemble of random matrices while the spectral properties of the integrable system follow a Poisson statistics. We find that whereas the Hubbard model without alternating potential is known to be integrable and its spectral properties follow a Poissonian statistics, the presence of an alternating potential causes a drastic change in the spectral properties, which resemble the one of a Gaussian ensemble of random matrices. However, to uncover this behavior one has to separately consider the blocks of all symmetries of the ionic Hubbard model.

Tunable transport in the mass-imbalanced Fermi-Hubbard model

The late-time dynamics of quantum many-body systems is organized in distinct dynamical universality classes, characterized by their conservation laws and thus by their emergent hydrodynamic transport. Here, we study transport in the one-dimensional Hubbard model with different masses of the two fermionic species. To this end, we develop a quantum Boltzmann approach valid in the limit of weak interactions. We explore the crossover from ballistic to diffusive transport, whose timescale strongly depends on the mass ratio of the two species. For timescales accessible with matrix product operators, we find excellent agreement between these numerically exact results and the quantum Boltzmann equation, even for intermediate interactions. We investigate two scenarios which have been recently studied with ultracold atom experiments. First, in the presence of a tilt, the quantum Boltzmann equation predicts that transport is significantly slowed down and becomes subdiffusive, consistent with previous studies. Second, we study transport probed by displacing a harmonic confinement potential and find good quantitative agreement with recent experimental data [N. Darkwah Oppong et al., arXiv:2011.12411]. Our results demonstrate that the quantum Boltzmann equation is a useful tool to study complex non-equilibrium states in inhomogeneous potentials, as often probed with synthetic quantum systems.

Glassy dynamics of the one-dimensional Mott insulator excited by a strong terahertz pulse

The elucidation of nonequilibrium states in strongly correlated systems holds the key to emergence of novel quantum phases. The nonequilibrium-induced insulator-to-metal transition is particularly interesting since it reflects the fundamental nature of competition between itinerancy and localization of the charge degrees of freedom. We investigate pulse-excited insulator-to-metal transition of the half-filled one-dimensional extended Hubbard model. Calculating the time-dependent optical conductivity with the time-dependent density-matrix renormalization group, we find that strong mono- and half-cycle pulses inducing quantum tunneling strongly suppress spectral weights contributing to the Drude weight σD, even if we introduce a large number of carriers Δnd. This is in contrast to a metallic behavior of σD∝Δnd induced by photon absorption and chemical doping. The strong suppression of σD in quantum tunneling is a result of the emergence of the Hilbert-space fragmentation, which makes pulse-excited states glassy.Received 9 February 2022Revised 15 May 2022Accepted 18 July 2022DOI:https://doi.org/10.1103/PhysRevResearch.4.L032019Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasLocalizationNonequilibrium statistical mechanicsOptical conductivityPhotoinduced effectQuantum statistical mechanicsPhysical Systems1-dimensional systemsGlassy systemsMott insulatorsCondensed Matter, Materials & Applied PhysicsStatistical PhysicsQuantum Information

Oscillating charge currents of one-dimensional Hubbard model in electric field

The time evolution properties of charge current for the one-dimensional Hubbard model in an electric field have been studied in a rigorous manner. We find that there is a complete and orthonormal set of time-evolution states for which the charge current can only keep zero or oscillate constantly, differing from the possible picture of damped or over-damped Bloch oscillations due to strong correlations. It is also found that, associated with these states, there is a set of constant phase factors, which are uniquely determined and are very useful on discussing the long-time evolution behaviors of the system.

Introducing screening in one-body density matrix functionals: Impact on charged excitations of model systems via the extended Koopmans' theorem

In this work we get insight into the impact of reduced density matrix functionals on the quality of removal/addition energies obtained using the extended Koopmans' theorem (EKT). Within the reduced density matrix functional theory (RDMFT) the EKT approach reduces to a matrix diagonalization, whose ingredients are the one- and two-body reduced density matrices. A striking feature of the EKT within RDMFT is that it opens a band gap, although it is too large, in strongly correlated materials, which are a challenge for state-of-the-art methods such as $GW$. Using the one-dimensional Hubbard model and the homogeneous electron gas as test cases, we find that (i) with exact or very accurate density matrices the EKT systematically overestimates the band gap in the Hubbard model and the bandwidth in the homogeneous electron gas and (ii) with approximate density matrices, instead, the EKT can benefit from error cancellation. In particular we test a new approximation which combines random-phase approximation screening with the power functional approximation to the two-body reduced density matrix introduced by Sharma et al. [Phys. Rev. B 78, 201103 (2008)]. An important feature of this approximation is that it reduces the EKT band gap in the studied models; it is hence a promising approximation for correcting the EKT band-gap overestimation in strongly correlated materials.