One-dimensional Hubbard Model Research Articles

Sampling Error Analysis in Quantum Krylov Subspace Diagonalization

Quantum Krylov subspace diagonalization (QKSD) is an emerging method used in place of quantum phase estimation in the early fault-tolerant era, where limited quantum circuit depth is available. In contrast to the classical Krylov subspace diagonalization (KSD) or the Lanczos method, QKSD exploits the quantum computer to efficiently estimate the eigenvalues of large-size Hamiltonians through a faster Krylov projection. However, unlike classical KSD, which is solely concerned with machine precision, QKSD is inherently accompanied by errors originating from a finite number of samples. Moreover, due to difficulty establishing an artificial orthogonal basis, ill-conditioning problems are often encountered, rendering the solution vulnerable to noise. In this work, we present a nonasymptotic theoretical framework to assess the relationship between sampling noise and its effects on eigenvalues. We also propose an optimal solution to cope with large condition numbers by eliminating the ill-conditioned bases. Numerical simulations of the one-dimensional Hubbard model demonstrate that the error bound of finite samplings accurately predicts the experimental errors in well-conditioned regions.

Fluctuation-dissipation theorem and entanglement dynamics in one-dimensional low-density Jaynes-Cummings Hubbard model after a quench.

In one-dimensional low-density Jaynes-Cummings Hubbard (JCH) models [Phys. Rev. E 106, 064107 (2022)2470-004510.1103/PhysRevE.106.064107], we proved that the eigenstate thermalization hypothesis (ETH) is valid when the tunneling strength and coupling strength are of the same order. Surprisingly, at the weak tunneling limit, we observed that the entanglement entropy and scaling law of kinetic energy operators also exhibit obvious quantum chaotic properties, this is an unexpected result. To substantiate these findings, we further discuss their nonequilibrium dynamics in this paper. Our analysis reveals that when the model is a weak tunneling limit after the quench and the initial state is an equilibrium state of chaos, the system reaches an equilibrium state. This observation supports the conclusion that the low-density JCH model at the weak tunneling limit is nonintegrable, corroborating our previous results [Phys. Rev. E 106, 064107 (2022)2470-004510.1103/PhysRevE.106.064107]. Additionally, by discussing the validity of the fluctuation-dissipation theorem (FDT) and the evolution behavior of entanglement entropy and fidelity, we numerically demonstrate the differences between the one-dimensional low-density JCH model and general nonintegrable systems. Specifically, in the low-density JCH model, when the Hamiltonian after the quench is integrable, the validity of FDT depends on the thermal behavior of the initial Hamiltonian, and a metastable state is observed during the evolution of entanglement entropy. Our research presents an an intriguing and unique nonintegrable model, enriching the current understanding of nonintegrable systems.

Unusual low-temperature behavior in the half-filled band of the one-dimensional extended Hubbard model in atomic limit.

Recently, a kind of finite-temperature pseudotransition was observed in several quasi-one-dimensional models. In this work, we consider a genuine one-dimensional extended Hubbard model in the atomic limit, influenced by an external magnetic field and with the arbitrary number of particles controlled by the chemical potential. The one-dimensional extended Hubbard model in the atomic limit was initially studied in the seventies and has been investigated over the past decades, but it still surprises us today with its fascinating properties. We rigorously analyze its low-temperature behavior using the transfer matrix technique and provide accurate numerical results. Our analysis confirms that there is an anomalous behavior in the half-filled band, specifically occurring between the alternating pair (AP) and paramagnetic (PM) phases at zero temperature. Previous investigations did not deeply identify this anomalous behavior, maybe due to the numerical simplicity of the model, but from an analytical point of view this is not so easy to manipulate algebraically because one needs to solve an algebraic cubic equation. In this study, we explore this behavior and clearly distinguish the pseudotransition, which could easily be mistaken with a real phase transition. This anomalous behavior mimics features of both first- and second-order phase transitions. However, due to its nature, we cannot expect a finite-temperature phase transition in this model.

Non-Hermitian linear-response theory for spin diffusion in quantum systems

Spin transport theory of non-Hermitian quantum systems, where the non-Hermiticity led to a series of exotic phenomena when the system experiences dissipation to an environment is proposed. The goal is a better understanding of the spin diffusion in the one-dimensional non-Hermitian XXZ model which allows the transport coefficients be computed analytically. Moreover, we analyzed the electric transport in the one-dimensional non-Hermitian Hubbard model which is a very important model of electrons strongly correlated, where we have investigated the effect of non-Hermitian parameters like the imaginary hopping on AC and DC conductivities of the system. We analyzed the large U limit of this model, where such non-Hermiticity contributes with a minus sign in the virtual exchange of quasiparticles and the behavior of the ground state energy and low-lying excitations is reversed.

Quantum phase transitions in one-dimensional nanostructures: a comparison between DFT and DMRG methodologies.

In the realm of quantum chemistry, the accurate prediction of electronic structure and properties of nanostructures remains a formidable challenge. Density functional theory (DFT) and density matrix renormalization group (DMRG) have emerged as two powerful computational methods for addressing electronic correlation effects in diverse molecular systems. We compare ground-state energies ( ), density profiles ( ), and average entanglement entropies ( ) in metals, insulators and at the transition from metal to insulator, in homogeneous, superlattices, and harmonically confined chains described by the fermionic one-dimensional Hubbard model. While for the homogeneous systems, there is a clear hierarchy between the deviations, , and all the deviations decrease with the chain size; for superlattices and harmonic confinement, the relation among the deviations is less trivial and strongly dependent on the superlattice structure and the confinement strength considered. For the superlattices, in general, increasing the number of impurities in the unit cell represents lower precision in the DFT calculations. For the confined chains, DFT performs better for metallic phases, while the highest deviations appear for the Mott and band-insulator phases. This work provides a comprehensive comparative analysis of these methodologies, shedding light on their respective strengths, limitations, and applications. The DFT calculations were performed using the standard Kohn-Sham scheme within the BALDA approach. It integrated the numerical Bethe-Ansatz (BA) solution of the Hubbard model as the homogeneous density functional within a local-density approximation (LDA) for the exchange-correlation energy. The DMRG algorithms were implemented using the ITensor library, which is based on the matrix product states (MPS) ansatz. The calculations were performed until the energy reaches convergence of at least .

Level statistics of the one-dimensional dimerized Hubbard model

The statistical properties of level spacings provide valuable insights into the dynamical properties of a many-body quantum systems. We investigate the level statistics of the Fermi–Hubbard model with dimerized hopping amplitude and find that after taking into account translation, reflection, spin and η pairing symmetries to isolate irreducible blocks of the Hamiltonian, the level spacings in the limit of large system sizes follow the distribution expected for hermitian random matrices from the Gaussian orthogonal ensemble. We show this by analyzing the distribution of the ratios of consecutive level spacings in this system, its cumulative distribution and quantify the deviations of the distributions using their mean, standard deviation and skewness.

Extracting off-diagonal order from diagonal basis measurements

Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multiscale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-dimensional extended Hubbard model with an off-diagonal order parameter, sandwiched between diagonal charge and spin density wave phases, using only diagonal descriptors. We study the model directly in the thermodynamic limit using the recently developed variational uniform matrix product states algorithm, and draw our samples from degenerate ground states related by global spin rotations, emulating the projective measurements that are accessible in experiments. Our results will have important implications for the study of exotic phases using optical lattice experiments. Published by the American Physical Society 2024

Proof of Completeness of the Local Conserved Quantities in the One-Dimensional Hubbard Model

Proof of Completeness of the Local Conserved Quantities in the One-Dimensional Hubbard Model

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.