We study extended Hubbard and Peierls-Hubbard models with up to N=10 sites and half-filled bands with modified boundary conditions. For the extended Hubbard model, the finite-size dependence of the transition from the charge-density-wave regime (CDW) to the spin-density-wave (SDW) regime and of the condensation transition is critically examined, clarifying some previous results. For the Peierls-Hubbard model and N\ensuremath{\le}10 with modified periodic boundary conditions, we always see a decrease of dimerization with U. Extrapolating these results to infinite N, we find that dimerization is enhanced by U for small electron-phonon coupling \ensuremath{\alpha}\ifmmode \tilde{}\else \~{}\fi{}0.75, whereas it decreases with U for 0.75\ensuremath{\le}1. By extending this model by a nearest-neighbor repulsion V, the amplitude of dimerization reaches a maximum at the same value U\ensuremath{\le}2V, where the CDW-to-SDW transition occurs in the extended Hubbard model.
We study extended Hubbard and Peierls-Hubbard models with up to N=10 sites and half-filled bands with modified boundary conditions. For the extended Hubbard model, the finite-size dependence of the transition from the charge-density-wave regime (CDW) to the spin-density-wave (SDW) regime and of the condensation transition is critically examined, clarifying some previous results. For the Peierls-Hubbard model and N\ensuremath{\le}10 with modified periodic boundary conditions, we always see a decrease of dimerization with U. Extrapolating these results to infinite N, we find that dimerization is enhanced by U for small electron-phonon coupling \ensuremath{\alpha}\ifmmode \tilde{}\else ~{}\fi{}0.75, whereas it decreases with U for 0.75\ensuremath{\le}1. By extending this model by a nearest-neighbor repulsion V, the amplitude of dimerization reaches a maximum at the same value U\ensuremath{\le}2V, where the CDW-to-SDW transition occurs in the extended Hubbard model.
"We have studied the extended Hubbard and the Peierls-extended Hubbard models in one dimension using primarily computational methods. We discuss the derivation of the models and the computational methods used to study them, focusing primarily on our extensions over existing numerical methods. In particular, we present the construction of ""operator-loop"" updates for spin-1/2 fermions and the inclusion of electron-phonon interactions in the Stochastic Series Expansion Quantum Monte Carlo method and the use of Quantum Parallel Tempering to improve the quality of the numerical data. For the extended Hubbard model at half-filling, we have confirmed the existence of a recently reported bond-order wave phase for the ground state in a small region of parameter space between the well-known charge- and spin-density wave phases. We have studied in detail the nature of the associated quantum phase transitions between all the phases. Upon doping away from half-filling, the doubly-degenerate, spin and charge gapped long-range-ordered ground states of the model are shown to evolve to a Luther-Emery liquid with no charge gap but a finite spin gap. For the Peierls-extended Hubbard model, which includes both electron-phonon and electron-electron interactions, we have studied the dimerization transition of the half-filled system in the presence of fully quantum mechanical phonons. We show that the transition to a dimerized state occurs only when the electron-phonon coupling exceeds a finite critical value that depends on the bare phonon frequency. This is in contrast to previous studies in similar models with adiabatic phonons that had predicted a dimerized state for arbitrarily small values of the electron-phonon coupling."
It has been well established that in both the weak- and strong-coupling limits, the phase transition from a spin-density wave (SDW) to a charge-density wave (CDW) in the half-filled extended-Hubbard model occurs precisely at the ratio U/Vc=2, where U and V are the on-site and intersite electron-electron interactions, respectively. However, a recent density-matrix renormalization group (DMRG) study and a previous Monte Carlo simulation jointly suggest a picture that challenges the well-established weak-coupling result. Here, a careful calculation shows such a picture is questionable. In particular, it clarifies the DMRG result, even with a moderate number of states kept in the density matrix, agrees with both the weak- and strong-coupling findings, and predicts a CDW/SDW phase boundary smoothly links those two coupling limits. Possible sources of these previous failures are discussed.
We study the effects of a bond-charge interaction (or correlated hopping) on the properties of the extended [i.e., with both on-site (U) and nearest-neighbor (V) repulsions] Hubbard model in one dimension at half-filling. Energy gaps and correlation functions are calculated by Lanczos diagonalization on finite systems. We find that, irrespective of the sign of the bond-charge interaction, X, the charge-density-wave (CDW) state is more robust than the spin-density-wave (SDW) state. A small bond-charge interaction term is enough to make the differences between the CDW and SDW correlation functions much less dramatic than when X=0. For X=t and fixed V2t (t is the uncorrelated hopping integral), there is an intermediate phase between a charge ordered phase and a phase corresponding to singly occupied sites, the nature of which we clarify: it is characterized by a succession of critical points, each of which corresponds to a different density of doubly occupied sites. We also find an unusual slowly decaying staggered spin-density correlation function, which is suggestive of some degree of ordering. No enhancement of pairing correlations was found for any X in the range examined. \textcopyright{} 1996 The American Physical Society.
By utilizing the twisted boundary conditions in the exact diagonalization method, we investigate the single-particle spectral function of the extended Peierls-Hubbard model at both half-filling and quarter filling. In one-dimensional (1D) interacting systems, the spin-charge separation can typically be identified in the single-particle spectral function by observing the distinct spinon and holon bands. At half filling, starting from the pure 1D Hubbard model with the on-site interaction U=10, we observe that the band structure indicative of the spin-charge separation gradually transitions to four individual bands as the Peierls instability δ increases. At U=10 and δ=0.2 where the spin-charge separation is still observable, increasing the nearest-neighbor interaction V can drive the system to a charge-density-wave (CDW) state when V≳U/2, without the obeservation of spinon and holon bands. At quarter-filling, on the other hand, the ground state of Peierls-Hubbard model manifests an antiferromagnetic Mott insulator in units of dimers. Increasing U results in only a very small gap in the single-particle spectrum because even for U=+∞, with the model transforming into a noninteracting half-filled dimerized tight-binding model, its gap determined by the Peierls instability δ remains small. Conversely, increasing V can effectively open the single-particle gap and make the spinon and holon bands more prominent.
To the conventional Peierls-Hubbard model, involving both on-site (U) and nearest-neighbor (V) Coulomb repulsions, we add ``off-diagonal'' terms, not expressible purely in terms of site densities, representing bond-bond (W) and mixed bond-site (X) electron-electron repulsive interactions involving nearest neighbors. We review earlier analyses of these interactions and discuss relative magnitudes of the parameters in applications to real materials. As a specific illustration, we investigate the effects of the off-diagonal W and X terms on dimerization in the one-dimensional, half-filled-band Peierls-Hubbard models, which have been widely applied to conjugated polymers (such as trans-polyacetylene) and to related quasi-one-dimensional charge-density wave (CDW) systems. Using both weak- and strong-coupling perturbation theory for large systems and exact diagonalizations of small systems, we investigate thoroughly the nature of the ground state of the model. For a broad range of the site-diagonal Hubbard parameters (U,V), including the values believed to be relevant to trans-polyacetylene, we find that the off-diagonal terms (W,X) initially enhance dimerization, thereby stabilizing the dimerized [or bond-order-wave (BOW)] ground state. For (unphysically) large values of W relative to U and V, dimerization is destroyed, and the BOW ground state goes over to a ferromagnetic ground state or a CDW ground state, depending on the relative sizes of U, V, and W. We conclude with a general discussion of the applicability of the Peierls-Hubbard models to quasi-one-dimensional materials, including the potential importance of the breaking of charge conjugation (``particle-hole'') symmetry by the X term.
We investigate two kinds of quantum phase transitions observed in the one-dimensional half-filled Peierls-Hubbard model with the next-nearest-neighbor hopping integral in the strong-coupling region $U\ensuremath{\gg}t$, ${t}^{\ensuremath{'}}$ [t ${(t}^{\ensuremath{'}})$, nearest- (next-nearest-) neighbor hopping; U, on-site Coulomb repulsion]. In the uniform case, with the help of the conformal field theory prediction, we numerically determine a phase boundary ${t}_{\mathrm{c}}^{\ensuremath{'}}(U/t)$ between the spin-fluid and the dimer states, where a bare coupling of the marginal operator vanishes and the low-energy and long-distance behaviors of the spin part are described by a free-boson model. To exhibit the conformal invariance of the systems on the phase boundary, a multiplet structure of the excitation spectrum of finite-size systems and a value of the central charge are also examined. The critical phenomenological aspect of the spin-Peierls transitions accompanied by the lattice dimerization is then argued for the systems on the phase boundary; the existence of logarithmic corrections to the power-law behaviors of the energy gain and the spin gap (i.e., the Cross-Fisher scaling law) are discussed.