We use a stochastic series-expansion quantum Monte Carlo method to study the phase diagram of the one-dimensional extended Hubbard model at half-filling for small to intermediate values of the on-site U and nearest-neighbor V repulsions. We confirm the existence of a novel, long-range-ordered bond-order-wave (BOW) phase recently predicted by Nakamura [J. Phys. Soc. Jpn. 68, 3123 (1999)] in a small region of the parameter space between the familiar charge-density-wave (CDW) state for $V\ensuremath{\gtrsim}U/2$ and the state with dominant spin-density-wave (SDW) fluctuations for $V\ensuremath{\lesssim}U/2.$ We discuss the nature of the transitions among these states and evaluate some of the critical exponents. Further, we determine accurately the position of the multicritical point, ${(U}_{m}{,V}_{m})=(4.7\ifmmode\pm\else\textpm\fi{}0.1,2.51\ifmmode\pm\else\textpm\fi{}0.04)$ (in energy units where the hopping integral is normalized to unity), above which the two continuous SDW-BOW-CDW transitions are replaced by one discontinuous (first-order) direct SDW-CDW transition. We also discuss the evolution of the CDW and BOW states upon hole doping. We find that in both cases the ground state is a Luther-Emery liquid, i.e., the spin gap remains but the charge gap existing at half-filling is immediately closed upon doping. The charge and bond-order correlations decay with distance r as ${r}^{\ensuremath{-}{K}_{\ensuremath{\rho}}},$ where ${K}_{\ensuremath{\rho}}$ is approximately
0.5 for the parameters we have considered. We also discuss advantages of using parallel tempering (or exchange Monte Carlo)---an extended ensemble method that we here combine with quantum Monte Carlo---in studies of quantum phase transitions.
We determine the ground-state phase diagram of the one-dimensional half-filled Hubbard model with on-site (nearest-neighbor) repulsive interaction U (V) and nearest-neighbor hopping t using the density-matrix renormalization group technique. Based on the results of the excitation gaps, Luttinger-liquid exponents, and bond-order-wave (BOW) order parameter, we confirm that the BOW phase appears in a substantial region between the charge-density-wave (CDW) and spin-density-wave phases. Each phase boundary is determined by multiple means and it allows us to make a cross-check on the validity of our estimations. We also find that the BOW-CDW transition changes from continuous to first order at the tricritical point (U(t),V(t)) approximately (5.89 t,3.10 t) and the BOW phase shrinks to zero at the critical end point (U(c),V(c)) approximately (9.25 t,4.76 t).
The quantum phases of one-dimensional spin
s=1/2 chains are discussed for models with two parameters, frustrating exchange
g=J2>0 between second neighbors and normalized nonfrustrating power-law exchange with exponent $\ensuremath{\alpha}$ and distance dependence ${r}^{\ensuremath{-}\ensuremath{\alpha}}$. The ground state (GS) at
g=0 has a long-range order (LRO) for $\ensuremath{\alpha}<2$ and long-range spin fluctuations for $\ensuremath{\alpha}>2$. The models conserve total spin
S=SA+SB, have singlet GS for any
g, $\ensuremath{\alpha}\ensuremath{\ge}0$ and decouple at
1/g=0 to linear Heisenberg antiferromagnets on sublattices
A and
B of odd- and even-numbered sites. Exact diagonalization of finite chains gives the sublattice spin $\ensuremath{\langle}{S}_{A}^{2}\ensuremath{\rangle}$, the magnetic gap
Em to the lowest triplet state, and the excitation ${E}_{\ensuremath{\sigma}}$ to the lowest singlet with opposite inversion symmetry to the GS. An analytical model that conserves sublattice spin has a first-order quantum transition at
gc=1/4ln2 from a GS with perfect LRO to a decoupled phase with
SA=SB=0 for $g\ensuremath{\ge}4/{\ensuremath{\pi}}^{2}$ and no correlation between spins in different sublattices. The model with $\ensuremath{\alpha}=1$ has a first-order transition to a decoupled phase that closely resembles the analytical model. The bond order wave (BOW) phase and continuous quantum phase transitions of finite models with $\ensuremath{\alpha}\ensuremath{\ge}2$ are discussed in terms of GS degeneracy where ${E}_{\ensuremath{\sigma}}(g)=0$, excited state degeneracy where ${E}_{\ensuremath{\sigma}}(g)={E}_{m}(g)$, and $\ensuremath{\langle}{S}_{A}^{2}\ensuremath{\rangle}$. The decoupled phase at large frustration has nondegenerate GS for any exponent $\ensuremath{\alpha}$ and excited states related to sublattice excitations.
The static structure factor S(q) of frustrated spin-1/2 chains with isotropic exchange and a singlet ground state (GS) diverges at wave vector q_m when the GS has quasi-long-range order (QLRO) with periodicity 2π/q_m but S(q_m) is finite in bond-order-wave (BOW) phases with finite-range spin correlations. Exact diagonalization and density matrix renormalization group (DMRG) calculations of S(q) indicate a decoupled phase with QLRO and q_m = π/2 in chains with large antiferromagnetic exchange between second neighbors. S(q_m) identifies quantum phase transitions based on GS spin correlations.
By using the density matrix renormalization group technique, the phase diagram of the half-filling extended Hubbard model is investigated. The conventional order parameter, the two-site entanglement entropy, and the block-block entanglement entropy are analyzed in detail. According to the numerical results, in the weak coupling region, an intermediate bond-order-wave (BOW) phase is shown to exist indeed between the charge-density-wave (CDW) and the spin-density-wave (SDW) phases. The critical phase transition points are determined by the singularity of the first order derivative of two-site entanglement entropy. In strong coupling region, a direct phase transition occurs from SDW phase to CDW phase, and shows discontinuous (first order) character accompanied with energy level crossing. The numerical results support the phase diagram proposed previously by some authors Sengupta et al., (2002); Sandvik et al., (2004); Zhang, (2004). Therefore, the quantum entanglement is a sensitive tool to describe quantum phase transitions in strongly correlated electron systems.
In this paper, by using the level spectroscopy method and bosonization theory, we discuss the evolution of the bond-order-wave (BOW) phase in a one-dimensional half-filled extended Hubbard model with the on-site Coulomb repulsion U as well as the inter-site Coulomb repulsion V and antiferromagnetic exchange J . After clarifying the generic phase diagrams in three limiting cases with one of the parameters being fixed at zero individually, we find that the BOW phase in the U-V phase diagram is initially enlarged as J increases from zero but is eventually suppressed as J increases further in the strong-coupling regime. A three-dimensional phase diagram is suggested where the BOW phase exists in an extended region separated from the spin-density-wave and charge-density-wave phases.