One-dimensional Hubbard Model Research Articles

Flat Band and η-Pairing States in a One-Dimensional Moiré Hubbard Model

A Moiré system is formed when two periodic structures have a slightly mismatched period, resulting in unusual strongly correlated states in the presence of particle-particle interactions. The periodic structures can arise from the intrinsic crystalline order and periodic external field. We investigate a one-dimensional Hubbard model with periodic on-site potential of period n 0, which is commensurate to the lattice constant. For large n 0, the exact solution demonstrates that there is a midgap flat band with zero energy in the absence of Hubbard interaction. Each Moiré unit cell contributes two zero energy levels to the flat band. In the presence of Hubbard interaction, the midgap physics is demonstrated to be well described by a uniform Hubbard chain in which the effective hopping and on-site interaction strength can be controlled by the amplitude and period of the external field. Numerical simulations are performed to demonstrate the correlated behaviors in the finite-sized Moiré Hubbard system, including the existence of an η-pairing state and bound pair oscillation. This finding provides a method to enhance the correlated effect by a spatially periodic external field.

Spectral properties of a one-dimensional extended Hubbard model from bosonization and time-dependent variational principle: Applications to one-dimensional cuprates

Spectral properties of a one-dimensional extended Hubbard model from bosonization and time-dependent variational principle: Applications to one-dimensional cuprates

All Local Conserved Quantities of the One-Dimensional Hubbard Model.

We present the exact expression for all local conserved quantities of the one-dimensional Hubbard model. We identify the operator basis constructing the local charges and find that nontrivial coefficients appear in the higher-order charges. We derive the recursion equation for these coefficients, and some of them are explicitly given. There are no other local charges independent of those we obtained.

Surface superconductor-insulator transition: Reduction of the critical electric field by Hartree-Fock potential

Recently, a surface superconductor-insulator transition has been predicted for a bulk superconductor in an electric field applied perpendicular to its surface. The related calculations were performed within a one-dimensional Hubbard model by numerically solving the Bogoliubov-de Gennes (BdG) equations without the Hartree-Fock (HF) interaction potential. The phase diagram of the surface superconducting, metallic, and insulating states was obtained as dependent on the electric field and temperature. This diagram was found to be in agreement with experimental results reported previously for (Li,Fe)OHFeSe thin flakes. In the present work, by taking into account the HF potential, we find that the latter acts as a kind of an extra electrostatic potential that enhances the electric-field effects on the surface states. The qualitative features of the phase diagram remain the same but the surface superconductor-insulator transition occurs at significantly lower electric fields, which supports prospects of its experimental observation in bulk samples.

Pump-probe spectroscopy of the one-dimensional extended Hubbard model at half filling

Pump-probe spectroscopy of the one-dimensional extended Hubbard model at half filling

On integrability of the one-dimensional Hubbard model

We find a family of solutions to Zamolodchikov's tetrahedral algebra corresponding to the fermionic R-operator for the free fermion model of the difference type in one of the spectral parameters, construct an extension of the R-operator for a system of two spins satisfying the Yang-Baxter equation, and find the local charges. We also construct a twisted monodromy operator, which leads to the one-dimensional Hubbard model.

Waterlike density anomaly in fermions

In this work we explore the one-dimensional extended Hubbard model as a fluid system modeling liquid phases of different densities. This model naturally displays two length scales of interaction, which are connected with waterlike anomalies. We analyze the density anomaly as a function of the model parameters, namely the hopping, on-site and first neighbor interactions. We show that this anomaly is present for a wide range of model parameters and is connected to a ground-state liquid–liquid critical point.

Umklapp scattering in the one-dimensional Hubbard model

Umklapp scattering in the one-dimensional Hubbard model

Classical Lie bialgebras for AdS/CFT integrability by contraction and reduction

Integrability of the one-dimensional Hubbard model and of the factorised scattering problem encountered on the worldsheet of AdS strings can be expressed in terms of a peculiar quantum algebra. In this article, we derive the classical limit of these algebraic integrable structures based on established results for the exceptional simple Lie superalgebra \mathfrak{d}(2,1;\epsilon)𝔡(2,1;ϵ) along with standard \mathfrak{sl}(2)𝔰𝔩(2) which form supersymmetric isometries on 3D AdS space. The two major steps in this construction consist in the contraction to a 3D Poincaré superalgebra and a certain reduction to a deformation of the \mathfrak{u}(2|2)𝔲(2|2) superalgebra. We apply these steps to the integrable structure and obtain the desired Lie bialgebras with suitable classical r-matrices of rational and trigonometric kind. We illustrate our findings in terms of representations for on-shell fields on AdS and flat space.

Disorder-induced spin-charge separation in the one-dimensional Hubbard model

Many-body localization is believed to be generically unstable in quantum systems with continuous non-Abelian symmetries, even in the presence of strong disorder. Breaking these symmetries can stabilize the localized phase, leading to the emergence of an extensive number of quasilocally conserved quantities known as local integrals of motion, or $l$ bits. Using a sophisticated nonperturbative technique based on continuous unitary transforms, we investigate the one-dimensional Hubbard model subject to both spin and charge disorder, compute the associated $l$ bits and demonstrate that the disorder gives rise to a novel form of spin-charge separation. We examine the role of symmetries in delocalizing the spin and charge degrees of freedom, and show that while symmetries generally lead to delocalization through multiparticle resonant processes, certain subsets of states appear stable.

Quench dynamics in the one-dimensional mass-imbalanced ionic Hubbard model

Using the time-dependent Lanczos method, we study the nonequilibrium dynamics of the one-dimensional ionic-mass imbalanced Hubbard chain driven by a quantum quench of the on-site Coulomb interaction, where the system is prepared in the ground state of the Hamiltonian with a different Hubbard interaction. The exact diagonalization method (Lanczos algorithm) is adopted to study the zero temperature phase diagram in equilibrium, which is shown to be in good agreement with previous studies using density matrix renormalization group (DMRG). We then study the nonequilibrium quench dynamics of the spin and charge order parameters by fixing the initial and final Coulomb interactions while changing the quenching time protocols. Our study shows that the time evolution of the charge and spin order parameters strongly depend on the quenching time protocols. In particular, the effective temperature of the system will decrease monotonically as the quenching time is increased. By taking the final Coulomb interaction strength to be in the strong coupling regime, we find that the oscillation frequency of the charge order parameter increases monotonically with the Coulomb interaction. By contrast, the frequency of the spin order parameter decreases monotonically with increasing Coulomb interaction. We explain this result using an effective spin model and a two-site Hubbard model in the strong coupling limit. Finally, we take the final Coulomb interaction strength to be in the weak coupling regime and find that the oscillation frequency of both the charge and spin order parameters increases monotonically with decreasing Coulomb interaction. Our study suggests strategies to engineer the relaxation behavior of interacting quantum many-particle systems.

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.