Reply to the Comment by Sandvik, Sengupta, and Campbell on ``Ground State Phase Diagram of a Half-Filled One-Dimensional Extended Hubbard Model''

In their Comment (see cond-mat/0301237), Sandvik, Sengupta, and Campbell present some numerical evidences to support the existence of an extended bond-order-wave (BOW) phase at couplings (U,V) weaker than a tricritical point (U_t,V_t) in the ground state phase diagram of the one-dimensional half-filled U-V Hubbard model. They claim that their results do not agree with the phase diagram proposed in my Letter (cond-mat/0204244), which shows a BOW phase for couplings stronger than the critical point only. However, I argue here that their results are not conclusive and do not refute the phase diagram described in the Letter.Comment: 1 page, published versio

Charge-density-wave formation in the Edwards fermion-boson model at one-third band filling

We examine the ground-state properties of the one-dimensional Edwards spinless fermion transport model by means of large-scale density-matrix renormalization-group calculations. Determining the single-particle gap and the Tomonaga-Luttinger liquid parameter ($K_\rho$) at zero temperature, we prove the existence of a metal-to-insulator quantum phase transition at one-third band filling. The insulator---established by strong correlation in the background medium---typifies a charge density wave (CDW) that is commensurate with the band filling. $K_\rho=2/9$ is very small at the quantum critical point, and becomes $K_\rho^{\rm CDW}=1/9$ in the infinitesimally doped three-period CDW, as predicted by the bosonization approach.Comment: 6 pages, 3 figures, contributions to SCES 201

Optical conductivity of the one-dimensional dimerized Hubbard model at quarter filling

We investigate the optical conductivity in the Mott insulating phase of the one-dimensional extended Hubbard model with alternating hopping terms (dimerization) at quarter band filling. Optical spectra are calculated for the various parameter regimes using the dynamical density-matrix renormalization group method. The study of limiting cases allows us to explain the various structures found numerically in the optical conductivity of this model. Our calculations show that the dimerization and the nearest-neighbor repulsion determine the main features of the spectrum. The on-site repulsion plays only a secondary role. We discuss the consequences of our results for the theory of the optical conductivity in the Bechgaard salts.Comment: 11 pages and 12 figure

Stripe formation in doped Hubbard ladders

We investigate the formation of stripes in 7\chunks \times 6 Hubbard ladders with 4\chunks holes doped away from half filling using the density-matrix renormalization group (DMRG) method. A parallelized code allows us to keep enough density-matrix eigenstates (up to $m=8000$) and to study sufficiently large systems (with up to 7\chunks = 21 rungs) to extrapolate the stripe amplitude to the limits of vanishing DMRG truncation error and infinitely long ladders. Our work gives strong evidence that stripes exist in the ground state for strong coupling ($U=12t$) but that the structures found in the hole density at weaker coupling ($U=3t$) are an artifact of the DMRG approach.Comment: 6 pages, 6 Figure

Parallelization Strategies for Density Matrix Renormalization Group Algorithms on Shared-Memory Systems

Shared-memory parallelization (SMP) strategies for density matrix renormalization group (DMRG) algorithms enable the treatment of complex systems in solid state physics. We present two different approaches by which parallelization of the standard DMRG algorithm can be accomplished in an efficient way. The methods are illustrated with DMRG calculations of the two-dimensional Hubbard model and the one-dimensional Holstein-Hubbard model on contemporary SMP architectures. The parallelized code shows good scalability up to at least eight processors and allows us to solve problems which exceed the capability of sequential DMRG calculations.Comment: 18 pages, 9 figure

Spectral function of the one-dimensional Hubbard model away from half filling

We calculate the photoemission spectral function of the one-dimensional Hubbard model away from half filling using the dynamical density matrix renormalization group method. An approach for calculating momentum-dependent quantities in finite open chains is presented. Comparison with exact Bethe Ansatz results demonstrates the unprecedented accuracy of our method. Our results show that the photoemission spectrum of the quasi-one-dimensional conductor TTF-TCNQ provides evidence for spin-charge separation on the scale of the conduction band width.Comment: REVTEX, 4 pages including 4 EPS figures (changed); correct chemical potential used to define excitation energies in figures and tex

Hole doped Hubbard ladders

The formation of stripes in six-leg Hubbard ladders with cylindrical boundary conditions is investigated for two different hole dopings, where the amplitude of the hole density modulation is determined in the limits of vanishing DMRG truncation errors and infinitely long ladders. The results give strong evidence that stripes exist in the ground state of these systems for strong but not for weak Hubbard couplings. The doping dependence of these findings is analysed.Comment: 2 pages, 2 figures, submitted to SCES0