We report quantum Monte Carlo (stochastic series expansion) results for the transition from a Mott insulator to a dimerized Peierls insulating state in a half-filled, one-dimensional extended Hubbard model coupled to optical bond phonons. Using electron-electron (e-e) interaction parameters corresponding approximately to polyacetylene, we show that the Mott-Peierls transition occurs at a finite value of the electron-phonon (e-ph) coupling. We discuss several different criteria for detecting the transition and show that they give consistent results. We calculate the critical e-ph coupling as a function of the bare phonon frequency and also investigate the sensitivity of the critical coupling to the strength of the e-e interaction. In the limit of strong e-e couplings, we map the model to a spin-Peierls chain and compare the phase boundary with previous results for the spin-Peierls transition. We point out effects of a nonlinear spin-phonon coupling neglected in the mapping to the spin-Peierls model.
We report quantum Monte Carlo (stochastic series expansion) results for the transition from a Mott insulator to a dimerized Peierls insulating state in a half-filled, one-dimensional extended Hubbard model coupled to optical bond phonons. Using electron-electron (e-e) interaction parameters corresponding approximately to polyacetylene, we show that the Mott-Peierls transition occurs at a finite value of the electron-phonon (e-ph) coupling. We discuss several different criteria for detecting the transition and show that they give consistent results. We calculate the critical e-ph coupling as a function of the bare phonon frequency and also investigate the sensitivity of the critical coupling to the strength of the e-e interaction. In the limit of strong e-e couplings, we map the model to a spin-Peierls chain and compare the phase boundary with previous results for the spin-Peierls transition. We point out effects of a nonlinear spin-phonon coupling neglected in the mapping to the spin-Peierls 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."
The zero and finite temperature spin-Peierls transitions in a quasi-one-dimensional spin-21 Heisenberg model coupled to adiabatic bond phonons is investigated using the stochastic series expansion (SSE) quantum Monte Carlo (QMC) method. The quantum phase transition from a gapless N\'eel state to a spin-gapped Peierls state is studied in the parameter space spanned by spatial anisotropy, interchain coupling strength, and spin-lattice coupling strength. It is found that for any finite interchain coupling, the transition to a dimerized Peierls ground state only occurs when the spin-lattice coupling exceeds a finite, nonzero critical value. This is in contrast to the pure 1D model (zero interchain coupling), where adiabatic/classical phonons lead to a dimerized ground state for any nonzero spin-phonon interaction. The phase diagram in the parameter space shows that for a strong interchain coupling, the relation between the interchain coupling and the critical value of the spin-phonon interaction is linear whereas for weak interchain coupling, this behavior is found to have a natural logarithmlike relation. No region was found to have a long range magnetic order and dimerization occurring simultaneously. Instead, the N\'eel state order vanishes simultaneously with the setting in of the spin-Peierls state. For the thermal phase transition, a continuous heat capacity with a peak at the critical temperature Tc shows a second order phase transition. The variation of the equilibrium bond length distortion ${\ensuremath{\delta}}_{eq}$ with temperature showed a power law relation which decayed to zero as the temperature was increased to Tc, indicating a continuous transition from the dimerized phase to a paramagnetic phase with uniform bond length and zero antiferromagnetic susceptibility.
The renormalization group technique is applied to one-dimensional electron-phonon Hubbard models at half filling and zero temperature. For the Holstein-Hubbard model, the results of one-loop calculations are congruent with the phase diagram obtained by quantum Monte Carlo simulations in the (U,gph) plane for the phonon-mediated interaction gph and the Coulomb interaction U. The incursion of an intermediate phase between a fully gapped charge-density-wave state and a Mott antiferromagnet is supported along with the growth of its size with the molecular phonon frequency ${\ensuremath{\omega}}_{0}$. We find additional phases enfolding the base boundary of the intermediate phase. A Luttinger liquid line is found below some critical ${U}^{*}\ensuremath{\approx}{g}_{\mathrm{ph}}^{*}$, followed at larger $U\ensuremath{\sim}{g}_{\mathrm{ph}}$ by a narrow region of bond-order-wave ordering which is either charge or spin gapped depending on U. For the Peierls-Hubbard model, the region of the (U,gph) plane with a fully gapped Peierls-bond-order-wave state shows a growing domination over the Mott gapped antiferromagnet as the Debye frequency ${\ensuremath{\omega}}_{D}$ decreases. A power-law dependence ${g}_{\mathrm{ph}}\ensuremath{\sim}{U}^{2\ensuremath{\eta}}$ is found to map out the boundary between the two phases, whose exponent is in good agreement with the existing quantum Monte Carlo simulations performed when a finite nearest-neighbor repulsion term V is added to the Hubbard interaction.
The one-dimensional extended Peierls-Hubbard model is studied at several band fillings using the density-matrix renormalization-group method. Results show that the ground state evolves from a Mott-Peierls insulator, with a correlation gap at half-filling, to a soliton lattice with a small band gap away from half-filling. It is also confirmed that the ground state of the Peierls-Hubbard model undergoes a transition to a metallic state at finite doping. These results show that electronic correlation effects should be taken into account in theoretical studies of doped polyacetylene. They also show that a Mott-Peierls theory could explain the insulator-metal transition observed in this material.
We present a numerical study of the Hubbard-Holstein model in one dimension at half filling, including finite-frequency quantum phonons. At half filling, the effects of the electron-phonon and electron-electron interactions compete with the Holstein phonon coupling acting as an effective negative Hubbard on-site interaction U that promotes on-site electron pairs and a Peierls charge-density wave state. Most previous work on this model has assumed that only Peierls or Mott phases are possible at half filling. However, there has been speculation that a third metallic phase exists between the Peierls and Mott phases. We confirm the intermediate phase, and show that the Luttinger liquid correlation exponent K(rho) >1 in this region, indicating dominant superconducting pair correlations. We explore the full phase diagram as a function of Hubbard U, phonon coupling constant, and phonon frequency.
The density matrix renormalization group method is used to investigate the Peierls transition for the extended Hubbard model coupled to quantized phonons. Following our earlier work on spin-Peierls systems, we use a phonon spectrum that interpolates between a gapped, dispersionless (Einstein) limit and a gapless, dispersive (Debye) limit to investigate the entire frequency range. A variety of theoretical probes are used to determine the quantum phase transition, including energy gap crossing, a finite-size scaling analysis, and bipartite quantum entanglement. All these probes indicate that a transition of Berezinskii-Kosterlitz-Thouless type is observed at a nonzero electron-phonon coupling gc for a nonvanishing electron-electron interaction. An extrapolation from the Einstein limit to the Debye limit is accompanied by an increase in gc for a fixed optical ($q=\ensuremath{\pi}$) phonon gap. We therefore conclude that the dimerized ground state is more unstable with respect to Debye phonons, with the introduction of phonon dispersion renormalizing the effective electron-lattice coupling for the Peierls-active mode. By varying the Coulomb interaction U, we observe a generalized Peierls transition, intermediate between the uncorrelated (U=0) and spin-Peierls ($U\ensuremath{\rightarrow}\ensuremath{\infty}$) limits, where U is the Hubbard Coulomb parameter. Using the extended Hubbard model with Debye phonons, we investigate the Peierls transition in trans-polyacetylene and show that the transition is close to the critical regime.